急需求助！如何在杆系模型中考虑剪切变形的影响？

以上是我建的模型命令流，还请大家指教，感激不尽！

FINISH
/CLEAR,ALL
/prep7
cyl4,0,0,0.0405,,0.0445
cyl4,0,0,0.0405
cyl4,0,0,0.022,,0.024
asel,all
APLOT
/PNUM,AREA,1
/REPLOT
*SET,es,2.06e11
*SET,ec,3.25e10
!*SET,fp,2.40e8
*SET,fy,3.00E8
*SET,fu,4.50e8
!*SET,fre,0.6e8
et,1,plane82
ET,2,PLANE82
et,3,plane82
ET,4,PLANE82
MP,EX,1,es
MP,PRXY,1,0.28
TB,MISO,1,1,6,
TBTEMP,0
TBPT ,, 0.0,0.0
TBPT ,, 0.8*fy/es , 0.8*fy,
TBPT ,, 1.2*fy/es,fy
TBPT ,, 12*fy/es, fy
TBPT ,, 120*fy/es ,fu
TBPT ,, 300*fy/es ,fu
MP,DENS,2,DENS2
MP,EX,2,ec
MP,PRXY,2,0.2
TB,MISO,2,1,30,
TBTEMP,0
TBPT,, 0 , 0
TBPT,, 0.000100000 , 3.25E+06
TBPT,, 0.000363291 , 9.27E+06
TBPT,, 0.000726581 , 1.76E+07
TBPT,, 0.001089872 , 2.49E+07
TBPT,, 0.001453163 , 3.12E+07
TBPT,, 0.001816453 , 3.66E+07
TBPT,, 0.002179744 , 4.10E+07
TBPT,, 0.002543034 , 4.44E+07
TBPT,, 0.002906325 , 4.69E+07
TBPT,, 0.003269616 , 4.83E+07
TBPT,, 0.003632906 , 4.88E+07
TBPT,, 0.005449359 , 5.04E+07
TBPT,, 0.007265813 , 5.15E+07
TBPT,, 0.009082266 , 5.25E+07
TBPT,, 0.010898719 , 5.33E+07
TBPT,, 0.012715172 , 5.40E+07
TBPT,, 0.014531625 , 5.46E+07
TBPT,, 0.016348078 , 5.52E+07
TBPT,, 0.018164532 , 5.57E+07
TBPT,, 0.019980985 , 5.62E+07
TBPT,, 0.021797438 , 5.66E+07
TBPT,, 0.025430344 , 5.74E+07
TBPT,, 0.029063251 , 5.81E+07
TBPT,, 0.032696157 , 5.88E+07
TBPT,, 0.036329063 , 5.94E+07
TBPT,, 0.054493595 , 6.17E+07
TBPT,, 0.072658127 , 6.35E+07
TBPT,, 0.090822658 , 6.50E+07
TBPT,, 0.10898719 , 6.62E+07
SAVE,ALL

mat,1
TYPE,1
SECNUM,1
esize,0.008,,
amesh,1 !!!!!在1#几何面上生成面单元
sectype,1,beam,mesh
secoffset,cent
secwrite,zzgg,sect,,1 !!!!!柱肢截面
secread,'zzgg','sect',,mesh
mat,2
TYPE,2
SECNUM,2 !定义截面号
esize,0.01,,
amesh,2
ESEL,S,TYPE,,2
sectype,2,beam,mesh
secoffset,cent
secwrite,hnt,sect,,2 !!!!!混凝土截面
secread,'hnt','sect',,mesh
mat,1
TYPE,3
SECNUM,3
esize,0.008,,
amesh,3
ESEL,S,TYPE,,3
sectype,3,beam,mesh
secoffset,cent
secwrite,qzh,sect,,3 !!!!!强轴截面
secread,'qzh','sect',,mesh
asel,all
SAVE,,,,ALL
finish
/clear,all
!!!!!!!!!!!!!以上为建立截面单元
finish
/clear,all
/filename,2.4CHP1-0.248
/prep7
ET,1,BEAM188
ET,2,BEAM188
!keyopt,1,6,1
!KEYOPT,2,6,2
!!!!!!!!!!!!注意这里改变长度
L=2.4
F0=L/1000
!!!!!!!!!!!!
*SET,dens1,7850
*SET,dens2,2500
*SET,es,2.06e11
*SET,ec,3.25e10
!*SET,Fp,2.4E8
*SET,fy,3.00e8
*SET,fu,4.50e8
!*SET,fre,0.6e8
MP,dens,1,den1
MP,EX,1,es
MP,PRXY,1,0.28
TB,MISO,1,1,6,
TBTEMP,0
TBPT ,, 0.0,0.0
TBPT ,, 0.8*fy/es , 0.8*fy,
TBPT ,, 1.2*fy/es,fy
TBPT ,, 12*fy/es, fy
TBPT ,, 120*fy/es ,fu
TBPT ,, 300*fy/es ,fu
MP,DENS,2,DENS2
MP,EX,2,ec
MP,PRXY,2,0.2
TB,MISO,2,1,30,
TBTEMP,0
TBPT,, 0 , 0
TBPT,, 0.000100000 , 3.25E+06
TBPT,, 0.000363291 , 9.27E+06
TBPT,, 0.000726581 , 1.76E+07
TBPT,, 0.001089872 , 2.49E+07
TBPT,, 0.001453163 , 3.12E+07
TBPT,, 0.001816453 , 3.66E+07
TBPT,, 0.002179744 , 4.10E+07
TBPT,, 0.002543034 , 4.44E+07
TBPT,, 0.002906325 , 4.69E+07
TBPT,, 0.003269616 , 4.83E+07
TBPT,, 0.003632906 , 4.88E+07
TBPT,, 0.005449359 , 5.04E+07
TBPT,, 0.007265813 , 5.15E+07
TBPT,, 0.009082266 , 5.25E+07
TBPT,, 0.010898719 , 5.33E+07
TBPT,, 0.012715172 , 5.40E+07
TBPT,, 0.014531625 , 5.46E+07
TBPT,, 0.016348078 , 5.52E+07
TBPT,, 0.018164532 , 5.57E+07
TBPT,, 0.019980985 , 5.62E+07
TBPT,, 0.021797438 , 5.66E+07
TBPT,, 0.025430344 , 5.74E+07
TBPT,, 0.029063251 , 5.81E+07
TBPT,, 0.032696157 , 5.88E+07
TBPT,, 0.036329063 , 5.94E+07
TBPT,, 0.054493595 , 6.17E+07
TBPT,, 0.072658127 , 6.35E+07
TBPT,, 0.090822658 , 6.50E+07
TBPT,, 0.10898719 , 6.62E+07
!钢管
*dim,aa,array,4
*dim,bb,array,35
*dim,cc,array,4
*dim,dd,array,7
aa(1)=0
aa(2)=0
aa(3)=0.320
aa(4)=0.320
bb(1)=0
bb(2)=0.060/2
bb(3)=0.060
*do, i,4,33,1
bb(i)=0.060+(i-3)*0.38 /5
*enddo
bb(34)=2.4-0.060/2
bb(35)=2.4
cc(1)=0
cc(2)=0.320
cc(3)=0.320
cc(4)=0
dd(1)=0.060
dd(2)=0.060+0.38
dd(3)=0.060+0.38*2
dd(4)=0.060+0.38*3
dd(5)=0.060+0.38*4
dd(6)=0.060+0.38*5
dd(7)=0.060+0.38*6
*do,j,1,4
*do,i,1,35,1
n,i+j*35-35,aa(j),bb(i),cc(j)
*enddo
*enddo

type,1
mat,1
SECNUM,1
SECTYPE,1,BEAM,MESH
SECOFFSET,CENT,,,
secread,'zzgg','sect',,MESH
*do,i,1,34,1
e,i,i+1
*enddo
*do,i,36,69,1
e,i,i+1
*enddo
*do,i,71,104,1
e,i,i+1
*enddo
*do,i,106,139,1
e,i,i+1
*enddo !柱肢钢管单元连接
TYPE,2
MAT,2
secnum,2
SECTYPE,2,BEAM,MESH,
SECOFFSET,CENT,,,
secread,'hnt','sect',,MESH
*do,i,1,34,1
e,i,i+1
*enddo
*do,i,36,69,1
e,i,i+1
*enddo
*do,i,71,104,1
e,i,i+1
*enddo
*do,i,106,139,1
e,i,i+1
*enddo !柱肢混凝土单元连接

*do,j,1,7
*do,i,1,4,1
n,200+i+10*(j-1),0,dd(j),i*0.320/5
n,300+i+10*(j-1),i*0.320/5,dd(j),0.320
n,400+i+10*(j-1),0.320,dd(j),0.320-i*0.320/5
n,500+i+10*(j-1),0.320-i*0.320/5,dd(j),0
*enddo
*enddo !!!强弱轴节点
*do,j,1,6
*do,i,1,4,1
n,600+i+10*(j-1),0.320-i*0.320/5,0.060+i*0.38/5+0.38*(j-1),0.320
n,700+i+10*(j-1),0.320-i*0.320/5,0.060+i*0.38/5+0.38*(j-1),0
n,800+i+10*(j-1),0,0.060+i*0.38/5+0.38*(j-1),0.320-i*0.320/5
n,900+i+10*(j-1),0.320,0.060+i*0.38/5+0.38*(j-1),0.320-i*0.320/5
*enddo
*enddo !!!强弱轴斜缀管节点
type,1
mat,1
SECNUM,3
SECTYPE,3,BEAM,MESH
SECOFFSET,CENT,,,
secread,'qzh','sect',,MESH
*do,h,1,4,1
*do,j,1,7,1
*do,i,201,203,1
e,i+10*(j-1)+(h-1)*100,i+1+10*(j-1)+(h-1)*100
*enddo
*enddo
*enddo
*do,i,1,7,1
e,3+(i-1)*5,201+(i-1)*10
e,204+(i-1)*10,38+(i-1)*5
e,38+(i-1)*5,301+(i-1)*10
e,304+(i-1)*10,73+(i-1)*5
e,73+(i-1)*5,401+(i-1)*10
e,404+(i-1)*10,108+(i-1)*5
e,108+(i-1)*5,501+(i-1)*10
e,504+(i-1)*10,3+(i-1)*5
*enddo
!!!强弱轴平缀管单元
*do,h,1,4,1
*do,j,1,6,1
*do,i,601,603,1
e,i+10*(j-1)+(h-1)*100,i+1+10*(j-1)+(h-1)*100
*enddo
*enddo
*enddo

*do,i,1,6,1
e,38+(i-1)*5,801+(i-1)*10
e,804+(i-1)*10,8+(i-1)*5
e,73+(i-1)*5,601+(i-1)*10
e,604+(i-1)*10,43+(i-1)*5
e,73+(i-1)*5,901+(i-1)*10
e,904+(i-1)*10,113+(i-1)*5
e,108+(i-1)*5,701+(i-1)*10
e,704+(i-1)*10,8+(i-1)*5
*enddo !!!强轴斜缀管单元
!改成偏压，间距为0.02m
*do,i,1,15
n,10100+i,0.0,0.0,0.020*i
n,10200+i,0.020*i,0.0,0.320
n,10300+i,0.320,0.0,0.020*i
n,10400+i,0.020*i,0.0,0.0 !构建下加载平面点四个边的
*enddo
*do,j,1,15
*do,i,1,15
n,20000+i+j*100,0.020*i,0.0,0.020*j !补充中间点 划分10*10分
*enddo
*enddo
*do,i,1,15
n,11100+i,0.0,2.4,0.020*i
n,11200+i,0.020*i,2.4,0.320
n,11300+i,0.320,2.4,0.020*i
n,11400+i,0.020*i,2.4,0.0 !构建上加载平面点四个边的
*enddo
*do,j,1,15
*do,i,1,15
n,30000+i+j*100,0.020*i,2.4,0.020*j !补充中间点 划分10*10分
*enddo
*enddo
et,5,beam188
MP,dens,5,0
MP,EX,5,2.1e15
MP,PRXY,5,0.0001
type,5
mat,5
SECNUM,5
sectype,5,beam,rect,duanban !截面号4 矩形端板
secdata,0.09,0.09
*do,i,1,14
e,10100+i,10101+i
e,10200+i,10201+i
e,10300+i,10301+i
e,10400+i,10401+i !下端板四条边连接成单元
*enddo
*do,i,1,15
e,10100+i,20001+i*100
*enddo
*do,i,1,15
e,10200+i,21500+i
*enddo
*do,i,1,15
e,10300+i,20015+i*100
*enddo
*do,i,1,15
e,10400+i,20100+i
*enddo
e,1,10101
e,1,10401
e,36,10115
e,36,10201
e,71,10215
e,71,10315
e,106,10301
e,106,10415
*do,j,1,15
*do,i,1,14
e,20000+i*100+j,20100+i*100+j
*enddo
*enddo
*do,i,1,15
*do,j,1,14
e,20000+i*100+j,20001+i*100+j
*enddo
*enddo
*do,i,1,14
e,11100+i,11101+i
e,11200+i,11201+i
e,11300+i,11301+i
e,11400+i,11401+i !上端板四条边连接成单元
*enddo
*do,i,1,15
e,11100+i,30001+i*100
*enddo
*do,i,1,15
e,11200+i,31500+i
*enddo
*do,i,1,15
e,11300+i,30015+i*100
*enddo
*do,i,1,15
e,11400+i,31500+i
*enddo
e,35,11101
e,35,11401
e,70,11115
e,70,11201
e,105,11215
e,105,11315
e,140,11301
e,140,11415
*do,j,1,15
*do,i,1,14
e,30000+i*100+j,30100+i*100+j
*enddo
*enddo
*do,i,1,15
*do,j,1,14
e,30000+i*100+j,30001+i*100+j
*enddo
*enddo
d,10208,Ux,0
d,10208,Uy,0
d,10208,Uz,0
d,10408,Ux,0
d,10408,Uy,0
d,10408,Uz,0
!*do,i,1,9
!d,2005+i*100,Ux,0
!d,2005+i*100,Uy,0
!d,2005+i*100,Uz,0
!*enddo
d,11208,Ux,0
d,11208,Uz,0
d,11408,Ux,0
d,11408,Uz,0
!*do,j,1,4
!*do,i,1,3
!d,8+(i-1)*5+(j-1)*20,uz,0
!*enddo
!*enddo


!d,825,Ux,0
!d,825,Uz,0
!d,845,Ux,0
!d,845,Uz,0
!d,825,Uy,-0.100
!d,845,Uy,-0.100
!DDELE,825,UY
!DDELE,845,UY
!d,825,Uy,-0.018
!d,845,Uy,-0.018
!静态分析
/solu
antype,static
OUTRES,ALL,ALL
PSTRES,ON
f,11208,fy,-1
f,11408,fy,-1
NSUBST,10,50,10
SOLVE
save
FINI
!施加初始缺陷
/SOLU
ANTYPE,BUCKLE !p415 特征值分析
BUCOPT,LANB,1
MXPAND,1
SOLVE
finish
/POST1
SET,FIRST
SET,LIST
PLDISP,1
FINI
/PREP7
UPGEOM,F0,1,1,2.4CHP1-0.248,RST
FINISH

/SOLU
TIME,1
ANTYPE,0
esel,s,type,,2
ekill,all
allsel
NLGEOM,ON !打开大变形开关
nropt,full
OUTRES,ALL,ALL
f,11208,fy,-135e3 !!!!!!!!
f,11408,fy,-135e3 !!!!!!!!
NSUBST,10,10,10
KBC,0
AUTOTS,1
SOLVE
/SOLU
time,2
ANTYPE,STATIC
nlgeom,on
nropt,full
fcum,add
outres,all,all
kbc,0
autots,1
SSTIF,ON
esel,s,TYPE,,2
EALIVE,all
allsel
nsubst,10,10,10
solve
/SOLU
time,3
ANTYPE,STATIC
nlgeom,on
nropt,full
upcoord,1,off
outres,all,all
kbc,0
autots,1
SSTIF,ON
d,11208,uy,-0.03
d,11408,uy,-0.03
nsubst,50,50,50
solve
finish
!提取结果
/post26
numvar,50
nsol,2,88,U,X
prod,7,2,,,,,,1000 !构件中点处荷载-挠度曲线
Rforce,3,10208,F,Y
Rforce,4,10408,F,Y
ADD,5,3,4, , , , ,1,1, ,
prod,6,5,,,,,,1/1000

/axlab,x,U(mm)
/axlab,y,P(kN)
xvar,7
plvar,6
PRVAR,6,7,, , , ,
ESOL,10,86,88,smisc,41, !近载2、4#管轴向应变
ESOL,11,86,88,smisc,42, !近载2、4#管弯矩应变
ADD,12,10,11, , , , ,-1e6,-1e6, ,
ESOL,15,78,80,smisc,41, !近载2、4#管轴向应变(L/4)
ESOL,16,78,80,smisc,42, !近载2、4#管弯矩应变(L/4)
ADD,17,15,16, , , , ,-1e6,1e6, ,
ESOL,18,94,96,smisc,41, !近载2、4#管轴向应变(L/4)
ESOL,19,94,96,smisc,42, !近载2、4#管弯矩应变(L/4)
ADD,14,18,19, , , , ,-1e6,1e6, ,
nsol,9,11208,U,y
nsol,20,35,U,y
nsol,21,70,U,y
nsol,22,105,U,y
nsol,23,140,U,y
nsol,26,1,U,y
nsol,27,36,U,y
nsol,28,71,U,y
nsol,29,106,U,y
prod,13,9,,,,,,-1000
nsol,30,88,U,Y !!!!!!!!中点纵向位移