老悬赏解答之函数位移荷载，一例，敬请版主们关注

flyingcow

为了使你更好的理解这个问题，请查看：http://www.simwe.com/forum/thread-741522-1-1.html 。这是我在ABAQUS版发了的帖子。今天在这里再发有两方面的原因：

1. 本来就期待问题得到解决，不管是用ANSYS 还是 ABAQUS
2. 今天很高兴看到版主们很久以前的“问题悬赏解答”课题之二：函数位移荷载的施加。这个帖子给我莫大的启发。使我想到，除了用原来所设计的方法（请阅读），我又看到了问题解决的曙光，即函数位移荷载。

于是我马不停蹄，在认真学习版主的帖子及随后的跟帖后，定义了一个一空间曲线（如后图）作为所要施加的函数位移，也即节点N2/N3的变形轨迹。

但是，不知为何，结构发生的变形十分诡异，好像发生了屈曲。请各位高手帮看看。并且，也不能实现100％荷载施加，即*DO,ii,1,NN
*ENDDO

这是按照版主们的帖子修改的APDL

finish
/clear
/FILNAME,triFrame,0

LL=2
pi=3.1415926535897932
H =LL*cos(pi/6)
R =LL*sin(pi/6)
h0=0.25/12
b0=1/12/12

/prep7
k,1
k,2, LL*sin(pi/6),0,LL*cos(pi/6)
k,3, 0,           0,LL*cos(pi/6)

l,1,2
l,2,3

et,1,beam4

mp,ex,1,200e9

mp,prxy,1,0.3

r,1,b0*h0,b0*h0*h0*h0/12,h0*b0*b0*b0/12,h0,b0

lesize,1,,,40
lesize,2,,,20
lmesh,all

d,1,ux,0
d,1,uy,0
d,1,uz,0
d,1,rotx,0
d,1,roty,0
d,1,rotz,0

d,42,ux,0
d,42,uy,0
! d,42,uy,0
d,42,rotx,0
d,42,roty,0
! d,42,rotz,0

/solu
antype,0
! antype, static

NN=200  ! 可以用来调节每一步荷载增加幅度的大小

*DO,ii,1,10
  time,ii

  i0=1  ! NODE＃2 的初始位置
  theta0= 2*pi - 2*pi/(NN-2)*(i0-1) !＝2*pi
  z0    = H - H/(NN-2)*(i0-1)
  x0    = z0/H*R*cos(theta0)
  y0    = z0/H*R*sin(theta0)

  i1=ii+1  ! NODE＃2 的新位置
  theta1= 2*pi - 2*pi/(NN-2)*(i1-1)
  z1    = H - H/(NN-2)*(i1-1)
  x1    = z1/H*R*cos(theta1)
  y1    = z1/H*R*sin(theta1)

  delz = z1-z0  ! 所需变形
  delx = x1-x0
  dely = y1-y0

  d,2,ux,delx
  d,2,uy,dely
  d,2,uz,delz

  nlgeom,on

  outres,all,all

  solve

*enddo

/post1
pldisp

以下是我所获得的变形

以及定义的函数位移－－实质是一绕圆锥表面的曲线

flyingcow

http://www.simwe.com/forum/viewt ... C%C9%CD%BD%E2%B4%F0