How to make UMATs for soils in Abaqus - Uniandes

81 Pages·2008·1.9 MB·English

HowtomakeUMATsforsoilsinABAQUS
How to make UMATs for soils in ABAQUS
GeotechnicalResearchGroup
UniversityofLosAndes
http://geotecnia.uniandes.edu.co
December21,2008
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
Introduction
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
Example:
MODIFIEDCAM-CLAYPlasticityconstitutiveequation:
T=Cep :ε (1)
withCep =

Ce, γ =0

1 (cid:104) (cid:104)−→ −→ (cid:105) −→ −→ (cid:105)
Cep =Ce− a 1⊗1+a T∗⊗1+1⊗T∗ +a T∗⊗T∗ , γ >1
 χ 1 2 3
wherea ,a ,a andχarefunctionofTandmaterialconstant,and
1 2 3
−→ T− 11
T∗ = 3 (2)
(cid:107)T− 11(cid:107)
3
YoumakeanUMATandbuildaFEMmodelinABAQUS!
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
FEM models
Step: Puestali Frame: 22
POR
(Avg: 75%)
−+++++++32468111........113570356875491356396321eeeeeeee−+++++++0000000071111222 (BSAoDvtVgto:2 m−−−−++++7321461245, ........%(97668891fr27293384)a85243692ceeeeeeeeti−−−−−−−−o00000000n33344333 = −1.0)
Simulación suelo−tuberia (c) Porepressureu SOiDmBu:l atecisóisn9 sBu.eoldob− t u bAebraiaqus/S(tadnda)rd VSersiohn 6e.7−a1 r Sats Detcr 1e3 11s:14s:08 TSA Pxacyific Standard Time 2008
Y SIOntDcerBpe::m Pteeusneitss 9t a B l2i.mo2d:o bSc ot e m pAp bTaaicmqtaueds =o/S s t ua 9pn0ed0raf.ir0cdia Vlersion 6.7−1 Sat Dec 13 11:14:08 SA Pacific StanFdardi gTimue 2r00e8:YZCXAM-SIPDntreCceimfrpoea:rm mrPyeLeu nVdeta s VA tr a: a lS2rim:Y2D U:oV Sc 2 otDemmpep fToaricmmtaoead t=ioo sd n u 9Sp0ece0raf.li0eci laFlactor: +1.000e+01
Z X PDreimfoarmrye Vda Vr:a Pr:O UR Deformation Scale Factor: +1.000e+01
1
1UMATfromFUENTES
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
FEMmodels
F(SA,E vSg2: −−−−−−−−−−−−−2711198765432125.............%21011111111211114556789905M)2332086421970eeeeeeeeeeeee+++++++++++++00000000000002221111111110 models (SA, vSg2: −−−−−−−−−−−−−2711198765432115.............%21011111111191117778889999)6671470470378eeeeeeeeeeeee+++++++++++++00000000000002221111111110
ODB: tabla.odb Abaqus/Standard Version 6.7−1 Sun Dec 21 15:36:51 SA Pacific Standard Time 2008
YZ2X34SIPDOntreDceimfrBpoea:r:m m rStyaeet benVdlptaa V− .r o3: a d,S1r :b 5, U A:S Su 2 A st2Debhpeau fqbTou,ri mmsS/eoaS htt=ialoe n n df r1Sae.cr0ida0 (lV0e eFarasciot)onr 6: .+71−.C10 0 0 Se+uA0n0 DeMc 21 -15C:42:15L SAA PacYific Stamndardo Timde YZ20e08l,XTySIPDntreceimfrpoea:rm mrSyeet enVdpta V− r 3: a ,S1r : 5, UA:S Su2 st2Dehpeu fbTo,ri mmSeoa ht=iloe n f r1Se.c0ia0l0e Factor: +1.0(00be+)00Hypoelasticmodel,Ty
2ModelfromFELLINinhttp://www.uibk.ac.at/geotechnik/res/hypopl.html
3UMATsfromFuentes.
4InHypoelasticity,theRICHART’sproposalistaken.
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
FEMmodels
FEM models
(SA, vSg2: −−−−−−−−−−−−+2711197654321815.............%21009999888290000853196429)7534186318579eeeeeeeeeeeee+++++++++++++00000000000002221111111100 (SA, vSg2: −−−−−−−−−−−−+2711197654321815.............%21009999888290000853196429)7533186308579eeeeeeeeeeeee+++++++++++++00000000000002221111111100
ODB: tabla.odb Abaqus/Standard Version 6.7−1 Wed Dec 10 17:34:46 SA Pacific Standard Time 2008 ODB: tabla_Fellin_Euler.odb Abaqus/Standard Version 6.7−1 Sun Dec 07 10:17:28 SA Pacific Standard Time 2008
YZ X(cSIPDntreceimfrpoea:rm m)rSyeet enVdpta V− r 3: a ,S1r : 5, UA:S Su2 st2DHehpeu fbTo,ri mmSeoa hty=iloe n f r1Se.cp0ia0l0e Faoctor: +p1.000el+0a0sticity model FUENTEYZSX’sSIPDntreceimfrpoea:rm mrSyeet enVdpta (V− r 3: a ,S9rd : 6, UA:S Su2 st2De)hpeu fbTo,ri mmSeoa ht=iloe n f r1Se.cH0ia0l0e Facytor: +p1.000oe+00plasticity model FELLIN’s
UMAT,Ty UMAT,Ty
5
5WOLFFERSDORFFequation
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
FEMmodels
FEM models
Step: Step−2 Frame: 505
(PAOvRg: 75%)
(SA, vSg2: −−−−−−−−−−−−−−27111118765421555..............%4321087531974486420901123496)66789531864112eeeeeeeeeeeeee++++++++++++++0000000000000022222111111100 −++++++++++++2611233455677.............653962952851870050505050500000000000000eeeeeeeeeeeee−++++++++++++00000000000006122222222222 Step: Step−2 Frame: 277
COODNB:S cOoLnIsDoAlidTaIOciNonHP.odb Abaqus/Standard Version 6.7−1 Fri Dec 12 14:37:38 SA Pacific Standard Time 2008
ZYX SIPDntreceimfrpoea:rm mrSyeet enVdpta V− r 2: a 2,Sr 7: , C7 US:o 2 Sn 2Dtseoeplfi odTraimmtioaent i=o n S2c7a.2le7 Factor: +1.000e+00
(e) Ty(kPa) COODNB:S cOoLnIsDoAlidTaIOciNonHP.odb Abaqus/Standard Version 6.7−1 Fri De(cf 1)2 14:37:38 SA Pacificu Standard Time 2008
YZ X SIPDntreceimfrpoea:rm mrSyeet enVdpta V− r 2: a 5,Pr 0: O C5U:Ro Sn Dtseoeplfi odTraimmtioaent i=o n S5c0a.0le0 Factor: +2.000e+02 (kPa)
Figure:Consolidationin3Delements,HypoplasticityFUENTES’sUMAT
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
INCREMENTALDRIVER,from(NIEMUNIS,2007)
INCREMENTAL DRIVER, from (NIEMUNIS, 2007)
600
pc=100kPa
qDeviatorstress[kPa]240000 =====1235150000000000kkkkkPPPPPaaaapa0=100kPa Volumetricstrain[-]εv−−−0000....000009633 pc=====123515000000000kkkkkPPPPPaaaaa
=1000kPa
0 0.12
0 0.03 0.06 0.09 0.12 0.15 − 0 0.03 0.06 0.09 0.12 0.15
Verticalstrainε1[-] Verticalstrainε1[-]
Figure:Drainedtriaxialpath,variationofpreconsolidatedpressurepc’,CAM-CLAY
model
6
6ForINCREMENTALDRIVERsee[Niemunis,2008]
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
INCREMENTALDRIVER,from(NIEMUNIS,2007)
INCREMENTAL DRIVER, from (NIEMUNIS, 2007)
2000 pc=200kPa 2000 M=1.2
Pa] ==12000000kkPPaa Pa]
[k1600 =3000kPa [k1600
q =4000kPa q
stress1200 p0=200kPa stress1200
Deviator 480000 Deviator 480000
0 0
0 0.01 0.02 0.03 0.04 0.05 0 400 800 1200
Verticalstrainε1[-] Meanstressp[kPa]
Figure:Undrainedtriaxialpath,variationofpreconsolidatedpressurepc,CAM-CLAY
model
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS
HowtomakeUMATsforsoilsinABAQUS
Introduction
INCREMENTALDRIVER,from(NIEMUNIS,2007)
Some errors in programming
45 45
40 40
Pa]35 Pa]35
k k
T[1230 T[1230
stress2205 stress2205
Shear1105 p0=100kPa Shear1105 p0=100kPa
5 5
Simulated Simulated
0 0
0 0.02 0.04 0.06 0.08 0.1 0 0.02 0.04 0.06 0.08 0.1
Shearstrainε12[-] Shearstrainε12[-]
(a) FELLIN’sUMAT (b) FUENTES’sUMAT
Figure:ProportionalpathinshearstresscontrolinHypoplasticity
GeotechnicalResearchGroupUniversityofLosAndeshttp://geotecnia.unHioanwdetos.medauk.ecUoMATsforsoilsinABAQUS