MES,CAD,CAE,FEA
News Search:
StartGaleriaForumDownloadArtykułyNewsyLinkiZasoby wiedzyKontaktZasady Stosujemy pliki cookies więcej...12/06/2019 12:36:03
Reklama
Nawigacja
Galeria
Forum
Download
Artykuły
FAQ
Linki
Zasoby wiedzy
Kontakt
Szukaj
Losowa Fotka
Wyniki obliczeń - naprężenia redukowane -konstrukcja odkształcona
Wyniki obliczeń - naprężenia redukowane -konstrukcja odkształcona
Kompostownik
Użytkownicy Online
CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne
WITAMY:
Konrad96
jako nowego użytkownika.

Zarejestrowanch Uzytkowników: 1365

Super Administratorzy: 1
Administratorzy: 1
Użytkownicy: 1363

Użytkownicy Online:

BE-FEA04:06:16
Konrad96 2 dni
piotrmadej 1 tydzień
piotr 3 tygodni
maq9633 3 tygodni

Gości Online: 11

Twoje IP to: 18.207.238.169

Kategorie Forum 53
Wątki na Forum 1243
Posty na Forum 5588
Komentarzy 83
Newsy 141
Artykuły 10
Ściągniętych plików 62
Kategorie Downloads 4
Ciekawe Strony 13
Postów w Shoutbox 404
Przykłady
deformacja konstrukcji

deformacja konstrukcji

Obliczenia trawersy
Klipsik

Klipsik

Fajne gify
Obciążenia

Obciążenia

Obliczenia trawersy
Przykładowe warunki kinematyczne brzegowe

Przykładowe warunki kinematyczne brzegowe

Obliczenia trawersy
Analiza źródeł awarii łyżki koparkowej -

Analiza źródeł awarii łyżki koparkowej -

Awaria łyżki koparki
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ół:

Darmowe arkusz Excela
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
Sprawy studenckie Praca Magisterska Ansys 2 piotr 11/14/2019 17:08
Jestem kompletnie zielony Współczynnik konwekcji 1 BE-FEA 09/30/2019 19:14
Wszystko o CAE, co nie pasuje do powyższych forum Bezpłatne webinarium CAE - jak działa program do symulacj... 0 MatLit 08/28/2019 11:20
Teoria MES Łączenie elementów skończonych różnego typu 1 BE-FEA 06/10/2019 21:59
Teoria MES Macierz pochodnych funkcji kształtu [B] 14 BE-FEA 06/05/2019 14:44
Jestem kompletnie zielony Jak połączyć wolne powierzchnie w Patranie 21 BE-FEA 05/31/2019 17:52
Teoria MES Rodzaje solverów MES 14 BE-FEA 05/24/2019 19:12
Jestem kompletnie zielony Problem z importem plików do Patrana 6 Solid-ny amator 04/23/2019 23:35
Jestem kompletnie zielony [ANSYS] Wytłaczanie uszczelek 1 Rogal 04/08/2019 11:26
Sprawy studenckie Odpłatna pomoc 2 piotr 03/28/2019 21:56
Sprawy studenckie Analiza strukturalna - skrzydło szybowca 4 BE-FEA 01/27/2019 15:04
Sprawy Abaqusa Sztywna cześć 2 BE-FEA 01/18/2019 12:46
Sprawy studenckie Pomoc w napisaniu prostego programu 2 pvalue 01/17/2019 17:51
Konkretne problemy Profil aluminiowy + wkładka termiczna 8 BE-FEA 12/13/2018 12:51
PRACA w CAE, oferty "biznesowe" POSZUKUJĘ STAŻU - CAE 0 mesuser94 11/07/2018 19:34
Zobacz Temat
CAD, MES -metoda elementów skończonych,obliczenia inżynierskie i metody numeryczne | Oprogramowanie CAE | Wszystko o CAE, co nie pasuje do powyższych forum
Autor MES w MATLAB - trochę nietypowo
popey
Użytkownik

Postów: 3
Data rejestracji: 26.06.08
Dodane dnia 06/26/2008 13:10
Szanowni Państwo

Na początku chciałbym się przywitać.
Założyłem ten wątek w dziale off-topic ponieważ wydaje mi się, że nie jest on zbytnio związany ani z programami CAE ani z prawdziwą inżynierią. Dlaczego? Otóż mój problem to prawdopodobnie " podstawy podstaw " MES - więc nie inżynieria. No i działam w Matlabie - więc nie specjalistyczne programy do MES
No, ale do rzeczy.

Wcześniej to samo pytanie zadałem na forum MATLAB, dopiero później odnalazłem Państwa stronę.

Otrzymałem zadanie stworzenia skryptu w Matlabie.
Problem wygląda następująco: istnieje wał (model Jeffcotta - wał bez masy, giętki, na nim umieszczony dysk o pewnej masie. Wał podparty na końcach z podporami o zmiennej sztywości - sztywność tę da się regulować).

Obliczenia matematyczne dały rozwiązanie postaci:

Mq" + Cq' + Kq = P

gdzie M to macierz mas, C - macierz tłumień i K - macierz sztywności
P zaś jest wektorem wymuszeń.

Jak widać - (chyba) typowy prosty problem.
Teraz jednak mam problem z napisaniem programu do tego. W związku z tym miałbym do Państwa prośbę: nie chodzi mi o podanie konkretnych rozwiązań, nie chcę, by ktokolwiek robił coś za mnie, lecz jedynie kilku wskazówek. Zapoznałem się z pracami dostępnymi w Internecie, lecz tamte rozwiązania nie korzystają z takiej postaci rozwiązania jakiej ja używam. Są tam (tak mi się wydaje, ponieważ materiały są po angielsku i to w bardzo technicznym jego wydaniu) bardzo skomplikowane przykłady z całkami itd., a mi się wydaje, że to wystarczą jakieś proste przekształcenia, żeby wykreślić przemieszczenia q.

Problem u mnie jest taki, że macierz K będzie się zmieniać, ale to później. Najpierw chciałbym napisać program rozwiązujący to równanie i wykreślający przesunięcia końców wału (zataczają one elipsę). Prosiłbym o kilka słów kluczowych - na co zwrócić uwagę lub od czego w ogóle zacząć? Z MES zetknąłem się dopiero " przed chwilą " - zrobiłem analizę tego wału, wyznaczyłem wsp. macierzy i teraz pozostaje już tylko program, ale zupełnie nie mam na niego pomysłu (zbyt mała wiedza, jak przypuszczam)

Byłbym bardzo wdzięczny za każdą, najmniejszą choćby wskazówkę.
Z góry dziękuję, pozdrawiam.
Edytowane przez popey dnia 06/26/2008 13:11
Autor RE: MES w MATLAB - trochę nietypowo
mnichv
Użytkownik

Postów: 90
Data rejestracji: 20.03.08
Dodane dnia 06/26/2008 16:39
function uklad3mas_mac1
%--------------------------------------------------------------------------
% k13
% |----------/\/\/\---------|
% | k12 k23 |
% O---/\/\/\---O---/\/\/\---O
% m1 m2 m3


m1=1; m2=1; m3=1;
M=diag([m1 m2 m3]) %<-- macierz bezwladnosci
k12=50; k23=5; k31=50;
k21=k12; k32=k23; k13=k31;
K=[(k12+k13) -k12 -k13
-k12 (k23+k12) -k23
-k13 -k23 (k23+k13)]; %<--macierz sztywnosci
B=[0 0 0
0 0 0
0 0 0]; %<-- macierz tlumienia
v0=[0 0 0]'; u0=[0 0 0]'; % warunki poczatkowe
y0=[v0; u0];
% czas=[0:50];
[w,k]=size(M);
Z=zeros(w,1);
Z1=zeros(w);
V=[Z1 M; M B];
W=[-M Z1; Z1 K];
[T,Y]=ode45(@licz,[0 5],y0,[],V,W) % rungy-kuta
hold on;
plot(T,Y(:,[2;5])); grid on; % wykres predkosci

function ddx=licz(t,y,V,W)
fw=[[0 0 -1]';[0 0 0]']; %wymuszenie
ddx=V\(fw-W*y);


to przyklad rozwiazania rownania roznikczkowego ukladu 3 masowego .... mysle ze na podstawie tego przyladu poradzisz sobie z Twoim problemem jak by co to pisz postaram sie pomoc
Edytowane przez mnichv dnia 06/26/2008 16:44
Autor RE: MES w MATLAB - trochę nietypowo
popey
Użytkownik

Postów: 3
Data rejestracji: 26.06.08
Dodane dnia 06/26/2008 23:40
Witam ponownie

Dziękuję za szybką odpowiedź - i to jaką.

Po wysłaniu mojego postu tutaj doszedłem do wniosku, że własnie może by zaprząc zwykłe rozwiązywanie równań różniczkowych w matlabie do tego problemu, poszperałem troszkę w internecie i wykombinowałem rozwiązanie także oparte na funkcji użytej przez Ciebie (jeżeli można) - ode45. Tyle, że mój kod był dwa razy dłuższy
Nie jestem dobrym programistą...

W każdym razie rozwiązuję teraz podane równanie, kreślę elipsy końcami wału i jestem szczęśliwy.

Teraz mam jeszcze jedno pytanie: ogólnie rozwiązanie miało być po to, by mieć model do " badań ". Teraz chciałbym zmieniać sztywności podpór, czyli wybrane współczynniki macierzy tak, aby zminimalizować drgania wału - czyli elipsy jakie kreślą końce wału.

Muszę obmyśleć w jaki sposób to zrobię - z jakich danych wyjściowych będę czerpał wskazówki, jakie metody optymalizacyjne użyć itp.

Z pewnością nie jest to już przedmiotem zainteresowania tego forum, zatem nie będę na siłę tutaj " smęcił "

Dziękuję raz jeszcze za pomoc!
Pozdrawiam
Autor RE: MES w MATLAB - trochę nietypowo
ADMIN
Administrator

Postów: 651
Miejscowość: Warszawa, PL
Data rejestracji: 18.10.07
Dodane dnia 06/27/2008 02:12
Problemy w MATLABie jak najbardziej nas interesują !!! To też się składa na CAE, może nie FEA, ale CAE ! Zresztą to forum (Wszystko o CAE, co nie pasuje do powyższych forum) jest poświęcone m.in. programom matematycznym


Pozdrawiam,
Admin

We live in an age when pizza gets to your home faster than police or an ambulance
Autor RE: MES w MATLAB - trochę nietypowo
popey
Użytkownik

Postów: 3
Data rejestracji: 26.06.08
Dodane dnia 06/27/2008 11:19
Witam

Trzymam Cię za słowo Pozdrawiam
Autor RE: MES w MATLAB - trochę nietypowo
wieczny_student
Użytkownik

Postów: 40
Miejscowość: Warszafka
Data rejestracji: 19.06.08
Dodane dnia 07/09/2008 01:16
Jak macie jakieś problemy z MATLAB to walcie do mnie Lubie rozwiązywać problemy w matlabie


Pozdrowienia,

Emalia:
wieczny_student-1976@o2.pl

student mechatroniki PW na 14-tym roku .
Autor RE: MES w MATLAB - trochę nietypowo
makman
Użytkownik

Postów: 3
Data rejestracji: 01.05.12
Dodane dnia 05/01/2012 14:26
Witam serdecznie,
Jak to się mówi odgrzebię kotleta i to dość starego jak widzę :P

Ale problem jest podobny więc nie ma sensu zakładać nowego wątku...

Za zadanie mam policzyć belkę jak na rysunku



próbowałem na kilka sposobów wklepać ode45 ale cały czas mi wali jakimiś errorami m.in jak ten :
Error using feval
Undefined function 'calkowanie' for input arguments of type 'double'.

Error in odearguments (line 88)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 114)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in belka2 (line 194)
wyn=ode45(@calkowanie,[0,5],y0)

Wiem, że jest długi weekend i najprawdopodobniej będę czekał długo na odpowiedź ale każda uwaga mile widziana ponieważ musze to zrobić jak najszybciej.

Docelowo muszę wyciągnąć drgania własne i przemieszczenia od czasu w poszczególnych węzłach aby móc wyrysować coś takiego ...


Edytowane przez makman dnia 05/02/2012 23:22
Autor RE: MES w MATLAB - trochę nietypowo
diadyn
Użytkownik

Postów: 34
Data rejestracji: 21.11.09
Dodane dnia 05/01/2012 19:02
Witam

po uruchomieniu Twojego programu pojawia się błąd w funkcji całkowanie, komunikat błędu:

Inner matrix dimensions must agree.
Error in ==> calkowanie at 188
yp=Vprim*W*y

w równaniu wektor y ma wymiar 24x1, a powinien mieć 48x1.
Porównując program mnichv z Twoim zauważyłem że zapomniałeś z wektorów y0 i y0w utworzyć jednego wektora
y0=[y0; y0w];

Pozdrawiam
Autor RE: MES w MATLAB - trochę nietypowo
makman
Użytkownik

Postów: 3
Data rejestracji: 01.05.12
Dodane dnia 05/01/2012 21:44
nie no wyniki wychodzą całkowicie absurdalne, coś jest nie tak w moim kodzie... a raczej w mojej metodologii rozwiązywania zadania.
Mam równanie Mq"+Kq=P
i teraz jak scałkować to równanie aby mieć współrzędne wszystkich węzłów ?

Szczerze powiedziawszy nie mam za dużo doświadczenia z mesami , matlabem i ogólnie zrozumienie tego zadania mnie przerasta a ksiązki z mes za wiele nie mówią o samym rozwiązywaniu problemów takich z jakimi siętu borykam czyli wyciągnięcie konkretnych wartości w przebiegu czasowym przy konkretnym wymuszeniu siłą /przesunięciem...

idę dalej z tym walczyć, jakby kogoś kopnęła wena jak można to zrobić to chętnie wysłucham bo powoli mi się pomysły kończą ;)
Edytowane przez makman dnia 05/02/2012 08:33
Autor RE: MES w MATLAB - trochę nietypowo
diadyn
Użytkownik

Postów: 34
Data rejestracji: 21.11.09
Dodane dnia 05/02/2012 09:57
Zanim zaczniesz rozwiązywać podane równanie poświęć trochę czasu na uporządkowaniem kodu
1.na początku definiujesz siłę p = 1000 potem do zmiennej p przypisujesz kolejno wartości od 0 do 20
2.w dwóch miejscach liczysz wartości własne [w,v]=eig(MP^-1*KP) i omega2 = eig(K,M), czy to konieczne ?
3.tak samo jest z uwzględnianiem warunków brzegowych, z macierzy Mw i Kw nie korzystasz więc możesz je usunąć

Na zdjęciach w węzłach masz po 2 stopnie swobody V i U, chyba powinno być V - przemieszczenie pionowe i jakiś kąt obrotu powiedzmy phi.
Wektorem P, przykładasz dwa momenty gnące, możesz to sprawdzić wyświetlając przemieszczenia węzłów uzyskanych z równania U=inv(K)*P.
Pamiętaj tylko że pionowe przemieszczenia węzłów są w nieparzystych wierszach, by je wyświetlić możesz użyć funkcji plot(U(1:2:end)),
a dla kątów plot(U(2:2:end))

A teraz równanie mq"+kq=P, możesz próbować dalej z funkcją ode45 lub użyć algorytmu Newmarka, kod poniżej.
Doklej go do swojego kodu, zaraz za poleceniem omega, wszystko poniżej usuń.

Najłatwiej jest założyć, że w chwili 0 układ znajduje się w spoczynku. Odpowiada to przypisaniu zerowych wartości wszystkim
przemieszczeniom - u0, prędkością - udot i przyspieszeniom - u2dot w układzie, oczywiście możesz nadać inne wartości.
Pobudzenie układu odbywa się przez przyłożenie sił w 11 stopniu swobody.

historię przemieszczeń węzłów znajdziesz w macierzy u, jeżeli chcesz obejrzeć stan układ w chwili t2 wywołujesz polecenie plot(u(1:2:end,2))



% -----------------------------------------------------------
[r,k]=size(K);

t=5;

% time step
dt=0.001;

% Przemieszczenie
u0=zeros(r,1);

% Prędkość
udot0=zeros(r,1);

% Przyspieszenie
u2dot0=zeros(r,1);

% Wymuszenie (Siła)
R0=zeros(r,0);

R0(12,1)=20/2;
R0(14,1)=20/2;

% parametry
beta=0.5;
alfa=0.25*(0.5+beta)^2;

% Stałe
a0=1/(alfa*dt^2);
a1=beta/(alfa*dt);
a2=1/(alfa*dt);
a3=1/(2*alfa)-1;
a4=beta/alfa-1;
a5=dt/2*(beta/alfa-2);
a6=dt*(1-beta);
a7=beta*dt;

% Efektywna macierz sztywności
K_e=K+a0*M;

u=horzcat(u0,zeros(r,length(dt:dt:t)-1));
udot=horzcat(udot0,zeros(r,length(dt:dt:t)-1));
u2dot=horzcat(u2dot0,zeros(r,length(dt:dt:t)-1));
R=horzcat(R0,zeros(r,length(dt:dt:t)-1));


for i=1:1:length(dt:dt:t)-1;

% effective load
R_e=R(:,i)+M*(a0*u(:,i)+a2*udot(:,i)+a3*u2dot(:,i));

u(:,i+1)=inv(K_e)*R_e;

u2dot(:,i+1)=a0*(u(:,i+1)-u(:,i))-a2*udot(:,i)-a3*u2dot(:,i);
udot(:,i+1)=udot(:,i)+a6*u2dot(:,i)+a7*u2dot(:,i+1);

end
clear a0 a1 a2 a3 a4 a5 a6 a7
clear K_e u0 udot0 u2dot0 R0
clear alfa beta



Autor RE: MES w MATLAB - trochę nietypowo
makman
Użytkownik

Postów: 3
Data rejestracji: 01.05.12
Dodane dnia 05/02/2012 23:25
Bardzo dziękuję za pomoc, bez metody newmarka raczej nie dałbym rady ;)

Poniżej cały kod, może komuś się przyda , rysowanie zrobione brute forcem - nie miałem już ochoty bawić się w pisanie pętli.
Pozdrawiam i życzę miłego wypoczynku przez resztę weekendu majowego ;)

l=4;
ro=7.85*10^3;
E=2.06*10^11;
%wprowadzanie wymuszen i war. początkowych
x=input('podaj wartość wychylenia środkowego bloczka od 0 do 1: ');
while x<0
disp(' !!! wychylenie musi być wartością dodatnią !!!');
while x>=1
disp(' !!! wychylenie musi być wartością mniejszą lub równą 1 !!!');
x=input('podaj wartość wychylenia środkowego bloczka od 0 do 1: ');
end
x=input('podaj wartość wychylenia środkowego bloczka od 0 do 1: ');
end

q=input('podaj wartość siły obciążającej środkowy bloczek (od -100 do 100): ');
while q>=100.1
disp(' !!! siła musi być wartością mniejszą lub równą 100 !!!');
q=input('podaj wartość siły obciążającej środkowy bloczek (od -100 do 100): ');
while q<=-100
disp(' !!! siła musi być wartością większą lub równą -100 !!!');
q=input('podaj wartość siły obciążającej środkowy bloczek (od -100 do 100): ');
end
end
%geometria el. skonczonego
L=l/5;

R=0.01;
R1=0.01;
A=3.14*R^2;
A1=3.14*R1^2;
L1=1.8;

J3=(3.14*R^4)/4;

J1=(3.14*R1^4)/4;
%boczne belki
MP=ro*L*[(13*A)/35+(6*J3)/(5*L^2) (11*A*L)/210+J3/(10*L) (9*A)/70-(6*J3)/(5*L^2) (-13*A*L)/420+J3/(10*L);
(11*A*L)/210+J3/(10*L) (A*L^2)/105+(2*J3)/15 (13*A*L)/420-J3/(10*L) (-A*L^2)/140-J3/30;
(9*A)/70-(6*J3)/(5*L^2) (13*A*L)/420-J3/(10*L) (13*A)/35+(6*J3)/(5*L^2) (-11*A*L)/210-J3/(10*L);
(-13*A*L)/420+J3/(10*L) (-A*L^2)/140-J3/30 (-11*A*L)/210-J3/(10*L) (A*L^2)/105+(2*J3)/15]

KP=(E/L)*[(12*J3)/(L^2) (6*J3)/L -(12*J3)/(L^2) (6*J3)/L;
(6*J3)/L 4*J3 -(6*J3)/L 2*J3;
-(12*J3)/(L^2) -(6*J3)/L (12*J3)/(L^2) -(6*J3)/L;
(6*J3)/L 2*J3 -(6*J3)/L 4*J3]
%środkowa belka
MP1=ro*L1*[(13*A1)/35+(6*J1)/(5*L1^2) (11*A1*L1)/210+J1/(10*L1) (9*A1)/70-(6*J1)/(5*L1^2) (-13*A1*L1)/420+J1/(10*L1);
(11*A1*L1)/210+J1/(10*L1) (A1*L1^2)/105+(2*J1)/15 (13*A1*L1)/420-J1/(10*L1) (-A1*L1^2)/140-J1/30;
(9*A1)/70-(6*J1)/(5*L1^2) (13*A1*L1)/420-J1/(10*L1) (13*A1)/35+(6*J1)/(5*L1^2) (-11*A1*L1)/210-J1/(10*L1);
(-13*A1*L1)/420+J1/(10*L1) (-A1*L1^2)/140-J1/30 (-11*A1*L1)/210-J1/(10*L1) (A1*L1^2)/105+(2*J1)/15]

KP1=(E/L1)*[(12*J1)/(L1^2) (6*J1)/L1 -(12*J1)/(L1^2) (6*J1)/L1;
(6*J1)/L1 4*J1 -(6*J1)/L1 2*J1;
-(12*J1)/(L1^2) -(6*J1)/L1 (12*J1)/(L1^2) -(6*J1)/L1;
(6*J1)/L1 2*J1 -(6*J1)/L1 4*J1]


M=zeros(24,24);
K=zeros(24,24);

p=0;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=2;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=4;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=6;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=8;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=10;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP1(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP1(iw,ik)+K(iw+p,ik+p);
end
end

p=12;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=14;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=16;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=18;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

p=20;
for iw=1:4
for ik=1:4
M(iw+p,ik+p)=MP(iw,ik)+M(iw+p,ik+p);
K(iw+p,ik+p)=KP(iw,ik)+K(iw+p,ik+p);
end
end

K(:,[1 2 23 24])=0;
K([1 2 23 24],:)=0;
K(1,1)=1;
K(2,2)=1;
K(23,23)=1;
K(24,24)=1;

M(:,[1 2 23 24])=0;
M([1 2 23 24],:)=0;
M(1,1)=1;
M(2,2)=1;
M(23,23)=1;
M(24,24)=1;

Q=zeros(1,24)';
Q(12,1)=q/2;
Q(14,1)=q/2;

U=inv(K)*Q;
omega2 = eig(K,M);
fmes = sqrt(omega2(2))/(2*pi);
omega = ( pi/(2*(L+1/2*L1)) ) * sqrt(E/ro);
f = omega/(2*pi) ;



% -----------------------------------------------------------
[r,k]=size(K);

t=5;

% time step
dt=0.1;

% Przemieszczenie
u0=zeros(r,1);
u0(11,1)=x;
u0(13,1)=x;

% Prędkość
udot0=zeros(r,1);

% Przyspieszenie
u2dot0=zeros(r,1);

% Wymuszenie (Siła)
R0=zeros(r,0);
R0(11,1)=q/2;
R0(13,1)=q/2;

% parametry
beta=0.5;
alfa=0.25*(0.5+beta)^2;

% Stałe
a0=1/(alfa*dt^2);
a1=beta/(alfa*dt);
a2=1/(alfa*dt);
a3=1/(2*alfa)-1;
a4=beta/alfa-1;
a5=dt/2*(beta/alfa-2);
a6=dt*(1-beta);
a7=beta*dt;

% Efektywna macierz sztywności
K_e=K+a0*M;

u=horzcat(u0,zeros(r,length(dt:dt:t)-1));
udot=horzcat(udot0,zeros(r,length(dt:dt:t)-1));
u2dot=horzcat(u2dot0,zeros(r,length(dt:dt:t)-1));
R=horzcat(R0,zeros(r,length(dt:dt:t)-1));


for i=1:1:length(dt:dt:t)-1;

% effective load
R_e=R(:,i)+M*(a0*u(:,i)+a2*udot(:,i)+a3*u2dot(:,i));

u(:,i+1)=inv(K_e)*R_e;

u2dot(:,i+1)=a0*(u(:,i+1)-u(:,i))-a2*udot(:,i)-a3*u2dot(:,i);
udot(:,i+1)=udot(:,i)+a6*u2dot(:,i)+a7*u2dot(:,i+1);

end
clear a0 a1 a2 a3 a4 a5 a6 a7
clear K_e u0 udot0 u2dot0 R0
clear alfa beta

u2=u(3,:)
rot2=u(4,:)
u3=u(5,:)
rot3=u(6,:)
u4=u(7,:)
rot4=u(8,:)
u5=u(9,:)
rot5=u(10,:)
u6=u(11,:)
rot6=u(12,:)
u7=u(13,:)
rot7=u(14,:)
u8=u(15,:)
rot8=u(16,:)
u9=u(17,:)
rot9=u(18,:)
u10=u(19,:)
rot10=u(20,:)
u11=u(21,:)
rot11=u(22,:)

jj_end=max(size(u));

for jj=1:1:jj_end


hold off
plot(0,0)
hold on

%el1
plot([0.2,1-(tan(rot2(jj))*0.5)],[2.8,2.8+u2(jj)])
plot([0.2,1+(tan(rot2(jj))*0.5)],[3.2,3.2+u2(jj)])
plot([1-(tan(rot2(jj))*0.5),1+(tan(rot2(jj))*0.5)],[2.8+u2(jj),3.2+u2(jj)])

%el2
plot([1-tan(rot2(jj))*0.5,1.8-(tan(rot3(jj))*0.5)],[2.8+u2(jj),2.8+u3(jj)])
plot([1+tan(rot2(jj))*0.5,1.8+(tan(rot3(jj))*0.5)],[3.2+u2(jj),3.2+u3(jj)])
plot([1.8-tan(rot3(jj))*0.5,1.8+tan(rot3(jj))*0.5],[2.8+u3(jj),3.2+u3(jj)])

%el3
plot([1.8-tan(rot3(jj)*0.5),2.6-(tan(rot4(jj))*0.5)],[2.8+u3(jj),2.8+u4(jj)])
plot([1.8+tan(rot3(jj)*0.5),2.6+(tan(rot4(jj))*0.5)],[3.2+u3(jj),3.2+u4(jj)])
plot([2.6-tan(rot4(jj)*0.5)*0.5,2.6+tan(rot4(jj))*0.5],[2.8+u4(jj),3.2+u4(jj)])

%el4
plot([2.6-tan(rot4(jj))*0.5,3.4-(tan(rot5(jj))*0.5)],[2.8+u4(jj),2.8+u5(jj)])
plot([2.6+tan(rot4(jj))*0.5,3.4+(tan(rot5(jj))*0.5)],[3.2+u4(jj),3.2+u5(jj)])
plot([3.4-tan(rot5(jj))*0.5,3.4+tan(rot5(jj))*0.5],[2.8+u5(jj),3.2+u5(jj)])

%el5
plot([3.4-tan(rot5(jj))*0.5,4.2],[2.8+u5(jj),2.8+u6(jj)])
plot([3.4+tan(rot5(jj))*0.5,4.2],[3.2+u5(jj),3.2+u6(jj)])

%el6 (uproszczony) brak rotacji ścian el 5 i 7 przyczepione na
%stałe do ścainek bocznych.
plot([4.2,6.2,6.2,4.2,4.2],[2+u6(jj),2+u6(jj),4+u7(jj),4+u7(jj),2+u6(jj)])

%el7
plot([7-tan(rot8(jj))*0.5,6.2],[2.8+u8(jj),2.8+u7(jj)])
plot([7+tan(rot8(jj))*0.5,6.2],[3.2+u8(jj),3.2+u7(jj)])
plot([7-tan(rot8(jj))*0.5,7+tan(rot8(jj))*0.5],[2.8+u8(jj),3.2+u8(jj)])

%el8
plot([7-tan(rot8(jj))*0.5,7.8-(tan(rot9(jj))*0.5)],[2.8+u8(jj),2.8+u9(jj)])
plot([7+tan(rot8(jj))*0.5,7.8+(tan(rot9(jj))*0.5)],[3.2+u8(jj),3.2+u9(jj)])
plot([7.8-tan(rot9(jj))*0.5,7.8+tan(rot9(jj))*0.5],[2.8+u9(jj),3.2+u9(jj)])

%el9
plot([7.8-tan(rot9(jj))*0.5,8.6-(tan(rot10(jj))*0.5)],[2.8+u9(jj),2.8+u10(jj)])
plot([7.8+tan(rot9(jj))*0.5,8.6+(tan(rot10(jj))*0.5)],[3.2+u9(jj),3.2+u10(jj)])
plot([8.6-tan(rot10(jj))*0.5,8.6+tan(rot10(jj))*0.5],[2.8+u10(jj),3.2+u10(jj)])

%el10
plot([8.6-tan(rot10(jj))*0.5,9.4-(tan(rot11(jj))*0.5)],[2.8+u10(jj),2.8+u11(jj)])
plot([8.6+tan(rot10(jj))*0.5,9.4+(tan(rot11(jj))*0.5)],[3.2+u10(jj),3.2+u11(jj)])
plot([9.4-tan(rot10(jj))*0.5,9.4+tan(rot10(jj))*0.5],[2.8+u11(jj),3.2+u11(jj)])

%el11
plot([9.4-tan(rot11(jj))*0.5,10.2],[2.8+u11(jj),2.8])
plot([9.4+tan(rot11(jj))*0.5,10.2],[3.2+u11(jj),3.2])

%sciana
plot(0,0,5,5, [0.2,0.2],[2.3,3.6])
plot(0,0,5,5, [10.2,10.2],[2.3,3.6])

%kreskowanie
for ii=9:12
plot([0.2 0.1], [4/15*ii 4/15*ii+0.1])
end
for ii=9:12
plot([10.3 10.2], [4/15*ii+0.1 4/15*ii])
end
pause(0.05)
end
disp('Siła wymuszenia=');
q
disp('Przesunięcie środkowego bloczka=');
x
disp('Częstość własna układu');
omega




Skocz do Forum:
Reklama
Najtańszy sklep z dobrymi oponami

Szukaj
Google
 
Logowanie
Nazwa Użytkownika

Hasło



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

Zapomniane hasło?
Wyślemy nowe, kliknij TUTAJ.
Ankieta dla Użytkowników
Jakiego używasz programu MES

Nastran (wszelkiej maści)

Ansys

Abaqus

Cosmos/M

Comsol

Algor

MARC

Inny

Musisz się zalogować, żeby móc głosować w tej Ankiecie.
Reklama
Shoutbox
Tylko zalogowani mogą dodawać posty w shoutboksie.

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

maq9633
10.11.2019
Witam, Zaczynam pisać pracę z płaskich den zbiorników ciśnieniowych i potrzebuje pomocy w programie Ansys workbench, proszę o kontakt

Vloiamwc
13.01.2019
Witam. Czy ktoś pomoże przygotować mnie do egzaminu z mes?

markos2101
30.01.2018
Dziękuję chociaż chodziło o znaczenie w MES - czyli zamiana obciążeń ciągłych na skupione

jasiu
29.01.2018
kongruentny: właściwy ; adekwatny ; zgodny ; poprawny ; godziwy ; dobry ; odpowiedni ; nadający się ; akuratny ; w sam raz ; należyty ; słuszny ; stosowny ; trafny ; przystający ;

markos2101
19.01.2018
Poszukuję definicji obciążenia kongruentnego - pomoże ktoś?

Adasiooo1
13.01.2018
witam szukam osób znających dobrze MES i abaqusa w zakresie żelbetu budownictwo

piotr
29.05.2017
Klachula nie jest do zadawania pytań od tego jest forum...

Geralt7586
28.05.2017
moźe ktoś mi pomoże ucze się femapa i mam problem jak użyję bar do rury to m=23kg a jak tube to 11,3 i wyrzuca mi blad SEKRRS przy tube a przy bar nie xd

kzurawski
08.04.2017
Jeśli chodzi o felgę to tu mamy jedną przeliczoną http://caelinux.org/wiki/
index.php/Contriblaws/C
ode_Aster/10_x_cases/plas
ticity


Archiwum
Wydarzenia
PonWtrCzwPiSoNie
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 31
Copyright Piotr© 2007- - 2019
Powered by Peter-Fusion © 2009 - Modified Aztec Theme