Программа 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;                # Время счета

gamma1 := .7462091399

beta := .2808129515

alpha := .2822393275

4.999

>    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);

>