|
某连续梁0#块空间分析
好久没发帖子了，把09年做的一个连续梁桥0#块的空间分析的项目分享给大家。
/COM, Structural
FINI
/CLE
/view,1,1,1,1
/Prep7
/TITLE,0# Analysis
!指定主梁单元类型c55混凝土 ------------------------------------------------------
et,1,45
!定义材料属性
mp,ex,1,3.55e10
mp,prxy,1,0.2
mp,dens,1,2600
!建立模型------------------------------------------------------------------------
!设置0#块顶板中心位置为坐标圆点，横桥向为x轴，竖向为y轴，不考虑横坡的影响
!建立-6~6米0#块单元,因为整个模型为对称结构，所以先建立四分之一模型
!建立顶板节点
k,1,0,0,0
k,2,0.5,0,0
k,3,2.775,0,0
k,4,3.575,0,0
k,5,4.575,0,0
k,6,6.875,0,0
k,7,0,-0.8,0
k,8,0.5,-0.8,0
k,9,2.775,-0.8,0
k,10,3.575,-1,0
k,11,4.575,-0.5,0
k,12,6.875,-0.18,0
*do,i,1,5,1
a,0+i,1+i,7+i,6+i
*enddo
!建立底板节点
k,13,0,-6.771,0
k,14,0.5,-6.771,0
k,15,2.775,-6.771,0
k,16,3.575,-6.771,0
k,17,0,-8.572,0
k,18,0.5,-8.572,0
k,19,2.775,-8.572,0
k,20,3.575,-8.572,0
*do,i,1,3,1
a,12+i,13+i,17+i,16+i
*enddo
!建立腹板
k,21,2.775,-5.071,0
k,22,3.575,-5.071,0
k,23,2.775,-6.572,0
k,24,3.575,-6.572,0
a,9,10,22,21
a,21,22,24,23
a,23,24,16,15
!沿Z轴向拉伸截面1.4米
VEXT,all, , ,0,0,1.4,,,,
!建立z=1.4米截面上的人洞关键点
k,107,0,-5.071,1.4
k,108,0.5,-5.071,1.4
k,117,0,-6.572,1.4
k,118,0.5,-6.572,1.4
a,kp(0,-0.8,1.4),kp(0.5,-0.8,1.4),108,107
a,kp(0.5,-0.8,1.4),kp(2.775,-0.8,1.4),kp(2.775,-5.071,1.4),108
a,108,kp(2.775,-5.071,1.4),kp(2.775,-6.572,1.4),118
a,117,118,kp(0.5,-6.771,1.4),kp(0,-6.771,1.4)
a,118,kp(2.775,-6.572,1.4),kp(2.775,-6.771,1.4),kp(0.5,-6.771,1.4)
!将z=1.4米截面拉伸1.2米，建立体
asel,s,loc,z,1.4
vext,all,,,0,0,1.2,,,,
kwpave,kp(2.77,0,2.6)
wprot,0,0,90
wpoff,0,0,-0.2
allsel
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!重新建立一次2.6米处关键点---------------------------------------------------------------------
k,101,0,0,2.6
k,102,0.5,0,2.6
k,103,2.775,0,2.6
k,104,3.575,0,2.6
k,105,4.575,0,2.6
k,106,6.875,0,2.6
k,107,0,-0.8,2.6
k,108,0.5,-0.8,2.6
k,109,2.775,-0.8,2.6
k,110,3.575,-1,2.6
k,111,4.575,-0.5,2.6
k,112,6.875,-0.18,2.6
!建立底板节点
k,113,0,-6.771,2.6
k,114,0.5,-6.771,2.6
k,115,2.775,-6.771,2.6
k,116,3.575,-6.771,2.6
k,117,0,-8.572,2.6
k,118,0.5,-8.572,2.6
k,119,2.775,-8.572,2.6
k,120,3.575,-8.572,2.6
!建立腹板
k,121,2.775,-5.071,2.6
k,122,3.575,-5.071,2.6
k,123,2.775,-6.572,2.6
k,124,3.575,-6.572,2.6
!复制节点-------------------------------------------------------------------------------------------
!Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imove
! Itime：拷贝份数
! Np1,Np2,Ninc：所选关键点
! Dx,Dy,Dz：偏移坐标
! Kinc：每份之间节点号增量
kgen,2,101,102,1,0,0,0.9,900
kgen,2,103,103,1,0.053,0,0.9,900
kgen,2,104,106,1,0,0,0.9,900
kgen,2,107,108,1,0,0.119,0.9,900
kgen,2,109,109,1,0.053,0.119,0.9,900
kgen,2,110,112,1,0,0,0.9,900
kgen,2,121,121,1,0.053,0,0.9,900
kgen,2,122,122,1,0,0,0.9,900
kgen,2,123,123,1,0.053,0,0.9,900
kgen,2,124,124,1,0,0,0.9,900
kgen,2,113,114,1,0,-0.159,0.9,900
kgen,2,115,115,1,0.053,-0.159,0.9,900
kgen,2,116,116,1,0,-0.159,0.9,900
kgen,2,117,118,1,0,0,0.9,900
kgen,2,119,119,1,0.053,0,0.9,900
kgen,2,120,120,1,0,0,0.9,900
allsel
*do,i,1,5,1
v,100+i,100+i+1,100+i+7,100+i+6,1000+i,1000+i+1,1000+i+7,1000+i+6
*enddo
*do,i,13,15,1
v,100+i,100+i+1,100+i+5,100+i+4,1000+i,1000+i+1,1000+i+5,1000+i+4
*enddo
v,100+9,100+10,100+22,100+21,1000+9,1000+10,1000+22,1000+21
v,100+21,100+22,100+24,100+23,1000+21,1000+22,1000+24,1000+23
v,100+23,100+24,100+16,100+15,1000+23,1000+24,1000+16,1000+15
allsel
*do,i,1,5,1
v,100+i,100+i+1,100+i+7,100+i+6,1000+i,1000+i+1,1000+i+7,1000+i+6
*enddo
*do,i,13,15,1
v,100+i,100+i+1,100+i+5,100+i+4,1000+i,1000+i+1,1000+i+5,1000+i+4
*enddo
v,100+9,100+10,100+22,100+21,1000+9,1000+10,1000+22,1000+21
v,100+21,100+22,100+24,100+23,1000+21,1000+22,1000+24,1000+23
v,100+23,100+24,100+16,100+15,1000+23,1000+24,1000+16,1000+15
!复制生成第0#块和1#块交界面的点----------------------------------------------------------
kgen,2,101,102,1,0,0,3.4,1900
kgen,2,103,103,1,0.2,0,3.4,1900
kgen,2,104,106,1,0,0,3.4,1900
kgen,2,107,108,1,0,0.45,3.4,1900
!kgen,2,109,109,1,0.2,0.45,3.4,1900
kgen,2,109,109,1,0.2,-0.2,3.4,1900
kgen,2,110,112,1,0,0,3.4,1900
kgen,2,121,121,1,0.2,-0.4,3.4,1900
kgen,2,122,122,1,0,-0.4,3.4,1900
kgen,2,123,123,1,0.2,-0.099,3.4,1900
kgen,2,124,124,1,0,-0.099,3.4,1900
kgen,2,113,114,1,0,-0.4,3.4,1900
kgen,2,115,115,1,0.2,-0.4,3.4,1900
kgen,2,116,116,1,0,-0.4,3.4,1900
kgen,2,117,118,1,0,0.247,3.4,1900
kgen,2,119,119,1,0.2,0.247,3.4,1900
kgen,2,120,120,1,0,0.247,3.4,1900
!生成体
allsel
*do,i,1,5,1
v,1000+i,1000+i+1,1000+i+7,1000+i+6,2000+i,2000+i+1,2000+i+7,2000+i+6
*enddo
*do,i,13,15,1
v,1000+i,1000+i+1,1000+i+5,1000+i+4,2000+i,2000+i+1,2000+i+5,2000+i+4
*enddo
v,1000+9,1000+10,1000+22,1000+21,2000+9,2000+10,2000+22,2000+21
v,1000+21,1000+22,1000+24,1000+23,2000+21,2000+22,2000+24,2000+23
v,1000+23,1000+24,1000+16,1000+15,2000+23,2000+24,2000+16,2000+15
!重新建立节点------------------------------------------------
k,3001,0,0,6
k,3002,0.5,0,6
k,3003,0.975,0,6
k,3004,2.475,0,6
k,3005,2.975,0,6
k,3006,3.575,0,6
k,3007,4.575,0,6
k,3008,6.875,0,6
k,3009,0,-0.35,6
k,3010,0.5,-0.35,6
k,3011,0.975,-0.35,6
k,3012,2.475,-0.65,6
k,3013,2.975,-1,6
k,3014,3.575,-1,6
k,3015,4.575,-0.5,6
k,3016,6.875,-0.18,6
!建立底板节点
k,3021,0,-7.171,6
k,3022,0.5,-7.171,6
k,3023,2.475,-7.171,6
k,3024,2.975,-7.171,6
k,3025,3.575,-7.171,6
k,3026,0,-8.324,6
k,3027,0.5,-8.324,6
k,3028,2.475,-8.324,6
k,3029,2.975,-8.324,6
k,3030,3.575,-8.324,6
!建立腹板
k,3031,2.975,-5.471,6
k,3032,3.575,-5.471,6
k,3033,2.975,-6.671,6
k,3034,3.575,-6.671,6
!复制生成1#和2#交界面的关键点---------------------------------------------------------------------
kgen,2,3001,3016,1,0,0,2,1000
kgen,2,3021,3025,1,0,0.158,2,1000
kgen,2,3031,3034,1,0,0.158,2,1000
kgen,2,3026,3030,1,0,0.196,2,1000
!生成体
allsel
*do,i,1,7,1
v,3000+i,3000+i+1,3000+i+9,3000+i+8,4000+i,4000+i+1,4000+i+9,4000+i+8
*enddo
*do,i,21,24,1
v,3000+i,3000+i+1,3000+i+6,3000+i+5,4000+i,4000+i+1,4000+i+6,4000+i+5
*enddo
v,3000+13,3000+14,3000+32,3000+31,4000+13,4000+14,4000+32,4000+31
v,3000+31,3000+32,3000+34,3000+33,4000+31,4000+32,4000+34,4000+33
v,3000+33,3000+34,3000+25,3000+24,4000+33,4000+34,4000+25,4000+24
v,3023,3024,3033,4023,4024,4033
!复制生成2#和3#交界面的关键点---------------------------------------------------------------------
kgen,2,4001,4016,1,0,0,3,1000
kgen,2,4021,4025,1,0,0.232,3,1000
kgen,2,4031,4034,1,0,0.232,3,1000
kgen,2,4026,4030,1,0,0.287,3,1000
!生成体
allsel
*do,i,1,7,1
v,4000+i,4000+i+1,4000+i+9,4000+i+8,5000+i,5000+i+1,5000+i+9,5000+i+8
*enddo
*do,i,21,24,1
v,4000+i,4000+i+1,4000+i+6,4000+i+5,5000+i,5000+i+1,5000+i+6,5000+i+5
*enddo
v,4000+13,4000+14,4000+32,4000+31,5000+13,5000+14,5000+32,5000+31
v,4000+31,4000+32,4000+34,4000+33,5000+31,5000+32,5000+34,5000+33
v,4000+33,4000+34,4000+25,4000+24,5000+33,5000+34,5000+25,5000+24
v,4023,4024,4033,5023,5024,5033
!复制生成3#和4#交界面的关键点---------------------------------------------------------------------
kgen,2,5001,5016,1,0,0,3,1000
kgen,2,5021,5025,1,0,0.227,3,1000
kgen,2,5031,5034,1,0,0.227,3,1000
kgen,2,5025,5030,1,0,0.28,3,1000
!生成体
allsel
*do,i,1,7,1
v,5000+i,5000+i+1,5000+i+9,5000+i+8,6000+i,6000+i+1,6000+i+9,6000+i+8
*enddo
*do,i,21,24,1
v,5000+i,5000+i+1,5000+i+6,5000+i+5,6000+i,6000+i+1,6000+i+6,6000+i+5
*enddo
v,5000+13,5000+14,5000+32,5000+31,6000+13,6000+14,6000+32,6000+31
v,5000+31,5000+32,5000+34,5000+33,6000+31,6000+32,6000+34,6000+33
v,5000+33,5000+34,5000+25,5000+24,6000+33,6000+34,6000+25,6000+24
v,5023,5024,5033,6023,6024,6033
!切割，并建立体，以有利于分网和建模------------------------------------------------------------
kwpave,kp(0.975,0,8)
vsbw,all
wpoff,0,0,1.5
vsbw,all
!删除部分体单元
vsel,s,loc,x,0.5,2.98
vsel,r,loc,z,3.5,6
vsel,r,loc,y,0,-1
VDELE,all, , ,1
!重新建立体
v,kp(2.48,0,3.5),kp(2.48,-1.05,3.5),kp(2.83,-1.05,3.5),kp(2.83,0,3.5),kp(2.47,0,6),kp(2.47,-0.65,6),kp(2.97,-1,6),kp(2.97,0,6)
v,kp(0.975,0,3.5),kp(0.975,-1.05,3.5),kp(2.48,-1.05,3.5),kp(2.48,0,3.5),kp(0.75,0,6),kp(0.975,-0.35,6),kp(2.47,-0.65,6),kp(2.47,0,6)
v,kp(2.48,-6.93,3.5),kp(2.83,-6.93,3.5),kp(2.83,-6.57,3.5),kp(2.48,-7.17,6),kp(2.98,-7.17,6),kp(2.98,-6.67,6)
v,1002,1008,100,99,2002,2008,3011,3003
v,163,1015,1023,151,65,72
!从中间分割腹板以便布设竖向预应力------------------------------------------------------------------
kwpave,kp(2.97,0,14)
vsel,s,loc,z,6,14
wpoff,0,0,0.3
vsbw,all
allsel
vsel,s,loc,z,0,2.6
wpoff,0,0,-0.1
vsbw,all
allsel
!分割顶板与腹板--------------------------------------------------------------------------------------------
kwpave,kp(2.77,-1.3,0)
wprot,0,90,0
wpoff,0,0,0.5
vsbw,all
wpoff,0,0,-0.5
allsel
!将需要布设钢束的地方网格加密
!lsel,s,loc,z,0
!lsel,r,loc,x,3.575
!lsel,r,loc,y,-1.3,-5.071
!LESIZE,all,0.2, , , , , , ,1
!切分出钢束位置--------------------------------------------------------
!切分出W01的位置，距离上缘，2.872米
wpoff,0,0,1.371
vsel,s,loc,z,0,3.5
vsbw,all
!切分出W02的位置，距离上缘，2.671米
wpoff,0,0,0.201
vsel,s,loc,z,0,3.5
vsbw,all
allsel
vglue,all
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!分网格------------------------------------------------------------------------------------------------------
vatt,1,1,1
esize,0.5
vsel,s,loc,z,0,1.4
vsweep,all
lsel,s,loc,z,1.4,8
lsel,r,loc,y,0
lsel,r,loc,x,6.875
LESIZE,ALL,0.5, , , ,1, , ,1,
LSEL,S, , , 172
LESIZE,all, , ,2, , , , ,1
lsel,s,loc,z,8,1.4
lsel,s,loc,z,8,14
lsel,r,loc,y,0
lsel,r,loc,x,6.875
LESIZE,ALL,1, , , ,1, , ,1,
lsel,s,loc,z,8,14
vsel,s,loc,z,1.4,2.6
vsweep,all
vsel,s,loc,z,2.6,3.5
vsweep,all
vsel,s,loc,z,3.5,6
vsweep,all
vsel,s,loc,z,6,8
vsweep,all
esize,1
vsel,s,loc,z,8,11
vsweep,all
et,2,45
mp,ex,2,3.55e10
mp,prxy,2,0.2
mp,dens,2,2600
vsel,s,loc,z,11,13
vatt,2,2,2
vsweep,all
!自由划分变截面处,将顶部与腹板切割位置移动下来0.5米后，就没有分网不良的问题，所以不用自由分网了--------
esize,0.25
VSEL, , , , 162
vatt,1,1,1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
!对称体————————————————————————————————————————————————————
allsel
VSYMM,X,all, , , ,0,0
allsel
VSYMM,Z,all, , , ,0,0
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!加载==============================================================================================================
hxy=78e4
hxs=0.7
!如果施加竖向预应力，要沿着腹板中心线切分------------------
!施加横向预应力,19.5吨/束
!定立节点选择集
*dim,AA,array,47
nsel,s,loc,y,-0.09
nsel,r,loc,x,6.875
aa(1)=3267
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fx,-hxy*hxs
*enddo
*do,i,3,20
f,aa(2*i+1),fx,-hxy*hxs
*enddo
*do,i,42,46
f,aa(i),fx,-hxy*hxs
*enddo
!*dim,AA,array,47
nsel,s,loc,y,-0.09
nsel,r,loc,x,-6.875
!nsel,u,loc,x,-13
aa(1)=4848
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fx,hxy*hxs
*enddo
*do,i,3,20
f,aa(2*i+1),fx,hxy*hxs
*enddo
*do,i,42,46
f,aa(i),fx,hxy*hxs
*enddo
allsel
!加1#and2#块预应力钢束===================================================================================================
!指定钢束单元类型 为link8单元,W01和W02为19-15.2，面积为2.64e-3m2------------------------------------------------------
et,3,link8
!定义材料属性
mp,ex,3,1.95e11
mp,prxy,3,0.3
mp,dens,3,7800
!直径15.24mm,标准强度1860MPa，单根张拉控制吨位193.9kN
areagjx=1.4e-4
*dim,gjarea,array,2
gjarea(1)=19*areagjx
!有效张拉吨位------------------------------------------------------------
F1=193.9e3 !1#钢绞线张拉控制吨位
F2=193.9e3 !2#钢绞线张拉控制吨位
E=1.95e11 !钢绞线弹模
L=1 !每根钢绞线长1m
!每根钢绞线单元的伸长量-----------------------------------------------------
detL1=F1*L/(E*areagjx) !1#钢绞线单元的伸长量
detL2=F2*L/(E*areagjx) !2#钢绞线单元的伸长量
!每根钢绞线单元的的拉应变-------------------------------------------------
istrn1=detL1/(detL1+L) !1#钢绞线单元的的拉应变
istrn2=detL2/(detL2+L) !2#钢绞线单元的的拉应变
!定义折减系数------------------------------------------------
zj1=0.85
zj2=0.85
!定义实常数------------------------------------------------
r,2,gjarea(1),zj1*istrn1 !定义1#钢束的实常数
r,3,gjarea(1),zj2*istrn2 !定义2#钢束的实常数
!建立W01,定制右侧钢束节点选择集，建立右侧钢束单元---------------------------------------------
*dim,xx,array,17
*dim,yy,array,1
*dim,zz,array,17
xx(1)=3.2015,3.1883,3.175,3.175,3.175,3.175,3.175,3.175,3.175,3.175
xx(11)=3.175,3.175,3.175,3.175,3.175,3.1883,3.2015
yy(1)=-2.872
zz(1)=-3.5,-3.05,-2.6,-2.2,-1.8,-1.4,-0.93333,-0.46667,0,0.46667
zz(11)=0.93333,1.4,1.8,2.2,2.6,3.05,3.5
TYPE, 3
MAT, 3
REAL, 3
ESYS, 0
SECNUM,
TSHAP,LINE
*do,i,1,16,1
e,node(xx(i),yy(1),zz(i)),node(xx(i+1),yy(1),zz(i+1))
*enddo
e,node(3.2750,-7.2130,-8.0000),node(xx(1),yy(1),zz(1))
e,node(xx(17),yy(1),zz(17)),node(3.2750,-7.2130,8.0000)
!建立左侧钢束单元
*do,i,1,16,1
e,node(-xx(i),yy(1),zz(i)),node(-xx(i+1),yy(1),zz(i+1))
*enddo
e,node(-3.2750,-7.2130,-8.0000),node(-xx(1),yy(1),zz(1))
e,node(-xx(17),yy(1),zz(17)),node(-3.2750,-7.2130,8.0000)
!建立W02,定制右侧钢束节点选择集，建立右侧钢束单元---------------------------------------------
!布设腹板钢束单元
*dim,xx,array,24
*dim,yy,array,24
*dim,zz,array,24
xx(1)=3.275,3.2603,3.2309,3.2015,3.1883,3.175,3.175,3.175,3.175,3.175
xx(11)=3.175,3.175,3.175,3.175,3.175,3.175,3.175,3.175,3.1883,3.2015
xx(21)=3.2162,3.2309,3.2603,3.275
yy(1)=-2.6772,-2.6267,-2.4993,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372
yy(11)=-2.372,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372,-2.372
yy(21)=-2.4357,-2.4993,-2.6267,-2.6772
zz(1)=-6.5,-5.5,-4.5,-3.5,-3.05,-2.6,-2.2,-1.8,-1.4,-0.93333
zz(11)=-0.46667,0,0.46667,0.93333,1.4,1.8,2.2,2.6,3.05,3.5
zz(21)=4,4.5,5.5,6.5
TYPE, 3
MAT, 3
REAL, 3
ESYS, 0
SECNUM,
TSHAP,LINE
*do,i,1,23,1
e,node(xx(i),yy(i),zz(i)),node(xx(i+1),yy(i+1),zz(i+1))
*enddo
e,node(3.2750,-7.1343,-11.0000),node(xx(1),yy(1),zz(1))
e,node(xx(24),yy(24),zz(24)),node(3.2750,-7.1343,11.0000)
!建立左侧钢束单元
*do,i,1,23,1
e,node(-xx(i),yy(i),zz(i)),node(-xx(i+1),yy(i+1),zz(i+1))
*enddo
e,node(-3.2750,-7.1343,-11.0000),node(-xx(1),yy(1),zz(1))
e,node(-xx(24),yy(24),zz(24)),node(-3.2750,-7.1343,11.0000)
!布设顶板钢束单元--------------------------------------------
!说明：顶板为直束，将其两根钢束的作用力合并为一根，作用在腹板顶缘下方0.7米处，采用集中力加载
!本节点力为2倍顶板钢束荷载的值
f,node(3.275,-0.7,-8),fz,1.6*19*193.9e3
f,node(-3.275,-0.7,-8),fz,1.6*19*193.9e3
f,node(3.275,-0.7,8),fz,-1.6*19*193.9e3
f,node(-3.275,-0.7,8),fz,-1.6*19*193.9e3
f,node(3.275,-0.7,-11),fz,1.6*19*193.9e3
f,node(-3.275,-0.7,-11),fz,1.6*19*193.9e3
f,node(3.275,-0.7,11),fz,-1.6*19*193.9e3
f,node(-3.275,-0.7,11),fz,-1.6*19*193.9e3
!钢束设置完毕-----------------------------------------------------------------
!钢束设置完毕-----------------------------------------------------------------
sxy=113.56e4 !竖向预应力
ss=0.7 !竖向预应力有效系数
!施加竖向预应力,113.56吨/2根，现在为两根合成一根---------------------------
!右侧顶板------------------------
*dim,AA,array,47
nsel,s,loc,y,0
nsel,r,loc,x,3.174,3.30
aa(1)=3353
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fy,-ss*sxy
*enddo
*do,i,3,20
f,aa(2*i+1),fy,-ss*sxy
*enddo
*do,i,42,46
f,aa(i),fy,-ss*sxy
*enddo
!左侧顶板---------------------------
!*dim,AA,array,47
nsel,s,loc,y,0
nsel,r,loc,x,-3.174,-3.30
aa(1)=5328
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fy,-ss*sxy
*enddo
*do,i,3,20
f,aa(2*i+1),fy,-ss*sxy
*enddo
*do,i,42,46
f,aa(i),fy,-ss*sxy
*enddo
!右侧底板-----------------------------------------------------
csys,0
allsel
nsel,s,loc,x,3.174,3.37
cm,sy1,node
nsel,s,loc,x,3.174,3.37
nsel,r,loc,y,-8.572
cm,sy11,node
!定义左侧局部坐标系，旋转坐标---------
!CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
k,1500,3.275,-7.562,14
k,2000,3.275,-7.6553,13
k,3000,3.275,-17.655,13.933
cskp,100,0,1500,2000,3000
!------------------------------------------------
csys,0
cmsel,s,sy1
nsel,r,loc,z,2.7,15
!在局部坐标选择左侧下边缘节点
csys,100
nsel,r,loc,y,-0.2,0.15
cm,sy12,node
!定义右侧局部坐标系，旋转坐标---------
csys,0
!CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
k,4000,3.275,-7.562,-14
k,5000,3.275,-7.6553,-13
k,6000,3.275,-17.655,-13.933
cskp,200,0,4000,5000,6000
!------------------------------------------------
csys,0
cmsel,s,sy1
nsel,r,loc,z,-2.7,-15
!在局部坐标选择又侧下边缘节点
csys,200
nsel,r,loc,y,-0.2,0.15
cm,sy13,node
!将点集汇总
cmsel,s,sy11
cmsel,a,sy12
cmsel,a,sy13
cm,sy1,node
FLST,5,2,1,ORDE,2
FITEM,5,2118
FITEM,5,9556
NSEL,U, , ,P51X
csys,0
!施加垂直向上的力---------
aa(1)=3380
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fy,ss*sxy
*enddo
*do,i,3,20
f,aa(2*i),fy,ss*sxy
*enddo
*do,i,42,46
f,aa(i),fy,ss*sxy
*enddo
!左侧底板------------------------------------------------------------
csys,0
allsel
nsel,s,loc,x,-3.174,-3.37
cm,sy2,node
nsel,s,loc,x,-3.174,-3.37
nsel,r,loc,y,-8.572
cm,sy21,node
csys,0
cmsel,s,sy2
nsel,r,loc,z,2.7,15
!在局部坐标选择左侧下边缘节点
csys,100
nsel,r,loc,y,-0.2,0.15
cm,sy22,node
!------------------------------------------------
csys,0
cmsel,s,sy2
nsel,r,loc,z,-2.7,-15
!在局部坐标选择又侧下边缘节点
csys,200
nsel,r,loc,y,-0.2,0.15
cm,sy23,node
!将点集汇总
cmsel,s,sy21
cmsel,a,sy22
cmsel,a,sy23
FLST,5,2,1,ORDE,2
FITEM,5,6156
FITEM,5,12779
NSEL,U, , ,P51X
csys,0
!施加垂直向上的力---------
aa(1)=5340
*do,i,1,46
aa(i+1)=nnear(aa(i))
nsel,u,,,aa(i)
*enddo
allsel
*do,i,2,6
f,aa(i),fy,ss*sxy
*enddo
*do,i,3,20
f,aa(2*i),fy,ss*sxy
*enddo
*do,i,42,46
f,aa(i),fy,ss*sxy
*enddo
csys,0
来自圈子: 桥梁结构分析有限元 |
-
o#块空间分析
-
|
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
