殼單元與實(shí)體單元 UMAT 有哪些差異？
瀏覽：1769
我在 Abaqus 中使用殼單元（如 S4R）調(diào)用 UMAT 時(shí)發(fā)現(xiàn)結(jié)果與實(shí)體單元不同，想確認(rèn)具體差異。
下面是一個(gè)簡(jiǎn)單的線彈性實(shí)體單元umat，如果我想改為用殼單元使用，應(yīng)該如何修改？
非常感謝大佬解答，解決了可以支付寶發(fā)您紅包，
我的微信：13635374241
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
1 RPL,DDSDDT,DRPLDE,DRPLDT,
2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),
1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3),
4 JSTEP(4)
C 線彈性例子
C 計(jì)算楊氏模量E泊松比XNU剪切模量G拉梅常數(shù)Lamda
E = PROPS(1)
Xmu = PROPS(2)
G = E/2/(1+Xmu)
Lamda = Xmu*E/(1+Xmu)/(1-2*Xmu)
C
C 雅克比矩陣初值為0
DDSDDE(i,j) = 0
C
C 填入雅克比矩陣左上角為lamda
DO i=1,NDI
DO j=1,NDI
DDSDDE(i,j)=Lamda
END DO
END DO
C
C 左上角主對(duì)角元素改為L(zhǎng)amda+2G
DO i=1,NDI
DDSDDE(i,i)=Lamda+2*G
END DO
C
C 填入雅克比矩陣右下角
DO i=NDI+1,NTENS
DDSDDE(i,i)=G
END DO
C
C 由應(yīng)變?cè)隽坑?jì)算應(yīng)力
DO i=1,NTENS
DO j=1,NTENS
STRESS(i)=STRESS(i)+DDSDDE(i,j)*DSTRAN(j)
END DO
END DO
RETURN
END
