!空間斜拉橋分析
finish
/clear
/PREP7 !進入前處理
et,1,Beam44 !定義梁單元(橋面系)
et,2,Pipe16 !定義管單元(拱)
et,3,Beam44 !定義梁單元(剛臂)
et,4,Pipe16 !定義管單元(拱腳加厚段)
et,5,Link8 !定義索單元(拉索)
mp,dens,1,7850 !定義質量密度
mp,prxy,1,1/3 !定義泊桑比
mp,ex,1,2.07e11 !定義彈性模量
acel,,,9.8 !定義加速度,Z正方向,大小為9.8
r,1,0.1523,0.1979,0.0599,0.70,1.61 !定義梁單元(橋面系)實常數
r,2,1.6,0.02 !定義管單元(拱)實常數
r,3,0.1546,0.0543,0.0543,0.70,0.75 !定義梁單元(剛臂)實常數
r,4,1.6,0.04 !定義管單元(拱腳加厚段)實常數
r,5,1.963e-3 !定義索單元實常數(直徑5cm)
n,1,0.00,-33.82,1.40 !生成橋面體系節點
n,2,1.76,-33.11,1.53
n,3,5.37,-31.34,1.78
n,4,8.74,-29.18,2.01
n,5,11.84,-26.66,2.23
n,6,14.65,-23.81,2.43
n,7,17.13,-20.67,2.60
n,8,19.78,-16.23,2.73
n,9,20.96,-13.66,2.84
n,10,22.27,-9.88,2.96
n,11,23.15,-5.98,3.02
n,12,23.59,-2.00,3.05
n,13,23.65,0.00,3.06 !對稱點
n,14,23.59,2.00,3.05
n,15,23.15,5.98,3.02
n,16,22.27,9.88,2.96
n,17,20.96,13.66,2.84
n,18,19.78,16.23,2.73
n,19,17.13,20.67,2.60
n,20,14.65,23.81,2.43
n,21,11.84,26.66,2.23
n,22,8.74,29.18,2.01
n,23,5.37,31.34,1.78
n,24,1.76,33.11,1.53
n,25,0.00,33.82,1.40
n,26,17.32,-16.64,-2.08 !生成拱節點
n,27,16.82,-16.60,-0.26
n,28,16.07,-16.23,2.53
n,29,13.51,-13.06,11.88
n,30,12.78,-11.51,14.62
n,31,12.11,-9.59,17.12
n,32,11.54,-7.28,19.25
n,33,11.11,-4.58,20.85
n,34,10.88,-1.57,21.74
n,35,10.84,0.00,21.74 !對稱點
n,36,10.88,1.57,21.74
n,37,11.11,4.58,20.85
n,38,11.54,7.28,19.25
n,39,12.11,9.59,17.12
n,40,12.78,11.51,14.62
n,41,13.51,13.06,11.88
n,42,16.07,16.23,2.53
n,43,16.82,16.60,-0.26
n,44,17.32,16.64,-2.08
allsel !全部選擇
!生成梁單元(橋面系)
type,1 !單元代號為1(前面定義梁單元)
mat,1 !材料特性代號為1
real,1 !實常數代號為1
*do,i,1,24,1 !進行循環操作開始
e,i,i+1 根據節點生成單元
*enddo !進行循環操作結束
!生成拱單元
type,2 !單元代號為2(前面定義管單元)
mat,1 !材料特性代號為1
real,2 !實常數代號為2
*do,i,29,40,1 !進行循環操作開始
e,i,i+1 根據節點生成單元
*enddo !進行循環操作結束
!生成剛臂單元
type,3 !單元代號為3(前面定義梁單元)
mat,1 !材料特性代號為1
real,3 !實常數代號為3
e,8,28 根據節點生成單元
e,18,42 根據節點生成單元
!生成拱(加厚段)單元
type,4 !單元代號為4(前面定義管單元)
mat,1 !材料特性代號為1
real,4 !實常數代號為4
e,26,27 !根據節點生成單元
e,27,28 !根據節點生成單元
e,28,29 !根據節點生成單元
e,41,42 !根據節點生成單元
e,42,43 !根據節點生成單元
e,43,44 !根據節點生成單元
!生成拉索單元
type,5 !單元代號為5(前面定義索單元)
mat,1 !材料特性代號為1
real,5 !實常數代號為5
e,3,33 !根據節點生成單元
e,4,32
e,5,31
e,6,30
e,7,29
e,9,31
e,10,32
e,11,33
e,12,34
e,14,36
e,15,37
e,16,38
e,19,41
e,20,40
e,21,39
e,22,38
e,23,37
allsel !全選
/eshape,1.0 !設置單元比例
eplot !繪制單元
d,26,all !約束拱腳(固定)
d,44,all !約束拱腳(固定)
d,1,uz !約束橋面體系(水平滑動)
d,25,uz !約束橋面體系(水平滑動)
!加載
*do,i,1,25,1
f,i,fz,-58000*4 !Z負方向58000N/m
*enddo
*do,i,1,7,1
f,i,mx, 58000*4*1.6*0.707 !X正方向的扭矩
f,i,my, 58000*4*1.6*0.707 !Y正方向的扭矩
*enddo
*do,i,8,18,1
f,i,my, 58000*4*1.6 !Y正方向的扭矩
*enddo
*do,i,19,25,1
f,i,mx, -58000*4*1.6*0.707 !X負方向的扭矩
f,i,my, 58000*4*1.6*0.707 !Y正方向的扭矩
*enddo
allsel !全選
/solu !進入求解器
solve !求解
/post1 !進入后處理
plnsol,s,1 !顯示第一主應力
prrsol !顯示支點反力