飞箭培训ppt课件

1、Simulate Geo-technical Problem(2nd part) Zhou Yongfa Project manager1. Analysis of Elastic-Plastic Problem and Yield Function2. Simulating method of anchor3. User Practice: Underground Digging in Mohr-Coulomb Rock Field Elastic-Plastic ProblemCalculation of the elastic-plastic stress around pinhole

2、Use yield function of Mohr-CoulombE-Plastic stress analysisA hole in the middle of the elastic plateR0: the radius of the hole P: the pressure around the plateR01.0 mP 1600 Pac 1000 Paf 0.2Problem DescriptionWith Mohr-Coulomb yield function, we have the theory result around the pinholeStress in elas

3、tic zone Problem descriptionStress in plastic zone Problem descriptionRadius of plastic zone The theory radius of plastic zone is 0.1293 m。 Input R01.0 m，p1600 Pa，c1000 Pa，f0.2 to above expressTheoretical solutionMohrCoulomb yield function Mohr-Coulomb1.Using solid mechanics formula library generate

4、 2D elastic-plastic program2.Modify the PRE parameter3.Add Mohr-Coulomb yield function4.Modeling and meshing5.Calculating and analysisSolving step FEPG Appwizard Solid E_Plastic 2DXY AWZ sdispe q4, ll2, ilu, outcore OKUsing formula lib1.Using solid mechanics formula library generate 2D elastic-plast

5、ic program2.Modify the PRE parameter3.Add Mohr-Coulomb yield function4.Modeling and meshing5.Calculating and analysisSolving stepModify PRE filematedata nespa aeq41.0e10;0.3;0;0;0.2;1.0e3;1.0;3000;0.6; nespa all20;100; nespb beq41.0e10;0.3;0;0;0.2;1.0e3;1.0;3000;0.6;matedata nespa aeq41.0e10;0.3;0;0

6、;0.2;1.0e3;0.0; nespa all20;1600; nespb beq41.0e10;0.3;0;0;0.2;1.0e3;0.0;1.Using solid mechanics formula library generate 2D elastic-plastic program2.Modify the PRE parameter3.Add Mohr-Coulomb yield function4.Modeling and meshing5.Calculating and analysisSolving stepdefidisp u vcoor x ycoef un,vnfun

7、c exx eyy exyshap %1 %2gaus %3$c6 dimension de(3),e(2,2),d(3,3),dv(3),dg(2,2)$c6 dimension dp(3,3),p(4)$c6 external prager Drucker-Prager function$c6 logical filestr$c6 common /gpstr/ gpstr(6,100000),str(6)nespa.pde $cv de(1) = e(1,1)$cv de(2) = e(2,2)$c6 de(3) = e(1,2)+e(2,1)$c6 call getstr(p,str,d

8、e,dv,d,dp,1,prag,prager)$c6 if (prag.gt.-p(2)*0.0001) then$c6 do i=1,3$c6 do j=1,3$c6 d(i,j)=d(i,j)-dp(i,j) nespa.pdeDrucker-Prager functionDrucker-PragerPlastic matrix CalculationStres.for & plas.forCalculate plastic matrix DpCalculate plastic stressCodes of yield function lies in plas.forPlas.for

9、real*8 FUNCTION PRAGER(X,U) implicit real*8 (a-h,o-z) DIMENSION U(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) - D S(2) = U(2) - D S(3) = U(3) Source codes ofDrucker-Prager functionIn plas.forDrucker-Prager functionnespa.pde UK = (u(1)+u(2)*pv-D SJ = (S(

10、1)*2+S(2)*2+UK*2)/2.+S(3)*2 alpha=tgo/dsqrt(9+12*tgo*2) ck=3*c/dsqrt(9+12*tgo*2) PRAGER = DSQRT(SJ) + ALPHA*D*3 - CK RETURN ENDSource codes ofDrucker-Prager functionIn plas.forAdd M-C function in plas.forMohrCoulomb yield function Add M-C function in plas.for real*8 FUNCTION MohrC(X,U) implicit real

11、*8 (a-h,o-z) DIMENSION U(4),X(4),EMSTR(3) emstr(1)=U(1)+U(2) coefb=-U(1)-U(2) coefc=U(1)*U(2)-U(3)*2 emstr(2)=(-coefb)+dsqrt(coefb*2-4*coefc)/2 emstr(3)=(-coefb)-dsqrt(coefb*2-4*coefc)/2 if (emstr(1).lt.emstr(2) then etemp=emstr(2) emstr(2)=emstr(1) emstr(1)=etemp end ifSource codes ofMohr-Coulomb f

12、unctionIn plas.forMohr-Coulomb functionAdd M-C function in plas.for if (emstr(1).lt.emstr(3) then etemp=emstr(3) emstr(3)=emstr(1) emstr(1)=etemp end if if (emstr(2).lt.emstr(3) then etemp=emstr(2) emstr(2)=emstr(3) emstr(3)=etemp end if tgo=x(1) c=x(2) cosf=1./dsqrt(1.+tgo*tgo)Source codes ofMohr-C

13、oulomb functionIn plas.forAdd M-C function in plas.for MohrC=emstr(1)-emstr(3)-2.0*c*cosf MohrC=MohrC+(emstr(1)+emstr(3)*(dsqrt(1-cosf*cosf) RETURN ENDSource codes ofMohr-Coulomb functionIn plas.forCompile plas.forfl32.exe /W0 /nologo /c plas.FORRight click mouse button on plas.for in FEPG interface

14、Then rundefidisp u vcoor x ycoef un,vnfunc exx eyy exyshap q 4gaus q$c6 dimension de(3),e(2,2),d(3,3),dv(3),dg(2,2)$c6 dimension dp(3,3),p(4)$c6 external MohrC$c6 logical filestr$c6 common /gpstr/ gpstr(6,100000),str(6)Modify nespa.pdeModify prager to MohrC $cv de(1) = e(1,1)$cv de(2) = e(2,2)$c6 de

15、(3) = e(1,2)+e(2,1)$c6 call getstr(p,str,de,dv,d,dp,1,prag, MohrC)$c6 if (prag.gt.-p(2)*0.0001) then$c6 do i=1,3$c6 do j=1,3$c6 d(i,j)=d(i,j)-dp(i,j) Modify nespa.pdeModify prager to MohrCCompile nespa.pdeClick run command button in FEPG interfaceCompile nespa.pdeIn command line runpde nespa aeq4Com

16、pile nespa.pdeLet.bat nespaRight click mouse button on nespa in FEPG interfaceThen rundisp sxx,syy,sxy,wp,vep,vyp,didcoor x ycoef un,vnshap q 4gaus qmass q vol$c6 dimension de(3),e(2,2),dv(3),dg(2,2)$c6 dimension dp(3,3),d(3,3),p(4)$c6 external MohrC$c6 logical filestr$c6 common /gpstr/ gpstr(6,1000

17、00),str(6)Modify nespb.pdeModify prager to MohrC $cv de(1) = e(1,1)$cv de(2) = e(2,2)$c6 de(3) = e(1,2)+e(2,1)$c6 call getstr(p,str,de,dv,d,dp,1,prag, MohrC)$c6 if (prag.gt.-p(2)*0.0001) then$c6 eid=0.0d0$c6 else$c6 eid=1.0d0$c6 endif Modify nespb.pdeModify prager to MohrCCompile nespb.pdeClick run

18、command button in FEPG interfaceCompile nespb.pdeIn command line runpde nespb beq4Compile nespb.pdeLet.bat nespbRight click mouse button on nespb in FEPG interfaceThen run1.Modify PRE file modify parameters2.Modify plas.for add Mohr-Coulomb function3.Modify nespa.pde modify yield function to M-C for

19、 displacement calculation4.Modify nespb.pde modify yield function to M-C for stress calculationSummarize5.Compile plas.for fl32.exe /W0 /nologo /c plas.FOR6.Compile nespa.pde pde nespa aeq47.Compile nespb.pde pde nespb beq48.Link to nespa.exe let.bat nespa9.Link to nespb.exe let.bat nespb Summarize1

20、.Using solid mechanics formula library generate 2D elastic-plastic program2.Modify the PRE parameter3.Add Mohr-Coulomb yield function4.Modeling and meshing5.Calculating and analysisSolving stepOne kind of material for the plateAssign materialDisplacement boundary conditionsAssign boundary conditions

21、Quadrilateral Meshing1.Using solid mechanics formula library generate 2D elastic-plastic program2.Modify the PRE parameter3.Add Mohr-Coulomb yield function4.Modeling and meshing5.Calculating and analysisSolving step1st Principal stress 3rd Principal stressResultTheory solution and calculate solution

22、CompareThe calculate Radius of plastic zone D0.130 mThe theory radius of plastic zone is 0.1293 m。 Compareplastic zoneelastic zone2. Simulating method of anchorSimulating anchorglobal coor. oxy & local coor. xEquations of anchorEquilibrium equations of anchor in local coordinate systemRecast to weak

23、 formGES & GLTPrepare 2 filesGES fileGLT fileDescript PDE in local coordinate systemTransform local coordinate system to global coordinate systemGES & GLTGESGLTPDE in local coor.Transform to global Coor.Assembling in global Coor.automaticGES filetrull2defidisp u,vvar u1,v1,u2,v2,refc rx,coor x,func

24、= ex,dord 1,1node 2$c6 pe=prmt(1)$c6 pa=prmt(2)$c6 fu=prmt(3)$c6 time=prmt(4)Trull2.gesGES file$c6 dt=prmt(5)$c6 imate=prmt(6)+0.5$c6 ielem=prmt(7)+0.5$c6 nelem=prmt(8)+0.5shapu=u1=(1.-rx)/2u2=(1.+rx)/2v=v1=(1.-rx)/2v2=(1.+rx)/2Trull2.gesGES filetranx=x(1)=(1.-rx)/2x(2)=(1.+rx)/2gaus 2-1. 1. 1. 1.fu

25、ncex=+u/xTrull2.gesGES filestifdist=+ex;ex*pe*pa+v;v*0.0load=+u*fuendTrull2.gesGLT filedefigsub = trugl2lsub = trull2gdim = 2ldim = 1gvar gu1,gv1,gu2,gv2,lvar lu1,lv1,lu2,lv2,node = 2vartlu1 = T1 gu1 gv1lu2 = T1 gu2 gv2lv1 = T2 gu1 gv1lv2 = T2 gu2 gv2Trugl2.gltExampleAnchor in rockElastic RockConsid

26、er self-weightErock=1.0E8Elastic anchorEanchor=1.0E10Anchor in rockGCN & GIO NFE file and PDE file name of rock , analysis flowPDE PDE description of rockPRE give considered element typeGCN & GIOdefia ell STARTc aSOLVc asdispe #elemtype q42dxyAnchor.gcnAnchor.gioPDE filedefidisp u vcoor x yfunc exx eyy exyshap %1 %2gaus %3mate pe pv fu fv 1.0e10;