基于vua的圆柱壳外爆损伤数值模拟.docxVIP

基于vua的圆柱壳外爆损伤数值模拟 由于abaqus的强大非线性功能，它可以模拟高速影响（例如，弹体破裂等）的大变形问题，这是一种越来越多的用户使用的。abaqus有一个相对完整的材料数据库，可以定义不同属性的材料。在材料损伤强度的学方面，abaqus也包含一些相关内容，但有限，使用起来不方便。作者利用abaqus的两个开发软件，将abaqus损坏规范引入两次开发软件，计算并显示解算方法能够有效模拟材料的损坏和损坏。 1 生成应力的计算 考虑含大量孔洞、体积为V的一块介质,假设实体的体积为Vs,孔洞的总体积为Vd,则有V=Vs+Vd.定义损伤D=Vd/V,可见损伤是无量纲的,本文假设材料损伤的主要形式是延性微空洞损伤,采用张伟推导出的损伤演化方程: ?D?t={34σ2s?σ20KsλclD(1?D)3/2?Dε˙pkk0σsσ0σs≤σ0?D?t={34σs2-σ02ΚsλclD(1-D)3/2-Dε˙kkpσsσ00σs≤σ0 式中,σs为实体应力,σs=σ/(1-D);σ0为产生损伤的阈值拉应力;Ks为体积模量;λ为单位面积表面能;cl为体波波速;ε˙ε˙pkk为塑性体应变率.在计算中用静水压力代替σs.同时,材料本构也考虑了损伤的影响,本文利用文献中提出的“应变等价原理”导出有关材料参数的损伤软化形式.应变等价原理可表达为“任何损伤材料的本构方程与无损伤材料均一致,只需将通常应力用实体应力替代”.下面给出2种参数含损伤的变化情况: 式中,E、K分别为考虑损伤影响后材料的弹性模量和体积模量. 2 自适应前应力指数自适应cog 在应力空间中材料的含损伤热塑性屈服准则为f(σ,εp,β,D,T)=0,式中,β为各向同性硬化参数,εp为塑性应变张量.本构关系采用流体弹塑性模型.应变增量为体应变增量和偏应变增量之和,体应变可用比容的变化表示.应力为静水压和偏应力之和.静水压p和比容v之间的关系满足Murnagham方程: p=Kγ[(v0v(1?D))γ?1p=Κγ[(v0v(1-D))γ-1 式中,γ为Murnagham指数. 假设材料满足各向同性硬化的Mises屈服准则,其屈服函数具体可写为f=(3/2)sijsij-k2=0,式中,sij为偏应力张量,k为屈服应力和塑性功的函数.则由广义Drucker公式可导出: dεpij=dλ?f?σij=3sijdλ(2)dεijp=dλ?f?σij=3sijdλ(2) 式中,dεpij为塑性应变增量,dλ为塑性流动因子,dsij为偏应力增量,deij为偏应变增量,depij为塑性偏应变增量,G为剪切模量. 3 材料及数的材料常数以及数的材料常数 ABAQUS为了使其具有更强的适用性和灵活性,在材料属性方面提供了可供用户定义的子程序,针对ABAQUS/Standard和ABAQUS/Explicit2个不同分析模块,ABAQUS有2个定义用户子程序方式:UMAT和VUMAT.本文分析的是材料动态损伤破坏问题,因此采用Explicit分析模块,用VUMAT对材料进行定义. 定义的过程:①在ABAQUS输入文件中定义材料所含有与求解相关变量的个数以及材料常数;②用FORTRAN或C语言编写材料本构程序,即在已知应变增量的基础上,给出材料的应力增量.ABAQUS在每个时间增量步的材料积分点均会调用VUMAT对材料进行求解计算.本文本构的特点,一是材料为流体弹塑性模型,二是要考虑损伤及温度对材料本构的影响.子程序实现过程如下(其中上标n+1和n分别表示新旧时刻): ①假设从n时刻到n+1时刻是弹性变形,塑性偏应变增量和塑性功增量都为0,计算n+1时刻的温度Tn+1,求新的材料参数En+1、Kn+1等. ②利用Murnagham方程,计算n+1时刻的单元静水压p. ③在已知应变增量dεij基础上,计算从n时刻到n+1时刻的偏应力增量dsij,求出n+1时刻的应力偏量(sij)n+1. ④计算fn+1,如果fn+10,则转步骤⑤,否则分下面2种情况处理: a.若fn= 0,则由n时刻到n+1时刻为塑性加载,计算出dλ、塑性偏应变增量(及塑性偏应变累积值)与塑性功累积值,重新求新时刻的温度Tn+1和新的材料参数Kn+1、Gn+1、Yn+1等.令fn+1=0,求出n+1时刻的偏应力(sij)n+1,然后转步骤⑤. b.如果fn0,则由n时刻到n+1时刻为部分塑性加载,调整n时刻的偏应力使fn=0,计算出从fn=0到n+1时刻的应变增量,转步骤④的a. ⑤给出此次时间步长增量后,计算应力张量各个分量, 退出子程序调用.图1是程序实现流程. 4 材料参数的选取 考虑一无限长圆柱壳受内部爆炸载荷的破裂情况.图2为圆柱壳的横截面图,内部为半径R1的圆柱体炸药,外层为厚度R2-R1的圆柱壳体.炸药从中心点起