Belka z odporem gruntu, dwie sily skupione i moment - CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne


Start :· Galeria :· Forum :· Download :· Artykuły :· Newsy :· Linki :· Zasoby wiedzy :· Kontakt :· Zasady	Stosujemy pliki cookies więcej...	11/23/2024 09:07:40

Reklama

Nawigacja

Galeria

Forum

Download

Artykuły

FAQ

Linki

Zasoby wiedzy

Kontakt

Szukaj

Losowa Fotka

Warunki brzegowe - siłowe
Kompostownik

Użytkownicy Online

CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne
WITAMY:
adrian24024
jako nowego użytkownika.

Zarejestrowanch Uzytkowników: 1400

Super Administratorzy: 1
Administratorzy: 1
Użytkownicy: 1398

Użytkownicy Online:

arepaj	25 tygodni
ndv	30 tygodni
piotr	46 tygodni
damian14100	79 tygodni
BE-FEA	81 tygodni

Gości Online: 31

Twoje IP to: 18.191.200.47

Kategorie Forum	53
Wątki na Forum	1268
Posty na Forum	5665
Komentarzy	83
Newsy	142
Artykuły	10
Ściągniętych plików	62
Kategorie Downloads	4
Ciekawe Strony	13
Postów w Shoutbox	413

Przykłady

Po drugiej iteracji
Problem Oli

model1
Model Felgi

Pola ciągnięń
Fajne gify

Skan tomograficzny zęba
Ząb

Przewijanie
Fajne gify

Partnerzy

Strony wspierającej nas firmy: Leasing,faktoring, kredyt gdy szukasz gotówki.
Opony to tanie opony letnie oraz katalog
zawierający firmy transportowe
Strony naszych przyjaciół:

Strony sponsorów:
Zawsze się znajdzie forum transportowe na giełdzie ładunków/ na stronie giełdy ładunków
ciekawe i miłe sercu rzeczy to: fajne fotki dzieci na stronie mojebaby.pl lub też warte obejrzenia zwierzęta domowe ze strony mojzwierz.pl

-->../lll

Ostatnie Komentarze

marek_ac dnia 11.01.2017 08:14

Demo download (nowy link): (poprawiono patrz: Linki-> Przyjaciele -> Vanadis)

Zobacz Komentarze Artykułów

marek_ac dnia 21.04.2015 05:53

... mała poprawka: double wprowadzono w wersji 1.3

Zobacz Komentarze Newsów

marek_ac dnia 20.04.2015 18:59

Jeżeli do oprogramowania inżynierskiego można podpiąć własne procedury to NVIDIA udostępnia darmowy driver NVCC za pomocą którego można kompilować kod napisany w C. Można to robić z Visual Studio (od wersji 10) pod Windows lub używając gcc pod Linuxem.
Wówczas oprogramowanie optymalizuje się pod posiadaną kartę graficzną. Warto mieć kartę o CUDA capability większym niż 2.0 ze względu na obsługę typu double.

Zobacz Komentarze Newsów

marek_ac dnia 14.04.2015 00:51

Witam,

Na dysku:

TUTAJ

zamieściłem demo (wersję windowsową) modelu rozprzestrzeniania się zanieczyszczeń.

W pliku vanadis.txt zawarty jest krótki opis modelu.

Marek

Zobacz Komentarze Artykułów

steleri dnia 17.03.2015 21:54

Czy ktoś z forum korzysta może z Elmera? Bardzo ciekawi mnie jak to działa. Bo z tego tutka wynika, że jest OK, a tymczasem ja próbowałem... bezskutecznie

Zobacz Komentarze Newsów

Ostatnie aktywne Tematy

rss1

rss2

Forum	Wątek	Odpowiedzi	Ostatni Post
POMOCNA DŁOń - pomagamy studentom rozwiązywać zadania	Analiza statyczna wózka wagonowego	3	BE-FEA	02/28/2022 18:26
Jestem kompletnie zielony	ABAQUS PROBLEM	2	Euklides_PL	01/25/2022 20:24
Ogłoszenia i newsy (branża CAE)	OBLICZENIA WYTRZYMAŁOŚCIOWE MES I RYSUNEK TECHNICZNY	0	1pkm	09/24/2021 01:50
Jestem kompletnie zielony	ABAQUS BŁĄD	1	BE-FEA	09/15/2021 19:44
PRACA w CAE, oferty "biznesowe"	Pomoc projektowa - Ansys fluent	3	BE-FEA	09/15/2021 19:32
Jestem kompletnie zielony	ABAQUS - ANALYTICAL RIGID	2	Euklides_PL	08/30/2021 21:20
Jestem kompletnie zielony	WSTĘPNE NAPRĘZENIA ŚCISKAJĄCE/SHRINK FIT ABAQUS	1	BE-FEA	06/03/2021 14:14
PRACA w CAE, oferty "biznesowe"	Szkolenia komercyjne	3	Adrian28	05/30/2021 11:57
Sprawy studenckie	OPTYMALIZACJA POWIERZCHNI ABAQUS	5	BE-FEA	03/23/2021 23:17
Sprawy studenckie	comsol heat transfer in porous media	1	BE-FEA	03/09/2021 17:55
Sprawy studenckie	ABAQUS-WARUNKI BRZEGOWE	1	BE-FEA	02/15/2021 00:01
Sprawy studenckie	Potrzebna pomoc - Femap	2	jasiu	01/27/2021 18:25
Sprawy studenckie	SolidWorks Simulation - podnoszenie zbiornika	3	BE-FEA	11/24/2020 14:33
Teoria MES	Kwantowy MES	1	BE-FEA	11/13/2020 14:32
Jestem kompletnie zielony	DEFORM	0	CwaniakzPKSU	11/05/2020 08:46

Zobacz Temat

CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne | Pozostałe wątki | POMOCNA DŁOń - pomagamy studentom rozwiązywać zadania

Autor

Belka z odporem gruntu, dwie sily skupione i moment

m10
Użytkownik

Postów: 1
Data rejestracji: 20.08.11

Dodane dnia 08/20/2011 15:55

Na poczatku chce wszystkich pozdrowic i przywitac.

Problem inzynieryjny do obliczenia.

Wymagany program do obliczenia SCILAB

Schemat belki

Po obliczeniu w programie belka na szybko wychodzi dobrze. to znaczy:

Momenty

TNACE I REAKCJE

UGIECIA

Jednak po wrzuceniu tego do SCILABA, program SCILAB wypuszcza bledy w obliczeniach na przyklad wiadomo ze w miejscu przylozenia momentu moment jest zerowy, a na wykresie wychodzi wartosc inna.

Wydruk z programu SCILAB

Katy obrotu

------------------------------

Momenty

-----------------------------

Odpor gruntu

----------------------------

Tnace

----------------------------

Ugiecia

-------------------------------

Wydruk z SCILABA

//belka sekcja dane
E=30000000;//mod. younGa [kpa]
BA=0.8//szerokosc belki na pierwszym odcinku
HA=1//wysokosc na pierszym przekroju belki
BB=1.2//szerokosc przekroju sprezystego
HB=0.6//wyskosc przekroju belki sprezystego
LP=245//liczba punktow na jaka dzielimy belke
ZAKRES=245//liczba punktow na jaka dzielimy belke
C=30000//liczba modul sprzystosci podloza
P1=5;//sila P1 o wartosci 5kN
P2=12;//sila P2 o wartosci 12kN
M=3;//wartosc momentu M wyrazona w wartosci 3kN
Q1=10;// wartosc obciazenia rozlozonego o wartosci Q1=10 kN
Q2=15;// wartosc obciazenia rozlozonego trojakatnego o wartosci Q2=15 kN
lambda=0.1;// dla warunku brzegowego stala Lambda

//obliczenia pomocnicze moment bezwladnosci
JA=BA*HA^3/12;
JB=BB*HB^3/12;

//obliczeniapomocnicze
QP1=P1/lambda;
QP2=P2/lambda;

//modul sprezystosci podloza
KA=C*BA;//przekroj pierwszy
KB=C*BB;//przekroj drugi

//stworzenie wektora momentu bezwladnosci
VJ=JA*ones(ZAKRES,1);
VJ(63)=2*JA*JB/(JA+JB);
VJ(123)=0;
VJ(63:ZAKRES)=JB;
//wydruk wartosci momentu bezwladnosci

//wektor sztywnosci podloza sprezystego
VK=KB*ones(ZAKRES,1);//wektor 245 elementowy
VK(1:2)=0;
VK(3)=0.5*KA;
VK(4:62)=KA;
VK(63)=0.5*(KA+KB);
VK(243)=0.5*KB;
VK(244:245)=0;

//wektor obciazen
VQ=zeros(ZAKRES,1);
VQ(33:123)=Q1;
VQ(63)=VQ(63)+P1;
VQ(123)=0.5*(Q1+Q2);
VQ(123:183)=(linspace(Q2,0,61))';
VQ(213)=QP2;

[(1:ZAKRES)'VQ]

//maciez glowna G;
G=zeros(ZAKRES, ZAKRES);

//wektor pomocniczy
WWJ=ones(ZAKRES, 1);

//macierz Glowna

for i=3:(ZAKRES-2) do
G(i, i-2)=VJ(i-1);
G(i, i-1)=-2*(VJ(i-1)+VJ(i));
G(i, i)=VJ(i-1)+4*VJ(i)+VJ(i+1)+VK(i)*lambda^4*WWJ(i)/E;
G(i, i+1)=-2*(VJ(i)+VJ(i+1));
G(i, i+2)=VJ(i+1);
end

//warunki brzegowe
Q=VQ*lambda^4/E;
G(1,2)=1;G(1,3)=-2;G(1,4)=1;Q(1)=0;
G(2,1)=-1;G(2,2)=2;G(2,4)=-2;G(2,5)=1;Q(2)=2*P1*lambda^3/(E*JA);
G($-1,243)=1;Q($-1)=0;
G($,242)=-1;G($,244)=1;Q($)=0;
G(63,:)=0;G(:,63)=0;G(63,63)=1;Q(63)=0;

//Rozwi¹zanie uk³adu równañ
//wektor ugiec
W=inv(G)*Q;

//Obliczenie pozosta³ych wielkoœci
fi=zeros(ZAKRES, 1);
M=zeros(fi);
T=zeros(fi);
R=zeros(fi);

for i=3:(ZAKRES-2) do
fi(i)=(-W(i-1)+W(i+1))/(2*lambda);
M(i)=-E*VJ(i)*(W(i-1)-2*W(i)+W(i+1))/(lambda^2);
T(i)=-E*(-VJ(i-1)*(W(i-2)-2*W(i-1)+W(i))+VJ(i+1)*(W(i)-2*W(i+1)+W(i+2)))/(2*lambda^3);
R(i)=VK(i)*W(i)*WWJ(i);
end

fi=fi(3:ZAKRES-2);
M=M(3:ZAKRES-2);
T=T(3:ZAKRES-2);
R=R(3:ZAKRES-2);

//Usuniecie wszystkich otwartych okien Graficznych:
xdel(winsid());

//Wykres ugiec:
xinit();
plot2d((0:(ZAKRES-5))*lambda, W(3:(ZAKRES-2)), [1], rect=[0, 1.1*min(W), (ZAKRES-5)*lambda, 1.1*max(W)]);
plot2d((0:(ZAKRES-5))*lambda, zeros(1,(ZAKRES-4)), [1]);
xname("Ugiecia");
xtitle('Ugiecia', 'x[m]', 'w[m]');
wy1=gcf();
wy1.children.axes_reverse = ["off","on","off"];
wy1.children.children.children(2).thickness=3;

//Wykres katów obrotu:
xinit();
plot2d((0:(ZAKRES-5))*lambda, fi, [1], rect=[0, 1.1*min(fi), (ZAKRES-5)*lambda, 1.1*max(fi)]);
plot2d((0:(ZAKRES-5))*lambda, zeros(1,(ZAKRES-4)), [1]);
xname("Katy obrotu");
xtitle('Katy obrotu', 'x[m]', 'fi[rad]');
wy2=gcf();
wy2.children.axes_reverse = ["off","on","off"];
wy2.children.children.children(2).thickness=3;

//Wykres momentów:
xinit();
plot2d((0:(ZAKRES-5))*lambda, M, [1], rect=[0, 1.1*min(M), (ZAKRES-5)*lambda, 1.1*max(M)]);
plot2d((0:(ZAKRES-5))*lambda, zeros(1,(ZAKRES-4)), [1]);
xname("Momenty");
xtitle('Momenty', 'x[m]', 'M[kNm]');
wy3=gcf();
wy3.children.axes_reverse = ["off","on","off"];
wy3.children.children.children(2).thickness=3;

//Wykres tnacych:
xinit();
plot2d((0:(ZAKRES-5))*lambda, T, [1], rect=[0, 1.1*min(T), (ZAKRES-5)*lambda, 1.1*max(T)]);
plot2d((0:(ZAKRES-5))*lambda, zeros(1,(ZAKRES-4)), [1]);
xname("Tnace");
xtitle('Tnace', 'x[m]', 'T[kN]');
wy4=gcf();
wy4.children.axes_reverse = ["off","on","off"];
wy4.children.children.children(2).thickness=3;

//Wykres odporu Gruntu:
xinit();
plot2d((0:(ZAKRES-5))*lambda, R, [1], rect=[0, 1.1*min(R), (ZAKRES-5)*lambda, 1.1*max(R)]);
plot2d((0:(ZAKRES-5))*lambda, zeros(1,(ZAKRES-4)), [1]);
xname("Odpór Gruntu");
xtitle('Odpór Gruntu', 'x[m]', 'r[kN/m]');
wy5=gcf();
wy5.children.axes_reverse = ["off","on","off"];
wy5.children.children.children(2).thickness=3;

//Zapisanie wyników do plików:
wyn=[W(3:(ZAKRES-2)), fi, M, T, R];
savematfile('wyniki.txt', 'wyn', '-ascii');

-----------------------------------------

Moje pytanie: czy ktos starszy z kolegow przejzalby kod programu SCILAB i ewentualnie skorygowal sekcje warunki brzegowe oraz maciez glowna.

Pozdrawiam

Edytowane przez m10 dnia 08/20/2011 15:55

Autor

RE: Belka z odporem gruntu, dwie sily skupione i moment

piotr
Super Administrator

Postów: 1432
Data rejestracji: 17.05.07

Dodane dnia 08/22/2011 20:44

Trochę długa prośba jak na końcówkę wakacji :f
Ale czy macierz sztywności nie powinna być symetryczna?

Skocz do Forum:

Reklama

Szukaj

Logowanie

Nie jesteś jeszcze naszym Użytkownikiem?
Kilknij TUTAJ żeby się zarejestrować.

Zapomniane hasło?
Wyślemy nowe, kliknij TUTAJ.

Ankieta dla Użytkowników

Reklama

Shoutbox

Tylko zalogowani mogą dodawać posty w shoutboksie.

Bulix
26.01.2022
Witam, szukam osoby która ogarnia program FEMM.

damian14100
25.01.2022
Witam, Czy znajduje się na forum osoba która ma doświadczenie w obliczeniach wytrzymałościowych wózków wagonowych ?

BE-FEA
31.03.2021
Nie pracuję w Deform, ale może coś podpowiem na zasadzie analogii do innych programów MES. Napisz jeśli nadal aktualne.

daniel8894
31.03.2021
Czy pomógłby ktoś w zrobieniu symulacji procesu skrawania w Deform 3 d

BE-FEA
09.03.2021
Postaram się pomóc, zapraszam do kontaktu jeśli nadal aktualne.

kinia22
09.03.2021
witam, potrzebuje kogoś kto pomoże w przygotowaniu projektu w programie COMSOL - przeplyw ciepła przez ośrodek porowaty!

BE-FEA
02.04.2020
Przerobiło mi jedną literę na emotkę. Chodzi o przycisk new thread

BE-FEA
02.04.2020
Żeby dodać nowy temat wystarczy wejść w odpowiedni dział i kliknąć mały przycisk "

ew thread" po prawej stronie.

KrzywaOHIO
30.03.2020
Jak dodać nowy wątek na forum?

Konrad96
02.12.2019
Witam, czy jest tu ktoś kto dobrze zna się na programie Robot struktural i pomógłby mi zamodelować belkę żelbetową podobną do belki użytej w badaniu? Proszę o kontakt jeśli ktoś chciałby pomóc

Archiwum

Wydarzenia

Listopad 2024

Pon	Wt	�r	Czw	Pi	So	Nie
				1	2	3
4	5	6	7	8	9	10
11	12	13	14	15	16	17
18	19	20	21	22	23	24
25	26	27	28	29	30

	115848793 Unikalnych wizyt

	Powered by Peter-Fusion © 2009 - Modified Aztec Theme