UMAT子程序在复合材料强度分析中的应用.docxVIP

UMAT 子程序在复合材料强度分析中的应用 讲师：孔祥宏版本：Abq。难度：关键词：强度分析，UMAT本例使用 UMAT 用户子程序进行复合材料单层板的应力分析和渐进损伤压缩强度分析，介绍UMAT 用户子程序编写方法及在Abaqus/CAE 中的设置。本章使用最大应变强度理论作为复合材料单层板的失效准 讲师：孔祥宏 版本：Abq 。 难度： 关键词：强度分 析，UMAT 知识要点： 强度分析 UMAT 用户子程序最大应变理论 刚度折减 .1 本章内容简介 本章通过两个实例介绍 UMAT 用户子程序在复合材料单层板的应力分析和强度分析中的应用。在第一个实例中，对一个简单的复合材料单层板进行应力分析，UMAT 子程序主要计算应力，不进行强度分析，本例用于验证UMAT 子程序的计算精度。在第二个实例中，对复合材料单层板进行渐进损伤强度分析，UMAT 子程序用于应力计算、强度分析和刚度折减。 本章所用复合材料为T700/BA9916，材料属性如表-1 所示。表-1 T700/BA9916 材料属性 参数 值 强度 值 E /GPa 1 E /GPa 2 E /GPa 3 μ 12 μ 13 μ 23 114 X /MPa T X /MPa C Y /MPa T Y /MPa C Z /MPa T Z /MPa C 2688 1458 236 175 G /GPa S /MPa 136 XY G /GPa S /MPa 136 XZ G /GPa S /MPa 23 YZ .2 实例一：UMAT 用户子程序应力分析 在使用 UMAT 用户子程序进行高级应用之前，应该先了解UMAT 子程序，熟悉UMAT 子程序的工作原理，了解UMAT 中的参数、变量的含义。为了便于读者快速了解和使用 UMAT，本例通过复合材料单层板的应力分析来介绍一个简单的UMAT 子程序。 读者可将本例中的单层板替换为层压板，进行对比分析。 复合材料单层板几何尺寸为 15mm×10mm 复合材料单层板几何尺寸为 15mm×10mm×，纤维方向为 45°，单 层板的 3D 实体模型如图-1 所示，X 轴方向为 0°方向，左侧面施加 X 轴向对称边界条件，下侧面施加Y 轴向对称边界条件，垂直于Z 轴且Z=0 的平面施加Z 轴向对称边界条件，右侧面施加100MPa 的拉力。 本例中单位 系统为 mm、MPa。 图-1 单层板边界条件及加载情况 . UMAT 用户子程序 本例使用的UMAT 用户子程序的全部代码如下，字母C 及“!”之后为注释内容。 1 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 2 1 RPL,DDSDDT,DRPLDE,DRPLDT, 3 2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 4 3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 5 4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC) 6 C 7 INCLUDE 8 C 9 CHARACTER*80 CMNAME 10 DIMENSION STRESS(NTENS),STATEV(NSTATV), 11 1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), 12 2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 13 3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3), 14 4 JSTEP(4) 15 16 DIMENSION EG(6), XNU(3,3), STRAND(6), C(6,6), STRESS0(6) 17 C**************************** 18CEG. E1,E2,E3,G12,G13,G2319 18 C EG. E1,E2,E3,G12,G13,G23 19 C XNU. NU12，NU21,NU13，NU31，NU23，NU32 20 C STRAND.....STRAINT AT THE END OF THE INCREMENT 21 C C.....6X6 STIFFNESS MATRIX 22 C STRESS0.....STRESS AT THE BEGINNING OF THE INCREMENT 23 C**************************** 24 C INITIALIZE X