用ansys模拟隧洞的分析实例

94117239 · 发表于 2006-5-2 19:43

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

谁有这方面的实例或者是电子书，麻烦上传一下，先谢了

AaronSpark · 发表于 2006-5-3 00:29

FINI /CLE /FILENAME,'TUNNEL' /UNIT,SI /CONFIG,NRES,1E6 !--------------------------------------------------------------------------------------------------------------------! !* TM_YS=16*1E9 MU_YS=0.24 C_YS=1.80*1E6 FI_YS=43 MD_YS=0.027*1E6/10 !* !*C20 TM_GXW=26*1E9 MU_GXW=0.21 C_GXW=2.4*1E6 FI_GXW=50 MD_GXW=0.02235*1E6/10 !* !* C10*! TM_HNT=22*1E9 MU_HNT=0.20 C_HNT=1.73*1E6 FI_HNT=49 MD_HNT=0.024*1E6/10 !* TM_MG=TM_YS MU_MG=MU_YS MD_MG=MD_YS NT_MG=90*1E3 C_MG=C_YS*1.10 FI_MG=FI_YS+5 ZH_MG=NT_MG/(1.2*1*1) !----------------------------------------------------------------------------------------------------------------------! !* R_HJ=6.230 H_ZQ=2.670 L_ZTZ=40 L_YTZ=40 Y_YJMIN=400 Y_YJMAX=415 L_YJ=Y_YJMAX-Y_YJMIN L_MG=3.5 JJ_MG=5.0 PJ_MG=5.0 H_CENDB=1.47 PD=2/1000 W_GXW=0.2 H_JS=831 DMJL_SD=40 H_CTOB=1.47 *DIM,CEN_ZSD,ARRAY,2,1 *DIM,CEN_YSD,ARRAY,2,1 CEN_ZSD(1,1)=-DMJL_SD/2,880 CEN_YSD(1,1)=DMJL_SD/2,880 !+H_CTOB H_ZQ=2.67 FLAG_RD=3 JL_RD=5 *DIM,Z_RDCC,ARRAY,3,2 Z_RDCC(1,1)=3,3,3 Z_RDCC(1,2)=CEN_ZSD(1,1),CEN_ZSD(2,1)+R_HJ+JL_RD+Z_RDCC(2,1)/2,L_ZTZ+7.5 *DIM,Y_RDCC,ARRAY,3,2 Y_RDCC(1,1)=3,3,3 Y_RDCC(1,2)=CEN_YSD(1,1),CEN_YSD(2,1)+R_HJ+JL_RD+Y_RDCC(2,1)/2,L_ZTZ+7.5 !* *DIM,DBKZDX,ARRAY,9,1 *DIM,DBKZDZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+3,1 DBKZDX(1,1)=-150,-64,-38,-20,0,20,42,64,150 DBKZDZ(1,1)=0,L_ZTZ DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)=L_ZTZ+Y_YJMAX-Y_YJMIN+L_YTZ *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG DBKZDZ(I+2,1)=L_ZTZ+PJ_MG*I *ENDDO *DIM,DBKZDGC,ARRAY,9,(Y_YJMAX-Y_YJMIN)/PJ_MG+3 DBKZDGC(1,1)=929,915,915,920,925,930,937,940,940 DBKZDGC(1,2)=929,915,915,920,925,930,937,940,940 DBKZDGC(1,3)=929,915,915,920,925,930,937,940,940 DBKZDGC(1,4)=929,915,915,920,925,930,937,940,940 DBKZDGC(1,5)=929,915,915,920,925,930,937,940,940 DBKZDGC(1,6)=929,915,915,920,925,930,937,940,940 !*************************************************************************************! !* !----------------------------------------------------------------------------------------------------------------------! /PREP7 ET,1,95 ET,2,43 KEYOPT,2,5,2 R,1,W_GXW 
MP,EX,7,TM_YS MP,PRXY,7,MU_YS MP,DENS,7,MD_YS MP,EX,1,TM_YS MP,PRXY,1,MU_YS MP,DENS,1,MD_YS MP,EX,2,TM_GXW MP,PRXY,2,MU_GXW MP,DENS,2,MD_GXW MP,EX,3,TM_MG MP,PRXY,3,MU_MG MP,DENS,3,MD_MG MP,EX,4,TM_HNT MP,PRXY,4,MU_HNT MP,DENS,4,MD_HNT MP,EX,5,TM_YS MP,PRXY,5,MU_YS MP,DENS,5,MD_YS MP,EX,6,TM_YS MP,PRXY,6,MU_YS TB,DP,1,,, TBMODIF,1,1,C_YS TBMODIF,1,2,FI_YS TB,DP,2,,, TBMODIF,1,1,C_GXW TBMODIF,1,2,FI_GXW TB,DP,3,,, TBMODIF,1,1,C_MG TBMODIF,1,2,FI_MG TB,DP,4,,, TBMODIF,1,1,C_HNT TBMODIF,1,2,FI_HNT TB,DP,5,,, TBMODIF,1,1,C_YS TBMODIF,1,2,FI_YS TB,DP,6,,, TBMODIF,1,1,C_YS TBMODIF,1,2,FI_YS !**************************************************! MAT_YS=1 MAT_GXW=2 MAT_MG=3 MAT_DB=4 MAT_SD=5 MAT_MGKW=6 MAT_TXYS=1 KBEGIN=0 *DO,I,1,9-1 *DO,J,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 K,KBEGIN+1,DBKZDX(I,1),H_JS,DBKZDZ(J,1) K,KBEGIN+2,DBKZDX(I,1),DBKZDGC(I,J),DBKZDZ(J,1) K,KBEGIN+3,DBKZDX(I,1),DBKZDGC(I,J+1),DBKZDZ(J+1,1) K,KBEGIN+4,DBKZDX(I,1),H_JS,DBKZDZ(J+1,1) K,KBEGIN+5,DBKZDX(I+1,1),H_JS,DBKZDZ(J,1) K,KBEGIN+6,DBKZDX(I+1,1),DBKZDGC(I+1,J),DBKZDZ(J,1) K,KBEGIN+7,DBKZDX(I+1,1),DBKZDGC(I+1,J+1),DBKZDZ(J+1,1) K,KBEGIN+8,DBKZDX(I+1,1),H_JS,DBKZDZ(J+1,1) V,KBEGIN+1,KBEGIN+2,KBEGIN+3,KBEGIN+4,KBEGIN+5,KBEGIN+6,KBEGIN+7,KBEGIN+8
 KBEGIN=KBEGIN+8 *ENDDO *ENDDO ALLS VGLUE,ALL NUMCMP,ALL VATT,MAT_YS !* CYL4,CEN_ZSD(1,1),CEN_ZSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN CYL4,CEN_YSD(1,1),CEN_YSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN BLOCK,CEN_ZSD(1,1)-R_HJ-L_MG,CEN_ZSD(1,1)+R_HJ+L_MG,CEN_ZSD(2,1)-H_ZQ,CEN_ZSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1) BLOCK,CEN_YSD(1,1)-R_HJ-L_MG,CEN_YSD(1,1)+R_HJ+L_MG,CEN_YSD(2,1)-H_ZQ,CEN_YSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1) !* VSEL,ALL VSEL,S,MAT,,MAT_YS CM,V_BSUB,VOLU VSEL,INVE CM,V_SUB,VOLU ALLS VSBV,V_BSUB,V_SUB CMDEL,V_BSUB CMDEL,V_SUB ALLS NUMCMP,ALL VATT,MAT_YS !* *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 WPAVE,CEN_ZSD(1,1),CEN_ZSD(2,1),DBKZDZ(I,1) CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1) CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) WPAVE,CEN_YSD(1,1),CEN_YSD(2,1),DBKZDZ(I,1) CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1) CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1) *ENDDO
ALLS VGLUE,ALL NUMCMP,ALL CM,ROCK,VOLU
WPAVE,0,0,0 *IF,FLAG_RD,EQ,-1,OR,FLAG_RD,EQ,2,THEN VSTART=1000 NUMSTR,VOLU,VSTART BLOCK,Z_RDCC(1,2)-Z_RDCC(1,1)/2,Z_RDCC(1,2)+Z_RDCC(1,1)/2,Z_RDCC(2,2)-Z_RDCC(2,1)/2,Z_RDCC(2,2)+Z_RDCC(2,1)/2,Z_RDCC(3,2)-Z_RDCC(3,1)/2,Z_RDCC(3,2)+Z_RDCC(3,1)/2 VSBV,ROCK,VSTART NUMCMP,VOLU ALLS CMDEL,ROCK CM,ROCK,VOLU *ENDIF *IF,FLAG_RD,EQ,1,OR,FLAG_RD,EQ,2,THEN VSTART=1000 NUMSTR,VOLU,VSTART BLOCK,Y_RDCC(1,2)-Y_RDCC(1,1)/2,Y_RDCC(1,2)+Y_RDCC(1,1)/2,Y_RDCC(2,2)-Y_RDCC(2,1)/2,Y_RDCC(2,2)+Y_RDCC(2,1)/2,Y_RDCC(3,2)-Y_RDCC(3,1)/2,Y_RDCC(3,2)+Y_RDCC(3,1)/2 VSBV,ROCK,VSTART NUMCMP,VOLU ALLS CMDEL,ROCK CM,ROCK,VOLU *ENDIF ALLS NUMCMP,ALL !* !----------------------------------------------------------------------------------------------------------------------! !* ALLS VATT,MAT_YS ALLS VSEL,S,LOC,X,CEN_ZSD(1,1)-1E-7,CEN_ZSD(1,1)+1E-7 CM,CM_TEMPV,VOLU ASLV,R ASEL,R,LOC,Y,CEN_ZSD(2,1)-H_ZQ-1E-7,CEN_ZSD(2,1)-H_ZQ+1E-7 VSLA,R CM,CM_ZDB,VOLU CMSEL,S,CM_TEMPV ASLV,S ASEL,R,LOC,Y,CEN_ZSD(2,1)-1E-7,CEN_ZSD(2,1)+1E-7 ASEL,U,LOC,X,CEN_ZSD(1,1)-R_HJ,CEN_ZSD(1,1)-R_HJ-L_MG ASEL,U,LOC,X,CEN_ZSD(1,1)+R_HJ,CEN_ZSD(1,1)+R_HJ+L_MG VSLA,R CM,CM_ZSD,VOLU CMSEL,S,CM_ZSD CMSEL,A,CM_ZDB VATT,MAT_SD CM,CM_ZSDKW,VOLU ALLS CMSEL,S,CM_TEMPV VSEL,U,MAT,,MAT_SD ASLV,R ASEL,R,LOC,Y,CEN_ZSD(2,1) VSEL,ALL VSLA,R VATT,MAT_MGKW CM,CM_ZMG,VOLU !右隧道 ALLS VSEL,S,LOC,X,CEN_YSD(1,1)-1E-7,CEN_YSD(1,1)+1E-7 CM,CM_TEMPV,VOLU ASLV,R ASEL,R,LOC,Y,CEN_YSD(2,1)-H_ZQ-1E-7,CEN_YSD(2,1)-H_ZQ+1E-7 VSLA,R CM,CM_YDB,VOLU CMSEL,S,CM_TEMPV ASLV,S ASEL,R,LOC,Y,CEN_YSD(2,1)-1E-7,CEN_YSD(2,1)+1E-7 ASEL,U,LOC,X,CEN_YSD(1,1)-R_HJ,CEN_YSD(1,1)-R_HJ-L_MG ASEL,U,LOC,X,CEN_YSD(1,1)+R_HJ,CEN_YSD(1,1)+R_HJ+L_MG VSLA,R CM,CM_YSD,VOLU CMSEL,S,CM_YSD CMSEL,A,CM_YDB VATT,MAT_SD CM,CM_YSDKW,VOLU ALLS CMSEL,S,CM_TEMPV VSEL,U,MAT,,MAT_SD ASLV,R ASEL,R,LOC,Y,CEN_YSD(2,1) VSEL,ALL VSLA,R VATT,MAT_MGKW CM,CM_YMG,VOLU 
ALLS
CMDEL,CM_TEMPV ALLS ASEL,S,LOC,X,DBKZDX(1,1) ASEL,A,LOC,X,DBKZDX(9,1) DA,ALL,UX,0 ALLS ASEL,S,LOC,Z,DBKZDZ(1,1) ASEL,A,LOC,Z,DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1) DA,ALL,UZ,0 ALLS ASEL,S,LOC,Y,H_JS DA,ALL,ALL ALLS
LTS_ZD=2 LTS_YD=2 *DIM,ZLT_HFX,ARRAY,LTS_ZD,1 ZLT_HFX(1,1)=5,3 !4,2 *DIM,YLT_HFX,ARRAY,LTS_YD,1 YLT_HFX(1,1)=5,3 !4,2 *DIM,LT_HFZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+2 LT_HFZ(1,1)=4 *DO,I,2,(Y_YJMAX-Y_YJMIN)/PJ_MG+1 LT_HFZ(I,1)=2 !PJ_MG *ENDDO LT_HFZ((Y_YJMAX-Y_YJMIN)/PJ_MG+2,1)=4 !10 
*DO,J,1,LTS_ZD *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 X_VCEN=(DBKZDX(J,1)+DBKZDX(J+1,1))/2 Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2 VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7 VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7 ASLV,R LSLA,R CM,CM_HFX,LINE LSEL,R,LOC,X,X_VCEN LESIZE,ALL,,,ZLT_HFX(J,1),,1 CMSEL,S,CM_HFX LSEL,R,LOC,Z,Z_VCEN LESIZE,ALL,,,LT_HFZ(I,1),,1 CMSEL,S,CM_HFX LSEL,U,LOC,X,X_VCEN LSEL,U,LOC,Z,Z_VCEN LESIZE,ALL,,,9,,1 CMSEL,S,CM_HFX MSHAPE,0,3D MSHKEY,1 VMESH,ALL ALLS *ENDDO !右端 *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 X_VCEN=(DBKZDX(9-J+1,1)+DBKZDX(9-J,1))/2 Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2 VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7 VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7 ASLV,R LSLA,R CM,CM_HFX,LINE LSEL,R,LOC,X,X_VCEN LESIZE,ALL,,,YLT_HFX(J,1),,1 CMSEL,S,CM_HFX LSEL,R,LOC,Z,Z_VCEN LESIZE,ALL,,,LT_HFZ(I,1),,1 CMSEL,S,CM_HFX LSEL,U,LOC,X,X_VCEN LSEL,U,LOC,Z,Z_VCEN LESIZE,ALL,,,10,,1 CMSEL,S,CM_HFX MSHAPE,0,3D MSHKEY,1 VMESH,ALL ALLS *ENDDO
*ENDDO *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 CMSEL,S,CM_ZSDKW ASLV,R LSLA,R
 X_VCEN=CEN_ZSD(1,1) Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2 VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7 CM,CM_VMESH,VOLU ASLV,R LSLA,R CM,CM_HFX,LINE ALLS CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7 LESIZE,ALL,,,2,,1 CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7 LESIZE,ALL,,,2,,1 ALLS
 CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7 LSEL,U,LOC,Z,Z_VCEN LESIZE,ALL,2 ALLS CMSEL,S,CM_HFX LSEL,R,LOC,Z,Z_VCEN LESIZE,ALL,,,LT_HFZ(I,1),,1 ALLS
 CMSEL,S,CM_HFX LSEL,U,LOC,X,X_VCEN LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7 LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7 LSEL,U,LOC,Z,Z_VCEN LESIZE,ALL,,,1,,1 ALLS
 ALLS CMSEL,S,CM_VMESH MSHAPE,0,3D VSWEEP,ALL ALLS *ENDDO
 CMSEL,S,CM_ZMG ESIZE,2 VSWEEP,ALL ALLS *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 CMSEL,S,CM_YSDKW ASLV,R LSLA,R
 X_VCEN=CEN_YSD(1,1) Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2 VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7 CM,CM_VMESH,VOLU ASLV,R LSLA,R CM,CM_HFX,LINE CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7 LESIZE,ALL,,,2,,1 CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7 LESIZE,ALL,,,2,,1 
 CMSEL,S,CM_HFX LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7 LSEL,U,LOC,Z,Z_VCEN LSEL,U,LOC,Y,CEN_YSD(2,1)-H_CTOB LESIZE,ALL,2 CMSEL,S,CM_HFX LSEL,R,LOC,Z,Z_VCEN LESIZE,ALL,,,LT_HFZ(I,1),,1
 CMSEL,S,CM_HFX LSEL,U,LOC,X,X_VCEN LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7 LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7 LSEL,U,LOC,Z,Z_VCEN LESIZE,ALL,,,1,,1 ALLS CMSEL,S,CM_VMESH MSHAPE,0,3D VSWEEP,ALL ALLS 
 CMSEL,S,CM_YMG ESIZE,2 VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7 VSWEEP,ALL ALLS *ENDDO

ALLS VSEL,U,MAT,,MAT_SD VSEL,U,MAT,,MAT_MGKW VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(9-LTS_YD,1) CM,CM_VMESH,VOLU
!* VSEL,S,,,CM_VMESH VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+2,1) ASLV,S LSLA,S LSEL,R,LOC,X,(DBKZDX(LTS_ZD+1,1)+DBKZDX(LTS_ZD+2,1))/2 LESIZE,ALL,,,5 !* VSEL,S,,,CM_VMESH VSEL,R,LOC,X,DBKZDX(LTS_ZD+2,1),DBKZDX(LTS_ZD+3,1) ASLV,S LSLA,S LSEL,R,LOC,X,(DBKZDX(LTS_ZD+2,1)+DBKZDX(LTS_ZD+3,1))/2 LESIZE,ALL,,,3 ALLS !* VSEL,S,,,CM_VMESH VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+4,1) ASLV,S LSLA,S LSEL,R,LOC,X,(DBKZDX(LTS_ZD+3,1)+DBKZDX(LTS_ZD+4,1))/2 LESIZE,ALL,,,3 !* VSEL,S,,,CM_VMESH VSEL,R,LOC,X,DBKZDX(LTS_ZD+4,1),DBKZDX(LTS_ZD+5,1) ASLV,S LSLA,S LSEL,R,LOC,X,(DBKZDX(LTS_ZD+4,1)+DBKZDX(LTS_ZD+5,1))/2 LESIZE,ALL,,,5 ALLS !* ALLS *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 !CMSEL,S,CM_VMESH VSEL,S,,,CM_VMESH VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1) VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+3,1) CM,CM_NEXTV,VOLU
 CMSEL,S,CM_NEXTV ESIZE,4 VSWEEP,ALL ALLS *ENDDO
ALLS LSEL,S,LOC,X,Y_RDCC(1,2)-Y_RDCC(1,1)/2-1E-7,Y_RDCC(1,2)+Y_RDCC(1,1)/2+1E-7 LSEL,R,LOC,Y,Y_RDCC(2,2)-Y_RDCC(2,1)/2-1E-7,Y_RDCC(2,2)+Y_RDCC(2,1)/2+1E-7 LSEL,R,LOC,Z,Y_RDCC(3,2)-Y_RDCC(3,1)/2-1E-7,Y_RDCC(3,2)+Y_RDCC(3,1)/2+1E-7 LESIZE,ALL,1 ALLS
 *DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1 VSEL,S,,,CM_VMESH VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1) CM,CM_NEXTV,VOLU
 CMSEL,S,CM_NEXTV VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+3+2,1) ESIZE,4 VSWEEP,ALL ALLS *ENDDO ALLS
VSEL,S,,,CM_ZSDKW ASLV,R VSEL,S,,,CM_ZMG ASLV,R CM,CM_ZGXW,AREA TYPE,2 MAT,MAT_GXW AMESH,ALL ALLS
VSEL,S,,,CM_YSDKW ASLV,R VSEL,S,,,CM_YMG ASLV,R CM,CM_YGXW,AREA TYPE,2 MAT,MAT_GXW AMESH,ALL ALLS
NUMCMP,ALL ALLS MPCHG,MAT_TXYS,ALL ALLS
SAVE,TNLMESH,DB

AaronSpark · 发表于 2006-5-3 00:31

!------------------------------------------------------- 开始求解 ---------------------------------------! /SOLU !* ACEL,,9.8 ANTYPE,STATIC,NEW TIME,1 NROPT,FULL AUTOTS,ON DELTIM,0.5,0.01,0.5 LNSRCH,ON PRED,ON ISWRITE,1 ALLS !* ALLS ESEL,S,TYPE,,2 EKILL,ALL !* ALLS SOLV FINI !*
/CLE /FILNAME,CALTUN /PREP7 RESUME,TNLMESH,DB FINI !* /SOLU !* ACEL,,9.8 !* ANTYPE,STATIC,NEW TIME,1 NROPT,FULL AUTOTS,ON DELTIM,0.5,0.01,0.5 LNSRCH,ON PRED,ON ESEL,U,TYPE,,2 ISFILE,READ,TUNNEL,IST,,0 !* ALLS ESEL,S,TYPE,,2 EKILL,ALL!* CNVTOL,F,0.05 ALLS !* SOLV
!---------------------------------------------------------------------- 开挖模拟 ---------------------------------------------------! !----------------------------------- 第一步开挖，这一步没有支护 ----------------------------------! ALLS VSEL,S,,,CM_ZSDKW VSEL,A,,,CM_YSDKW VSEL,R,LOC,Z,DBKZDZ(1,1),DBKZDZ(2,1) ESLV,R !ESEL,A,TYPE,,2 EKILL,ALL TIME,2 NROPT,FULL AUTOTS,ON DELTIM,0.3,0.001,0.5 LNSRCH,ON PRED,ON NCNV,2 ALLS SOLV SAVE,KW-1,DB !* !---------------------------- 循环做第二步到第五步的开挖，第一步到第四步的支护 ----------------------------! !* ALLS ESEL,S,TYPE,,2 EKILL,ALL !---------------------------- 第I步开挖,第I-1步支护 ---------------------! *DIM,KWD_TOU,ARRAY,4 *DIM,KWD_WEI,ARRAY,4 *DIM,ZHD_TOU,ARRAY,4 *DIM,ZHD_WEI,ARRAY,4 KWD_TOU(1,1)=3,4,5,6 KWD_WEI(1,1)=2,3,4,5 ZHD_TOU(1,1)=2,3,4,5 ZHD_WEI(1,1)=1,2,3,4 *DO,I,1,4 !* TOU=KWD_TOU(I,1) WEI=KWD_WEI(I,1) !* ALLS VSEL,S,,,CM_ZSDKW VSEL,A,,,CM_YSDKW VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLV,R EKILL,ALL !* !* TOU=ZHD_TOU(I,1) WEI=ZHD_WEI(I,1) !* ALLS VSEL,S,,,CM_ZDB VSEL,A,,,CM_YDB VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLV,R EALIVE,ALL MPCHG,MAT_DB,ALL !* ALLS VSEL,S,,,CM_ZMG VSEL,A,,,CM_YMG VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLV,R MPCHG,MAT_MG,ALL !* ALLS ASEL,S,,,CM_ZGXW ASEL,A,,,CM_YGXW ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLA,R EALIVE,ALL MPCHG,MAT_GXW,ALL !* ALLS ASEL,S,,,CM_ZGXW ASEL,A,,,CM_YGXW ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLA,S SFE,ALL,,PRES,,ZH_MG ALLS !* TIME,I+2 NROPT,FULL AUTOTS,ON DELTIM,0.3,0.001,0.5 LNSRCH,ON PRED,ON NCNV,2 !OUTRES,ALL,ALL ALLS SOLV SAVE,KW-%I+1%,DB !* *ENDDO !----------------------------------- 第5步支护，这一步没有开挖 ----------------------------------! TOU=6 WEI=5 !* ALLS VSEL,S,,,CM_ZDB VSEL,A,,,CM_YDB VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLV,R EALIVE,ALL MPCHG,MAT_DB,ALL !* ALLS VSEL,S,,,CM_ZMG VSEL,A,,,CM_YMG VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLV,R MPCHG,MAT_MG,ALL !* ALLS ASEL,S,,,CM_ZGXW ASEL,A,,,CM_YGXW ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) ESLA,R EALIVE,ALL MPCHG,MAT_GXW,ALL ALLS ASEL,S,,,CM_ZGXW ASEL,A,,,CM_YGXW ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1) SFA,ALL,,PRES,ZH_MG ALLS !* ALLS TIME,7 NROPT,FULL AUTOTS,ON DELTIM,0.3,0.01,0.5 LNSRCH,ON PRED,ON !OUTRES,ALL,ALL ALLS SOLV SAVE,KW-6,DB /EXIT,NOSAVE

AaronSpark · 发表于 2006-5-3 00:34

隧道上部挖方或填方的宏 NBpre.mac !!=================================== !!=================================== !!=================================== !!=================================== /com,The analysis of nearby construction fini /clear,nostart /filname,NYcon !!=================================== !!=================================== !!=================================== !!=================================== !!=================================== !填写衬砌的参数 multipro,'start',7 *cset,1,3,E_num1,'Lining element number',3 *cset,4,6,Ex_mat1,'Youngs modulus(Pa)',1e10 *cset,7,9,Nxy_mat1,'Poissons ratio',0.3 *cset,10,12,Den_mat1,'Density(kg/m^3)',3000 *cset,13,15,Dep_mat1,'Height of Lining(m)',0.5 *cset,16,18,e_sizex,'Number size of Lining(x)',15 *cset,19,21,e_sizey,'Number size of Lining(y)',10 *cset,61,62,'Enter the attibutes of','Lining' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写混凝土参数，包括单轴抗拉极限强度和单轴抗压极限强度 multipro,'start',2 *cset,1,3,ConcrRl,'Tensile cracking stress(MPa)',1.3 *cset,4,6,ConcrRa,'crushing stress(MPa)',10.5 *cset,61,62,'Enter the attibutes of','Concrete' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写围岩的参数 multipro,'start',7 *cset,1,3,E_num2,'Lining element number',42 *cset,4,6,Ex_mat2,'Youngs modulus(Pa)',1e9 *cset,7,9,Nxy_mat2,'Poissons ratio',0.3 *cset,10,12,Co_mat2,'Cohesion(Pa)',1e4 *cset,13,15,FrA_mat2,'Friction Angle(0)',20 *cset,16,18,FlA_mat2,'Flow Angle(0)',0 *cset,19,21,Den_mat2,'Density(kg/m^3)',3000 *cset,61,62,'Enter the attibutes of',' Block' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写注浆参数 multipro,'start',9 *cset,1,3,E_num3,'Zhujiang element number',E_num2 *cset,4,6,Ex_mat3,'Youngs modulus(Pa)',Ex_mat2 *cset,7,9,Nxy_mat3,'Poissons ratio',Nxy_mat2 *cset,10,12,Co_mat3,'Cohesion(Pa)',Co_mat2 *cset,13,15,FrA_mat3,'Friction Angle(0)',FrA_mat2 *cset,16,18,FlA_mat3,'Flow Angle(0)',FlA_mat2 *cset,19,21,Den_mat3,'Density(kg/m^3)',Den_mat2 *cset,22,24,DH,'Zhujiang houdu(m)',3 *cset,25,27,Zj_step,'Zhujiang step',0 *cset,61,62,'Enter the attibutes of','Zhujiang' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写衬砌的半径以及围岩的取值范围 !以及地表到隧道中心的距离 multipro,'start',4 *cset,1,3,R_Lining,'Radius of the Lining',5 *cset,4,6,Width,'Width of the Rock range',100 *cset,7,9,Height,'Height of the Rock range',60 *cset,10,12,H_R,'Distance of surface and Lining',20 *cset,61,62,'Enter the attibutes of',' Location' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写开挖的高度及分几层采挖 multipro,'start',2 *cset,1,3,Cave,'Height of the cavation',20 *cset,4,6,Step,'How many cavationg step',4 *cset,61,62,'Enter the attibutes of',' cavation' *cset,63,64,'Note:Unit of Each Parameter!' multipro,'end' !!=================================== !!=================================== !填写每层加载或卸载高度 *do,i,1,Step *ask,H_step%i%,Input %i% step height(m) *enddo !!=================================== !!=================================== !指定是开挖还是卸载 *ask,Code,Cavation or Load(Y or N) !!=================================== !!=================================== !确定土层释放系数数组 *do,i,1,Step CaveCoff%i%=H_step%i%/Cave *enddo !!=================================== !!=================================== !!=================================== !!=================================== !!=================================== !进入前处理器 /prep7 shpp,off !!=================================== !定义衬砌支护参数 et,1,E_num1 mp,ex,1,Ex_mat1 mp,nuxy,1,Nxy_mat1 mp,dens,1,Den_mat1 sectype,1,beam,rect,chenqi,0 secoffset,cent secdata,1,Dep_mat1 !!=================================== !!=================================== !定义实常数 r,1,Dep_mat1,Dep_mat1/12,Dep_mat1 !!=================================== !!=================================== !定义围岩参数 et,2,E_num2 keyopt,2,3,2 mp,ex,2,Ex_mat2 mp,nuxy,2,Nxy_mat2 mp,dens,2,Den_mat2 tb,dp,2 tbdata,1,Co_mat2,FrA_mat2,FlA_mat2 !!=================================== !!=================================== !定义注浆参数 et,3,E_num3 keyopt,3,3,2 mp,ex,3,Ex_mat3 mp,nuxy,3,Nxy_mat3 mp,dens,3,Den_mat3 tb,dp,3 tbdata,1,Co_mat3,FrA_mat3,FlA_mat3

AaronSpark · 发表于 2006-5-3 00:35

!!=================================== !建模 !!=================================== !!=================================== !!给衬砌建模 !将衬砌分成相等的四个部分 !建立衬砌的四个关键点 *afun,deg !改变系统角度的单位 csys,1 k,1,R_Lining,40 k,2,R_Lining,140 k,3,R_Lining,220 k,4,R_Lining,320 csys,0 larc,1,2,3,R_Lining !线1为顶拱 larc,2,3,1,R_Lining !线2为左边墙 larc,3,4,1,R_Lining !线3为仰拱 larc,4,1,2,R_Lining !线4为右边墙 type,1 mat,1 real,1 lsel,s,,,1,3,2 lesize,all,,,e_sizex allsel,all lsel,s,,,2,4,2 lesize,all,,,e_sizey lsel,all lmesh,all cm,Chenqi,elem cm,Chenqind,node lsel,none esel,none nsel,none !!=================================== !!=================================== !!=================================== !!=================================== !!给注浆层建模 csys,1 k,5,R_Lining+DH,40 k,6,R_Lining+DH,140 k,7,R_Lining+DH,220 k,8,R_Lining+DH,320 csys,0 larc,5,6,7,R_Lining+DH !线5为注浆顶 larc,6,7,5,R_Lining+DH !线6为注浆左侧 larc,7,8,5,R_Lining+DH !线7为注浆底 larc,8,5,6,R_Lining+DH !线8为注浆右侧 lsel,s,,,5,7,2 lesize,all,,,e_sizex lsel,all lsel,s,,,6,8,2 lesize,all,,,e_sizey lsel,none lstr,1,5 lstr,2,6 lstr,3,7 lstr,4,8 lesize,all,,,3,1.5 lsel,all al,1,9,5,10 al,2,10,6,11 al,3,11,7,12 al,4,12,8,9 type,2 mat,2 mshape,0,2d mshkey,1 allsel,all esel,none amesh,all cm,Zhujiang,elem lsel,none asel,none esel,none nsel,none !!=================================== !!=================================== !!=================================== !!=================================== !给围岩建模 k,9,-Width/2,H_R k,10,Width/2,H_R lsel,none lstr,9,10 !l=13 lesize,all,,,e_sizex,-0.6 lsel,all !!卸载 *if,lwcase(Code),eq,'y',then !此处判断是否为卸载 *do,i,1,Step,1 *if,i,eq,1,then H=CaveCoff1*Cave *else H=H+CaveCoff%i%*Cave *endif k,10+(i-1)*2+1,-Width/2,H_R-H k,10+(i-1)*2+2,Width/2,H_R-H lsel,none *if,i,eq,1,then lstr,9,10+(i-1)*2+1 lstr,10,10+(i-1)*2+2 *else lstr,10+(i-1)*2+1,10+(i-2)*2+1 lstr,10+(i-1)*2+2,10+(i-2)*2+2 *endif *if,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,2,then lesize,all,,,1 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,2,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,4 lesize,all,,,2 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,4,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,8 lesize,all,,,3 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,8,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,15 lesize,all,,,4 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,15,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,30 lesize,all,,,6 *else lesize,all,,,8 *endif lsel,none lstr,10+(i-1)*2+1,10+(i-1)*2+2 lesize,all,,,e_sizex,-0.6 lsel,all *enddo lsel,all *do,i,1,Step,1 *if,i,eq,1,then H=CaveCoff1*Cave lsel,s,loc,x,-Width/2,Width/2 lsel,r,loc,y,H_R-H,H_R al,all lsel,all *else H=H+CaveCoff%i%*Cave lsel,s,loc,x,-Width/2,Width/2 lsel,r,loc,y,H_R-H,H_R-H+CaveCoff%i%*Cave al,all lsel,all *endif *enddo *get,kmax,kp,,num,max k,kmax+1,-Width/2,-(Height-H_R) k,kmax+2,Width/2,-(Height-H_R) lsel,none lstr,kmax+1,kmax+2 lesize,all,,,e_sizex,-0.6 lsel,none lstr,kmax+1,kp(-Width/2,H_R-Cave,0) lesize,all,,,e_sizey lsel,none lstr,kmax+2,kp(Width/2,H_R-Cave,0) lesize,all,,,e_sizey lsel,all *get,lmax,line,,num,max lsel,none lstr,5,kp(Width/2,H_R-Cave,0) lstr,6,kp(-Width/2,H_R-Cave,0) lstr,7,kmax+1 lstr,8,kmax+2 lesize,all,,,e_sizey,3.5 lsel,all al,5,lmax+1,13+3*Step,lmax+2 al,6,lmax+2,15+3*Step,lmax+3 al,7,lmax+3,14+3*Step,lmax+4 al,8,lmax+4,16+3*Step,lmax+1 !!以下将划分围岩网格 type,2 mat,2 mshape,0,2d mshkey,1 allsel,all cmsel,s,Zhujiang cmsel,u,Chenqind asel,s,,,4+Step+1,4+Step+4,1 amesh,all *do,i,1,Step *if,i,eq,1,then H=CaveCoff1*Cave !用于储存开挖深度 *else H=H+CaveCoff%i%*Cave *endif allsel,all esel,none nsel,none asel,s,,,4+i,4+i amesh,all cm,Sub%i%,elem nsel,s,loc,y,R_H-H,R_H-H cm,Subnd%i%,node asel,all *enddo allsel,all cm,Rock0,elem *do,i,1,Step,1 cmsel,u,Sub%i% cm,Rock%i%,elem *enddo eplot !!加载 *elseif,lwcase(Code),eq,'n' !此处判断是否为加载 *do,i,1,Step,1 *if,i,eq,1,then H=CaveCoff1*Cave *else H=H+CaveCoff%i%*Cave *endif k,10+(i-1)*2+1,-Width/2,H_R+H k,10+(i-1)*2+2,Width/2,H_R+H lsel,none *if,i,eq,1,then lstr,9,10+(i-1)*2+1 lstr,10,10+(i-1)*2+2 *else lstr,10+(i-1)*2+1,10+(i-2)*2+1 lstr,10+(i-1)*2+2,10+(i-2)*2+2 *endif *if,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,1,then lesize,all,,,1 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,1,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,3 lesize,all,,,1 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,3,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,10 lesize,all,,,2 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,10,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,18 lesize,all,,,3 *elseif,distkp(10+(i-1)*2+1,10+(i-2)*2+1),gt,18,and,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,30 lesize,all,,,6 *else lesize,all,,,8 *endif lsel,none lstr,10+(i-1)*2+1,10+(i-1)*2+2 lesize,all,,,e_sizex,-0.6 lsel,all *enddo allsel,all *do,i,1,Step,1 *if,i,eq,1,then H=CaveCoff1*Cave lsel,s,loc,x,-Width/2,Width/2 lsel,r,loc,y,H_R,H_R+H al,all lsel,all *else H=H+CaveCoff%i%*Cave lsel,s,loc,x,-Width/2,Width/2 lsel,r,loc,y,H_R+H-CaveCoff%i%*Cave,H_R+H al,all lsel,all *endif *enddo *get,kmax,kp,,num,max k,kmax+1,-Width/2,-(Height-H_R) k,kmax+2,Width/2,-(Height-H_R) lsel,none lstr,kmax+1,kmax+2 lesize,all,,,e_sizex,-0.6 lsel,none lstr,kmax+1,9 lesize,all,,,e_sizey lsel,none lstr,kmax+2,10 lesize,all,,,e_sizey lsel,all *get,lmax,line,,num,max lsel,none lstr,5,10 lstr,6,9 lstr,7,kmax+1 lstr,8,kmax+2 lesize,all,,,e_sizey,3.5 lsel,all al,5,lmax+1,13,lmax+2 al,6,lmax+2,15+3*Step,lmax+3 al,7,lmax+3,14+3*Step,lmax+4 al,8,lmax+4,16+3*Step,lmax+1 lsel,all !!以下将划分围岩网格 type,2 mat,2 mshape,0,2d mshkey,1 *do,i,1,Step,1 *if,i,eq,1,then H=0 !用于储存开挖深度 *else H=H+CaveCoff%i-1%*Cave *endif allsel,all esel,none nsel,none asel,s,,,4+i,4+i amesh,all cm,Sub%i%,elem nsel,s,loc,y,H_R+H cm,Subnd%i%,node asel,all *enddo allsel,all *do,i,1,Step,1 cmsel,u,Sub%i% *enddo cm,Rock0,elem cmsel,s,Rock0 *do,i,1,Step,1 cmsel,a,Sub%i% cm,Rock%i%,elem *enddo allsel,all asel,s,,,4+Step+1,4+Step+4,1 amesh,all *else *endif allsel,all !check,esel,err eplot fini

AaronSpark · 发表于 2006-5-3 00:35

!!=================================== !!=================================== !!=================================== !!=================================== /eof NBsol.mac !!=================================== !!=================================== !!=================================== !!=================================== !!以下将进入求解器计算 /solu allsel,all !!=================================== !!=================================== !以下步骤将提取开挖或加载交界面处的节点号 !!提取节点号 !!!!!!!!!!!!每条线上的节点数应等于e_sizex+1!!!!!!!!!!!!!!!!!!!!! !!!!!!!交界线的总数应与加载和卸载的步骤数Step相同!!!!!!!!!!!!!!!!!! !!!!!!!定义一数组Node_num(e_sizex+1,Step)来存储所有的交界线节点!!!!!!!!! allsel,all *dim,Node_num,array,e_sizex+1,Step *if,lwcase(Code),eq,'y',then *do,i,1,Step,1 *if,i,eq,1,then H=CaveCoff1*Cave *else H=H+CaveCoff%i%*Cave *endif nsel,s,loc,y,H_R-H,H_R-H *do,j,1,e_sizex+1,1 *get,nmin,node,,num,min Node_num(j,i)=nmin nsel,u,,,Node_num(j,i) *enddo allsel,all *enddo *elseif,lwcase(Code),eq,'n' *do,i,1,Step,1 *if,i,eq,1,then H=0 *else H=H+CaveCoff%i-1%*Cave *endif nsel,s,loc,y,H_R+H,H_R+H *do,j,1,e_sizex+1,1 *get,nmin,node,,num,min Node_num(j,i)=nmin nsel,u,,,Node_num(j,i) *enddo allsel,all *enddo *else *endif allsel,all eplot !!=================================== !!=================================== !全局求解控制 !!求解器选择 antype,static,new eqslv,,, nlgeom,on nropt,full,,off !激活自适应下降 cnvtol,f cnvtol,m kbc,0 ncnv,2 !当计算不收敛时，中断计算，但不结束程序 autots,on rescontrol,define,all,-1,2 outres,all,1 outpr,all,1 allsel,all ddele,all,all ealive,all !!用一个三维数组NRf来存储等效节点力（释放荷载），维数 !!为(e_sizex+1)×4×Step,第一列是节点号，第二至第四列分别 !!存储Fx,Fy,Mz;面号对应施工步。 !!****************************************************** !!用一个三维数组NRf来存储等效节点力（释放荷载），维数 !!为(e_sizex+1)×4×Step,第一列是节点号，第二至第四列分别 !!存储Fx,Fy,Mz;面号对应施工步。 *if,lwcase(Code),eq,'y',then *dim,NRf,array,e_sizex+1,4,Step !该数组仅在开挖卸载有用 *endif !!****************************************************** !正式求解 *do,i,1,Step+1,1 /solu allsel,all *if,i,eq,1,then time,1 *else antype,static,rest parres,new,parafile,dat time,i deltim,0.1,0.01,0.2,1 *endif !!****改变注浆层材料***** *if,Zj_step,eq,i-1,then esel,s,,,Zhujiang mpchg,3,all allsel,all *endif !!*********************** *if,lwcase(Code),eq,'n',then !!加载计算 *if,i,eq,1,then esel,none *do,j,1,Step,1 esel,a,,,Sub%j% *enddo ekill,all esel,s,live nsle,s nsel,inve d,all,ux,0 d,all,uy,0 nsel,s,loc,y,-(Height-H_R),-(Height-H_R) d,all,uy,0 nsel,s,loc,x,-Width/2,-Width/2 nsel,r,loc,y,-(Height-H_R),H_R d,all,ux,0 nsel,s,loc,x,Width/2,Width/2 nsel,r,loc,y,-(Height-H_R),H_R d,all,ux,0 allsel,all acel,,9.8 *else allsel,all esel,none *do,j,1,i-1,1 esel,a,,,Sub%j% *enddo ealive,all nsle,s ddele,all,all nsel,s,loc,x,-Width/2,-Width/2 nsel,a,loc,x,Width/2,Width/2 d,all,ux,0 *endif allsel,all *endif *if,lwcase(Code),eq,'y',then !!卸载计算 !!****第一步求解原始地貌***** *if,i,eq,1,then allsel,all nsel,s,loc,x,-Width/2,-Width/2 nsel,a,loc,x,Width/2,Width/2 d,all,ux,0 nsel,all nsel,s,loc,y,-(Height-H_R),-(Height-H_R) d,all,uy,0 allsel,all acel,,9.8 !!*************************** !!****求解开挖步************* *elseif,i,gt,1 esel,none *do,j,1,i-1,1 esel,a,,,Sub%j% *enddo ekill,all esel,s,live nsle,s nsel,inve d,all,ux,0 d,all,uy,0 allsel,all ! *do,j,1,e_sizex+1,1 !在以下可以考虑分级释放，但 ! *if,i,ge,3,then !考虑到挖方本身的支护效应没有 ! fdele,NRf(j,1,i-2),fx !故省略 ! fdele,NRf(j,1,i-2),fy ! fdele,NRf(j,1,i-2),mz ! *else ! *endif ! f,NRf(j,1,i-1),fx,-NRf(j,2,i-1) ! f,NRf(j,1,i-1),fy,-NRf(j,3,i-1) ! f,NRf(j,1,i-1),mz,-NRf(j,4,i-1) ! *enddo *else *endif allsel,all *endif allsel,all /pbc,all,,1 solve

AaronSpark · 发表于 2006-5-3 00:37

!!************************************************* !!进入后处理器得到该步计算得到的开挖释放荷载,命令仅对开挖卸载有效 *if,lwcase(Code),eq,'y',and,i,le,Step,then fini /post1 allsel,all set,i esel,s,,,Sub%i% *do,j,1,e_sizex+1,1 NRf(j,1,i)=Node_num(j,i) nsel,s,,,NRf(j,1,i) spoint,NRf(j,1,i) fsum,,, *get,NRf(j,2,i),fsum,0,item,fx *get,NRf(j,3,i),fsum,0,item,fy *get,NRf(j,4,i),fsum,0,item,mz *enddo fini *else *endif !!************************************************* !!************************************************* !以下两步极其重要 i=i+1 parsav,all,parafile,dat !!************************************************* *enddo fini /eof NBpost.mac !!=================================== !!=================================== !!=================================== !!=================================== !在处理中将得到衬砌单元的节点的每一步的力，施工位移，以及单元的安全系数 !以下进入后处理模块 fini /post1 !首先定义各数组 *dim,ElemNum,array,2*(e_sizex+e_sizey),3,1 !数组ElemNum存储衬砌单元和每个单元的I,J节点号 !维数i=2×(e_sizex+e_sizey),j=3.第一列是节点号，第二列是总位移，第三列是位移的方向 *dim,NdData,array,2*(e_sizex+e_sizey),5,Step !数组NdData存储衬砌所有节点号，坐标以及每个施工步位移和方向 !维数是i=2×(e_sizex+e_sizey),j=6,k=Step.第1列是节点号，2，3列是坐标，第四列是位移，第五列是位移的方向，面号对 !应施工步 *dim,ElemDatI,array,2*(e_sizex+e_sizey),7,Step+1 !数组ElemDatI存储单元I节点的数据 !维数为i=2×(e_sizex+e_sizey),j=7,k=Step+1,第一列是单元号，第二列是轴力，第三列是弯矩，第四列是偏心距，第五列是 !安全系数，第六列是安全与否，第七列是受拉还是受压 *dim,ElemDatJ,array,2*(e_sizex+e_sizey),7,Step+1 !数组ElemDatJ存储单元I节点的数据 !维数为i=2×(e_sizex+e_sizey),j=7,k=Step+1,第一列是单元号，第二列是轴力，第三列是弯矩，第四列是偏心距，第五列是 !安全系数，第六列是安全与否，第七列是受拉还是受压 !*************************************** !初始化数组ElemNum esel,s,,,Chenqi nsle,s *do,i,1,2*(e_sizex+e_sizey),1 *get,EMin,elem,,num,min ElemNum(i,1)=EMin ElemNum(i,2)=nelem(ElemNum(i,1),1) ElemNum(i,3)=nelem(ElemNum(i,1),2) esel,u,,,EMin *enddo allsel,all !数组ElemNum的数据已经全部填充完毕 !*************************************** !*************************************** !初始化数组NdData esel,s,,,Chenqi nsle,s csys,1 *do,i,1,2*(e_sizex+e_sizey),1 *get,NMin,node,,num,min *do,j,1,Step,1 NdData(i,1,j)=NMin NdData(i,2,j)=nx(NdData(i,1,j)) NdData(i,3,j)=ny(NdData(i,1,j)) *enddo nsel,u,,,NMin *enddo csys,0 allsel,all !*************************************** !*************************************** !初始化数组ElemDatI和ElemDatJ *do,i,1,2*(e_sizex+e_sizey),1 *do,j,1,Step+1,1 ElemDatI(i,1,j)=ElemNum(i,1) ElemDatJ(i,1,j)=ElemNum(i,1) *enddo *enddo !*************************************** !*************************************** !填充NdData(节点号与坐标已经初始化），得到衬砌所有节点在每个施工步以及累加之前施工步的位移位移 lcdef,1,1 *do,i,2,Step+1,1 lcdef,i,i lcase,i lcoper,sub,1 *do,j,1,2*(e_sizex+e_sizey),1 dispx=ux(NdData(j,1)) dispy=uy(NdData(j,1)) Alldisp=sqrt(dispx**2+dispy**2) DTheta=atan(dispy/dispx) NdData(j,4,i-1)=Alldisp NdData(j,5,i-1)=DTheta *enddo *enddo !数组NdData的数据已经全部填充完毕 !*************************************** !*************************************** !填充ElemDatI和ElemDatJ(单元号已经初始化），得到衬砌单元的所有参数 allsel,all *do,i,1,Step+1,1 set,i *do,j,1,2*(e_sizex+e_sizey),1 *get,forxI,elem,ElemDatI(j,1,i),smisc,1 *get,formI,elem,ElemDatI(j,1,i),smisc,6 *get,forxJ,elem,ElemDatJ(j,1,i),smisc,7 *get,formJ,elem,ElemDatJ(j,1,i),smisc,12 ElemDatI(j,2,i)=forxI !I轴力 ElemDatI(j,3,i)=formI !I弯矩 ElemDatI(j,4,i)=abs(formI/forxI) !I偏心距 ElemDatJ(j,2,i)=forxJ !J轴力 ElemDatJ(j,3,i)=formJ !J弯矩 ElemDatJ(j,4,i)=abs(formJ/forxJ) !J偏心距 !****以下得到安全系数****************** !ElemDatI a=1-1.5*abs(formI/forxI)/Dep_mat1 !轴力偏心影响系数 *if,abs(formI/forxI),le,0.2,then safecoff=a*ConcrRa*1e6*1*Dep_mat1/abs(forxI) !安全系数 ElemDatI(j,5,i)=safecoff *if,safecoff,ge,2.4,then ElemDatI(j,6,i)=1 !1表示安全 *elseif,safecoff,lt,2.4 ElemDatI(j,6,i)=0 !0表示不安全 *else *endif ElemDatI(j,7,i)=1 !1表示衬砌受压控制 *elseif,abs(formI/forxI),gt,0.2 safecoff=1.75*ConcrRl*1e6*1*Dep_mat1**2/abs(forxI*(6*abs(formI/forxI)-Dep_mat1)) ElemDatI(j,5,i)=safecoff *if,safecoff,ge,3.6,then ElemDatI(j,6,i)=1 !1表示安全 *elseif,safecoff,lt,3.6 ElemDatI(j,6,i)=0 !0表示不安全 *else *endif ElemDatI(j,7,i)=0 !0表示衬砌受拉控制 *endif !ElemDatJ a=1-1.5*abs(formJ/forxJ)/Dep_mat1 !轴力偏心影响系数 *if,abs(formJ/forxJ),le,0.2,then safecoff=a*ConcrRa*1e6*1*Dep_mat1/abs(forxJ) !安全系数 ElemDatJ(j,5,i)=safecoff *if,safecoff,ge,2.4,then ElemDatJ(j,6,i)=1 !1表示安全 *elseif,safecoff,lt,2.4 ElemDatJ(j,6,i)=0 !0表示不安全 *else *endif ElemDatJ(j,7,i)=1 !1表示衬砌受压控制 *elseif,abs(formJ/forxJ),gt,0.2 safecoff=1.75*ConcrRl*1e6*1*Dep_mat1**2/abs(forxI*(6*abs(formI/forxI)-Dep_mat1)) ElemDatJ(j,5,i)=safecoff *if,safecoff,ge,3.6,then ElemDatJ(j,6,i)=1 !1表示安全 *elseif,safecoff,lt,3.6 ElemDatJ(j,6,i)=0 !0表示不安全 *else *endif ElemDatJ(j,7,i)=0 !0表示衬砌受拉控制 *endif *enddo *enddo !!*****************以下得到优化状态变量***********************!! !总的施工位移状态变量SigDisp,是开挖或卸载导致的总的变形量,由数组NdData得到 SigDisp=NdData(1,4,Step) *do,i,1,2*(e_sizex+e_sizey),1 *if,NdData(i,4,Step),gt,DispTota,then SigDisp=NdData(i,4,Step) *else *endif *enddo !单元I端节点受拉或受压时的安全系数状态变量，分别共有的数量为Step+1个，用SigYaI%i% !表示受压，用SigLaI%i%表示受拉，i代表施工步。由数组ElemDatI得到 *do,i,1,Step+1,1 !首先在第i步的数组ElemDatI中得到一个受压，一个受拉单元数据,该步骤非常重要，它能够 !搜索该施工步中是否有受压或受拉单元，或者两者都有 *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatI(j,7,i),eq,1,then SigYaI%i%=2.4 *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatI(j,7,i),eq,0,then SigLaI%i%=3.6 *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatI(j,7,i),eq,1,and,ElemDatI(j,5,i),Lt,SigYaI%i%,then SigYaI%i%=ElemDatI(j,5,i) *else *endif *if,ElemDatI(j,7,i),eq,0,and,ElemDatI(j,5,i),Lt,SigLaI%i%,then SigLaI%i%=ElemDatI(j,5,i) *else *endif *enddo *enddo !单元J端节点受拉或受压时的安全系数状态变量，分别共有的数量为Step+1个，用SigYaJ%i% !表示受压，用SigLaJ%i%表示受拉，i代表施工步。由数组ElemDatJ得到 *do,i,1,Step+1,1 !首先在第i步的数组ElemDatI中得到一个受压，一个受拉单元数据,该步骤非常重要，它能够 !搜索该施工步中是否有受压或受拉单元，或者两者都有 *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatJ(j,7,i),eq,1,then SigYaJ%i%=2.4 *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatJ(j,7,i),eq,0,then SigLaJ%i%=3.6 *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatJ(j,7,i),eq,1,and,ElemDatJ(j,5,i),Lt,SigYaJ%i%,then SigYaJ%i%=ElemDatJ(j,5,i) *else *endif *if,ElemDatJ(j,7,i),eq,0,and,ElemDatJ(j,5,i),Lt,SigLaJ%i%,then SigLaJ%i%=ElemDatJ(j,5,i) *else *endif *enddo *enddo fini /eof

AaronSpark · 发表于 2006-5-3 00:38

NBopt.mac !!=================================== !!=================================== !!=================================== !!=================================== !进入优化处理器 /opt opclr opanl,'NBOptDesign','mac' opvar,Dep_mat1,dv,0.01,1 opvar,Cave,dv,0.1,H_R-2*R_Lining *do,i,1,Step+1,1 *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatI(j,7,i),eq,1,then opvar,SigYaI%i%,sv,2.4,40 !opvar,SigYaI%i%,obj *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatI(j,7,i),eq,0,then opvar,SigLaI%i%,sv,3.6,40 !opvar,SigLaI%i%,obj *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatJ(j,7,i),eq,1,then opvar,SigYaJ%i%,sv,2.4,40 !opvar,SigYaJ%i%,obj *exit *else *endif *enddo *do,j,1,2*(e_sizex+e_sizey),1 *if,ElemDatJ(j,7,i),eq,0,then opvar,SigLaJ%i%,sv,3.6,40 !opvar,SigLaJ%i%,obj *exit *else *endif *enddo *enddo opdata,,, oploop,prep,proc,all opprnt,on opkeep,on !optype,run !opfrst,45 !opexe /eof

AaronSpark · 发表于 2006-5-3 00:40

!又一个隧道模拟开挖命令流 /com fini /cle *set,x1,-30 *set,y1,-40 *set,w,64.9 *set,h,58.3 *set,th,0.5 /prep7 k,,0,0 $k,,0,3.85 k,,0.88,5.5 $k,,2.45,6.15 k,,4.02,5.5 $k,,4.9,3.85 $k,,4.9,0 larc,1,2,6,8.13 $larc,2,3,6,3.21 larc,3,4,6,2.22 $larc,4,5,2,2.22 larc,5,6,2,3.21 $larc,6,7,2,8.13 $larc,7,1,4,6 a,1,2,3,4,5,6,7 blc4,x1,y1,w,h aovl,2,1 lsel,s,line,,1,6, lcomb,all,, lsel,all, et,1,beam3 et,2,plane42 keyopt,2,3,2 r,1,th,th*th*th/12,th, , mp,ex,1,2.6e7 mp,prxy,1,0.2 mp,dens,1,25 mp,ex,2,1.3e6 mp,prxy,2,0.32 tb,dp,2 tbdata,1,200,30, mp,dens,2,22 mp,ex,3,1.3e6 mp,prxy,3,0.32 tb,dp,3 tbdata,1,200,30, mp,dens,3,22 lsel,s,line,,1 lesize,all,,,50, latt,1,1,1 lmesh,all alls lsel,s,line,,7 lesize,all,,,17, latt,1,1,1 lmesh,all alls lsel,s,line,,8,11 lesize,all,,,40 alls mopt,split,on mopt,split,err mopt,split,warn mopt,qmesh,main asel,s,area,,3 aatt,2,,2 amesh,all alls asel,s,area,,1 aatt,3,,2 amesh,all alls lsel,s,loc,x,x1 lsel,a,loc,x,x1+w dl,all,,ux,0 alls lsel,s,loc,y,y1 dl,all,,uy,0 alls acel,,1 fini /solu antype,static deltim,0.1,0.05,0.2 autots,on !使用自动时间步 pred,on !打开时间步长预测器 lnsrch,on !打开线性搜索 nlgeom,on !打开大位移效果 nropt,full !设定牛顿-拉普森选项 cnvtol,f,,0.02,0 esel,s,type,,1 ekill,all esel,all esel,s,live nsle,s nsel,inve d,all,all,0 nsel,all esel,all solve esel,s,mat,,3 ekill,all esel,s,type,,1 ealive,all nsle,s ddele,all,all esel,all esel,s,live nsle,s nsel,inve d,all,all,0 nsel,all esel,all solve fini /post1 etable,if,smisc,1 etable,jf,smisc,7 etable,im,smisc,6 etable,jm,smisc,12

94117239 · 发表于 2006-5-3 17:54

真是谢谢你了，我现在遇到的情况这样的，我们先对四类围岩进行固结灌浆，然后达到三类的标准，在进行喷锚支护，这中间钢筋，喷射混凝土，锚杆，改造后的围岩应该采用那种单元来模拟呢，谢谢指导

94117239 · 发表于 2006-5-3 17:56

还有接触面的约束＆围岩地影响范围应该怎么取

AaronSpark · 发表于 2006-5-4 07:44

围岩可以solid单元 锚杆用杆单元 拱架可用梁单元

94117239 · 发表于 2006-5-4 09:49

谢谢指导，我在仔细琢磨一下

94117239 · 发表于 2006-5-9 10:49

锚杆怎么加荷载啊，面荷载怎么加，锚杆的参数怎么取

追风 · 发表于 2010-7-28 22:24

非常感谢啊，救了我的命啊！

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 用ansys模拟隧洞的分析实例

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例

回复：（94117239）用ansys模拟隧洞的分析实例