资源描述：
《大位移引起的收敛问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、【讨论】问题悬赏解答——课题之三（大位移引起的收敛问题）出了点问题，重新编辑了以下，如下：fini/clear/config,nres,5000/prep7et,1,beam23,,1,,1,,2r,1,0.025mp,ex,1,2790mp,nuxy,,0.42tb,biso,1tbdata,1,23,0K,1,0,0,0 K,2,0,0.4,0k,3,0,3.6,0l,1,2l,2,3lesize,1,,,200lesize,2,,,400lmesh,all/solud,1,uxd,1,uyd,1,rotzd,2,rotznlg
2、eom,onnsubst,120outres,all,alltime,1f,2,fx,-0.0006solvetime,2f,2,fx,0solvetime,3d,202,rotzddele,2,rotzsolveDELTAX=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、,-3.6+3.6*SQRT(YY) SOLV *ENDDO /post1pldisp加载过程是：固定节点1位移（ux,uy,rotz）,固定节点2转动自由度（rotz）,在节点2施加水平载荷；卸掉水平载荷；固定节点202（顶点）转动自由度，解除对节点2转动自由度的限制；对节点202施加一定轨迹的位移载荷。 1看看程序中的加载存在什么问题？？？2程序中有什么不合理的地方吗？这是一个实际的工程问题，就是求在这个加载过程后金属线的形状。虽然结果收敛，但是还有点问题。prvar(ux,uy)的结果（只显示了有问题的一部分）：
4、 TIME 202UX 202UY 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...
5、 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，U
6、Y都先归零，然后再加载位移，而实际情况应该是：加载的位移载荷在以前的力载荷引起的位移的基础上进行。看看问题出在哪里，如何解决？？？在第四步之前加一句：DCUM,ADD试试，不过后面所加位移量要做相应的修改，应该加上前后两步的位移增量。我现在没时间验算。另外可以看看命令手册中DCUM的说明，也许会有些启发。恩，我把程序该了一下，设置了DCUM，ADD，在课题二中仅仅加位移载荷的情况下运行成功。但在这里还是不行，加载位移载荷之前，UX，UY还是归零。不知还有什么办法，会不会和预应力之类的设置有关？？？ 改动部分为：DELTAX=0.1
7、 TNUM=3.4/DELTAX*DIM,YY,,TNUMDCUM,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,UX,XX *IF,I,EQ,1,THEN D,202,UY,YY(I)-3.49 *ELSE D,202,UY,YY(I)-YY(I-1) *E
8、NDIF NLGEOM,ON OUTRES,ALL,ALL SOLV *ENDDO也就是说现在的问题是力边界产生的位移那部分无法累计。有个想法：先加力边界，算出作用点的位移；将该位移作为边界替代原来的力边
