Response of single DOF spring-mass-dashpot model to half-sine pulse
Response of single DOF spring-mass-dashpot model to half-sine pulse
(OP)
Folks-
The problem is to find the solution to mx''+cx'+kx=f where f is a half-sine pulse. The initial condition are x(0)=a positive constant, say B and x'(0)=0. Years ago I attempted (over many hours!) a closed form solution to this problem because it is a practical problem in electronics packaging design (e.g. MIL-STD-810 Functional Shock). The forcing function is easy to construct with unit step functions and a sine function. I tried several methods concluding with attempt with Laplace transorms. Each solution got so long and convoluted and lengthy that I tried a new approach. Now I'm asking you folks. BTW I've no luck finding a solution in the literature.
The problem is to find the solution to mx''+cx'+kx=f where f is a half-sine pulse. The initial condition are x(0)=a positive constant, say B and x'(0)=0. Years ago I attempted (over many hours!) a closed form solution to this problem because it is a practical problem in electronics packaging design (e.g. MIL-STD-810 Functional Shock). The forcing function is easy to construct with unit step functions and a sine function. I tried several methods concluding with attempt with Laplace transorms. Each solution got so long and convoluted and lengthy that I tried a new approach. Now I'm asking you folks. BTW I've no luck finding a solution in the literature.
Tunalover





RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
F(t)
=(t<0)* (B*k)
+( t>=0 * t=<.5/f)*(sin(2*pi*f*t))
+(t>.5/f)*????
is that right? what is the forcing function after the end of the pulse?
I think Laplace is the right approach, but I bet a crafty bugger could use a superimposition of two systems with a time delay to do it classically.
Cheers
Greg Locock
Please see FAQ731-376 for tips on how to make the best use of Eng-Tips.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
That's a bit hard to understand physically.
Cheers
Greg Locock
Please see FAQ731-376 for tips on how to make the best use of Eng-Tips.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Cyril M. Harris, Shock and Vibration Handbook, McGraw-Hill, Chapter 8. My copy is the 3rd Edition, 1988. Amazon says the 5th ed. is now out; you should be able to find a copy in a good uni. library.
While the book doesn't give an explicit closed-form solution, it does present a plot of the response function versus reduced frequency. Also, this chapter provides two references for that problem (good luck finding them, you'll have to probably dig in the university library archives):
Criner, H. E., G. D. McCann, and C. E. Warren: Journal of Applied Mechanics, 12:135 (1945).
and
Mindlin, R. D., F. W. Stubner, and H. L. Cooper: Proceedings of the Society of Experimental Stress Analysis, 5:2, 69 (1948)
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Regards
Dave
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
He solves this problem for initial conditions of x(0) = 0 and x'(0) = 0. You should be able to use the same methodology to solve for x(0) = constant.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
The waveform replicates a drop shock. It's really not that unreasonable. The duration is typically 10-12ms.
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Cheers
Greg Locock
Please see FAQ731-376 for tips on how to make the best use of Eng-Tips.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Your forcing function can easily be represented in laplace transforms (produce of sin and window function which can be created as combination of step functions or "heaviside" functions).
When you transform the differential equation into Laplace domain, it becomes an algebraic equation. Solve for X(s) and then inverse transform to find x(t).
Attached is solution for the complete problem including initial conditions using Maple to help with the algebra. Without Maple it would have been ugly (especially the partial fraction expansion to find the inverse laplace transform), but do-able (if you had several hours).
ht
If for some reason that doens't work, go here:
http://home.houston.rr.com/electricpete/link3.htm
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
I'm looking for the closed form solution but thanks anyway.
electricpete-
Thanks for the input! As I said earlier, I tried the Laplace transform solution but that led to excessive algebra. I looked at your solution though and it looks like you managed it. What software generated the formulae?
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
http://www.maplesoft.com/
If you qualify for student version, it's a helluva lot cheaper
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
I have no problem with partial fraction expansions. What I have to do is find a quiet place for two hours to try to come up with your result. I'm one of those geeks who has to do the math in order to believe it!
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
If the duration of the forcing pulse is less than 1/3 of the period of the natural frequency of the single DOF system, then the response will behave as if the pulse is 1/3 duration. IE if the pulse duration gets shorter, then the systems just behaves the same.
As I write this, I have a little gremlin going off that says it might be less than 1/10 of the natural period.
Does this help?
Cheers
Solider
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
I also learned a rule of thumb that, for manually portable (you don't need a truck, forklift, or dolly to move it) electronic equipment, if the pulse duration is greater than 13ms then the shock pulse can be fairly replaced with a static load on the CG with magnitude AW where A is the amplitude of the pulse (in g's) and W is the weight of the assembly. Has anyone else heard this one?
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Here is the derivation for a single-degree-of-freedom system subjected to a half-sine pulse base excitation.
http://www.vibrationdata.com/sbase.pdf
Tom Irvine
http://www.vibrationdata.com
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Tunalover
RE: Response of single DOF spring-mass-dashpot model to half-sine pulse
Star for you Tom, amazing.
Cheers
Greg Locock
Please see FAQ731-376 for tips on how to make the best use of Eng-Tips.