Моделирование частичного закрытия трещиновидных полостей в горящем твердом топливе при действии объемных сил Текст научной статьи по специальности «Физика»

16. Ковалев, В.А. Волновые задачи теории поля и термомеханика / В.А. Ковалев, Ю.Н. Радаев. Саратов: Изд-во Сарат. ун-та, 2010. 328 с.

17. Ковалев, В.А. Распространение связанных 0№11-термоупругих волн в длинном цилиндрическом волноводе / В.А. Ковалев, Ю.Н. Радаев // Вестн. ЧГПУ им. И.Я. Яковлева. Сер. Механика предельного состояния. 2010. № 2(8). С. 207-255.

18. Ковалев, В.А. Прохождение теплового 0№11-волнового сигнала с высокой окружной гармоникой через цилиндрический волновод / В.А. Ковалев, Ю.Н. Радаев, А.Е. Романов // Актуальные проблемы приклад-

УДК 539.375

ной математики, информатики и механики: сб. трудов междунар. конф., посвящ. 80-летию проф. Д.Д. Ивле-ва. Воронеж: Изд. центр Воронеж. гос. ун-та, 2010. С. 173-180.

19. Ковалев, В.А. Элементы теории поля: вариационные симметрии и геометрические инварианты / В.А. Ковалев, Ю.Н. Радаев. М.: Физматлит, 2009. 156 с.

20. Ковалев, В.А. Волновые задачи теории поля и термомеханика / В.А. Ковалев, Ю.Н. Радаев // Математическая физика и ее приложения: материалы междунар. конф. / под ред. чл.-кор. РАН И.В. Воловича и проф. Ю.Н. Радаева. Самара: Книга, 2010. С. 165-166.

МОДЕЛИРОВАНИЕ ЧАСТИЧНОГО ЗАКРЫТИЯ ТРЕЩИНОВИДНЫХ ПОЛОСТЕЙ В ГОРЯЩЕМ ТВЕРДОМ ТОПЛИВЕ ПРИ ДЕЙСТВИИ ОБЪЕМНЫХ СИЛ

В.М. Мирсалимов, Б.Э. Рустамов

Азербайджанский технический университет, Баку, кафедра сопротивления материалов E-mail: mir-vagif@mail.ru

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

Ключевые слова: горящее твердое топливо, трещиновидные полости с концевыми зонами, силы сцепления, объемные силы, контактные напряжения.

Simulation of Partial Closure Crack-Visible Cavities in Burning Solid Fuel Under the Influence of Body Forces

V.M. Mirsalimov, B.E. Rustamov

Azerbaijan Technical University, Baku, Chair of Strength of Materials E-mail: mir-vagif@mail.ru

On the basis of the theory of elasticity mathematical description of the model for the covering crack-visible cavities with end zones in which cohesive forces of material act, in the burning solid-fuel has been performed. It is accepted that the interaction of crack-visible cavity surfaces under the influence of body and surface loads leads to the appearance of overlap zones of their surfaces. The determination of the unknown parameters that characterize the closure crack-visible cavities, is reduced to the solving of the system of singular integrodifferential equations. Using the procedure of algebraization integral equations are reduced to the system of nonlinear algebraic equations which is solved by successive approximations. Normal and tangential contact stresses, the tractions in the bonds, the values of the size of the end contact zones, where the faces of the crack-visible cavities are closed, have been found.

Key words: burning solid fuel, crack-visible cavity with end zones, cohesive forces, body forces, contact stresses.

1. ПОСТАНОВКА ЗАДАЧИ

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

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

Рассматривается напряженно-деформированное состояние среды, ослабленной системой произвольно размещенных прямолинейных трещиновидных полостей, поверхности которых горят. Пусть в среде, занимающей плоскость хОу, имеется N прямолинейных трещиновидных полостей длиной 2¡к (к = 1, 2,..., N). Принято, что трещиновидные полости имеют концевыми зоны со связями между берегами. В центрах трещиновидных полостей разместим начала локальных систем координат О&у к, оси которых совпадают с линиями трещин и

образуют углы ак с осью Ох (рис. 1).

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

В материалах при наличии вблизи трещины зон с нарушенной структурой в процесс разрушения вовлекается достаточно большая часть трещины. В этом случае область процесса разрушения рассматривают как некоторый слой (концевую зону), примыкающий к трещине и содержащий материал с частично нарушенными связями. Одна из возможностей моделирования такого слоя состоит в рассмотрении его как части (продолжения) трещины и в явном приложении к концевой зоне к берегам трещиновидной полости сил сцепления, сдерживающих ее раскрытие. Используется модель трещиновидной полости со связями между берегами в концевой зоне, получившая экспериментальное подтверждение [3-6] для композиционных материалов с полимерным связующим, когда адгезионная прочность меньше прочности полимеров. Рассматривается модель трещиновидной полости с силами сцепления (связями), непрерывно распределенными в узких концевых зонах полости и имеющими заданную диаграмму деформирования. Считается, что на частицы среды действуют объемные силы ^ = X + %У, где X(х,у) и У(х,у) — заданные функции. При х, у ^ го компоненты вектора перемещений и компоненты тензора напряжений в топливе стремятся к нулю.

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

При нагружении топлива в связях, соединяющих берега трещиновидных полостей, будут возникать в общем случае нормальные (хк) и касательные дХк Ук (хк) усилия (к = 1, 2, ..., N). Эти напряжения заранее неизвестны и подлежат определению в процессе решения задачи механики разрушения. В процессе нагружения топлива при некотором соотношении физических и геометрических параметров среды и действующих нагрузок будут возникать области сжимающих напряжений, в которых берега трещиновидных полостей (или их часть) могут войти в контакт, что приведет к появлению контактных напряжений на данном участке берегов трещиновидных полостей.

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

У /^

0 *

уД

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

где берега трещин вошли в контакт, будут возникать нормальные pyk (Xk) и касательные pxk yk (Xk) напряжения. Величины этих контактных напряжений заранее неизвестны и подлежит определению в процессе решения краевой задачи механики разрушения. Напомним, что в рассматриваемом случае каждая трещиновидная полость состоит из трех областей: внутренней и двух концевых зон. Внутренняя область трещиновидной полости — это противоположные берега, нагруженные давлением pk(xk) газа. Вторая область трещиновидной полости — это концевые зоны (Ai)k,di¡k) и (d2jk, A2jk) со связями между берегами. Наконец, третья область трещиновидной области — это концевые контактные зоны (-lk, A1)k) и (A2)k,lk), где берега трещиновидной полости вошли в контакт. Краевые условия в рассматриваемой задаче имеют вид

°Ук = -Pk (xk), Txk yk =0 на L' (к = 1, 21...1N),

аУк = qyk (xk), TXkyk = qxkyk (xk) на L", (1)

ayk = pyk (xk), TXk yk = pXk yk (xk) на L ■

Здесь L' — совокупность берегов трещиновидных полостей, поверхности которых горят, L" — совокупность концевых зон трещиновидных полостей, в которых действуют силы сцепления; L''' — совокупность концевых зон, в которых берега трещиновидных полостей вошли в контакт.

Основные соотношения рассматриваемой задачи должны быть дополнены уравнениями газовой динамики и уравнениями для перемещений. Уравнения течения газа в к-й полости для плоского стационарного случая имеют вид [7]

Pk = Cok РХ, ^ + = Cik, (2)

k 2 х - 1 pk

d

dx (hkpkVxk ) = Psf (Pk), hk = hok(xk)+ f (pk) dpk + Vk,

Cok = Pokp-Х, Cik = 1 2 f2(Pik) + ■

2 \PokJ X - 1 Pik

Здесь pok, pok и pik, pik — давление и плотность газа в камере сгорания при х^; =0 ив конце полости соответственно (величины с индексом 1к подлежат определению), hok(xk) — заданная толщина начальной трещиновидной полости; Vxk — скорость газа, — — показатель адиабаты; ps — плотность твердой фазы, f (p) — некоторая функция. Течение газа считается локально-изэнтропическим и безвихревым, а газ — идеальным.

Дополнительные уравнения для перемещений имеют вид

(v+ (xk, 0) - v-(xk, 0)) - i(u+ (xk, 0) - u- (xk, 0))=0 на L''', (3)

(v+ (xk, 0) - v-(xk, 0)) - i(u+ (xk, 0) - u-(xk, 0)) = C(xk,Ck)[qyk (xk) - iqxkyk (xk)] на L'', (4)

7к = V < + <Ук ■

Функции С(хк) можно рассматривать как эффективные податливости соответствующих связей, зависящие от натяжения связей, 7к — модуль вектора усилий в соответствующих связях, — V-), (п+ — и-) — соответственно нормальная и касательная составляющие раскрытия берегов к-й трещи-новидной полости.

2. МЕТОД РЕШЕНИЯ

При наличии объемных сил решение задачи представим в виде суммы

7х = 70 + 7 у, 7У = 7у + 7у, тху = тху + ТХ'У,

где 7у, 70, туу — любое частное решение уравнений плоской теории упругости при наличии объемных сил, 7у, 7Уу, туу — решение уравнений плоской теории упругости при отсутствии объемных сил.

Используя комплексные представления Колосова — Мусхелишвили [8], для компонент тензора напряжений имеем

ах + ау — 4Бе

Ф(г) -

1

дд

ау ах + 2гтху — 2

гФ'(г) + Ф(г) +

2(1 + ко) дг_ 1 д

(5)

2(1 + ко) дг

(коО - А)

в которые входят две аналитические функции Ф(г) и Ф(г) комплексной переменной г — х + ¿у и две функции д(г,г) и А1!(г,г), представляющие собой любые частные решения уравнений:

д 2д

Г — А (г, г),

д 2 А!

= А(г, г)

(6)

дгдг ч~' дг2

Здесь ко — 3 — 4г>, V — коэффициент Пуассона материала топлива.

С помощью формул (5) краевые условия задачи (1) можно записать в виде граничной задачи для отыскания двух комплексных потенциалов Ф(г) и Ф(г):

Ф(Хк ) + Ф(хк) + хк Ф'(хк) + Ф(Хк ) — / (Хк ),

(7)

где / (хк) — — Рк (хк) + /о (хк) на I/, / (хк) — ОУк (хк) — ¿ОХ ук (хк) + /о (хк) на I'

/(хк) — Рук (хк) — ¿Рхкук (хк) + /о (хк) на I"',

/о (хк ) —

-Бе

дд

1 + ко дг 2(1 + ко)

дг

дА!

дг

при у к — 0.

Комплексные потенциалы Ф(г) и Ф(г), дающие решение граничной задачи (7), ищем в виде [9]

N к

1 ^ дк (¿) ¿Ь

Ф(г) —2П £/

к=1 ,

— 1'к

Ь — гк

N ч

Ф! (г) — ¿5>-2гак / к=!

—II

дк (ь) ег

Ь — гк (Ь — гк )2

дк(Ь)

Здесь Тк — Ьегак + гк — е гак (г — ), дк(хк) — искомые функции, характеризующие раскрытия смещений при переходе через линию соответствующей трещиновидной полости с концевыми зонами

дк (хк) —

2^

д

¿(1 + ко) дхк

К (хк,0) — и_ (хк,0) + ¿^+(хк,0) — v— (хк, 0))]-

Определив по формулам (5) напряжения на оси хп (п — 1, 2,..., N) и подставив их в граничные условия (7), получим систему N сингулярных интегральных уравнений

N

^ [Кпк(Ь,х)дк(Ь) + Ьпк(Ь,х)дк(Ь)] <И — п/(х) |х| <1

(8)

к = !

Здесь Кпк(Ь,х), 1пк(Ь,х) определяются по известным соотношениям ([9], формулы (11.8)).

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

дк — 0 (к — 1, 2, ..., N)

(9)

— 1к

С помощью процедуры алгебраизации [9, 10] система комплексных сингулярных интегральных уравнений (8) при отмеченных выше условиях (9) сводится к системе N х М комплексных алгебраических уравнений для определения N х М неизвестных д°°(¿т) — vn(¿т) — гпП(Ьт) (п — 1, 2,..., N; т — 1, 2, ...,М):

М N

'У ] ^ ] 1 к |^до (Ьт )Кп к (1 к ¿то 1п хг) + д° (Ьт)1п к (1 к ¿то ¿пхг ^ — /п(хг) (г — 1 2, . . . , М — 1),

1

М

т=! к=1

. 2 т^П _ 1 ПГ

где ¿т — соб 2'2М_! п, хг — соб '

М

Х>П (Ьт)—0 (п — 1, 2,..., N),

(10)

т=!

М

1

1

к

о

к

п

к

Если в (10) перейти к комплексно сопряженным значениям, получим еще N х М алгебраических уравнений. В правые части (10) входят неизвестные значения усилий дУк (хк) и дХкУк (хк) в связях, контактные напряжения рУк (хк), рХкУк (хк) в узловых точках, принадлежащих соответствующим концевым зонам, а также неизвестные параметры А1;к и Л2,к концевых контактных зон и значения давления газа в узловых точках полостей, поверхности которых горят.

Условиями, определяющими неизвестные контактные напряжения и усилия в связях, возникающие на берегах трещиновидных полостей в соответствующих концевых зонах, являются дополнительные уравнения (3), (4).

В рассматриваемой задаче эти дополнительные условия удобнее записать для производной раскрытия перемещений берегов трещиновидной полости:

^(хк) = ¡(ТГм ^ [и+(хк'0) - (хк, 0)+

(хк, 0) - ^(х к, 0))]=0 (к = 1, 2,..., N) на (11)

(х к) = ^(1+Цко) £[С(хк,ак)(9ук(хк) - «О^ук(хк))] на Ь", (12)

где х к — аффикс точек берегов соответствующих концевых зон к-й трещиновидной полости.

Требуя выполнения условий (11) в узловых точках, содержащихся в контактных концевых зонах (—1 к, А^к) и (Л2,к,1 к), получим уравнения для определения приближенных значений контактных напряжений ) и ук (¿ш^) в узловых точках

£0 (*ш1,к ) = 0 (к = 1, 2,..., N; Ш1,к = 1, 2,...,М1,к),

где М1;к — число узловых точек, принадлежащих концевым контактным зонам к-й трещиновидной полости.

Для построения недостающих алгебраических уравнений для нахождения приближенных значений усилий дук (£ш2,к), <ЭХкук (¿ш2,к) (т2, к = 1, 2,..., М2,к) потребуем выполнения условий (12) в узловых точках, содержащихся в концевой зоне, где действуют силы сцепления (связи). При этом используем метод конечных разностей. В результате получена комплексная алгебраическая система из М2,к уравнений для определения приближенных значений дУк (¿ш2 к), дХкУк (¿ш2 к) в узловых точках концевой зоны (А1>к,^1)к) и (¿2>к, Л2,к) со связями между берегами.

Для определения размеров концевых контактных зон необходимо еще 2N уравнений. Записывая условия конечности напряжений, находим еще 2N недостающих уравнений в виде

м 2т — 1

(—1)ш£0(*ш)<*6 = 0 (к = 1, 2,..., N),

ш=1

Л 2т — 1

Е(—:1)м+ш ^о (¿ш^2^ п=0.

ш=1

Из-за неизвестных размеров концевых контактных зон система алгебраических уравнений оказалась нелинейной.

3. МЕТОДИКА ЧИСЛЕННОГО РЕШЕНИЯ И ЕГО АНАЛИЗ

Полученные системы уравнений позволяют при заданных характеристиках связей материала определить давление газа в каждой полости, усилия в связях дУк(х к) и дХкУк(х к), контактные напряжения рУк (х к), рХкУк (х к), размеры концевых контактных зон, напряженно-деформированное состояние среды.

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

Для нахождения напряжений и размеров концевых контактных зон необходимо задать закон изменения объемных сил. Функции X(х,у) и У(х,у) для упрощения разложим в ряды Тейлора в окрестности начала координат и в этом разложении ограничимся несколькими первыми членами разложения. С помощью интегрирования уравнений (6) находим

Q(z,z) = J dz J F(z,z) dz, Fi(z,z) = J dz J F(z,z) dz.

По найденным функциям Q(z,z) и Fi(z,z) согласно (6) находим функцию fo(xk).

Для решения системы уравнений использовался метод последовательных приближений [10]. В каждом приближении система уравнений решалась методом Гаусса с выбором главного элемента для различных значений M.

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

Нелинейная часть кривой деформирования связей представлялась в форме билинейной зависимости [12], восходящий участок которой соответствовал упругому деформированию связей (0 < Vk < V*) с максимальным натяжением связей. При Vk = \(u+ — u-) + i(v++ — v-)\ > V* закон деформирования описывался нелинейной зависимостью, определяемой двумя точками (V*,а*) и (5c,ac). Причем при <jc > а* имеем возрастающую линейную зависимость (линейное упрочнение, соответствующее упругопластической деформации связей).

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

\(v+ — v -) — i(u+ — u-)\ = 5c, (13)

где 5c — характеристика сопротивления материала топлива развитию в нем трещины.

Считается, что разрыв связей на краю концевой области при xk = di:k или xk = d2ik происходит при выполнении условия (13). На основании полученного решения, учитывая соотношение (4), предельное условие запишется в виде

для левого конца концевой зоны C(dijk,ak(dijk))ak(dijk) = 5c,

для правого конца концевой зоны C(d2jk,<k(d2jk))<k = 5c.

Совместное решение полученной системы уравнений и условия (13) позволяет при заданных объемных силах определить критическую интенсивность давления pk(xk) газа в каждой полости, усилия в связях, контактные напряжения, размеры контактных концевых зон для состояния предельного равновесия.

Рассмотрен приближенный расчет в предположении, что

dk

hok (x) = hko = const, f (p) = A + Bp, hok >> vk, hok >> hck, J Pk (xk) dxk =1 dk (pik + Pok )■

o

Здесь А и В — эмпирические постоянные, 2vk — взаимное смещение противоположных стенок k-й трещиновидной полости, hck — толщина сгоревшего слоя на одной стенке.

При вычислении смещений, усилий в связях, контактных напряжений и размеров концевых контактных зон считали, что на стенки k-й полости действует постоянное давление (pik + pok)/2.

На основании этих упрощений из уравнений системы (2) находим

pik pok w ,r Ps

PX = PX, pik Vxik — pok Vx0k = pik Pok hok

Adk + 2 Bdk (pik + pok)

Ps l A i P^ \ VXlk i X pik _ Vxok | X pok — (A + Bpik ), —p.--1--7- = ^--1--7-

Pok 2 X— 1 Pik 2 X— 1 Pok

Из этих соотношений определяются р1к, рхк, УХ1к, УХ0к.

Проведен параметрический анализ зависимости контактных напряжений рУк и рХкУк от размера трещиновидной полости, ее расположения и физико-механических параметров связей, относительного

z

z

z

размера концевой области. Результаты расчета контактных напряжений рудля случая одной трещины при «1 = 30°, г? = 0, 2^1 егпв/12 вдоль правой концевой контактной зоны для разных значений относительного размера концевой области е1 = (11 — й1>2)/11 представлены в виде графиков на рис. 2 (кривая 1 соответствует линейному закону деформирования связей, а кривая 2 — билинейному (нелинейному) закону деформирования связей). Здесь — силовой фактор. Сплошные линии для е1 = 0, 25, а штриховые для е1 = 0, 5. При расчетах были использованы безразмерные координаты х'к (к = 1):

0,3

-1,0

-0,5

ру1 Р

// // ____ У 0,2 А —'ч \ Ч Ч

Л // Л > // / А/' /' /* Л // / ^ / г г / > 0,1 _ \ \ ч \ ч \ ч > ч ч ч\ \\ \\ V \\ \ 4 \ А \\

// // II /' А А 1// С7/ \Л АЛ

0,5

1,0

Рис. 2. Зависимости нормальных контактных напряжений при а1 = 30° вдоль правой концевой зоны для разных значений относительного размера концевой области е1 = (11 — й1,2)/11

= ^2,к + 1к + 1к — ^2,к ' хк — + хк •

2

2

Результаты расчетов нормальных контактных напряжений для случая двух трещин а1 = 45°, г? = 0, 25^^/10 и «2 = 30°, ¿0 = 0,15^2егп0/18 вдоль правых концевых контактных зон представлены в виде графиков при относительных размерах концевых областей (11 — й1)2)/11 = 0,25 (рис. 3, а), (12 — ^2>2)/12 = 0, 25 (рис. 3, б). Здесь кривые 1 для линейного закона деформирования связей, а кривые 2 — для билинейного закона деформирования связей.

х

0

Рис. 3. Зависимости нормальных контактных напряжений для случая двух трещин а1 = 45° и а2 = 30° вдоль правой контактной зоны первой трещины (а) и второй трещины (б)

Максимальные значения контактных напряжений рУк и рХкУк находятся в средней части контактных зон, где берега трещиновидных полостей смыкаются. Подобные зависимости контактных напряжений имеют место и в левых концевых контактных зонах. При увеличении размера концевой зоны сил сцепления материала уровень усилий и qXkук в связях снижается, соответственно уменьшаются и контактные напряжения рук и рХкук. Значительное влияние на контактные напряжения оказывает место расположения трещиновидной полости. Когда полости расположены близко друг к другу, то расчеты показывают увеличение размеров концевых контактных зон. Различие видов взаимного влияния трещиновидных полостей объясняется различиями их расположения.

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

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

Библиографический список

1. Соркин, Р.Е. Газотермодинамика ракетных двигателей на твердом топливе / Р.Е. Соркин. М.: Наука, 1967. 319 с.

2. Шапиро, Я.М. Основы проектирования ракет на твердом топливе/ Я.М. Шапиро, Г.Ю. Мазинг, Н.Е. Прудников. М.: Военное изд-во, 1978. 621 с.

3. Budianscky, B. Fiber-matrix de bonding effects on cracking in aligned fiber ceramic composites / B. Budianscky, A.G. Evans, J.W. Hutchinson // Intern. J. Solid Structures. 1995. Vol. 32, № 3-4. P. 315-328.

4. Ji, H. Adhesion via Connector Molecules: The Many-stitch Problem / H. Ji, P.G. de Gennes // Macromolecules. 1993. Vol. 26. P. 520-525.

5. Cox, B.N. Concepts for bridged cracks fracture and fatigue / B.N. Cox, D.B. Marshall // Acta Met. Mat. 1994. Vol. 42, № 2. P. 341-363.

6. Goldstein, R.V. Modelling of the adhesion strength and fracture kinetics of the microelectronic package polymer-polymer joints / R.V. Goldstein, V.F. Bakirov, M.N. Perelmuter // Proc. Inst. Phys. Technol.

Russian Ac. of Sci. 1997. Vol. 13. P. 115-125.

7. Черепанов, Г.П. Механика хрупкого разрушения / Г.П. Черепанов. М.: Наука, 1974. 640 с.

8. Мусхелишвили, Н.И. Некоторые основные задачи математической теории упругости / Н.И. Мусхелишвили. М.: Наука, 1966.

9. Панасюк, В.В. Распределение напряжений около трещин в пластинках, оболочках / В.В. Панасюк, М.П. Саврук, А.П. Дацышин. Киев: Наук. думка, 1976. 443 с.

10. Мирсалимов, В.М. Неодномерные упругопласти-ческие задачи / В.М. Мирсалимов. М.: Наука, 1987. 256 с.

11. Ильюшин, А.А. Пластичность / А.А. Ильюшин. М.; Л.: Гостехиздат, 1948.

12. Гольдштейн, Р.В. Рост трещин по границе соединения материалов / Р.В. Гольдштейн, М.Н. Перельмутер // Проблемы механики: сб. статей к 90-летию со дня рождения А.Ю. Ишлинского / под ред. Д.М. Климова. М.: Физматлит, 2003. С. 221-238.

УДК 539.319

ОБ ОДНОМ ПОДХОДЕ К РЕШЕНИЮ НЕКОТОРЫХ ОДНОМЕРНЫХ КОНТАКТНЫХ ЗАДАЧ

М.А. Осипенко, Ю.И. Няшин

Пермский государственный технический университет, кафедра теоретической механики E-mail: oma@theormech.pstu.ac.ru

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

Ключевые слова: струна, балка, мембрана, пластина, односторонний контакт, одномерная задача, аналитическое решение.

ВВЕДЕНИЕ

A Certain Approach to Solving of Some One-Dimensional Contact Problems

M.A. Osipenko, Yu.I. Nyashin

Perm State Technical University, Chair of Theoretical Mechanics E-mail: oma@theormech.pstu.ac.ru

The paper deals with the problems of the unbonded contact of beams, strings, circular membranes and plates. A new approach to solving of such problems is suggested. This approach includes the rigorous problem statement, the elementary proof of the uniqueness of solution and the analytical solution construction method. The method is based on the iterative correction of the contact region. A number of examples of this method application are given.

Keywords: string, beam, membrane, plate, unbonded contact, one-dimensional problem, analytical solution.

Существует ряд одномерных контактных задач, решения которых могут быть получены в явном аналитическом виде. Это задачи об одностороннем контакте двух струн, балок, осесимметричных мембран или пластин под заданной нагрузкой (рис. 1, a, Ь, c, d соответственно) и о таком же контакте струны, балки, мембраны или пластины с абсолютно твердым телом (рис. 2, a, Ь, c, d соответственно).

Физические предположения в этих задачах примем самыми простыми: все объекты являются однородными, толщины их геометрически равны нулю, контакт происходит без трения, перемещения считаются малыми (на рис. 1 и 2 перемещения для наглядности показаны сильно увеличенными). Теории изгиба балок и пластин используем также простейшие: теории Бернулли и Кирхгофа соответственно [1].