第1页 / 共21页
第2页 / 共21页
第3页 / 共21页
第4页 / 共21页
第5页 / 共21页
第6页 / 共21页
第7页 / 共21页
第8页 / 共21页
abaqus二次开发.pdf
Power Law UMAT
December
Mechanics
No summmation in the following expressions
ij Sij p ij
p k tr
For shear components introduce engineering strain
ij ij
Sii
m
and Sij
m
ii
ij
ij
ii X
ij ij
X
X
ij
X
ii
ii
jj
ij
and
ij
kl
ikjl
I ijkl
ikjl iljk
J
ijt t
klt t
ijt t
klt t
tr
kl
kl
kk
kk
kl
kl
Now we can compute S
ii
kk
P
ij
P
kl
kk
kl
Sii
kk
Sij
kk
Sii
kl
Sij
kl
m
m
m
m
kk
ii
ii
kk
kkii ik
m
kk
ij
ij
kk
m
kk ij
m
m
m
m
m
kl
ii
ii
kl
m
klii
m
m
m
m
kl
ij
ij
kl
klij ikjl
p
kk
ij k
tr
kk
ij kij
and
p
kl
ij k
tr
kl
ij
Now we can compute Jacobian
Sii
kk
p
kk
m
m
kkii ik k
J iikk
J ijkk
J iikl
J ijkl
ii
kk
ij
kk
ii
kl
ij
kl
Sij
kk
m
p
kl
m
m
Sii
kl
m
Sij
kl
kk ij
klii
m
m
klij ikjl
and Stress
ii Sii p
ii k tr
m
m
ij
ij S ij
Coding
P N LT k
T ERM
m
T ERM
T ERM
For Jacobian
m
T ERM T ERM
J iikk T ERM T ERM kkii ik P N LT
J iikl J klii T ERM ii kl
J ijkl
T ERM
T ERM ij kl ikjl
For Stress
ii T ERM ii P N LT tr
ij
T ERM ij
ABAQUS input le Uniaxial Tension
HEADING
UMAT POWER LAW INCOMPRESSIBLE MATERIAL CD UMATPLT
WAVEFRONT MINIMIZATIONSUPPRESS
NODENSETALLN
ELEMENTTYPECDELSETALLE
SOLID SECTIONELSETALLEMATERIALALLE
MATERIALNAMEALLE
USER MATERIALCONSTANTS
E v POWER sig eps StTol Pnlt
EEE
USER SUBROUTINE
SUBROUTINE UMATSTRESSSTATEVDDSDDESSESPDSCD
RPLDDSDDTDRPLDEDRPLDTSTRANDSTRAN
TIMEDTIMETEMPDTEMPPREDEFDPREDMATERLNDINSHRNTENS
NSTATVPROPSNPROPSCOORDSDROTPNEWDTCELENT
DFGRD DFGRDNOELNPTKSLAYKSPTKSTEPKINC
INCLUDE ABAPARAMINC
CHARACTER CMNAME
DIMENSION STRESSNTENSSTATEVNSTATV
DDSDDENTENSNTENSDDSDDTNTENSDRPLDENTENS
STRANNTENSDSTRANNTENSTIMEPREDEFDPRED
PROPSNPROPSCOORDSDROT
DFGRD DFGRD
DIMENSION STRANTDELTA
C
C
C
C
PARAMETER ONE D TWO D THREE D SIX D
HALFD ZERO D
DATA NEWTONTOLER D
KRONECKERS DELTA
C
C
C
DATA DELTA D D D
D D D
D D D
UMAT FOR ISOTROPIC ELASTICITY AND ISOTROPIC PLASTICITY
J FLOW THEORY
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
PROPS E
PROPS NU
PROPS POWER POWER LAW EXPONENT
PROPS SIG
PROPS EPS
PROPS STTOL ELASTIC STRAINYIELD STRAIN
PROPS PNLT PENALTY FOR INCOMPRESSIBILTY
COMPUTE TOTAL STRAIN
DO KNTENS
STRANTK STRANKDSTRANK
ENDDO
C
C
C
COMPUTE MEAN STRAIN
STNMN ZERO
DO K NDI
STNMN STNMN STRANTKSTRANTK
ENDDO
DO K NDINTENS
STNMN STNMN HALFSTRANTKSTRANTK
ENDDO
C
C
C
C
C
C
C
C
C
C
C
C
C
STNMN TWOTHREESTNMNHALF
ELASTIC PROPERTIES
EMODPROPS
ENUPROPS
IFENUGT ANDENULT ENU
EBULKEMODONETWOENU
EGEMODONEENU
EGEGTWO
EGTHREEEG
ELAMEBULKEGTHREE
MATERIAL PROPERTIES AND PENALTIES
POWER PROPS
SIG PROPS
EPS PROPS
STTOL PROPSEPS
PNLT PROPSEMOD
SWITCH FOR LINEARPOWER LAW BEHAVIOR
IF STNMNLESTTOL THEN
ELASTIC STIFFNESS
DO KNTENS
DO KNTENS
DDSDDEKKZERO
CONTINUE
CONTINUE
DO KNDI
DO KNDI
DDSDDEKKELAM
CONTINUE
DDSDDEKKEGELAM
CONTINUE
DO KNDINTENS
DDSDDEKKEG
CONTINUE
C
C
C
CALCULATE STRESS FROM ELASTIC STRAINS
DO KNTENS
DO KNTENS
STRESSKSTRESSKDDSDDEKKDSTRANK
CONTINUE
CONTINUE
C
C
C
C
C
C
C
C
C
C
C
C
ELSE
NONLINEAR BEHAVIOR
PRESSURE
VOLDEF
DO K NDI
VOLDEF VOLDEF STRANTK
ENDDO
PRESS PNLTVOLDEF
PRECOMPUTED TERMS
TERM TWOTHREESTNMNEPS POWERONESIG EPS
TERM POWERONETWOTHREESTNMNSTNMN
TERM HALFTERMTERM
STRESS
DO K NDI
STRESSK TERMSTRANTK PRESS
ENDDO
DO K NDI NTENS
STRESSK HALFTERMSTRANTK
ENDDO
相关推荐
2023年江西萍乡中考道德与法治真题及答案.doc
2012年重庆南川中考生物真题及答案.doc
2013年江西师范大学地理学综合及文艺理论基础考研真题.doc
2020年四川甘孜小升初语文真题及答案I卷.doc
2020年注册岩土工程师专业基础考试真题及答案.doc
2023-2024学年福建省厦门市九年级上学期数学月考试题及答案.doc
2021-2022学年辽宁省沈阳市大东区九年级上学期语文期末试题及答案.doc
2022-2023学年北京东城区初三第一学期物理期末试卷及答案.doc
2018上半年江西教师资格初中地理学科知识与教学能力真题及答案.doc
2012年河北国家公务员申论考试真题及答案-省级.doc
2020-2021学年江苏省扬州市江都区邵樊片九年级上学期数学第一次质量检测试题及答案.doc
2022下半年黑龙江教师资格证中学综合素质真题及答案.doc
