> | restart; |
> | with(linalg):with(plots):with(plottools): |
Warning,
the protected names norm and trace have been redefined and unprotected
Warning,
the name changecoords has been redefined
Warning,
the name arrow has been redefined
> | libname:="C:/",libname: |
> | with(mecapac3d): |
Paso 1. Definimos las coordenadas generalizadas del sistema en una lista que se denominará cg.
> | cg:=[alpha,theta]: |
> | rotdisc:=evalm(rota(d/R*alpha,1) &*rota(Pi/2,2)): |
> | d1:=[disco,[0,0,0],rotdisc,M,R]: |
> | m1:=[punto,0,R*cos(theta),R*sin(theta),m]: |
Paso 3. Definición de los elementos gráficos que definiran nuestro sistema de ejes.
> | ejeX:=[vector,[0,0,0],[10,0,0],red]: |
> | ejeY:=[vector,[0,0,0],[0,10,0],green]: |
> | ejeZ:=[vector,[0,0,0],[0,0,10],blue]: |
> | TO := [texto,[0,0,-1],"O"]: |
> | TX := [texto,[10,0,1],"X"]: |
> | TY := [texto,[0,10,1],"Y"]: |
> | TZ := [texto,[0,0,11],"Z"]: |
Paso 4. Definición de la variable sistema que agrupa en una lista todos los elementos anteriores.
> | cgsub:=[d*cos(alpha),d*sin(alpha),R]: |
> | rotsub:=rota(alpha,3): |
> | s1:=[subsistema2,cgsub,rotsub,[d1,m1]]: |
> | sistema:=[s1,ejeX,ejeY,ejeZ,TO,TX,TY,TZ]: |
Paso 5. Obtención de la energía cinética del sistema mediante fT asignándola a la variable T.
> | T:=fT(sistema); |
Paso 6. Obtención de la energía potencial del sistema mediante fV asignándola a la variable V.
> | V:=fV(sistema); |
> | L:=T-V: |
> | ecua:=ec_lag(); |
> | g:=9.8:M:=15:m:=5:R:=2.5:d:=3: |
> | res:=fint([0,1,0,0]): |
> | dibu3(2,50); |
> |