薄膜破坏过程数值模拟的MPM方法  

NUMERICAL SIMULATION OF THIN FILM FAILURE WITH MPM

作　　者：郑勇刚[1] 顾元宪[1] 陈震[2] 

机构地区：[1]大连理工大学工程力学系 [2]美国密苏里大学哥伦比亚分校土木工程系

出　　处：《力学学报》2006年第3期347-355,共9页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家自然科学基金(10421002;10032030)海外青年学者合作研究基金(10228206)资助项目.

摘　　要：将连续与不连续本构模型相结合用于薄膜破坏过程的模拟计算．采用von Mises关联塑性本构模型描述材料的弹塑性变形过程,在塑性硬化和软化阶段通过检查分岔是否发生,从而判断是否有破坏出现,即材料出现应变局部化现象．一旦出现破坏,则在发生破坏的区域采用位移不连续(decohesion)本构模型进行模拟, 直到破坏面两侧材料完全分离．分离阶段采用相应的分离算法同时运用黏性边界条件以提高计算效率．结果表明结合连续与不连续本构模型的方法,能很好地模拟材料从局部化到完全破坏的过程,体现了材料破坏连续建模的思想,可用于预测材料发生破坏的环境和形式,从而通过改变材料工作环境提高其使用寿命;同时数值算例也显示了物质点法(material-point-method,简称MPM)的健壮性和有效性．Thin metallic films under compressive loads often fail due to delaminating away from the substrates, which has been observed in many experiments and arouses many researchers＇ interests. Recently, much research has been conduced to model and simulate this physical phenomenon from theoretical and computational aspects and most of them are based on the conventional elastic stability theory and interracial fracture mechanics （Shen 2004）. In fact, the thin film failure often involves both localization and decohesion processes and the evolution of localization in the film is quite possible is the reason of delamination （Chen 2003）, which is characterized by the interaction between geometrical and material instabilities. To simulate such a multi-physical phenomenon, a careful designed model should be constructed. Before localization occurs, thin films undergo usually elastoplastic deformations, and after that displacement jump appears in some local areas, which can be modeled by a decohesion-based model （Schreyer et al. 2002, Chen et al. 2003）. It was pointed out that the initiation of localization could be regarded as the transition from continuous to discontinuous failure modes （Chen 1996）, which can be determined via a bifurcation analysis of the acoustic tensor. The purpose of this paper is to simulate thin film failure phenomenon with the combination of continuous and discontinuous constitutive models and the Material Point Method （MPM, Sulsky et al. 1994）. The MPM is chosen as a basic method in this paper because of its simplicity and high performance in modeling large deformation, impact/contact, blast, penetration and other challenging physical processes. At the beginning of the deformations of the segment, associated yon Mises elastoplasticity constitutive model is used. Whether a local failure will appear in the segment can be identified through the bifurcation analysis in the plastic deformation stage, namely, localization. If localization occurs, the decohesion-based model will be ado

关 键 词：物质点法 局部化 分岔 位移不连续 薄膜 

分 类 号：O344.4[理学—固体力学] O347.4[理学—力学]