>

restart;

Размеры

>	c1:=12:c2:=20:

>	h1:=3:h2:=5: h3:=7.:

Бесконечность

>	infin:=1e20:

Координаты центров тяжести

>	y1:= h1/2; y2:=h1+h2+h3/2; y3:=(h1+h2+h3)/2;

Площади

>	F1:=-h1*c1;

>	F2:=-h3*c2;

>	F3:=(h1+h2+h3)*c2;

Общая площадь

>	F:=add(F||i,i=1..3);

Координата ц.т.

>	yc:=add(F||i*y||i,i=1..3)/F;

Моменты инерции

>	Jy1:=-h1*c1^3/12:

>	Jy2:=-h2*c1^3/12:

>	Jy3:=(h1+h2+h3)*c2^3/12:

>	Jy:=add(Jy||i,i=1..3);

>	Jx1:=(-h1^3*c1/12+F1*(y1-yc)^2);

>	Jx2:=(-h3^3*c1/12+F2*(y2-yc)^2);

>	Jx3:=(h1+h2+h3)^3*c2/12+F3*(y3-yc)^2:

>	Jx:=add(Jx||i,i=1..3);

Радиусы инерции в квадрате

>	ix2:=Jx/F; iy2:=Jy/F;

>

Уравнения контура фигуры в отрезках

>	b[1]:=h1+h2+h3-yc:  a[1]:=infin:

>	b[2]:=infin: a[2]:=c2/2:

>	b[3]:=-yc:   a[3]:=infin:

>	b[4]:=infin: a[4]:=-c2/2:

Координаты угл точек ядра сечения

>	for i to 4 do

>	x0[i]:=-iy2/a[i]; y0[i]:=-ix2/b[i];

>	printf("%8.3f,%8.3f\n",x0[i],y0[i]);

>

od:

  -0.000,  -0.600

  -7.135,  -0.000

  -0.000,   1.305

   7.135,  -0.000

>