Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1200

.pdf
Скачиваний:
6
Добавлен:
15.11.2022
Размер:
11.69 Mб
Скачать

Тогда

интеграл

примет вид

 

 

 

 

*

 

 

 

 

*2

 

К —,Г

+ 2 . + I

и2

 

 

 

J

j

Р

 

2 du9

“ - ^ е х р

 

 

 

2 а ?

 

 

 

 

 

- 2 . + I

 

где обозначено

 

 

 

 

 

 

1

( ±v*2 — a2

у * — а , ^

 

 

а± 2.+1

 

 

 

Замечая, что

е 2 d u ~ Ф(ы),

— оо

где Ф(и) — интеграл Лапласа, получим

- у2

После вычисления остальных интегралов приходим к формуле

М Г ) = -£Мехр

 

(v* ~ ai)

 

 

2 л

 

2 а 2

 

+ ехр

( » ;+ « |) 8

[О*

+

 

2 а ?

 

 

 

 

 

 

 

 

+ £ | е * Р

 

(oj—а,) 2

Ф (“—1. +г)] +

 

2а?

[Ф («+,<+2)

 

 

 

 

 

+

ехр Г

 

( р 2 ~ Ь а 2

[«* («+,. _а) —«»

 

 

L

 

2°!

 

 

Здесь использованы

 

обозначения

 

 

1

 

J i z f f — р ± 1 * Г ^ ) ; aj = IL

U±j,±k

' у j — ^2

°у

ff/t /

' <v

 

 

 

(140)

(141)

Для определения функции надежности достаточно подставить выра­ жение (140) в формулы (132) или (133)*.

* О приложении этих формул к стационарным процессам см. подстрочное примечание на стр. 200.

При некоторых дополнительных предположениях формула (140) может быть заменена соотношениями, более удобными для расчетов. Пусть, например, процессы v-^t) и v2(t) — центрированные и стохас­ тически независимые. Тогда аг = а2 р = 0. Формула (140) прини­ мает вид

(Г) =

— ехр [

1— ] Фх

 

 

 

2®?

\

'7777777?

'//////, W

/ / / / A

 

 

1

 

 

\\

я» X

W

------- ^

0

t

 

*

^

 

 

777777?.'//////‘ Х<у’/ / / / /

$- — 2i 1

Рис. 67

+

ШО

 

 

 

 

(142)

ехр

 

 

 

JI

 

 

 

<*1

Здесь

 

Фх(и) — интеграл

 

Лапласа,

представленный в форме*

 

 

 

 

 

1

X2

 

Фх(и) =

2

dx.

 

/ 2 л

 

 

 

Заметим, что имеет место нера­ венство

ф (« +2. + 1) - ф (« _ 2>+1) < 1

и т. д. Отсюда вытекает следующее неравенство для среднего числа по­ ложительных пересечений границы Г:

 

 

 

1

»

1 3о

 

(Г )< - ^ - ехр------

>

 

 

 

 

-с;

 

Ч

 

 

 

2ат

 

 

 

 

 

 

I

ш 2

ехр

2~ а2)2

 

^

2я

2а2

 

 

го1

+ ех р

J

+ ехр

К + 'ч Н >+

2о2

(а2+ аг)2 1 (143)

2а?

Очевидно, что в правой части стоит сумма среднего числа пересе­ чений линий v x = ± v \ и v2 = ±v*2, т. е. сумма среднего числа выбросов

процесса за пределы полос —v \ < v x < v \ и —v? < v2 < V2 (рис. 67). Если выбросы — редкие события, а корреляция между компонентами Vi(/) и v2(t) — достаточно слабая, то правая часть неравенства (143) может быть использована для приближенной оценки среднего числа выбросов из области й 0- Таким образом, приходим к формуле

v+ ( D » | l ехр

(»;-«!)* ] + ехр

К + Д.)Ч

 

 

2а?

 

2о\

 

 

 

 

/

*

\2

 

(i»3 + g2)2

2 -

д2)

+ ехр

 

2а|

(144)

 

 

2а2

* См. примечание на стр. 155.

Формулы (140), (142) и (144) могут быть обобщены на случай про­ извольного числа измерений. Ограничимся тем, что вычислим среднее число положительных пересечений поверхности /z-мерного паралле­ лепипеда стационарным гауссовским процессом v(/), полагая компо­ ненты этого процесса стохастически независимыми. Итак, пусть до­ пустимая область Q0 задана в виде

vV < v k < vk (к==1' 2> л)«

а совместная плотность вероятности процесса v(/) и его первой произ­ водной — в виде

P(V’

Здесь использованы обозначения

 

 

Ph (vh) = " т = — ехР

(\ - ° к )2

к/

Y2n<3,

У

Ч

 

 

 

 

 

 

vl

 

K - 5 5 .,.,“

P V

где ah — математические ожидания компонентов vh(t)\ а\ — их диспер­ сии; щ — эффективные частоты. Подставляя эти выражения в форму­ лу (130), получим

М Г)= 2

М °1) + М °Г)] $ Pn-\-k(vk)vkdvk J llr

 

 

k= 1

 

 

‘ ^ ^

 

 

Отсюда после вычисления интегралов найдем

 

 

 

м г>= 2

-S-iexp

К -» * ) 1

 

 

 

k= 1

 

Ч

 

 

 

 

П ф |

V,

— а,

 

+ ехр

 

vi ~ ai \ _ ф / ui

~ ui

(145)

Ч

 

 

 

 

 

 

i*k

 

 

 

Заметим, что для высоконадежных систем можно приближенно принять, что

п ф

V,

— а,

- ф

1.

k

Тогда придем к приближенной формуле

 

п

К - « * ) 2

 

(1>

exp

+ exp

' *= I 2я

2 а f

н

дающей, как и формула (144), оценку сверху для среднего числа пере­ сечений. Каждый из членов суммы, стоящей в правой части формулы (146), представляет собой среднее число выбросов в единицу времени

 

из слоя v*k

<Vk<Vk. Будем в дальнейшем

 

называть формулы типа (144) и (146) фор­

 

мулами полосового приближения.

 

 

Аналогично могут быть

получены фор­

 

мулы для

 

случая, когда область Q0 яв­

 

ляется кругом и, вообще, n-мерным ша­

 

ром.

Некоторые

вычисления,

основанные

 

на асимптотических формулах

типа

Кра­

 

мера, были проделаны А. М. Шепутисом

 

1119].

приближенного вычисления на­

 

Для

 

дежности

 

в

случае областей сложной кон­

 

фигурации

 

можно

рекомендовать

метод

мажорантных оценок. Допустим, что

в

пространстве

качества взяты

такие две области Qx и й 2, что &i cz

 

cz й 2. Например,

граница 1\

области Q1 вписана в границу

Г, а

граница

Г2 области

Q2 описана

вокруг Г. Тогда для функции

надежности

имеем

 

 

 

где Px(t)

и P 2(t) — функции надежности для областей

и Q2 соответ­

ственно.

В качестве областей

и й 2 целесообразно взять такие, для

которых функции надежности и средние числа выбросов легко вычис­ ляются (рис. 68).

§ ШЛО. Применение теории надежности к расчету оптимальной виброзащиты оборудования

Вопросам защиты приборов и оборудования от вибрационных воз­ действий посвящена обширная литература (см. например, [54]). Если внешнее воздействие является гармоническим, то задача о выборе па­ раметров виброзащитного устройства решается сравнительно элемен­ тарными средствами. Но в реальных условиях вибрационные воздей­ ствия обычно носят случайный характер. Решение задач об оптималь­ ной виброзащите при случайных внешних воздействиях встречает ряд аналитических и принципиальных трудностей.

В теории оптимальной виброзащиты используются критерии, ана­ логичные критерию минимума среднеквадратической ошибки в теории автоматического управления. Например, ставится условие, чтобы средний квадрат перемещения защищаемого объекта относительно

основания, средний квадрат абсолютного ускорения объекта и т. п. при­ нимали минимальные значения. Оставляя пока в стороне вопрос о тех­ ническом истолковании подобных условий, заметим, что они являются источниками некоторых далеко идущих затруднений. Одним из первых встречается следующее затруднение. Если подбирать параметры виброзащитного устройства из условия минимума среднеквадратического абсолютного ускорения, то придем к тривиальному решению: жест­ кость подвешивания и потери в демпфере должны быть минимальны. Однако при этом получаются недопустимо большие относительные пе­ ремещения объекта. Если же минимизировать среднеквадратическое относительное перемещение, то становятся недопустимо большими пе­ регрузки на объекте. Из этого затруднения выходят, дополняя усло­ вие минимума одного параметра ограничением, которое накладывается на другой параметр. При этом в зависимости от того, что минимизирует­ ся — средний квадрат ускорения или средний квадрат относительного перемещения, получаются различные решения.

Следующая трудность встречается, когда оказывается, что опти­ мальные параметры линейной виброзащитной системы зависят от ин­ тенсивности воздействия. В реальных конструкциях интенсивность вибрационного воздействия редко остается постоянной в течение сро­ ка эксплуатации конструкции. Линейные системы, параметры кото­ рых подобраны применительно к некоторому уровню воздействия, пе­ рестают быть оптимальными при изменении этого уровня. Для выхода из создавшегося затруднения приходится интерпретировать линейную систему как результат статистической линеаризации некоторой не­ линейной системы.

Эти трудности усугубляются, если объект обладает несколькими (тем более — бесконечным числом) степенями свободы и если внешнее воздействие является нестационарным случайным процессом. Пред­ ставляется, что эти трудности не присущи внутренне теории вибро­ защиты, а скорее являются следствием неудачного выбора критерия для оптимизации. Инженер требует от системы виброзащиты, чтобы она обеспечивала^надежное функционирование объекта. Конечно, сред­ ние квадраты перемещений и ускорений на объекте в некоторой степе ­ ни характеризуют условия надежного функционирования. Но с тех­ нической точки зрения было бы правильнее минимизировать вероят­ ность того, что за время эксплуатации объекта его параметры хотя бы раз выйдут из области допустимых значений. Это эквивалентно тре­ бованию, чтобы мера надежности системы принимала максимальное значение. Оптимизация по надежности, будучи более естественной и обоснованной, в то же время снимает трудности, возникающие при применении более частных критериев. Этот метод оптимизации приме­ ним как к линейным, так и нелинейным системам с произвольным чис­ лом степеней свободы; не накладывается также ограничений на стохас­ тическую природу внешних воздействий.

Общая постановка задачи о проектировании виброзащиты форму­ лируется следующим образом. Пусть некоторая механическая систе­ ма с конечным или бесконечным числом степеней свободы прикрепля­

ется к основанию при помощи конечного числа опор. В общем случае свойства этих опор неизвестны; тогда говорят о выборе структуры виб­ розащиты. Однако чаще система виброзащиты ищется в классе прос­ тых линейных связей, содержащих упругость и вязкое трение. В этом случае неизвестными параметрами являются координаты опор, коэф­ фициенты их жесткости и вязкости. Под действием внешних сил или ускорений, сообщаемых основанию, в системе возникает некоторое виб­ рационное поле. Из технических соображений выбирается система па­ раметров качества и допустимая область в пространстве параметров ка­ чества. Далее вычисляется надежность системы как функция неизвест­ ных параметров виброзащиты. Последние находятся из условия, что­ бы надежность, достигаемая к некоторому моменту времени, была мак­ симальна. Как правило, надежность как функция параметров вибро­ защиты не имеет изолированного максимума. Кроме того, указанные параметры обычно могут принимать значения лишь из некоторой огра­ ниченной области. Поэтому проектирование виброзащиты сводится к неклассической задаче оптимизации, которая может быть разрешена лишь численными методами.

Одним из центральных остается вопрос о выборе параметров ка­ чества и области их допустимых значений. Чаще всего требуется, что­ бы максимальные виброперегрузки оборудования не превышали не­ которых предельных значений. Пусть а(г, /) — поле абсолютных уско­ рений в системе; а* — предельно допустимое значение вибрационного ускорения. При этих условиях функция надежности вводится как

P (0 = P [ sup макс || а (г, т ) ||< а „ 1.

(147)

Если изолируемая система является абсолютно твердым телом, то максимум абсолютного ускорения достигается в одной из точек, при­ надлежащих поверхности тела. Если тело — многогранник, то сле­ дует проверить каждую из его вершин. Во всяком случае соотношение (147) точно или хотя бы приближенно может быть заменено следую­ щим:

P(t) = P\ sup макс||a (r;-, x ) ||< a j .

(U 8)

Таким образом, пространство параметров качества становится много­ мерным эвклидовым пространством, и метод, изложенный в предыду­ щем параграфе, может быть применен для вычисления функциц на_ дежности.

Наряду с ограничениями, накладываемыми на абсолютные уско­ рения, могут быть также поставлены ограничения для относительных перемещений различных точек системы. В соотношения типа (148) мо­ гут входить также проекции ускорений на некоторые направления и т . п. В этих случаях число измерений пространства параметров ка­ чества должно быть увеличено.

Поясним постановку и метод решения задачи на простейшем при­ мере. Будем трактовать объект как систему с одной степенью свободы.

Пусть объект массой М при помощи упругой связи с жесткостью с и вязкой связи с коэффициентом трения k прикреплен к основанию, которое совершает колебания с переносным ускорением а0(/) (рис. 69). Перемещение и(() объекта относительно основания удовлетворяет уравнению

и + 2ги + а1и= —а0(0.

(149)

где е = kl2My азо = с/М. Пусть, далее, для надеж­ ного функционирования объекта требуется, чтобы относительное перемещение и и абсолютное уско­

рение а = а0 + и не превышали по модулю пре­ дельно допустимых значений и* и а* соответственно:

М < « , . |а |< в ..

(150)

oft)

м

»t(t)

*7777777777777777.

Рис. 69

Система виброзащиты будет оптимальной по надежности, если ве­ роятность пребывания системы в допустимой области (150) за время 0 < т < 7 будет максимальной. Таким образом, приходим к крите­ рию оптимизации в виде

РГ sup \и(т )\< и^

sup |а (т )|< а Л =макс.

(151)

|_0<т<Г

J

 

В дальнейшем ограничимся случаем, когда ускорение основания a0(t) представляет собой стационарный гауссовский процесс. Среднее число v_|_(r) выбросов из допустимой области в единицу времени не будет зависеть от времени. Тогда, как это следует из формул (132) и (133), вместо критерия оптимизации (151) может быть принят крите­ рий

V+ (T) = MHH.

(152)

Поскольку допустимая область (150) представляет собой прямоуголь­ ник, то среднее число выбросов г+(Г) находится по формулам (140)

и (144).

Возьмем для процесса а0(0 спектральное представление

00

ао (0 = <ао> + $ А (а>) eml d(0.

— оо

Здесь Л(со) — обобщенная случайная функция (спектр процесса). Представляя в аналогичной форме процесс u(t)

00

ы (*) = <«> + 5 и((й)еш йш

— оо

 

и используя уравнение (149), получим

 

« ■ > = - - ^

о 53>

“ о

 

8В. Зак. 1481

Далее стандартным методом (§ 1.7) найдем связь между спектральной плотностью <Sa0(co) входного процесса a0(t) и моментами второго поряд­ ка выходных процессов:

Г

12>

J

|(0ц + 2('е(й —со2!2

—гш'

'

* | й)% + 2ш о|2 S0o(<o)

 

| COQ-f-2teco —со2 I2

oo

S Go(со) со2 day

'“l,>= I

 

(154)

| COQ-J-2/80)— Cl)2 |2

 

„2

 

 

cog + 2/eco |2 5flo (со) со2 rfco

< M 2 > =

| со2 +

2/eco— со2 |2

 

<u* n) — —

(со2 -f- 2/eco) S a^ (со) day

Icog+

2/eco —со2 12

 

Чтобы провести вычисления далее, необходимо задаться спект­ ральной плотностью S (J0( CQ) . Пусть, например, ускорение основания a0(t) является узкополосным случайным процессом с несущей часто­ той 0. Его спектральную плотность аппроксимируем при помощи дель­ та-функции, т. е. примем, что

5ПоИ ^ б ( М - е ) ,

(155)

где сто — дисперсия ускорения a0(t). Подставляя (155) в формулы (154), получим

~

°о

1

ю

а

а,

 

j

1+ И '

 

02 \2

 

 

 

У

\

0,0)

 

 

 

е '

(156)

 

)!

1

СОо >

е \2

У

[— •

\

 

Wo /

При этом через у обозначен декремент затухания собственных коле­ баний, т. е.

2ле

На рис. 70 представлен график зависимости среднего квадратичес­ кого смещения ои от собственной частоты со0 при различных значениях у. При этом принято, что 0 = 1000 секг1>сго= 10 0 м-сек . Здесь же

нанесена кривая для среднего смещения ( и), построенная по формуле (153) (принято, что <я0> = 1° М’Секг2). На рис. 71 даны кривые для сред­ него квадратического ускорения оа. Эти кривые аналогичны зависи­ мостям динамического коэффициента для линейной системы с одной

степенью свободы при гармоническом возбуждении. Оптимальные ус­ ловия виброзащиты следует искать в области малых собственных час­ тот (о0, имеющих порядок

со0

I < Д р > 1

ы*

 

На рис. 72 приведены результаты вычислений среднего числа по­ ложительных выбросов v^_(F) из области (150) при и* = 10~3 л*, а* = =-20 м-секг2. Сплошные линии соответствуют вычислениям по формуле

полосового приближения (144). При этом = u(t)yv2 =

a{t)yосц = о) =

=

0. Нисходящая кривая дает число выбросов из полосы —и* <; и <

<

а*; восходящие кривые — из полосы—а* <1 а <

а*. Здесь же

показано суммарное число выбросов. Расчеты по точной формуле (140) показывают, что ошибка формулы полосового приближения становит­ ся заметной лишь при достаточно больших числах v+(r).

Наилучшая виброзащита по \Ц-(Г) обеспечивается при исчезающе малом демпфировании и при собственной частоте со0 ж 140 сек-1. Ошибка в выборе частоты со0 на ± 1 0 сект1 дает увеличение числа вы­ бросов примерно на порядок. Поэтому при расчете технических систем необходимо учитывать разброс параметров виброзащитного устройст­ ва, трактуя функции надежности типа (147), (148) и (151) как условные надежности.

Заметим, что ввиду узкополосности процесса надежность Р(^) ЛуЧ. ше оценивать через среднее число выбросов огибающей, т. е. пользо­ ваться формулами типа (105).

Выше был рассмотрен узкополосный случайный процесс на входе системы. Многие реальные воздействия представляют собой широкопо­

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]