哪位用ansys的link8算过问题啊？ link8模拟的横向连接杆
提取的轴力老是为零，杆变形后明显被拉长了
命令流如下：
/cle
*set,x5,1800000000
*set,y5,170000000
*set,x4,1800000000
*set,y4,180000000
*set,x3,1800000000
*set,y3,200000000
*set,x2,1800000000
*set,y2,180000000
*set,x1,1800000000
*set,y1,170000000
*set,k4,108000000
/PREP7
ET,1,BEAM4
et,2,link8 !定义link8单元
SECTYPE, 1, BEAM, I, ganggong, 0 !*定义1号工字形截面*
SECOFFSET, CENT !*截面质心不偏移*
SECDATA,0.1,0.1,0.2,0.0114,0.0114,0.007,0,0,0,0 !*1号截面参数*
R,1,0.0035578,0.0000237,0.00000158,0.1,0.2,0, !设置实常数，编号，截面面积，IZZ（水平轴）,IYY，Z轴方向厚度，Y轴方向厚度，转动
r,2,379.94e-6 !定义link8单元的面积
et,3,combin14
r,3,k4 !设置地层反力弹簧的实常数K，单位Pa
!定义弹簧单元
et,4,combin14
et,5,combin14
et,6,combin14
et,7,combin14
et,8,combin14
!设置X方向约束弹簧的实常数K，单位Pa
r,4,x1
r,5,x2
r,6,x3
r,7,x4
r,8,x5
et,9,combin14 !y方向弹簧
et,10,combin14
et,11,combin14
et,12,combin14
et,13,combin14
!设置y方向约束弹簧的实常数K，单位Pa
r,9,480e6
r,10,470e6
r,11,460e6
r,12,450e6
r,13,500e6
mp,EX,1,200e9 !材料特性（钢拱、杆件、弹簧）
mp,PRXY,1,0.3
mp,DENS,1,7850
!建立几何建模
k,1,100,100,0 !创建第一榀各关键点
k,2,93.544,97.621,0
k,3,100,106.88,0
k,4,106.456,97.621,0
larc,2,3,1,6.88
larc,3,4,1,6.88
k,101,100,100,1 !创建第二榀各关键点
k,102,93.544,97.621,1
k,103,100,106.88,1
k,104,106.456,97.621,1
larc,102,103,101,6.88
larc,103,104,101,6.88
k,201,100,100,2 !创建第三榀各关键点
k,202,93.544,97.621,2
k,203,100,106.88,2
k,204,106.456,97.621,2
larc,202,203,201,6.88
larc,203,204,201,6.88
k,301,100,100,-1
k,302,93.544,97.621,-1
k,303,100,106.88,-1
k,304,106.456,97.621,-1
larc,302,303,301,6.88
larc,303,304,301,6.88
k,401,100,100,-2
k,402,93.544,97.621,-2
k,403,100,106.88,-2
k,404,106.456,97.621,-2
larc,402,403,401,6.88
larc,403,404,401,6.88
!设置单元网格
lesize,1,,,22
lesize,2,,,22
lesize,3,,,22
lesize,4,,,22
lesize,5,,,22
lesize,6,,,22
lesize,7,,,22
lesize,8,,,22
lesize,9,,,22
lesize,10,,,22
!划分网格
type,1
mat,1
real,1
lmesh,all
type,2 !形成杆件单元
mat,1
real,2
e,93,48
*repeat,21,1,1
e,48,3
*repeat,21,1,1
e,3,138
*repeat,21,1,1
e,138,183
*repeat,21,1,1
e,135,90
*repeat,21,-1,-1
e,90,45
*repeat,21,-1,-1
e,45,180
*repeat,21,-1,-1
e,180,225
*repeat,21,-1,-1
e,92,47
e,47,2
e,2,137
e,137,182
!添加地层弹簧单元
!psprng,节点号，tran(径向施加弹簧),弹簧K值（需要根据地层抗力折算）,1（一个长度）
type,3
mat,1
real,3
local,1000,1,100,100,0
csys,1000
*do,i,1,45
psprng,i,tran,k4,1
*enddo
local,1001,1,100,100,1
csys,1001
*do,i,46,90
psprng,i,tran,k4,1
*enddo
local,1002,1,100,100,2
csys,1002
*do,i,91,135
psprng,i,tran,k4,1
*enddo
local,1003,1,100,100,-1
csys,1003
*do,i,136,180
psprng,i,tran,k4,1
*enddo
local,1004,1,100,100,-2
csys,1004
*do,i,181,225
psprng,i,tran,k4,1
*enddo
csys,0
!设施X向支座弹簧
type,4
mat,1
real,4
psprng,1,tran,x1,-1,0,0,,
psprng,24,tran,x1,1,0,0,,
type,5
mat,1
real,5
psprng,46,tran,x2,-1,0,0,,
psprng,69,tran,x2,1,0,0,,
type,6
mat,1
real,6
psprng,91,tran,x3,-1,0,0,,
psprng,114,tran,x3,1,0,0,,
type,7
mat,1
real,7
psprng,136,tran,x4,-1,0,0,,
psprng,159,tran,x4,1,0,0,,
type,8
mat,1
real,8
psprng,181,tran,x5,-1,0,0,,
psprng,204,tran,x5,1,0,0,,
!设施Y向支座弹簧
type,9
mat,1
real,9
psprng,1,tran,y1,0,-1,0,,
psprng,24,tran,y1,0,-1,0,,
type,10
mat,1
real,10
psprng,46,tran,y2,0,-1,0,,
psprng,69,tran,y2,0,-1,0,,
type,11
mat,1
real,11
psprng,91,tran,y3,0,-1,0,,
psprng,114,tran,y3,0,-1,0,,
type,12
mat,1
real,12
psprng,136,tran,y4,0,-1,0,,
psprng,159,tran,y4,0,-1,0,,
type,13
mat,1
real,13
psprng,181,tran,y5,0,-1,0,,
psprng,204,tran,y5,0,-1,0,,
!加Z向约束，加载，求等效节点荷载
/sol
nsel,s,node,,1,181,45 !加Z向节点约束
nsel,a,node,,24,204,45
d,all,uz,0
allsel
acel,0,10,0 !施加重力加速度
allsel !施加节点等效荷载（根据围岩压力用EXCEL计算出来的）
f, 1 ,fx, 20052.295893/3
*repeat,5,45,0
f, 3 ,fx, 40587.718035/3
*repeat,5,45,0
f, 4 ,fx, 41397.02104/3
*repeat,5,45,0
f, 5 ,fx, 41889.932123/3
*repeat,5,45,0
f, 6 ,fx, 42062.684072/3
*repeat,5,45,0
f, 7 ,fx, 41913.956588/3
*repeat,5,45,0
f, 8 ,fx, 41444.886315/3
*repeat,5,45,0
f, 9 ,fx, 40659.058335/3
*repeat,5,45,0
f, 10 ,fx, 39562.47862/3
*repeat,5,45,0
f, 11 ,fx, 38163.52813/3
*repeat,5,45,0
f, 12 ,fx, 36472.89891/3
*repeat,5,45,0
f, 13 ,fx, 34503.51226/3
*repeat,5,45,0
f, 14 ,fx, 32270.419825/3
*repeat,5,45,0
f, 15 ,fx, 29790.688865/3
*repeat,5,45,0
f, 16 ,fx, 27083.271705/3
*repeat,5,45,0
f, 17 ,fx, 24168.860695/3
*repeat,5,45,0
f, 18 ,fx, 21069.730325/3
*repeat,5,45,0
f, 19 ,fx, 17809.56681/3
*repeat,5,45,0
f, 20 ,fx, 14413.28721/3
*repeat,5,45,0
f, 21 ,fx, 10906.84875/3
*repeat,5,45,0
f, 22 ,fx, 7317.050685/3
*repeat,5,45,0
f, 23 ,fx, 3671.32941/3
*repeat,5,45,0
f, 2 ,fx, 0
*repeat,5,45,0
f, 25 ,fx, -3673.544105/3
*repeat,5,45,0
f, 26 ,fx, -7319.019365/3
*repeat,5,45,0
f, 27 ,fx, -10908.559655/3
*repeat,5,45,0
f, 28 ,fx, -14414.73236/3
*repeat,5,45,0
f, 29 ,fx, -17810.7419/3
*repeat,5,45,0
f, 30 ,fx, -21070.63448/3
*repeat,5,45,0
f, 31 ,fx, -24169.49675/3
*repeat,5,45,0
f, 32 ,fx, -27083.64596/3
*repeat,5,45,0
f, 33 ,fx, -29790.81098/3
*repeat,5,45,0
f, 34 ,fx, -32270.30261/3
*repeat,5,45,0
f, 35 ,fx, -34503.1715/3
*repeat,5,45,0
f, 36 ,fx, -36472.35319/3
*repeat,5,45,0
f, 37 ,fx, -38162.798415/3
*repeat,5,45,0
f, 38 ,fx, -39561.58801/3
*repeat,5,45,0
f, 39 ,fx, -40658.03203/3
*repeat,5,45,0
f, 40 ,fx, -41443.75095/3
*repeat,5,45,0
f, 41 ,fx, -41912.739848/3
*repeat,5,45,0
f, 42 ,fx, -42061.414707/3
*repeat,5,45,0
f, 43 ,fx, -41888.639244/3
*repeat,5,45,0
f, 44 ,fx, -41395.733806/3
*repeat,5,45,0
f, 45 ,fx, -40586.465574/3
*repeat,5,45,0
f, 24 ,fx, -20051.681983/3
*repeat,5,45,0
f, 1 ,fy, 21057.88552/2
*repeat,5,45,0
f, 3 ,fy, 36280.811236/2
*repeat,5,45,0
f, 4 ,fy, 24494.544842/2
*repeat,5,45,0
f, 5 ,fy, 12521.069935/2
*repeat,5,45,0
f, 6 ,fy, 451.8981695/2
*repeat,5,45,0
f, 7 ,fy, -11620.727391/2
*repeat,5,45,0
f, 8 ,fy, -23604.537256/2
*repeat,5,45,0
f, 9 ,fy, -35407.940847/2
*repeat,5,45,0
f, 10 ,fy, -46940.7263/2
*repeat,5,45,0
f, 11 ,fy, -58114.750026/2
*repeat,5,45,0
f, 12 ,fy, -68844.610505/2
*repeat,5,45,0
f, 13 ,fy, -79048.300847/2
*repeat,5,45,0
f, 14 ,fy, -88647.835583/2
*repeat,5,45,0
f, 15 ,fy, -97569.846761/2
*repeat,5,45,0
f, 16 ,fy, -105746.14459/2
*repeat,5,45,0
f, 17 ,fy, -113114.23882/2
*repeat,5,45,0
f, 18 ,fy, -119617.81615/2
*repeat,5,45,0
f, 19 ,fy, -125207.17034/2
*repeat,5,45,0
f, 20 ,fy, -129839.58298/2
*repeat,5,45,0
f, 21 ,fy, -133479.64916/2
*repeat,5,45,0
f, 22 ,fy, -136099.54813/2
*repeat,5,45,0
f, 23 ,fy, -137679.25697/2
*repeat,5,45,0
f, 2 ,fy, -138204.79489/2
*repeat,5,45,0
f, 25 ,fy, -137674.43673/2
*repeat,5,45,0
f, 26 ,fy, -136094.17481/2
*repeat,5,45,0
f, 27 ,fy, -133473.82598/2
*repeat,5,45,0
f, 28 ,fy, -129833.41633/2
*repeat,5,45,0
f, 29 ,fy, -125200.76738/3
*repeat,5,45,0
f, 30 ,fy, -119611.28323/2
*repeat,5,45,0
f, 31 ,fy, -113107.68133/2
*repeat,5,45,0
f, 32 ,fy, -105739.66491/2
*repeat,5,45,0
f, 33 ,fy, -97563.543025/2
*repeat,5,45,0
f, 34 ,fy, -88641.80105/2
*repeat,5,45,0
f, 35 ,fy, -79042.622595/2
*repeat,5,45,0
f, 36 ,fy, -68839.36846/2
*repeat,5,45,0
f, 37 ,fy, -58110.01612/2
*repeat,5,45,0
f, 38 ,fy, -46936.56345/2
*repeat,5,45,0
f, 39 ,fy, -35404.402435/2
*repeat,5,45,0
f, 40 ,fy, -23601.666775/2
*repeat,5,45,0
f, 41 ,fy, -11618.557755/2
*repeat,5,45,0
f, 42 ,fy, 453.34518/2
*repeat,5,45,0
f, 43 ,fy, 12521.783395/2
*repeat,5,45,0
f, 44 ,fy, 24494.52505/2
*repeat,5,45,0
f, 45 ,fy, 36280.069935/2
*repeat,5,45,0
f, 24 ,fy, 21057.336625/2
*repeat,5,45,0
/sol
allsel,all
SOLVE
/post1
ALLSEL
esel,s,type,,2 !选取杆单元
eplot
plnsol,u,sum
prnsol,dof
etable,s-axis,ls,1 !显示轴力
smult,n-axis,s-axis,,2e-4,1
pretan,n-axis
plls,n-axis,n-axis,1,0
