Содержание к диссертации
Введение
1 Разрешающие уравнения и вывод дисперсионного уравнения 11
1.1 Уравнения состояния 12
1.2 Уравнения движения 15
1.3 Распространение гармонических волн 15
1.4 Дисперсия гармонических волн 17
1.4.1 Изгибные волны 19
1.4.2 Волны расширения 21
1.5 Численный анализ дисперсионных уравнений 22
2 Анализ дисперсионного уравнения 36
2.1 Длинноволновые низкочастотные приближения 37
2.1.1 Изгибные волны 37
2.1.2 Волны растяжения 39
2.2 Длинноволновые высокочастотные приближения 41
2.2.1 Изгибные волны 43
2.2.2 Волны расширения 46
2.3 Коротковолновые высокочастотные приближения 48
2.3.1 Случай 1(a) 48
2.3.2 Случай 1(b) 51
2.3.3 Случай 2 54
2.4 Коротковолновый предел фундаментальной моды 56
3 Длинноволновые низкочастотные модели 59
3.1 Симметричные движения 60
3.1.1 Относительные порядки перемещений 60
3.1.2 Асимптотически приближенные уравнения 62
3.1.3 Модельная задача: торцевая нагрузка 71
3.2 Антисимметричные движения 74
3.2.1 Относительные порядки перемещений 75
3.2.2 Асимптотически приближенные уравнения 76
3.2.3 Модельная задача: торцевая нагрузка 86
4 Длинноволновые высокочастотные модели 92
4.1 Изгибные волны 93
4.1.1 Асимптотическое интегрирование:первое семейство частот запирания 93
4.1.2 Второе семейство частот запирания 102
4.2 Волны расширения 110
4.2.1 Асимптотическое интегрирование:первое семейство частот запирания 110
4.2.2 Второе семейство частот запирания 115
4.3 Иллюстративные численные результаты 115
5 Длинноволновые модели для идеализированных волокнистых упругих конструкций 121
5.1 Разрешающие уравнения и дисперсионное соотношение 122
5.2 Численный анализ 126
5.3 Анализ дисперсионного уравнения 127
5.3.1 Длинноволновые высокочастотные приближения 128
5.3.2 Длинноволновые низкочастотные приближения 129
5.3.3 Коротковолновые высокочастотные приближения 132
5.3.4 Коротковолновый предел фундаментальной моды 134
5.4 Длинноволновые асимптотические модели 135
5.4.1 Длинноволновые высокочастотные движения 135
5.4.2 Длинноволновые низкочастотные движения 140
Литература 145
- Уравнения движения
- Длинноволновые высокочастотные приближения
- Антисимметричные движения
- Волны расширения
Введение к работе
Существует множество новых технологий, например в аэрокосмической промышленности или при проектировании сложных механических конструкций, которые требуют наличия материалов с конкретными динамическими характеристиками, некоторые из которых не встречаются в обычных, монолитных системах. Композитные материалы являются одним из решений этой проблемы. Такие материалы обычно являются разнородной смесью двух или более составляющих, которые объединены вместе. Среди составляющих композитов, одна фаза, так называемая армирующая фаза или арматура, часто в форме волокон или листов, внедрена в другой материал, матрицу. Арматура - это твердые, прочные материалы, состоящие обычно из углерода, стекла, керамики или металла. В аэрокосмической индустрии широкое применение находят матрицы на основе эпоксидных или резиноподобных материалов, см. Daniel & Ishai (1994) и Kelly et al. (1993). Композитные материалы, включая ламинированные и волокнистые, играют значительную роль во многих реальных научных и инженерных задачах, в первую очередь благодаря их уникальным физическим и механическим характеристикам, таким, как особая прочность, стойкость к разрушению и износостойкость. В настоящее время исследование и разработка композитов принадлежит к одному из наиболее интенсивно изучаемому разделу в науке о материалах и инженерной науке. В диссертации мы будем исследовать композиты, армированные волокнами.
Рассматривая механические свойства армированных композитов, мы сосредоточимся только на макроскопическом поведении композитов. Таким образом, мы не будем рассматривать задачи, включающие локальные взаимодей-
Введение ствия между волокнами и матрицей. С этой точки зрения принято использовать непрерывные теории при описании механических или динамических характеристик композита, как единого целого. В рамках данных теорий предполагается, что волокна являются скорее неотъемлемым свойством материала, чем некоторой формой включения. Идея разработки непрерывной теории для описания армированных материалов, по-видимому, берет начало в серии статей Адкинса, см. например Adkins (1958) и Green & Adkins (1960) в теории больших деформаций резиноподобных материалов, армированных нерастяжимыми струнами.
Хотя мы ограничиваемся непрерывными теориями и макроскопическими моделями, свойства рассматриваемых нами композитных материалов происходят от свойств и геометрического расположения их составляющих. Мы полагаем, что волокна представляют собой одно семейство параллельных прямых линий, представленных в композите в одном направлении. Таким образом, в макроскопическом масштабе композит будет вести себя как локально трансвер-сально изотропный материал с осью трансверсальной изотропии совпадающей с направлением волокон. Трансверсально изотропное линейно-упругое твердое тело характеризуется пятью константами, см. например Green (1982).
Однонаправленные упорядоченные волокна могут иметь большую прочность если нагрузка направлена параллельно волокнам. Такие композиты имеют сильно выраженное анизотропное поведение и их модули зависят от направления нагрузки и волокон. Для того, чтобы решить эту проблему, волокна зачастую упорядочивают в слои, которые располагают определенным образом, получая так-называемые ламинаты. Теоретическое и экспериментальное исследование распространения волн и колебаний в таких структурах является областью, к которой проявляется огромный интерес в последние несколько лет, благодаря большому количеству промышленных приложений таких структур. Подробный список ссылок может быть найден в обзорных статьях Chimenti (1997) and Lowe (1995). Однако, постановка трехмерных динамических задач и построение асимптотически приближенных теорий применительно к многослойным структурам является достаточно сложной задачей, когда возможная
Введение физическая интерпретация часто теряется в алгебраических деталях. Следовательно, представляется оправданным использование некоторых идеализированных теорий для выявления особенностей таких материалов и обеспечения улучшенного физического понимания их поведения.
Для некоторых классов композитов, например резины, армированной углеродными или нейлоновыми волокнами, их свойства могут быть идеализированы при предположении, что композит является несжимаемым и (или) нерастяжимым в направлении волокон. Механический эффект введения данных ограничений обсуждается в контексте конечной упругости в Green & Adkins (1960) и Truesdell & Noll (1965). Заметим, что такие идеализированные теории следует применять с осторожностью, поскольку они могут быть менее реалистичны в некоторых обстоятельствах, например по отношению к теориям бесконечно малых деформаций. Однако, потенциальное упрощение достигаемое при использовании данных ограничений обычно оправдывает использование таких приближенных упрощенных теорий.
В диссертации мы будем накладывать ограничение несжимаемости в обычной трансверсально изотропной теории, получая уравнения состояния, имеющие только три константы материала. Такая непрерывная модель была детально исследована и ее упрощения использовались в Musgrave (1970) и позже, в серии статей Чадвика, см. Chadwick (1993), Chadwick (1994а) и Chadwick (1994b) .Для многих реальных композитов, волокна являются прочными и жесткими по отношению к матрице. Таким образом, сопротивление композита к растяжению в направлении волокон может быть на несколько порядков больше по величине, чем их сопротивление к поперечному растяжению или к осевому или поперечному сдвигу. Это свойство позволяет рассматривать композит в качестве сильно анизотропного твердого тела, с большим модулем Юнга в направлении волокон по отношению к модулю Юнга в других направлениях. В предельном случае, это может быть идеализировано с помощью предположения нерастяжимости в направлении волокон. Подробный перечень механических свойств идеализированных материалов и точные ограничения, накладываемые
Введение на непрерывную модель, могут быть найдены в Spencer (1972).
Несколько статей было посвящено сравнению результатов полученных с помощью идеализированной теории для несжимаемого материала с нерастяжимыми волокнами с результатами получаемыми из линейной анизотропной теории упругости (Lekhnitskii (1963)). Результаты показывают, что методы теории по-гранслоя или теории сильного возмущения, могут быть применены для анализа задач, включающих сильно анизотропные материалы. Rogers & Pipkin (1971) успешно применяли такие теории к задаче плоского изгиба консольной балки. Этот анализ показывает, что идеализированная теория дает корректное решение первого порядка во внутренней области балки. Однако, на верхней и нижней гранях балки идеализированная теория предсказывает большие концентрации напряжений, которые связаны с погранслоями напряжения сильного растяжения или сжатия.
Одной из основных характеристик некоторых типов сред является дисперсия волн, указывающая, что фазовая скорость волны (и частота) являются функциями волнового числа. Хотя дисперсия и не является атрибутом распространения гармонических упругих волн в бесконечной области, она является свойством упругих слоистых структур. Дисперсионное уравнение, или частотное уравнение, получается подстановкой решения в форме распространяющейся волны в уравнения движения и граничные условия и одновременного их удовлетворения. Бесконечное число решений дисперсионного уравнения для каждого волнового числа описывает частотный спектр бесконечного количества мод. Это уравнение обеспечивает нас фундаментальными свойствами распространения волн. Дисперсионная природа распространения волн в бесконечной упругой плите была впервые обнаружена Lord Rayleigh (1888-1889) и Lamb (1889-1890), где были выведены частотные уравнения для волн сжатия и изгибных волн плоской деформации. Природа низких (фундаментальных) мод передачи волн и физика волн которую они формируют были позднее установлены Lamb (1917) и Davies (1948). Более поздние работы Holden (1951), Tolstoy к Usdin (1957), Mindlin (I960) и другие открыли природу более высоких действитель-
Введение ных мод, называемых гармониками. Позднее эти исследования были расширены на более сложные случаи, включая анализ трансверсально изотропного упругого слоя в Green (1982), многослойных армированных волокнами композитов в Achenbach &; Herrmann (1968) и Baylis к. Green (1986) и структур периодически слоистых моноклиническими слоями в Nayfeh (1991). Из-за сложности основных дисперсионных уравнений, большинство результатов в этих статьях были получены с помощью численного анализа дисперсионных уравнений. В недавних статьях показана ценность и потенциально высокая полезность асимптотического анализа дисперсионных уравнений, особенно по отношению к определению динамической реакции на высокоскоростные ударные воздействия или построению приближенных теорий, см. например Kossovitch к, Rogerson (1999), Kaplunov et al. (2000). В частности, подробное исследование природы дисперсионного уравнения требуется до численного обращения интегральных преобразований при определении реакции на ударные воздействия, см. например Rogerson (1992) и Green (1993).
Аналитические и численные методы для анализа дисперсионных уравнений были также обобщены для предварительно напряженной несжимаемой плиты, см. Ogden & Roxburgh (1993) и Rogerson & Fu (1995), для некоторых типов многослойных сред, см. Rogerson & Sandiford (1996) и Rogerson & Sandiford (2002), для почти несжимаемой среды, см. Sandiford & Rogerson (2000) и для трансверсально изотропных плит, см. Kossovitch et al. (2002а). Однако, алгебраическая сложность дисперсионных уравнений для некоторых структур делает привлекательными приближенные теории, где упрощенные модели позволяют предсказать динамическое поведение для некоторых режимов распространения волн. Существуют хорошо известные структурные теории более низкого порядка, в основном основанные на физических гипотезах, например теории плит Кирхгофа и Миндлина, теория оболочек Кирхгофа-Лява и улучшенные теории Тимошенко-Рейсснера, см. например Lord Rayleigh (1877), Love (1927), Timoshenko (1938),Mindlin (1955), Reissner (1964). В этих теориях делаются попытки совместить обычную физику точной теории (высокой размерности)
Введение с упрощениями элементарных теорий. Диапазон применимости этих теорий обычно ограничен, например низкочастотная и длинноволновая область фундаментальной моды в частотном спектре. Корректно предсказывая распространение волн в ведущем порядке, эти теории не всегда полностью согласуются с точной трехмерной теорией.
Позже были многочисленные попытки доказать корректность приближенных теорий без каких-либо физических предположений. Для этого применялся метод асимптотического интегрирования точных трехмерных уравнений движения, см. например Goldenveiser (1966). Совсем недавно этот метод был расширен для динамических задач, также позволяя улучшить приближенные теории, учитывая длинноволновый и коротковолновый вклад мод, см. Achenbach (1969) в контексте стационарных колебаний изотропного упругого слоя и Berdichevskii (1977), см. также Kaplunov (1990) для нестационарных задач. Полный и подробный перечень асимптотических методов, требуемых для определения динамической реакции тонкостенных упругих конструкций может быть найден в Kaplunov et al. (1998). Недавно были разработаны некоторые предварительные приближенные теории для трансверсально изотропной и несжимаемой предварительно нагруженной плиты, см. Kaplunov et al. (2000), Pichugin &; Rogerson (2002).
Целью данной диссертации является вывод серии длинноволновых моделей для выявления динамической реакции несжимаемых трансверсально изотропных плит. Диссертация организована следующим образом. Вывод основных разрешающих уравнений и анализ дисперсии гармонических волн в упругой плите из несжимаемого трансверсально изотропного материала приводится в главе 1. Рассматривается несжимаемое трансверсально изотропное упругое тело, находящееся в естественном состоянии равновесия при отсутствии деформаций. Относительно этого состояния, бесконечно малые, зависящие от времени движения, задаются в форме распространяющихся волн. Сначала выводятся линеаризованные уравнения движения и линеаризованное условие несжимаемости. Характеристическое уравнение выводится из условия существования
Введение нетривиальных решений для плоской волны из разрешающих уравнений и условия несжимаемости. Мы рассматриваем распространение гармонических волн в плите со свободными лицевыми поверхностями, имеющей конечную толщину и бесконечную в продольном направлении. Решения характеристического уравнения могут быть использованы для представления соответствующих волновых решений в перемещениях и напряжений как суперпозиции линейно независимых функций. После подстановки этих решений в граничные условия, получаем однородную систему из шести уравнений с шестью неизвестными. Симметрия задачи позволяет разделить систему на симметричные и антисимметричные решения, которые приводят к, так называемым, дисперсионным уравнениям расширения и изгиба. Данные уравнения затем анализируются численно для выявления их принципиальных свойств, а также для того, чтобы выделить различные асимптотические режимы распространения волн.
После проведенного в главе 1 численного анализа, во второй главе мы сосредоточиваемся на асимптотическом анализе дисперсионных уравнений, соответствующих распространению гармонических волн в несжимаемой транс-версально изотропной плите при отсутствии напряжений на лицевых поверхностях. Исследование включает вывод асимптотических разложений для малых и больших значений волнового числа, задающих фазовую скорость (или частоту) как функцию волнового числа для каждой ветки. Сначала выводятся длинноволновые асимптотические разложения для фундаментальных волн изгибной и продольной составляющих. Полученные асимптотические приближения второго порядка согласуются с численным решением дисперсионного уравнения в приемлемом диапазоне волновых чисел. В длинноволновом высокочастотном режиме получены разложения третьего порядка для двух семейств частот запирания для изгибной и продольной составляющих. Проведенное сравнение этих разложений с численным решением показывает совпадение в окрестности частот запирания. Коротковолновые приближения дисперсионного уравнения выводятся для многих возможных комбинаций параметров материала, описывая различную реакцию материала, включая осцилляции гармоник
Введение при kh —> со. Для коротковолнового низкочастотного типа движения показано, что предел при большом значении волнового числа для фундаментальных мод изгибной и продольной составляющих является скоростью соответствующей поверхностной волны Релея. Материал второй главы послужил основой для статьи Kossovitch et al. (2002а).
Вывод моделей более низкой размерности мотивируется стремлением сократить количество существенных параметров таким образом, чтобы было возможно выявлять качественные особенности процесса. В главе 3 мы начинаем с двух длинноволновых низкочастотных моделей, одна для симметричных, а другая для антисимметричных движений. Численный и асимптотический анализ соответствующих дисперсионных уравнений был выполнен в первых двух главах, показывая, что существуют два конечных длинноволновых предела дисперсионного уравнения, связанного с продольной составляющей и один предел в случая изгибного движения. Полученные асимптотические разложения позволяют нам проанализировать асимптотическую структуру распространения волн в соответствующем длинноволновом низкочастотном режиме. Затем определяются относительные порядки перемещений и вводится соответствующее масштабирование пространственных и временных переменных для представления уравнений движения, условия несжимаемости и граничных условий в форме, подходящей для асимптотического интегрирования. Это позволяет получить иерархическую систему краевых задач для различных асимптотических порядков. Решая данную систему для каждого порядка, мы получаем решения для компонент перемещения и давления в терминах прогиба срединной поверхности для изгибных волн и в терминах разрешающих расширений для продольных перемещений. Данные функции могут быть получены как решения соответствующих разрешающих уравнений. Корректность данного подхода демонстрируется тем фактом, что дисперсионное уравнение, связанное с асимптотическим разрешающим уравнением более низкого порядка, совпадает с соответствующим приближением точного дисперсионного уравнения, выведенного во второй главе. Для иллюстрации потенциальных приложений данных длинноволновых низко-
Введение частотных моделей более низкого порядка сформулированы и решены некоторые модельные задачи на ударную нагрузку.
Закончив вывод асимптотически корректных моделей для длинноволнового низкочастотного движения, в главе 4 мы обращаемся к длинноволновым высокочастотным моделям в трансверсально изотропной несжимаемой плите. Асимптотический анализ, проведенный в главе 2 показал, что существует два различных асимптотических случая, связанных с двумя семействами частот запирания. Разложения второго порядка, полученные в главе 2 позволяют определить относительные порядки перемещений и давления для данного типа движения. Также вводится соответствующее масштабирование пространственных координат и времени. Асимптотически приближенные уравнения выводятся для построения асимптотической модели, позволяющей нам представить перемещения и давление в терминах одной функции, называемой длинноволновой амплитудой. Разрешающие уравнения ведущего и уточненного второго порядка получаются из условия существования решения в процессе вывода модели. Корректность полученных двумерных моделей демонстрируется сравнением каждого точного дисперсионного уравнения связанного с разрешающими уравнениями с соответствующим асимптотическим разложением дисперсионного уравнения, связанного с полной трехмерной теорией, полученным в главе 2. В конце данной главы коротко обсуждаются некоторые интересные особенности асимптотических моделей, связанные с возможным существованием отрицательной групповой скорости.
Заключительная пятая глава посвящена построению длинноволновых асимптотических моделей более низкой размерности применительно к нерастяжимому в направлении волокон материалу. Мы начинаем с короткого обзора разрешающих уравнений и вывода дисперсионных уравнений для симметричных и антисимметричных движений в плите, состоящей из данного материала. Затем выполняется численный и асимптотический анализ, позволяющий выявить, что, в отличи от растяжимого случая, длинноволновый предел антисимметричной фундаментальной моды не равен нулю. Это объясняется тем, что нало-
Введение женное условие нерастяжимости препятствует изгибу в классическом смысле. Для полноты, приведены также коротковолновые асимптотические приближения. Затем выводится асимптотическая модель более низкой размерности для длинноволнового движения в окрестности частот запирания. Способом, похожим на описанный в главе 4, выводится иерархическая система, приводящая к системе краевых задач, которая решается способом асимптотического интегрирования. НДС определяется в терминах одной функции, называемой длинноволновой амплитудой, в терминах которой и получены разрешающие уравнения ведущего и уточненного второго порядка. На каждом порядке дисперсионное уравнение, связанное с разрешающим уравнением для длинноволновой амплитуды, совпадает с приближениями точного дисперсионного уравнения. Похожий анализ проводится для длинноволнового низкочастотного движения, связанного с фундаментальной модой. Асимптотическая корректность снова показывается по отношению к длинноволновому низкочастотному движению.
Уравнения движения
В данной диссертации мы имеем дело с гармоническими волнами, в частности с плоскими волнами. Таким образом будем искать решения для уравнений дви Глава 1. Разрешающие уравнения жения и условия несжимаемости в форме гармонической распространяющейся волны где v - фазовая скорость, к - волновое число, (c0,se) = (cos 0, sin в) обозначают проекцию вектора нормали к волне на плоскость Ох\х2 и число q - должно быть определено из уравнений движения и условия несжимаемости. Подставляя (1.22) в уравнения (1.17)-(1.20), получаем однородную систему уравнений относительно Ux, U2, Uz и Р в виде Условие существование нетривиального решения системы (1.23) позволяет заключить, что q удовлетворяет где масштабированная скорость волны v2 = pv2 и Мы будем обозначать два возможных корня квадратного уравнения для q2 (1.24)i, как qf и q2. Для последующего использования в асимптотическом анализе дисперсионного уравнения выпишем следующие стандартные результаты Глава 1. Разрешающие уравнения Для вообще говоря отличных и ненулевых корней характеристического уравнения (1.24), любое решение для щ, и?,, щ или р может быть представлено как суперпозиция шести линейно независимых функцийexp(kqmx3) и ехр(—kqmXz), т Є {1,2,3}. Эти комбинации могут быть выражены в следующем виде в которой Uj и p(m\ j є {1,2,3}, m Є {1,2,3,4,5,6} - двадцать четыре произвольные постоянные. Эти константы не все являются независимыми, т.к. они связаны через уравнения движения (1.17)-(1.19) и условие несжимаемости (1.20).
Эти уравнения могут быть использованы для получения решений для перемещений и давления в терминах только шести независимых констант, а именно Глава 1. Разрешающие уравнения Мы будем рассматривать только трехмерный случай и, следовательно, можем по умолчанию подразумевать, что с ф 0 и s ф 0. Можно легко показать, что что вместе с предыдущим соглашением 0 0 0 ф 0 к Sg ф 0 значит, что обнуление b(qm,v) происходит всегда одновременно с обнулением f{qm,v) и h(qm,v). Данная ситуация связана с вырожденным решением, т.к. граничные условия удовлетворяются автоматически. Решения (1.29) могут быть теперь подставлены в соответствующее подмножество (1.12)-(1.14), позволяя найти соответствующие представления длясгіз, Ст2з и зз- Для свободных граничных условий стіз, 023 и зз все обнуляются при Хз = ±/І, и могут быть использованы для получения следующей системы однородных уравнений Особенностью симметричной плиты с одинаковыми граничными условиями на каждой грани является возможность разделения задачи на антисимметричное и симметричное решения с заменой системы шесть на шесть (1.31) на две системы три на три. Первая из этих систем может быть получена вычитанием (1.31)i из (1.31)г, вычитанием уравнения (1.31)з из (1.31)4 и сложением (1.31)s с (1.31)б, получая в результате получается похожим образом, добавляя (1.31)і к (1.31)2, складывая (1.31)з с (1.31)4 и вычитая уравнение (1.31)s из (1.31)б, получаем Заметим, что в изгибном случае Щт = Щ , т Є {1,2,3}, и, следовательно «і, и2 и р являются антисимметричными функциями относительно Хз и из - симметричная, см. (1.29). Данный тип движения обычно называется из-гибпым или антисимметричным движением. Необходимо отметить, что когда Щт = Щ т система (1.33) удовлетворяется тривиально, следовательно дисперсионное уравнение для антисимметричных (изгибных) волн получается требованием, чтобы система уравнений (1.32) имела нетривиальное решение, и Глава 1. Разрешающие уравнения может быть представлено в следующей форме где бездисперсионный множитель опущен. Для последующего асимптотического анализа потребуется представление для функций G(q, v), Ф(д, v) и Ф(г;) с более прозрачной зависимостью от v и q. Следовательно, на этом этапе запишем следующее представление Глава 1. Разрешающие уравнения Заметим, что компоненты перемещения и давление (1.29), связанные с изгиб-ным движением могут быть записаны в виде Условие подразумевает, что все компоненты пе ремещений и давление являются симметричными функциями по отношению к х3, за исключением и3, которая антисимметрична. Данный тип колебаний обычно называется расширяющим или симметричным. Граничные условия (1.32) удовлетворяются тождественно и нам нужно только рассмотреть однородную систему линейных уравнений (1.33). Эта система имеет нетривиальное решение, когда определитель соответствующей матрицы равен нулю. Таким образом, дисперсионное уравнение, связанное с волнами расширения выражается в виде
Длинноволновые высокочастотные приближения
Предметом нашего интереса теперь будет длинноволновые высокочастотные колебания для которых v —» со при kh — 0. Как было замечено ранее, это лава 2. Анализ дисперсионного уравнения длинноволновое поведение гармоник в окрестности частот запирания. Анализ характеристического уравнения (1.24) предполагает, чтод2 - величина порядка 0(v2)y один из корней (1.24)i, скажем д2, является величиной порядка 0(v2) и второй д2 - порядка 0(1). Раскладывая эти корни в степенной ряд по v2 и подставляя их в характеристическое уравнение (1.24), может быть выведено следующее приближение Две скорости волны, связанные с q\ являются величинами порядка 0(1), и, следовательно, не соответствуют высокочастотному колебанию. Заметим, что отсутствие третьей волны, связанной с #з является следствием наложения условия несжимаемости, которое делает невозможным распространение продольной волны.
В обычном материале существует третья возможная большая скорость, связанная с дз см- например Kaplunov et al. (2000) для случая сжимаемой транс-версально изотропной плиты. Подставляя разложения (2.14) в дисперсионное уравнение (1.34), получаем следующее приближение Так как все коэффициенты Fx \ F2 , F$ , Fx , F2 , F% являются величинами порядка 0(1), а дз - 0(1), и, следовательно, tanh(g3 ) является величиной порядка 0(г -1), дисперсионное уравнение (2.16) может быть асимптотически Глава 2. Анализ дисперсионного уравнения сбалансировано только если T\(h) 0(v2) или T2(h) 0(v2). Эти два случая связаны с двумя возможными семействами частот запирания. В первом случае T\(h) 0(v2), следовательно ведущий порядок аргумента гиперболического тангенса равен і (п + ) 7г. Таким образом, мы можем разложить этот аргумент по степеням малого параметра kh, а именно в котором параметры ф2 и ф{ являются неизвестными величинами вообще говоря порядка 0(1). Соответствующее разложение для Ті(h) задается Введем параметры Л{ = (n + \)ircz, п = 1,2,3..., которые являются соответствующими частотами запирания. Теперь мы можем использовать (2.14), чтобы выписать
Данные приближения позволяют нам аппроксимировать все необходимые функции в (2.16), таким образом Подставляя разложения (2.18)-(2.24) назад в приближение дисперсионного уравнения второго порядка (2.16), мы получаем выражения для ф{ и ф{ъ форме Глава 2. Анализ дисперсионного уравнения Теперь мы можем использовать (2.15) для получения приближения третьего порядка для масштабированной частоты в форме изображен результат сравнения асимптотических разложений (2.27) и (2.28) с численным решением в окрестности частот запирания для первых семи гармоник. В частности, при увеличении номера гармоники п совпадение становится частично лучше, с неразличимыми асимптотическим и численным решениями для п = 7. Рис. 2.3: Сравнение численных и асимптотических решений дисперсионного уравнения для первых ненулевых частот запирания для изгибных волн, ЦЕ Ць №Т 1: 2 : 3, в = 75, corresponding to щ = 1.3228, ц3 = 1.4377, 4 = 1.2433. Подставляя приближение (2.14) в дисперсионное уравнение (1.40), мы получаем ведущий порядок приближения дисперсионного уравнения, а именно Далее подразумевается, что функции материальных параметров Лі, А2 являются в общем случае величинами порядкаО(І). Так как q3 — 1+0(гГ2), очевидно, что T3(h) = 0(v l), указывая на то, что Tx{h) = 0(v 4) или T2(h) = 0{v 4). Каждое из эти асимптотических соотношений соответствует гармоникам, связанным с одним из двух возможных семейств частот запирания. Сначала рассмотрим случай, в котором Ті (/i) = 0(v 4) и, следовательно, асимптотический порядок Ті (h) подразумевает следующие формы разложений
Антисимметричные движения
Вывод асимптотически приближенных уравнений включает в себя масштабирование пространственных и временной переменных. Масштабирование пространственных переменных имеет целью сбалансировать полутолщину плиты h и нормальную координату с длиной типичной волны I с координатами в плоскости плиты, показывая, что подходящим выбором будет где i, 2 и
С являются новыми безразмерными пространственными координатами. Для изгибного движения типичная волна рассматриваемого типа распространяется с масштабированной скоростью rj y/FW, где F = F± cj+F Cgsl+ F3 sj, см. разложение (2.8), и, таким образом, проходит расстояние длины ти пичной волны I за время l\fp/{r]\jFW) и соответсвующий масштаб времени может быть выбран в виде с г обозначающим безразмерное время. В соответствии с распределением относительных порядков перемещений в (3.69), компоненты перемещения и давления должны иметь следующую асимптотическую структуру с u\, «2, p являющимися нечетными функциями и Из четной по . Все величины с верхним индексом обозначают величины одного асимптотического порядка. Теперь мы можем переписать уравнения движения (1.17)-(1.19) в терминах безразмерных переменных Глава 3. Длинноволновые низкочастотные модели в которых запятая в нижнем индексе обозначает дифференцирование по отношению к указанной масштабированной (пространство или время) переменной. Уравнения (3.73)-(3.75) должны быть решены вместе с соответствующим образом масштабированным условием несжимаемости и при соответствующих граничных условиях свободных от напряжений лицевых поверхностях Вся информация, касающаяся длинноволновых низкочастотных вкладов различных слагаемых в уравнения (3.73)-(3.77) теперь выражается непосредственно через наличие 772. Это позволяет нам организовать схему для получения разрешающих уравнений меньшего порядка для этого типа движения. Будем искать решения в виде степенных рядов
Подставляя решение (3.78) в задачу (3.73)-(3.77) мы получаем иерархические системы различного порядка. Ведущий порядок Для ведущего порядка уравнения движения (3.73)-(3.75) принимают вид Глава 3. Длинноволновые низкочастотные модели которые должны быть решены вместе с ведущим порядком условия несжимаемости и соответствующим граничным условием Общее решение краевой задачи (3.81), (3.83)з может быть представлено как в котором Щт = и;т(Ъ,&,т). Напомним, что мы используем двойной верхний индекс для обозначения функций независимых от , с первым индексом, обозначающим порядок приближения и вторым, обозначающим степень любого возможного множителя С"1. После подстановки решения (3.84) в краевую задачу (3.79), (3.83\ можем (о) получить решение ДЛЯ Ui в виде Данная процедура может быть повторена для задачи (3.80), (3.83)г, получая Таким образом, мы получили ведущий порядок решений для компонент перемещений в терминах функции U3 , или, более точно, линейной комбинации этой функции и ее производных. Мы будем называть ведущий порядок нормального перемещения ведущим порядком изгиба срединной поверхности. Знание этой функции определяет длинноволновое низкочастотное движение в ведущем порядке. Однако, разрешающее уравнение для этой функции не может быть получено без введения в рассмотрение более высоких порядков уравнений.
Волны расширения
Закончив вывод модели для изгибных волн, обратимся к модели для длинноволновых высокочастотных движений типа расширения-сжатия. Снова можем наблюдать два различных асимптотических случая, соответствующих двум семействам частот запирания. Мы начинаем с сравнения относительных порядков перемещений и давления. Это дает нам необходимую информацию об асимптотической структуре данного вида перемещений. В случае волн расширения (симметричном) относительные порядки перемещений для обоих случаев снова задаются (4.1) и (4.2), соответственно. Для того, чтобы изучить движение в окрестности частот запирания, полезно ввести масштабированные переменные пространства и времени i, 2, С и т определенные следующим образом Cm где m = 3 для первого набора частот запирания и m = 2 для второго. Асимптотическая структура перемещений и давления является такой же, как для изгибного типа движения, а именно Глава 4. Длинноволновые высокочастотные модели в которой верхний индекс обозначает безразмерные величины одного асимптотического порядка. Система уравнений движения (1.17)-(1.19) теперь может быть переписана в терминах безразмерных переменных Уравнения движения (4.76)-(4.78) должны быть решены вместе с соответствующим образом масштабированным условием несжимаемости и при соответствующих граничных условиях Разложения в степенной ряд в виде (4.12) позволяют получить иерархические системы обыкновенных дифференциальных уравнений. Результаты для каждого порядка представлены ниже. Ведущий порядок Уравнения движения ведущего порядка должны быть решены вместе с условием несжимаемости ведущего порядка Глава 4. Длинноволновые высокочастотные модели при граничных условиях ведущего порядка Следуя схеме, похожей на использованную для антисимметричного случая, получаем следующее решение системы (4.81)-(4.85)
Таким образом, решения ведущего порядка для компонент перемещения и , г Є {1,2,3} и давления р получены в терминах функции щ; 0 — и{;(fi, 2, т) и ее производных. Заметим, что щ; определяет длинноволновое высокочастотное движение в ведущем порядке длинноволновую амплитуду. Эта функция не может быть определена без рассмотрения более высоких порядков. Второй порядок На втором порядке, уравнения движения принимают вид Глава 4. Длинноволновые высокочастотные модели которые должны быть решены вместе с условием несжимаемости второго порядка и граничными условиями Первое уравнение движения (4.87), рассмотренное вместе с решениями ведущего порядка (4.86) с применением граничного условия (4.91)i, показывает, что имеет место при условии с с и д определенными после разложения (2.33). Т.к. компонента перемещения Из является линейной функцией «і и ее производных, см. (4.86), уравнение (4.93) также удовлетворяется. Используя этот факт мы можем заменить слагаемое, умножаемое на г) 2 в третьем уравнении движения (4.89) и получить решения для «з и Р 2\ а именно Глава 4. Длинноволновые высокочастотные модели Общее решение для щ также получено, однако, т.к. оно не требуется для вывода уточненного разрешающего уравнения для длинноволновой амплитуды, мы его не приводим. Третий порядок Для уточнения разрешающего уравнения для длинноволновой амплитуды, нам необходимо рассмотреть только третье уравнение движения третьего порядка и соответствующее граничное условие Как и в задаче второго порядка, разрешающее уравнение для щ} и ttf теперь может быть получено как необходимое условие для разложения «і из задачи (4.95)-(4.96), принимая вид Уравнение, аналогичное (4.38), задается следующим образом Это уравнение может быть выражено в терминах исходных переменных, а именно Дисперсионное соотношение, связанное с данным уравнением соответствует точному дисперсионному уравнению (1.40) в том смысле, что все три порядка разложения для частоты (2.33) совпадают полностью (первые два порядка в случае разрешающего уравнения ведущего порядка (4.93)). Глава 4.
Длинноволновые высокочастотные модели 4.2.2 Второе семейство частот запирания Анализ асимптотического поведения плиты в окрестности второго семейства частот сдвигового резонанса очень похож на рассмотренный выше случай. Соответствующий вывод приводит к получению разрешающих уравнений ведущего и второго порядка для длинноволновой амплитуды, а именно Уравнение, аналогичное (4.72), задается с S\c И E\S определенными после разложения (2.34). Асимптотическая корректность проверяется сравнением с разложением (2.34).