两种材料,USDFLD子程序中,区别在于两种材料的强度参数不一样,其他都一样,如何编辑让他们放在一起,计算机能够识别?

浏览:83967 回答:3
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,TIME,DTIME,
1 CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,
2 NDI,nshr,coord,jmac,jmtyp,matlayo,laccflg)
C
INCLUDE 'ABA_PARAM.INC'
C
C STRENGTH PARAMETERS
PARAMETER(Xc=600,Xt=700,Yc=120,Yt=70,
$ S12=45)
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),T(3,3),TIME(2)
DIMENSION ARRAY(15),JARRAY(15),coord(*),jmac(*),jmtyp(*)
C
C CALCULATE COEFFICIENT OF TSAI-WU
F1 = 1/Xt-1/Xc
F2 = 1/Yt-1/Yc
F11 = 1/(Xt*Xc)
F22 = 1/(Yt*Yc)
F66 = 1/S12**2
F12 = -0.5/(Xt*Xc*Yt*Yc)**0.5
C
C INITIALIZE FAILURE FLAGS FROM STATEV.
CR11 = STATEV(1)
CR22 = STATEV(2)
CR12 = STATEV(3)
C
C GET STRESSES FROM PREVIOUS INCREMENT
CALL GETVRM('S',ARRAY,JARRAY,FLGRAY,jrcd,
$ jmac, jmtyp, matlayo, laccflg)
ST11 = ARRAY(1)
ST22 = ARRAY(2)
ST33 = ARRAY(3)
ST12 = ARRAY(4)
ST13 = ARRAY(5)
ST23 = ARRAY(6)
C
C TSAI-WU CRITERION
CR = F1*ST11 + F2*ST22 + F11*ST11**2 + F22*ST22**2 +
$ 2*F12*ST11*ST22 + F66*ST12**2
C
C FAILURE TYPE
IF (CR>=1) THEN
H1 = F1*ST11 + F11*ST11**2
H2 = F2*ST22 + F22*ST22**2
H6 = F66*ST12**2
IF (CR11<1) THEN
IF ((H1 .GE. H2).AND.(H1.GE.H6)) THEN
CR11=1
ELSE
CR11=0
ENDIF
ENDIF
IF (CR22<1) THEN
IF ((H2 .GE. H1).AND.(H2.GE.H6)) THEN
CR22=1
ELSE
CR22=0
ENDIF
ENDIF
IF (CR12<1) THEN
IF ((H6 .GE. H1).AND.(H6.GE.H2)) THEN
CR12=1
ELSE
CR12=0
ENDIF
ENDIF
ENDIF
C
C UPDATE STATEV
STATEV(1) = CR11
STATEV(2) = CR22
STATEV(3) = CR12
C
C UPDATE FIELD VARIABLES
IF (CR11>0) FIELD(1) = 1
IF (CR22>0) FIELD(2) = 1
IF (CR12>0) FIELD(3) = 1
RETURN
END
以上就是USDFLD子程序,我在poperty项,定义两种材料,即两种User Defined Field,这两种材料的不同之处就是把上边程序中Xc=600,Xt=700,Yc=120,Yt=70,S12=45 数据修改下,两种材料两种数据,其他都一样一样,在计算时,如何把他们联系起来,计算机识别,求操作,用M1和M2 分别定义两种材料名..不胜感激!
邀请回答 我来回答

全部回答

(2)
默认 最新
oneturtle

子程序里面有个变量CMNAME,这个代表相应材料的名字,也就是你定义的材料名字。子程序里面添加if判断CMNAME对应哪个,不同材料对应不同参数取值。

2021年4月20日
评论 点赞
whbhust
楼主这个问题解决了吗?
2015年4月23日
评论 1 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP