0
Research Papers

# A Contact Force Model With Nonlinear Compliance and Residual Indentation1

[+] Author and Article Information
Xiaogang Xiong

Ryo Kikuuwe

e-mail: kikuuwe@mech.kyushu-u.ac.jp

Motoji Yamamoto

e-mail: yama@mech.kyushu-u.ac.jp
Department of Mechanical Engineering,
Kyushu University,
Fukuoka, 819-0395, Japan

This paper is an extended version of the authors’ conference paper [26]. Simulation results regarding the parameter $γ$ are newly included here. More detailed results and discussions are also presented.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received August 4, 2012; final manuscript received March 22, 2013; accepted manuscript posted May 7, 2013; published online September 16, 2013. Assoc. Editor: Wei-Chau Xie.

J. Appl. Mech 81(2), 021003 (Sep 16, 2013) (8 pages) Paper No: JAM-12-1370; doi: 10.1115/1.4024403 History: Received August 04, 2012; Revised March 22, 2013; Accepted May 07, 2013

## Abstract

The contact between two bodies is a complicated phenomenon in which the force and the relative position have nonlinear relations. Empirical results in the literature show that, in some mechanical systems such as biological tissues, the relation between the contact force and the indentation is characterized by the following three features: (i) continuity of the force at the time of collision, (ii) a Hertz-like nonlinear force-indentation curve, and (iii) nonzero indentation at the time of loss of contact force. The conventional Hunt–Crossley (HC) model does not capture the feature (iii) as the model makes the contact force and the indentation reach zero simultaneously. This paper proposes a compliant contact model based on a differential-algebraic equation that satisfies all three features. The behaviors of the model and the effect of the parameters in the model are investigated through numerical simulations.

<>
Your Session has timed out. Please sign back in to continue.

## References

Yamamoto, T., Vagvolgyi, B., Balaji, K., Whitcomb, L. L., and Okamura, A. M., 2009, “Tissue Property Estimation and Graphical Display for Teleoperated Robot-Assisted Surgery,” Proceedings of IEEE International Conference on Robotics and Automation (ICRA'09), Kobe, Japan, May 12–17, pp. 4239–4245.
Yamamoto, T., 2011, “Applying Tissue Models in Teleoperated Robot-Assisted Surgery,” Ph.D. thesis, Johns Hopkins University, Baltimore, MD.
Mu, X., and Wu, Q., 2006, “On Impact Dynamics and Contact Events for Biped Robots Via Impact Effects,” IEEE Trans. Syst., Man, Cybern., Part B: Cybern., 36(6), pp. 1364–1372.
Flores, P., Ambrósio, J., Claro, J. C. P., and Lankarani, H. M., 2006, “Influence of the Contact-Impact Force Model on the Dynamic Response of Multi-Body Systems,” Proc. Inst. Mech. Eng. Part K: J. Multi-Body Dyn.220(1), pp. 21–34.
Cross, R., 1999, “The Bounce of a Ball,” Am. J. Phys., 67(3), pp. 222–227.
Cross, R., 2011, Physics of Baseball and Softball, Springer, New York.
Ismail, K. A., and Stronge, W. J., 2012, “Viscoplastic Analysis for Direct Impact of Sports Balls,” Int. J. Non-Linear Mech., 47(2), pp. 16–21.
Hanley, K., Collins, F., Cronin, K., Byrne, E., Moran, K., and Brabazon, D., 2012, “Simulation of the Impact Response of a Sliotar Core With Linear and Non-Linear Contact Models,” Int. J. Impact Eng.50, pp.113–122.
Schwager, T., and Pöschel, T., 2007, “Coefficient of Restitution and Linear-Dashpot Model Revisited,” Granular Matter, 6(9), pp. 465–469.
Ramírez, R., Pöschel, T., Brilliantov, N. V., and Schwager, T., 1999, “Coefficient of Restitution of Colliding Viscoleastic Spheres,” Phys. Rev. E, 60(4), pp. 4465–4472.
Gilardi, G., and Sharf, I., 2002, “Literature Survey of Contact Dynamics Modelling,” Mech. Mach. Theory, 37(10), pp. 1213–1239.
Hunt, K. H., and Crossley, F. R. E., 1975, “Coefficient of Restitution Interpreted as Damping in Vibroimpact,” ASME J. Appl. Mech., 42(2), pp. 440–445.
Machado, M., Moreira, P., Flores, P., and Lankarani, H. M., 2012, “Compliant Contact Force Models in Multibody Dynamics: Evolution of the Hertz Contact Theory,” Mech. Mach. Theory, 53, pp. 99–121.
Lankarani, H. M., and Nikravesh, P. E., 1990, “A Contact Force Model With Hysteresis Damping for Impact Analysis of Multibody Systems,” ASME J. Mech. Des., 112(3), pp. 369–376.
Marhefka, D. W., and Orin, D. E., 1999, “A Compliant Contact Model With Nonlinear Damping for Simulation of Robotic Systems,” IEEE Trans. Syst. Man Cybern., Part A: Syst. Humans, 29(6), pp. 566–572.
Gonthier, Y., Mcphee, J., Lange, C., and Piedbœuf, J. C., 2004, “A Regularized Contact Model With Asymmetric Damping and Dwell-Time Dependent Friction,” Multibody Syst. Dyn., 11(3), pp. 209–233.
Flores, P., Machado, M., Silva, M. T., and Martins, J. M., 2011, “On the Continuous Contact Force Models for Soft Materials in Multibody Dynamics,” Multibody Syst. Dyn., 25(3), pp. 357–375.
Roy, A., and Carretero, J. A., 2012, “A Damping Term Based on Material Properties for the Volume-Based Contact Dynamics Model,” Int. J. Non-Linear Mech., 47(3), pp. 103–112.
Haddadi, A., and Hashtrudi-Zaad, K., 2008, “A New Method for Online Parameter Estimation of Hunt–Crossley Environment Dynamic Models,” Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2008), Nice, France, September 22–26, pp. 981–986.
Zhang, Y., and Sharf, I., 2009, “Validation of Nonlinear Viscoelastic Contact Force Models for Low Speed Impact,” ASME J. Appl. Mech., 76(5), p. 051002.
Sørensen, S., Hansen, M., Ebbesen, M., and Mouritsen, O., 2011, “Implicit Identification of Contact Parameters in a Continuous Chain Model,” Model. Identif. Control, 32(1), pp. 1–15.
Gharib, M., and Hurmuzlu, Y., 2012, “A New Force Model for Low Coefficient of Restitution Impact,” ASME J. Appl. Mech., 79(6), p. 064506.
Shivaswamy, S., and Lankarani, H. M., 1997, “Impact Analysis of Plates Using Quasi-Static Approach,” ASME J. Mech. Des., 119(3), pp. 376–381.
Minamoto, H., and Kawamura, S., 2009, “Effects of Material Strain Rate Sensitivity in Low Speed Impact Between Two Identical Spheres,” Int. J. Impact Eng., 36(5), pp. 680–686.
Adams, M., Lawrence, C. J., Urso, M. E. D., and Rance, J., 2004, “Modelling Collisions of Soft Agglomerates at the Continuum Length Scale,” Powder Technol., 140(3), pp. 268–279.
Xiong, X., Kikuuwe, R., and Yamamoto, M., 2012, “A Differential-Algebraic Contact Model With Nonlinear Compliant,” Proceedings of ASME Dynamic Systems and Control Conference, Fr. Lauderdale, FL, October 17–19, ASME Paper No. DSCC/MOVIC2012-8642.
Lankarani, H. M., and Nikaravesh, P. E., 1994, “Continuous Contact Force Models for Impact Analysis on Multibody Systems,” Nonlinear Dyn., 5(2), pp. 193–207.
Ismail, K. A., and Stronge, W. J., 2008, “Impact of Viscoplastic Bodies: Dissipation and Restitution,” ASME J. Appl. Mech., 75(6), p. 061011.
Yigit, A. S., Christoforou, A. P., and Majeed, M. A., 2011, “A Nonlinear Visco-Elastoplastic Impact Model and the Coefficient of Restitution,” Nonlinear Dyn., 66(4), pp. 509–521.
Xiong, X., Kikuuwe, R., and Yamamoto, M., 2013, “A Differential-Algebraic Method to Approximate Nonsmooth Mechanical Systems by Ordinary Differential Equations,” J. Appl. Math.2013, p. 320276.
Flores, P., Leine, R., and Glocker, C., 2010, “Modeling and Analysis of Planar Rigid Multibody Systems With Translational Clearance Joints Based on the Non-Smooth Dynamics Approach,” Multibody Syst. Dyn., 23(2), pp. 165–190.
Leine, R. I., and Nijmeijer, H., 2004, Dynamics and Bifurcations of Non-Smooth Mechanical Systems, (Lecture Notes in Applied and Computational Mechanics, Vol. 18), Springer-Verlag, Berlin.
Brogliato, B., Daniilidis, A., Lemaréchal, C., and Acary, V., 2006, “On the Equivalence Between Complementarity Systems, Projected Systems and Differential Inclusions,” Syst. Control Lett., 55(1), pp. 45–51.

## Figures

Fig. 1

The graphs of dio(x) and max(0,-x)

Fig. 2

Force-indentation curves: (a) typical empirical result adopted from, e.g., Fig. 4 in Ref. [1], Fig. 2.6 in Ref. [2], and Fig. 2 in Ref. [5], (b) the KV model (5), (c) the HC model (6) without any external force (solid line) and with an external pulling force (dashed line), and (d) the authors’ previous contact model (7)

Fig. 5

Influence of β1 on the behaviors of the new model (10) integrated by the RK4 with a time step size of 0.001 s. The parameters are set as: fe = 0 N, λ = 1.5, and K = 104 N/m1.5. The initial conditions are set as: p(0) = -0.1 m, p·(0) = 2 m/s, and a(0) = 0m1.5.

Fig. 6

Influence of β2 on the behaviors of the new model (10) integrated by the RK4 with a time step size of 0.001 s. The parameters are set as: fe = 0 N, λ = 1.5, and K = 104 N/m1.5. The initial conditions are set as: p(0) = -0.1 m, p·(0) = 2 m/s, and a(0) = 0m1.5.

Fig. 8

(a), (b) The COR as a function of β1, β2 and the impact velocity obtained from the new model (10). (c), (d) The COR as a function of the impact velocity. The parameters are set as: fe = 0 N, λ = 1.5, and K = 104 N/m1.5. The initial conditions are set as: p(0) = -0.1 m and a(0) = 0m1.5.

Fig. 7

Influence of γ on the behaviors of the new model (10) integrated by the RK4 with a time step size of 0.001 s. The parameters are set as: fe = 0 N, λ = 1.5, and K = 104 N/m1.5. The initial conditions are set as: p(0) = -0.1 m and a(0) = 0m1.5.

Fig. 9

(a) The COR as a function of γ and the impact velocity obtained from the new model (10). (b) The COR as a function of the impact velocity. The parameters are set as: fe = 0 N, λ = 1.5, and K = 104 N/m1.5. The initial conditions are set as: p(0) =-0.1 m and a(0) = 0m1.5.

Fig. 3

The force-indentation curves of the new model (10) and the HC model (6) integrated by the RK4 with a time step size of 0.001 s. The parameters are chosen as: fe = 0 N, λ = 1.5, K = 104 N/m1.5, γ = 2×103s-1, β1 = 3×10-3 s, β2 = 0.1 s/m1.5, and b2 = 0.35 s/m. The initial conditions are set as: p(0) = -0.1 m, p·(0) = 2 m/s, and a(0) = 0m1.5.

Fig. 4

Simulation of the bouncing motion by using the new model (10) and the HC model (6) integrated by the RK4 with a time step size of 0.001 s. The parameters are set as: fe = Mg, λ = 1.5, K = 104 N/m1.5, γ = 2×103s-1, β1 = 1.45×10-3 s, β2 = 0.2 s/m1.5, and b2 = 0.2 s/m. The initial conditions are set as: p(0) = -0.5 m, p·(0) = 0 m/s, and a(0) = 0m1.5.

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections