Hi
I would like to include a periodic external force into a multiple mass
damper-system. My first thought was to include it in the forcelist (using
KanesMethod)
which would then become something like in the code below. (where b is the
frequency of the periodic force)
This however, returns an error whenever i want to integrate the dynamical
equations.
Does anyone know how to include the external force f(t) = sin(b*t), so i
can integrate the dynamical equations?
Thanks
q1,u1,q2,u2 = dynamicsymbols('q1 u1 q2 u2')
q1d, q2d = dynamicsymbols('q1 q2',1)
m1,c1,k1, m2,c2,k2, g,t,b = symbols('m1 c1 k1 m2 c2 k2 g t b')
N = ReferenceFrame('N')
O=Point('O')
O.set_vel(N,0)
P = Point('P')
P.set_vel(N,-u1*N.y)
Q=Point('Q')
Q.set_vel(N,-(u2+u1)*N.y)
kd=[q1d-u1,q2d-u2]
FL =
[(P,(2*k1*q1-k1*q2+c1*u1)*N.y),(Q,((k1+k2)*q2-k1*q1+c2*u2+sin(b*t))*N.y)]
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sympy.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/b1b1fd3f-196f-4151-a3c1-c94e82ef18c4%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
