CC BY
82
ской модели нестационарных волновых процессов в оболочках вращения при ударных торцевых воздействиях, а также предназначенные для решения краевых задач для компонент напряженно-деформированного состояния (НДС) с различными показателями изменяемости и динамичности. Приведена классификация асимптотических приближений. Изложено построение низкочастотного и коротковолнового приближений, погранслоев в окрестностях квазифронта, фронтов волн расширения и сдвига и условного фронта поверхностных волн Рэлея. Описаны схемы расчленения нестационарного НДС в случаях ударных торцевых продольных воздействий тангенциального, изгибающего, а также воздействия нормального типа.
Asymptotic Methods in Dynamics of Shells under Shock Loading L.Yu. Kossovich
The paper deals with the asymptotic methods, developed for creating a mathematic model of non-stationary wave propagation in shells of revolution under shock impacts of tangential, bending types and shock impacts of normal type; the methods are also aimed at solving the boundary value problems for the strain-stress state (SSS) components with different values of variability and dynamicity indices. Classification of asymptotic approximations is also presented. This classification defines three different types of separation scheme of non-stationary SSS. This scheme uses the following asymptotic approximations: short-wave and low-frequency ones, boundary layers in the vicinities of the quasi-front, the dilatation and shear wave fronts, and the front of Rayleigh surface waves. The schemes of ranges of applicability of approximate theories and schemes for the longitudinal stress resultant, bending moment and transverse shear force
тизация результатов разработки асимптотических методов исследования нестационарных волн для тонких упругих оболочек вращения.
последнего характерного типа НДС (возникающего при торцевых ударных воздействиях нормального типа и при нормальных поверхностных нагрузках), изученного в последние несколько лет, а также более подробным изложением разработанных асимптотических методов решения краевых задач для приближенных составляющих.
Представленные исследования являются обобщением асимптотического подхода, разработанного А.Л. Гольденвейзером [2,3] в задачах статики и стационарной динамики. Основы схемы расчленения нестационарного НДС на составляющие с различными показателями
МЕХАНИКА
УДК 539.3
АСИМПТОТИЧЕСКИЕ МЕТОДЫ В ДИНАМИКЕ ОБОЛОЧЕК ПРИ УДАРНЫХ ВОЗДЕЙСТВИЯХ
У \
Л.Ю. Коссович
Саратовский государственный университет,
кафедра математической теории упругости и биомеханики
E-mail: KossovichLYu@info.sgu.ru
В работе описаны асимптотические методы, разработанные для построения математиче-
are represented.
ВВЕДЕНИЕ
Данный обзор обобщает материалы, изложенные в [1], добавлением
Целью представленной работы является окончательная система-
© Л.Ю. Коссович, 2008
изменяемости впервые представлены в монографии [4]. Решения для нестационарных волн строятся с помощью безмоментной и моментной составляющих Кирхгофа - Лява, квазистатического погранслоя типа Сен-Венана и введенного в [4] нового типа составляющей — быстроизменяющегося по продольной координате и времени динамического погранслоя, описывающего НДС в некоторый начальный промежуток времени и в достаточно широкой окрестности фронтов волн. Эта схема расчленения нестационарного НДС явилась принципиально новой, отличающейся в корне от схем расчленения, использующихся в статике и стационарной динамике [2, 3]. В отличие от метода положения решений с различными показателями изменяемости (в частности, метода погранслоя), используется вариант метода сращиваемых разложений. При этом доказано существование областей согласования составляющих Кирхгова - Лява и динамических погранслоев различного типа, где эти решения совпадают с заданной асимптотической точностью. Наличие областей согласования показывает достаточную точность постановки краевых задач для всех составляющих, а также корректность представленной схемы расчленения. Отметим, что применение этой схемы основывается на существенной неоднородности показателей изменяемости и динамичности по координатам и времени.
Дальнейший анализ решений для динамического погранслоя показал, что динамический погранс-лой сам является достаточно неоднородным по изменяемости, вместо которого можно использовать три составляющие различного типа, которые и были введены в [5, 6]. Среди них — погранслой параболического типа в окрестности квазифронта, коротковолновые приближения для квазиплоской и квазиантиплоской задач теории упругости, а также погранслой гиперболического типа в окрестностях фронтов волн расширения и сдвига. Причем если в [5, 6] уравнения параболического погранслоя выводились на базе теории оболочек высокого порядка приближения, то в [7, 8] была разработана методика непосредственного их вывода прямо из точных трехмерных уравнений (в этих работах решались аналогичные задачи для случаев анизотропных и вязкоупругих оболочек).
Отметим, что предыдущие результаты относятся к случаям, когда нестационарные волны в оболочке возникают при ударных торцевых воздействиях, относящихся по классификации У.К. Нигула [9] к случаям продольного воздействия тангенциального типа ЬТ и продольного воздействия изгибающего типа ЬМ. В последние годы был рассмотрен третий и завершающий, с точки зрения типа схемы расчленения нестационарного НДС, вид воздействия: нормальное воздействие вида NW [9]. В отличие от предшествующих случаев, появляется новый тип НДС в окрестности условного фронта поверхностных волн Рэлея. Рассмотрение этого типа НДС началось с анализа нестационарного НДС в случае задачи Лэмба для полуплоскости [10, 11]. Выведено гиперболическое уравнение, описывающее распространение волны Рэлея вдоль границы полуплоскости, а затухание решения в глубь полуплоскости определяется путем последовательного решения двух однотипных задач Неймана для потенциалов Ламе. Полученные методы были использованы для построения погранслоя в пластинах и цилиндрических оболочках в работах [12-14]. В целом нестационарное НДС при таких воздействиях исследовано в работах [15-19].
Указанными работами полностью завершилось построение асимптотического подхода к исследованию нестационарных волн в тонких упругих оболочках. Корректность окончательной схемы расчленения определяется наличием областей согласования между различными приближениями.
но продольное воздействие тангенциального типа,
продольное воздействие изгибающего типа и нормальное воздействие, соответствующие классификации У.К. Нигула [9]. В первом случае (ЬТ-нагружение) мы имеем ненулевое значение продольного
1. ПОСТАНОВКА ЗАДАЧИ
Изучается действие ударной нагрузки, приложенной к краю и зависящей от времени как единичная функция Хевисайда.
Рассмотрим оболочку вращения (рис.1), отнесенную к криволинейным координатам (в, в, г), где 5 — длина дуги вдоль образующей, в — угол в
окружном направлении, г — координата внешней £
нормали срединной поверхности.
Рассматриваются три типа воздействия, а имен-
Рис. 1
нормального усилия. Если мы обозначим напряжения и перемещения ст, и V. = 1,2,3), то
следующие граничные условия соответствуют этой нагрузке:
стц = 1Н(£)/(0), ?2 = ?3 =0 при 5 = 0 ,
(1.1)
/ (0) = сов(т0), или / (0) = — вт(ш0).
Здесь I — амплитуда нагрузки, Н(£) — ступенчатая функция Хевисайда, £ — время. Во втором случае (ЬМ-нагружение) мы имеем на торце ненулевой изгибающий момент. Пример соответствующих граничных условий имеет вид
стц = 1гН(£)/(0), ?2 = ?3 = 0 при 5 = 0 . (1.2)
В третьем же случае (ЫШ-нагружение) мы имеем на торце ненулевую перерезывающую силу. Этой нагрузке соответствуют следующие граничные условия:
ст^ = 1Н (£)/(0), ?(1) = =0 при 5 = 0 (1.3)
Рассматриваем только однородные начальные условия
д?^ ■
V. = =0 (г = 1, 2, 3) при 1 = 0 (1.4)
и примем, что лицевые поверхности оболочки свободны от напряжений
ст 1 з = ст 23 = стзз =0 при г = ±Н, (1.5)
где Н — полутолщина оболочки. Мы также предположим, что упругие волны не достигают второго края оболочки и ограничим себя случаем, когда расстояние, пройденное фронтом, соизмеримо с характерным значением радиусов кривизны.
Приведем точные трехмерные уравнения теории упругости. Поскольку длинноволновые составляющее выводятся для самого общего случая тонких оболочек, выпишем сначала трехмерные уравнения для этого общего случая геометрии.
Пусть срединная поверхность оболочки отнесена к ортогональным криволинейным координатам (а15а2); г — как и прежде, нормальная координата. Примем для механических параметров следующие обозначения: р — плотность, Е — модуль Юнга, V — коэффициент Пуассона. Тогда уравнения движения и закона Гука запишутся в виде
1 дст.. 1 дст,. дст3. 1 дН, 1 дН.
------------1---------,--1-----------1------------- (ст.. — ст,, ) +-----(ст.,- + ст-п) +
Н. да. Н, да, дг Н.Н, да/ 331 Н.Н, да/ у ^
1 дН. Н3 1 дН. д2?
+ Н.Н, дг ст3г + Н ^7стг3 — рд2 = 0
1 дст.3 + 1 дст,3 + дст33 1 дН. 1 дН, + 1 дН. Н, +
Н. да. + Н, да, + дг Н. дг ^ Н, дг ^ + Н.Н, дг ст33+ , 1 дН, 1 дН. д2 ?3 _ .
+ Н.Н, да. стг3 + Н.Н, да, ст'3 Р д£2 ’
(1.6)
Е
2(1 + V )к2
1 — V
диз 1 дги 1 дН— 1 дН—
3 +-------------— +----------------- V- + —
Н— да— Н. Н— да. г
дг 1 дН.
Н— дг
vз +
1 дVг
+ Н дО~ + Я.#— да—
Н. дг
Е
Озз —
2(1 + V )к2
+
1 дН.
V /1 д^ 1 дv—
1 — V дО~ + Н— да— + Я.#— да— ^
1 дН—
Н Н— да.
V. +
Е
оз. —
2(1 + V)
Е
2(1 + V)
_1_ дУз + д^ Н да,- дг
1 дН.
Н. дг
1 дН
vз +
+
1 , д^" Н— дг Щ) + дг
1 дН—
1 дvi Н— да—
+
Н. дг
1 дv- 1
Н да!
дН.
Н. Н— да—
■V;. —
1 дН
Н Н— да.
где г — ^ — 1, 2 и
(1.7)
К = — — С1
1 — 2v
С1 —
Е (1 — V)
Е
С2 —
,2(1 — V)’ 1 у(1 + V)(1 — 2v)р ’ " у 2(1 + V)р •
Здесь с1 — скорость волны расширения, с2 — скорость волны сдвига, Н — параметры Ламе:
(1.8)
Н — ( 1 + "5“
Я.
(1.9)
где А. — коэффициенты первой квадратичной формы срединной поверхности, Д. — ее главные радиусы кривизны.
Перепишем уравнения (1.6)-(1.7) для оболочки вращения. Обозначим через А — А(з) расстояние от точки срединной поверхности до оси вращения, деленное на Я, где Я — характерное значение радиусов кривизны срединной поверхности. Тогда уравнения движения и закон Гука запишутся в виде
до-1
дз 1 ЯА д$ + дг
1 до 22 . ЯА д0 + Са 1 2 + доз2
дз дг
доз1 . 1 2 з Ь с д + дозз
доз1 + ( ) д2-У1 _
+ к(Оц — О 22) — Р-Е^РГ — 0,
дз
ЯА
Е
дг
+ 2ко 12 — р
д Ч д£2
д£2
— о,
+ к д2 vз — 0
+ коз1 — р“д^ — 0
О11 —
О22 —
озз —
Оз1 —
О32 —
О12 —
2(1 + V )к2 1 — V
Е V
2(1 + V )к2 1 — V
Е V
2(1 + V )к2 Е
2(1 + V)
Е
2(1 + V)
Е
2(1 + V)
дv3 1 дv2 \ дv1
а7 + ая Ж + *''7 + а7
дv3 дv1 \ 1 дv2
а7 + аг) + Зда^т + ^
дv1 1 дv2 \ дv3
1 — V V аГ + АйаТ + ь‘] + а7
дv3 дv1
дз + дг
1 дv3 дv2
АЯ~Ш + д7
1 дv1 дv2
Айа» + а7 — Ь2
где к —
1 ¿А А ¿з '
(1.10)
(1.11)
Ом —
V
—
—
V
2. КЛАССИФИКАЦИЯ АСИМПТОТИЧЕСКИХ ПРИБЛИЖЕНИИ
В теории упругих тонкостенных конструкций асимптотические методы имеют важнейшее значение как при выводе приближенных теорий, так и при решении краевых задач для составляющих НДС. Однако данная работа посвящена, в основном, общим вопросам построения НДС с помощью приближенных теорий, построенных асимптотическим методом на базе точной трехмерной теории.
Вывод асимптотически оптимальных уравнений для приближенных теорий производится с помощью метода асимптотического интегрирования трехмерных уравнений теории упругости, который основан на малости геометрического параметра n = h/R.
Динамические процессы в тонких упругих оболочках характеризуются следующими двумя физическими параметрами: l — отношение длины волны к R и T — отношение временного масштаба к Rc-1. Удобно их выразить в степенях малого параметра n
Показатели д, а являются основными понятиями. Они называются показателями изменяемости и динамичности.
В случае стационарной динамики параметры I и Т (и соответственно д и а) являются для каждой стационарной волны в области, занятой телом, как правило, постоянными. В нестационарной динамике, наоборот, они зависят от координат и времени, то есть показатели д и а имеют локальный характер.
Растянем масштаб независимых переменных, положив
и предположим, что дифференцирование по безразмерным переменным ^ и т не изменяет асимптотический порядок неизвестных величин. Зависимость неизвестных величин от толщинной переменной £ часто определяется явно при выводе приближенных уравнений асимптотическим методом. Таким образом, введение независимых переменных (2.2) с учетом изменяемости НДС позволяет методом асимптотического интегрирования трехмерных уравнений теории упругости вывести асимптотически приближенные уравнения для составляющих НДС при различных значениях показателей изменяемости и динамичности.
При построении нестационарного НДС с граничными условиями типа (1.1)-(1.3) используются в основном составляющие с показателями изменяемости и динамичности, удовлетворяющими следующим условиям:
Введем классификацию рассматриваемых асимптотических приближений. Приближения, для которых д < 1 и а < 1, отнесены к так называемым длинноволновым низкочастотным; приближения, для которых д < 1 и а — 1, отнесены к длинноволновым высокочастотным; приближения, для которых д — 1 и а < 1, отнесены к коротковолновым низкочастотным и, наконец, приближения, для которых д — 1 и а — 1, отнесены к коротковолновым высокочастотным. Эта классификация представлена на рис.2. Отметим, что длинноволновые низкочастотные приближения определяют безмоментную и изгибную составляющие, описываемые классической двумерной теорией Кирхгофа - Лява.
Для длинноволновых приближений длина волны много больше толщины тела (I >> п), в то время как для низкочастотных приближений временной масштаб много больше времени прохождения фронтами волн расширения и сдвига расстояния, равного толщине оболочки (Т >> п). Для коротковолновых приближений эти длины и для высокочастотных приближений эти времена одного порядка, то есть I ~ п и Т ~ п соответственно.
l = nq, T = n“.
(2.1)
(2.2)
0 < q < 1, 0 < а < 1.
(2.3)
Рис. 2
Разделим длинноволновые приближения на тангенциальное и поперечное (рис. 3). Для тангенциальных приближений тангенциальные компоненты вектора перемещений велики по сравнению с его нормальной компонентой (V. >> v3). Наоборот, для поперечных приближений нормальная компонента вектора перемещений велика по сравнению с его тангенциальными компонентами: v3 >> V..
Рис. 3
Коротковолновые низкочастотные приближения не представляют большого интереса в динамике. Показано [6], что они являются частным случаем (квазистатическим) коротковолновых высокочастотных приближений. В нестационарной динамике они дают хорошо известный квазистатический погранслой типа Сен-Венана, описанный для статики в работе [2]. Поэтому мы часто называем длинноволновые низкочастотные приближения низкочастотными и коротковолновые высокочастотные приближения — коротковолновыми.
Асимптотические приближения, отвечающие условию
д > 1, а > 1,
(2.4)
являются обобщением коротковолновых приближений. Они используются в нестационарной динамике при построении погранслоев в окрестностях фронтов волн. Выделены следующие типы погранслоев (рис. 4), которые названы погранслоями гиперболического типа (в окрестностях фронтов волн), параболического типа (в окрестности квазифронта двумерной волны растяжения) и эллиптического типа (в окрестности условного фронта поверхностных волн Релея).
Механика
17
3. НИЗКОЧАСТОТНЫЕ ПРИБЛИЖЕНИЯ
Рассмотрим длинноволновые низкочастотные приближения трехмерных уравнений теории упругости в общем случае тонких оболочек произвольной формы.
В случае плоской пластины симметричный и антисимметричный по нормальной координате типы граничных условий порождают симметричный и антисимметричный типы НДС. При этом длинноволновые низкочастотные приближения разделяются на два типа: тангенциальные и поперечные, соответствующие теориям растяжения и изгиба тонких пластин. Естественно предположить, что также имеют место два разных рассматриваемых приближения в общем случае тонких оболочек. Начнем рассмотрение с тангенциального приближения.
При построении тангенциального приближения показатели изменяемости и динамичности связаны соотношением q = a (0 < q < 1) и компоненты НДС имеют следующую асимптотику:
V = R(nq v° + n1+q ), v3 = R(nv0 + n2q ),
0ii = E (a“, + 4CT1(), CTij = -E(CT° + n01), (3.1)
03i = E(n3-3q<7° + n2-q), 033 = E(n2-2q0°3 + n3-2q013), (i = j = 1,2).
Здесь предполагается, что величины с индексами имеют один и тот же асимптотический порядок.
Формулы (3.1) обеспечивают разделение искомого НДС на основную и дополнительную компоненты. Величины, входящие в основное НДС, отмечены индексом «0», в то время как величины, входящие в дополнительное НДС, отмечены индексом «1».
Основное НДС аналогично НДС плоского слоя, тогда как дополнительное — результат действия поверхностной кривизны. Если основное НДС симметрично относительно срединной поверхности, то дополнительное антисимметрично, и наоборот. Для симметричного НДС v, иц, 0j, 033 являются четными функциями толщинной координаты, и v3, 03i — нечетные функции. Для антисимметричного НДС все величины имеют противоположную четность, а именно V, 0 ц, 0j, о33 являются нечетными функциями толщинной координаты, а v3, 03i — четными.
В рассматриваемом случае основное НДС симметрично относительно срединной поверхности. С точки зрения использования рассматриваемых приближений в нестационарных задачах произведенный набор требований, накладываемых на показатели изменяемости, динамичности и интенсивности, позволяет удовлетворить этим приближениям (по обычной терминологии безмоментной составляющей) тангенциальной части граничных условий в двумерной форме, а также описать распространение волн сжатия и сдвига (асимптотически главная часть разрешающей системы для приближения содержит старшие производные от тангенциальных перемещений по продольной координате и времени).
Описание процесса асимптотического интегрирования трехмерных уравнений теории упругости приведено подробно в [10, 6]. Здесь мы укажем только основные его вехи. Подставим асимптотики
(3.1) в уравнения (1.5)-(1.6) и проведем замену переменных согласно (2.2) при a = q. Разделяя каждое из уравнений на две части, принимая во внимание четность членов по Z и пренебрегая членами порядка O(n2-2q) (асимптотическая погрешность классической теории оболочек Кирхгова -Лява), приходим к замкнутой системе относительно асимптотически главных компонент V3, v°, 0°, и 0°j. Выведена также система, определяющая асимптотически второстепенные компоненты через асимптотически главные. Эти две системы дают возможность установить следующую зависимость от нормальной координаты:
V° = v(0) > v3 = v30), 0°i = 0LO) > 0°j = 0ljO) >
v1 = Cv(1), v° = с43), 4 = C0i1) > 0i = C01, (3.2)
01 = 0(°) + z 2 0(2) 0O = 0(°) + z 2 0(2)
03i = 03i + z 03i , 033 = 033 + z 033 •
Выведены уравнения для величин с индексами в скобках, не зависящих от Z. Запишем разрешающую систему относительно следующих общепринятых в теории оболочек двумерных величин: Ti —
продольные усилия, Sj — сдвигающие усилия, u, — тангенциальные перемещения точек срединной поверхности, w — прогиб:
Т =
2ЕЛ, Ґ 1 ди* , т / 1 ди, , т
Т-^1 А да + ^ + Ж + Ча дО“ + к'г + д
Ж Г А* д
<?• • =
1 + V I А, да.
и*
А*
А, д / и,
А* да* V А
2_ дЛ
А* да*
+
1 д^*“ А“ да“
+ к, (Т* — Т“) + 2к*^*“ — 2рЛ,
д2и*
дї2
= 0,
где к* =
1 дА*
А* А, да,
Ті + Т2 +2 ,д2ю =0
Ж + Д2 +2р,і зї2 =0 (і = І = 1,2).
(3.3)
Рассмотрим теперь поперечное приближение. Показатели изменяемости и динамичности связаны соотношением q = (а + 1)/2 (0 < а < 1), а асимптотика НДС определяется следующим образом:
V = Я(пи° + п2д V,1), = Я(п9 ^0 + п9+1 VI),
Стгг = Е (П1 9 /°* + /г), ^ = Е (/ 9 /°/ + ),
(3.4)
аэг = Е(П2 2903г + ПЭ 2903г), 0ЭЭ = Е(/ ^0ЗЭ + П2 90ЗЭ) (* = 3 = 1 2)
Основное НДС антисимметрично относительно срединной поверхности. С точки зрения использования рассматриваемого приближения в нестационарных задачах, указанный выбор показателей изменяемости, динамичности и интенсивности позволяет поперечному приближению (изгибной составляющей по обычной терминологии) удовлетворить нетангенциальной части граничных условий, снять невязку, появившуюся в этих условиях после построения безмоментной составляющей, а также обеспечить частичное моделирование поперечной изгибной волны (показатели изменяемости и динамичности выбираются так, чтобы асимптотически главная часть разрешающей системы содержала старшие производные от прогиба по продольной координате и времени). Выведены системы для асимптотически главных и второстепенных компонент НДС [6, 10], дающие возможность установить следующую зависимость от нормальной координаты:
і _ я,(°) „,° _ я,(°)
Юо = Vg
1 (°)
1 (°) ^ = /У ,
/°* = с/*0,
= С/У,
(3.5)
г° _ Л°)
г3*
'3* + С °3* ’ ^33 _ ^33 "Т" Ч и33 І ^33 _ 4^33 "Т" Ч ^33 •
Выведены уравнения для величин с индексами в скобках, не зависящих от £. Запишем разрешающую систему относительно усилий, перемещений точек срединной поверхности, а также изгибающих моментов С*, скручивающих моментов Н, и перерезывающих сил N:
2 _(2)
Т1 -~(°) 33
2 (2)
° (1)
3 (3)
с* = —
2ЕЛ3
1д
1 дт \ к 1 дт 3(1 — V2) [ А* да* ^ А* да* у + * А, да, + ^
2Е^3
н = —
* 3(1 + V) \ А, да, ^А* да*
1д
1 дт
к 1 дт
— , а да;
т* =
2Ж
1 ди* , т / 1 ди, , т
1 — V 2 1 А* да + к*и + д* + Ча д0” + ^ + д
_ = ГА* д ; 1 — ^А, да,
+
А,- д
А* да* А
и,
Л дТ
А* да* +
+ к, (Т* — Т) + 2к* S*J — 0,
Т± + ^2 — ± дЖ1
^1 ^2 А1 да1
1_ т1
а, да,
А* да*
+
1 дЖ2 д2т
I- 75------к2Ж1 — к1 Ж2 + 2Р^ о,2 = 0?
А2 да 2 дї2
+ к, (С* — С,) + 2к*Н*, — Ж* — 0 .
1 дт А, да,
+ к
1 дт 3 А* да*
(3.6)
Отметим, что полученные приближения являются, по сути, асимптотическим обобщением соответствующих составляющих теории Кирхгофа - Лява. Одним из принципиальных положений этого обобщения является возможность (а во многих случаях и необходимость) учета изменения длины
нормального элемента:
4. КОРОТКОВОЛНОВЫЕ ПРИБЛИЖЕНИЯ
£з —
2ЕЛ,
(Т + Т2 )
(3.7)
В случае коротковолновых высокочастотных приближений (д = а = 1) легко выписать асимптотику НДС. Дело в том, что при д = 1 длина волны имеет порядок толщины, и все направления становятся эквивалентными. Очевидно, что искомое НДС симметрично или антисимметрично относительно срединной поверхности с асимптотической погрешностью О(п). Следовательно, асимптотика НДС имеет вид
аы — Е(акг + паУ , уш — Дп (УШ + пуШ) (к, 1, т — 1, 2, 3)
(4.1)
Здесь, как в предыдущем пункте, мы разделяем искомое НДС на основную (величины с индексом «0») и дополнительную компоненты.
Выпишем уравнения в размерной форме для асимптотически главных компонент НДС:
1 даіі 1 да, , даз
А* да
+
А,- да,
+
1 дст*з 1 да,з дазз , . , д2^з
+—;—~-------1—~-+ к, а*з + к*а,з — Р 0^2 — О,
А* да і
А, да,
дг
Е
азз —
2(1 + V )к2 Е
2(1 + V )к2
V / дуз 1 ду
і I я + Т” + к, ) + л я ,Ъг ^
1 — V \ дг А, да, ) А* да*
д£2 1 ду*
к* у,
(4.2)
1 ду.
1 — V V А* да*
1 ду + А, да,
+
дуз
дг
Е
аз* —
2(1 + V)
_! дуз + дУ*
А* да* дг
Е
2(1 + V)
1 ду* 1 ду,
I- я--------+ ------к*у* — к,
А, да, А* да*
(і — І — 1, 2)
Уравнения (4.2) аналогичны уравнениям трехмерной упругости при условии, что метрика срединной плоскости слоя совпадает с метрикой срединной поверхности оболочки. Коэффициенты этих уравнений не зависят от толщинной координаты а3. Более того, они являются функциями коэффициентов первой квадратичной формы срединной поверхности.
В случае оболочки вращения разрешающие уравнения (4.2) могут быть значительно упрощены. Рассмотрим в деталях случай, когда колебания являются коротковолновыми только в направлении меридиана. Предположим, что изменяемость по окружной координате имеет показатель изменяемости р < 1. Тогда выделяются два типа асимптотик НДС. Асимптотика первого типа имеет вид
ап — Еа*, сті2 — Еп
— Еп1-Ра
і2,
у1 — Лпу*, у2 — Дп
— Лп2-рУ*
азі — Еа*і, аз2 — Еп уз — Лпу* (I — 1, 2, 3).
— Еп1-ра*2,
(4.3)
При этом асимптотически главными компонентами НДС являются ац (I = 1, 2, 3), 031, VI и у3. Заметим, что эти напряжения и перемещения определяют плоскую задачу теории упругости, а асимптотически главные части разрешающих уравнений совпадают с уравнениями плоской задачи теории упругости. Асимптотически второстепенные напряжения и перемещения появляются в этих уравнениях только с множителем п2-2р.
Анализ показал [4, 6], что для построения решения для колебаний с погрешностью
5 = О (п + п2-2р) разрешающие уравнения коротковолнового приближения должны быть построены
V
а—
*
с асимптотической точностью <5, то есть члены порядка О (п + п2-2р) должны быть в разрешающих уравнениях сохранены.
Исследования также показали [4, 6], что второе разрешающее уравнение движения может быть максимально упрощено (оно заменяется уравнением первого порядка), и системы разрешающих уравнений в размерном виде записываются следующим образом:
да
11
+
1 да
12
да
+ ~7Г1 + к(ац - а22) - р
дз ДА д6>
1 ду1 ду2
ДА~Ш - "дЗ"
дг
д2 у1 д*2
= 0,
= 0,
да
31
дз а11 =
а22 =
а33 =
а31 =
а32 =
а12 =
+
1
ДА д6> Е
да32 да33
+ —-------+ каз1
дг
V
2(1 + V )к2 1 — V
Е V
2(1 + V )к2 1 — V
Е V
+
ду1
дз
ду3 ду1
дг + дз
+
1 дг>2 да д#
2(1 + V )к2 Е
2(1 + V)
Е
дг>1
+
1 — V V дз ДА
ду3 ду1 дз + дг 1 ду3 ду2 АД Ж + д/
+ кгч +
+ кг>1 дУз
дг
2(1 + V)
Е /1 ду1 ду2 2(1 + V) ^АД д^ + д7
Отметим, что модифицированное второе уравнение движения имеет очевидный физический смысл. Оно соответствует отсутствию вращения вокруг нормали к срединной поверхности. Заметим также, что с введением этого модифицированного уравнения система описывает только волны расширения и поперечные волны, определяющие деформацию изгиба.
Перейдем к рассмотрению асимптотики второго типа:
аи = Еп1 р аш У1 = Дп2-^*,
а12 = Еа*2, аз1 = Еп1 ра*ь аз2 = Еа*2, у2 = Дпу*, уз = Дп2-РV* (1 = 1, 2, 3).
(4.4)
При этом величины, характеризующие антиплоскую задачу теории упругости (напряжения а12, аз2 и перемещение у2), являются асимптотически главными, а асимптотически главные части разрешающих уравнений совпадают с уравнениями антиплоской задачи теории упругости.
Анализ показал, что в рассматриваемом случае можно модифицировать первое и третье уравнения движения, и система разрешающих уравнений в размерном виде записывается следующим образом:
ду1 1 ду2
“дз + Да "д^ = ’
Уз = 0,
1 да22 да12 даз2 д2у2
22 + + 2ка12 — р—/ = 0,
ДА
а11 = азз = а12 =
дз
Е
дг
V
д*2
1
ду2 ду1
+ дз
2(1 + V)к2 \1 — V ДА д$
Ev / ду1 1 ду2
2(1 + V2)к2 V"дз + ДА д^
Е / 1 ду1 ду2 к 2(1 + V) ^ДАд^ + д7 — У2
Е
а22 =
V д^1
2(1 + V)к2 \1 — V дз
+
аз1 =
Е дУ1
2(1 + V) д/’
Е
аз2 =
1 дг>2 ДА~Ш
дУ2
(4.5)
2(1 + V) дг ’
Физический смысл модифицированных уравнений движения ясен: они соответствуют отсутствию объемного расширения и нормального перемещения.
5. ПОГРАНСЛОЙ В ОКРЕСТНОСТИ КВАЗИФРОНТА
Хорошо известно, что как классическая теория обобщенного плоского напряженного состояния, так и теория оболочек Кирхгофа - Лява искажают скорость волны расширения. Вместо с1 они задают скорость сз = у/Е/р(1 — V2). Фронт волны расширения по классическим двумерным теориям пластин и оболочек является фактически квазифронтом: в его окрестности трехмерное решение хотя и быстро изменяется, но является непрерывным.
Предположим, что в окрестности квазифронта преобладает длинноволновый характер решения и показатель изменяемости принадлежит полуинтервалу
3 < ?< 1. (5.1)
Классическая теория оболочек не работает на этом полуинтервале, и при построении асимпто-
тически оптимальных уравнений следует учитывать члены более высокого порядка малости. Введем безразмерные переменные х и *з [7, 8]
1
х = (*Сз — з), *з = Д (5.2)
и предположим, что дифференцирование по х и *з не изменяет асимптотического порядка неизвестных
величин. Считаем, что д = а = 2/3 и толщина погранслоя имеет порядок О (п2/з).
Пусть асимптотика НДС в этой области определяется выражениями:
,Л ,,2 = Д (п47з-р„° + пг/з-Ру2
«1 = Д (п2/зV? + п5/з«!), ^2 = Д (п4/з-р«0 + п7/з-М),
Уз = Д (п«з + п4/зV1) , а11 = Е(а?1 + па^), а22 = Е (а°2 + па^) ,
азз = Е (п2/за°з + п5/за1з) , а1з = Е ^Псг]°з + п4/заЬ) ,
а2з = Е (п5/з-ра?з + П2-Мз) , а12 = Е (п2/з-ра°2 + п5^-^) .
(5.3)
Пренебрегая асимптотически малыми членами в разрешающих уравнениях (1.9)-(1.10) с учетом
(5.2)-(5.3) и интегрируя эти уравнения по £, устанавливаем зависимость неизвестных величин от нормальной координаты:
2„(2)\ „,°_„,(°) „Д_„,(°)
У? = +п2/з +с2 <42;), У? = 4з),
Уз = гч
а?1 = а(1,з+п2/з (а1°1)1+с, а?2 = а(2), а° = 4°, (5 4)
/ = СУг(1), У° = С^з1), 4 , ^ ,
а1 = а(з) + А 2 а(2) а° = а(0) + Л 2 а(2)
°3г = °3г + ^ °3г , азз = азз + Ч азз .
Здесь учитывается тот факт, что для построения уравнений погранслоя в окрестности квазифронта компоненты НДС и° и а°1 должны определяться с большей, чем остальные компоненты, точностью
О (п2/з).
Введем усилия и перемещения срединной поверхности по формулам:
Т1 = 2Е^а(°)°, Т2 = 2Е^а(2), = 2Е^2/з-р а/,
’ ^ (5 5)
и1 = Дп2/зУ1°°, и2 = Дп4/з-ру2°) , 'ш = Дп4/зУз°) .
Тогда для асимптотически главных двумерных компонент НДС и1, 71, Т2 получаем следующие уравнения:
д4и1 к дзи1 V2 д6и1
+ Д о о з— ТЛ--72 п 6 = 0,
дхзд*з 2 дхз 6(1 — V)2 дх6 /г
(5.6)
2ЕЛ, 1 ди1 2Е^ 1 ди1
7 =-------------------------1 Т2 =---------------------------1.
1 — V2 Дп2/з дх , 1 — V2 Дп2/з дх
Рассмотрим погранслой для волн, распространяющихся в направлении возрастания координаты 5. В этом случае можно получить разрешающее уравнение относительно продольного перемещения, допускающее простые решения, выражающиеся через функцию Эйри. Получим искомое уравнение, проинтегрировав (5.6) по х:
ди і к
~дЦ + Я 2 “ 1
д 3 и 1
6(1 — V)2 дх3
= 0.
(5.7)
Перейдем к исходным координатам 5, 0, £. Уравнения для асимптотически главных компонент примут вид:
ди
21,2
V 2 Л
д3и1
1 ди1 к
дз + с3 д£ +2 и1 6(1 — V)2еЦ д£3
= 0,
2ЕЛ ди
1 — V2 д5
Т2 =
а для асимптотически второстепенных компонент —
ди?
ди2 1 ди1
+ дАЖ = ’
д5
$12 =
д5
+
2ЕЛ^ ди1
1 — V2 д5
1 V
Я + ^ Іиі = 0
Я1 Я2
(5.8)
ЕЛ 1 + V
1
ди1 ди2
ЯАЖ + “д7
(5.9)
2
V
6. ПОГРАНСЛОЙ В ОКРЕСТНОСТИ ФРОНТОВ ВОЛН РАСШИРЕНИЯ И ВОЛНЫ СДВИГА
В этой части мы описываем вывод уравнений погранслоя в окрестности фронта на примере волны расширения. Рассмотрим сначала погранслой в случае продольных воздействий тангенциального и изгибающего типов вида (1.1) и (1.2), когда основной разрыв решения, порожденный ударной нагрузкой, переносится нормальным напряжением а11. Уравнения выводятся для оболочек вращения нулевой гауссовой кривизны (цилиндрические и конические оболочки, а также пластины), когда фронт волны образуется совокупностью нормалей к срединной поверхности.
Введем безразмерные переменные
1
х =
^е1
~Я
(6.1)
и предположим, что дифференцирование по х, ¿4 не изменяет порядка неизвестных величин. Считаем, что толщина погранслоя имеет порядок О (п2) и д = а = 2. Определим асимптотику НДС формулами
у1 = , ^2 = Яг/
= Рп3-Р-Л*
У3 = ДпМ,
стіі = Еп 1а1і а13 = Еа-
О 22 = Еп
-1 і О22,
о 33 = Еп 1_ 1
(6.2)
г*3,
О23 = Еп2 Р О?3,
о 12 = Еп1 рст*~
12,
где все величины со звездочкой имеют одинаковый асимптотический порядок.
Переходя к переменным (6.1) и к величинам со звездочкой и удерживая асимптотически главные члены, получим из уравнений (1. 10)-(1.11) замкнутую подсистему для асимптотически главных компонент «2, а^1, а22, а|з:
д2
дС2
'и
-2
д2V* А' д-и*
дхд^1 1
— т- =0,
А дх
д-и*
2(1 + V)к2 дх ’
*
О22 =
д-и*
(1 — 2v)(1 + V) дх ’
*
О33 =
д-и*
(1 — 2v)(1 + V) дх
Вторая подсистема служит для определения асимптотически второстепенных компонент а2з, а22 через асимптотически главную компоненту г 2:
д^*
1 д-и *
дх А
= 0,
д-и* д-и*
—1 +--------------3 = 0,
дС + дх ’
*
01 3 =
1
*
О23 =
1
2(1 + V) V дС А
*
01 2 =
1
1 + V А
(6.3)
у2 , ^3 , О 1 3,
(6.4)
*
і
V
V
*
ди*
Здесь, как и в предшествующем случае, мы имеем право пренебрегать в уравнении асимптотически второстепенными членами, оставляя в них только асимптотически главную часть. Данный вывод следует из того, что безразмерные переменные (6.1) имеют порядок единицы в рассматриваемой окрестности фронта.
В исходных координатах рассматриваемые подсистемы запишутся в виде:
d2 v1 д2г>і
ds2 + dz2 ci dt2
1 d2 v1 A' dv1
1 + — = 0,
E
dv1
Ev
an =
2(1 + v)k2 ds
^22 =
A ds
dv1
Ev
(1 — 2v)(1 + v) ds ’
азз =
dv1
(6.5)
(1 — 2v)(1 + v) ds
dv2 . 1 dv1 =o dV1 — dV3 =o
ds + AR d0 = ’ dz ds = ’
a13 =
E dv1
E
a23 =
dv2 1 dv3
~d7 + AR"d^
a 12 =
E 1 dv1
(6.6)
1 + v dz 2(1 + v) ^ dz ' AR d0 ) 1 " ^ 1 + v AR
Можно проверить, что решения должны удовлетворять только одному граничному условию на лицевых поверхностях
а33 = 0 при z = ±h (6.7)
и одному граничному условию на торце s = 0:
а11 = IH (t)f (0) (LT — нагружение),
а11 = IzH (t)s(0) (LM — нагружение).
(6.8)
Отметим, что если исходные трехмерные уравнения теории упругости описывают три типа волн (продольные волны — волны расширения и поперечные волны — волны сдвига), то полученные уравнения погранслоя в окрестности фронта волны описывают только волны одного типа, задающие рассматриваемый фронт.
При воздействии вида ЫШ передний фронт волны сдвига несет главный разрыв напряжения ст13, отражающий скачок напряжения в начальный момент времени. В рассматриваемом случае оболочек вращения нулевой гауссовой кривизны фронт образован нормалями к срединной поверхности и определяется уравнением 5 = е2£.
Асимптотику НДС погранслоя определим следующим образом:
v1 = Rn2 v3, v2 = Rn
= Rn4-pvi
V3 = Rnv3,
a11 = Ea
<T13 = En a
11,
-1 3
13,
a22 = Ea
22,
a33 = Ea
33,
(6.9)
a23 = En1 p a^,
a12 = En2 pa3-
12
В рассматриваемом случае асимптотически главной компонентой вектора перемещений является ■и3, а главным напряжением — а13. Тогда аналогично предыдущему случаю выводим систему для асимптотически главных компонент НДС
d2V3 d2V3 A' dv3 1 d2V3
+-----^"o--1------ -----------2 ^.O -- 0,
dz2 ds2 A ds c2 dt2
и для асимптотически второстепенных компонент
a13 =
E dv3 2(1 + v) ds
(6.10)
_d_
ds
d2 v2 d2 v2 1 d2v2
dz2 + ds2 c2 dt2
A' dv2 A ds
dv1 + dV3 = 0 ds + dz = ,
1 d / d2v1 d2v1
+ ARd^ v’dz2 + "ds2
1
d2 v1 \
A' d2v1 + ^^7^ =0,
c2 dt2 J A2Rdsd0
a11 =
E dv3 1 + v dz ,
2v E
dv3
E
a22 =
1 dv3
(1 — 2v)(1 + v) dz ’
a33 =
E dv3 1 + v dz
(6.11)
E
2(1 + v) AR
dv2 + _1 2(1 + v) V ds + AR
dv1
3
3
3
7. ПОГРАНСЛОИ В ОКРЕСТНОСТИ УСЛОВНОГО ФРОНТА ПОВЕРХНОСТНЫХ ВОЛН РЕЛЕЯ
Построим погранслой в окрестности условного фронта поверхностных волн Релея для цилиндрической оболочки [13, 19].
Анализ разрешающих уравнений показывает, что искомое НДС, строящееся с показателями изменяемости и динамичности, равными единице, относится к классу коротковолновых и является симметричным или антисимметричным относительно срединной поверхности с асимптотической погрешностью О(п):
сты = Е(ст° + п^У , Ут = Лп/У + П^У) , м,т = 1, 2, 3,
(7.1)
где индексом «0» обозначена основная составляющая НДС (симметричная или антисимметричная относительно срединной поверхности), а индексом «1» обозначена второстепенная составляющая (соответственно антисимметричная или симметричная). Тогда разрешающие уравнения для основной составляющей будут переписаны в виде, совпадающем с уравнениями для плоского слоя:
-2 д2г>1 д2г>1 1 д2 у3 д2
к 2^гг + ^ттт + -------------------— тттгтт----------~ ~ = 0,
дС2 д£2 1 — 2^ дСд£ дт2
1 д2и1 , д2у3 , -2 д2у3 д2у3
(7.2)
+
+ к-
1 — 2^д£д£ дС2
Е (дУ1
СТ11 =
ст22 =
стзз =
+
д£2 дт2
V ду3
= 0,
2(1 + V)к2Л, \ дС 1 — V д£ у ’ Е
V ду1 + V ду3
2(1 + V)к2Л, \1 — vдС 1 — V д£ у ’ Е
V ду1 + ду3
(7.3)
СТ13 =
2(1 + V)к2Л, \1 — V дС д^ У Е (ду1 ду3 2(1+ V)Л. V9( + 9С
где С = 5/^, С = ¿/Л., т = ¿С2/с1.
Аналогично случаю выведения уравнений погранслоя для полуполосы [14] выделим частное решение уравнений (7.2)-(7.3), удовлетворяющее только граничным условиям (1.3) при /(0) = 1. Обозначим это решение верхним индексом «0»:
,(°)
= 0,
,(°)
= — I (т — С )Н (т — С),
(0) _ „(») _ „(0) _
СТ11 =
22
= ст33 = о,
г(0) = 13
Е1
2(1 + V)
Н (т — С).
(7.4)
Решение (7.4), построенное с малой асимптотической погрешностью О(п), представляет собой поперечную волну, распространяющуюся вдоль оси цилиндрической оболочки в плоскостях, проходящих через эту ось. Построенное таким образом решение не удовлетворяет граничным условиям (1.5) на лицевых поверхностях.
Представим решение рассматриваемой задачи в виде следующей суммы:
НДС = НДС(0) + НДС(1),
(7.5)
где для решения с верхним индексом «1» на торце и лицевых поверхностях ставятся следующие граничные условия:
С = 0:
ду
(1)
дС
= 0,
с = ±1 : ст(3) = о,
(1)
= 0,
ст(1) = — °13 =
Е1
(7.6)
2(1 + V)
Н (т — С).
1
3
1
Граничные условия (7.5) показывают, что деформированное торцевое сечение остается прямолинейным и перпендикулярным к оси цилиндра. Следовательно, если взять бесконечную цилиндрическую оболочку с симметричным нагружением, показанным на рис.6, то деформирование каждой из ее симметричных частей эквивалентно деформированию рассматриваемой полубесконечной оболочки. Поэтому можно перейти к исследованию эквивалентной задачи для бесконечной оболочки, представленной на рис. 6.
Рис. 6
Тогда вывод уравнений погранслоя можно провести аналогично случаю бесконечной цилиндрической оболочки [13, 14] с помощью символического метода А.И. Лурье. Запишем систему уравнений
(7.2)-(7.3) в операторном виде
¿2у1 , { -2^2 ^2 А , д5 ¿У3 п
и? +(к д«- ^У1 + т-2їїіс ’
¿у1 -2 ¿2у3 . 2 Я2\ п
1-2?Ц + к 1с2 + (д - 9т) Ю3 = 0’
СТ11 =
733 =
<713 =
Е / V ¿у3 2(1 + V)к2Н у1 — V Е
+ У1 ,
где операторы , дт имеют вид Граничные условия имеют вид
С = 0:
2(1 + V)к2Н \ ¿С 1 _ V
2(ТТ?)Н (д{У3 + Ю '
я д я д
д{ = д?’ дт = д?
¿У3
¿С
= 0 , у1 = 0 ,
£ = ±1 : 733 =0 , 713 = —
Е/
2(1 + V )Н
Н (т — ?).
(7.7)
(7.8)
(7.9)
(7.10)
Решаем систему (7.7), рассматриваемую как систему обыкновенных дифференциальных уравнений, и удовлетворяем граничным условиям (7.10). Пренебрегая членами малого порядка [13,14] и рассматривая большие значения времени т >> 1 (когда фронты волн проходят расстояние, много
большее толщины полосы): т = 0(Т), Т >> 1, введем малый параметр е = 1/Т и масштабированные
переменные у, то:
Со — кдто * *
У =---------------, Со = то = ет,
(7.11)
где коэффициент кД = СД/е2 < 1, ед — скорость поверхностных волн Релея. Коэффициент кд является корнем известного трансцендентного уравнения:
(2 - кД )2 - 4\Л - КД\/1 - к2кД = 0*
(7.12)
Введем новые операторы ду, дт
Я Я - А
дУ п , дТ0 о *
ду 0 дто
(7.13)
Отбрасывая в операторном решении асимптотически второстепенные члены и возвращаясь к исходным операторам, введем объемный и сдвиговой потенциалы Ламе ^ и 0. Разделим их на части, обозначая индексом «1» отрицательное по С направление распространения возмущения, а индексом «2» — положительное:
У = уі + У2, Ф = Фі + 02*
(7.14)
Эти потенциальные функции определяются эллиптическими уравнениями:
1 - СД 1 ,2
1-----------Д
1 ,2
д2Уі + д 2 у і = 0
с2 у дз2 дг2
еД У д2фі д2фі
+
с2 У дз2 дг2
=О
(7.15)
при следующих граничных условиях:
д2ф1 2 д2ф1 , 2 Вд2ф2
— СД о 2---+ 2кДВ‘
Кд в
д£2
дз2
2 ид 2 ф1 + д2 ф2 2кдВ^^- +
дз2
д£2
дз2 2 д2 ф1 Кд^
1 — КД/2
Кд В
1 — КД/2
$, г = Л,
$, г = —Л,
(7.16)
дуі дг ду2 _ дг
1 — КД/2
дфі дф2
1 — кД/2 V дз
дз дз дфі дф2
дз
г = Л,
г = —Л,
(7.17)
где постоянная В и функция $ определяются выражениями:
В = 2
кд
22
,
к к
4кД
1 — кД 1 — к2кД 2 — кД
-і
$ =
2
2/я(т — С), С > о,
—1 /я(т — С), С < 0.
Перемещения и напряжения выражаются через потенциальные функции по формулам
= дуі , , дФі + дф2
дз дз дг дг ’
дуі ду2
!’1 = 77 + 77
дфі
дз
дф2 дз ’
Е
733 =
1 + V
К
д2 ф2 д2ф2
дздг дздг
Е
1 + V
— | 1 — :ія) ^
V 2 у у <9з2 ду2
д2уі + д2у2 — Л — кДу д2Фі + д2ф2
дздг дздг V 2 / V дз2 дз2
(7.18)
1
1
д
8. ВОЛНЫ ПРИ иГ НАГРУЖЕНИИ
Проиллюстрируем расчленение НДС на составляющие с различными показателями изменяемости при продольных воздействиях тангенциального типа схемой решения для продольного усилия в цилиндрической оболочке (осесимметричный случай) в некоторый момент времени, который больше времени пробега передним фронтом волны толщины оболочки (рис.7) и схемой областей применимости приближенных теорий (рис.8).
Рис. 7 Рис. 8
Применяются [4, 6] следующие приближенные теории: К1 — гиперболический погранслой в окрестности фронта волны расширения, К2 — квазиплоская квазисимметричная задача теории упругости, К3 — параболический погранслой в окрестности квазифронта, К4 — тангенциальное длинноволновое приближение (безмоментная составляющая теории Кирхгофа - Лява), К5 — наложение двумерной составляющей и квазистатического погранслоя типа Сен-Венана; В1, В2, В3 — области согласования.
Исследования работы [4] с учетом выделения вместо динамического погранслоя двух типов по-гранслоев и квазиплоской задачи, выполенного в [6], показали, что область согласования между без-моментной составляющей и погранслоем в окрестности квазифронта определяется асимптотическим неравенством
П1-р + П4/3-2р + П2/3 << 1 (сзЬ - в) << 1. (8.1)
К
В работе [20] найдено расположение области согласования между коротковолновой составляющей (квазиплоская задача теории упругости) и погранслоем в окрестности фронта волны расширения:
П2 << "1 (с1 Ь - в) <<п- (8.2)
К
Отметим, что проведенные исследования сняли дискутирующийся долгое время вопрос о возможности использования параболических уравнений теории Кирхгофа - Лява для описания волновых процессов в оболочках, указав с помощью асимптотических оценок границы области применимости безмоментных и изгибных составляющих. Результаты исследования полностью согласуются с выводами В.В. Новожилова и Л.И. Слепяна [21, 22] о действии принципа Сен-Венана в динамике пластин
и оболочек.
В соответствии с (1.1) граничные условия для безмоментной составляющей записываются в виде
Т1 = 2^/Я (*)/(0), и2 = 0. (8.3)
Рассмотрим случай быстроизменяющей по окружной координате неосесимметричной нагрузки, когда р > 0. Для решения краевой задачи для безмоментной составляющей в [4]
использовано интегральное преобразование Лапласа по времени с последующим обращением изображения различными асимптотическими методами в зависимости от диапазона изменения продольной пространственной переменной и времени. Так, для волны растяжения при
Со(то — Со) = О (п2р) (Со = в/К, то = с3Ь/К )метод прифронтовой асимптотики дал решение
де разложения по функциям Бесселя
Т1р) = 2Л/
¿о ^\j2ml п-2рСот^ + 2(1 со^)Тр ¿2 ^\j2rn2 п-2рСо тр^ + ...
/ (*)# (тр),
(8.4)
тр = то — Со, т* = тпр ,
а при Со(то — Со) >> П2р метод перевала позволяет представить решение для волны растяжения в виде
Т1(р) = 2М
1 — (1 — ^)т* Со + (1 — V)т2 с2
+
пр/2Со
п1/2 т*/2то Тр/4
1-
2 пр
(1 — V )тр
2р
4/2
I Со .
ехр —т* — +
т*Тр . т* Тр
сое-------------+ эш---------
1/2
1+ О
Тр
1/2
/ (»),
(8.5)
Т = т 2 — С2
Тр т0 С0 •
Решение для погранслоя в окрестности квазифронтов легко получается из уравнения (5.8) и имеет вид [4,6]
Т1 = 2М
+ Аг(—у')Ф'
/ (»), У =
2(1 — V )2
V2п2Со
1/3
(то - Со)
(8.6)
где Аг(-) — функция Эйри.
Решение для погранслоя в окрестности фронта волны расширения строится с помощью двойного интегрального преобразования Лапласа по времени и Фурье — по пространственной продольной координате. Обращение ищется с помощью разложения изображения в бесконечный ряд по вычетам и метода профронтовой асимптотики (разложение по обратным степеням параметра преобразования Лапласа):
1
п2 ^ (п — 1/2)2
П=1 4 1 /
21/2 (п — 1/2)пп ^С^СКЬ — Со
/ (*)# (^ — Со
(8.7)
в ви
п
2
п
п
0
У
Уравнения для квазистатического погранслоя в рассматриваемом случае торцевого воздействия совпадают с уравнениями плоской задачи теории упругости для полуполосы [2]. Методы решения соответствующей краевой задачи представлены в [23, 24].
Отметим также, что коротковолновая составляющая имеет малоамплитудный характер и практического расчетного интереса не представляет. Ее анализ необходим только для доказательства корректности предложенной схемы расчленения.
В случае р = 0 изменяется только схема построения безмоментной составляющей [4]. Здесь применяется только метод прифронтовой асимптотики.
Методы обращения изображения по Лапласу, разработанные для цилиндрической оболочки, могут быть применены для общего случая оболочек вращения после нахождения решения в изображениях асимптотическим методом экспоненциального представления.
9. ВОЛНЫ ПРИ LM НАГРУЖЕНИИ
Рассмотрим схему расчленения НДС на составляющие при продольных воздействиях изгибающего типа (граничные условия (1.2)) на примере цилиндрической оболочки. Схема решения для изгибающего момента С1 в некоторый момент времени, больший времени пробега передним фронтом волны толщины оболочки, изображена на рис.9, а схема областей применимости приближенных теорий изображена на рис.10.
Рис. 9
Рис. lO
R1 — гиперболический погранслой в окрестности фронта волны расширения, R2 — квазиплоская квазиантисимметричная задача теории упругости, R3 — поперечное длинноволновое приближение (изгибная составляющая теории Кирхгофа - Лява), R4 — наложение двумерной составляющей и квазистатического погранслоя типа Сен-Венана, B1, B2 — области согласования.
Исследования работы [4] с учетом выделения в рассматриваемом типе НДС вместо динамического погранслоя — погранслоя в окрестности фронта и квазиплоской задачи, выполненного в [б], показали, что область согласования между изгибной и коротковолновой составляющими определяется асимптотическим неравенством
1/2 c_t / c_tX1/4
■ (9.l)
n1/2 (f) + n1-pÏ <<’w(ï)
В соответствии с (1.2) граничные условия для изгибной составляющей записываются в виде
2
G1
з h3 IH (t)/(Ö), w
О.
(9.2)
Рассмотрим сначала случай неосесимметричной нагрузки при р > 1/2. В этом случае при использовании интегрального преобразования Лапласа по временной переменной удается найти простую асимптотику для корней характеристического уравнения, содержащего малый параметр [4]. После использования метода перевала получаем решение для изгибающего момента при то >> п2р-1 в виде
G1 = з h31\
1 —
(1 — v)m* Со
2
np
exp
&\ з3/4n2Р-3/2Со Vm2 + 4n2-20pт02
—m.- I
n,) n1/223/2T« (m2 + 4^)
x (cos Tm — sin Tm )
1 + O
пто
С2 0
(9.3)
Tm
з1/2Со
+
2
m2 то
4^т0 ' 31/2п2р-1'
При т0 = О (п2р-1) (а также для всех значений времени в диапазоне т0 = О (1) при р < 1/2) решение строится при помощи разложения изображений в ряды по обратным степеням корня квадратного из параметра преобразования Лапласа, когда оригинал решения представляется в виде ряда по функциям £т,с и , введенным и описанным в [4]:
X
2
G1 = з h3 Jd1,c (y, то ) — n2r-2p n r [D2,c (y, то ) — D2,s (y, то )] —
,2r-2p (1 — V)m2 + n4r-4p_ m4С02
з1/2
8 ■ 31/2n2r.
D3,s(y, то) + ... >,
(9.4)
y=
3V4 Со 21/2 nr’
r =
1/2, p < 1/2, p, p > 1/2.
Решение для погранслоя в окрестности фронта волны расширения строится, как и в предыдущем случае воздействия типа ЬТ, с помощью двойных интегральных преобразований, теории вычетов и метода прифронтовой асимптотики [6]:
П=1
(9.5)
Решение для коротковолновой составляющей не строится вследствие его малоамплитудного характера.
10. ВОЛНЫ ПРИ Ш НАГРУЖЕНИИ
Рассмотрим схему расчленения НДС на составляющие при торцевом воздействии нормального типа (граничные условия (1.3)) также на примере цилиндрической оболочки. Схема расчленения для перерезывающей силы N изображена на рис. 11, а схема областей применимости приближенных теорий представлена на рис. 12.
к2
ІІІІІІБ,
ШБ Б2
1^4 1*5
Рис. 11
Применяются [19] следующие приближенные теории: — квазиплоская квазиантисимметричная
задача теории упругости (коротковолновое приближение), ^2 — погранслой в окрестности фронта волны сдвига, ^3 — погранслой в окрестности условного фронта поверхностной волны Релея, ^4 — поперечное длинноволновое приближение (изгибная составляющая теории Кирхгофа - Лява), ^5 — наложение двумерной составляющей и квазистатического погранслоя типа Сен-Венана. Здесь же обозначениями В, В2, В3, В4, В5 изображены соответствующие области согласования приближенных теорий. При этом анализ схемы решения показывает, что разрывной скачок наблюдается на фронте волны сдвига, а непрерывный скачок — на фронте условного фронта волны Релея, неверно воспринимаемого двумерной теорией типа Тимошенко как дополнительный фронт сдвиговой волны. Таким образом, неточность теории типа Тимошенко заключается здесь в том, что она неверно воспринимает условный фронт поверхностной волны Релея, фактически являющийся квазифронтом, как фронт с соответствующим разрывом решения, что полностью противоречит точной картине явления по трехмерной теории.
Для построения решения изгибной составляющей применяются методы, указанные в предыдущем параграфе. Решение в окрестности фронта волны сдвига строится методом прифронтовой асимптотики [6]. Для решения задачи о погранслое в окрестности условного фронта поверхностной волны Релея используется метод двойных интегральных преобразований: Лапласа — по времени и Фурье — по продольной координате [25].
ЗАКЛЮЧЕНИЕ
В работах У.К. Нигула [9] и других авторов показано, что все случаи нестационарных волн в тонкостенных телах при ударных торцевых воздействиях сводятся к трем основным схемам, определяемым тремя типами нагрузки.
Описанные в статье асимптотические методы исследования нестационарных волн в тонких упругих оболочках вращения составляют в совокупности замкнутую математическую модель исследования простейшего класса нестационарных процессов в тонкостенных конструкциях. Замкнутую как в смысле полноты разработанных математических методов, так и в смысле влияния особенностей рассмотренных в классе материалов и геометрии конструкции на свойства НДС и соответственно на принципы расчленения НДС на составляющие с различными показателями изменяемости. Полнота математических методов заключается как в построении полной системы асимптотически оптимальных уравнений, описывающих все необходимые составляющие нестационарного НДС, так и в разработке аналитических методов решения краевых задач для составляющих. Последние представляют как самостоятельный интерес, так и интерес для построения решений модельных задач при разработке численных методов исследования или для разработки гибридных аналитически-численных алгоритмов.
В рамках разработанной математической модели проведены дальнейшие исследования, углубившие понимание рассматриваемых физических явлений с точки зрения формирования установившихся процессов при различных видах воздействий. Так, в работе [15] исследуется, в частности, стремление решения для изгибной составляющей в оболочках вращения к стационарному с ростом времени и принципиальные отличия со временем свойств этих решений от соответствующих решений для пластин. В работах [16, 17] исследовано нестационарное НДС при ударно приложенных гармонических торцевых воздействиях на оболочки вращения и выявлено расположение зон перехода нестационарной части решения в гармонические колебания.
Разработка математической модели указанного класса определила принципы подхода к исследованию самых разнообразных динамических задач для тонкостенных конструкций, которые базируются на широчайших возможностях асимптотического метода. Приведем примеры последних разработок, обобщающих рассмотренный асимптотический подход.
В плане изучения нового класса ударных воздействий и изменения общей схемы расчленения характерна работа [18], которая указывает, в частности, что при высокочастотном краевом нагружении НДС около торца будет описываться вместо длинноволновой низкочастотной составляющей (теория Кирхгофа - Лява) длинноволновой высокочастотной. Локальное нормальное ударное воздействие сжимающего типа на упругую плиту рассматривалось в [19]: для решения этой очень характерной задачи использовалась уточненная асимптотическая 2D модель. Отметим также, что изменение геометрического класса задач — переход к оболочкам общего очертания (или к воздействиям, локально действующим на торцы оболочки) — требует кардинально нового подхода к выбору систем координат, связанного со сложной геометрией фронтов и квазифронтов. В [18] предложено срединную поверхность произвольных оболочек относить не к линиям кривизны, а к полугеодезической системе координат и изучать процесс распространения волн вдоль геодезических. Данный подход реализован в работах [20-21] на примере скошенной цилиндрической оболочки.
Интересными теоретически и практически важными явились исследования распространения нестационарных волн в подкрепленных оболочках вращения [22, 23]. Здесь необходимым оказался учет квазистатического погранслоя типа Сен-Венана при исследовании волнового процесса в подкреплениях.
Рассмотрены нестационарные процессы в оболочках из других материалов со свойствами, принципиально влияющими на характер НДС. Так, в работах [7, 24] выведены асимптотически оптимальные уравнения для различного типа погранслоев в трансверсально изотропных тонких пластинах и оболочках.
В работах [7, 25] рассмотрено нестационарное НДС в оболочках из вязкоупругих материалов. Введение нового параметра — показателя интенсивности времени релаксации — позволило описать новую классификацию нестационарного НДС конструкций в случае вязкоупругих материалов с дифференциальной формой уравнений состояния. Выведены асимптотически оптимальные уравнения всех составляющих.
Анализ математической модели, описанной в рассматриваемой статье, и ее обобщений на различные случаи геометрии тонкостенных конструкций и материалов показывает универсальный характер разработанной математической теории.
Библиографический список
1. Коссович Л.Ю., Каплунов Ю.Д. Асимптотический анализ нестационарных упругих волн в тонких оболочках вращения при ударных торцевых воздействиях // Изв. Сарат. ун-та. Новая сер. 2001. Т. 1, вып. 2. С. 115128
2. Гольденвейзер А.Л. Теория упруги тонких оболочек. М.: Наука, 1976. 512 с.
3. Гольденвейзер А.Л., Лидский В.Б., Товстик П.Е. Свободные колебания тонких упругих оболочек. М.: Наука, 1979. 384 с.
4. Коссович Л.Ю. Нестационарные задачи теории упругих тонких оболочек. Саратов: Изд-во Сарат. ун-та, 1986. 176 с.
5. Kaplunov J.D. On the quasi-front in two-dimensional shell theories // C.R. Acad. Sci. Paris. Ser. 2. 1991. V. 313. P. 731-736.
6. Kaplunov J.D., Kossovich L.Yu., Nolde E.V. Dynamics of thin walled elastic bodies. San Diego Academic Press, 1998.
7. Шевцова Ю.В. Погранслой в окрестности квазифронта в транвсерсально изотропной цилиндрической оболочке // Проблемы прочности элементов конструкций под действием нагрузок и рабочих сред: Межвуз. науч. сб. Саратов: Изд-во СГТУ, 2000. С. 114-117.
8. Бажанова Н.С., Коссович Л.Ю., Сухоловская М.С. Нестационарные волны в вязкоупругих оболочках: модель Максвелла // Известия вузов. Северокавказский регион. Естественные науки. 2000. №2. С. 17-24.
9. Nigul U. Regions of effective aplication of the methods of three-dimensional and two-dimensional analysis of transient stres waves in shells and plates // Intern. J. Solids and Structures. 1969. V. 54. P. 607-627.
10. Каплунов Ю.Д., Коссович Л.Ю. Асимптотическая модель для вычисления дальнего поля волны Рэлея в случае упругой полуплоскости // Докл. АН. 2004. Т. 395, № 4. С. 482-484.
11. Kaplunov J., Kossovich L., Zakharov A. An explicit asymptotic model for the Bleustein-Gulyaev wave. C.R. Mecanique. 2004. № 332. P. 487-492.
12. Коссович Л.Ю., Кушеккалиев А.Н. Поле Рэлея в бесконечном упругом слое // Математика, механика: Сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2004. Вып.5. С. 159-161.
13. Ковалев В.А., Коссович Л.Ю., Таранов О.В. Поле Рэлея в задаче Лэмба для цилиндрической оболочки // Известия вузов. Северо-Кавказский регион. Естественные науки. Спецвыпуск. 2004. С. 52-54.
14. Ковалев В.А., Коссович Л.Ю., Таранов О.В. Дальнее поле волны Рэлея для упругой полуполосы при дей-
ствии торцевой нагрузки // Известия РАН. МТТ. 2005. №5. С. 89-96.
15. Кушеккалиев А.Н. Решение задач о распространении волн в трансверсально-изотропной цилиндрической оболочке при нормальных воздействиях // Механика деформируемых сред. Саратов: Изд-во Сарат. ун-та, 2002. Вып. 14. С. 106-115.
16. Кушеккалиев А.Н. Волны типа Рэлея в полубеско-нечной пластине при нормальном воздействии поперечного типа // Механика деформируемых сред. Саратов: Изд-во Сарат. ун-та, 2004. Вып. 15. С.66-73.
17. Коссович Л.Ю., Кушеккалиев А.Н. Анализ приближений в задаче Лэмба для бесконечного упругого слоя // Известия вузов. Северо-Кавказский регион. Естественные науки. 2003. №2. С. 10-22.
18. Коссович Л.Ю., Кушеккалиев А.Н. Расчленение нестационарного НДС в задаче Лэмба для бесконечного слоя на составляющие с разными показателями изменяемости // Тр. III Всероссийской конференции по теории упругости с международным участием. Ростов н/Д, 2003. С. 232-234.
19. Ковалев В.А., Таранов О.В. Расчленение нестационарного НДС цилиндрических оболочек при ударных торцевых воздействиях нормального типа // Смешанные задачи механики деформированного тела: Материалы V Рос. конф. с международным участием. Саратов: Изд-во Сарат. ун-та, 2005. С. 191-193.
20. Кириллова И.В. Области согласования погранслоя и коротковолнового высокочастотного приближения // Математическое моделирование и управление в технических системах: Сб. трудов / СГАУ. Саратов, 1998. С. 3-11.
21. Новожилов В.В., Слепян Л.И. О принципе Сен-Венана в динамике стержней // ПММ. 1965. Т. 29, вып. 2. С. 261-281.
22. Слепян Л.И. Нестационарные упругие волны. Л.:Судостроение, 1972. 374 с.
23. Гусейн-Заде М.И. Об условиях существования затухающих решений плоской задачи теории упругости для полуполосы // Прикл. мат. и мех. 1965. Т.29, вып. 2. С. 393-399.
24. Гусейн-Заде М.И. О необходимых и достаточных условиях существования затухающих решений плоской задачи теории упругости для полуполосы // Прикл. мат. и мех. 1965. Т.29, вып. 4. С. 752-760.
25. Ковалев В.А., Таранов О.В. Анализ точного и приближенного решения для погранслоя в окрестности условного фронта поверхностной волны Рэлея в упругой полуполосе // Вестник СамГУ. Естественнонаучная сер. 2007. № 6(56). С. 57-60
