大位移引起的收斂問題

1、討論】問題懸賞解答課題之三（大位移引起的收斂問題）出了點(diǎn)問題，重新編輯了以下，如下：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é)點(diǎn)1位移(ux,uy,rotz ),固定節(jié)點(diǎn)2轉(zhuǎn)動自由度(rotz), 在節(jié)點(diǎn) 2 施加

3、水平載荷；卸掉水平載荷；固定節(jié)點(diǎn) 202(頂點(diǎn))轉(zhuǎn)動自由 度，解除對節(jié)點(diǎn) 2 轉(zhuǎn)動自由度的限制； 對節(jié)點(diǎn) 202 施加一定軌跡的位移載 荷。1 看看程序中的加載存在什么問題？2 程序中有什么不合理的地方嗎？ 這是一個(gè)實(shí)際的工程問題，就是求在這個(gè)加載過程后金屬線的形狀。雖然 結(jié)果收斂，但是還有點(diǎn)問題。prvar(ux,uy) 的結(jié)果(只顯示了有問題的一部分) ：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可以看到，第三個(gè)載荷步完成后，施加的力載荷引起頂部（節(jié)點(diǎn)202 ）的位移分別是 UX=-0.266964 和 UY=-0.110036 。然后從第四個(gè)載荷步開始 施加位移載荷，但是第四步施加的位移載荷并沒有累計(jì)前三步力載荷所引起的位移。這與實(shí)際情況不符，實(shí)際情況下，t=4時(shí)，UX=0.05+ （-0.266 ）,UY=-0.389+ （ -0.11 ）。從位移圖中也可以看出，加載位移載荷之

5、前， UX， UY 都先歸零，然后再 加載位移，而實(shí)際情況應(yīng)該是：加載的位移載荷在以前的力載荷引起的位 移的基礎(chǔ)上進(jìn)行?？纯磫栴}出在哪里，如何解決？在第四步之前加一句：DCUM,ADD試試，不過后面所加位移量要做相應(yīng)的修改，應(yīng)該加上前后兩步的位移增量。我現(xiàn)在沒時(shí)間驗(yàn)算。另外可以看看命令手冊中DCUM的說明，也許會有些啟發(fā)。恩，我把程序該了一下，設(shè)置了 DCUM，ADD，在課題二中僅僅加位移載荷的情況下運(yùn)行成功。但在這里還是不行，加載位移載荷之前，UX，UY還是歸零。不知還有什么辦法，會不會和預(yù)應(yīng)力之類的設(shè)置有關(guān)？ ？改動部分為：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)在的問題是力邊界產(chǎn)生的位移那部分無法累計(jì)。 有個(gè)想法：先加力邊界，算出作用點(diǎn)的位移；將該位移作為邊

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

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位移歸零，錯(cuò)誤，如圖所示:改進(jìn)一下，把節(jié)點(diǎn)2位移載荷引起節(jié)點(diǎn)202的位移（0.025 ）給彌補(bǔ)上去。time,1d,2,ux,0.01solve time,2d,202,ux,0.01+ 0.025solve雖然最終點(diǎn)加載大小相同，但加載過程不對，還是歸零了，錯(cuò)誤！如圖所 示:time,1d,2,ux,0.01d,202,ux,0.025!（先計(jì)算然后再加到這里，可能用 *GET或其他命令啟動重新分析會更好，先不管這些）solve time,2d,202,ux,0.01天啊，這不是我

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