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

1、在任意面施加任意方向任意变化的压力在某些特殊的应用场合,可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷,当然,这在一定程度上可以通过ANSYS表面效应单元实现。如果利用ANSYS的参数化设计语言,也可以非常完美地实现此功能,下面通过一个小例子描述此方法。!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!finish/prep7et,500,shell63press=100e6amesh, allesla, snsla,s,1! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将! 所有节点坐标系旋转到局部坐标系下.*get,enm

2、ax,elem,num,maxdofsel,s,fx,fy,fzfcum,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

3、!单元中心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,nnf,nelem(i,j,fx,f_n !压力的作用方向为X方向! f,nelem(i,j,fy,f_n !压力的作用方向为Y方向! f,nelem(i,j,fz,f_n !压力的作用方向

4、为Z方向*enddo*endifesla,s*enddoaclear,allfcum,repl !将力的施加方式还原为缺省的"替代"dofsel,allallsel说明:本信息在任意面施加任意方向任意变化的压力在某些特殊的应用场合,可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷,当然,这在一定程度上可以通过ANSYS表面效应单元实现。如果利用ANSYS的参数化设计语言,也可以非常完美地实现此功能,下面通过一个小例子描述此方法。!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!finish/prep7et,500,shell6

5、3press=100e6amesh,allesla,snsla,s,1! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将! 所有节点坐标系旋转到局部坐标系下.*get,enmax,elem,num,maxdofsel,s,fx,fy,fzfcum,add !将力的施加方式设置为"累加",而不是缺省的"替代"*do,i,1,enmax*if,esel,eq,1,then*get,ae,elem,i,area !此命令用单元真实面积,如用投影面积,请用下几条命令! *get,ae,elem,i,aproj,x !此命令用单元X投影面积,如用

6、真实面积,请用上一条命令! *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

7、,nnf,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,allfcum,repl !将力的施加方式还原为缺省的"替代"dofsel,all-ANSYS5.7版本具有函数加载功能,可以很方便地在模型表面施加函数变化的各种载荷,在ANSYS5.6版本中,也可以通过变通的方式来实现此功能,其思路是:首先选定所要施加函数变化表面载荷的表面上的节点,利用ANSYS的参数数组和嵌入函

8、数知识写一简单的命令流,定义好相应节点位置的面载荷值,然后通过在节点上施加面载荷来完成。下面以在一圆柱表面施加函数变化载荷为例:/prep7et,1,45cyl4,0.5,3vsweep,allasel,s,loc,y,0.01,1nsla!*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(1sf,

9、all,pres,0利用ANSYS随机振动分析功能实现随机疲劳分析-ANSYS随机振动分析功能可以获得结构随机振动响应过程的各种统计参数(如:均值、均方根和平均频率等,根据各种随机疲劳寿命预测理论就可以成功地预测结构的随机疲劳寿命。本文介绍了ANSYS随机振动分析功能,以及利用该功能,按照Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法进行ANSYS随机疲劳计算的具体过程。在工程应用中,汽车、飞行器、船舶以及其它各种机械或零部件,大多是在随机载荷作用下工作,当它们承受的应力水平较高,工作达到一定时间后,经常会突然发生随机疲劳破坏,往往造成灾难性的后果。因此,预测结构或

10、零部件的随机疲劳寿命是非常有必要的。ANSYS随机振动分析功能十分强大,主要表现在以下方面:具有位移、速度、加速度、力和压力等PSD类型;能够考虑a阻尼、b阻尼、恒定阻尼比和频率相关阻尼比;能够定义基础和节点PSD激励;能够考虑多个PSD激励之间的相关程度:共谱值、二次谱值、空间关系和波传播关系等;能够得到位移、应力、应变和力的三种结果数据:1s位移解,1s速度解和1s加速度解;在工程界,疲劳计算广泛采用名义应力法,即以S-N曲线为依据进行寿命估算的方法,可以直接得到总寿命。下面围绕该方法举例说明ANSYS随机疲劳分析的一般原理。当应力历程是随机过程时,疲劳计算相对比较复杂。但已经有许多种分析

11、方法,这里仅介绍一种比较简单的方法,即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线性累计损伤