异型密封圈计算泄漏量与参数化优化过程仿真.pdf

异型密封圈计算泄漏量与参数化优化过程

1.%2.%3后处理

图325接触压力云图

从接触压⼒云图可以得出，异型密封圈主密封接触压⼒达到18MPa，并且接触部位的接触状态稳定，

没有发⽣断层和突变，说明不会发⽣泄漏。

图326轴承支撑反作用力

根据泄漏量计算公式，需要知道密封圈所受载荷的⼤⼩，此处便使⽤⽀撑轴承的反作⽤⼒间接得到密

封圈所受载荷，⽽密封间隙宽度则需要通过计算接触面积来求得，这⾥便需要使⽤命令来实现：

1set,last!选择计算迭代的最终步

2cmsel,s,pressure!通过名称选择pressure的下半部分的所有节点

!通过名称增选pressure2的下半部分的所有节

3cmsel,a,pressure2

点

4esln,s,1!选择附加到节点的element

!只重新选择类型为ldquo;cid_1rdquo;的接

5esel,r,type,,cid1

触element

6etable,estat,cont,stat!保存已选择单元列表

7esel,s,etab,estat,3!选择接触部分

8etable,c_area,VOLU!选择接触元件面积

9ssum!面积相加

10*get,t_area,ssum,0,item,c_area!将结果存储到t_area

11my_areat_area/4.84!求取密封间隙宽度并将结果输出到my_area

2.%2.%3泄漏量计算

经过上述的仿真流程，得到了密封间隙的宽度和密封圈的径向载荷，然后再使⽤MATLAB平台将公式

编⼊功能函数进⾏泄漏量的计算。以下便是泄漏量计算功能函数的代码：

function[s]LLJS(L,W)%L为密封间隙宽度，W为受到的径向载荷

LL/1000;%换算单位为M

P113.16e6;%内压

P211.76e6;%外压

k0.5;%暴露损失所占比重

A0.2e-3;%振幅M

v220e-6;%运动粘度M2/S

rho930;%密度KG/M3

etav*rho;%动力粘度N*S/M2

nRotation80;%转速rpm

D55e-3;%直径M

N36;%周期个数

VnRotation*pi()*D/60;%线速度M/S

WW*N/D/pi();

hsqrt(0.1602*eta*V*L.^2/W);

f@(x)sin(2*N*x/D)./(A.^2*

(sin(2*N*x/D)).^2+1).^0.5;

s(1-k)/2*N*A*h*V*integral(@(x