Оценка минимального коэффициента безопасности с использованием метода интерполяционных полиномов Текст научной статьи по специальности «Строительство и архитектура»

УДК 621.3.08

ОЦЕНКА МИНИМАЛЬНОГО КОЭФФИЦИЕНТА БЕЗОПАСНОСТИ С ИСПОЛЬЗОВАНИЕМ МЕТОДА ИНТЕРПОЛЯЦИОННЫХ ПОЛИНОМОВ

© 2014 Ю.Л. Тарасов, С.Н. Перов, С.А. Чернякин

Самарский государственный аэрокосмический университет имени академика С.П. Королёва

Поступила в редакцию 09.12.2014

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

Ключевые слова: коэффициент безопасности, метод интерполяционных полиномов, несущая способность, шпангоут, обечайка, вероятность, безотказная работа

Из-за недостаточной информации об эксплуатационных нагрузках и приближенного представления о значениях механических характеристик материала конструкции, определяющих его сопротивление внешним нагрузкам, основным методом оценки прочностной надежности до настоящего времени является назначение запасов прочности [1, 2]. Значение этих запасов принимаются в зависимости от стабильности условий нагружения, справочных данных о механических характеристиках, уровня технологии и ряда других факторов. Допустимые значения запасов прочности назначают с учётом инженерного опыта создания подобных конструкций. До настоящего времени отсутствуют теоретическое и экспериментальное обоснования составляющих запаса прочности, не учитывается стохастическая природа действующих квазистатических и циклических нагрузок и характеристик используемых конструкционных материалов, что приводит к существенному увеличению металлоемкости конструкций и назначению неоправданно высоких значений коэффициентов безопасности.

Тарасов Юрий Леонидович, доктор технических наук, профессор кафедры космического машиностроения. E-mail: proch@ssau.ru

Перов Сергей Николаевич, доктор технических наук, доцент кафедры космического машиностроения. Email: perov@imi-samara.ru Чернякин Сергей Алексеевич, аспирант

Нормативный подход имеет и недостатки. Во-первых, наблюдается традиционность в проектировании. Из-за этого нередко нормативные материалы превращаются в своего рода тормоз на пути совершенствования показателей металлоемкости конструкции. Во-вторых, наблюдается нарушение системного подхода к проблеме. При этом нормативные материалы не охватывают весь комплекс вопросов, связанных с оценкой потребных значений несущей способности, а ограничиваются регламентацией лишь некоторых из них. В существующих материалах основных отраслей техники отсутствуют явные количественные связи нормативных нагрузок и соответствующих поправочных коэффициентов с потребной надежностью конструкций в целом и их частей. В сущности, при таком подходе уровень надежности конструкций оценивается лишь качественно. В-третьих, при реализации детерминистического нормативного подхода не может быть в полном объеме учтено разнообразие условий эксплуатации конструкций, сочетание различных факторов, статистический разброс механических свойств материала, геометрических параметров, начальная дефектность элементов конструкции. Указанные обстоятельства обусловливают повышение удельного веса вероятностных моделей при обеспечении прочности и надёжности.

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

Эти модели вполне применимы, если нагруже-ние представляет собой единичный дискретный акт или последовательность таких актов и можно исключить из рассмотрения временные эффекты, процессы накопления повреждений и т.п. С некоторыми оговорками эти модели могут быть использованы также для нагрузок, непрерывно развертывающихся во времени, если в расчеты ввести распределение максимальных значений нагрузок на всем рассматриваемом отрезке времени.

Установим связь между показателями уровня прочностной надежности, т.е. вероятности безотказной работы, с традиционным для детерминированных прочностных расчетов коэффициентом безопасности / с учетом рассеивания значений тех параметров, которые используются в вычислениях. Установим связь между уровнем прочностной надёжности (вероятности безотказной работы) и традиционными для прочностных детерминистических расчётов коэффициентами безопасности, а также запасами прочности при учёте рассеяния параметров, входящих в расчёт. Условие прочности для элемента конструкции записывается в виде:

1 =

Я

N1.

(1)

где Ятт - уинимальная несущая способность, ^ртах - максимальная расчётная нагрузка.

Для элементов конструкций летательных аппаратов максимальная расчётная нагрузка представляется в виде

Np =

шах ^ г

тах

(2)

здесь N'max - максимальная эксплуатационная нагрузка,/- коэффициент безопасности.

В статической постановке связь расчётных (нормативных) величин с учётом рассеяния со статическими характеристиками устанавливается с помощью зависимостей:

^Шах = N + ^,

Яшт = Я -

(3)

(4)

Коэффициенты вариации величин записываются следующим образом:

Q - SN

= N'

а - ^

= Я •

(5)

(6)

Рис. 1. Плотности распределения нагрузки/К(х) и несущей способности/^)

Учитывая (3)-(6) получим новое выражение для коэффициента избытка прочности (1):

1 =

Яшт _ Я

_ NШax / {N + aNSN)_

Я -а.^к Я _ 1 Я 1 -ак я

/ (N + aN&NN) fN 1 + aNN^

(7)

Введём условный коэффициент избытка прочности, соответствующий отношению средних величин N и Я :

я

%сп = N •

(8)

Запишем выражение коэффициента избытка прочности п с учётом (7) и (8)

1 -аЯ

1 =/ 1 Х+О^Ш ■

(9)

Здесь N, Я - средние значения нагрузки N и несущей способности Я, aN, ая - доли отклонения расчётных величин ^^ Ятт от соответствующих значений N и Я, выраженные в долях среднего квадратического отклонения, ¿к, -среднеквадратические отклонения величин N и Я . На рис. 1 представлены плотности распределения нагрузки N и несущей способности Я.

Гауссова мера надёжности находится следующим образом:

г- — ~ — SY

—=Я - N S2=s2R+s2N.

С учётом (5), (6) и (8) имеем:

7 = -

Я - N

Я -1

N

N

л -1

I усл

л -1

I усл

2

А N2

N2

Л усл - 1

Лл +

У

(10)

Установим вероятность безотказной работы. При нормальном законе распределения величин N и R вероятность безотказной работы конструкции H можно определить по таблицам функций нормального распределения:

Н = Ф

\ у

= Ф (/)■

(11)

Формулы (9) и (10) дают возможность связать традиционный коэффициент безопасности f с вероятностью безотказной работы Н. Для этого выражение (10) разрешается относительно цусл:

Лусл ~ N ~

1 + .

1 -(1 -удк )2 1 -(Г3„)

1 )2

(12)

По заданным статистическим характеристикам нагрузки N и несущей способности R, учитывающим рассеяние этих параметров, можно решать две задачи: определить вероятность безотказной работы конструкции H (11) и условный коэффициент избытка прочности цусл (12); на основании формулы (9) найти минимальное значение коэффициента безопасности, обеспечивающее заданное значение вероятности безотказной работы H конструкции:

„1 1 - аЯЯ

Г = -л

Л усл 1 + аNN^■

(13)

статистической динамики. Одним из таких методов является метод интерполяционных полиномов [3].

Рассмотрим применение изложенного выше подхода совместно с использованием метода интерполяционных полиномов к задаче об оценке минимального коэффициента безопасности при анализе напряжённо-деформированного состояния стыка днища топливного бака ракеты-носителя с обечайкой, подкреплённого шпангоутом. На рис. 2 представлено поперечное сечение шпангоута и примыкающие к нему участки днища и обечайки.

приняв минимально допустимое значение коэффициента избытка прочности.

Для сравнительно простых систем, для которых параметр поведения системы получается оценить с помощью простых формул, для отыскания минимального значения f используется метод линеаризации. Для сложных вырожденных систем, где связь между характеристиками системы, параметром внешней нагрузки и параметром поведения системы не получается выразить в явном виде, используются методы

Рис. 2. Схема рассматриваемой конструкции

Первым этапом решения этой задачи является выполнение детерминированного расчёта. В качестве исходных данных для анализа напряжённо-деформированного состояния элемента бака ракеты-носителя берутся значения мат. Ожиданий, приведенных в табл. 1 величин. Для расчёта целесообразно воспользоваться методом конечных элементов. В качестве программного пакета можно использовать пакет А№У8, позволяющим в дальнейшем реализовать метод интерполяционных полиномов с помощью удобного встроенного языка программирования АРОЬ. Данный программный пакет позволяет выявить опасное сечение, в котором достигаются наибольшие напряжения (рис. 2).

Далее проводится решение задачи статистической динамики для вырожденных систем методом интерполяционных полиномов. В качестве параметра внешнего воздействия берётся давление наддува, действующего в баке ракеты-носителя. В качестве свойств системы принимаются модуль упругости материала E, толщина обечайки ¿об., толщина днища ¿дн., радиус обечайки Rоб., геометрические характеристики шпангоута А и /2. За параметр, характеризующий предельное состояние конструкции при котором не происходит её разрушения принимается условный предел текучести а0,2. Законы

распределения и числовые характеристики указанных случайных величин приведены в табл. 1.

25.7204 73.739Е 131.751 Ж 131.763 237.771 _52.2341_105.246_153^57_211.265_264.28

Рис. 2. Опасное сечение элемента конструкции бака ракеты-носителя

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

здесь о1>о2>о3 - главные напряжения.

Для вычисления минимального значения коэффициента безопасности, обеспечивающего

Следует иметь ввиду, что это значение коэффициента безопасности обеспечивает тот уровень вероятности безотказной работы конструкции, который соответствует гауссовой мере надёжности у. В рассматриваемом примере она взята равной 3,72, что соответствует уровню вероятности безотказной работы Н=0,9999.

заданный уровень надёжности, методом интерполяционных полиномов следует выполнить следующие шаги:

1. Определяются законы распределения случайных входных величин и их вероятностные характеристики, а также задаётся потребный уровень надёжности (задаётся Гауссова мера надёжности у, соответствующая потребному уровню надёжности Н);

2. Для каждой конкретной реализации входных случайных величин с использованием программного пакета А№У8 определяется напряжённое состояние в наиболее опасном сечении. Выполнив q решений, находятся q реализаций выходного параметра;

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

4. Определяется условный коэффициент безопасности в соответствии с формулой:

^ = 11 + [1"(1" -

5. Определяется коэффициент безопасности, полагая ц=1 по формуле:

Коэффициенты „ и ао1ЬЕе!:, характеризующие степень рассеяния случайных величин в долях среднего квадратического отклонения взяты по единице.

Изложенный алгоритм достаточно просто реализуется в том же программном пакете А№У8, в котором производится расчёт

Таблица 1. Параметры распределения

Случайная Закон рас- Математическое Среднее квадратиче-

величина пределения ожидание ское отклонение

V нормальный = 0,45 МПа 5р = 0,0067 мш

Е нормальный (Е) = 7 ■ 104 МПа 5е = 0,1 (Е) МПа

&0.2 нормальный (Оол) = 300 МПа .;_-.. = : : МПа

8*6. равномерный (6о5) = 5,5 мм Ч* ^ мм

Ят равномерный <5дн.) = 5,3 мм = 0,2 --■- мм

¿1 равномерный = 163,8 мм 5'1 = °'3ММ

к равномерный {У = 53,6 мм 5!2 = 0'3 ММ

^оВ. равномерный <Яс5.> = 1324,7мм = 0,3 ■ ■: : мм

напряжённо-деформированного состояния конструкции. Это возможно с помощью средств встроенного языка программирования APDL. Представленный выше алгоритм можно изобразить в виде блок-схемы, показанной на рис. 3. В табл. 2 приведены результаты, полученные

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

ИСХОДНЫЕ ДАННЫЕ К РАСЧЕТУ

Задание Гауссовой меры надёжности, обеспечивающей требуемый уровень вероятности безотказной работы

—► Оценка условного коэффициента запаса прочности

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

Рис. 3. Блок-схема Таблица 2. Результаты расчёта

Число узлов Напряжения по Мизесу Коэффициент безопасности

математическое ожидание, МПа среднее квадратиче-ское отлонение, МПа

2x2x2x2x2x2x2 264,1513 18,8945 1,05895

3x3x3x3x3x3x3 264,3782 19,0878 1,05905

Выводы:

1. Необходимо иметь достоверную и наиболее полную информацию по вероятностным характеристикам свойств системы, параметров

предельного состояния и параметров внешних воздействий, как на этапе изготовления изделия, так и на этапе его эксплуатации.

2. Обеспечивать заданную вероятность безотказной работы можно двумя путями: а) добиваться меньшего рассеяния значений параметров свойств системы, параметров предельного состояния и параметров внешних воздействий (т.е. значения случайных величин должны лежать в меньшем диапазоне коэффициентов а), что в свою очередь приводит к ужесточению условий эксплуатации и технологии изготовления изделия; второй путь заключается в снижении параметра поведения системы, либо в повышении параметра предельного состояния. Добиться этого можно изменяя конструктивные параметры, физико-механические свойства и/или параметры внешнего воздействия.

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

удобством численной реализации данного метода.

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

СПИСОК ЛИТЕРАТУРЫ:

1. Надёжность и эффективность в технике: Справочник в 10 т. Т. 1: Методология. Организация. Терминология; под ред. А.И. Рембезы. - М.: Машиностроение, 1987. 224 с.

2. Надёжность и эффективность в технике: Справочник в 10 т. Т. 2: Математические методы в теории надёжности и эффективности; под ред. Б.В. Гне-денко. - М.: Машиностроение, 1987. 280 с.

3. Перов, С.Н. Обеспечение надёжности трубопроводных систем / С.Н. Перов, С.И. Аграфенин, Ю.В. Скворцов, Ю.Л. Тарасов. Самара: СНЦ РАН, 2008. 246 с.

ESTIMATION THE MINIMUM SAFETY COEFFICIENT WITH USING THE INTERPOLATION POLYNOMS METHOD

© 2014 Yu.L. Tarasov, S.N. Perov, S.A. Chernyakrn

Samara State Aerospace University named after the academician S.P. Korolyov

The algorithm allowing to connect the set level of reliability of a design with safety coefficient is given in article. The algorithm is stated methodically and is easily realized at the solution of practical tasks. At the solution of practical task on stress-strain state of interface of cylindrical shell with spherical bottom through the elastic bulkhead the method of interpolation polynoms was used. The minimum value of safety coefficient for the considered task is received. The analysis of influence of dispersion the parameters of behavior, properties of system and parameter of limit state on value of minimum safety coefficient is carried out. Ways of ensuring the set level of probability of failure-free operation are specified.

Keywords: safety coefficient, interpolation polynoms method, bearing capacity, bulkhead, shell, probability, failure-free operation

Yuriy Tarasov, Doctor of Technical Sciences, Professor at the Space Mechanical Engineering Department. E-mail: proch@ssau.ru

Sergey Perov, Doctor of Technical Sciences, Associate Professor at the Space Mechanical Engineering Department. E-mail: perov@imi-samara.ru Sergey Chernyakin, Post-graduate Student