基于abaqus的折弯片断裂问题仿真.docxVIP

基于abaqus的折弯片断裂问题仿真 0 v形槽折弯片特性 在有限框架下，xferm在不连续的区域执行有限函数，并添加描述不连续边界的方法。有限软件abaqus扩展了基于预重叠方法的有限元法功能。假设虚拟节点法，并不具有表示非连续性位移场的性能，裂纹的传播仅限于网格。abaqus有限功能的动态扩展过程不需要重建网格。 V形槽折弯片为两面开槽的金属薄片, 金属片之间的过渡区域用圆弧连接, 端部施加局部压力, 上面的金属板施加转角位移.典型的V形槽折弯片见图1, 折弯片结构的参数示意见图2.图2中, Rr为折弯片连接处圆弧转角的半径, 是影响断裂的主要设计参数之一. 对下面的金属板固支约束, 上面的金属板施加向上35°转角.不同的Rr值, 折弯片的开裂时间和断裂扩展路径有所不同.采用Abaqus的XFEM功能对该问题进行分析, 忽略金属板上开孔对开裂的影响. 1 单元类型的影响 影响裂纹扩展的因素有过渡区圆弧半径、网格尺寸和单元类型等.根据裂纹的扩展方向和裂纹处的应力, 裂纹为I型和Ⅱ型裂纹叠加的混合裂纹.在轴力和位移载荷作用下, 不同的圆弧半径和结构的应力场会有差别, 裂纹开裂通过一个单元的时间会有所差别, 圆弧半径的敏感性分析见表1. 圆弧半径为Rr=0.2, 0.3和0.5 mm时裂纹的开裂与扩展云图见图3~5. Abaqus中的XFEM使用1阶单元, 单元类型有CPE4, CPS4, C3D4, C3D8以及它们的缩减积分和非协调单元.1阶单元和1阶缩减积分单元的裂纹路径相同, 折弯片采用二维平面应变 (缩减积分) 单元CPE4和CPE4R, 2种单元类型的计算结果比较见图6. 在一个增量步中, 裂纹扩展方向的改变角度不能超过90°.断裂区网格的细化和裂纹的扩展方向不稳定.裂纹在一个单元内部的改方向接近90°, 裂纹扩展超出断裂区, 导致结果不收敛.网格对裂纹扩展方向的影响见图7. 2 lefm方法 Abaqus的XFEM分为LEFM和cohesive element等2种算法, 其中, LEFM方法适用于弹性材料的分析, cohesive element适用于弹塑性材料的分析.计算时没有考虑材料的塑性, 采用二维模型对裂纹的开裂和扩展进行分析. LEFM方法需要在interaction模块中定义断裂的接触属性, 赋给断裂区域;而cohesive element方法是在材料中定义破坏准则, 赋给断裂区域. 如图8所示, cohesive element裂纹的扩展方向较LEFM方法平滑一些, 但扩展路径容易发生较大转角. 3 静态通用分析的计算 Abaqus的XFEM功能可以使用static general和dynamic implicit等2种分析步类型.隐式动态分析比通用静态分析所需的分析步多, 考虑质量的影响, 结果较准确.静态分析步中位移载荷和压力的加载方式分为线性加载与平滑加载, 线性加载容易造成不收敛或超出迭代次数, 而平滑加载容易收敛.在相同网格的条件下, cohesive element计算所需要的增量步在1 000步以内, 较LEFM方法少一些, 不同算法的增量步数目见表2.静态通用分析的计算成本较隐式动态分析的低一些, 但是会出现收敛问题. 4 拉格力的影响 折弯片的断裂仿真分析表明, 单元类型对裂纹的路径影响很小.裂纹相邻处单元的应力会出现棋格效应.在同样的载荷下, 折弯片连接处圆弧的半径影响裂纹的开裂和完全断开的时间, 细化网格和裂纹扩展的方向有任意性, 超出断裂区.cohesive element方法的计算结果较LEFM平滑一些, 所需要的增量步也少一些.