求解带摩擦接触问题的拟高斯迭代法.docVIP

求解带摩擦接触问题的拟高斯迭代法 （标题2、三号、黑体、加粗） 金朝 作者2 （宋体、五号字） 作 者 （宋体、五号字） 中国铁道科学院 （宋体、五号字） 单 位2 （宋体、五号字） （正文字体为宋体、五号字，行间距为18磅） 摘 要：（小四号字体加粗）本文提出了一种新的求解带摩擦的接触问题的数值算法，即拟高斯迭代法。它对法向接触力和切向接触力进行交替迭代，并利用高斯迭代法求解法向接触力，利用分块高斯迭代法求解切向接触力。同其它的数值算法相比，该算法保留了关于接触力的柔度矩阵的稀疏性和对称性，利用矩阵乘向量可以分步进行的技巧，该算法只存储关于接触力的柔度矩阵的下三角形矩阵的非零元和对角矩阵。根据可能接触边界的分布特点，将区域分解成不同的子区域，引进拉格朗日乘子表示接触力，保证了各子区域的网格剖分和位移求解是完全独立的。基于上述算法和有限元程序自动生成系统开发了相应的求解带摩擦接触问题的软件，数值实验表明，程序是正确的，算法是高效的。（字体楷体_GB2312、五号字） 关键词：摩擦，接触问题，拟高斯迭代法，拉格朗日乘子，区域分解（请选出4至5个关键字，如：有限元、接触、FEPG、并行） Abstract（小四号字体加粗） A numerical code RFPA3D based on linear elastic damage mechanics on mesoscopic scale and FEM are developed to simulate the complete progressive 3-D failure and macroscopic mechanical behaviours of rock-like materials. Rock properties such as elastic constants, peak strength, and Poisson ratio are randomly distributed to reflect the initial random distributed weakness in mesoscopic scale. By introducing heterogeneity of rock properties into the model, the RFPA3D code can simulate non-linear deformability of a quasi-brittle behaviour with an ideal brittle constitutive law. Each element in rock samples is considered to be homogenous and obeys a simple constituent law. A self-adaptive loading method is adopted to control the number of elements as less as possible that reach the failure threshold. A specimen was subjected to uniaxial tension. Complete stress-strain curves are obtained. Different from 2-D failure, 3-D samples show more ductile features. Crack initiation, propagation, coalescence and stress distribution can all be observed obviously during progressive failure process. （字体楷体_GB2312、五号字） Key words: （小四号字体加粗） impact，thin-walled tube，dynamic progressive buckling，non-linear FEM（字体楷体_GB2312、五号字） （注意：本模板仅作为对论文字体、段落、图片、参考文献和标题等的示范。除中英文的摘要和关键字必须编写之外，论文文章结构由作者自行决定。本文的文章结构不属于我们对论文编写的规定之一） 1、引言（大标题顶格、标题三、四号字、加粗） 接触问题在工程和自然界中大量存在，如高速火车的弓网接触问题，带断层或裂缝的拱坝稳定性问题，地震震源破裂和滑坡问题等。其特点和难点是接触边界和接触力事先是未知的，当考虑摩擦滑动时，能量变分原理不再成立［1］。目前求解接触问题的方法有解析法[2][3]，半解析法[4][5]