MES,CAD,CAE,FEA
News Search:
StartGaleriaForumDownloadArtykułyNewsyLinkiZasoby wiedzyKontaktZasady Stosujemy pliki cookies więcej...03/28/2024 20:13:36
Reklama
Nawigacja
Galeria
Forum
Download
Artykuły
FAQ
Linki
Zasoby wiedzy
Kontakt
Szukaj
Losowa Fotka
Tradycyjny zawias
Tradycyjny zawias
Fajne gify
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:

piotr12 tygodni
ndv33 tygodni
damian1410045 tygodni
BE-FEA47 tygodni
adrian2402451 tygodni

Gości Online: 23

Twoje IP to: 54.160.243.44

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
model 3

model 3

Model Felgi
Felga model 1 mesh

Felga model 1 mesh

Model Felgi
Wazon

Wazon

FEM-art, ukryte piękno obliczeń
Ssanie sprężanie... kupa pracy

Ssanie sprężanie... kupa pracy

Fajne gify
Naprężenia redukowane w kompozycie

Naprężenia redukowane w kompozycie

Ząb
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
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 | 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.

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
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
Copyright Piotr© 2007- - 2024
Powered by Peter-Fusion © 2009 - Modified Aztec Theme