Elliptic boundary layer in shells of revolution under surface shock loading of normal type
- 作者: Kirillova I.V.1
-
隶属关系:
- Saratov State University
- 期: 编号 5 (2024)
- 页面: 48–59
- 栏目: Articles
- URL: https://journals.rcsi.science/1026-3519/article/view/277074
- DOI: https://doi.org/10.31857/S1026351924050045
- EDN: https://elibrary.ru/UBBQQG
- ID: 277074
如何引用文章
全文:
详细
This article presents a method for solving the boundary-value problem for an elliptic boundary layer occurring in thin-walled shells of revolution under impact loads of normal type applied to the face surfaces. The elliptic boundary layer is formed in the vicinity of the conditional front of Rayleigh surface waves and is described by elliptic equations with boundary conditions determined by hyperbolic equations. In the general case of shells of revolution, methods for solving elliptic boundary layer equations developed for shells of revolution with zero Gaussian curvature cannot be applied. The previously considered approach using Laplace and Fourier integral transforms fails because the governing equations become equations with variable coefficients. The method proposed in this article for solving the equations of the elliptic boundary layer is based on the use of asymptotic representations of the Laplace-transformed solutions (in time) in exponential form. Numerical calculations of normal stresses based on the obtained analytical solutions are provided for the case of a spherical shell.
全文:
1. Введение. При ударных воздействиях нормального типа NW (в соответствии с классификацией У.К. Нигула [1, 2]), возникает эллиптический погранслой, определяющий особый тип нестационарного напряженно-деформированного состояния (НДС) в тонких пластинах и оболочках в окрестности условного фронта поверхностных волн Рэлея. Данный погранслой и выделяет нормальный тип ударного воздействия NW среди других видов ударных воздействий – продольного, изгибающего типа (LM) и продольного, тангенциального типа (LT).
Асимптотический подход к выводу уравнений теории оболочек был впервые разработан в трудах А.Л. Гольденвейзера [3, 4]. Концепция изменяемости напряженно-деформированного состояния (НДС) по координатам и её количественная оценка через показатели изменяемости позволила сформулировать принципы асимптотического интегрирования точных трёхмерных уравнений теории упругости в статике и получить этим методом уточнённые уравнения теории оболочек. Также удалось выделить составляющие НДС (безмоментная составляющая, простой краевой эффект и др.) как в статике, так и стационарной динамике, а также построить асимптотические методы их определения.
Разработанный для задач статики и стационарной динамики асимптотический подход дал возможность создать принципиально новую асимптотическую теорию нестационарных процессов деформации тонких оболочек [5–8], основанную на варианте метода сращиваемых разложений, позволяющего разделить фазовую область на различные участки, обладающие теми специфическими значениями показателей изменяемости и динамичности, которые соответствуют известным особенностям поведения нестационарных волн. В качестве таких особенностей, отличающих нестационарные задачи от задач статики и стационарной динамики, являются различные виды фронтов и квазифронтов, отражающих ударный характер приложения нагрузки.
В рамках разрабатываемой асимптотической теории нестационарных волн появилась необходимость нахождения новых асимптотических методов вывода приближенных теорий из точных трёхмерных уравнений теории упругости. Асимптотические уравнения двумерных составляющих, соответствующих классической теории Кирхгофа–Лява, получены в работе [5]. Общий подход к выводу уравнений для НДС в окрестностях фронтов волн и волновых квазифронтов изложен наиболее полно в работах [6–8]. В итоге были выделены безмоментная и изгибная составляющая двумерной теории Кирхгофа–Лява, гиперболические погранслои в окрестностях фронтов волн расширения и сдвига, параболический погранслой в окрестности квазифронта, связанного с ложным фронтом волны растяжения-сжатия по двумерной теории Кирхгофа–Лява и эллиптический погранслой в окрестности условного фронта поверхностных волн Рэлея. При различных типах ударных воздействий эти компоненты применяются в разных комбинациях, но полнота таких представлений доказывается наличием и определением расположения так называемых областей согласования соседних составляющих.
Проведённые исследования показали, что ни одна из известных двумерных теорий не может описать полностью рассматриваемую волновую картину вследствие высокой изменяемости и принципиально разного поведения нестационарного НДС в малых окрестностях фронтов и квазифронтов.
В работах [8, 9] окончательно описано построение асимптотически оптимальных уравнений эллиптического погранслоя в окрестности условного фронта поверхностных волн Рэлея. Выявлены уникальные свойства этого погранслоя: поведение решения вне лицевых поверхностей оболочки описывается эллиптическими уравнениями, а на поверхности для потенциальных функций, определяющих решение, описывается гиперболическими уравнениями. Отметим, что такие качественные свойства решения базируются на анализе дальнего поля волны Рэлея в случае упругой полуплоскости в классической задаче Лэмба [10], объясняющего эффект так называемых поверхностных волн. В работе [9] также получено решение для эллиптического погранслоя в случае цилиндрической оболочки: используются интегральные преобразования Лапласа по времени и Фурье по продольной координате.
В общем случае оболочек вращения, когда разрешающие уравнения содержат переменные коэффициенты, метод их решения, разработанный в [9], неприменим. В данной статье описывается новый разработанный асимптотический метод решения искомых уравнений в случае нормальной ударной нагрузки на лицевые поверхности произвольных оболочек вращения. Используются асимптотические методы в пространстве интегрального преобразования Лапласа по времени. В качестве примера приведены результаты расчёта нормального напряжения для сферической оболочки.
2. Постановка задачи. Рассмотрим оболочку вращения, изображенную на рис.1, где система координат (, , z) привязана к срединной поверхности: – длина дуги вдоль образующей, – угол в окружном направлении, z – внешняя нормаль к срединной поверхности. Напряжения и перемещения оболочки обозначим как ij и vi (i, j = 1, 2).
Рис. 1. Полубесконечная оболочка вращения.
Для определенности рассмотрим осесимметричное НДС.
Граничные условия на лицевых поверхностях оболочки выберем в форме, аналогичной [9]:
где h – полутолщина оболочки, t – время, L – ширина кольца нагрузки на поверхности. Ударное воздействие P(a, t) определяется следующим образом
Здесь p – амплитуда нагрузки, а H(t) – единичная функция Хевисайда. При этом рассматриваются однородные начальные условия
Разрешающие уравнения эллиптического погранслоя при рассматриваемых воздействиях описаны в [8]. В исходной форме они записываются относительно объёмных и сдвиговых потенциалов i, i(i = 1, 2), когда потенциальные функции 1, 1 определяют распространение возмущений в положительном направлении оси z, а 2 и 2 – в отрицательном. Ограничиваясь рассмотрением только волны, инициируемой лицевой поверхностью z = -h, выпишем, соответственно, разрешающую систему уравнений только для потенциалов 1, 1 (индекс «1» опускаем) в безразмерной форме в координатах = a/h, z = z/h, = c2t/h (c2 – скорость волны сдвига). Тогда уравнения, описывающие НДС рассматриваемого погранслоя, имеют форму
(2.1)
и являются уравнениями эллиптического типа, где c1 и cR – скорости волн расширения и поверхностных волн Рэлея; v – коэффициент Пуассона. Граничные условия на лицевой поверхности описываются уже гиперболическими уравнениями
(2.2)
где E – модуль Юнга, = h/R – малый параметр тонкостенности оболочки, R – характерное значение радиусов кривизны, B – расстояние до оси вращения, а постоянные g, kc, определяются выражениями
Отметим, что поведение НДС в области рассматриваемого погранслоя можно наглядно анализировать на примере нормального напряжения 33, выражающегося через потенциалы следующим образом:
3. Асимптотическое решение для эллиптического погранслоя в оболочках вращения при ударном воздействии нормального типа на лицевые поверхности. В цитируемой работе [8] приведено решение для рассматриваемого эллиптического погранслоя в случае цилиндрической оболочки. В этом случае уравнения (2.1)–(2.2) становятся уравнениями с постоянными коэффициентами, и для их решения были применены интегральные преобразования Лапласа по временной, а Фурье – по продольной переменной.
Разрешающие уравнения эллиптического погранслоя для рассматриваемого общего случая оболочек вращения имеют переменные коэффициенты, и поэтому вышеизложенная методика для нашего общего случая не подходит. Однако переменные коэффициенты являются при этом медленно изменяющимися по продольной координате. Это позволяет использовать в нашей работе, в качестве ключевого, метод экспоненциального представления в пространстве преобразования Лапласа по временной переменной. Главной особенностью предлагаемого здесь асимптотического подхода является определение выражения для изображения по Лапласу искомого решения в рассматриваемом случае оболочек вращения через известное решение для цилиндрической оболочки.
Рассмотрим сначала поведение потенциала на границе = -1, определяемое первым уравнением системы (2.2), которое для его изображения по Лапласу примет вид
(3.1)
где s – параметр преобразования Лапласа.
Проанализируем свойства решения для двойного интегрального преобразования (Лапласа по времени и Фурье по продольной координате) граничного значения потенциальной функции 0 в случае цилиндрической оболочки. Оно легко получается из (3.1) при B(0) = const:
(3.2)
Обращение преобразования Фурье для изображения (3.2) даёт, в соответствии с [11], следующее выражение для 0L:
(3.3)
В решении (3.3) слагаемые с определяют волны, возбуждаемые границами приложения ударной поверхностной нагрузки 0 = ±L0. Так, в области 0 > 0 граница 0 = L0 возбуждает прифронтовое поле в окрестности условного фронта поверхностных волн Рэлея 0 = kRt0 + L0 (слагаемое exp{-[(0 - L0)/(kRe)]s} ) и прифронтовое поле в окрестности фронта 0 = L0 – kRt0 (слагаемое exp{[(0 - L0)/(kRe)]s} ). В этой же области x0 > 0 граница 0 = -L0 возбуждает прифронтовое поле в окрестности фронта 0 = kRt0 - L0 (слагаемое exp{-[(0 + L0)/(kRe)]s} ).
Отметим, что решение (3.3) можно непосредственно получить из (3.1) как линейную комбинацию частного решения неоднородного уравнения и общих решений уравнения однородного, определяющих вышеописанные волновые поля. Проведенный анализ позволяет построить асимптотическое решение уравнения (3.1) для общего случая оболочек вращения. Так, решения для волновых составляющих можно получить из однородного уравнения
(3.4)
методом экспоненциальных представлений с погрешностью O(e) в виде
(3.5)
где (s)0 /e – специальный вид функции изменяемости, C(s) – постоянная интегрирования, а (0, s) – медленно изменяющаяся функция интенсивности, представляемая разложением в степенной ряд по малому параметру .
Подставляя (3.5) в (3.4), получаем с погрешностью O(e2) следующее уравнение
которое даёт выражение для функции изменяемости
и нулевое приближение для функции интенсивности
Поскольку частное решение неоднородного уравнения (3.1) с погрешностью O(e) записывается в виде
полностью соответствующему такому же частному решению в (3.3), то общее решение уравнения (3.1) может быть получено при использовании волновых решений (3.5) за счёт выбора постоянных интегрирования C(s) при условии его непрерывности в точках 0 = ±L0. В итоге получаем следующее решение:
(3.6)
Полученное решение (3.6) полностью обобщает соответствующее решение (3.3) для случая цилиндрической оболочки. Поскольку нас интересуют асимптотики только в малых окрестностях условных фронтов волн Рэлея, то решение (3.6) можно переписать в форме
(3.7)
Таким образом, из вида асимптотики решения (3.7) следует, что она связана с соответствующим решением для цилиндрической оболочки следующим образом:
(3.8)
где L0,c является соответствующим решением для цилиндрической оболочки.
Перейдём к определению потенциалов , в основной области. Анализируем систему (2.1)–(2.2) и выражение (2.3), приходим к выводу, что только первое уравнение в граничных условиях (2.2) содержит в асимптотически главной части волновой оператор, определяющий НДС в окрестностях условных фронтов поверхностных волн Рэлея. Следовательно, только в этом уравнении нужно учитывать величины нулевого и первого порядка малости; в остальных же уравнениях достаточно удерживать только асимптотически главные величины. Следовательно, аналогично зависимости (3.8), все компоненты НДС могут быть представлены в виде
где ij,c, vLi,c соответствующие решения для случая цилиндрической оболочки. Следовательно, аналогичные зависимости имеют место и для оригиналов:
Окончательно, используя результаты решения задачи об эллиптическом погранслое для цилиндрической оболочки [8], получаем следующее решение для напряжения 33:
(3.9)
Асимптотический анализ решения (3.9) показывает, что вне малых окрестностей условного фронта волны Рэлея его порядок уменьшается до O(). Отметим также следующую особенность рассматриваемого решения: вне узкой прифронтовой зоны эллиптического погранслоя [8, 9] порядка O(2) решение для 33 имеет порядок O(1) только в малой приторцевой зоне + 1 = O(). Следовательно, важнейшее предположение о достаточности отдельного рассмотрения волн, инициируемых каждой лицевой поверхностью, полностью оправдывается при построении решения в этой области:
где, как и показывают численные расчёты, имеет место первый скачок напряжения 33.
4. Анализ решения для эллиптического погранслоя в случае сферической оболочки. Рассмотрим свойства решения для нормального напряжения 33 в зависимости от пространственных переменных и времени на примере сферической оболочки. На рис. 2 изображена схема сечения срединной поверхности такой оболочки плоскостью, проходящей через центр сферы. Для наглядности определения зависимости B (расстояния от текущей точки срединной поверхности до оси вращения) введены следующие обозначения: точка T – пересечение торцевой поверхности со срединной линией; точка A соответствует координате a0; точка O2 – граничная точка оси вращения; a0 – текущая координата; углы T и a соответствуют радиусам до точек T и A. Поскольку a = a/R, получаем выражение для B(a0) через a0 и угол T
(4.1)
Рис. 2. Геометрия сечения сферической оболочки.
Выражение (4.1) даёт возможность полностью определить закон распределения напряжения 33 по координатам и времени в области действия эллиптического погранслоя.
На рис. 3. изображены графики нормального напряжения s33 в малой окрестности условного фронта волны Рэлея 0 = kRt0 + L0, возбуждаемого границей приложения нагрузки 0 = L0. Здесь кривые 1, 2, 3 соответствуют значениям нормальной координаты z = -0.99, -0.96, -0.92. Расчёты проведены для момента безразмерного времени t0 = 1. Геометрические и механические параметры принимались следующими: L0 = 0.2, T = 0.93, n = 0.3, kR = 0.93.
Рис. 3. График нормального напряжения 33 в малой окрестности условного фронта волны Рэлея 0 = kRt0 + L0 в момент времени t0 = 1 для значений нормальной координаты = -0.99,-0.96,-0.92.
Расчёты подтверждают выводы предшествующего параграфа: первый скачок напряжения сосредоточен в узкой прифронтовой зоне и полностью сосредоточен в малой окрестности лицевой поверхности.
5. Заключение. В данной статье впервые представлен специфический асимптотический метод решения задачи для эллиптического погранслоя в произвольных оболочках вращения при ударных поверхностных нагрузках нормального типа. Он основан на решении на первом этапе базовой задачи для граничного значения потенциальной функции методом экспоненциальных представлений в пространстве интегрального преобразования Лапласа по временной переменной с последующим представлением компонент НДС через соответствующие компоненты эталонной задачи для цилиндрической оболочки. Проведённые численные расчёты полностью подтверждают правильность предложенного асимптотического подхода.
参考
- Nigul U. Regions of effective application of the methods of three-dimensional and two-dimensional analysis of transient stress waves in shells and plates // Int. J. Solid and Structures. Vol. 5. № 6. 1969. P. 607–627.
- Nigul U. Comparison of the results of the analysis of transient wave processes in shells and plates according to the theory of elasticity and approximate theories // Applied Mathematics and Mechanics. 1969. Vol. 33. Issue 2. Pp. 308–332.
- Gol’Denveizer A.L. Theory of Elastic Thin Shells // Pergamon. 1961, p. 680.
- Gol’Denveizer A.L., Lidskiy V.B., Tovstik P.E. Free vibrations of thin elastic shells. 1979. 384 p.
- Kaplunov Yu.D., Kirillova I.V., Kossovich L.Yu. Asymptotic integration of the dynamic equations of the theory of elasticity for the case of thin shells. J. Appl. Math. Mech. Volume 57, Issue 1, pp. 95–103.
- Kossovich L.Yu. Non-stationary problems of the theory of elastic thin shells. Saratov, 1986, 176 p.
- Kaplunov J.D., Nolde E. V., Kossovich L.Y. Dynamics of thin walled elastic bodies // Academic Press. 1998. P. 226.
- Kirillova I.V., Kossovich L.Y. Asymptotic theory of wave processes in shells of revolution under surface impact and normal end actions Mechanics of Solids. 2022. Т. 57. № 2. P. 232–243.
- Kirillova I.V., Kossovich L.Y. refined equations of elliptic boundary layer in shells of revolution under normal shock surface loading. Vestnik of the St. Petersburg University: Mathematics. 2017. Т. 50. № 1. С. 68–73.
- Kaplunov Yu.D., Kossovich L.Yu. Asymptotic model for calculating the far field of a Rayleigh wave in the case of an elastic half-plane. Reports of the Academy of Sciences. 2004. Vol. 395. No. 4. Pp. 482–484.
- Polyanin A.D., Manzhirov A.V. Handbook of integral equations. 2nd ed. Chapman & Hall/CRC Press. 2008. P. 1108.
补充文件
