Расчет основных характеристик в радиальном ГСП во FLEX.PDE

Предложенный метод расчета реализован в среде Flex.PDE в виде программы BIARING1R.PDE. Программа включает в себя разделы:

- variables (переменная) с переменной temp (безразмерный квадрат давления в зазоре)

- select в данном разделе задается число шагов для вычисления функций задаваемых оператором stages

- definitions в разделе задаются:

Pa- давление на кромке подшипника,

Ps- давление наддува,

Rn- наружный радиус подшипника,

Rv- внутренний радиус подшипника,

Rd- радиус окружности расположения питателей,

dp- диаметр питателей,

Н- величина зазора подшипника,

k1- покозатель адиабаты газа,

ass- скорость звука,

А- площадь питателей,

Gamma- конструктивный параметр,

N- число питателей.

.

- equations в разделе приводится уравнения Рейнольдца при этом правая часть уравнения определяет влияние вращения на распределение давления в зазоре при sigma = 0 вращение отсутствует.

- boundaries в разделе строятся две области Region 1 и Region 2 определяющие смазываемую поверхность опоры УГСП, и задаем граничные условия на кромках подшипника Pa^2,.

- monitors в разделе происходит отслеживание результатов расчетов в процессе работы программы.

- plots в разделе показаны конечные результаты численных значений в виде поля давлений (contour (p)).

- summary показаны значений в определенных точках в численном виде (report (p)).

- history команда связанная с stages.

{BIARING1R.PDE }

{ ************ }

title " Radial GSP "

Variables

temp { dimensionless pressure square in the gap }

Select

! stages=14

Definitions

Par=0.98e5 {davlenie na kromke podshipnika, n/m2 }

! Psr=6*Par! 4.9e5! 4.9e5 {davlenie nadduva, n/m2 }

Rr1=0.075/2! 0.112 { radius podshipnika, m}

dpr=0.9e-3 { diametr pitateley, m}

Lr=0.13/2

Lr1=0.3*Lr

L=Lr/Rr1

L1=Lr1/(Rr1)

cr=26e-6

p=abs(temp)^(1/2)

! L1=0.25

e=0.5! staged(0.01, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45)

! tempi=array(0.72, 0.38, 0.28, 0.25)

er=e*cr

ALF=0.8

H1(f)=cr-er*cos(f)

k1=1.4 { pokazatel adiabaty gaza }

mu=179.2e-7{nc/m2}

ass=341.7 {skorost zvuka, m/c }

Psr=staged (5*Par, 6*Par, 8*Par, 10*Par)

A(f)=if (pi*dpr*dpr/4< =pi*dpr*H1(f)) then pi*dpr*dpr/4 else pi*dpr*H1(f) {plochad pitately }

gamma(f)=(24*mu*ass*A(f)*Rr1*ALF)/(pi*dpr*H1(f)^3*Psr)*((2/(k1+1))^((k1+1)/(2*(k1-1)))) {konstruktivniy parametr }

N=6{chislo pitateley}

Pa=Par/Psr

dp=dpr/(2*Rr1)

sig(f)=if (f< N/2) then(-1) else 1

sig1(f)=if (f< N/2) then(0) else -1

omega=0! staged(100, 1000, 4000, 10000, 30000)

LAMBDA=(6*mu*omega*Rr1*Rr1)/(Psr*cr^2)

nOM=omega*30/pi

Px=(p-Pa)*cos(x)

Preal=p-Pa

h=1-e*cos(x)

F1=(2/(k1+1))^(2*k1/(k1-1))

F2=sqrt(2/(k1-1))*(2/(k1+1))^(-(k1+1)/(k1-1)/2)

q(f)=if (f< =F1) then(1) else (((abs(f))^(1/k1)-(abs(f))^((k1+1)/(2*k1)))^(1/2))

Initial values

temp =0.999

equations { define the heatflow equation }

dy(h^3*dy(temp))+dx(h^3*dx(temp))=2*LAMBDA*dx(h*p)

Boundaries

Region " POVERX" { the outer boundary defines the copper region }

Start (-pi, L)

value (temp) =Pa ^2 {convection boundary}

Line to (pi, L)

periodic(x-2*pi, y) { walk the boundary Counter-Clockwise }

Line to (pi, -L)

value (temp) =Pa ^2

Line to (-pi, -L)

Nobc(temp)