![{v2 = -2*g*(m[A]+m[B]+m[C])/(6*m[A]+3*m[B]+2*m[C]), v1 = -2*g*(3*m[A]-m[C])/(6*m[A]+3*m[B]+2*m[C])}](images/P221.gif)
Программа 22. Задача 41
| > | restart; |
| > | JC:=m[C]*Rc^2/2: # Момент инерции |
| > | JB:=m[B]*Rb^2/2: # Момент инерции |
| > | wB:=v1/Rb: # Угловая скорость |
| > | wC:=(v1+v2)/Rc: # Угловая скорость |
Кинетическая энергия
| > | T:=m[A]*v1^2/2+JB*wB^2/2+JC*wC^2/2+m[C]*v2^2/2: |
| > | N:=-m[A]*g*v1-m[C]*g*v2: # Мощность |
| > | Q1:=diff(N,v1): #Обобщенная сила 1 |
| > | Q2:=diff(N,v2): #Обобщенная сила 2 |
Уравнения Лагранжа
| > | u1:=collect(diff(T,v1),v1)=Q1: |
| > | u2:=collect(diff(T,v2),v1)=Q2: |
| > | solve({u1,u2},{v1,v2}); |
| > |