Нестационарные волны в пластинах при нормальных ударных воздействиях

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

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

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

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

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

В случае простых геометрических объектов (напр., плит) алгоритм степенных рядов может быть успешно применен в форме символического метода А. И. Лурье [54,55] или в форме метода начальных функций В.З. ^ Власова [11,12]. Символические формулы трехмерной теории динамики плит представлены в работе [62] для деформации, антисимметричной относительно срединной поверхности. В символической теории применяются бесконечные, во всех других вариантах — усеченные ряды.

Состояние исследований по изучению переходных волновых процессов деформации изгибного типа, вызванные действующей нагрузкой в бесконечных и полубесконечных плитах и балках, охарактеризовано в обзорах [1,8,10,14,30,40,57,61,64,70,79,86,91,100].

На основе элементарной теории изгиба решены многочисленные задачи расчета плит и оболочек. В динамике балок первая поправка к элементарной теории была внесена еще Релеем, предложившим учесть инерцию вращения. С. П. Тимошенко показал в 1921 г., что учет сдвига дает поправку такого же порядка. Он вывел уравнение балки, учитывающее поправки от сдвига и от инерции вращения [104]. Эти поправки приводят к # отказу от элементарной теории (прямолинейная нормаль к срединной поверхности не остается нормалью после деформации).

Я.С. Уфлянд был первым, применившим [79] теорию типа Тимошенко к анализу переходных волновых процессов, вызванных сосредоточенной импульсной нагрузкой в бесконечной плите и балке. Он использовал преобразование Лапласа. Поскольку полученные интегралы Меллина (в отличие от теории Кирхгоффа) не поддались обращению в замкнутом виде, ц то Я. С. Уфлянд предложил приводить их к суммам определенных интегралов по берегам срезов плоскости s, обеспечивающих однозначность подынтегральных функций. На основе этих интегралов в [79] были получены некоторые численные данные.

Позже различные аспекты применения теории плит и балок типа Тимошенко рассматривались в многочисленных работах [8,30,58,79,84,85,86,91,99], в том числе задача расчета бесконечной плиты ^ под действием сосредоточенной силы — в работах [30,99].

В работах [8,30,58,79,84,85,86,91,99] были решены конкретные задачи при помощи операционного исчисления. В основном использовалось преобразование Лапласа [8,30,58,79,84,86,99], но нашло применение и преобразование Фурье по координате [86,91], при построении формальных решений в виде контурных интегралов, которые были обращены методом перевала. Методом интегрирования по берегам срезов плоскости s получен ряд диаграмм в работе [84].

В работах [8,9,58] приведены результаты, полученные методом сеток с выделением частных решений, переносящих разрывы первых и вторых производных нормального прогиба и угла поворота.

Проблема применимости теории типа Тимошенко рассматривалась на основе трехмерной теории в работе [61]. Было установлено, что теория типа Тимошенко неприменима в области от фронта волны сжатия до условного # фронта поверхностных волн Релея, а также в непосредственной близости от последней. В [61] приведен также ряд численных данных о точности теории типа Тимошенко в области, где её можно считать практически применимой.

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

Ю. Микловиц [98,100] исследовал осесимметричный волновой процесс, возбужденный в бесконечной плите действием двух сосредоточенных сил, приложенных к свободным поверхностям плиты и направленных друг к другу. Были применены преобразование Лапласа по времени и преобразование Ханкеля по координате. Первое обращение выполнялось при помощи вычетов, второе — методом перевала (с точностью трех первых контурных интегралов). В данном случае возбуждаются только моды, симметричные относительно срединной поверхности. При использовании метода перевала был использован малочастотный участок диаграммы групповых скоростей первых трех симметричных мод, где максимальное значение групповой скорости соответствует первой моде в точке со = 0, и равен величине с^ (4.6). В работе [100] линия г = cj называется «фронтом» и указывается, что вблизи «фронта» доминируют малочастотные волны первой моды. Фактически фронт перемещается со скоростью с^, и в области cj > г > с J движутся высокочастотные волновые группы более высоких мод.

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

В работе [61] рассматривался волновой процесс изгиба полубесконечной плиты, возбужденной напряжениями, внезапно приложенными по толщине плиты. Были применены преобразование Лапласа по времени и синуси косинуспреобразования Фурье по координате. Преобразование Фурье было обращено при помощи вычетов, а преобразование Лапласа — методом перевала с учетом шести первых контурных интегралов, которые в данном случае связаны с антисимметричными модами. Роль отдельных мод зависит от рассматриваемой величины, но, грубо говоря, вклад первых двух мод доминирует при л < с2/ (с2 — скорость распространения волн сдвига), а при x>c2t быстро растет вклад следующих мод, что ограничивало (при учете шести первых контурных интегралов) область исследования условием x<.2c2t. Численные результаты [61] позволили сделать выводы о применимости двумерных теорий.

В работе [92] рассматривается плоская задача для плоскопараллельного упругого слоя со свободными границами. Источником колебаний является внезапно приложенная поперечная сила (рассматривается несколько вариантов распределения заданных напряжений сдвига). Силы распределены по сечению так, что возбуждаются только изгибные волны, смещения в которых антисимметричны относительно срединной линии слоя. Напряжения и смещение в слое представляются в виде суперпозиции нормальных волн. Для четырех первых нормальных волн рассчитаны кривые зависимости фазовой и групповой скоростей от волнового числа. Для оценки встречающихся интегралов используется метод стационарной фазы. Показано, что при 0 0.37^/ отмечается существенное различие в результатах точной теории, элементарной теории и теории Тимошенко. Рассматривается поле напряжений при больших t в районе фронта волны сдвига и фронта поверхностных волн Релея. При воздействии силы вдоль поверхности полосы в промежутке с, напряжения на поверхности y=h испытывает скачок в районе фронта волны Релея. В случае, распределение силы вдоль нормальной координаты изменяется по параболическому типу.

• Динамике стержней и пластин посвящены работы В. В. Новожилова и Л. И. Слепяна [68,74,75]. В них изучались свойства аналитических решений задач при рассмотрении переходных волновых процессов в стержнях и пластинах, анализировалась работа принципа Сен-Венана. Было показано, что главная часть деформаций, которая соответствует внезапно приложенным самоуравновешенным по сечению нагрузкам, локализируется вблизи фронтов волн и того сечения, где приложена нагрузка. Проведенный анализ областей действия различных теорий позволил использовать для решения задач в начале движения и в окрестностях фронтов волн разные методы.

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

Направление исследований, связанное с асимптотическим интегрированием уравнений трехмерной теории упругости, развивалось в работах А. Л. Гольденвейзера, А. Грина, Б. Новотны и др.

В работах Г. И. Петрашеня и Л. А. Молоткова [72,73] методом асимптотического анализа был исследован вопрос о границах применимостей уравнений классической теории изгиба плит в динамических задачах. Основная идея, использованная в [70], состоит в сопоставлении формальных решений (контурных интегралов), полученных на основе трехмерной теории щ и двумерных теорий. В ходе этого сопоставления пользуются разложением точного изображения по степеням нормальной координаты z.

Основополагающие понятия показателя изменяемости напряженно-деформированного состояния (НДС) по пространственным координатам и операции растяжения масштаба в уравнениях теории упругости связаны, в первую очередь, с работами АЛ. Гольденвейзера [14−24]. При рассмотрении статических задач, посвященных построению двумерной теории оболочек, ф вводился малый безразмерный параметр, равный отношению толщины оболочки к характерному радиусу.

Введение

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

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

Динамике тонких упругих оболочек посвящены работы.

A.Л.Гольденвейзера [17,21]. В работе [17] для интегралов двумерных динамических уравнений теории оболочек была введена классификация, аналогичная статическому случаю. Было показано, что при построении классификации в динамике необходимо учитывать изменяемость напряженного состояния по времени. В работе [22] рассмотрены динамические трехмерные уравнения теории упругости и свойства их интегралов в случае, когда тело тонкое и его лицевые поверхности не закреплены. Получена связь данных интегралов и интегралов двумерных уравнений теории оболочек и теории погранслоя. В работе [24] сформулирован модифицированный принцип Сен-Венана, обуславливающий затухание асимптотически главной части НДС, вызванной системой сил, приложенных к торцу тонкого упругого тела. Получены условия выполнения модифицированного принципа Сен-Венана и изучена возможность их использования при построении итерационных процессов интегрирования общих уравнений теории упругости.

Асимптотический метод был применен также для изучения свободных колебаний оболочек на основе двумерной классической теории оболочек. Им посвящены работы А. Л. Гольденвейзера, Ю. Д. Каплунова, Л. Ю. Коссовича,.

B.Б.Лидского, Е. П. Товстика и др.

Колебаниям оболочек посвящены работы А. Л. Гольденвейзера и Ю. Д. Каплунова. В [21] в рамках трехмерной теории упругости рассматривалась задача об установившихся колебаниях тонкой упругой осесимметричной оболочки вращения произвольного очертания под действием краевой нагрузки, меняющейся по гармоническому закону. Исследовался вопрос об использовании равноизменяющихся решений уравнений теории упругости для приближенного использования вынужденных колебаний оболочек при частотах, исключающих применение теории Кирхгофа-Лява. В [20] было установлено, что (при достаточно больших значениях частоты вынужденных колебаний) двумерная теория оболочек становится неприемлемой, т. е. существуют критические значения частот, при превышении которых изменяемость в меридиональном направлении оказывается существенно большей, чем в направлении параллелей.

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

Одними из первых работ, связанных с изучением динамического НДС оболочек и использованием метода расчленения, были работы Н. А. Алумяэ, JI. Поверуса, У. К. Нигула [3−6,59,61,66,67,103]. В [3] исследовался осесимметричный переходный процесс в полубесконечной круговой цилиндрической оболочке, вызванный действием краевой нагрузки, меняющейся во времени по синусоидальному закону. В [4] краевая нагрузка на оболочку задавалась во времени функцией Хевисайда, а по дуговой координате менялась по закону косинуса. В [3] асимптотическое обращение контурных интегралов от изображений решений по Лапласу позволило разложить НДС на безмоментное решение и краевые эффекты. В [4] расчленение напряженных состояний проводилось с учетом показателя изменяемости по времени.

Работа У. К. Нигула [66] посвящена изучению начального этапа переходных процессов деформации упругих круговых цилиндрических оболочек, вызванных торцевой нагрузкой, которая действует или возрастает до максимального значения за время, которое соизмеримо со временем пробега волнами деформаций пути, равного характерному размеру срединной поверхности оболочки или меньше этого времени. На основании анализа переходных процессов деформации плит была получена следующая качественная картина деформации: в момент приложения ударной нагрузки к торцу в оболочке возникает сложная система первичных волн, как продольных, так и поперечных, которые начинают распространяться вглубь оболочки. При этом первичные волны, распространяющиеся в продольном направлении, взаимодействуя с лицевыми поверхностями, отражаются от них и, в свою очередь, порождают вторичные волны, причем каждая из волн порождает и продольную, и поперечную отраженную волну. Таким образом, в оболочке возникает сложная система волн. Поскольку с течением времени количество волн катастрофически растет, отдельное построение и анализ числа элементарных волн практически неосуществимы. Поэтому при расчете тонкостенных конструкций следует изучать вклад не отдельно взятой волны, а суммировать вклад всех волн пакета. В данной работе выявлены области и условия применимости приближенных теорий для аппроксимации суммарного вклада отдельных волн, проведено численное сопоставление решений по теории упругости и по приближенным теориям. В [66,67,103] приведены результаты исследований по изучению областей применимости теории Кирхгофа-Лява и теории типа Тимошенко при осесимметричной деформации оболочек вращения, вызванной локальной нагрузкой, проанализированы результаты использования методов перевала, стационарной фазы, конечных разностей и прифронтовой асимптотики. Один из главных выводов, сформулированных в данных работах, заключается в том, что характер НДС существенно зависит от изменяемости воздействия во времени.

В работе Berkowitz Н.М. [82] решалась задача об определении осесимметричных напряжений, возникающих в результате удара полубесконечной упругой цилиндрической оболочки, двигающейся с осевой скоростью, по направлению к жесткой преграде. С помощью безмоментной теории были получены не только главный член асимптотического разложения для больших значений времени, но и дополнительные члены. Это позволило оценить точность решения для точек, находящихся далеко позади фронта волны. Кроме того, было получено решение, имеющее силу вблизи фронта волны.

В работе [56] А. П. Малышева и В. И. Паничкина рассматривались одномерные волновые процессы в оболочках вращения, возникающие при действии быстроизменяющихся нагрузок. В работе показано, что величина разрыва на фронте продольной волны остается постоянной. Полученные результаты демонстрируют разделение НДС оболочки на безмоментное состояние и краевой эффект.

При решении задач нестационарной динамики большинство используемых методов основывается на применении интегральных преобразований Лапласа и Фурье. Для получения решения в начальные промежутки времени Н. Д. Векслер в работах [8,9] использовал метод прифронтовой асимптотики, основанный на разложении изображений по Лапласу в ряды по отрицательным степеням параметра преобразования. С удалением от фронтов использовались методы перевала и стационарной фазы. Разнообразные методы обращения решения двумерных задач для стержней и пластин и безмоментных задач для цилиндрических оболочек, основанные на аналитических свойствах интегральных преобразований, освящены в [74,75].

Значительный вклад в развитие асимптотических методов решения нестационарных задач теории оболочек внесли работы Л. Ю. Коссовича [4146]. Исследования нестационарного волнового НДС оболочек вращения проводились в [42] с использованием понятия показателя изменяемости, введенного А. Л. Гольденвейзером. Рассматривался один из важных классов нестационарных задач — класс задач о распространении волн деформаций в оболочках вращения под действием ударных нагрузок, приложенных к торцу оболочки. В работе [42] исследовались осесимметричные волновые процессы в оболочках вращения. Были выделены тангенциальная составляющая, описывающая распространение волны растяжения-сжатия. И нетангенциальная, описывающая распространение изгибной и сдвиговой волн. В случае неосесимметричных воздействий в [43] проведено исследование областей согласования интегралов теории Кирхгофа-Лява, описывающих НДС в областях оболочки, примыкающих к торцу и быстроизменяющихся составляющих, описывающи напряженное состояние в прифронтовых областях. Установление областей согласования динамического краевого эффекта, описываемого теорией Киргофа-Лява, и антисимметричной составляющей быстроизменяющегося прифронтового поля было проведено также в [44]. Здесь же определены зоны действия приближенных теорий и расположения областей согласования этих теорий.

В монографии [45] были разработаны асимптотические подходы к нестационарным задачам для оболочки вращения с меридианом произвольной формы, основанные на расчленении НДС на безмоментную, моментную составляющие, и (здесь понимались быстроизменяющиеся составляющие) погранслой с различными показателями изменяемости. Были установлены области действия в фазовой плоскости всех составляющих и доказано сращивание безмоментной составляющей и плоского обратносимметричного погранслоя. Была показана связь различного типа асимптотик двумерных решений со свойствами областей их применимости. На основании такого анализа был разработан подход к получению таких асимптотик без использования интегральных преобразований.

В дальнейшем концепция расчленения нестационарного НДС, разработанная Л. Ю. Коссовичем в [45], получила свое дальнейшее развитие. Так, с показателями изменяемости и динамичности равными 2, были выявлены в работах [37,38,95,96] погранслой в малых окрестностях фронтов волн расширения и сдвига, построены их асимптотически оптимальные уравнения и доказано сращивание указанных погранслоёв с решениями для квазиплоской задачи. В работах [39,46] исследовано нестационарное НДС при ударно приложенных гармонических торцевых воздействиях на оболочки вращения и выявлено расположение зон перехода нестационарной части решения в гармонические колебания. В работах [7,81] рассмотрено нестационарное НДС в оболочках из вязкоупругих материалов.

Введение

нового параметра — показателя интенсивности времени релаксации позволило описать новую классификацию нестационарного НДС конструкций в случае вязкоупругих материалов с дифференциальной формой уравнений состояния. Выведены асимптотически оптимальные уравнения всех составляющих.

Несмотря на значительное количество работ, посвященных применению асимптотических методов к задачам динамики оболочек в двумерной постановке, асимптотический вывод двумерных уравнений из трехмерных уравнений теории упругости ранее не был осуществлен. В работе Ю. Д. Каплунова, И. В. Кирилловой, Л. Ю. Коссовича [35] проведено асимптотическое интегрирование трехмерных динамических уравнений теории упругости для случая тонких оболочек. Были получены предельные двумерные системы.

Исследования Ю. Д. Каплунова [31−34,88] внесли существенный вклад в изучение нестационарных волновых процессов.

В [31] к построению двумерных уравнений теории оболочек, описывающих высокочастотные НДС малой изменяемости, применялся метод асимптотического интегрирования трехмерных динамических уравнений теории упругости. Были установлены области применимости и погрешность полученных уравнений. Исследования высокочастотных колебаний оболочки вращения проводились в [32].

В [33] Ю. Д. Каплуновым и Е. В. Нольде проводился асимптотический анализ трехмерных динамических уравнений теории упругости для случая изгиба пластин. В отличие от большинства работ, посвященных асимптотическому построению двумерной динамической теории пластин, два безразмерных асимптотических параметра (показатель изменяемости и показатель динамичности) полагались независимыми.

В работах [23,88] обсуждаются линейные TP теории (по имени С. П. Тимошенко и Е. Рейснера) пластин и оболочек, т. е. теории, учитывающие деформацию поперечного сдвига и инерцию вращения. Ставится вопрос об их построении асимптотическим методом и о вытекающих из этого оценках погрешностей для задач статики и динамики. Предлагается метод расширения применимости TP теорий для динамических задач.

Принципиально новым является результат, приведенный в [94]. Было показано, что для определения НДС оболочки вращения в окрестности квазифронта вместо общих уравнений коротковолновой высокочастотной составляющей можно использовать уравнения теории Кирхгофа-Лява с учетом оператора приведенной инерции. Уточнение этих уравнений позволяет расширить область применимости классической двумерной теории.

Значительное время сфера применения асимптотических методов в нестационарной динамике тонких упругих тел ограничивалась случаями нестационарных волн в пластине, круговой цилиндрической оболочке и оболочках вращения. Результаты исследований в области асимптотической теории тонких упругих тел обобщены Ю. Д. Каплуновым, Л. Ю. Коссовичем, Е. В. Нольде в монографии [93]. На основе трехмерных уравнений теории упругости получены асимптотически оптимальные уравнения низкочастотных, высокочастотных и длинноволновых высокочастотных приближений, позволяющие в совокупности описать как стационарные, так и нестационарные динамические процессы. Разработаны двумерные теории высшего порядка пластин и оболочек, рассмотрены задачи колебания оболочек вращения, тонких тел в среде, излучения тонкими телами. Выведены асимптотически оптимальные уравнения динамических погранслоев в окрестностях фронтов волн расширения и сдвига, в окрестности квазифронта.

Окончательная форма расчленения нестационарного НДС оболочек вращения описана в [47] при торцевых ударных воздействиях видов LT и LM.

Использование предложенных методов для новых классов задач было проведено в следующих направлениях:

— изучения установления стационарных процессов (Березин В.Л., Парфенова Я. А., Никонов А.В.),.

— анализ волн в случаях новых классов ударных воздействий (Каплунов Ю.Д., Мухомодияров P.P., Шевцова Ю.В.),.

— исследования отраженных и прошедших волн в подкрепленных оболочек вращения (Копнина А.Ю.),.

— исследование влияния трансверсально изотропного упругого материала на нестационарные волны (Шевцова Ю.В., Мухомодияров P.P.,),.

— исследования влияния вязкоупругости материалов на нестационарные волны (Бажанова Н.С., Анофрикова Н. С., Сухоловская М.С.), что подтверждает универсальность разработанных методов.

Данная диссертационная работа посвящена исследованию нестационарного волнового НДС при ударных воздействиях нормального типа на лицевые поверхности и на торец пластин. В работе развивается асимптотический подход к решению задач о переходных волновых процессах в пластинах.

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

Исследования выполнены на базе линейной теории упругости. В первом пункте первой главы приведена постановка трехмерных задач теории упругости для бесконечной и полу бесконечной пластины. Далее в главе приведена постановка задачи по двумерной теории, а именно, по теории типа Тимошенко. Поставлены граничные условия на торце и на лицевых поверхностях пластин, моделирующие ударную нагрузку нормального типа. Уравнения записаны в декартовой системе координат.

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

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

В третьей главе исследуется задача распространения волн в полубесконечной пластине, к торцу которой приложена ударная нагрузка типа NW, постоянно распределенная по толщине и не изменяющаяся со временем. В первом параграфе поставленная задача изучается на основе трехмерной теории упругости. Формальное решение уравнения находится с помощью двукратных интегральных преобразований: Лапласа — по времени и синуси косинуспреобразования Фурье — по продольной координате. Обращение двукратного интегрального преобразования произведено в следующем порядке: сначала обращение преобразования Фурье по теории вычетов, а затем преобразование Лапласа. Анализируется решение на примере перерезывающего усилия А^, напряжения оп и перемещения и3. Во втором параграфе данная задача изучается по теории типа Тимошенко. Решение уравнения проводится с помощью интегрального преобразования Лапласа по времени. Обращение преобразования Лапласа осуществляется с помощью метода перевала. При малых значениях времени обращение преобразования проводится с помощью разложения по отрицательным степеням корня из параметра преобразования, а при больших значениях времени — метода прифронтовой асимптотики.

Четвертая глава посвящена построению уравнения погранслоя в окрестности фронта поверхностной волны Релея символическим методом А. И. Лурье. В первом пункте приводится вывод уравнения погранслоя для упругого полупространства при нормальном воздействии [36]. Во втором параграфе выводится приближенные уравнений погранслоя в окрестности фронта поверхностных волн Релея для пластины, при сосредоточенном ударном нормальном воздействии на лицевые поверхности. Применение интегральных преобразований Лапласа и Фурье к этой задаче дает решение, полностью совпадающее с решением, полученным во второй главе по трехмерной теории упругости.

В заключении диссертации сформулированы основные выводы.

Основные результаты диссертации опубликованы в работах [48,49,50,52,53].

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

Научная новизна диссертации: впервые разработан асимптотический подход к решению нестационарных задач для пластин при нормальных воздействиях, основанный на расчленении НДС, выделены границы применимости различных составляющих решения. Кроме того, впервые асимптотическими методами проанализировано и приближенно аналитически определено решение в окрестности условного фронта поверхностных волн Релея. Также построены уравнения погранслоя для поля Релея в пластине символическим методом А. И. Лурье на основе анализа трехмерных уравнений теории упругости.

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

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

Апробация работы. Основные результаты исследований, выполненных в диссертации, доложены: на III Всероссийской конференции по теории упругости с международным участием (г. Ростов-на-Дону, 2003 г.), научных семинарах кафедры математической теории упругости и биомеханики Саратовского государственного университета.

ЗАКЛЮЧЕНИЕ

.

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

1. На основе общей схемы расчленения НДС пластин на составляющие с различными показателями изменяемости разработана методика асимптотического исследования нестационарного НДС при ударных лицевых и торцевых воздействиях нормального типа.

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

3. Предложен новый метод построения решения в окрестности условного фронта поверхностных волн Релея методом стационарной фазы.

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