ANSYS在任意面施加任意载荷

  • 格式:doc
  • 大小:22.50 KB
  • 文档页数:2

【讨论】在任意面施加任意载荷(总结)

一:关于压力的加载(感谢和aggie网友的探讨)

加载过程在循环后进行,利用SFFUN命令把节点与节点压力相匹配。

缺点:每个节点需要匹配,导致数组维数太大。

如,下面指定节点

1、2,3,4处的表面载荷值。

*DIM,ABC,ARRAY,4

ABC(1)=10,20,30,40

SFFUN,PRES,ABC(1)

SF,ALL,PRES,100 !匹配后,还需加载

则节点1处的压力值为110,节点2处压力载荷为120,依次类推。

Main Menu>Solution>Settings>Node Function命令用于把数组变量为节点加载。

在圆柱外表面0~180度之间加余弦分布载荷。

全部程序如下:

/prep7

*afun,deg

et,1,45

mp,ex,1,2e11

mp,prxy,1,0.3

cyl4,,,0.5,,,,1

smrtsize,1

vsweep,all

*get,etotnum,node,,count

*dim,t1,array,etotnum,1,1

csys,1

*do,i,1,etotnum

*if,ny(i),gt,0,then

*if,ny(i),lt,180,then

*if,nx(i),eq,0.5,then

t1(i)=100*cos(ny(i))

*else

t1(i)=0

*endif

*endif

*endif

*enddo

sffun,pres,t1(1) !建立节点号与表面载荷的对应关系

sf,all,pres,0 !加载

二:关于力和位移的加载(这里转David2003网友的帖子,原贴由于设置积分限制,很多网友看不见,这里一并转来,做一个总结)

在内孔120度范围内的面上加载与位置有关的面力,还有一部分施加了位移约束。

/prep7

*afun,deg

et,1,45

mp,ex,1,2e11

mp,prxy,1,0.3

wpave,0,0,0

cyl4,0,0,1,360,1.5,360,0.5

smrtsize,1

vsweep,all

*get,etotnum,node,,count

*status

CSWPLA,11,1,1,1

nrotat,all

cscir,11,1

*do,i,1,etotnum

*if,ny(i),gt,0,then

*if,ny(i),lt,60,then

*if,nx(i),eq,1,then

d,i,all,0

*endif

*endif

*endif

*enddo

*do,i,1,etotnum

*if,ny(i),gt,120,then

*if,ny(i),lt,240,then

*if,nx(i),eq,1,then

p=cos(ny(i))

f,i,fx,p

*endif

*endif

*endif

*enddo

eplot

几点说明

1:加载任意函数载荷时,由于有限元计算采取离散加载到结点的方法,使得加载数值在相邻结点之间以直线近似代替曲线函数。这就要求网格划分的足够细,才能提高近似精度。

2:对于施加压力,结点上压力(df/ds)的单位是力/面积,对每个节点压力数值积分后的值就是总力的大小,对结点施加压力不能象加力那样(可以在单个结点上加),而应该选用数组匹配的方式加载。