延性单晶非均匀变形的三维有限元模拟  

Three Dimensional FEA Simulation of Nonuniform Deformation in Ductile Single Crystals

作　　者：冯露[1] 张克实[1] 张光[1] 

机构地区：[1]西北工业大学,西安710072

出　　处：《应用力学学报》2002年第4期5-9,共5页Chinese Journal of Applied Mechanics

基　　金：国家自然科学基金 19972 0 5 5 ;航空基础科学基金 0 0C5 3 0 2 2;西北工业大学博士创新基金

摘　　要：对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6 .1商用有限元软件中。该方法的特点是 ,利用晶体塑性的动力学方程 ,获得一个关于晶体弹性变形梯度的演化方程 ,采用半隐式积分方案进行求解。本文推导出一种新的应力应变本构矩阵。按此方式更新本构矩阵 ,计算速度和计算稳定性大大提高。加载方式 ,边界条件和变形程度等因素影响着滑移系的启动状况 ,这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况 ,全面地考虑了FCC单晶材料 12个可能滑移系在变形过程中的启动状况 ,合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。Nonuniform deformations of ductile single crystals subject to tensile loading are simulated by using three dimensional finite element methods. Rate dependent crystal plasticity model together with a new time integration algorithm has been implemented into ABAQUS6.1 finite element software. The feature of this algorithm is that an evolution equation for the crystal elastic deformation gradient is obtained from crystal kinematics, and an implicit scheme integrates this equation. A new constitutive matrix is presented. It is shown that constitutive matrix updating this way yields significant reduction in computational time. In simulation twelve possible slip systems for FCC single crystals are considered during practical deformations. Nonuniform deformations such as strain localization and necking are simulated rationally and really.

关 键 词：延性单晶 非均匀变形 三维有限元模拟 晶体粘塑性 数值积分方法 

分 类 号：O733[理学—晶体学]