Материал: Прочность, устойчивость, колебания. Т. 3

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

516 Статистические задачи колебаний и устойчивости

Квазистатические методы основаны на применении хорошо извест­ ных формул теории вероятностей. Эти методы с успехом могут быть применены к тем задачам, в которых случайные факторы описываются при помощи конечного (практически не слишком большого) количества случайных величин. Для того чтобы использовать квазистатические методы, необходимо располагать решениями соответствующих детерми*

Рис. 3. Классификация методов решения задач статистической дин

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

Корреляционные методы основаны на использовании связи между корреляционными (или моментными) функциями «входных» параметров (например, нагрузок) и «выходных» параметров (прогибов, внутренних усилий, напряжений). Эти связи могут выражаться как при помощи дифференциальных и интегро-дифференциальных уравнений, так и — в простейших случаях — при помощи конечных соотношений. Спек­ тральный метод и метод канонических разложений занимают промежу­ точное место между корреляционными и квазистатическими методами. Область применения корреляционных методов — задачи, в которых

Квазистатические методы статистической динамики

517

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

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

КВАЗИСТАТИЧЕСКИЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧ СТАТИСТИЧЕСКОЙ ДИНАМИКИ

Предварительные замечания. Квазистатические методы основаны на формулах элементарной теории вероятностей н не требуют привле­ чения теории случайных процессов. Применение квазнстатических методов предполагает предварительное решение соответствующих детерминистических (статических или динамических) задач [3, 4, 25, 28, 501. Приложения квазнстатических методов к различным задачам устойчивости оболочек даны в работах (10, 12, 21—23]. Аналогичные методы применяют в теории автоматического управления и радиотех­ нике под названием «методы безынерционных преобразований» |27].

Условия применения квазнстатических методов. Эти условия со­ стоят в следующем [41:

1) входные параметры, влияющие на поведение конструкции, могут быть описаны при помощи конечного числа случайных величин

и1р иг> . ., ит, <71, <7г,

•* <7г с известной совместной плотностью ве­

роятности р (и1ши2,

., «ш, <7г»

• •* <7/0; здесь иа — внутренние

параметры, характеризующие форму и механические свойства конструк­ ции; <7а — параметры, характеризующие нагрузку и другие внешние условия эксплуатации или эксперимента;

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

чии Ух, о2. • •> ^л; 3) указанные две группы параметров связаны между собой извест­

ными детерминистическими зависимостями

 

 

«а = ^а(«1 .

<71. 9г.

Чг). (а = 1, 2,

л),

(1)

которые предполагают однозначными и кусочно-непрерывными. Определение совместной плотности вероятности для выходных пара­

метров. Плотность вероятности р (ь>1, о2, уя) определяют по из­ вестным формулам для плотностей вероятности функций случайных

518 Статистические задачи колебаний и устойчивости

величин. В частности, если п ^ т и если соотношения (1) допускают решение

иа = и а (Щ, о2,

оп, ип+1, .

ит, Ях, Яг,

. .. Яг)

= 1 , 2 ,

л),

(2)

где Уа — однозначные и непрерывные функции переменных »1, у2, . . .

I>л, то решение дается формулой [4]

Р (»1, ^2.

»л) =

= 1

 

Уп, «л+1.

Ях.......... ?г) х

 

 

 

„ I

д ( У 1, Г/,

Ц»)

 

А I

^(^1. ^21

дип+1-'-(1ит (1ях'--с1Яг.

•. ”п)

 

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

 

 

О/!

дУ2

 

дУп

 

 

дох

дох

' * ^1»!

б {Ух. У*

Уп)

дУх дУ2

 

 

до2

до.,

*

до..

д(о1г о2,

0„)

 

 

 

 

 

 

дУх дУ2

 

дУп

 

 

доп

до,,

*‘"

доп

(3)

(4)

Формула (3) может быть легко обобщена на случаи, когда функции У1, и 2, . . ., Уп не являются однозначными (тогда область интегриро­ вания разбивается на ряд подобластей), на случай, когда эти функции

являются

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

для

и1э у2>

оп будет содержать дельтаобразные выбросы) и т. д.

Если

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

усреднением уравнений (1).

п < т -{- г.

Нетрудно переписать формулу (3) для случая, когда т <

Если же п^> т + г, то часть из параметров у1, о2, .

оп будет, оче­

видно, функционально зависимой. Тогда можно принять, что п = т + г

и искать

совместную

плотность вероятности для

параметров о1У

у2* • • •»

^т+г-

 

 

 

 

 

' Вместо формул (2) имеем соотношения

 

 

 

иа = Уа(их,

о2ш

от+г) (а = 1,

2, . .

т)\

 

7а=<?а(ы!»

о2,

от+г)

(а = 1,

2,

г).

Для совместной плотности р (цх, о2,

;;|+г) имеем формулу

 

 

Р&Х,

о2,

от+г) =

 

 

 

= Р(<Л. . .

Утш

 

<?,) X

 

 

д(Ух, . .

Ущ,

 

фг)

(5)

 

X

д (уь

о2, ..

ут+г)

 

 

 

 

 

Квазистатические методы статистической динамики

519

Определение вероятности отказа. Допустим, что области нарушения нормальных условий эксплуатации принадлежат точки в пространстве

Ох, у2» ■•

•»

Для которых V («!, о2,

у„) >

0. Тогда вероятность

отказа (см. также гл. 8, т. 1) определяют по формуле

 

 

 

р = [ • • • [ Р (и1.

“п)

4ог

а«п

(6)

 

 

 

^>0

 

 

 

 

 

Эту вероятность можно найти также иначе

[4]. Если области V (г^,

у2,

•.

Уп) !> 0 соответствуют

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

и1, и2, . .

ит, щ1%

Яг*

•. Яг, область 'К, К , ы2, .

ит,

Я1. Яг*

•. Яг) >

0, то

 

 

р ==1 ‘ “ ]

 

Яи

Яг’

•• ^ *

 

 

 

 

ЧГ1>0

 

 

 

 

 

 

 

 

Х4щ(1и2.

.й и т йЯ1.

. &Яг-

(7)

Задачи, связанные с обработкой и интерпретацией опытных данных. Пусть независимыми случайными величинами являются внутренние параметры системы и1г и2, ., ит. Требуется найти совместную плот­ ность вероятности для предельных значений параметров 01э ^2, . . щг.

Будем обозначать эти предельные значения через </2* • » Яг» Допустим, что известны детерминистические соотношения между пара­

метрами <7^ <72,

<7Г и внутренними

 

параметрами

конструкции иь

м2,

-1 ит) а также совместная плотность вероятности

р (мх, и2,

 

.

 

ит). Пусть далее г<^т

и, кроме того,

существуют однозначные

и

непрерывные

зависимости

 

 

 

 

 

 

 

“а

=

и а (?’>

.......... я\, иг+1............ ит)

(о = 1, 2,

г);

(8)

тогда совместная плотность вероятностир (<7*,

?*) может быть

найдена по формуле типа

(3):

 

 

 

 

 

 

 

 

 

Р (я\,

Ч1

Яг) =

 

 

 

 

 

 

 

\р (1 > I. У*.

 

 

 

Пт) X

 

 

 

 

а(Уь

у г,

.. •• Уг)

 

диГ4-1. .. дищ-

 

 

(9)

 

 

 

а(<7*, я\.

«’)

 

 

 

 

 

 

 

Некоторые ограничения, накладываемые

на зависимости (9), могут

быть легко сняты. Случай, когда однозначная зависимость отсутствует, рассмотрена книге [4].

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

520

Статистические задачи колебаний и устойчивости

по формулам типа (3), (5) и (9). Пусть,

например, г = т, а зависимости

=

<?а(“1. “2-

■’ ит) однозначны

и непрерывны. Тогда

Р («].

V

А

ит) = р (<?;.

<?;.

<?ш) X

^ 2'

Ф/л

(10)

д(ии «2,

 

ит)

 

 

Распределение вероятностей для предельных сил р

д2,*

д*)

можно получить, исходя из формулы (7). Допустим,

что параметры

Я\*

•» Яг фиксированы и равны дх, <72,

цг

соответственно.

Тогда получим следующее выражение для условной вероятности отказа:

Р ( * | ^1» 72» 7Г) —

 

- И

Р ( « 1 .

ит) йих с1ий. . йит\ (11)

 

 

 

(«Г

и2‘

т V '2-

о >

 

 

Формула (11) дает по существу совместную функцию распределения

для

д2,

Яг■ Для

совместной

плотности

вероятности имеем

формулу

 

 

 

 

 

р (ч1

я\,

дгР ( * \д \ ,д \,

■■,<?*)

 

 

 

( 12)

 

 

 

 

дЧ\дя:2...дч,

Примеры применения квазистатических методов. Ряд работ [3, 12, 22, 23] посвящен следующей задаче: пластинку или оболочку нагружают внешними силами, заданными с точностью до одного общего множи­ теля — параметра д. Этот параметр весьма медленно (квазистати­ чески) и монотонно возрастает от нуля до некоторого конечного значения. Требуется найти распределение параметров деформации

(обычно — обобщенных

координат,

характеризующих

нормальный

прогиб), достигаемое к

концу процесса нагружения.

В статье [3]

рассмотрена задача о распределении

вероятностей полного прогиба

упругой пологой цилиндрической панели со смещающимися кромками, сжатой осевыми силами интенсивностью д. Параметр начального прогиба считают случайным, параметр нагрузки — детерминированным. Вероятность «хлопка» для той же задачи вычислена в статье [3] в пред­ положении, что начальные прогибы подчиняются симметричному нор­ мальному распределению со стандартом аи. Эта вероятность показана на рис. 4 как функция нагрузки. Здесь Р (* ) — вероятность «хлопка»; д — величина осевого усилия.

Задача о распределении критических сил подверглась наиболее подробному исследованию (3, 12, 22, 23]. В простейшем случае, когда

Источник: https://studfile.net/preview/19992756/