请教一下,为什么下面三段程序计算的挠度不同?谢谢!
程序1为梁保持不变,仅改变重力的方向
程序2 为梁旋转60.36度,重力沿y方向,改变了节点坐标方向
程序3,没有改变节点坐标方向,其他同程序2.
程序1
************创建截面***************
fini
/clear
/prep7
et,1,82
blc5,0,0,1,1.2
blc5,0,0,0.8,1
absa,1,2
aesize,all,0.1
amesh,all
secwrite,sect1,sect,,1
fini
/clear
/prep7
et,1,82
pcirc,0.3,0,0,360
aesize,all,0.1
amesh,all
secwrite,sect2,sect,,1
fini
************创建单元求解***************
/clear
/prep7
!clocal,11,0,,,,60.36
!csys,11
et,1,beam188
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,ex,2,2.1e20
sectype,1,beam,mesh
secoffset,origin
secread,'sect1',,,mesh
sectype,2,beam,mesh
secoffset,cent
secread,'sect2',,,mesh
k,1,
k,2,70
k,3,70,0.1
lstr,1,2
lsel,s,,,1
latt,1,,1,,3,,1
lesize,1,,,70
lmesh,1
n,800,70,0.449
n,801,70,-0.559
type,1
mat,2
secnum,2
n,1001,,,0.1
e,2,800,1001
e,2,801,1001
!csys,0
d,1,all
d,1001,all
acel,8.52,4.85 !重力加速度做了分解
fini
/solu
solve
程序2
************创建截面***************
fini
/clear
/prep7
et,1,82
blc5,0,0,1,1.2
blc5,0,0,0.8,1
absa,1,2
aesize,all,0.1
amesh,all
secwrite,sect1,sect,,1
fini
/clear
/prep7
et,1,82
pcirc,0.3,0,0,360
aesize,all,0.1
amesh,all
secwrite,sect2,sect,,1
fini
************创建单元求解***************
/clear
/prep7
clocal,11,0,,,,60.36
csys,11
et,1,beam188
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,ex,2,2.1e20
sectype,1,beam,mesh
secoffset,origin
secread,'sect1',,,mesh
sectype,2,beam,mesh
secoffset,cent
secread,'sect2',,,mesh
k,1,
k,2,70
k,3,70,0.1
lstr,1,2
lsel,s,,,1
latt,1,,1,,3,,1
lesize,1,,,70
lmesh,1
n,800,70,0.449
n,801,70,-0.559
type,1
mat,2
secnum,2
n,1001,,,0.1
e,2,800,1001
e,2,801,1001
csys,0
nsel,s,,,all
nrotat,all
d,1,all
d,1001,all
acel,,9.8
fini
/solu
solve
程序3
************创建截面***************
fini
/clear
/prep7
et,1,82
blc5,0,0,1,1.2
blc5,0,0,0.8,1
absa,1,2
aesize,all,0.1
amesh,all
secwrite,sect1,sect,,1
fini
/clear
/prep7
et,1,82
pcirc,0.3,0,0,360
aesize,all,0.1
amesh,all
secwrite,sect2,sect,,1
fini
************创建单元求解***************
/clear
/prep7
clocal,11,0,,,,60.36
csys,11
et,1,beam188
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,ex,2,2.1e20
sectype,1,beam,mesh
secoffset,origin
secread,'sect1',,,mesh
sectype,2,beam,mesh
secoffset,cent
secread,'sect2',,,mesh
k,1,
k,2,70
k,3,70,0.1
lstr,1,2
lsel,s,,,1
latt,1,,1,,3,,1
lesize,1,,,70
lmesh,1
n,800,70,0.449
n,801,70,-0.559
type,1
mat,2
secnum,2
n,1001,,,0.1
e,2,800,1001
e,2,801,1001
d,1,all
d,1001,all
acel,,9.8
fini
/solu
solve
|