大位移引起的收斂問題

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,0K,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,rotz nlgeom,on nsubst,120 outres,all,all time,1 f

2、,2,fx,-0.0006 solve time,2f,2,fx,0 solve time,3 d,202,rotz ddele,2,rotz solve DELTAX=0.05 TNUM=3.4/DELTAX *DO,I,1,TNUMTIME,I+3XX=DELTAX*IYY=1-XX*XX/(3.4*3.4)*IF ,YY,LT,1E-6,THENYY=0*ENDIFD,202,UX,XXD,202,UY ,-3.6+3.6*SQRT(YY)SOLV*ENDDO/post1pldisp加載過程是：固定節(jié)點1位移(ux,uy,rotz ),固定節(jié)點2轉動自由度(rotz), 在節(jié)點 2 施加

3、水平載荷；卸掉水平載荷；固定節(jié)點 202(頂點)轉動自由 度，解除對節(jié)點 2 轉動自由度的限制； 對節(jié)點 202 施加一定軌跡的位移載 荷。1 看看程序中的加載存在什么問題？2 程序中有什么不合理的地方嗎？ 這是一個實際的工程問題，就是求在這個加載過程后金屬線的形狀。雖然 結果收斂，但是還有點問題。prvar(ux,uy) 的結果(只顯示了有問題的一部分) ：TIME202 UX202 UY2.9315-0.266964-0.1100363.0000-0.266964-0.1100363.00050.260417E-04-0.202758E-063.00100.520833E-04-0.405

4、515E-063.00180.911458E-04-0.709651E-063.95220.476104E-01-0.370689E-034.00000.500000E-01-0.389294E-03UXUY可以看到，第三個載荷步完成后，施加的力載荷引起頂部（節(jié)點202 ）的位移分別是 UX=-0.266964 和 UY=-0.110036 。然后從第四個載荷步開始 施加位移載荷，但是第四步施加的位移載荷并沒有累計前三步力載荷所引起的位移。這與實際情況不符，實際情況下，t=4時，UX=0.05+ （-0.266 ）,UY=-0.389+ （ -0.11 ）。從位移圖中也可以看出，加載位移載荷之

5、前， UX， UY 都先歸零，然后再 加載位移，而實際情況應該是：加載的位移載荷在以前的力載荷引起的位 移的基礎上進行?？纯磫栴}出在哪里，如何解決？在第四步之前加一句：DCUM,ADD試試，不過后面所加位移量要做相應的修改，應該加上前后兩步的位移增量。我現(xiàn)在沒時間驗算。另外可以看看命令手冊中DCUM的說明，也許會有些啟發(fā)。恩，我把程序該了一下，設置了 DCUM，ADD，在課題二中僅僅加位移載荷的情況下運行成功。但在這里還是不行，加載位移載荷之前，UX，UY還是歸零。不知還有什么辦法，會不會和預應力之類的設置有關？ ？改動部分為：DELTAX=0.1TNUM=3.4/DELTAX*DIM, YY

6、 ”TNUMDCUM,ADD *DO,I,1,TNUMTIME,I+3XX=DELTAXYY(I)=1-XX*XX*I*I/(3.466*3.466)*IF ,YY(I),LT,1E-6,THENYY(I)=0*ENDIFYY(I)=3.49*SQRT(YY(I)D,202,UX,XX*IF ,I,EQ,1,THEND,202,UY ,YY(I)-3.49*ELSED,202,UY ,YY(I)-YY(I-1)*ENDIFNLGEOM,ONOUTRES,ALL,ALLSOLV*ENDDO也就是說現(xiàn)在的問題是力邊界產生的位移那部分無法累計。 有個想法：先加力邊界，算出作用點的位移；將該位移作為邊

7、界替代原來 的力邊界(這需要新建一個與原來完全一樣的模型) ，這樣后續(xù)位移加載 就能累計了。不過這種方法有個問題就是力的加載點和后續(xù)位移加載點不 同的話可能無法累計(因為邊界不同) 。你可以調試一下。另外需要比較 一下用位移代替原來力加載后結構的應力分布，各點的位移情況是否一致。通過 yymavrerick 的啟發(fā)，我終于搞定了這個問題（我真想獎勵我 10 分， 哈哈），有兩種方案：一、重新建立模型：在力載荷作用完后，記錄每個節(jié)點的坐標值，然后運 用這些節(jié)點重新建立模型，在新模型基礎上施加位移載荷。二、等效冗余加載。 （下面通過一個例子來說明） 這個思路很好，運用等效冗余加載的方法可以解決所提

8、到的問題。 看以下命令流：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,0K,2,0,0.3,0K,3,0,0.6,0l,1,2l,2,3 lesize,ALL,200 lmesh,all/solu dcum,add d,1,ux d,1,uy d,1,rotz nlgeom,on nsubst,120 outres,all,all time,1 d,2,ux,0.01 solve time,2 d,202,ux,0

9、.01 solve /post26 nsol,2,202,u,x nsol,3,202,u,yplvar,2,3位移歸零，錯誤，如圖所示:改進一下，把節(jié)點2位移載荷引起節(jié)點202的位移（0.025 ）給彌補上去。time,1d,2,ux,0.01solve time,2d,202,ux,0.01+ 0.025solve雖然最終點加載大小相同，但加載過程不對，還是歸零了，錯誤！如圖所 示:time,1d,2,ux,0.01d,202,ux,0.025!（先計算然后再加到這里，可能用 *GET或其他命令啟動重新分析會更好，先不管這些）solve time,2d,202,ux,0.01天啊，這不是我

10、夢寐以求的過程曲線嗎？我激動得要哭了，哈哈!重復以下等效冗余加載得要點：節(jié)點M加載位移U1引起節(jié)點N位移U2D, M, UX, U1和D, M，UX，U1D，N，UX，U2的加載效果是一樣的，但是D，M，UX，U1D，N，UX，U2的加載方式可以通過設置 DCUM, ADD實現(xiàn)載荷步施加在不同節(jié)點時節(jié)點位移的正確累加。恩，別高興太早，可能還有別的問題，大家看看吧！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

11、,0,0,0 K,2,0,0.4,0 k,3,0,3.6,0 l,1,2l,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 /post1 set,last *get,x,node,202,u,x *get,y,node,202,u,y /solu

12、 lsclear,all 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 kbc,1 d,202,ux,x d,202,uy,y solve DELTAX=0.05 TNUM=3.4/DELTAX kbc,0 *DO,I,1,TNUMTIME,I+3XX二DELTAX*IYY=3.6*sqrt(abs(1-XX*XX/(3.4*3.4)D,202,UX,x+xxD,202,U Y, y+yy-3.6SOLV