ABAQUS初学者用户子程序小例子.pdfVIP

实例：最简单的杆件单轴拉伸，材料本构模型采用线弹性模型， E ＝206000MPa, ν=0.3 。左端约束，右端施加均布力q=1 。 Author: xueweek 希望这个小例子能帮助初学者减少一两天的摸索时间 1 材料本构模型：σ=Eε 在 Property 中定义材料时，在 General 下选中 User Material ，输入 206000 ，0.3 。 以上两个数值代表E 、 ν，在用户子程序中代表着PROPS(1) 、PROPS(2) 。然后在General 下选中 Depvar ，由于该例子中没有使用状态变量，因此在第一项中输入大于 0 的数值即 可。 建模大家都会，故省略 2 ABAQUS 调用用户子程序 大家在论坛上也可以看到很多方法，我这里指介绍最简单的一种。 在 Job Manager 中点击Edit 选项，在 General 选项的最后一项中选择自己建立好的用户 Fortran 子程序文件。（注：用户子程序文件可以使用文本编辑器进行编辑，当然也可以用 编译器，如果对自己的用户子程序文件的语法不放心，可以先用 Fortan 编译器进行编译， 不过编译前要先建立 Project ，关于Fortran 编译，这里不再介绍）。 完成后，点击 submit 即可进行分析。 3 结果 以下两张图分布是用户子程序和 ABAQUS 自带的材料模型得到的应力 云图，可以看出两种图形完全相同。 4 UMAT 子程序 对于初学者来说，需要注意的是，FORTRAN 对于程序语言格式上的要求。例如，对于 FORTRAN 语言，前六个字符必须空出来，等等。检查语法最好的方法就是在 FORTRAN 编 译器上进行编译。 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT, 1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED, 2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT, 3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) include aba_param.inc CHARACTER*8 CMNAME DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS), 1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS), 2 TIME(2),PREDEF(1),DPRED( 1),PROPS(NPROPS),COORDS(3),DROT(3,3), 3 DFGRD0(3,3),DFGRD1(3,3) C UMAT FOR ISOTROPIC ELASTICITY C CANNOT BE USED FOR PLANE STRESS C C PROPS(1) - E C PROPS(2) - NU C C IF (NDI.NE.3) THEN WRITE (6,*) THIS UMAT MAY ONLY BE USED FOR ELEMENTS 1 WITH THREE DIRECT STRESS COMPONENTS CALL XIT ENDIF C C ELASTIC PROPERTIES EMOD=PROPS(1) ENU=PROPS(2) EBULK3=EMOD/(1-2*ENU) EG2=EMOD/(1+ENU) EG=EG2/2 EG3=3*EG ELAM=(EBULK3-EG2)/3 C C ELASTIC STIFFNESS C DO K1=1, NDI DO K2=1, NDI