声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5038|回复: 15

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

[复制链接]
发表于 2006-5-2 19:43 | 显示全部楼层 |阅读模式

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

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

x
谁有这方面的实例或者是电子书,麻烦上传一下,先谢了
回复
分享到:

使用道具 举报

发表于 2006-5-3 00:29 | 显示全部楼层

回复:(94117239)用ansys模拟隧洞的分析实例

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

<P>ALLS<BR>VSEL,U,MAT,,MAT_SD<BR>VSEL,U,MAT,,MAT_MGKW<BR>VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(9-LTS_YD,1)<BR>CM,CM_VMESH,VOLU</P>
<P>!*<BR>VSEL,S,,,CM_VMESH<BR>VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+2,1)   <BR>ASLV,S<BR>LSLA,S<BR>LSEL,R,LOC,X,(DBKZDX(LTS_ZD+1,1)+DBKZDX(LTS_ZD+2,1))/2<BR>LESIZE,ALL,,,5   <BR>!*<BR>VSEL,S,,,CM_VMESH<BR>VSEL,R,LOC,X,DBKZDX(LTS_ZD+2,1),DBKZDX(LTS_ZD+3,1) <BR>ASLV,S<BR>LSLA,S<BR>LSEL,R,LOC,X,(DBKZDX(LTS_ZD+2,1)+DBKZDX(LTS_ZD+3,1))/2<BR>LESIZE,ALL,,,3   <BR>ALLS<BR>!*<BR>VSEL,S,,,CM_VMESH<BR>VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+4,1)  <BR>ASLV,S<BR>LSLA,S<BR>LSEL,R,LOC,X,(DBKZDX(LTS_ZD+3,1)+DBKZDX(LTS_ZD+4,1))/2<BR>LESIZE,ALL,,,3   <BR>!*<BR>VSEL,S,,,CM_VMESH<BR>VSEL,R,LOC,X,DBKZDX(LTS_ZD+4,1),DBKZDX(LTS_ZD+5,1)  <BR>ASLV,S<BR>LSLA,S<BR>LSEL,R,LOC,X,(DBKZDX(LTS_ZD+4,1)+DBKZDX(LTS_ZD+5,1))/2<BR>LESIZE,ALL,,,5   <BR>ALLS<BR>!*<BR>ALLS<BR>*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1<BR> !CMSEL,S,CM_VMESH<BR> VSEL,S,,,CM_VMESH<BR> VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)<BR> VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+3,1)<BR> CM,CM_NEXTV,VOLU</P>
<P> CMSEL,S,CM_NEXTV<BR> ESIZE,4<BR> VSWEEP,ALL<BR> ALLS<BR>*ENDDO</P>
<P>ALLS<BR>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<BR>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<BR>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<BR>LESIZE,ALL,1   <BR>ALLS</P>
<P><BR>*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1<BR> VSEL,S,,,CM_VMESH<BR> VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)<BR> CM,CM_NEXTV,VOLU</P>
<P> CMSEL,S,CM_NEXTV<BR> VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+3+2,1)<BR> ESIZE,4<BR> VSWEEP,ALL<BR> ALLS<BR>*ENDDO<BR>ALLS</P>
<P>VSEL,S,,,CM_ZSDKW   <BR>ASLV,R<BR>VSEL,S,,,CM_ZMG<BR>ASLV,R<BR>CM,CM_ZGXW,AREA<BR>TYPE,2<BR>MAT,MAT_GXW<BR>AMESH,ALL<BR>ALLS</P>
<P>VSEL,S,,,CM_YSDKW   <BR>ASLV,R<BR>VSEL,S,,,CM_YMG<BR>ASLV,R<BR>CM,CM_YGXW,AREA<BR>TYPE,2<BR>MAT,MAT_GXW<BR>AMESH,ALL<BR>ALLS</P>
<P>NUMCMP,ALL<BR>ALLS<BR>MPCHG,MAT_TXYS,ALL <BR>ALLS</P>
<P>SAVE,TNLMESH,DB  <BR></P>
发表于 2006-5-3 00:31 | 显示全部楼层

回复:(94117239)用ansys模拟隧洞的分析实例

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

回复:(94117239)用ansys模拟隧洞的分析实例

<STRONG>隧道上部挖方或填方的宏<BR><BR></STRONG>NBpre.mac<BR>!!===================================<BR>!!===================================<BR>!!===================================<BR>!!===================================<BR>/com,The analysis of nearby construction<BR>fini<BR>/clear,nostart<BR>/filname,NYcon<BR>!!===================================<BR>!!===================================<BR>!!===================================<BR>!!===================================<BR><BR>!!===================================<BR>!填写衬砌的参数<BR>multipro,'start',7<BR>  *cset,1,3,E_num1,'Lining element number',3<BR>  *cset,4,6,Ex_mat1,'Youngs modulus(Pa)',1e10<BR>  *cset,7,9,Nxy_mat1,'Poissons ratio',0.3<BR>  *cset,10,12,Den_mat1,'Density(kg/m^3)',3000<BR>  *cset,13,15,Dep_mat1,'Height of Lining(m)',0.5<BR>  *cset,16,18,e_sizex,'Number size of Lining(x)',15<BR>  *cset,19,21,e_sizey,'Number size of Lining(y)',10<BR>  *cset,61,62,'Enter the attibutes of','Lining'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写混凝土参数,包括单轴抗拉极限强度和单轴抗压极限强度<BR>multipro,'start',2<BR>  *cset,1,3,ConcrRl,'Tensile cracking stress(MPa)',1.3<BR>  *cset,4,6,ConcrRa,'crushing stress(MPa)',10.5<BR>  *cset,61,62,'Enter the attibutes of','Concrete'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写围岩的参数<BR>multipro,'start',7<BR>  *cset,1,3,E_num2,'Lining element number',42<BR>  *cset,4,6,Ex_mat2,'Youngs modulus(Pa)',1e9<BR>  *cset,7,9,Nxy_mat2,'Poissons ratio',0.3<BR>  *cset,10,12,Co_mat2,'Cohesion(Pa)',1e4<BR>  *cset,13,15,FrA_mat2,'Friction Angle(0)',20<BR>  *cset,16,18,FlA_mat2,'Flow Angle(0)',0<BR>  *cset,19,21,Den_mat2,'Density(kg/m^3)',3000<BR>  *cset,61,62,'Enter the attibutes of',' Block'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写注浆参数<BR>multipro,'start',9<BR>  *cset,1,3,E_num3,'Zhujiang element number',E_num2<BR>  *cset,4,6,Ex_mat3,'Youngs modulus(Pa)',Ex_mat2<BR>  *cset,7,9,Nxy_mat3,'Poissons ratio',Nxy_mat2<BR>  *cset,10,12,Co_mat3,'Cohesion(Pa)',Co_mat2<BR>  *cset,13,15,FrA_mat3,'Friction Angle(0)',FrA_mat2<BR>  *cset,16,18,FlA_mat3,'Flow Angle(0)',FlA_mat2<BR>  *cset,19,21,Den_mat3,'Density(kg/m^3)',Den_mat2<BR>  *cset,22,24,DH,'Zhujiang houdu(m)',3<BR>  *cset,25,27,Zj_step,'Zhujiang step',0<BR>  *cset,61,62,'Enter the attibutes of','Zhujiang'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写衬砌的半径以及围岩的取值范围<BR>!以及地表到隧道中心的距离<BR>multipro,'start',4<BR>  *cset,1,3,R_Lining,'Radius of the Lining',5<BR>  *cset,4,6,Width,'Width of the Rock range',100<BR>  *cset,7,9,Height,'Height of the Rock range',60<BR>  *cset,10,12,H_R,'Distance of surface and Lining',20<BR>  *cset,61,62,'Enter the attibutes of',' Location'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写开挖的高度及分几层采挖<BR>multipro,'start',2<BR>  *cset,1,3,Cave,'Height of the cavation',20<BR>  *cset,4,6,Step,'How many cavationg step',4<BR>  *cset,61,62,'Enter the attibutes of',' cavation'<BR>  *cset,63,64,'Note:Unit of Each Parameter!'<BR>multipro,'end'<BR>!!===================================<BR><BR>!!===================================<BR>!填写每层加载或卸载高度<BR><BR>*do,i,1,Step<BR>   *ask,H_step%i%,Input %i% step height(m) <BR>*enddo<BR>!!=================================== <BR><BR>!!===================================<BR>!指定是开挖还是卸载<BR>*ask,Code,Cavation or Load(Y or N)<BR>!!===================================<BR><BR>!!===================================<BR>!确定土层释放系数数组<BR>*do,i,1,Step <BR>   CaveCoff%i%=H_step%i%/Cave<BR>*enddo<BR>!!===================================<BR><BR>!!===================================<BR>!!===================================<BR>!!===================================<BR>!!===================================<BR>!进入前处理器<BR>/prep7<BR><BR>shpp,off<BR>!!===================================<BR>!定义衬砌支护参数<BR>et,1,E_num1<BR>mp,ex,1,Ex_mat1<BR>mp,nuxy,1,Nxy_mat1<BR>mp,dens,1,Den_mat1<BR>sectype,1,beam,rect,chenqi,0<BR>secoffset,cent<BR>secdata,1,Dep_mat1<BR>!!===================================<BR><BR>!!===================================<BR>!定义实常数<BR>r,1,Dep_mat1,Dep_mat1/12,Dep_mat1<BR>!!===================================<BR><BR>!!===================================<BR>!定义围岩参数<BR>et,2,E_num2<BR>keyopt,2,3,2<BR>mp,ex,2,Ex_mat2<BR>mp,nuxy,2,Nxy_mat2<BR>mp,dens,2,Den_mat2<BR>tb,dp,2<BR>tbdata,1,Co_mat2,FrA_mat2,FlA_mat2<BR>!!===================================<BR><BR>!!===================================<BR>!定义注浆参数<BR>et,3,E_num3<BR>keyopt,3,3,2<BR>mp,ex,3,Ex_mat3<BR>mp,nuxy,3,Nxy_mat3<BR>mp,dens,3,Den_mat3<BR>tb,dp,3<BR>tbdata,1,Co_mat3,FrA_mat3,FlA_mat3<BR>
发表于 2006-5-3 00:35 | 显示全部楼层
!!===================================<BR><BR>!建模<BR>!!===================================<BR>!!===================================<BR>!!给衬砌建模<BR>!将衬砌分成相等的四个部分<BR>!建立衬砌的四个关键点<BR><BR>*afun,deg  !改变系统角度的单位<BR>csys,1<BR>k,1,R_Lining,40<BR>k,2,R_Lining,140<BR>k,3,R_Lining,220<BR>k,4,R_Lining,320<BR>csys,0<BR>larc,1,2,3,R_Lining  !线1为顶拱<BR>larc,2,3,1,R_Lining  !线2为左边墙<BR>larc,3,4,1,R_Lining  !线3为仰拱<BR>larc,4,1,2,R_Lining  !线4为右边墙<BR>type,1<BR>mat,1<BR>real,1<BR>lsel,s,,,1,3,2<BR>lesize,all,,,e_sizex<BR>allsel,all<BR>lsel,s,,,2,4,2<BR>lesize,all,,,e_sizey<BR>lsel,all<BR>lmesh,all<BR>cm,Chenqi,elem<BR>cm,Chenqind,node<BR>lsel,none<BR>esel,none<BR>nsel,none<BR>!!===================================<BR>!!===================================<BR><BR>!!===================================<BR>!!===================================<BR>!!给注浆层建模<BR>csys,1<BR>k,5,R_Lining+DH,40<BR>k,6,R_Lining+DH,140<BR>k,7,R_Lining+DH,220<BR>k,8,R_Lining+DH,320<BR>csys,0<BR>larc,5,6,7,R_Lining+DH  !线5为注浆顶<BR>larc,6,7,5,R_Lining+DH  !线6为注浆左侧<BR>larc,7,8,5,R_Lining+DH  !线7为注浆底<BR>larc,8,5,6,R_Lining+DH  !线8为注浆右侧<BR>lsel,s,,,5,7,2<BR>lesize,all,,,e_sizex<BR>lsel,all<BR>lsel,s,,,6,8,2<BR>lesize,all,,,e_sizey<BR>lsel,none<BR>lstr,1,5<BR>lstr,2,6<BR>lstr,3,7<BR>lstr,4,8<BR>lesize,all,,,3,1.5<BR>lsel,all<BR>al,1,9,5,10<BR>al,2,10,6,11<BR>al,3,11,7,12<BR>al,4,12,8,9<BR>type,2<BR>mat,2<BR>mshape,0,2d<BR>mshkey,1<BR>allsel,all<BR>esel,none<BR>amesh,all<BR>cm,Zhujiang,elem<BR>lsel,none<BR>asel,none<BR>esel,none<BR>nsel,none<BR>!!===================================<BR>!!===================================<BR><BR>!!===================================<BR>!!===================================<BR>!给围岩建模<BR>k,9,-Width/2,H_R<BR>k,10,Width/2,H_R<BR>lsel,none<BR>lstr,9,10    !l=13<BR>lesize,all,,,e_sizex,-0.6<BR>lsel,all     <BR><BR>!!卸载<BR>*if,lwcase(Code),eq,'y',then  !此处判断是否为卸载<BR>   *do,i,1,Step,1<BR>     *if,i,eq,1,then<BR>        H=CaveCoff1*Cave<BR>  *else         <BR>    H=H+CaveCoff%i%*Cave<BR>     *endif<BR>     k,10+(i-1)*2+1,-Width/2,H_R-H <BR>     k,10+(i-1)*2+2,Width/2,H_R-H <BR>     lsel,none<BR>     *if,i,eq,1,then<BR>        lstr,9,10+(i-1)*2+1<BR>  lstr,10,10+(i-1)*2+2<BR>  *else<BR>    lstr,10+(i-1)*2+1,10+(i-2)*2+1<BR>          lstr,10+(i-1)*2+2,10+(i-2)*2+2<BR>     *endif<BR>     *if,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,2,then<BR>        lesize,all,,,1<BR>  *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<BR>          lesize,all,,,2<BR>  *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<BR>          lesize,all,,,3<BR>  *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<BR>          lesize,all,,,4<BR>  *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<BR>          lesize,all,,,6<BR>  *else<BR>    lesize,all,,,8<BR>     *endif<BR>     lsel,none<BR>     lstr,10+(i-1)*2+1,10+(i-1)*2+2<BR>     lesize,all,,,e_sizex,-0.6<BR>     lsel,all<BR>   *enddo<BR>   lsel,all<BR>   *do,i,1,Step,1<BR>     *if,i,eq,1,then<BR>       H=CaveCoff1*Cave   <BR>       lsel,s,loc,x,-Width/2,Width/2<BR>       lsel,r,loc,y,H_R-H,H_R<BR>       al,all<BR>       lsel,all<BR>     *else<BR>       H=H+CaveCoff%i%*Cave<BR>       lsel,s,loc,x,-Width/2,Width/2<BR>       lsel,r,loc,y,H_R-H,H_R-H+CaveCoff%i%*Cave<BR>       al,all<BR>       lsel,all<BR>     *endif<BR>   *enddo<BR>   *get,kmax,kp,,num,max<BR>   k,kmax+1,-Width/2,-(Height-H_R)<BR>   k,kmax+2,Width/2,-(Height-H_R)<BR>   lsel,none<BR>   lstr,kmax+1,kmax+2<BR>   lesize,all,,,e_sizex,-0.6<BR>   lsel,none<BR>   lstr,kmax+1,kp(-Width/2,H_R-Cave,0)<BR>   lesize,all,,,e_sizey<BR>   lsel,none<BR>   lstr,kmax+2,kp(Width/2,H_R-Cave,0)<BR>   lesize,all,,,e_sizey<BR>   lsel,all<BR>   *get,lmax,line,,num,max<BR>   lsel,none<BR>   lstr,5,kp(Width/2,H_R-Cave,0)<BR>   lstr,6,kp(-Width/2,H_R-Cave,0)<BR>   lstr,7,kmax+1<BR>   lstr,8,kmax+2<BR>   lesize,all,,,e_sizey,3.5<BR>   lsel,all<BR>   al,5,lmax+1,13+3*Step,lmax+2<BR>   al,6,lmax+2,15+3*Step,lmax+3<BR>   al,7,lmax+3,14+3*Step,lmax+4<BR>   al,8,lmax+4,16+3*Step,lmax+1<BR>   <BR>   !!以下将划分围岩网格<BR>   type,2<BR>   mat,2<BR>   mshape,0,2d<BR>   mshkey,1<BR>   allsel,all<BR>   cmsel,s,Zhujiang<BR>   cmsel,u,Chenqind<BR>   asel,s,,,4+Step+1,4+Step+4,1<BR>   amesh,all<BR>   *do,i,1,Step   <BR>      *if,i,eq,1,then<BR>           H=CaveCoff1*Cave        !用于储存开挖深度<BR>      *else<BR>           H=H+CaveCoff%i%*Cave<BR>      *endif<BR>      allsel,all    <BR>      esel,none<BR>      nsel,none<BR>      asel,s,,,4+i,4+i<BR>      amesh,all<BR>      cm,Sub%i%,elem <BR>      nsel,s,loc,y,R_H-H,R_H-H<BR>      cm,Subnd%i%,node<BR>      asel,all<BR>   *enddo<BR>   allsel,all<BR>   cm,Rock0,elem<BR>   *do,i,1,Step,1<BR>      cmsel,u,Sub%i%<BR>      cm,Rock%i%,elem<BR>   *enddo<BR>   eplot<BR><BR>!!加载<BR>*elseif,lwcase(Code),eq,'n'  !此处判断是否为加载<BR>   *do,i,1,Step,1<BR>     *if,i,eq,1,then<BR>        H=CaveCoff1*Cave<BR>  *else<BR>    H=H+CaveCoff%i%*Cave<BR>     *endif<BR>     k,10+(i-1)*2+1,-Width/2,H_R+H <BR>     k,10+(i-1)*2+2,Width/2,H_R+H <BR>     lsel,none<BR>     *if,i,eq,1,then<BR>        lstr,9,10+(i-1)*2+1<BR>  lstr,10,10+(i-1)*2+2<BR>  *else<BR>    lstr,10+(i-1)*2+1,10+(i-2)*2+1<BR>          lstr,10+(i-1)*2+2,10+(i-2)*2+2<BR>     *endif<BR>     *if,distkp(10+(i-1)*2+1,10+(i-2)*2+1),le,1,then<BR>        lesize,all,,,1<BR>  *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<BR>          lesize,all,,,1<BR>  *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<BR>          lesize,all,,,2<BR>  *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<BR>          lesize,all,,,3<BR>  *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<BR>          lesize,all,,,6<BR>  *else<BR>    lesize,all,,,8<BR>     *endif<BR>     lsel,none<BR>     lstr,10+(i-1)*2+1,10+(i-1)*2+2<BR>     lesize,all,,,e_sizex,-0.6<BR>     lsel,all<BR>   *enddo<BR>   allsel,all<BR>   *do,i,1,Step,1<BR>     *if,i,eq,1,then<BR>       H=CaveCoff1*Cave<BR>       lsel,s,loc,x,-Width/2,Width/2<BR>       lsel,r,loc,y,H_R,H_R+H <BR>       al,all<BR>       lsel,all<BR>     *else<BR>       H=H+CaveCoff%i%*Cave<BR>       lsel,s,loc,x,-Width/2,Width/2<BR>       lsel,r,loc,y,H_R+H-CaveCoff%i%*Cave,H_R+H <BR>       al,all<BR>       lsel,all<BR>     *endif<BR>   *enddo<BR>   *get,kmax,kp,,num,max<BR>   k,kmax+1,-Width/2,-(Height-H_R)<BR>   k,kmax+2,Width/2,-(Height-H_R)<BR>   lsel,none<BR>   lstr,kmax+1,kmax+2<BR>   lesize,all,,,e_sizex,-0.6<BR>   lsel,none<BR>   lstr,kmax+1,9<BR>   lesize,all,,,e_sizey<BR>   lsel,none<BR>   lstr,kmax+2,10<BR>   lesize,all,,,e_sizey<BR>   lsel,all<BR>   *get,lmax,line,,num,max<BR>   lsel,none<BR>   lstr,5,10<BR>   lstr,6,9<BR>   lstr,7,kmax+1<BR>   lstr,8,kmax+2<BR>   lesize,all,,,e_sizey,3.5<BR>   lsel,all<BR>   al,5,lmax+1,13,lmax+2<BR>   al,6,lmax+2,15+3*Step,lmax+3<BR>   al,7,lmax+3,14+3*Step,lmax+4<BR>   al,8,lmax+4,16+3*Step,lmax+1<BR>   lsel,all<BR>  !!以下将划分围岩网格<BR>   type,2<BR>   mat,2<BR>   mshape,0,2d<BR>   mshkey,1<BR>   *do,i,1,Step,1   <BR>      *if,i,eq,1,then<BR>           H=0 !用于储存开挖深度<BR>      *else<BR>           H=H+CaveCoff%i-1%*Cave<BR>      *endif<BR>      allsel,all<BR>      esel,none<BR>      nsel,none<BR>      asel,s,,,4+i,4+i<BR>      amesh,all<BR>      cm,Sub%i%,elem  <BR>      nsel,s,loc,y,H_R+H<BR>      cm,Subnd%i%,node  <BR>      asel,all<BR>   *enddo<BR>   allsel,all<BR>   *do,i,1,Step,1<BR>      cmsel,u,Sub%i%<BR>   *enddo<BR>   cm,Rock0,elem<BR>   cmsel,s,Rock0<BR>   *do,i,1,Step,1<BR>      cmsel,a,Sub%i%<BR>      cm,Rock%i%,elem<BR>   *enddo<BR>   allsel,all<BR>   asel,s,,,4+Step+1,4+Step+4,1<BR>   amesh,all  <BR>*else<BR>*endif<BR>allsel,all<BR>!check,esel,err<BR>eplot<BR>fini<BR>
发表于 2006-5-3 00:35 | 显示全部楼层

回复:(94117239)用ansys模拟隧洞的分析实例

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

回复:(94117239)用ansys模拟隧洞的分析实例

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

回复:(94117239)用ansys模拟隧洞的分析实例

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

回复:(94117239)用ansys模拟隧洞的分析实例

!又一个隧道模拟开挖命令流<BR><BR>/com  <BR>fini  <BR>/cle  <BR>*set,x1,-30  <BR>*set,y1,-40  <BR>*set,w,64.9  <BR>*set,h,58.3  <BR>*set,th,0.5  <BR>  <BR>/prep7  <BR>  <BR>k,,0,0 $k,,0,3.85  <BR>k,,0.88,5.5 $k,,2.45,6.15  <BR>k,,4.02,5.5 $k,,4.9,3.85 $k,,4.9,0  <BR>  <BR>larc,1,2,6,8.13 $larc,2,3,6,3.21  <BR>larc,3,4,6,2.22 $larc,4,5,2,2.22  <BR>larc,5,6,2,3.21 $larc,6,7,2,8.13 $larc,7,1,4,6  <BR>  <BR>a,1,2,3,4,5,6,7  <BR>blc4,x1,y1,w,h  <BR>aovl,2,1  <BR>  <BR>lsel,s,line,,1,6,  <BR>lcomb,all,,  <BR>lsel,all,  <BR>  <BR>et,1,beam3  <BR>et,2,plane42  <BR>keyopt,2,3,2  <BR>r,1,th,th*th*th/12,th, ,  <BR>  <BR>mp,ex,1,2.6e7  <BR>mp,prxy,1,0.2  <BR>mp,dens,1,25  <BR>  <BR>mp,ex,2,1.3e6  <BR>mp,prxy,2,0.32  <BR>tb,dp,2  <BR>tbdata,1,200,30,  <BR>mp,dens,2,22  <BR>  <BR>mp,ex,3,1.3e6  <BR>mp,prxy,3,0.32  <BR>tb,dp,3  <BR>tbdata,1,200,30,  <BR>mp,dens,3,22  <BR>  <BR>lsel,s,line,,1  <BR>lesize,all,,,50,  <BR>latt,1,1,1  <BR>lmesh,all  <BR>alls  <BR>  <BR>lsel,s,line,,7  <BR>lesize,all,,,17,  <BR>latt,1,1,1  <BR>lmesh,all  <BR>alls  <BR>  <BR>lsel,s,line,,8,11  <BR>lesize,all,,,40  <BR>alls  <BR>  <BR>mopt,split,on  <BR>mopt,split,err  <BR>mopt,split,warn  <BR>mopt,qmesh,main  <BR>  <BR>asel,s,area,,3  <BR>aatt,2,,2  <BR>amesh,all  <BR>alls  <BR>  <BR>asel,s,area,,1  <BR>aatt,3,,2  <BR>amesh,all  <BR>alls  <BR>  <BR>lsel,s,loc,x,x1  <BR>lsel,a,loc,x,x1+w  <BR>dl,all,,ux,0  <BR>alls  <BR>  <BR>lsel,s,loc,y,y1  <BR>dl,all,,uy,0  <BR>alls  <BR>  <BR>acel,,1  <BR>                              <BR>fini  <BR>  <BR>/solu  <BR>antype,static  <BR>  <BR>deltim,0.1,0.05,0.2  <BR>autots,on !使用自动时间步  <BR>pred,on !打开时间步长预测器  <BR>lnsrch,on !打开线性搜索  <BR>nlgeom,on !打开大位移效果  <BR>nropt,full !设定牛顿-拉普森选项  <BR>cnvtol,f,,0.02,0  <BR>  <BR>esel,s,type,,1  <BR>ekill,all  <BR>esel,all  <BR>esel,s,live  <BR>nsle,s  <BR>nsel,inve  <BR>d,all,all,0  <BR>nsel,all  <BR>esel,all  <BR>solve  <BR>  <BR>esel,s,mat,,3  <BR>ekill,all  <BR>  <BR>esel,s,type,,1  <BR>ealive,all  <BR>nsle,s  <BR>ddele,all,all  <BR>esel,all  <BR>esel,s,live  <BR>nsle,s  <BR>nsel,inve  <BR>d,all,all,0  <BR>nsel,all  <BR>esel,all  <BR>  <BR>solve  <BR>fini  <BR>/post1  <BR>etable,if,smisc,1  <BR>etable,jf,smisc,7  <BR>etable,im,smisc,6  <BR>etable,jm,smisc,12  <BR>
 楼主| 发表于 2006-5-3 17:54 | 显示全部楼层
真是谢谢你了,我现在遇到的情况这样的,我们先对四类围岩进行固结灌浆,然后达到三类的标准,在进行喷锚支护,这中间钢筋,喷射混凝土,锚杆,改造后的围岩应该采用那种单元来模拟呢,谢谢指导
 楼主| 发表于 2006-5-3 17:56 | 显示全部楼层
还有接触面的约束&围岩地影响范围应该怎么取
发表于 2006-5-4 07:44 | 显示全部楼层

回复:(94117239)用ansys模拟隧洞的分析实例

围岩可以solid单元<BR>锚杆用杆单元<BR>拱架可用梁单元
 楼主| 发表于 2006-5-4 09:49 | 显示全部楼层
谢谢指导,我在仔细琢磨一下
 楼主| 发表于 2006-5-9 10:49 | 显示全部楼层
<P>锚杆怎么加荷载啊,面荷载怎么加,锚杆的参数怎么取</P>
发表于 2010-7-28 22:24 | 显示全部楼层
非常感谢啊,救了我的命啊!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 07:25 , Processed in 0.075620 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表