No to urcite, chlastani nebude, v patek jsem sice byl, ale jenom slusne, protoze pres vikend sem se musel ucit. Uz mam za sebou dve zkousky, tu posledni jsem psal dneska odpoledne. Myslim ze obe dopadli docela dobre, ale krutou pravdu se teprve dozvim. Treba kruta nebude - sebevedomej odhat oboje kolem 90, pesimistickey 60-70 ... tak uvidime. Zjistil jsem ze zkousky nejsou tak hrozny, ze me nejvic pijou krev ty ukoly.
Kuprikladu dneska - dynamika, musim v matlabu (kterej se musime naucit uplne samotny) napsat funkci na odezvu elastoplastickeho systemu na zemetreseni. Uz jsem nad tim stravil hodiny, linearni systemy jsou primitivni, ale jakmile me to zacne jit nonlinearne aka elastoplasticky tak jsem v cajzlu. Vim jak to mam napsat teoreticky, ale ten zasranej matlab me proste nejde ;-). Jo jo, sorry, ale nekde se ventilovat musim.
skriptik na metodu "central difference" co me funguje jenom na pul, ma to vytahnou acceleraci ze zemetresni, "centralni differenci" na deformaci v case i, a pak non-linearni funkci preji na i+1, a vyflusnout deformace, casove intervaly non-linearniho chovani a konecne vztah sila-deformace
datat=data';
dt=0.02;
m=10;
k=157.7536*m;
t(1)=0;
u(1)=0;
pi(1)=0;
pi(2:length(data(:))+1)=-1000*datat(:);
a=m/(dt)^2;
b=(2*a);
khat=a;
phat(1)=pi(1)-fs(1)+b*u(1);
fy=.171*386
for i=1:750
u(i+1)=phat(i)/khat;
t(i+1)=i*0.02;
fs(i+1)=fs(i)+k*(u(i+1)-u(i));%problem with fs
if fs(i+1) > fy
fs(i+1) = fy
end
if fs(i+1) < -fy
fs(i+1)= -fy
end
phat(i+1)=pi(i+1)-a*u(i)-fs(i+1)+b*u(i+1);%since fs is wrong ....
end
plot(t,u)
display(max(u))
xlabel('T (s)')
ylabel('u (in)')
title('Displacement Response')
grid on
grid minor
no jo, machrovat se musi, prece si to vsechno nenecham jenom pro sebe ;-), tak, a je to venku, a uz me ani tak nesere, ale do zitra ho uz asi nespravim, ale verejne ventilovani na blogu je ocividne pomerne zabavne ;-), zejtra az me bude pit krev geotechnika tak mozna napisu neco z tamtoho soudku, tam uz taky potrebujem matlab :-(
Neco optimistictejsiho - mam se dobre, posledni dobou se preziram, protoze se ted rad varim, k snidani pravidelne vajicka se slaninou, k veceri bud grill a nebo neco na panvi, pripadne cokoliv uvari alex. Pocasi se nam taky dari, a v berkeley se me libi cim dal tim vic.
Na kole jezdim pravidelne, uz ho mam konecne dany poradne do kupy, brdzi, prehazuje, nevzrze, a jezdi pomerne rychle. O vikendu se asi nekam pojedu podivat, protoze dalsi zkouska je az pristi stredu. Teda jestli si nekdo nevpomene ze budem muset prijit na prednasky i o vikendu, s cim uz teda strasi minimalne dva profesori.
Tesim se dom, uz mam skoro letenku!
Přihlásit se k odběru:
Komentáře k příspěvku (Atom)
1 komentář:
Ja nerozumim ani tomu, co pises cesky, natoz tem matematickym velbludam!
Okomentovat