Genspeeds2

# Genspeeds2 -...

This preview shows pages 1–2. Sign up to view the full content.

Sheet1 Page 1 % Blocks problem from homework problem 4.2 % used to demonstrate generalized speeds overwrite on % Set up problem newtonian a frames b,c bodies d points o mass d=md inertia d,i1,i2,i3,0,0,0 variables q{6}',u{6}' % Note that q is for generalized coordinates % and u is for generalized speeds % Rotation matrices simprot(a,b,3,q1) simprot(b,c,1,q2) simprot(c,d,3,q3) a_d=a_b*b_c*c_d % Kinematic relationships and rotational kinematics for second approach w_d_a>=u1*d1>+u2*d2>+u3*d3> w_d_a_alt>=q1'*a3>+q2'*b1>+q3'*c3> express(w_d_a_alt>,d) diff>=w_d_a>-w_d_a_alt> eqn[1]=dot(diff>,d1>) eqn[2]=dot(diff>,d2>) eqn[3]=dot(diff>,d3>) solve(eqn,q1',q2',q3') expand(q1') expand(q2') expand(q3') q4'=u4 q5'=u5 q6'=u6 alf_d_a>=dt(w_d_a>,d) % Via million dollar formula % Translational kinematics for second approach

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: p_o_do&gt;=q4*a1&gt;+q5*a2&gt;+q6*a3&gt; v_do_a&gt;=dt(p_o_do&gt;,a) a_do_a&gt;=dt(v_do_a&gt;,a) % Dynamics equations for second approach zero=fr()+frstar() expand(zero) % Note that to solve these equations via numerical integration, % you need to start with initial conditions for the following % variables: q1-&gt;q6 and u1-&gt;u6. Then two sets of first order % (rather than one set of second order) differential equations % must be solved via numerical integration: q' = f(q,u) and Sheet1 Page 2 % u' = f(u). The dependence of the dynamical (i.e., u') % differential equations on q has disappeared, and instead % the kinematical (i.e., q') differential equations now % depend on q, but in a much simpler form....
View Full Document

## Genspeeds2 -...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online