Tibbitts Ira B, Kakarla Deepika, Siskey Stephanie, Ochoa Jorge A, Ong Kevin L, Brannon Rebecca M
University of Utah, Dept. of Mechanical Engineering, 2134 MEB, 50 S Central Campus Dr., Salt Lake City, UT 84112.
Ortho Development Corp., 12187 Business Park Dr., Draper, UT 84020.
Exp Mech. 2014 Feb 1;54(2):187-198. doi: 10.1007/s11340-013-9778-6.
To validate models of contact mechanics in low speed structural impact, slender rods were impacted in a drop tower, and measurements of the contact and vibration were compared to analytical and finite element (FE) models. The contact area was recorded using a novel thin-film transfer technique, and the contact duration was measured using electrical continuity. Strain gages recorded the vibratory strain in one rod, and a laser Doppler vibrometer measured speed. The experiment was modeled analytically on a one-dimensional spatial domain using a quasi-static Hertzian contact law and a system of delay differential equations. The three-dimensional FE model used hexahedral elements, a penalty contact algorithm, and explicit time integration. A small submodel taken from the initial global FE model economically refined the analysis in the small contact region. Measured contact areas were within 6% of both models' predictions, peak speeds within 2%, cyclic strains within 12 με (RMS value), and contact durations within 2 μs. The global FE model and the measurements revealed small disturbances, not predicted by the analytical model, believed to be caused by interactions of the non-planar stress wavefront with the rod's ends. The accuracy of the predictions for this simple test, as well as the versatility of the diagnostic tools, validates the theoretical and computational models, corroborates instrument calibration, and establishes confidence that the same methods may be used in experimental and computational study of contact mechanics during impact of more complicated structures. Recommendations are made for applying the methods to a particular biomechanical problem: the edge-loading of a loose prosthetic hip joint which can lead to premature wear and prosthesis failure.
为验证低速结构冲击下的接触力学模型,在落塔中对细长杆进行冲击,并将接触和振动测量结果与解析模型和有限元(FE)模型进行比较。使用一种新型薄膜转移技术记录接触面积,利用电连续性测量接触持续时间。应变片记录一根杆中的振动应变,激光多普勒振动计测量速度。该实验在一维空间域上采用准静态赫兹接触定律和延迟微分方程组进行解析建模。三维有限元模型使用六面体单元、罚函数接触算法和显式时间积分。从初始全局有限元模型中提取的一个小子模型在小接触区域经济地细化了分析。测量的接触面积在两个模型预测值的6%以内,峰值速度在2%以内,循环应变在12 με(均方根值)以内,接触持续时间在2 μs以内。全局有限元模型和测量结果揭示了解析模型未预测到的小扰动,认为是由非平面应力波前与杆端的相互作用引起的。这个简单测试的预测准确性以及诊断工具的通用性验证了理论和计算模型,证实了仪器校准,并建立了信心,即相同的方法可用于更复杂结构冲击过程中接触力学的实验和计算研究。针对将这些方法应用于一个特定生物力学问题提出了建议:松动的人工髋关节边缘加载,这可能导致过早磨损和假体失效。