|
就是这个啦：
/prep7
/title,analysis of a cftb
!主梁节点
/prep7
z=-3.0124 !(=-6.0248/2)
y=-22.9207 ! (=-45.3653/2-0.2380)
b=0
n,b+1,-84,y,z $ n,b+2,-69,y,z
n,b+3,-57,y,z $ n,b+4,-45,y,z
*do,i,b+5,b+14,1
n,i,-45+9*(i-4),y,z
*enddo
n,b+15,56,y,z
ngen,2,15,all,,,,,5.72
!主梁左侧节点由对称产生，这里不写
!吊杆底部节点放在后面
！拱肋关键点
c=30 !c的值可变
n,c+1,0,1.3454,-0.2130 $ n,c+26,0,-1.3454,0.2130
n,c+2,-3,1.2736,-0.2016 $ n,c+27,-3,-1.4182,0.2245
n,c+3,-6,1.0591,-0.1677 $ n,c+28,-6,-1.6374,0.2592
n,c+4,-9,0.7014,-0.1110 $ n,c+29,-9,-2.0027,0.3170
n,c+5,-12,0.2006,-0.0316 $ n,c+30,-12,-2.5146,0.3981
n,c+6,-15,-0.4438,0.0702 $ n,c+31,-15,-3.1728,0.5023
n,c+7,-18,-1.2311,0.1949 $ n,c+32,-18,-3.9785,0.6298
n,c+8,-21,-2.1621,0.3423 $ n,c+33,-21,-4.9315,0.7807
n,c+9,-24,-3.2372,0.5125 $ n,c+34,-24,-6.0321,0.9549
n,c+10,-27,-4.4564,0.7055 $ n,c+35,-27,-7.2811,1.1526
n,c+11,-30,-5.8205,0.9214 $ n,c+36,-30,-8.6785,1.3738
n,c+12,-33,-7.3287,1.1602 $ n,c+37,-33,-10.2262,1.6188
n,c+13,-36,-8.9825,1.4219 $ n,c+38,-36,-11.9238,1.8876
n,c+14,-39,-10.7820,1.7068 $ n,c+39,-39,-13.7722,2.1802
n,c+15,-42,-12.7303,2.0152 $ n,c+40,-42,-15.7699,2.4964
n,c+16,-45,-14.8188,2.3459 $ n,c+41,-45,-17.9269,2.8379
n,c+17,-48,-17.0570,2.7002 $ n,c+42,-48,-20.2351,3.2033
n,c+18,-51,-19.4426,3.0778 $ n,c+43,-51,-22.6984,3.5932
n,c+19,-54,-21.9766,3.4790 $ n,c+44,-54,-25.3178,4.0079
n,c+20,-57,-24.6590,3.9036 $ n,c+45,-57,-28.0949,4.4475
n,c+21,-60,-27.4903,4.3518 $ n,c+46,-60,-31.0313,4.9123
n,c+22,-63,-30.4717,4.8238 $ n,c+47,-63,-34.1279,5.4025
n,c+23,-66,-33.6029,5.3194 $ n,c+48,-66,-37.3871,5.9185
n,c+24,-69,-36.8825,5.8386 $ n,c+49,-69,-40.8127,6.4608
n,c+25,-71,-39.1646,6.1998
nsel,s,,,31,79,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
!（49＝79+1-31）
nsel,s,,,80,128,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
nsel,s,,,80,128,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
nsel,all
nsel,s,,,32,49,1
nsym,x,146 ,all,,,
nsel,s,,,57,74,1
nsym,x,139 ,all,,,
nsel,s,,,81,98,1
nsym,x,133 ,all,,,
nsel,s,,,106,123,1
nsym,x,126 ,all,,,
nsel,s,,,130,147,1
nsym,x,120 ,all,,,
nsel,s,,,155,172,1
nsym,x,113,all,,,
nsel,all
nplot
!补右拱脚节点
n,286,55,-26.221,4.1509
n,287,55,-26.1038,4.8917
n,288,55,-25.9865,5.6324
n,289,57,-24.6590,3.9036
n,290,57,-24.4245,4.64435
n,291,57,-24.4245,5.3851
！横撑中心节点
nsel,s,,,5
ngen,2,287,all,,,0,10.9969,0
nsel,s,,,5
ngen,2,288,all,,,0,13.9382,0
nsel,s,,,7
ngen,2,287,all,,,0,18.9422,0
nsel,s,,,7
ngen,2,288,all,,,0,21.6896,0
nsel,s,,,9
ngen,2,287,all,,,0,21.5753,0
nsel,s,,,9
ngen,2,288,all,,,0,24.2661,0
nsel,s,,,11
ngen,2,287,all,,,0,18.9422,0
nsel,s,,,11
ngen,2,288,all,,,0,21.6896,0
nsel,s,,,13
ngen,2,287,all,,,0,10.9969,0
nsel,s,,,13
ngen,2,288,all,,,0,13.9382,0
nsel,all
nplot
*do,i,1,11,1
n,301+i,-45+9*(i-1),-22.6827,2.6876
*enddo
*do,i,1,6,1
nsel,s,,,105+3*(i-1)
ngen,2,208-2*(i-1),all,,,0,-0.4981,0.01358
*enddo
*do,i,1,5,1
nsel,s,,,234+3*(i-1)
ngen,2,85-2*(i-1),all,,,0,-0.4981,0.01358
*enddo
nsel,all
nplot
!立柱
!324，2要耦合
!325，3要耦合
!326，17要耦合
!327，18要耦合
nsel,s,,,2,3,1
ngen,2,322,all,,,0,-(1.6-0.238),0
nsel,s,,,17,18,1
ngen,2,309,all,,,0,-(1.6-0.238),0
nsel,s,,,326
ngen,2,2,all,,,0,0,3.131
nsel,s,,,327
ngen,2,2,all,,,0,0,1.196
nsel,all
nplot
!定义单元类型
et,1,beam4
et,2,link10
keyopt,2,2,2
keyopt,2,3,0
et,3,shell63
!定义材料
!拱肋弦杆
MP,EX,1,0.527e11
MP,nuxy,1,0.25
MP,dens,1,3188
!拱肋腹杆，横撑强、弱构件
MP,EX,2,2.06e11
MP,nuxy,2,0.3
MP,dens,2,7850
!拱肋上下平联
MP,EX,3,0.355e11
MP,nuxy,3,0.2
MP,dens,3,2600
!吊杆
MP,EX,4,1.95e11
MP,nuxy,4,0.3
MP,dens,4,7850
!立柱、立柱横梁
MP,EX,5,0.3e11
MP,nuxy,5,0.2
MP,dens,5,2500
!鱼骨梁主梁
MP,EX,6,0.355e11
MP,nuxy,6,0.2
MP,dens,6,4119
!鱼骨梁刺梁
MP,EX,7,3.55e11 !弹性模量输入主梁的十倍
MP,nuxy,7,0.2
MP,dens,7,0.1
!定义实常数（面积除了拱肋弦杆外，其他输入实际面积以和密度对应）
!拱肋弦杆
r,1,0.37,0.03956,0.03956,0.7,0.7
rmore,,0.02357
!拱肋腹杆、横撑强构件
r,2,0.01225,0.000233,0.000233,0.4,0.4
!拱肋横撑弱构件
r,3,0.004825,0.0000233,0.0000233,0.2,0.2
!立柱
r,4,1,0.083333,0.083333,1,1
rmore,,0.1426
!立柱横梁
r,5,1.44,0.1728,0.1728,1.2,1.2
rmore,,0.2957
!鱼骨梁主梁
r,6,6.718,2.061,103.4,2.5,1.6 !注意输入
rmore,0.02357
!鱼骨梁刺梁
r,7,0.01,1030.4,20.016,0.1,0.1 !抗弯惯性矩为主梁的十倍
!吊杆
r,8,0.00535,1e-8
!上下平联板
r,9,0.5，0.5，0.5，0.5，
!主梁
！刺梁
type,1
mat,7
real,7
*do,i,1,15,1
e,i,i+15
*enddo
！主梁
type,1
mat,6
real,6
*do,i,1,14,1
e,i,i+1
!e,e+15,i+16 !三主梁的边梁
*enddo
!拱肋
！弦杆
type,1
mat,1
real,1
*do,i,1,24,1
e,i+30,i+31
e,i+128,i+129
*enddo
*do,i,1,23,1
e,i+55,i+56
e,i+153,i+154
*enddo
*do,i,1,17,1
e,i+177,i+178
e,i+195,i+196
e,i+249,i+250
e,i+267,i+268
*enddo
e,31,178 $ e,129,250
e,56,196 $ e,154,268
e,213,286 $ e,285,288
e,195,289 $ e,267,291
！竖腹杆
type,1
mat,2
real,2
*do,i,1,24,1
e,i+30,i+55
e,i+128,i+153
*enddo
*do,i,1,18,1
e,i+177,i+195
e,i+249,i+267
*enddo
！斜腹杆
type,1
mat,2
real,2
*do,i,1,23,1
e,i+31,i+55
e,i+129,i+153
*enddo
*do,i,1,17,1
e,i+178,i+195
e,i+250,i+267
*enddo
e,178,56 $ e,250,154
e,213,289 $ e,285,291
!横撑
！横撑强构件
type,1
mat,2
real,2
e,292,293 $ e,292,43
e,294,295 $ e,294,37
e,296,297 $ e,296,31
e,298,299 $ e,298,183
e,300,301 $ e,300,189
！横撑弱构件
type,1
mat,2
real,3
e,293,43 $ e,292,68
e,295,37 $ e,294,62
e,297,31 $ e,296,56
e,299,183 $ e,298,201
e,301,189 $ e,300,207
!吊杆
type,2
mat,4
real,8
e,302,318 $ e,303,317 $ e,304,316
e,305,315 $ e,306,314 $ e,307,313
e,308,319 $ e,309,320 $ e,310,321
e,311,322 $ e,312,323
e,95,120 $ e,92,117 $ e,89,114
e,86,111 $ e,83,108 $ e,80,105
e,216,234 $ e,219,237 $ e,222,240
e,225,243 $ e,228,246
!立柱
！立柱横梁
type,1
mat,5
real,5
e,324,326 $ e,326,328
e,325,327 $ e,327,329
！竖立柱
type,1
mat,5
real,4
e,328,54 $ e,329,50
!上下平联板
type,3
mat,3
real,9
！左平联
*do,i,1,24,1
e,i+30,i+31,i+80,i+79
e,i+79,i+80,i+129,i+128
*enddo
*do,i,1,23,1
e,i+55,i+56,i+105,i+104
e,i+104,i+105,i+154,i+153
*enddo
e,31,178,214,80 $ e,80,214,250,129
e,56,196,232,105 $ e,105,232,268,154
！右平联
*do,i,1,17,1
e,i+177,i+178,i+214,i+213
e,i+213,i+214,i+250,i+249
e,i+195,i+196,i+232,i+231
e,i+231,i+232,i+268,i+267
*enddo
e,195,289,290,231 $ e,231,290,291,267
e,213,286,287,249 $ e,249,287,288,285
/eshape,1
eplot
cs,11,0,297,296,298,1,1
csys,11
dsys,11
nsym,z,329,all !前面一共three two nine 个节点
ensym,585,,329,all !前面一共five eight five个单元
eplot
csys,0
dsys,0
/replo
nummrg,all, , , ,low
numcmp,all
*do,i,1,6,1
cp,next,ux,312+i,102+3*i
cp,next,uy,312+i,102+3*i
cp,next,uz,312+i,102+3*i
*enddo
*do,i,1,5,1
cp,next,ux,318+i,231+3*i
cp,next,uy,318+i,231+3*i
cp,next,uz,318+i,231+3*i
*enddo
*do,i,1,11,1
cp,next,ux,301+i,18+i
cp,next,uy,301+i,18+i
cp,next,uz,301+i,18+i
*enddo
*do,i,1,6,1
cp,next,ux,616+i,416+3*i
cp,next,uy,616+i,416+3*i
cp,next,uz,616+i,416+3*i
*enddo
*do,i,1,5,1
cp,next,ux,622+i,545+3*i
cp,next,uy,622+i,545+3*i
cp,next,uz,622+i,545+3*i
*enddo
*do,i,1,11,1
cp,next,ux,605+i,332+i
cp,next,uy,605+i,332+i
cp,next,uz,605+i,332+i
*enddo
*do,i,1,2,1
cp,next,ux,16+i,325+i
cp,next,uy,16+i,325+i
cp,next,uz,16+i,325+i
cp,next,rotx,16+i,325+i
cp,next,roty,16+i,325+i
cp,next,rotz,16+i,325+i
*enddo
*do,i,1,2,1
cp,next,ux,1+i,323+i
cp,next,uy,1+i,323+i
cp,next,uz,1+i,323+i
cp,next,rotx,1+i,323+i
cp,next,roty,1+i,323+i
cp,next,rotz,1+i,323+i
*enddo
*do,i,1,2,1
cp,next,ux,330+i,627+i
cp,next,uy,330+i,627+i
cp,next,uz,330+i,627+i
cp,next,rotx,330+i,627+i
cp,next,roty,330+i,627+i
cp,next,rotz,330+i,627+i
*enddo
/solu
antype,2
modopt,lanb,10
mxpand,10, , ,0
lumpm,0
pstres,0
modopt,lanb,10,0,100, ,off
d,286, , , , 288,2 ,all
d,289, , , , 291,2 ,all
d,600, , , , 602,2 ,all
d,603, , , , 605,2 ,all
d,467, , , , , ,all
d,369, , , , , ,all
d,491, , , , , ,all
d,393, , , , , ,all
d,55, , , , , ,all
d,153, , , , , ,all
d,79, , , , , ,all
d,177, , , , , ,all
d,344,,,,,,ux,uy,uz
d,30,,,,,,ux,uy,uz
d,15,,,,,,ux,uy,uz
d,330,,,,,,uy,uz
d,1,,,,,,uy,uz
d,16,,,,,,uy,uz
solve
finish |
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
