交叉支撑梁的屈曲分析
@王培齐
FINISH
/CLEAR
/FILENAME,Schwedler
/TITLE,The buckling analyze of HFB keep section
Bf=0.0801 !上翼缘宽度
H1=0.0403 !上翼缘高度
TF1=0.0034 !上翼缘厚度
H=0.2634 !腹板截面高度
TW=0.0056 !腹板厚度
TF2=0.0054 !下翼缘厚度
bf2=Bf !下翼缘宽度
L=1.8 !跨度
L1=0.054
force=-1
/PREP7
ET,1,solid65
mp,ex,1,3.25e10
mp,prxy,1,0.2
mp,ex,2,2.06e11
mp,prxy,2,0.3
ET,2,BEAM188
MP,EX,3,2.06e11 !弹性模量
MP,PRXY,3,0.3 !泊松比
MP,DENS,3,7850
SECTYPE,1,BEAM,CSOLID
SECDATA,0.003 !杆的实常数
*do,I,0,tw,0
k,,i-L1,,bf2/2
k,,i-L1,,tw/2
k,,i-L1,,-tw/2
k,,i-L1,,-bf2/2
k,,i-L1,tf2,bf2/2
k,,i-L1,tf2,tw/2
k,,i-L1,tf2,-tw/2
k,,i-L1,tf2,-bf2/2
k,,i-L1,tf2+h,bf/2
k,,i-L1,tf2+h,(bf)/2-tf1
k,,i-L1,tf2+h,tw/2
k,,i-L1,tf2+h,-tw/2
k,,i-L1,tf2+h,-(bf)/2+tf1
k,,i-L1,tf2+h,-bf/2
k,,i-L1,tf2+h+tf1,bf/2
k,,i-L1,tf2+h+tf1,(bf)/2-tf1
k,,i-L1,tf2+h+tf1,tw/2
k,,i-L1,tf2+h+tf1,-tw/2
k,,i-L1,tf2+h+tf1,-(bf)/2+tf1
k,,i-L1,tf2+h+tf1,-bf/2
k,,i-L1,tf2+h+h1-tf1,bf/2
k,,i-L1,tf2+h+h1-tf1,(bf)/2-tf1
k,,i-L1,tf2+h+h1-tf1,tw/2
k,,i-L1,tf2+h+h1-tf1,-tw/2
k,,i-L1,tf2+h+h1-tf1,-(bf)/2+tf1
k,,i-L1,tf2+h+h1-tf1,-bf/2
k,,i-L1,tf2+h+h1,bf/2
k,,i-L1,tf2+h+h1,(bf)/2-tf1
k,,i-L1,tf2+h+h1,tw/2
k,,i-L1,tf2+h+h1,-tw/2
k,,i-L1,tf2+h+h1,-(bf)/2+tf1
k,,i-L1,tf2+h+h1,-bf/2
*enddo
a,1,2,6,5
a,2,3,7,6
a,3,4,8,7
a,6,7,12,11
*do,j,0,12,6
*do,i,9,13,1
a,i+j,i+j+1,i+j+7,i+j+6
*enddo
*enddo
k,999,-L1
k,1000
k,1001,tw,
k,1002,L/2-tw/2
k,1003,L/2+tw/2
k,1004,L-tw
k,1005,L
k,1006,L+L1
*do,i,999,1005,1
l,i,i+1
*enddo
vdrag,all,,,,,,51,52,53,54,55,56
asel,s,loc,x,L-0.001,L+0.001
vdrag,all,,,,,,57
allsel,all
v,43,75,74,42,36,68,67,35
v,51,83,73,41,40,72,70,38
v,107,139,138,106,100,132,131,99
v,115,147,137,105,104,136,134,102
v,171,203,202,170,164,196,195,163
v,179,211,201,169,168,200,198,166
vsel,s,loc,y,tf2+h+h1/2
vsel,u,loc,z,(bf)/2-tf1/2
vsel,u,loc,z,-(bf)/2+tf1/2
vatt,1,,1
allsel,all
vsel,u,loc,y,tf2+h+h1/2
vsel,a,loc,z,(bf)/2-tf1/2
vsel,a,loc,z,-(bf)/2+tf1/2
vatt,2,,1
allsel,all
vgen,2,all,,,,,-0.300,,, !偏移0.3m
K,,0,tf2,-tw/2
K,,0,tf2,-0.3+tw/2
K,,0,tf2+h,-0.3+tw/2
K,,0,tf2+h,-tw/2
K,,0,tf2+h/2,-0.15
K,,L/2+tw/2,tf2,-tw/2
K,,L/2+tw/2,tf2,-0.3+tw/2
K,,L/2+tw/2,tf2+h,-0.3+tw/2
K,,L/2+tw/2,tf2+h,-tw/2
K,,L/2+tw/2,tf2+h/2,-0.15
K,,L,tf2,-tw/2
K,,L,tf2,-0.3+tw/2
K,,L,tf2+h,-0.3+tw/2
K,,L,tf2+h,-tw/2
K,,L,tf2+h/2,-0.15
L,517,513
L,517,514
L,517,515
L,517,516
L,522,518
L,522,519
L,522,520
L,522,521
L,527,523
L,527,524
L,527,525
L,527,526
LSEL,S,LINE,,1304,1315,1
LATT,3,3,2
LESIZE,ALL,,,10
LMESH,ALL
vsel,s,loc,y,tf2+h+h1/2
vsel,u,loc,z,(bf)/2-tf1/2
vsel,u,loc,z,-(bf)/2+tf1/2
vsel,u,loc,z,-0.3+(bf)/2-tf1/2
vsel,u,loc,z,-0.3-(bf)/2+tf1/2
vatt,1,,1
allsel,all
vsel,u,loc,y,tf2+h+h1/2
vsel,a,loc,z,(bf)/2-tf1/2
vsel,a,loc,z,-(bf)/2+tf1/2
vsel,a,loc,z,-0.3+(bf)/2-tf1/2
vsel,a,loc,z,-0.3-(bf)/2+tf1/2
vatt,2,,1
allsel,all
esize,0.027
mshape,0
mshkey,1
vmeshh,all
nsel,s,loc,x,0
nsel,s,loc,x,L/2+tw/2
nsel,s,loc,x,L+0.000001,L-0.000001
allsel,all
cerig,32,170
cerig,22,5739
cerig,2,146
cerig,12,5682
cerig,73,2796
cerig,63,8354
cerig,43,2791
cerig,53,8335
cerig,114,5246
cerig,104,10804
cerig,84,5241
cerig,94,10785
allsel,all
nsel,s,loc,x,0-0.001,0+0.001
nsel,a,loc,x,L-0.001,L+0.001
nsel,u,loc,y,0
nsel,u,loc,y,tf2+h
nsel,u,loc,y,tf2+h+tf1
nsel,u,loc,y,tf2+h+h1/2
nsel,u,loc,y,tf2+h+tf1
nsel,u,loc,y,tf2+h+h1-tf1
nsel,u,loc,y,tf2+h+h1
nsel,r,loc,z,-tw/2,tw/2
d,all,uz
nsel,r,loc,y,tf2
d,all,uy
nsel,r,loc,x,0
d,all,ux
allsel,all
allsel,all
nsel,s,loc,x,0-0.001,0+0.001
nsel,a,loc,x,L-0.001,L+0.001
nsel,u,loc,y,0
nsel,u,loc,y,tf2+h
nsel,u,loc,y,tf2+h+tf1
nsel,u,loc,y,tf2+h+h1/2
nsel,u,loc,y,tf2+h+tf1
nsel,u,loc,y,tf2+h+h1-tf1
nsel,u,loc,y,tf2+h+h1
nsel,r,loc,z,-0.3-tw/2,-0.3+tw/2
d,all,uz
nsel,r,loc,y,tf2
d,all,uy
nsel,r,loc,x,0
d,all,ux
allsel,all
!施加荷载
nsel,s,loc,y,tf2+h+h1
nsel,r,loc,x,L/2
F,all,FY,-1/16
ALLSEL,ALL
/solu
pstres,on
solve
finish
/solu
antype,buckle
bucopt,lanb,5
mxpand,5
outres,all,all
solve
finish
问题1，交叉支撑用什么单元比较好。
问题2，本命令流中，交叉支撑采用梁单元，由于约束不对导致屈曲分析中跨中交叉支撑没有协同工作，这种情况应该怎样定义梁单元与实体单元的接触
邹正刚 
