基于能量极值原理的单晶体塑性滑移的有限变形分析Ξ.PDFVIP

Vol. 23 No. 3 　　 　　　　　　 　 　　　　　　　　　 　 第23卷 第3期 固 体 力 学 学 报 2002年 　9 月 ACTA MECHANICA SOLIDA SINICA September 2002 基于能量极值原理的单晶体塑性 滑移的有限变形分析 冯 　露 　　张克实 　　张 　光 ( 西北工业大学工程力学系 ,西安 ,710072) 摘 　要 　对 FCC 单晶体的率无关弹塑性力学响应的本构关系进行了数值模拟. 用一个基于 能量极值原理的数值计算方法来处理复杂的多面塑性问题. 这种算法可以有效地模拟单晶体多 滑移系的启动. 在这一理论框架下 ,增量的应力应变关系可以从所构造的能量函数中推导出来. 对于滑移系激活情况的判定则可转化为求解活动约束的非线性数学规划问题. 通过对时间的离 散 ,此问题又可细化为逐步二次规划问题 ,并采用有效集法来搜索启动滑移系 ,进而求得弹塑性 本构关系. 数值结果表明该方法具有稳定、收敛、可行的特点. 在数值计算的基础上研究了单轴 拉伸下晶体的不同取向对单晶体硬化程度和滑移系激活情况的影响. 关键词 　晶体塑性 ,弹塑性有限变形 ,多滑移系启动 ,非线性规划 ,有效集法 1 　引言 对于不均匀材料例如多晶体的宏观热力学性质的描述通常建立在至少两个唯象尺度上 即宏观描述和微观描述. 宏观尺度的研究与均匀的连续介质相关联 ,而微观尺度的研究往往 是选取具有统计特性的代表性体积单元 ,它包含着微观成分 ,例如多晶体中具有代表性的微 元就是单晶体. 为了预测晶体材料的组织演化过程 ,宏观模型具有难以克服的局限性 ,例如 ( ) 它们难以预测由于晶体结构变形 如晶格旋转 导致的各向异性性质. 较完善的模型通常考 虑了由于随动强化所导致的 Bauschinger 效应 ,却忽略了变形中屈服面形状的改变. 与此相 反 ,具有代表性的微观体积单元 ———晶体的本构模型得到了较为完善的发展. 近年来所谓的 连续滑移理论得到了一致的认可[1～5 ] . 建立在位错动力学模型上的晶体滑移硬化规律 ,在文 献中多有报导[6 ] . 这些理论为将连续滑移理论应用于描述晶体的非弹性变形奠定了良好的 基础. 在这一研究背景下 ,有必要发展一种适用于晶体塑性的稳定的应力计算方法. 对单晶或多晶体的力学行为响应的描述大多是在粘塑性框架下进行的. 在这种情况下 , 对晶体滑移系滑移应变的计算是通过求解一组非线性方程来完成的. 事实上 ,所有的滑移系 并不是同时启动的 ,而是根据作用在每一滑移系上的分切剪应力来判定滑移系的激活情况. 在弹塑性率无关框架下 ,确定滑移系的启动情况非常困难. 这一问题可以归属于多面塑性这 一数学框架. 在这种情况下 ,这些约束构成给定滑移面上的屈服函数. 与标准的多面塑性 ( ) multi - surface plasticity 所不同的是 ,这些约束又是线性相关的. 这将导致在给定变形模式 ( ) ( ) 国家自然科学基金 、航空科学基金 00C53022 和西北工业大学博士论文创新基金资助. 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 3 期 　 　冯 　露等 : 　基于能量极值原理的单晶体塑性滑移的有限变形分析