fini
/clear
/filname,xiaohengliang
/title,yingxiangxian
/prep7
et,1,beam3
mp,ex,1,206e9
mp,prxy,1,0.3
r,1,4.851e-3,5.017e-5,0.25
k,1,1.5
k,2,1.5+1.625
k,3,1.5+1.625*2
k,4,1.5+1.625*2+3
k,5,1.5+1.625*2+3+2
k,6,1.5+1.625*2+3+2+2.75
k,7,1.5+1.625*2+3+2+2.75+1.2
k,8,1.5+1.625*2+3+2+2.75+1.2+0.3
*do,i,1,7
l,i,i+1
*enddo
lsymm,x,all
l,9,1
lesize,all,0.05
lmesh,all
/solu
dk,all,uy
ksel,s,loc,x,1.5
dk,all,ux
nsel,s,loc,x,-(1.5+1.625*2+3+2+2.75+1.2),-(1.5+1.625*2+3+2)
n0=node(-13.7,0,0)
nmax=ndinqr(0,13)
*dim,p1node,,nmax
*dim,nodex,,nmax
p1node(1)=n0
nodex(1)=nx(n0)
!n00=node(-13.75,0,0)
!nsel,u,,,n00
*do,i,2,nmax
ni=nnear(n0)
p1node(i)=ni
nodex(i)=nx(ni)
nsel,u,,,n0
n0=ni
*enddo
allsel,all
*do,i,1,nmax
time,i
fdele,all,all
f,p1node(i),fy,-1
solve
*enddo
/post26
nsol,3,72,u,y,uy_72
vput,nodex,2
/xrange,nx(p1node(1)),nx(p1node(nmax)) !定义影响线长度范围
/axlab,x,length(m) !定义X轴名称
!/gropt,divx,(nx(p1node(nmax))-nx(p1node(1)))/10!定义X轴刻度数
/gropt,view,on !将所绘制影响线图设为可缩放模式(缺省不缩放)
nodefan=node(-12.5,0,0)
rforce,4,nodefan,f,y,ry_2
/axlab,y,ry_2
plvar,4
|