- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
复杂加载模式下非线性混合强化材料模型的应用
非线性比例加载模拟的可行性
塑料形成的模拟加工精度主要取决于材料性能参数、结构关系的选择和边界条件的调整。随着计算机技术的发展, 模拟计算的效率问题已基本解决, 但模拟的精度问题仍然是相关学者们关注的热点。
在经典塑性理论中, 屈服面的作用是限定材料的弹性区域, 对其初始位置及后继变化的描述是塑性本构方程的重要内容。因为屈服面及其变化很难测定, 所以在理论塑性力学中经常使用等向强化法则 (Isotropic hardening law, 简称Iso) , 即假设初始屈服面在塑性变形过程中仅半径增大, 而圆心不变, 这种强化法则仅能够准确预测线性比例加载情况。对处理非线性比例加载情况时, 文献提出了线性随动强化理论 (Kinematic hardening law, 简称KIN) 。该理论假设塑性变形过程中初始屈服平面仅改变位置, 而不改变形状大小。此外, 为了同时描述塑性变形中应力屈服面的形状和位置变化, 常将等向强化法则与随动强化法则结合使用, 从而组成各向同性动态强化法则 (Isotropic-Kinematic hardening law, 简称ISO-KIN) 。为准确描述反向加载过程中的包辛格效应和快速应力应变强化现象, 文献提出了多种模型, 其中, Chaboche模型因使用单屈服面而比使用双屈服面的模型计算效率高。将这些模型和经典的ISO-KIN强化模型结合起来可以更好的描述反向加载过程中的包辛格效应和快速应力应变强化现象。
合理的硬化模型可以更好的描述变形路径对材料塑性变形行为的影响, 因此在虚拟加工过程中使用正确的应力应变本构关系, 是准确预测成形过程的前提。飞机蒙皮拉形过程中板料承受拉、弯、扭等多种复合运动, 拉形加载轨迹复杂多样。不同的空间加载轨迹决定了板料不同的应变路径。成形过程中零件内存在不同的应力、应变状态, 其复杂程度远非线性比例加载情况所能描述。但等向强化法则和线性随动强化法则因其使用方便仍被广泛应用于拉伸成形的有限元模拟研究中, 影响了模拟的准确性。为提高复杂加载路径下拉形工艺的预测精度, 本文分别对单凸蒙皮拉形以及复杂S形蒙皮拉形工艺过程进行研究, 从而得到Chaboche强化模型的适用条件及其能否准确使用的判定法则。
1 ziegerr模型
屈服面可描述为:
式中f———屈服方程表达式, 为m阶齐次方程
σ———屈服应力, 为柯西应力/MPa
α———背应力/MPa, 在背应力的作用下屈服表面由初始位置移动到当前的位置
———表示屈服面大小的等效应力, 是等效应变的函数, 表示等效应力把屈服面作为一阶齐次方程进行计算
由等效塑性功增量dwiso= (σ-α) dεp=σisodε, 对式 (1) 求导可得:
由Ziegler模型可知, dα~ (σ-α) ;另外, 通常假设dα~dεp。在Chaboche模型中, 背应力增量由两部分组成, 即dα=dα1-dα2, 其中dα1表示加载过程的背应力增量, dα2表示卸载过程的背应力增量, 为了区分加载和卸载 (反向加载) 过程中的快速应力应变强化行为, 则式 (2) 可改写为:
其中
将式 (4) 代入屈服函数f中, 可得:
则,
由式 (3) ~式 (6) 可得:
其中
则线弹性各向同性本构关系可以表达为:
式中C———ABAQUS/STANDARD用户子材料程序UMAT的单元刚度矩阵
对于每一个迭代时间步长内的应变增量dε, 由式 (7) 和式 (8) 可得塑性应变增量dεp, 背应力和屈服应力增量按式 (5) 、式 (6) 和式 (9) 计算, 依照Jaumann应力增量法则均为分别增加。
2 应力更新过程
在求解板料成形过程等非线性边值问题时, 加工过程模拟的每一步数值迭代通常包含3个步骤。第一, 通过尝试改变位移增量得到离散的应变增量;第二, 按照弹塑性本构方程由离散的位移增量求得应力及其他变量, 如塑性应变、背应力等的增量;第三, 检验需要求解的要素是否在误差允许范围内, 如果满足要求, 迭代过程结束, 否则就尝试产生新的应变增量使迭代过程循环。
在第二步应力更新过程中, 首先假定对于给定的应变增量Δε, 使更新后的应力满足弹性条件:
其中‘T’表示尝试状态, 下标表示应力更新时间步, 同时将上一步的应变增量和背应力保存到本步参加运算, 即:
评价应力状态是否为弹性状态的误差表示为Tole, 如果式 (12) 所示的屈服关系在该误差范围内, 则第n+1步迭代被认为是弹性状态。
如果式 (12) 所示应力状态不在误差范围内, 则该应力状态被认为是弹-塑性状态, 上一个弹性应力状态被认为是求解该弹塑性问题的初始状态。通过获取塑性应变增量得到式 (13) 所示的位于新屈服面上的塑性应力状态。
文档评论(0)