Программа 37
> | restart;#Кручение. вычисление поправочных коэфф. |
> | N:=100: |
> | t1:=time(): # Засекаем время счета |
> | k:=2*n-1: |
Функция напряжений
> | F:=evalf(8*a^2/Pi^3*sum(k^(-3)*(-1)^(n-1)*(1-cosh(k*y*Pi/a)/cosh(k*b*Pi/2/a))*cos(k*x*Pi/a),n=1..N)): |
> | b:=4*a: # Соотношение сторон |
> | tau_zy:=-G*teta*diff(F,x):# Касательные напряжения |
> | tau_zx:=G*teta*diff(F,y): |
> | tau_max:=subs(y=0,x=a/2,tau_zy): |
> | tau_min:=subs(y=-b/2,x=0,tau_zx): |
> | gamma1:=tau_min/tau_max; |
Крутящий момент
> | M:=2*G*teta*int(int(F,x=-a/2..a/2),y=-b/2..b/2): |
> | beta:=M/(G*teta*a^3*b); |
> | alpha:=M/(tau_max*a^2*b); |
> | time()-t1; # Время счета |
> | a:=0.1:teta:=0.1:G:=80: |
> | plot3d(tau_zy,x=-a/2..a/2,y=-b/2..b/2,axes=BOX,lightmodel=light2,orientation=[75,60],shading=NONE); |
> |