飞箭培训

北京飞箭软件公司BeijingFEGENsoftwareCo.Ltd.(2ndpart)ZhouYongfaProjectmanagerFEGENSoft飞箭培训1.AnalysisofElastic-PlasticProblemandYieldFunction2.Simulatingmethodofanchor3.UserPractice:UndergroundDigginginMohr-CoulombRockFieldElastic-PlasticProblemFEGENSoft飞箭培训1.Calculationoftheelastic-plasticstressaroundpinholeUseyieldfunctionofMohr-CoulombE-PlasticstressanalysisFEGENSoft飞箭培训AholeinthemiddleoftheelasticplateR0:theradiusoftheholeP:thepressurearoundtheplateR0＝1.0mP＝1600Pac＝1000Paf＝0.2ProblemDescriptionFEGENSoft飞箭培训WithMohr-Coulombyieldfunction,wehavethetheoryresultaroundthepinholeStressinelasticzoneprRr2201prR22010rProblemdescriptionFEGENSoft飞箭培训Stressinplasticzone1sin1sin20Rrctgcr1sin1sin1sin1sin20RrctgcProblemdescriptionFEGENSoft飞箭培训Radiusofplasticzonesin2sin10sin1ctgcctgcpRRPThetheoryradiusofplasticzoneis0.1293m。InputR0＝1.0m，p＝1600Pa，c＝1000Pa，f＝0.2toaboveexpressTheoreticalsolutionFEGENSoft飞箭培训cos2sin3131cFMohr－CoulombyieldfunctionMohr-CoulombFEGENSoft飞箭培训1.Usingsolidmechanicsformulalibrarygenerate2Delastic-plasticprogram2.ModifythePREparameter3.AddMohr-Coulombyieldfunction4.Modelingandmeshing5.CalculatingandanalysisSolvingstepFEGENSoft飞箭培训FEPGAppwizardSolidE_Plastic2DXYAWZsdispeq4,ll2,ilu,outcoreOKUsingformulalibFEGENSoft飞箭培训1.Usingsolidmechanicsformulalibrarygenerate2Delastic-plasticprogram2.ModifythePREparameter3.AddMohr-Coulombyieldfunction4.Modelingandmeshing5.CalculatingandanalysisSolvingstepFEGENSoft飞箭培训ModifyPREfilematedatanespaaeq41.0e10;0.3;0;0;0.2;1.0e3;1.0;3000;0.6;nespaall20;100;nespbbeq41.0e10;0.3;0;0;0.2;1.0e3;1.0;3000;0.6;matedatanespaaeq41.0e10;0.3;0;0;0.2;1.0e3;0.0;nespaall20;1600;nespbbeq41.0e10;0.3;0;0;0.2;1.0e3;0.0;FEGENSoft飞箭培训1.Usingsolidmechanicsformulalibrarygenerate2Delastic-plasticprogram2.ModifythePREparameter3.AddMohr-Coulombyieldfunction4.Modelingandmeshing5.CalculatingandanalysisSolvingstepFEGENSoft飞箭培训defidispuvcoorxycoefun,vnfuncexxeyyexyshap%1%2gaus%3$c6dimensionde(3),e(2,2),d(3,3),dv(3),dg(2,2)$c6dimensiondp(3,3),p(4)$c6externalpragerDrucker-Pragerfunction$c6logicalfilestr$c6common/gpstr/gpstr(6,100000),str(6)nespa.pdeFEGENSoft飞箭培训……$cvde(1)=e(1,1)$cvde(2)=e(2,2)$c6de(3)=e(1,2)+e(2,1)$c6callgetstr(p,str,de,dv,d,dp,1,prag,prager)$c6if(prag.gt.-p(2)*0.0001)then$c6doi=1,3$c6doj=1,3$c6d(i,j)=d(i,j)-dp(i,j)……nespa.pdeDrucker-PragerfunctionFEGENSoft飞箭培训Drucker-Prager0)(212121HJIJIF，sinsin233＝cossincH633＝FEGENSoft飞箭培训PlasticmatrixCalculationStres.for&plas.forCalculateplasticmatrixDpCalculateplasticstressCodesofyieldfunctionliesinplas.forFEGENSoft飞箭培训Plas.forreal*8FUNCTIONPRAGER(X,U)implicitreal*8(a-h,o-z)DIMENSIONU(4),S(3),x(4)tgo=x(1)q0=x(4)c=x(2)+q0*U(4)PV=X(3)D=(U(1)+U(2))*(1+pv)/3.S(1)=U(1)-DS(2)=U(2)-DS(3)=U(3)SourcecodesofDrucker-PragerfunctionInplas.forDrucker-PragerfunctionFEGENSoft飞箭培训nespa.pdeUK=(u(1)+u(2))*pv-DSJ=(S(1)**2+S(2)**2+UK**2)/2.+S(3)**2alpha=tgo/dsqrt(9+12*tgo**2)ck=3*c/dsqrt(9+12*tgo**2)PRAGER=DSQRT(SJ)+ALPHA*D*3-CKRETURNENDSourcecodesofDrucker-PragerfunctionInplas.forFEGENSoft飞箭培训AddM-Cfunctioninplas.forcos2sin3131cFMohr－CoulombyieldfunctionFEGENSoft飞箭培训AddM-Cfunctioninplas.forreal*8FUNCTIONMohrC(X,U)implicitreal*8(a-h,o-z)DIMENSIONU(4),X(4),EMSTR(3)emstr(1)=U(1)+U(2)coefb=-U(1)-U(2)coefc=U(1)*U(2)-U(3)**2emstr(2)=((-coefb)+dsqrt(coefb**2-4*coefc))/2emstr(3)=((-coefb)-dsqrt(coefb**2-4*coefc))/2if(emstr(1).lt.emstr(2))thenetemp=emstr(2)emstr(2)=emstr(1)emstr(1)=etempendifSourcecodesofMohr-CoulombfunctionInplas.forMohr-CoulombfunctionFEGENSoft飞箭培训AddM-Cfunctioninplas.forif(emstr(1).lt.emstr(3))thenetemp=emstr(3)emstr(3)=emstr(1)emstr(1)=etempendifif(emstr(2).lt.emstr(3))thenetemp=emstr(2)emstr(2)=emstr(3)emstr(3)=etempendiftgo=x(1)c=x(2)cosf=1./dsqrt(1.+tgo*tgo)SourcecodesofMohr-CoulombfunctionInplas.forFEGENSoft飞箭培训AddM-Cfunctioninplas.forMohrC=emstr(1)-emstr(3)-2.0*c*cosfMohrC=MohrC+(emstr(1)+emstr(3))*(dsqrt(1-cosf*cosf))RETURNENDcos2sin3131cFSourcecodesofMohr-CoulombfunctionInplas.forFEGENSoft飞箭培训Compileplas.forfl32.exe/W0/nologo/cplas.FORRightclickmousebuttononplas.forinFEPGinterfaceThenrunFEGENSoft飞箭培训defidispuvcoorxycoefun,vnfuncexxeyyexyshapq4gausq$c6dimensionde(3),e(2,2),d(3,3),dv(3),dg(2,2)$c6dimensiondp(3,3),p(4)$c6externalMohrC$c6logicalfilestr$c6common/gpstr/gpstr(6,100000),str(6)Modifynespa.pdeModifypragertoMohrCFEGENSoft飞箭培训……$cvde(1)=e(1,1)$cvde(2)=e(2,2)$c6de(3)=e(1,2)+e(2,1)$c6callgetstr(p,str,de,dv,d,dp,1,prag,MohrC)$c6if(prag.gt.-p(2)*0.0001)then$c6doi=1,3$c6doj=1,3$c6d(i,j)=d(i,j)-dp(i,j)……Modifynespa.pdeModifypragertoMohrCFEGENSoft飞箭培训Compilenespa.pdeClickruncommandbuttoninFEPGinterfaceFEGENSoft飞箭培训Compilenespa.pdeIncommandlinerunpdenespaaeq4FEGENSoft飞箭培训Compilenespa.pdeLet.batnespaRightclickmousebuttononnespainFEPGinterfaceThenrunFEGENSoft飞箭培训dispsxx,syy,sxy,wp,vep,vyp,didcoorxycoefun,vnshapq4gausqmassq