壳单元与实体单元 UMAT 有哪些差异？
浏览：655
我在 Abaqus 中使用壳单元（如 S4R）调用 UMAT 时发现结果与实体单元不同，想确认具体差异。
下面是一个简单的线弹性实体单元umat，如果我想改为用壳单元使用，应该如何修改？
非常感谢大佬解答，解决了可以支付宝发您红包，
我的微信：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 计算杨氏模量E泊松比XNU剪切模量G拉梅常数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 左上角主对角元素改为Lamda+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 由应变增量计算应力
DO i=1,NTENS
DO j=1,NTENS
STRESS(i)=STRESS(i)+DDSDDE(i,j)*DSTRAN(j)
END DO
END DO
RETURN
END
![[免费案例]Ensight案例教程分享](https://img.jishulink.com/cimage/245b3ca9e2c939e40491a25edae94515.jpeg?image_process=resize,fw_576,fh_320,)
