Расчет составной конструкции
Пример с.54 (Кирсанов М.Н. Решебник.Теоретическая механика)
Ввод координат опор А, Е и шарнира С
restart;
| > | al:=Pi/3: Угол альфа |
| > | xa:=0: xc:=4*cos(al)+6: xe:=xc+4*cos(al)+2: |
| > | ya:=0: yc:=-4*sin(al): ye:=yc-4*sin(al): |
Момент
| > | M:=100: |
Действующие силы. Порядок записи - [Величина, x, y , угол(от полож.напр. оси х против час.стр.)]
| > | FL:=[[20,2*cos(al),-2*sin(al),0], P [10,xc-2,yc,Pi/180*50], Q [16,2*cos(al),-2*sin(al),-Pi/2], Gab [24,xc-3,yc,-Pi/2]]: Gbc |
Силы, приложенные к правой части (СЕ)
| > | FR:=[[16,xc+1,(yc-ye)/2,-Pi/2], Gcd [8,xe-1,ye,-Pi/2]]: Gde |
Процедура вычисления суммы моментов относительно точки с кординатами (x,y)
| > | Mom:=proc(x,y,F) evalf(add(F[i][1]*(cos(F[i][4])*(y-F[i][3]) -sin(F[i][4])*(x-F[i][2])),i=1..nops(F))): end: |
Процедуры вычисления суммы проекций на оси
| > | Prx:=proc(F) evalf(add(F[i][1]*cos(F[i][4]),i=1..nops(F))):end: |
| > | Pry:=proc(F) evalf(add(F[i][1]*sin(F[i][4]),i=1..nops(F))):end: |
Уравнение равновесия.Сумма моментов относительно опоры
| > | eq1:=Xa*(ye-ya)-Ya*(xe-xa)+Mom(xe,ye,FR)+Mom(xe,ye,FL)+M; |
Уравнение равновесия.Сумма моментов относительно шарнира
| > | eq2:=Xa*(yc-ya)-Ya*(xc-xa)+Mom(xc,yc,FL); |
| > | sol:=solve({eq1,eq2},{Xa,Ya}); assign(sol): |
Вычисление реакции опоры из уравнений проекций
| > | Xe:=-(Xa+Prx(FL)+Prx(FR)); |
| > | Ye:=-(Ya+Pry(FL)+Pry(FR)); |
Проверка. Сумма моментов относительно опоры А
| > | evalf(Mom(xa,ya,FR)+Mom(xa,ya,FL)+M+Xe*(ya-ye)-Ye*(xa-xe)); |
| > |