资源描述：
《ABAQUS用户子程序小例子.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实例：最简单的杆件单轴拉伸，材料本构模型采用线弹性模型，E＝206000MPa,ν=0.3。左端约束，右端施加均布力q=1。Author:xueweek希望这个小例子能帮助初学者减少一两天的摸索时间1材料本构模型：σ=Eε在Property中定义材料时，在General下选中UserMaterial，输入206000，0.3。以上两个数值代表E、ν，在用户子程序中代表着PROPS(1)、PROPS(2)。然后在General下选中Depvar，由于该例子中没有使用状态变量，因此在第一项中输入大于0的数值即可。建模大家都会，故省略2ABAQUS调用用户子程序大家在论坛上也可以看到很
2、多方法，我这里指介绍最简单的一种。在JobManager中点击Edit选项，在General选项的最后一项中选择自己建立好的用户子程序文件。（注：用户子程序文件可以使用文本编辑器进行编辑，当然也可以用Fortran编译器，如果对自己的用户子程序文件的语法不放心，可以先用Fortan编译器进行编译，不过编译前要先建立Project，关于Fortran编译，这里不再介绍）。完成后，点击submit即可进行分析。3结果以下两张图分布是用户子程序和ABAQUS自带的材料模型得到的应力云图，可以看出两种图形完全相同。4UMAT子程序对于初学者来说，需要注意的是，FORTRAN对于程序语言格
3、式上的要求。例如，对于FORTRAN语言，前六个字符必须空出来，等等。检查语法最好的方法就是在FORTRAN编译器上进行编译。SUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,1DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,2CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,3PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTE
4、P,KINC)include'aba_param.inc'CHARACTER*8CMNAMEDIMENSIONSTRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),1DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),2TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),3DFGRD0(3,3),DFGRD1(3,3)CUMATFORISOTROPICELASTICITYCCANNOTBEUSEDFORPLAN
5、ESTRESSC----------------------------------------------------------------CPROPS(1)-ECPROPS(2)-NUC----------------------------------------------------------------CIF(NDI.NE.3)THENWRITE(6,*)'THISUMATMAYONLYBEUSEDFORELEMENTS1WITHTHREEDIRECTSTRESSCOMPONENTS'CALLXITENDIFCCELASTICPROPERTIESEMOD=PROP
6、S(1)ENU=PROPS(2)EBULK3=EMOD/(1-2*ENU)EG2=EMOD/(1+ENU)EG=EG2/2EG3=3*EGELAM=(EBULK3-EG2)/3CCELASTICSTIFFNESSCDOK1=1,NDIDOK2=1,NDIDDSDDE(K2,K1)=ELAMENDDODDSDDE(K1,K1)=EG2+ELAMENDDODOK1=NDI+1,NTENSDDSDDE(K1,K1)=EGENDDOCCCALCULATESTRESSCDOK1=1,NTENSDOK2=1,NTENSSTRESS(K2)=STRESS(K2)+DDSDDE(K2,K1)*D
7、STRAN(K1)ENDDOENDDOCRETURNEND5INP输入文件*Heading**Jobname:Job-1Modelname:Model-1**Generatedby:Abaqus/CAE6.10-1*Preprint,echo=NO,model=NO,history=NO,contact=NO****PARTS***Part,name=Part-1*EndPart******ASSEMBLY***Assembly,name=Assembly***Instance,
