maple_matt_2

0 Reputation

One Badge

12 years, 341 days

MaplePrimes Activity


These are replies submitted by maple_matt_2

Thank you!  I am still confused as to why this happened.  I literally copied and pasted the elements of the solution, as Maple displays them, ((1) in blue) from EQN when forming this expression.  One would think that Maple would know its own code!

So I guess the rule is, whether copied and pasted or not, always go through and put a multiplication in anything whatever that is supposed to be a multiplication.  This is a little bit of a PITA, but now I know. 

 

Thank you!  I am still confused as to why this happened.  I literally copied and pasted the elements of the solution, as Maple displays them, ((1) in blue) from EQN when forming this expression.  One would think that Maple would know its own code!

So I guess the rule is, whether copied and pasted or not, always go through and put a multiplication in anything whatever that is supposed to be a multiplication.  This is a little bit of a PITA, but now I know. 

 

I don't see my worksheet showing up, so I will try again.

 

Kinematics_2.mw

 

  restart

NULL

E[before] := E[Typesetting:-delayDotProduct(pot, collar)*initial]+E[Typesetting:-delayDotProduct(pot, beam)*initial]:

NULL

NULL

NULL

E[Typesetting:-delayDotProduct(pot, collar)*initial] := M[c]*g*L*sin(beta):

E[Typesetting:-delayDotProduct(pot, beam)*initial] := (1/2)*M[b]*g*L*sin(beta):

 

E[Typesetting:-delayDotProduct(kin, wheel)] := (1/2)*Iota[wheel]*omega[wheel]^2+(1/2)*M[w]*V[Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(c, m), wheel)]^2:

NULL

NULL

E[Typesetting:-delayDotProduct(kin, collar)] := (1/2)*M[c]*V[c]^2:

NULL

E[Typesetting:-delayDotProduct(kin, beam)] := (1/2)*Iota[Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(c, m), beam)]*omega[beam]^2+(1/2)*M[b]*V[Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(c, m), beam)]^2:

NULL

E[after] := E[Typesetting:-delayDotProduct(kin, wheel)]+E[Typesetting:-delayDotProduct(kin, collar)]+E[Typesetting:-delayDotProduct(pot, collar)]+E[Typesetting:-delayDotProduct(kin, beam)]+E[Typesetting:-delayDotProduct(pot, beam)]:

NULL

NULL

NULL

Now with all of the components in terms of either θ or Vc we can solve.

 

 

EQN := E[before] = E[after];

M[c]*g*L*sin(beta)+(1/2)*M[b]*g*L*sin(beta) = (3/4)*M[w]*V[c]^2/tan(theta)^2+(1/2)*M[c]*V[c]^2+M[c]*g*L*sin(theta)+(1/6)*M[b]*V[c]^2*cos(theta)^2+(1/8)*M[b]*V[c]^2*(1+1/tan(theta)^2)+(1/2)*M[b]*g*L*sin(theta)

(1)

 

 

S := solve(EQN, V[c])

(2)

NULL

M[w] := 21:

``

(3)

evalf[20](S)

-2.5935955955772686654, 2.5935955955772686654

(4)

evalf[20](sqrt((M[c]*g*L*sin(beta)+(1/2)*M[b]*g*L*sin(beta)-M[c]*g*L*sin(theta)-(1/2)*M[b]*g*L*sin(theta))/((3/4)*M[w]/tan(theta)^2+(1/2)*M[c]+(1/6)*M[b]*cos(theta)^2+(1/8)*M[b](1+1/tan(theta)^2))))

2.6692045646476235434

(5)

NULL

NULL

NULL

``


Download Kinematics_2.mw

Page 1 of 1