I have a body made of cylinder with sphere at one base. It's constrained to move in z-direction though prismatic joint. Spherical base of the body is in contact with a plane as shown in the figure.
A time varying external force F(t) pushes the body towards the plane. F(t) is a perodic in nature and is described as:
where Tp is time period of F(t)
Ideally it's expected that normal reaction should eventually converge to applied external force. This holds true when Tp is greater than 0.12 sec. Figure below shows the case of 0.12 sec:
Pressure force is essentially applied external force. When Tp gets smaller, normal reaction is not able to rise & fall in given amount of time. Figure below shows the case of Tp = 0.01 sec:
Values of contact parameters are as follow:
- Mass: 0.031 lbf (0.014 kg)
- Contact stiffness: 5.7e06 lbf/in (1e09 N/m)
- Contact damping: 5.7lbf/(in/s) (1e03 N/m/s)
- Penetration exponent: 2.2 (standard for steel to steel contact)
To check whether contact parameters are suitable or not, I solved the same problem in 1D using simulink. There,the normal reaction was converging to applied external force in steady state even for low value of Tp. I'm not sure if this is a solver issue or something else. Currently I'm using daessc solver.