在ANSYS中施加任意面方向变化载荷的方法

本文格式为Word版，下载可任意编辑——在ANSYS中施加任意面方向变化载荷的方法在任意面施加任意方向任意变化的压力在某些特别的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。假使利用ANSYS的参数化设计语言，也可以十分完美地实现此功能，下面通过一个小例子描述此方法。

!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!!!finish/prep7

et,500,shell63press=100e6amesh,allesla,snsla,s,1

!假使载荷的反向是一个特别坐标系的方向,可在此建立局部坐标系,并将!所有节点坐标系旋转到局部坐标系下.*get,enmax,elem,,num,maxdofsel,s,fx,fy,fz

fcum,add!!!将力的施加方式设置为\累加\而不是缺省的\替代\*do,i,1,enmax*if,esel,eq,1,then

*get,ae,elem,i,area!此命令用单元真实面积，如用投影面积，请用下几条命令

!*get,ae,elem,i,aproj,x!此命令用单元X投影面积，如用真实面积，请用上一条命令

!*get,ae,elem,i,aproj,y!此命令用单元Y投影面积!*get,ae,elem,i,aproj,z!此命令用单元Z投影面积xe=centrx!单元中心X坐标(用于求解压力值)ye=centry!单元中心Y坐标(用于求解压力值)ze=centrz!单元中心Z坐标(用于求解压力值)

!下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.p_e=(xe-10)*press+(ye-5)*pressf_tot=p_e*aeesel,s,elem,,insle,s,corner

*get,nn,node,,countf_n=f_tot/nn*do,j,1,nn

f,nelem(i,j),fx,f_n!压力的作用方向为X方向!f,nelem(i,j),fy,f_n!压力的作用方向为Y方向!f,nelem(i,j),fz,f_n!压力的作用方向为Z方向*enddo*endif

esla,s*enddoaclear,all

fcum,repl!!!将力的施加方式还原为缺省的\替代\dofsel,allallsel说明：本信息

在任意面施加任意方向任意变化的压力在某些特别的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。假使利用ANSYS的参数化设计语言，也可以十分完美地实现此功能，下面通过一个小例子描述此方法。

!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!!!finish/prep7

et,500,shell63press=100e6amesh,allesla,snsla,s,1

!假使载荷的反向是一个特别坐标系的方向,可在此建立局部坐标系,并将!所有节点坐标系旋转到局部坐标系下.*get,enmax,elem,,num,maxdofsel,s,fx,fy,fz

fcum,add!!!将力的施加方式设置为\累加\而不是缺省的\替代\*do,i,1,enmax*if,esel,eq,1,then

*get,ae,elem,i,area!此命令用单元真实面积，如用投影面积，请用下几条命令

!*get,ae,elem,i,aproj,x!此命令用单元X投影面积，如用真实面积，请用上一条命令

!*get,ae,elem,i,aproj,y!此命令用单元Y投影面积!*get,ae,elem,i,aproj,z!此命令用单元Z投影面积xe=centrx!单元中心X坐标(用于求解压力值)ye=centry!单元中心Y坐标(用于求解压力值)ze=centrz!单元中心Z坐标(用于求解压力值)

!下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.p_e=(xe-10)*press+(ye-5)*pressf_tot=p_e*aeesel,s,elem,,insle,s,corner

*get,nn,node,,count

f_n=f_tot/nn*do,j,1,nn

f,nelem(i,j),fx,f_n!压力的作用方向为X方向!f,nelem(i,j),fy,f_n!压力的作用方向为Y方向!f,nelem(i,j),fz,f_n!压力的作用方向为Z方向*enddo*endifesla,s*enddoaclear,all

fcum,repl!!!将力的施加方式还原为缺省的\替代\dofsel,all

allsel本信息来源：CAD教育网.

ANSYS5.7版本具有函数加载功能，可以很便利地在模型表面施加函数变化的各种载荷，在ANSYS5.6版本中，也可以通过变通的方式来实现此功能，其思路是：首先选定所要施加函数变化表面载荷的表面上的节点，利用ANSYS的参数数组和嵌入函数知识写一简单的命令流，定义好相应节点位置的面载荷值，然后通过在节点上施加面载荷来完成。

下面以在一圆柱表面施加函数变化载荷为例：/prep7et,1,45

cyl4,,,0.5,,,,3vsweep,all

asel,s,loc,y,0.01,1

nsla!

*get,nmax,node,,num,max,*get,nmin,node,,num,min,*afun,deg

*dim,t1,array,nmax,1,1,csys,1

*do,k,nmin,nmax*if,nsel(k),eq,1,thent1(k)=1000*sin(ny(k))*elset1(k)=0*endif*enddo!

sffun,pres,t1(1)sf,all,pres,0

利用ANSYS随机振动分析功能实现随机疲乏分析

ANSYS随机振动分析功能可以获得结构随机振动响应过程的各种统计参数（如：均值、均方根和平均频率等），根据各种随机疲乏寿命预计理论就可以成功地预计结构的随机疲乏寿命。本文介绍了ANSYS随机振动分析功能，以及利用该功能，依照Steinberg提出的基于高斯

分布和Miner线性累计损伤定律的三区间法进行ANSYS随机疲乏计算的具体过程。

１．随机疲乏现象普遍存在

在工程应用中，汽车、飞行器、船舶以及其它各种机械或零部件，大多是在随机载荷作用下工作，当它们承受的应力水平较高，工作达到一定时间后，经常会突然发生随机疲乏破坏，往往造成灾难性的后果。因此，预计结构或零部件的随机疲乏寿命是十分有必要的。

2．ANSYS随机振动分析功能介绍

ANSYS随机振动分析功能十分强大，主要表现在以下方面：

具有位移、速度、加速度、力和压力等PSD类型；

能够考虑a阻尼、b阻尼、恒定阻尼比和频率相关阻尼比；

能够定义基础和节点PSD鼓舞；

能够考虑多个PSD鼓舞之间的相关程度：共谱值、二次谱值、空间关系和波传播关系等；

能够得到位移、应力、应变和力的三种结果数据：1s位移解，1s速度解和1s加速度解；

3．利用ANSYS随机振动分析功能进行疲乏分析的一般原理

在工程界，疲乏计算广泛采用名义应力法，即以S-N曲线为依据进行寿命估算的方法，可以直接得到总寿命。下面围绕该方法举例说明ANSYS随机疲乏分析的一般原理。

当应力历程是随机过程时，疲乏计算相对比较繁杂。但已经有大量种分析方法，这里仅介绍一种比较简单的方法，即Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区

间法（应力区间如图1所示）：

应力区间发生的时间

-1s~+1s68.3%的时间

-2s~+2s27.1%的时间-3s~+3s4.33%的时间99.73%

大于3s的应力仅仅发生在0.27%的时间内，假定其不造成任何损伤。在利用Miner定律进行

疲乏计算时，将应力处理成上述3个水平，总体损伤的计算公式就可以写成：其中：

：等于或低于1s水平的实际循环数目（0.6831）；：等于或低于2s水平的实际循环数目（0.271）；

：等于或低于3s水平的实际循环数目（0.0433）；

,,：根据疲乏曲线查得的1s、2s和3s应力水平分别对应许可循环的次数。

综上所述，针对Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法的ANSYS随机疲乏分析