Руководство к решению задач по механике материалов и конструкций. Егодуров Г.С - 232 стр.

Const kl=50;
Var i,i1,i2,n,n1:integer;
         a,l,h,bb,b,b1,c1,c2,q,E,F,z1,z,R1,R2,Ix:real;
         T,V,S,K,M,Wx,ZZ:array[1..k1] of real;
         Otlet:char;
Begin
Clrscr;
Writeln(‘Размер балки’);
Writeln(‘l(mm),h(mm),bb(mm),b(mm)’);
Readln(l,h,bb,b);
Writeln(‘Нагрузка q(H/mm),F(H), Модуль упругости Е(MПа’);
Readln(q,F,E);
Writeln(‘Число шагов интегрирования n’);
Readln(n);
Writeln(‘Результаты расчета’);
Writeln;
Writeln(‘z mm MH*m S MПа                  V mm T рад’);
Writeln;
a:=4*l;
z1:=a/n; {Шаг интегрирования }
n1:=n+1; {Количество узлов}
R1:=3*F/2+3*q*1/4; R2:=9*q*1/4-F/2;
For i:=1 to n1 do
Beqin
z:=z1*(i-1);
ZZ[I]:=z;
{Уравнение изгибающих моментов по участкам}
If (z>=0) and (z<=l) then M[i]:=-F*z;
If (z>l) and (z<=3*l) then M[i]:=-F*z+R1*(z-l)-q*(z-l)*(z-l)/2;
If (z>3*l) and (z<=4*l) then M[i]:=-F*z+R1*(z-l)-q*(z-l)*(z-l)/2+R2*(z-3*l);
{Геометрические характеристики}
If (z>=0) and (z<=l) then b1:=bb;
If (z>1) and (z<=4*l) then b1:=bb-(bb-b)*(z-l)/(3*l);
Ix:=b1*h*h*h/12;
Wx[I]:=2*Ix/h;
S[i]:=M[i]/Wx[i]; {Нормальные напряжения}
K[i]:=M[i]/(E*Ix); {Кривизна}
end;
T[1]:=0;
For i:=1 to n do T[i+1]:=T[i]+z1*(K[i]+K[i+1])/2;
V[1]:=0;
For i:=1 to n do V[i+1]:=V[i]+z1*(T[i]+T[i+1])/2;
c1:=(V[6]-V[16])/(2*l);
c2:=V[6]+c1*l;
For I:=1 to n1 do
Beqin
T[i]:=T[i]+c1; {Угол поворота}
V[i]:=-V[i]+c1*ZZ[i]+c2; {Прогиб}
Writerln(ZZ[i]:6:1,M[i]/1E3:9,S[i]:9:2,V[i]:8:2,T[i]:11:5);
End;
Readln; End.




                                               232