Skip to main content
1-Visitor
January 28, 2011
Question

mathcad: help with Numerical Soln to Differential Eqn

  • January 28, 2011
  • 12 replies
  • 6554 views

Please help to implement following differential equation numerical solution in mathcad

Numerical solution for single degree of freedom elastic system.
Spring stiffness = k = 2000 lbf/ft
Weight = W = 64 lbf

Force function:
F(t) = 500t+50 if t <=0.1
=-1250t+225 if 0.1 <=t<= 0.14
= 50 if t>=0.14
Acceleration function:
A(t) = ½ F(t) – 1000 *Y

Say time step = dt = 0.02
Initial Condition:
Displacement:
Y0 = 0 that is displacement is zero when t = 0
Y1 = ½ * A0 *dt^2 that is first time step which t = 0.02 this gives displacement but in terms of Acceleration at t =
Acceleration:
A0 = 0 when t = 0

Solution:
Y (t+dt) = 2 * Y (t) – Y(t-dt) + A(t) * dt^2
For t – 0 to 0.4 with dt = 0.02

Return:
Y, A

12 replies

1-Visitor
January 28, 2011

Please provide a worksheet.


Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

1-Visitor
January 28, 2011

attached is the mathcad sheet. This is all i could do

19-Tanzanite
January 28, 2011

Please post the worksheet in the other thread you started, where everyone has been commenting.

This is why we don't like cross posts.

1-Visitor
January 30, 2011

Took a bit of work, but:

1-Visitor
January 31, 2011

Thanks for your help.