首页/文章/ 详情

子弹打靶—复合材料受高速冲击的应变率效应VUAMT

2小时前浏览34


对于静力线弹性问题,应力应变增量呈现线性关系,简单理解就是:

DStress=K*DSrain

然而试验中发现,高速冲击作用下,材料的模量和强度有了一定程度的“增强”,并且增强效果与应变率有关。所谓应变率就是应变的变化率:DSrain/Dt。这时,材料的刚度K会受到应变率的影响,见下图。

   

材料应变率效应

图来源:杨东, 姜紫薇, 郑志军. 高温高应变率下钛合金Ti6Al4V的动态力学行为及本构关系[J]. 高压物理学报, 2024, 38(1): 014101. doi: 10.11858/gywlxb.20230743


造成这个现象的原因比较复杂,应力在材料内部是以波的形式传播的,高速问题中还要考虑能量守恒,且不同的材料阻抗等参数也不同,这些因素共同导致冲击速度影响了材料的本构。


本文就以弹丸穿透靶板的案例介绍一下,如何基于VUMAT实现应变率本构和冲击失效效果。


VUMAT和UMAT

VUMAT整体功能与UMAT类似,主要区别是前者适用于动力学问题,且可以获取到上一时间步的应力、应变。比如stressOld就代表前一步应力,需要更新的当前应力则为stateNew。


还有一点,VUMAT需要完成能量密度的更新,否则可能不收敛。


本构模型

应力应变基本关系如下:


   
   

采用Hashin准则进行失效判断:

   

对于失效单元,进行刚度折减,实现结构整体渐进失效效果

   

采用如下形式描述应变率效应,其中X代表模量或者强度参数,下标o代表准静态参数值,下标rate代表当前应变率下的材料参数值。为系数由试验测定。


   

将上述模型嵌入到VUMAT中。


有限元模型

模型中包含复合材料靶板、弹丸,其中弹丸设置为刚体。靶板采用多层实体单元,注意不要使用ABAQUS的铺层定义模块,这会导致子程序无法使用。

   

网格部分,需要打开靶板的单元删除设置,否则将无法模拟出单元删除效果。另外需要设置接触关系,适应冲击接触。


仿真效果

失效区域扩展过程:

   

弹丸穿透过程:

   


01


学术道德规范


我方仅提供技术方面的实现与咨询服务,不涉及具体论文、专利等知识产权方面的工作。

02


不涉密承诺


我方不承接涉密技术服务。


来源:静界有限元
Abaqus复合材料二次开发UM材料试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-12
最近编辑:2小时前
静界有限元
博士 签名征集中
获赞 23粉丝 0文章 7课程 0
点赞
收藏
作者推荐

『静界有限元』工作室介绍

静界有限元工作室介绍part1工作室简介静界有限元1面向客户工作室主要面向在校学生、科研院所老师提供结构有限元仿真(含二次开发)、流体力学仿真、算法开发、软件开发服务。2服务内容结构有限元:复合材料失效分析、机织复合材料多尺度仿真、高速冲击仿真、疲劳分析、UMAT、VUMAT流体力学仿真:外流场仿真、单向流固耦合仿真、结冰冰形计算算法开发:数值计算、仿真结果二次处理软件开发:Matlab语言、Python语言part2工作室特色01深厚的FEM求解器技术积淀工作室长期致力于自研有限元求解器研制,对有限元方法有独到理解。可以更好的帮助客户利用现有的商用软件完成本构模型的二次开发,实现复杂问题的仿真。自研FEM求解器02大量UMAT/VUMAT开发经验针对材料静力失效、疲劳失效、低速冲击、高速冲击相关问题,工作室开发了大量的UMAT、VUMAT子程序,熟悉各种强度准则、渐进失效的实现方法。UMAT失效仿真03坚实的工业软件开发基础工作室在各类工业软件开发中,掌握了三维渲染、算法、交互、设备通信技术,为商用软件二次开发奠定了基础,可以帮助客户提高仿真结果的展示度。自研3D模块对有限元结果二次处理part3业务能力介绍1方向一:结构有限元复合材料失效分析:开发UMAT子程序,嵌入Hahin、蔡-吴等强度准则,采用渐进损伤方法,实现材料强度预测、失效模式预测、失效扩展过程预测、失效区域预测。复合材料接头损失区域应力应变曲线机织复合材料多尺度仿真:掌握机织复合材料单胞建模方法,开展单胞模型的弹性常数仿真计算、失效强度预测(基于UMAT子程序),并进一步实现宏观模型的工况仿真。单胞、全尺寸模型仿真高速冲击仿真:开发VUMAT子程序,考虑应变率效应,嵌入Hahin、蔡-吴等强度准则,采用渐进损伤方法,通过单元删除技术,实现结构高速冲击失效过程模拟。子弹穿透复合材料靶板疲劳分析:开发UMAT子程序,嵌入Hahin、蔡-吴等强度准则以及材料疲劳寿命模型,采用渐进损伤方法,实现结构S-N曲线预测、失效区域扩展过程预测。S-N曲线2方向二:流体力学仿真外流场仿真:掌握FLUENT、CFX、Fensap等CFD软件技术,精通ANSA等CFD网格划分软件,实现外流场气动仿真计算。单向流固耦合仿真:掌握单向流固耦合分析技术,将CFD获取的压力、温度等数据加载到ABAQUS中,进一步完成结构变形分析。结构变形结果自研结冰计算软件:自研Spring2D翼型结冰软件,免费开放使用。软件可以在数秒内完成翼型冰型计算,精度比肩Lewice软件。自研结冰软件Fensap结冰计算:精通Fensap结冰计算软件,可以实现复杂三维模型的流场计算、水滴撞击系数、多步冰形计算推进。某探头冰形进气口水滴撞击系数计算3方向三:算法开发数值计算方法:掌握FEM、CFD面元法、复合材料力学算法、常规数值计算方法,可以支持用户算法层面的定制化需求。面元法数值计算方法:掌握FEM、CFD面元法、复合材料力学算法、常规数值计算方法,可以支持用户算法层面的定制化需求。参数化网格生成三维云图渲染4方向四:软件开发工业软件技术:掌握基于PyQT、Matlabapp界面开发技术,熟悉Matlab、Python语言,精通设备通信采集技术、数据库技术、模块集成技术,支持多个工业软件项目落地。Matlab界面软件数据采集软件(基于Python)part4合作模式01技术服务针对客户具体仿真、计算、开发需求,完成相关工作,并撰写技术报告。02教学针对用户具体技术需求,开展一对一教学服务。03企业服务针对企业的重要项目需求,以外协单位方式承担技术任务,通过合同和技术协议进行工作内容、周期的强约束。来源:静界有限元

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈