24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tytccn

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3
帖子: 6
在线: 1.5小时
虫号: 3201166
注册: 2014-05-12
专业: 固体力学

[求助] 计算过程中出现了很大位移，但算出来的结果显示却是正常的已有1人参与

计算中出现

但结果中并没有异常的位移。
例如UZ，非常的小，而上图中显示在122节点有1080000的位移。不知道是不是我理解错了。

最后附上代码。一个很简单的模型，一根梁分成若干单元，与地面之间用弹簧支撑，约束住整根梁的ux，roty，rotz和梁两端的uz。荷载是一个集中力。

finish
/clear,start
/prep7

E2=2e4 !弹性模量
v2=0.167 !泊松比
dens2=2500 !密度

mp,ex,2,E2
mp,prxy,2,v2
mp,dens,2,dens2

et,201,188
keyopt,201,3,3
sectype,1,beam,rect,ganggui
secdata,90,90,5,5

et,301,14
keyopt,301,2,2
r,301,7.5e4

k,9,0,500,0
k,10,0,500,6000
l,9,10
k,11,0,0,0

lesize,all,25

lsel,s,,,1
latt,2,,201,,11,,1
lmesh,all

allsel,all
nsel,s,loc,y,500-0.001,500+0.001
cm,allnode,node

cmsel,s,allnode,node
ngen,2,1000,all,,,0,-100,0
nsel,u,,,1,1000
cm,diceng_node,node

cmsel,s,allnode,node
*get,nodenum,node,0,COUNT
*dim,d_nodelist,array,nodenum,2
*get,d_nodelist(1,1),node,0,num,min
d_nodelist(1,2)=d_nodelist(1,1)+1000
*do,in,2,nodenum,1
d_nodelist(in,1)=ndnext(d_nodelist(in-1,1))
d_nodelist(in,2)=ndnext(d_nodelist(in-1,1))+1000
*enddo

*do,in,1,nodenum,1
type,301
real,301
e,d_nodelist(in,1),d_nodelist(in,2)
*enddo

cmsel,s,diceng_node,node
d,all,all,0
cmsel,s,allnode,node
d,all,ux,0
d,all,roty,0
d,all,rotz,0
nsel,r,node,,1,2
d,all,uz,0

allsel,all
nsel,s,node,,117
f,all,fy,-10000

allsel,all

回复此楼