我不是专门研究裂纹扩展的,前不久一个偶然的机会涉及到这一方面,查阅了一些资料。先将我查阅的资料整理如下,希望对大家有帮助,请大家批评指正,加以完善。
首先必须把单元间共节点的节点离散,可以采用 LS-Prepost 或 femp 实现。然后在通过 MATLAB 或者其他语言编写小程序,对位于同一个位置的节点建立节点集,添加 *CONSTRAINED_TIED_NODES_FAILURE 关键字。采用此方法来实现裂纹模拟的缺点是前处理太麻烦。应用实例可参考白金泽《LSDYNA3d 基础理论与实例分析》。
关于这个关键字本版内有很多讨论,可以搜索一下。需要注意的是,在 LSDYNA 971R4 之前的版本中,这个材料模型所带的失效模式均只适用于单点积分的二维和三维实体单元。但是在 R4 之后的版本中,这个关键字有了很大的改进:
1、去除了单点积分的限制,同时还支持 3 维壳单元和厚壳单元中的 type1 和 type2。
2、可以定义初始损伤值,增加了几种损伤模型,具体可以参考 LSDYNA 971R5 版的关键字。
有些材料模型本身就带有失效的,可以定义单元的失效来模拟裂纹的拓展。如 *MAT_PLASTIC_KINEMATIC 等。如果某些材料模型不带失效模式,可以采用方法 2,或者通过自定义材料本构来实现裂纹的模拟。
这个方法个人觉得有些牵强,但是在有些文献中也见过。在定义裂纹前必须已知可能出现裂纹的区域,通过带有失效模式的面对面的绑定接触 CONTACT_TIED_SURFACE_TO_SURFACE_FAILURE 或者用弹簧单元来模拟裂纹面。
某些材料模型如 *MAT_120(*MAT_GURSON),*MAT_120_JC(*MAT_GURSON_JC),*MAT_120_RCDC(*MAT_GURSON_RCDC),还有一些 damage 模型,如 *MAT_96(*MAT_BRITTLE_DAMAGE)等,用损伤值来代替裂纹,通过观察损伤云图来判断裂纹的扩展。
[版权声明] :本文文字、代码及图片版权归原作者所有,任何媒体、网站或个人未经本网协议授权不得采集、整理、转载或以其他方式复制发表。已经本站协议授权的媒体、网站,在使用时必须注明“稿件来源:学研谷”。