隧道初期支护与二次衬砌接触模拟求助（附命令流）

xyghzzj · 发表于 2009-6-1 22:10:48

最近做了一个公路隧道开挖的二维数值模拟，用梁单元模拟初期支护和二次衬砌，初期支护和二次衬砌采用接触单元，结果初期支护和二次衬砌接触不上，求高人指点迷津，在下不省感激！

xyghzzj · 发表于 2009-6-1 22:11:37

! 一建模

! ANSYS环境设定
finish
/clear
/filename,sd

! 1. 参数定义
DefineParameter

! 2. 实体建模
GenModel

! (6) 生成部件与组件
GenCom

! (7) 改变梁单元圆弧段的法线方向
ReverseArc

! 3. 划分网格
GenGrid

! 4. 生成二次衬砌
profile2

! 5. 生成接触单元
GenContact

! 二加载与求解
Calculate

! 三计算结果分析
!Post

xyghzzj · 发表于 2009-6-1 22:11:48

! 1. 参数定义

! 应力释放系数
! 开挖应力释放系数
Factor_1=0.1
! 支护应力释放系数
Factor_2=0.3
! 二次衬砌应力释放系数
Factor_3=0.6

! 所需计算时步总数
Number_Totle=6
! 目标面Z坐标
Loc_Z=0

! (1)几何参数设定

! 设置角度单位
! 设置角度单位为弧度
*afun,rad
! 弧度转角度因子
PI=acos(-1)
RadToDeg=180/PI
DegToRad=PI/180
! 设置角度单位为度
*afun,deg

g=10

! 隧道衬砌内轮廓拱部圆弧圆心
x_1=0
y_1=0
! 隧道衬砌内轮廓拱部圆弧半径
R_1=5.75
! 隧道衬砌内轮廓拱部对应圆心角
st_1=180

! 隧道衬砌内轮廓边墙部位圆弧圆心
x_2=-2.300
y_2=0
! 隧道衬砌内轮廓边墙部位圆弧半径
R_2=8.05
! 隧道衬砌内轮廓边墙部位对应圆心角
st_2_s=0
st_2_e=10.4311
st_2=st_2_s+st_2_e

! 隧道衬砌内轮廓拱脚部位圆弧圆心
x_3=4.388
y_3=-1.231
! 隧道衬砌内轮廓拱脚部位圆弧半径
R_3=1.25
! 隧道衬砌内轮廓拱脚部位对应圆心角
st_3_s=10.4311
st_3_e=71.3914
st_3=st_3_e-st_3_s

! 隧道衬砌内轮廓仰拱圆心
x_4=0
y_4=11.800
! 隧道衬砌内轮廓仰拱半径
R_4=15.000
! 隧道衬砌内轮廓仰拱对应圆心角
st_4=37.2172

! 临时仰拱圆心
x_h=0
y_h=30
! 临时仰拱半径
R_h=30
! 临时仰拱圆心角
st_h=15

! 初期支护厚度
T_IniLin=0.26
! 初期支护面积
A_IniLin=1*T_IniLin

! 二次衬砌厚度
T_SecLin=0.55
! 二次衬砌面积
A_SecLin=1*T_SecLin

! 隧道初期支护拱部圆弧半径
R_1=R_1+T_SecLin+T_IniLin/2
! 隧道初期支护边墙圆弧半径
R_2=R_2+T_SecLin+T_IniLin/2
! 隧道初期支护拱脚圆弧半径
R_3=R_3+T_SecLin+T_IniLin/2
! 隧道初期支护仰拱圆弧半径
R_4=R_4+T_SecLin+T_IniLin/2

! 隧道单洞开挖最大宽度
B=2*(R_2+x_2)
! 隧道单洞开挖最大高度
H=(R_1+y_1)+(R_4-y_4)

! 计算区域:
! 考虑工程的需要和有限元离散误差以及计算误差，一般选计算范围均不小于3~4倍洞室直径。
! 地层右侧计算边界与隧道内轮廓右边界距离
W_right=60
! 地层下侧计算边界与隧道内轮廓拱底距离
H_bot=45
! 隧道埋深(地层上侧计算边界与隧道内轮廓拱顶距离)
H_top=32

! 隧道周围网格加密范围
l_grid=10
! 上下台阶分界线半径
R_bench=30

! (2)材料参数设定

! 初期支护材料参数
! 衬砌类型:Z5q
! 围岩类别:5级围岩
! 喷射混凝土:260mm厚C25钢筋砼

! 围岩材料参数：5级围岩
! 围岩容重,N/m^3
r_Rock=18.5e3
! 容重转换为密度,单位由N/m^3转换为国际单位制kg/m^3
DENS_Rock=r_Rock/g
! 围岩弹性模量，Pa
E_Rock=1.5e9
! 围岩泊松比
U_Rock=0.4
! 围岩粘聚力，Pa
C_Rock=0.125e6
! 围岩内摩擦角，度
FI_Rock=23.5

! 初期支护材料参数：C25砼
! 初期支护容重,kN/m^3
r_IniLin_C25=22e3
! 容重转换为密度,单位由N/m^3转换为国际单位制kg/m^3
DENS_IniLin_C25=r_IniLin_C25/g
! 初期支护弹性模量，Pa
E_IniLin_C25=23e9
! 初期支护泊松比
U_IniLin_C25=0.2

! 二次衬砌材料参数：C30钢筋砼
! 二次衬砌容重,N/m^3
r_SecLin_C30=25e3
! 容重转换为密度,单位由N/m^3转换为国际单位制kg/m^3
DENS_SecLin_C30=r_SecLin_C30/g
! 二次衬砌弹性模量，Pa
E_SecLin_C30=31e9
! 二次衬砌泊松比
U_SecLin_C30=0.2

! 约束弹簧刚度
k_spring=1e3

! (3)定义单元类型

! 进入前处理器
/prep7
! 定义围岩单元
et,1,plane42
! 定义为平面应变问题
keyopt,1,3,2

! 定义梁单元
et,2,beam3
! 打开输出内力选项
keyopt,2,6,1

! 定义目标单元
et,3,TARGE169
! 定义接触单元
et,4,conta171
! 闭合间隙或减少初始穿透
!KEYOPT,4,5,3
! 壳的厚度影响
KEYOPT,4,11,1

! (4)定义实常数

! 定义初期支护的面积、惯性矩和梁高
r,1,T_IniLin,T_IniLin**3/12,T_IniLin
! 定义二次衬砌的面积、惯性矩和梁高
r,2,T_SecLin,T_SecLin**3/12,T_SecLin
! 定义接触单元实常数，法向接触刚度因子1，最大渗透范围0.001
r,3,,,1,0.001,,
! 定义接触单元实常数，最大接触摩擦0.5MPa
RMORE,,,0.5e6

! (5)定义材料属性

! 剩余围岩的材料属性
mp,ex,1,E_Rock
mp,prxy,1,U_Rock
mp,dens,1,DENS_Rock
tb,dp,1,,,
tbmodif,1,1,C_Rock
tbmodif,1,2,FI_Rock

! 围岩材料性能是一样的，但是在这里以不同的材料编号加以区分“剩余围岩”和“被挖去围岩”
! 的目的是为了后面方便地选择被挖去的有效围岩单元。

! 隧道内轮廓以内被挖去围岩的材料属性
mp,ex,2,E_Rock
mp,prxy,2,U_Rock
mp,dens,2,DENS_Rock+0.1
tb,dp,2,,,
tbmodif,1,1,C_Rock
tbmodif,1,2,FI_Rock

! 初期支护材料属性
mp,ex,3,E_IniLin_C25
mp,prxy,3,U_IniLin_C25
mp,dens,3,DENS_IniLin_C25

! 二次衬砌材料属性
mp,ex,4,E_SecLin_C30
mp,prxy,4,U_SecLin_C30
mp,dens,4,DENS_SecLin_C30

! 接触单元材料属性
mp,mu,5,0.2

xyghzzj · 发表于 2009-6-1 22:12:21

! (1) 划分拱部与侧墙接触面单元

! 设置单元类型
type,3
! 设置实常数
real,3
! 设置材料属性
mat,5

! 生成目标面接触单元节点部件
cmsel,s,ASSM_IniLin_line
lplot
nsll,s,1
nplot
esln,s,1
CM,CM_Con1_IniLin_node,node
! 生成目标面接触单元
esurf,all,top
! 改变接触面接触单元外法向
!esurf,,reverse
allsel,all

! 检查目标面接触单元外法向
!ESEL,S,TYPE,,3
!/PSYMB,ESYS,1
!allsel,all

! 设置单元类型
type,4
! 设置实常数
real,3
! 设置材料属性
mat,5

! 生成接触面接触单元节点部件
cmsel,s,CM_SecLin_line
lplot
nsll,s,1
nplot
esln,s,1
CM,CM_Con1_SecLin_node,node
! 生成接触面接触单元
esurf,all,Top
! 改变接触面接触单元外法向
esurf,,reverse
allsel,all

! 检查接触面接触单元外法向
!ESEL,S,TYPE,,4
!/PSYMB,ESYS,1

! 施加拱顶弹簧约束
psprng,NODE(0,R_1-T_SecLin/2-T_IniLin/2,0),tran,k_spring,1
psprng,NODE(0,R_1-T_SecLin/2-T_IniLin/2,0),tran,k_spring,,-1
! 施加拱底弹簧约束
psprng,NODE(0,y_4-R_4+T_SecLin/2+T_IniLin/2,0),tran,k_spring,1
psprng,NODE(0,y_4-R_4+T_SecLin/2+T_IniLin/2,0),tran,k_spring,,1

zbhualala · 发表于 2009-6-2 15:18:38

谢谢楼主，辛苦啦

10202 · 发表于 2010-3-28 06:28:45

我也需要做这个，但必须要用接触单元？初衬和二衬必然不产生拉应力...何不1个plane421个beam3...

xiazuguo2008 · 发表于 2010-6-9 15:16:41

好资料啊！

qingw8 · 发表于 2010-7-27 16:18:09

楼主这么辛苦，发了那么多，咋没有高手帮忙解答啊

zhuchh0720 · 发表于 2010-10-10 08:55:54

我是初学者，太好了，正需要这个

江beyond · 发表于 2011-1-4 12:17:44

太棒了！！

zhangmeng123 · 发表于 2011-1-5 12:07:01

太好了，有没有ABABQUS的

apoleo · 发表于 2011-1-8 11:04:03

我也是隧道专业的

mobeimobei · 发表于 2011-4-22 21:39:45

二次衬砌不产生拉应力？

banana5117 · 发表于 2011-4-23 10:38:24

lz看来回复，估计对仿真论坛彻底失望了！

zhuhaha · 发表于 2011-6-7 20:47:13

我看很多人也没有怎么施加接触啊

zdpjll · 发表于 2011-8-6 22:00:18

接触没算出来啊

小浩瀚宇宙 · 发表于 2011-8-14 20:24:53

回复 1# xyghzzj

很不错啊

残梦断忆 · 发表于 2011-8-26 23:14:49

我做的初支换个二衬是用link1单元模拟的，因为初支和二衬之间几乎没有剪应力。其实用link10更好，但是总是不收敛我也不知道什么原因，你的比我高级多了。你看看我的命令流：
finish
/clear
/prep7

*set,shuxiang,106502.4 !竖向荷载
*set,zuoce,23963.04 !水平荷载
*set,youce,-23963.04

*set,uxishu,350e6 !水平弹性抗力系数
*set,vxishu,350e6 !竖向弹性抗力系数

!定义单元类型、实常数及材料性质
et,1,beam3
et,2,link1
r,1,0.12,0.12*0.12*0.12/12,0.12 !初支截面特征
mp,ex,1,2.95e10 !初支弹性模量
mp,dens,1,2300 !初支密度
mp,prxy,1,0.167  !初支泊松比

r,2,0.3,0.3*0.3*0.3/12,0.3 !二衬截面特征
r , 3 , 0.412200194 , 0.005836377 , 0.412200194
r , 4 , 0.439410241 , 0.007070161 , 0.439410241
r , 5 , 0.483074787 , 0.009394245 , 0.483074787
r , 6 , 0.542887965 , 0.013333661 , 0.542887965
r , 7 , 0.618992878 , 0.019764039 , 0.618992878
r , 8 , 0.711171137 , 0.029973753 , 0.711171137
r , 9 , 0.819535363 , 0.045869272 , 0.819535363
r , 10 , 0.943818868 , 0.070062353 , 0.943818868
r , 11 , 1.015203881 , 0.087192386 , 1.015203881
r , 12 , 0.8968 , 0.060104301 , 0.8968

mp,ex,2,3.23e9 !二衬弹性模量
mp,dens,2,2200 !二衬密度
mp,prxy,2,0.167  !二衬泊松比

mp,ex,3,2.95e10 !接触弹性模量
mp,prxy,3,0.167  !接触泊松比
!接触单元实常数
r , 3 , 0.118244315 , 0
r , 4 , 0.117132524 , 0
r , 5 , 0.121586893 , 0
r , 6 , 0.123461705 , 0
r , 7 , 0.125601396 , 0
r , 8 , 0.127971493 , 0
r , 9 , 0.130561719 , 0
r , 10 , 0.133378805 , 0
r , 11 , 0.136404282 , 0
r , 12 , 0.136854871 , 0
r , 13 , 0.137930673 , 0
r , 14 , 0.142507634 , 0
r , 15 , 0.147716851 , 0
r , 16 , 0.153403131 , 0
r , 17 , 0.159414197 , 0
r , 18 , 0.165767935 , 0
r , 19 , 0.172421049 , 0
r , 20 , 0.179240648 , 0
r , 21 , 0.18614813 , 0
r , 22 , 0.19318184 , 0
r , 23 , 0.200273639 , 0
r , 24 , 0.207347138 , 0
r , 25 , 0.214399086 , 0
r , 26 , 0.221348273 , 0
r , 27 , 0.228180527 , 0
r , 28 , 0.234954277 , 0
r , 29 , 0.244017905 , 0
r , 30 , 0.257277042 , 0
r , 31 , 0.270519955 , 0
r , 32 , 0.281456855 , 0
r , 33 , 0.291122681 , 0
r , 34 , 0.299459716 , 0
r , 35 , 0.30637926 , 0
r , 36 , 0.311815102 , 0
r , 37 , 0.315751941 , 0
r , 38 , 0.318070316 , 0
r , 39 , 0.318820343 , 0
r , 40 , 0.31797671 , 0
r , 41 , 0.315556029 , 0
r , 42 , 0.311564001 , 0
r , 43 , 0.306046787 , 0
r , 44 , 0.299055291 , 0
r , 45 , 0.290618109 , 0
r , 46 , 0.280883616 , 0
r , 47 , 0.269941916 , 0
r , 48 , 0.257821 , 0
r , 49 , 0.244799002 , 0
r , 50 , 0.230969871 , 0
r , 51 , 0.216550983 , 0
r , 52 , 0.20182712 , 0
r , 53 , 0.187079906 , 0
r , 54 , 0.171812416 , 0
r , 55 , 0.157769066 , 0
r , 56 , 0.148973731 , 0
r , 57 , 0.143676514 , 0
r , 58 , 0.138974739 , 0
r , 59 , 0.134895351 , 0
r , 60 , 0.131534123 , 0
r , 61 , 0.128937119 , 0
r , 62 , 0.127189664 , 0
r , 63 , 0.126353794 , 0
r , 64 , 0.126442437 , 0
r , 65 , 0.127416946 , 0
r , 66 , 0.129288082 , 0
r , 67 , 0.132045281 , 0
r , 68 , 0.135513479 , 0
r , 69 , 0.139682231 , 0
r , 70 , 0.144507528 , 0
r , 71 , 0.14984035 , 0
r , 72 , 0.160737676 , 0
r , 73 , 0.171207559 , 0
r , 74 , 0.175946315 , 0
r , 75 , 0.180720122 , 0
r , 76 , 0.185448114 , 0
r , 77 , 0.190199397 , 0
r , 78 , 0.194970748 , 0
r , 79 , 0.199623076 , 0
r , 80 , 0.214308955 , 0
r , 81 , 0.261335834 , 0
r , 82 , 2.012357112 , 0
r , 83 , 1.907137798 , 0

!初支坐标
k , 1 , -4.06 , -0.82
k , 2 , -4.0827 , -0.518
k , 3 , -4.1053 , -0.216
k , 4 , -4.128 , 8.60E-02
k , 5 , -4.1507 , 0.388
k , 6 , -4.1734 , 0.69
k , 7 , -4.196 , 0.992
k , 8 , -4.2187 , 1.294
k , 9 , -4.2414 , 1.596
k , 10 , -4.264 , 1.898
k , 11 , -4.2867 , 2.2
k , 12 , -4.3037 , 2.4866
k , 13 , -4.31 , 2.7736
k , 14 , -4.3056 , 3.0606
k , 15 , -4.2906 , 3.3473
k , 16 , -4.265 , 3.6333
k , 17 , -4.2288 , 3.9181
k , 18 , -4.1821 , 4.2013
k , 19 , -4.1249 , 4.4827
k , 20 , -4.0573 , 4.7617
k , 21 , -3.9794 , 5.038
k , 22 , -3.8913 , 5.3112
k , 23 , -3.7931 , 5.581
k , 24 , -3.685 , 5.847
k , 25 , -3.5671 , 6.1087
k , 26 , -3.4396 , 6.3659
k , 27 , -3.3027 , 6.6182
k , 28 , -3.1564 , 6.8653
k , 29 , -2.99 , 7.099
k , 30 , -2.8051 , 7.3186
k , 31 , -2.6042 , 7.5236
k , 32 , -2.3883 , 7.7128
k , 33 , -2.1588 , 7.8851
k , 34 , -1.9168 , 8.0395
k , 35 , -1.6639 , 8.1752
k , 36 , -1.4014 , 8.2914
k , 37 , -1.1309 , 8.3873
k , 38 , -0.85391 , 8.4626
k , 39 , -0.57201 , 8.5166
k , 40 , -0.28683 , 8.5491
k , 41 , 0 , 8.56
k , 42 , 0.28683 , 8.5491
k , 43 , 0.57201 , 8.5166
k , 44 , 0.85391 , 8.4626
k , 45 , 1.1309 , 8.3873
k , 46 , 1.4014 , 8.2914
k , 47 , 1.6639 , 8.1752
k , 48 , 1.9168 , 8.0395
k , 49 , 2.1588 , 7.8851
k , 50 , 2.3883 , 7.7128
k , 51 , 2.6042 , 7.5236
k , 52 , 2.8051 , 7.3186
k , 53 , 2.99 , 7.099
k , 54 , 3.1564 , 6.8653
k , 55 , 3.3027 , 6.6182
k , 56 , 3.4396 , 6.3659
k , 57 , 3.5671 , 6.1087
k , 58 , 3.685 , 5.847
k , 59 , 3.7931 , 5.581
k , 60 , 3.8913 , 5.3112
k , 61 , 3.9794 , 5.038
k , 62 , 4.0573 , 4.7617
k , 63 , 4.1249 , 4.4827
k , 64 , 4.1821 , 4.2013
k , 65 , 4.2288 , 3.9181
k , 66 , 4.265 , 3.6333
k , 67 , 4.2906 , 3.3473
k , 68 , 4.3056 , 3.0606
k , 69 , 4.31 , 2.7736
k , 70 , 4.3037 , 2.4866
k , 71 , 4.2867 , 2.2
k , 72 , 4.264 , 1.898
k , 73 , 4.2414 , 1.596
k , 74 , 4.2187 , 1.294
k , 75 , 4.196 , 0.992
k , 76 , 4.1734 , 0.69
k , 77 , 4.1507 , 0.388
k , 78 , 4.128 , 8.60E-02
k , 79 , 4.1053 , -0.216
k , 80 , 4.0827 , -0.518
k , 81 , 4.06 , -0.82

!二衬坐标
k , 82 , -2.37E-16 , -1.8971
k , 83 , -0.33338 , -1.8885
k , 84 , -0.66588 , -1.8629
k , 85 , -0.99662 , -1.8202
k , 86 , -1.3247 , -1.7606
k , 87 , -1.6494 , -1.6843
k , 88 , -1.9697 , -1.5914
k , 89 , -2.2848 , -1.4823
k , 90 , -2.5939 , -1.3571
k , 91 , -2.8962 , -1.2162
k , 92 , -3.1908 , -1.06
k , 93 , -3.4239 , -0.81591
k , 94 , -3.5505 , -0.50855
k , 95 , -3.6372 , -0.17906
k , 96 , -3.7165 , 0.15376
k , 97 , -3.788 , 0.48956
k , 98 , -3.8516 , 0.82796
k , 99 , -3.9073 , 1.1685
k , 100 , -3.955 , 1.5111
k , 101 , -3.9945 , 1.855
k , 102 , -4.0259 , 2.2
k , 103 , -4.0428 , 2.4729
k , 104 , -4.0498 , 2.7463
k , 105 , -4.0468 , 3.0198
k , 106 , -4.0337 , 3.293
k , 107 , -4.0107 , 3.5655
k , 108 , -3.9778 , 3.8369
k , 109 , -3.9349 , 4.107
k , 110 , -3.8822 , 4.3754
k , 111 , -3.8197 , 4.6416
k , 112 , -3.7476 , 4.9054
k , 113 , -3.6658 , 5.1664
k , 114 , -3.5746 , 5.4242
k , 115 , -3.474 , 5.6785
k , 116 , -3.3641 , 5.9289
k , 117 , -3.2452 , 6.1752
k , 118 , -3.1174 , **169
k , 119 , -2.9808 , 6.6538
k , 120 , -2.8291 , 6.8812
k , 121 , -2.6572 , 7.0938
k , 122 , -2.4694 , 7.2925
k , 123 , -2.2667 , 7.4761
k , 124 , -2.0505 , 7.6434
k , 125 , -1.8219 , 7.7935
k , 126 , -1.5825 , 7.9254
k , 127 , -1.3335 , 8.0384
k , 128 , -1.0765 , 8.1318
k , 129 , -0.8131 , 8.2051
k , 130 , -0.54479 , 8.2577
k , 131 , -0.27321 , 8.2894
k , 132 , 0 , 8.3
k , 133 , 0.27321 , 8.2894
k , 134 , 0.54479 , 8.2577
k , 135 , 0.8131 , 8.2051
k , 136 , 1.0765 , 8.1318
k , 137 , 1.3335 , 8.0384
k , 138 , 1.5825 , 7.9254
k , 139 , 1.8219 , 7.7935
k , 140 , 2.0505 , 7.6434
k , 141 , 2.2667 , 7.4761
k , 142 , 2.4694 , 7.2925
k , 143 , 2.6572 , 7.0938
k , 144 , 2.8291 , 6.8812
k , 145 , 2.9808 , 6.6538
k , 146 , 3.1174 , **169
k , 147 , 3.2452 , 6.1752
k , 148 , 3.3641 , 5.9289
k , 149 , 3.474 , 5.6785
k , 150 , 3.5746 , 5.4242
k , 151 , 3.6658 , 5.1664
k , 152 , 3.7476 , 4.9054
k , 153 , 3.8197 , 4.6416
k , 154 , 3.8822 , 4.3754
k , 155 , 3.9349 , 4.107
k , 156 , 3.9778 , 3.8369
k , 157 , 4.0107 , 3.5655
k , 158 , 4.0337 , 3.293
k , 159 , 4.0468 , 3.0198
k , 160 , 4.0498 , 2.7463
k , 161 , 4.0428 , 2.4729
k , 162 , 4.0259 , 2.2
k , 163 , 3.9945 , 1.855
k , 164 , 3.955 , 1.5111
k , 165 , 3.9073 , 1.1685
k , 166 , 3.8516 , 0.82796
k , 167 , 3.788 , 0.48956
k , 168 , 3.7165 , 0.15376
k , 169 , 3.6372 , -0.17906
k , 170 , 3.5505 , -0.50855
k , 171 , 3.4239 , -0.81591
k , 172 , 3.1908 , -1.06
k , 173 , 2.8962 , -1.2162
k , 174 , 2.5939 , -1.3571
k , 175 , 2.2848 , -1.4823
k , 176 , 1.9697 , -1.5914
k , 177 , 1.6494 , -1.6843
k , 178 , 1.3247 , -1.7606
k , 179 , 0.99662 , -1.8202
k , 180 , 0.66588 , -1.8629
k , 181 , 0.33338 , -1.8885

!创建初支线
*do,i,1,80
l,i,i+1
*enddo

!创建二衬线
*do,i,82,180
l,i,i+1
*enddo
l,181,82

!创建连接线
*do,i,1,81
l,i,i+91
*enddo

!选择初衬并定义属性
lsel,s,,,1,80,1
latt,1,1,1

!选择二衬并定义属性
*do,i,3,12
lsel,s,,,103-i
latt,2,i,1
*enddo
*do,i,3,12
lsel,s,,,158+i
latt,2,i,1
*enddo
lsel,s,,,101,160,1
latt,2,2,1
lsel,s,,,171,180,1
latt,2,2,1
lsel,s,,,81,90,1
latt,2,2,1
!选择接触并定义属性
*do,i,3,83,1
lsel,s,,,178+i
latt,3,i,2
*enddo
allsel,all
lesize,all,,,1             !定义各个线所划分的单元个数
lmesh,all                   !对所有线进行单元划分

!加弹簧
!y方向
*do,i,83,92
a=350000000
k=a*abs(nx(i+1)-nx(i-1))/2
psprng,i,tran,k, , -0.3,
*enddo
*do,i,172,180
a=350000000
k=a*abs(nx(i+1)-nx(i-1))/2
psprng,i,tran,k, , -0.3,
*enddo
a=350000000
k=a*abs(nx(83)-nx(181))/2
psprng,82,tran,k, , -0.3,

a=350000000
k=a*abs(nx(180)-nx(82))/2
psprng,181,tran,k, , -0.3,

!x方向

*do,i,83,92
a=350000000
k=a*abs(ny(i+1)-ny(i-1))/2
psprng,i,tran,k, -0.3, ,
*enddo
*do,i,172,180
a=350000000
k=a*abs(ny(i+1)-ny(i-1))/2
psprng,i,tran,k, 0.3, ,
*enddo
a=350000000
k=a*abs(ny(180)-ny(82))/2
psprng,181,tran,k, 0.3, ,

*do,i,2,18
a=350000000
k=a*abs(ny(i+1)-ny(i-1))/2
psprng,i,tran,k, -0.3, ,
*enddo
*do,i,64,80
a=350000000
k=a*abs(ny(i+1)-ny(i-1))/2
psprng,i,tran,k, 0.3, ,
*enddo
a=350000000
k=a*abs(ny(2)-ny(1))
psprng,1,tran,k, -0.3, ,
a=350000000
k=a*abs(ny(81)-ny(80))
psprng,81,tran,k, 0.3, ,
allsel
finish

!进入求解层，施加荷载，定义荷载步等
/solu

allsel
fcum,add,,  !一定要，使荷载能叠加

!竖向荷载
*do,i,13,68
a=-shuxiang*0.5*abs(nx(i)-nx(i+1))
f,i+1,fy,a
f,i,fy,a
*enddo

!水平荷载
*do,i,82,91
a=zuoce*0.5*abs(ny(i)-ny(i+1))
f,i+1,fx,a
f,i,fx,a
*enddo
*do,i,1,40
a=zuoce*0.5*abs(ny(i)-ny(i+1))
f,i+1,fx,a
f,i,fx,a
*enddo
*do,i,41,80
a=youce*0.5*abs(ny(i)-ny(i+1))
f,i+1,fx,a
f,i,fx,a
*enddo
*do,i,172,180
a=youce*0.5*abs(ny(i)-ny(i+1))
f,i+1,fx,a
f,i,fx,a
*enddo
a=youce*0.5*abs(ny(181)-ny(82))
f,82,fx,a
f,181,fx,a

!施加重力
acel,,9.8

solve !求解
finish

!进入后处理
/post1
!显示变形图
pldisp,1

!仅显示单元类型1
esel,r,type,,1

!显示弯矩图
etable,mi,smisc,6
etable,mj,smisc,12
plls,mi,mj,-1

!显示轴力图
etable,fi,smisc,1
etable,fj,smisc,7
plls,fi,fj,1

!更新单元表数据
etable,refl

!打印单元表数据
pretab,mi,mj,fi,fj

guol1989 · 发表于 2011-11-29 21:14:19

先下看看，顶一下了

zjwcq · 发表于 2012-8-5 14:16:48

:hug::hug::hug::hug:

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[ansys隧洞分析专区] 隧道初期支护与二次衬砌接触模拟求助（附命令流）

本帖子中包含更多资源

相关帖子