【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
模型和边界条件程序已经给出：/prep7
et,1,beam23,,1,,1,,2
r,1,0.05
mp,ex,1,2790
mp,nuxy,,0.42
tb,biso,1
tbdata,1,23,0
K,1,0,0,0,
K,2,0,3.6,0,
l,1,2
lesize,all,,,100
lmesh,all
/solu
nlgeom,on
!pred,on,,on
d,1,ux
d,1,uy
d,1,rotz
d,2,rotz
------
------
要求顶点按照图示的椭圆（椭圆的两个半轴分别为3.6和3.4）弧线移动，请写出位移加载的命令流（用循环方式）。
回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
草草做了一下，估计还没有吃透斑竹的愿意，算是抛砖引玉吧......FINI
/CLE
/prep7
et,1,beam23,,1,,1,,2
r,1,0.05
mp,ex,1,2790
mp,nuxy,,0.42
tb,biso,1
tbdata,1,23,0
K,1,0,0,0,
K,2,0,3.6,0,
l,1,2
lesize,all,,,100
lmesh,all
D,1,UX,0
D,1,UY,0
D,1,ROTZ,0
D,2,ROTZ,0
/SOLU
ANTYPE,0
DELTAX=0.2
TNUM=3.4/DELTAX
*DO,I,1,TNUM
TIME,I
XX=DELTAX*I
YY=1-XX*XX/(3.4*3.4)
*IF,YY,LT,1E-6,THEN
YY=0
*ENDIF
D,2,UX,XX
D,2,UY,-3.6+3.6*SQRT(YY)
NLGEOM,ON
OUTRES,ALL,ALL
SOLV
*ENDDO
回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
pldisp:回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
太棒了，说不定老兄一帖就把这个课题画上句号了，呵呵！！！回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
老兄过奖了，由于时间关系，草草做了一下。不知道老兄出这个课题的愿意是什么，主要想要解决哪方面的问题。回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
恩，这个问题比较单一，主要就是练练结构非线形分析中大位移载荷的施加。还有就是练练循环了。简单的加载就是直接写上载荷的值，比较复杂的载荷可以通过数组参数或者函数边界条件的方法来施加。函数边界条件方法的变量一般是时间、温度或者位置。该题中，施加的两个位移载荷之间具有一定的函数关系。
马上就有课题三了，比较复杂一些，到时候要参与啊！！！
回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
不要考虑曲屈么?回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
我想，如果施加的是力载荷，屈曲发生时可能会导致计算不收敛，施加的是位移载荷时，应该不会导致结构的不稳定，所以在整个过程中都可以得到收敛的解。回复： 【讨论】问题悬赏解答——课题之二（函数位移载荷的施加）
上面的程序是位移载荷替代形式的加载，以下是累加形式的加载：FINI
/CLEAR
/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,3.6,0,
l,1,2
lesize,all,,,100
lmesh,all
/solu
D,1,UX,0
D,1,UY,0
D,1,ROTZ,0
D,2,ROTZ,0
ANTYPE,0
DELTAX=0.1
TNUM=3.4/DELTAX
*DIM,YY,,TNUM
DCUM,ADD
*DO,I,1,TNUM
TIME,I
XX=DELTAX
YY(I)=1-XX*XX*I*I/(3.4*3.4)
*IF,YY(I),LT,1E-6,THEN
YY(I)=0
*ENDIF
YY(I)=3.6*SQRT(YY(I))
D,2,UX,XX
*IF,I,EQ,1,THEN
D,2,UY,YY(I)-3.6
*ELSE
D,2,UY,YY(I)-YY(I-1)
*ENDIF
NLGEOM,ON
OUTRES,ALL,ALL
SOLV
*ENDDO
/post1
pldisp 高啊，实在是高。小弟PH！！！！！！！ 真高呀，佩服佩服！！！！！！！！！！！！！！！！！1
页:
[1]