粘弹性接触半解析算法研究.docVIP

  粘弹性接触半解析算法研究# 彭向征，王战江** （重庆大学机械传动国家重点实验室重庆 400030） 5 10 15 20 25 摘要：提出了一种基于半解析法的粘弹性接触快速算法。在接触求解中运用最小余能原理， 并使用共轭梯度法(CGM)和快速傅里叶变换(FFT)。由于蠕变效应具有路径相关性，需计算 时考虑前步求解结果。通过分析三参数固体材料以及 Maxwell 材料的接触问题，验证了模 型的正确性。在 Maxwell 材料的接触分析中，最大接触面积与最大载荷不一致，最大接触 面积滞后于最大载荷，载荷卸载时压入面积快速下降为零。该方法可以用来求解橡胶等粘 弹性材料的接触问题，并可进一步分析橡胶纹理的接触问题。 关键词：接触力学；半解析法；粘弹性接触；共轭梯度法 中图分类号：O357 Visco-elastic contact analysis based on Semi-analytical method PENG Xiangzheng, WANG Zhanjiang (State Key Laboratory of Mechanical Transmission,Chongqing University,Chongqing 400030) Abstract:In this paper, a numerical approach based on Semi-analytical method is proposed for the visco-elastic. The contact solver is based on minimization of complementary energy, which integrated with the conjugate gradient method (CGM) algorithm and the fast Fourier transform (FFT) technique for improving numerical efficiency. The effect of creep is considered by using results from previous step. Good agreements can be observed by comparing the results from present method and analytical solutions for the 3-parameter solid materials and the Maxwell materials. Keywords: contact mechanics, semi-analytical method，visco-elastic contact，conjugate cradient method 0 引言 在粘弹性接触力学研究方面，Lee 和 Radok[1]建议用比拟的方法，将粘弹性应力应变关 30 35 40 系式中相应的积分运算代替弹性解中的弹性常数实现粘弹性接触的计算。他们进一步指出， 这种方法可以应用于接触面积增加的加载过程。但是这种方法应用于接触面积减少的情况 时，接触区外部也将有接触压力，这样积分运算难以完成卸载过程的求解。Ting[2]提出一种 方法来求解卸载过程，假设在时刻 t 时，接触尺寸 a(t)是减小的，而前面必定有一个时刻 t1， 其接触尺寸 a(t1)是增加并等于 a(t)。通过将时间分为 0-t1 和 t1-t 两部分进行积分计算得到相 关压力和压入量。但是这种方法公式复杂。最近，Greenwood[3]提出通过叠加一系列刚性平 底压头和粘弹性半无限大体解的方法得到任意轴对称形貌粘弹性接触的解。进一步，他用这 种方法分析刚性圆锥形压头、球形压头与具有三参数固体模型特性的半无限大体的压入和卸 载过程。这种解析方法较 Ting 的方法简单而能求解多种轴对称问题接触的解。Goryacheava[4] 对层状材料的粘弹性问题进行了研究，并得出了解析解，但是所使用粘弹性模型及层状材料 边界条件的假设过于简单。 基金项目：基金项目：教育部博士点科研基金（20110191120009） 作者简介：彭向征（1987-）男，研究生，非线性接触力学 通信联系人：王战江（1981-）男，副教授，粘弹性力学. E-mail: pxztop2@163.com -1-   Naghie 等[5]研究了不同泊松比下的粘弹性接触问题，但是泊松比不随时间发生变化。上 述研究者均集中在研究解析解。简单形貌的粘弹性接触可以得到解析解，但是对于任意形貌 45 50 55 的接触问题，特别是具有粗糙度的表面，需要借助数值的方法求解。此外，对于变泊松比问