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

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

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

et,500,shell63 press=100e6 amesh, all esla, s nsla,s,1

! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将 ! 所有节点坐标系旋转到局部坐标系下. *get,enmax,elem,,num,max dofsel,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)*press f_tot=p_e*ae esel,s,elem,,i nsle,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 *endif

esla,s *enddo aclear,all

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

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

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

et,500,shell63 press=100e6 amesh,all esla,s nsla,s,1

! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将 ! 所有节点坐标系旋转到局部坐标系下. *get,enmax,elem,,num,max dofsel,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)*press f_tot=p_e*ae esel,s,elem,,i nsle,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 *endif esla,s *enddo aclear,all

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

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

--------------------------------------------------------------------------------

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

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

cyl4,,,0.5,,,,3 vsweep,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,then t1(k)=1000*sin(ny(k)) *else t1(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 ~+1s 68.3%的时间

-2s ~+2s 27.1%的时间 -3s ~+3s 4.33%的时间 99.73%

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

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

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

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

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

综上所述，针对Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法的A NSYS随机疲劳分析的一般过程是：

(1) 计算感兴趣的应力分量的统计平均频率(应力速度/应力)；

(2) 基于期望(工作)寿命和统计平均频率，计算1 ，2 和3 水平下的循环次数 、 和 ；

(3) 基于S-N曲线查表得到 、 和 ； (4) 计算疲劳寿命使用系数。

显然，根据其他随机疲劳分析方法和ANSYS随机振动分析结果，我们还可以进行许多类似的疲

劳分析计算。