大位移引起的收敛问题.docVIP

【讨论】问题悬赏解答——课题之三（大位移引起的收敛问题）出了点问题，重新编辑了以下，如下： fini /clear /config,nres,5000 /prep7 et,1,beam23,,1,,1,,2 r,1,0.025 mp,ex,1,2790 mp,nuxy,,0.42 tb,biso,1 tbdata,1,23,0 K,1,0,0,0? ?K,2,0,0.4,0 k,3,0,3.6,0 l,1,2 l,2,3 lesize,1,,,200 lesize,2,,,400 lmesh,all /solu d,1,ux d,1,uy d,1,rotz d,2,rotz nlgeom,on nsubst,120 outres,all,all time,1 f,2,fx,-0.0006 solve time,2 f,2,fx,0 solve time,3 d,202,rotz ddele,2,rotz solve DELTAX=0.05??TNUM=3.4/DELTAX??*DO,I,1,TNUM????TIME,I+3????XX=DELTAX*I????YY=1-XX*XX/(3.4*3.4)????*IF,YY,LT,1E-6,THEN??????YY=0????*ENDIF????D,202,UX,XX????D,202,UY,-3.6+3.6*SQRT(YY)????SOLV??*ENDDO??/post1 pldisp 加载过程是：固定节点1位移（ux,uy,rotz）,固定节点2转动自由度（rotz）,在节点2施加水平载荷；卸掉水平载荷；固定节点202（顶点）转动自由度，解除对节点2转动自由度的限制；对节点202施加一定轨迹的位移载荷。 ??1看看程序中的加载存在什么问题？？？ 2程序中有什么不合理的地方吗？ 这是一个实际的工程问题，就是求在这个加载过程后金属线的形状。虽然结果收敛，但是还有点问题。 prvar(ux,uy)的结果（只显示了有问题的一部分）： ??????TIME? ?? ???202 UX? ?? ???202 UY? ?? ? ??????????????????UX? ?? ?? ?? ?UY? ?? ? ???2.9315? ?? ? -0.266964? ???-0.110036? ??????3.0000? ?? ? -0.266964? ???-0.110036? ??????3.0005? ?? ???0.260417E-04 -0.202758E-06 ???3.0010? ?? ???0.520833E-04 -0.405515E-06 ???3.0018? ?? ???0.911458E-04 -0.709651E-06 ... ???3.9522? ?? ???0.476104E-01 -0.370689E-03 ???4.0000? ?? ???0.500000E-01 -0.389294E-03 可以看到，第三个载荷步完成后，施加的力载荷引起顶部（节点202）的位移分别是UX=-0.266964和UY=-0.110036。然后从第四个载荷步开始施加位移载荷，但是第四步施加的位移载荷并没有累计前三步力载荷所引起的位移。这与实际情况不符，实际情况下，t=4时，UX=0.05+（-0.266），UY=-0.389+（-0.11）。 从位移图中也可以看出，加载位移载荷之前，UX，UY都先归零，然后再加载位移，而实际情况应该是：加载的位移载荷在以前的力载荷引起的位移的基础上进行。 看看问题出在哪里，如何解决？？？ 在第四步之前加一句：DCUM,ADD试试，不过后面所加位移量要做相应的修改，应该加上前后两步的位移增量。我现在没时间验算。另外可以看看命令手册中DCUM的说明，也许会有些启发。 恩，我把程序该了一下，设置了DCUM，ADD，在课题二中仅仅加位移载荷的情况下运行成功。但在这里还是不行，加载位移载荷之前，UX，UY还是归零。不知还有什么办法，会不会和预应力之类的设置有关？？？ ??改动部分为： DELTAX=0.1??TNUM=3.4/DELTAX *DIM,YY,,TNUM DCUM,ADD *DO,I,1,TNUM????TIME,I+3????XX=DELTAX????YY(I)=1-XX*XX*I*I/(3.466*3.466)????*IF,YY(I),LT,1E-6,THEN??????YY(I)=0????*ENDIF ??YY(I)=3.49*SQRT(YY(I)) ??D,202,