ABAQUS用户材料子程序开发及应用

华中科技大学硕士学位论文ABAQUS用户材料子程序开发及应用姓名：杨曼娟申请学位级别：硕士专业：桥梁与隧道工程指导教师：卫军20050417IRankineMohr-CoulombABAQUSABAQUSABAQUSUMATUMATUMATFORTRANRankineMohr-CoulombEulerUMATFLAC-3DUMAT()ss=fABAQUSABAQUSABAQUSRankineMohr-CoulombABAQUSUMATEulerIIAbstractTunnel-typeanchorage,akindofgeotechnicalengineeringstructure,isthemainsupportingsystemofasuspensionbridge.Sinceitishardtoperformafull-scalemodeltest,numericalanalysisbecomesanimportantbasisandalsooneofthemostdifficultandpivotalproblemsindesign.Duetothecomplexityofgeotechnicalproblems,asingleconstitutivemodelcan'tsatisfypracticalrequestssometimes.Inthispaper,anewmodelbycombiningRankineandMohr-Coulombcriteriaisadopted,andtheanalysisisdoneatthesimulatingplatformofABAQUS,whichisfamousfiniteelementanalysissoftware.ABAQUSisthemostcapablesoftwareespeciallyforsimulatingnonlinearproblems.Itincludesmanydifferentkindsofmaterialmodelsandfailurecriteria,andprovidesanopeninterfaceandsomeusersubroutinesthatallowuserstoexpandthemainfunctionintheformofcode.However,theconstitutivemodelneededinthisresearchisnotprovided,soitisnecessarytoredevelopitbyuser-definedmaterialsubroutineinterface.EstablishingaconstitutivemodelandchoosinganintegralalgorithmicmethodarethecrucialmattersfordevelopingUMATsubroutine.Inthispaper,basedontheredevelopingtechniqueofUMAT,themodelbycombiningRankineandMohr-Coulombcriteriawasestablished,andcorrespondinginterfaceprogramwasdeveloped.ThealgorithmicmethodforupdatingstressadoptedtheimplicitEulerone,whichisdifferentfromcommonexplicitmethods.TheUMATsubroutineshowscorrectnessandaccuracybyelementtests.Basedonthismodel,athree-dimensional(3D)elasto-plasticanalysisonthetunnel-typeanchorageofSiduhesuspensionbridgewasperformed.TheresultshadgoodconsistencywiththeonefromFLAC-3DanalysisprovidedbyYangtzeRiverScientificResearchInstitute,whichshowsthefeasibilityandvalidityoftheinterfacesubroutine.TheUMATsubroutinedevelopedinthispaperisapplicableforallyieldfunctionsinIIItheformof()ss=fonlybysolving“fiveparameters”mentionedinthispaper.ThisworkexpandstheapplicationrangeofABAQUSinelasto-plasticfiniteelementanalysis,andsetanexamplefordevelopingotherfunctionsofABAQUSandfiniteelementsoftware.Keywords:RankinecriterionMohr-CoulombmodelABAQUSsoftwareUMATsubroutineimplicitEulermethodelasto-plasticanalysis______111.112Mohr-CoulombMohr-Coulomb3RankineMohr-CoulombMohr-CoulombABAQUS2ABAQUSABAQUSABAQUSDrucker-PragerMohr-CoulombABAQUSRankineMohr-CoulombABAQUSUMATRankineMohr-CoulombABAQUS1.2ANSYSMSCABAQUS4ANSYSANSYS53ANSYS6MSCMARCMARC7MSCMARCABAQUSUserSubroutinesABAQUSUser-definedMaterialMechanicalBehaviorUMATABAQUS8UMATGyoo-DongJung9UMATKyunghoonLeeSomnathGhosh10UMATVoronoi11Ducan-ChangUMATABAQUSUMAT41.3ABAQUSUMATUMATFORTRAN(1)ABAQUSUMATABAQUSUMAT(2)RankineMohr-Coulomb(3)UMATRankineMohr-CoulombEulerEuler(4)UMAT(5)FLAC-3DRankineMohr-Coulomb52ABAQUS2.1ABAQUSABAQUS12ABAQUSABAQUS1ABAQUS2ABAQUSABAQUS3ABAQUS4338ABAQUS6Druker-PragerMohr-CoulombABAQUS4ABAQUSABAQUS2.2ABAQUSABAQUSABAQUS/StandardABAQUS/ExplicitABAQUS/StandardABAQUS/AquaABAQUS/USAABAQUS/PreABAQUS/PostABAQUSABAQUS/CAEABAQUS/Design13ABAQUSpre-processingsimulationpost-processing2-1ABAQUS/PreABAQUSABAQUS/StandardABAQUS/Standard7ABAQUS/PostX-Y2-1ABAQUS2.3ABAQUSABAQUSNewton-Raphson132-2PΔABAQUSABAQUS/Prejob.inpABAQUS/StandardABAQUS/Explicitjob.resjob.filABAQUS/Post80KPΔacauaKABAQUSaIPaIaRaaIPR-=aR2-2aaRABAQUSaRABAQUS0.52-2aRABAQUSac0uuuaa-=Δ1ABAQUSABAQUSaKaRABAQUS2-3ABAQUSABAQUS9ABAQUS16ABAQUS25ABAQUS5ABAQUS502-32.42.4.1ABAQUSABAQUSUtilityRoutineABAQUS6.4421510ABAQUSABAQUS2.4.2ABAQUSABAQUSFortranSUBROUTINES(x1,x2,,xn)INCLUDE‘ABA_PARAM.INC’ABAQUS/StandardORINCLUDE‘VABA_PARAM.INC’)ABAQUS/ExplicitRETURNENDx1x2xnABAQUSABAQUSSUBROUTINEDLOADKSTEPKINCCOORDSFaba_param.incvaba_param.incABAQUSABAQUSRETURNEND11.forCAEEDITJOBGENERALUSERSUBROUTINEFILEABAQUS.COMMANDabaqusjob=job-nameuser={source-file|object-file}14(1)FortranABAQUSABAQUSFortranKABAQUS(2)OPEN1518100ABAQUS(3)ABAQUSSTEPSTEPINCREMENTABAQUSSTEPINCREMENTKSTEPKINCABAQUSABAQUS2.5UMATUser-definedMaterialMechanicalBehaviorUMATABAQUSFortranABAQUSUMATABAQUSABAQUS1512*USERMATERIALUMATUMAT1ABAQUS2ABAQUS3UMATJacobianUMATUMATUMATSUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,1RPL,DDSDDT,DRPLDE,DRPLDT,2STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,3NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,4CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)INCLUDE‘ABA_PARAM.INC’CHARACTER*80CMNAMEDIMENSIONSTRESS(NTENS),STATEV(NSTATV),1DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),2STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),3PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)usercodingtodefineDDSDDE,STRESS,STATEV,SSE,SPD,SCDand,ifnecessary,RPL,DDSDDT,DRPLDE,DRPLDT,PNEWDT13RETURNENDUMATDDSDDEDDSDDTDRPLDENDINSHRUMATUMATDDSDDE(NTENSNTENS)NTENSNTENSJacobianes∂∂sΔeΔDDSDDE(ij)jiJacobian*USERMATERIALUNSYMMSTRESS(NTENS)NDINSHRUMATUMATUMATUMATUMATUMATSTATEV(NSTATEV)UMATUSDFLDUEXPANUMATABAQUS*DEPVARPROPS(NPROPS)*USERMATERIALCONSTANTS14USERMATERIALSSESPDSCDSTRAN(NTENS)DSTRAN(NTENS)DTIMENDINSHRNTENSNTENSNDINSH