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)的单位是力/面积,对每个节点压力数值积分后的值就是总力的大小,对结点施加压力不能象加力那样(可以在单个结点上加),而应该选用数组匹配的方式加载。