哪位好心人能给一个abaqus子程序应用的实例

逍遥生1314

求一个fortran编写的子程序代码和对应的cae模型，老师叫做本构的二次开发，从百度上找的例子根本算不了，提示子程序错误就30处（pis：完全复制的人家的程序，验证就有错误）

ddongxupeng

1. SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
   
2.      1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
   
3.      2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,
   
4.      3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
   
5.        include 'aba_param.inc'
   
6.        CHARACTER*80 CMNAME
   
7.        DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
   
8.      1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
   
9.      2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
   
10.      3 DFGRD0(3,3),DFGRD1(3,3)
    
11.       PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0)
    
12. C UMAT FOR ISOTROPIC ELASTICITY
    
13. C CANNOT BE USED FOR PLANE STRESS
    
14. C ----------------------------------------------------------------
    
15. C PROPS(1) - E
    
16. C PROPS(2) - NU
    
17. C ----------------------------------------------------------------
    
18. C
    
19.       IF (NDI.NE.3) THEN
    
20.       WRITE (6,*) 'THIS UMAT MAY ONLY BE USED FOR ELEMENTS
    
21.      1 WITH THREE DIRECT STRESS COMPONENTS'
    
22.       CALL XIT
    
23.       ENDIF
    
24.       EMOD=PROPS(1)
    
25.       ENU=PROPS(2)
    
26.       EBULK3=EMOD/(1-2*ENU)
    
27.       EG2=EMOD/(1+ENU)
    
28.       EG=EG2/2
    
29.       EG3=3*EG
    
30.       ELAM=(EBULK3-EG2)/3
    
31. C
    
32. C ELASTIC STIFFNESS
    
33. C
    
34.       DO K1=1, NDI
    
35.       DO K2=1, NDI
    
36.       DDSDDE(K2, K1)=ELAM
    
37.       END DO
    
38.       DDSDDE(K1, K1)=EG2+ELAM
    
39.       END DO
    
40.       DO K1=NDI+1, NTENS
    
41.       DDSDDE(K1 ,K1)=EG
    
42.       END DO
    
43. C
    
44. C CALCULATE STRESS
    
45. C
    
46.       DO K1=1, NTENS
    
47.       DO K2=1, NTENS
    
48.       STRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1)
    
49.       END DO
    
50.       END DO
    
51.       RETURN
    
52.        END

复制代码

逍遥生1314

:Q:Q什么都不会，有没有来指点下

ddongxupeng

模型代码

icanlax

ddongxupeng 发表于 2014-12-18 20:46

谢谢层主，程序没有问题，运行成功

平行线不相交

:):):):):):)

164177231

楼主，还在不？能不能说一下，这个程序带入到那里试呢？还请楼主不吝赐教啊

断水流

ddongxupeng 发表于 2014-12-18 20:46

谢谢，完全好用

tml2012

ddongxupeng 发表于 2014-12-18 20:46

不知能不能发点弹塑性本构的子程序，十分感谢

江南2016

楼主有做纤维模型的子程序没，我下了个纤维模型的子程序，但是不懂怎么建了，有小实例可以共享下吗

dengbo110

学习了！

liyushu

谢谢分享