|

楼主 |
发表于 2006-12-17 22:36
|
显示全部楼层
syms x y a b w0 k n E1 E A H h v D q w f1 f2 f3 f4 f5 f6 f7
%D=E*h^3/12/(1-v^2)
w=w0/(pi/2-1)*(cos((x^2/a^2+y^2/b^2)*pi/2)+pi/2*(x^2/a^2+y^2/b^2-1))
f1=D*(diff(w,x,4)+2*diff(diff(w,x,2),y,2)+diff(w,y,4))
f1=simple(f1)
f2=k*w-q
f3=f1+f2
f4=f3*w/w0*(pi/2-1)
f5=int(int(f4,x,-a*sqrt(1-y^2/b^2),a*sqrt(1-y^2/b^2)),y,-b,b)
%f5=4*int(int(f4,x,0,a*sqrt(1-y^2/b^2)),y,0,b)
%f5=4*int(int(f3,x,0,a*sqrt(1-y^2/b^2)),y,0,b)
%simple(f5)
%pretty(ans)
%f6=solve(f3,w0) |
|