基于互补理论和扩展有限元法的摩擦接触裂纹应力强度因子计算  

Stress Intensity Factor Calculation for Cracks with Frictional Contact Based on Complementary Theory and Extended Finite Element Method

作　　者：郑安兴 罗先启[2] 钱镜林 Anxing Zheng;Xianqi Luo;Jinglin Qian(School of Water Conserancy and Environmental Engineering,Zhejiang University of Water Resources and Electric Porer,Hangzhou,310018;School of Naval Architecture,Ocean and Crvil Engineering,Shanghai Jiaotong Universitry,Shanghai,200240)

机构地区：[1]浙江水利水电学院水利与环境工程学院,杭州310018 [2]上海交通大学船舶海洋与建筑工程学院,上海200240

出　　处：《固体力学学报》2020年第1期50-58,共9页Chinese Journal of Solid Mechanics

基　　金：国家自然科学基金项目(51279100);浙江省自然科学基金项目(LQ18E090003)资助。

摘　　要：从虚功方程出发,结合扩展有限元离散技术与接触条件的非线性互补表述,建立了摩擦接触裂纹问题的扩展有限元非线性互补模型,将不等式接触条件转化为非线性互补类的非光滑方程组,并采用基于广义导数的非光滑阻尼牛顿法求解方程组,无需引入任何额外人工变量以及迭代求解.以含中心倾斜裂纹平板和边裂纹平板为例,运用相互作用积分法计算摩擦接触裂纹的应力强度因子,将其结果与理论解进行对比分析,该方法都能给出精确的计算结果;基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟,计算结果表明,受压裂纹数值结果与实验结果比较吻合,从而验证了论文方法的有效性与正确性.The extended finite element method is a numerical method for modeling discontinuities within the classical finite element framework. This method can treat arbitrary cracks independently of mesh and crack growth with minimal remeshing. Based on the principle of virtual work equation, and combined with the nonlinear complementary expression for contact conditions and the extended finite element discrete technique, a nonlinear complementary model for frictional contact crack problem with the extended finite element method is presented. First, the conditions that describe frictional contact are formulated as a system of non-smooth equations based on variational inequality theory, and the non-smooth damped Newton method is given based on the definitions of generalized derivative to directly solve the system of equations without any extra parameters and iteration. Then, taking the finite plane with a single inclined crack and the edge-cracked plate as examples, the stress intensity factor of finite plane with a closed crack is calculated by the interaction integral approach. The numerical results calculated by the extended finite element method are basically consistent with the theoretical solution. Finally, the propagation of an inclined crack under uniaxial compression is simulated by the extended finite element method. The numerical results are in agreement with the experimental ones, indicating that the proposed method can accurately simulate the crack growth under uniaxial compression. Numerical examples are presented to demonstrate the correctness and effectiveness of this method.

关 键 词：扩展有限元法 裂纹 摩擦接触 应力强度因子 非线性互补法 

分 类 号：O346.1[理学—固体力学] TB115[理学—力学] O313.5[理学—数学]