Использование квазистатического подхода в динамике стержневых систем тема диссертации и автореферата по ВАК РФ 01.02.04, кандидат физико-математических наук Синильщикова, Галина Александровна

  • Синильщикова, Галина Александровна
  • кандидат физико-математических науккандидат физико-математических наук
  • 2008, Санкт-Петербург
  • Специальность ВАК РФ01.02.04
  • Количество страниц 86
Синильщикова, Галина Александровна. Использование квазистатического подхода в динамике стержневых систем: дис. кандидат физико-математических наук: 01.02.04 - Механика деформируемого твердого тела. Санкт-Петербург. 2008. 86 с.

Оглавление диссертации кандидат физико-математических наук Синильщикова, Галина Александровна

стр:

ВВЕДЕНИЕ 2

1 . ИСПОЛЬЗОВАНИЕ РЕАКЦИЙ СВЯЗЕЙ КАК ОБОБЩЕННЫХ

КООРДИНАТ ПРИ МОДЕЛИРОВАНИИ КОЛЕБАНИЙ

СТЕРЖНЕВЫХ СИСТЕМ 10

1.1 Применение квазистатического подхода к задаче об изгибных. колебаниях консоли 11

1.2 Результаты расчетов 17

2. ИСПОЛЬЗОВАНИЕ ВНУТРЕННИХ СИЛОВЫХ ФАКТОРОВ КАК

ОБОБЩЕННЫХ КООРДИНАТ ПРИ МОДЕЛИРОВАНИИ!

КОЛЕБАНИЙ СТЕРЖНЕВЫХ СИСТЕМ^' 18

2.1 Об использовании внутренних силовых факторов как обобщенных координат при моделировании колебаний стержневых систем 18

2.2. Задача о продольных колебаниях свободного стержня 19

2.3 Применение квазистатического подхода к задаче об^ изгибных колебаниях свободного стержня, симметричного относительно серед ины -■; . 25

2.4 Применение квазистатического подхода к задаче об изгибных колебаниях несимметричного свободного стержням1 - 28

2.5. Задача об изгибных колебаниях свободного стержня с сосредоточенными массами 31

3. РАЗВИТИЕ ТРЕЩИНЫ В ТОНКОМ БРУСЕ ПРИ

ИМПУЛЬСНОМ НАГРУЖЕНИИ. 34

3:1. Постановка задачи 34

3.2. Определение прогиба консоли в квазистатическом приближении 36

3.3. Представление прогиба консоли и определение кинетической и потенциальной энергии в72-ом приближении 40

3.4. Определение обобщенных сил, соответствующих введенным координатам 45

3.5; Уравнения Лагранжа II рода и условия начала и прекращения развития трещиньг 47

3.6. Переход к,безразмерным переменным^ 53

3.7 Алгоритм решения ^ 57

3.8 Определение характерных значений импульса I и приближенной зависимости атах(Т) при импульсном нагружении 61

3.9 Результаты расчетов 68

Рекомендованный список диссертаций по специальности «Механика деформируемого твердого тела», 01.02.04 шифр ВАК

Введение диссертации (часть автореферата) на тему «Использование квазистатического подхода в динамике стержневых систем»

Задачи динамики балок и балочных систем можно отнести к «классическим» задачам механики деформируемого твердого тела (МДТТ). К настоящему времени разработан целый ряд подходов к решению этих задач. Эти подходы можно условно отнести к одной из трех групп: точные методы (основанные, как правило, на разложении колебаний по собственным формам) [1, 6, 34, 41-43], приближенные методы, использующие априорно задаваемые формы колебаний, (например, основанные на предположении о подобии формы колебаний форме статического прогиба) [9, 41-42, 53, 57] и численные методы (методы конечных разностей и конечных элементов) [2, 46]. Тем не менее, эти задачи по-прежнему остаются актуальными. Сравнительная простота таких задач делает их идеальным полем для отработки новых, перспективных методов. После такой отработки эти методы могут быть применены для решения более сложных задач МДТТ, как это было с численными методами.

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

Такой подход впервые был применен Герцем [10, 12, 52] при решении задачи о соударении шаров (данной задаче эквивалентна также задача соударения шара с твердой стенкой). Шары при ударе деформируются, в пятне контакта возникает сила соударения. Эта сила уравновешивается силамиинерции поступательного движения соударяющихся тел. При этом деформации считаются такими, какими они были бы в статике под действием силы соударения и сил инерции.

В работах С.А. Зегжды и В.Н. Вернигора [6, 19] квазистатический подход Герца был обобщен и применен к широкому кругу задач динамики упругих систем. В [11, 15-16, 39, 57] данный подход был применен к различным задачам динамики балок. В работах [16, 39] было предложено консольную балку рассматривать как свободный стержень под действием изгибающего момента и перерезывающей силы в заделке, которые должны быть такими, чтобы элемент балки, соответствующий заделке, оставался неподвижным. Таким образом, в рассмотрение вводятся две голономные связи. Свободный стержень может поступательно перемещаться, поворачиваться и изгибаться. И деформация стержня под действием сил инерции поступательного и вращательного движения считается квазистатической.

Данная работа является развитием этого направления. Ее особенностью является рассмотрение внутренних силовых факторов и реакций связей в качестве обобщенных лагранжевых координат. Такой подход позволяет представить балки с разными схемами закрепления в виде систем с конечным (24) числом степеней свободы. Для проверки корректности такой замены был решен ряд задач динамики собственных продольных и изгибных колебаний балок переменного сечения. Сравнение полученных результатов (первой частоты и формы - для продольных колебаний и первых двух частот и форм -для изгибных) с имеющимися аналитическими решениями показало высокую точность предлагаемого подхода (относительная погрешность определения первой частоты, как правило, не превосходит 1%).

Разработанный подход стал основой для решения задачи о динамике раскрытия трещины в тонком брусе. Задачи о развитии трещины под действием динамической нагрузки, являются актуальными и представляют большой научный и практический интерес. В настоящее время механика развития трещин является одной из самых актуальных и быстро развивающихся областей механики деформируемого твердого тела. Вместе с тем, в большинстве работ [20-21, 28-30, 32, 38, 40, 44-45, 47-51, 55-56] развитие трещины рассматривается как квазистатический процесс, который определяется, с одной стороны, накоплением повреждений в материале (работой сил разрушения), а с другой стороны - напряженно-деформированным состоянием (НДС) в окрестности вершины трещины. При решении динамических задач обычно учитывается нестационарность НДС, обусловленная изменением внешней нагрузки, а также - колебаниями рассматриваемой системы. При этом динамические эффекты, связанные собственно с раскрытием трещины, как правило, не учитываются или учитываются приближенно.

Специфической областью теории разрушения, является модель «балочного» приближения [7-8, 13-17, 26-27, 38-39, 49]. В этих работах рассматриваются задачи о развитии трещины в плоскости симметрии балки или пластины. В экспериментах для выполнения этого, несколько искусственного условия применялись специальные меры, например, использовались технологические надрезы, либо использовались кристаллы из анизотропных материалов. В этих условиях при больших (по сравнению с толщиной) длинах трещины каждый из образующихся берегов трещины можно рассматривать как балку-консоль, для которой вершина трещины играет роль заделки. В реальных случаях условие симметрии трещины не выполняется, однако решение задачи в такой упрощенной постановке является весьма важным и актуальным, поскольку оно позволяет исследовать общие закономерности динамики разрушения.

Впервые задача в такой постановке решалась Дж. Гилманом [8]. Он проводил как экспериментальные, так и теоретические исследования. Для раскрытия трещины в кристалле использовался клин, вводимый в трещину. При проведении расчетов форма прогиба консоли (берега трещины) считалась совпадающей с формой статического прогиба.

В [13] решается задача раскрытия трещины в брусе под действием сосредоточенных сил, приложенных к свободным краям берегов трещины. Каждый из берегов трещинь1 рассматривается как балка-консоль переменной длины Ь, которая принимается за обобщенную лагранжеву координату. Как показано в [13], поскольку при развитии трещины к рассматриваемой механической системе (консольной балке) новые точки присоединяются при нулевой скорости, к такой системе могут применяться уравнения Лагранжа второго рода [25, 31]. В [13] решение (функция прогиба консоли) представлено в виде суммы бесконечного ряда по балочным функциям. В отличие от консолей постоянной длины, для которых решение по каждой форме может быть получено независимо от остальных; несмотря на ортогональность балочных функций, данная задача сводится к системе дифференциальных уравнений, число которых стремится к бесконечности. В случае, если характерное время нарастания нагрузки, вызывающее раскрытие трещины, соизмеримо с периодами собственных колебаний консоли или больше их, можно ограничиться несколькими первыми формами или квазистатическим приближением, как это и было сделано в [13]. В противном случае, например, при импульсном нагружении, задача становится трудноразрешимой.

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

В работах В.Сысика [16, 39] задача решалась в двух постановках. В* первой, рассматривалось внедрение в трещину клина с заданной начальной скоростью. Движение берега (консоли) описывалось совокупностью балочных функций (совпадающих с функциями формы для консоли фиксированной длины). При этом первые п форм учитывались динамически, а остальные квазистатически (учитывалась только их потенциальная энергия, а кинетическая считалась равной нулю). Расчеты проводились для п= 1.

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

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

Экспериментальные исследования для случая симметричной трещины проводились на кафедре теории упругости Санкт-Петербургского государственного университета под руководством Ю.В. Петрова [4]. В этих экспериментах давление возрастает до своего максимума и падает до нуля за несколько микросекунд. Это время существенно меньше, чем время пробега волны по размеру берега трещины.

В данной работе рассматривается динамика раскрытия трещины при кратковременном (импульсном) нагружении ее берегов. При импульсном на-гружении возникают колебания по большому числу форм и неучет кинетической энергии высших форм приводит к значительным погрешностям. Вследствие малого времени воздействия к моменту его окончания кинетическая энергия консоли на несколько порядков превышает ее потенциальную энергию. Динамический учет большого числа форм связан с вычислительными трудностями (с увеличением номера формы резко увеличивается число зна)гчащих цифр, которые необходимо учитывать при вычислениях). Решить эту проблему позволяет обобщение квазистатического подхода на случай, когда деформация происходит под действием сил инерции, уравновешенных реакциями связей. Благодаря эффективному способу обезразмеривания переменных получена безразмерная система уравнений, не зависящая от числовых параметров (за исключением величин, характеризующих внешнюю нагрузку). Построенная математическая модель позволяет проанализировать результаты уникальных экспериментов, проводимых на кафедре теории упругости Санкт-Петербургского государственного университета под руководством Ю.В. Петрова [4]. Отмечен ступенчатый характер раскрытия трещины и задержка ее раскрытия после окончания действия импульсной нагрузки, что согласуется с результатами экспериментов.

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

Основные результаты, выносимые на защиту: • Показано, что, если свободный стержень переменного поперечного сечения мысленно разбить на две части и предположить, что каждая из частей под действием сил, приложенных к ней со стороны другой части, деформируется квазистатически, то при изгибных колебаниях система имеет четыре степени свободы, а при продольных - две. Полагая, что сечение, которым стержень разделяется на две части, является неподвижным, придем к рассмотрению консоли. Система имеет две степени свободы. Обобщенными координатами являются величины, пропорциональные изгибающе-5 му моменту и перерезывающей силе в заделке. Найденные приближенныезначения первой частоты для консолей в виде клина и конуса больше точных соответственно на 0,1% и 0,2%. При продольных колебаниях погрешность приближенного значения больше, чем при изгибных. В частности,для стержня в виде усеченного f клина с.заделанным концом приближенное значение:первой частоты выше точного на 1,1%.• Построена квазистатическая модель динамики раскрытия трещины в тонком брусе при импульсном нагружении ее исходных берегов. В* нулевом приближении модель имеет три степени свободы. Обобщенными* лагран-жевыми координатами являются изгибающий момент и перерезывающая сила в вершине трещины, а также длина трещины.• Предложен алгоритмпостроения последующих приближений, основанный; на добавлении к этим трем координатам новых координат, позволяющих динамически учесть несколько первых форм колебаний свободного стержня. Показано, что можно ограничиться моделью, имеющей шесть степеней свободы.• Показано, что раскрытие трещины носит ступенчатый характер.• - Определена минимальная, величина импульса, при котором начинает развиваться трещина.• Построена феноменологическая модель, позволяющая определить зависимость длины раскрытия трещины от импульса приложенной нагрузки;Для всех. рассматриваемых задач получены системы алгебраических и дифференциальных уравнений,, для решения которых использовался пакет «Mathematica-5.2», а для решения задачи о развитии трещины была составлена программа на языке Fortran 95.

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

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

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

Апробация работы. Полученные в работе результаты были представлены автором на следующих конференциях: Международная конференция «Четвертые Окуневские чтения», Санкт-Петербург, 2004 г. [23-24]; Международная научно-техническая конференция «Вычислительная- механика деформируемого твердого тела», Москва, 2006 г. [36]; Международная конференция по механике «Четвертые Поляховские чтения», Санкт-Петербург, 2006 г [37].

Результаты докладывались на семинарах кафедры теоретической и прикладной механики СПбГУ (2007 - 2008 гг.), а также на семинаре Института проблем машиноведения РАН «Нестационарные задачи механики и физики» (2007 г.), на секции теоретической механики в Доме Ученых им. М.Горького (2007 г.).

Публикации. По теме диссертации имеется 6 публикаций [17, 23-24, 3537], в том числе 1 статья в журнале, рекомендованном ВАК. В совместной работе [17] соавтору принадлежит постановка задачи и метод решения, автору принадлежит реализация предложенного метода и результаты расчетов.!Iч1. ИСПОЛЬЗОВАНИЕ РЕАКЦИЙ СВЯЗЕЙ КАК ОБОБЩЕННЫХ > КООРДИНАТ ПРИ МОДЕЛИРОВАНИИ КОЛЕБАНИЙСТЕРЖНЕВЫХ СИСТЕМВ данной главе показывается целесообразность использования реакций связей упругих механических систем как обобщенных лагранжевых коордиIV нат. Сделать это удается благодаря применению квазистатического подхода крешению динамических задач теории упругости [12, 57]. В настоящей главе такой подход используется для расчета собственных колебаний балок переменного сечения. В соответствии с определением обобщенными лагранже-выми координатами называются параметры, однозначно определяющие положение механической системы [33]. Следовательно, в случае если переме-, щение любой точки тела может быть однозначно выражено через реакции' связи, последние могут рассматриваться в качестве обобщенных лагранжевых координат.

Одним из первых квазистатический подход к решению динамических задач теории упругости применил Г.Герц. Он рассмотрел задачу о соударении упругих тел [12, 52]. Силу соударения Р в этой задаче можно рассматривать как реакцию связи, не позволяющей телам проникать друг в друга. В каждый момент времени сила Р уравновешивается силами инерции поступательного движения соударяющихся тел. Если считать соударение продолжающимся сколь угодно долго по отношению к периоду первой основной формы колебаний каждого из соударяющихся тел, то можно принять, что тела при соударении деформируются квазистатически. Это предположение позволило Герцу выразить силу соударения Р через относительное смещение а центров масс упругих тел (местное смятие). Таким образом, в модели, построенной Герцем, и кинетическая, и потенциальная энергии деформации соударяющихся тел выражаются как через Р и Р, так и через а и ос. Поэтому сила Р в модели Герца, имеющей одну степень свободы, может быть выбрана в качестве лагранжевой координаты. Идея Герца о целесообразности квазистатического подхода к решению динамических задач теории упругости нашла свое развитие в работах В.Н. Вернигора, С.А. Зегжды и М.П. Юшкова [6, 12, 18,57].

Использование реакций опорных элементов в качестве обобщенных ла-гранжевых координат в сочетании с использованием квазистатического подхода позволило с достаточно высокой точностью определить первые две собственные частоты и формы изгибных колебаний консолей переменного поперечного сечения.

1.1 Применение квазистатического подхода к задаче об изгибных колебаниях консолиРассмотрим задачу о собственных изгибных колебаниях консоли переменного сечения - рис. 1.1. Координата х отсчитывается от свободного конца балки, / - длина консоли, А(рс) и Дх) - площадь и момент инерции поперечного сечения х соответственно, М и <2 — реакции связей, равные перерезывающей силе и изгибающему моменту в заделке, у(х, /) - прогиб сечения х консоли.

А о —j—;—►Г fym х Рис. 1.1. Консоль переменного поперечного сечения.

Условия заделки консоли эквивалентны наличию двух голономных связейдх „,Используя принцип освобождаемости от связей, консоль будем рассматривать как свободную балку, к правому концу которой приложен изгибающий момент M{t) и перерезывающая сила Q{t). Свободная балка может, во-первых, перемещаться поступательно вдоль оси у и поворачиваться вокругоси, перпендикулярной плоскости ху, и, во-вторых, изгибаться. Будем сначала считать, что изгиб балки носит квазистатический характер, т.е. происходит под действием сил инерции поступательного и вращательного движения балки как абсолютно твердого тела, уравновешенных реакциями связей. Интенсивность сил инерции выражается через значения Mit) и Q(i), и потому прогиб любого сечения консоли в квазистатике однозначно определяется заданием величин Mit) и Q{t). Следовательно, в квазистатике изгибающий момент Mit) и перерезывающую силу Qit) можно рассматривать как обобщенные лагранжевы координаты. Предположение о квазистатической деформации данной балки под действием сил реакций Q(t) и M(t) и сил инерции позволяет выразить и кинетическую, и потенциальную энергию деформацииконсоли через величины Q(t), Q, M{t),M, т.е. позволяет построить модель с двумя степенями свободы. Для удобства при решении рассматриваемой задачи за обобщенные координаты принимаются величины:Л ) Ej(l) ■ - ' EJ(1) ' Здесь Е - модуль упругости, J{J) - момент инерции поперечного сечения балки на заделанном конце, т.е. при х=1 (см. рис. 1.1).

Введем безразмерные переменные, полагая:40' 40где А(х) и /(х) - соответственно площадь и момент инерции поперечного сечения л:.

Положение центра масс консоли переменного поперечного сечения определяется по формуле1А(г)с1гг, = 1с = 1—,J A(z)dzУскорение wc поступательного движения балки таково:„.(,)- = МОМО (1Л)р / 40 ¡A{z)dz р /4Л(/) pi(z) dz о оЗдесь р - плотность материала.

Отметим, что построенные функции h\{z) и h2(z), удовлетворяющие краевым условиям (1.7), можно рассматривать как удачно выбранные функции для консоли переменного сечения, позволяющие при рассмотрении ее колебаний использовать метод Рэлея-Ритца.

Как следует из выражений (1.7) и (1.8), первая форма собственных колебаний, соответствующая частоте соь имеет видЗдесь C,=l, С2 —с22 ®1 а221.2 Результаты расчетовВ качестве апробации изложенного метода использовались точные решения для консоли в виде клина единичной толщины (A(z) = z, j(z) = z3) иконуса ( A{z) = z2, j{z) = z4 ). Значение первой собственной частоты для клина, определенное Кирхгофом [42], таково: (Di=5,315.

Приближенное значение, вычисленное по описанной выше методике при помощи пакета «Mathematica»: cöi=5,318.

Таким образом, данное значение величины coj всего на 0,07% выше точного.

Точное значение первой собственной частоты для конуса [42]: coi=8,719.

Приближенное значение, вычисленное по описанной выше методике при помощи пакета «Mathematica»: (öi=8,735.

Таким образом, расчетное значение величины coi для конуса выше точного на 0,2%. Более высокое, в сравнении с клином, значение погрешности, связано с тем, что у клина масса убывает пропорционально z, а у конуса — пропорционально z2. Для консоли постоянного сечения будем поэтому иметь достаточно большую точность. Действительно, в этом случае точное значение: cöi=3}51602 [34], а приближенное: ©1=3,51604, т.е. погрешность составляет 5,7-10"4 %.

2. ИСПОЛЬЗОВАНИЕ ВНУТРЕННИХ СИЛОВЫХ ФАКТОРОВ КАК ОБОБЩЕННЫХ КООРДИНАТ ПРИ МОДЕЛИРОВАНИИ КОЛЕБАНИЙ СТЕРЖНЕВЫХ СИСТЕМ 2.1. Об использовании внутренних силовых факторов как обобщенных координат при моделировании колебаний стержневых системВ данной главе показывается целесообразность использования внутренних силовых факторов в характерных сечениях стержневых систем как обобщенных лагранжевых координат.

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

Следовательно, эти величины могут рассматриваться как обобщенные ла-гранжевы координаты.

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

Введем следующие обозначения: г, = хх/1, = х2/1 — безразмерные координаты (гх=2- г2); Ах(х,), А2(лг2) — площади поперечных сечений стержней 1 и 2 соответственно, = Ах{хх)/Ах{1) ; А2(г2)= А2(х2)/Ах(/) безразмерные площади поперечных сечений.

2.3 Применение квазистатического подхода к задаче об изгибных колебаниях свободного стержня, симметричного относительно серединыРассмотрим теперь изгибные колебания свободного симметричного стержня длины 21 - рис. 2.2. Как и в предыдущем примере разбиваем его на два стержня длины /.

2.4 Применение квазистатического подхода к задаче об изгнбных колебаниях несимметричного свободного стержняРассмотрим теперь несимметричный относительно середины стержень длины 21. Как и в предыдущих случаях, разбиваем его на два стержня длины / — рис.2.3. Снова будем отсчитывать продольную координату для каждой из половин от ее свободного конца. Теперь из-за несимметричности целого стержня придется рассматривать каждую его половину в отдельности, поэтому добавим к прежним следующие обозначения: ^(х,), У2 (х2) — моменты инерции поперечного сечения стержней 1 и 2 соответственно, J^(zl) = Jl{xl)/ J(l) J1(z1>)= J2{x■2)/ J(l) — безразмерные моменты инерции поперечного сечения.ООх.км м.

21Фfrм1 Sjyiw Х\ гп<—1--> стержень 1Х2<Ртстержень 2ОгРис. 2.3. Разбиение исходного несимметричного стержня на стержни 1 и 2 при расчете изгибных колебаний.

2.5. Задача об изгибных колебаниях свободного стержня с сосредоточенными массамиТеперь рассмотрим случай, когда на свободном стержне имеются сосредоточенные массы, жестко связанные с ним. Пусть в общем случае имеется щ масс ту, расположенных на стержне 1 в точках с координатами ху и п2 масс т2р расположенных на стержне 2 в точках с координатами Ху. — рис. 2.4.

М Мти Щп т21 Ш2пОт\\ Щп ОТЪПО 21 т2п ^21.-¿гРис. 2.4. Разбиение исходного стержня с сосредоточенными массами на стержни 1 и 2Будем рассматривать функции ^ 1(х() и Л2(х2) как обобщенные (содержащие 8-функции Дирака). Тогда имеем:>1 РЗдесь А0;(х;) - площадь поперечного сечения стержня без учета сосредоточенных масс, /-1,2. Тогда выражение для безразмерной площади поперечного сечения примет вид:.Мх,) & тч з( \Далее можно использовать формулы (2.15 - 2.20). В случае, если и сам целый стержень и сосредоточенные массы симметричны относительно центра, можно использовать формулы (2.9 - 2.13).

3. РАЗВИТИЕ ТРЕЩИНЫ В ТОНКОМ БРУСЕ ПРИ ИМПУЛЬСНОМ НАГРУЖЕНИИ3.1. Постановка задачиРассматривается динамика развития трещины в тонком брусе ширины b и толщины 2h. В брусе имеется технологический надрез длины L*, к берегам которого прикладывается импульсная нагрузка p(t)1 в виде давления, равномерно распределенная по всей площади берега. Размеры b и h считаем соизмеримыми друг с другом, а Ц » h. Толщиной технологического надреза посравнению с толщиной бруса можно пренебречь. С конца надреза начинает развиваться трещина. Расчетная схема нагружения приведена на рис.3.1.

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

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

3.2.). Исходное положение балки изображено сплошной линией, а пунктиром показано, как она изогнется под действием сил инерции ее поступательного и вращательного движения. Будем сначала считать, что изгиб балки носит квазистатический характер, т.е. происходит под действием сил инерции поступательного и вращательного движения балки как абсолютно твердого тела, уравновешенных реакциями связей. Интенсивность сил инерции выражается через значения M{t) и <2(0 > и потому прогиб любого сечения консоли в квазистатике однозначно определяется заданием величин M{t) и Q{t). Следовательно, в квазистатике перерезывающую силу 0{t), изгибающий моментM(i) и переменную длину балки L(t) можно рассматривать как обобщенные лагранжевы координаты. Заметим, что данная модель ни в квазистатическом, ни в последующих приближениях не позволяет описать смыкание берегов трещины, т.е. полагаем L > О. Также полагаем, что расчеты ведутся, пока 0 (см. рис.3.2).

Рис. 3.2. Расчетная схема распространения трещины в брусе: а - основные обозначения; б- принцип освобождаемости от связей.

При приложении статической нагрузки трещина начинает раскрываться тогда, когда изгибающий момент M(t) в ее вершине становится равен критическому значению [22]:M* = J2y.ЕЛ (3.1)Q Т-1 Г Ьк3Здесь Е — модуль упругости, J =--момент инерции поперечного сечения12балки, у - постоянная Гриффитса.

3.2. Определение прогиба консоли в квазистатическом приближенииПрогиб свободной балки длины b(t) может быть представлен в видеу(х, t) = д0 (г) + cp(iX* - L(t)) + £ q( {t)Xi( х ^/=1(3.2)где Xif x ^sMÙJ— балочные функции свободной балки [1, 34]:хА, =ShX; -SinX;L\t) chXi-cosXiЗдесь = l,.,oo) — возрастающая последовательность корней уравнения chXj cos Xi = l. Первое слагаемое в выражении (3.2) соответствует поступательному перемещению свободной балки, а второе — ее повороту вокруг конца х = L{t).

Отсюдадд{ дд1М = рЪЫ} — —рЫг1}д0, 0 =3 22 22 ф з:М(со2#(. = МЬ(3.6)(3.7)Здесь — значение обобщенной координаты в квазистатике.

Отметим, что, как показано в магистерской диссертации И.В.Мещерского [25] и (для более общего случая) в [31] для систем переменной массы, в случае, если новые массы присоединяются к системе при нулевой скорости, уравнения Лагранжа сохраняют обычный вид.

3.3. Представление прогиба консоли и определение кинетической и потенциальной энергии в «-ом приближенииИспользуя выражения (3.2), (3.3), (3.9), прогиб консоли в л-ом приближении представим в видеп< х ^/=1( х \(ш иг \X(3.13)Здесь У;.(г) = Х1.(г)-Х).(1)-Х/(1)(г-1). Величины Ь, Ак,к = 1,2, г = входящие в выражение (3.13), будем рассматривать как обобщенные лагранжевы координаты.

Похожие диссертационные работы по специальности «Механика деформируемого твердого тела», 01.02.04 шифр ВАК

Заключение диссертации по теме «Механика деформируемого твердого тела», Синильщикова, Галина Александровна

ЗАКЛЮЧЕНИЕ

I 1

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Синильщикова, Галина Александровна, 2008 год

1. Бабаков И.М. Теория колебаний. М.: Наука, 1968. 559 с.

2. Бате К., Вилсон Е. Численные методы анализа и метод конечных элементов. Пер: с англ. М.: Стройиздат, 1982.

3. Беляев Н.М. Труды по теории упругости и пластичности. М. 1957.

4. Березкин А.Н., Кривошеев С.И., Петров Ю.В., Уткин A.A. Эффект запаздывания старта трещины при пороговых импульсных нагрузках // Доклады Академии Наук. 2000; т.375. №3'. с. 328-330.

5. Бидерман B.JI. Теория удара. М. 1952.6: Вернигор В.Н, Михайлов A.JI. Модальный анализ механических колебаний упругих систем. Рыбинск 2001.

6. Георгиев И.Г. Динамика развития трещины в тонком брусе // Вестник молодых ученых. Серия: Прикладная математика и механика. СПб.: изд-во СПбГТУ, №4, 2000; С.92-101.

7. Гилман Дж.Дж. Скол, пластичность и вязкость кристаллов // Атомный механизм разрушения. М: Металлургиздат, 1963. С.220-253.

8. Дарков А.В:, Шпиро Г.С. Сопротивление материалов. М.: Высшая школа, 1975. 654 с.

9. Зегжда, С.А. Задача о соударении деформируемых тел//Прикладная механика. Вып.4. Л., 1979.

10. Зегжда С.А. Соударение кольца и балки // Колебания и устойчивость механических систем. Л., 1988.

11. Зегжда С.А. Соударение упругих тел. Изд-во СПбГУ. 1997, 316 с.

12. Зегжда С.А., Морозов Н.Ф., Семенов Б.Н. О «балочном» подходе в задачах распространения трещин // Механика твердого тела. №3.1999.

13. Зегжда С.А., Синилыцикова Г.А. Развитие трещины в тонком брусе при импульсном нагружении // Вестн. С.-Петербург, ун-та. Сер. 1. Вып.З. — 2007.-С. 15-23.

14. Зегжда С.А., Солтаханов Ш.Х., Юшков М.П. Уравнения движения не-голономных систем и вариационные принципы механики. Новый класс задач управления. М.: Наука, 2005. 268 с.

15. Зегжда С.А., Юшков М.П. Применение уравнений Лагранжа первого рода при исследовании собственных колебаний вала с дисками // Механика твердого тела. N4,1999. с.31-35.

16. Карзов Г.П., Марголин Б.З., Швецова В.А. Физико-механическое моделирование процессов разрушения. СПб.: Политехника, 1993. 391 с.

17. Косторов Б.В., Никитин Л.В., Флитман Л.М. Механика хрупкого разрушения. // Изв. АН СССР. МТТ 1972. №2 с.60-68.

18. Ландау Л.Д., Лившиц Е.М., Теоретическая физика т.7, М.: Наука. 1965. 204 с.

19. Мещерский И.В. Динамика точки переменной массы / Сб. И.В.Мещерский. Работы по механике тел переменной массы. М-Л. Гос-техиздат, 1949.

20. Михайлов А.М. Динамические задачи теории трещин в балочном при82ближении IIПМТФ. 1966. №5. С. 167-172.

21. Михайлов A.M. Обобщение балочного подхода к задачам теории развития трещин // ПМТФ.' 1969. №3. С. 171-175.

22. Молчанов А.Е. Никитин JI.B. Динамика трещин продольного сдвига после потери устойчивости // Изв. АН СССР. МТТ. 1972. №2. С.60-68.

23. Морозов Н.М., Никишков Г.П. Метод конечных элементов в механике , разрушения. М.: Наука, 1980. 254 с.30: Морозов Н.Ф. Паукшто М.В. Дискретные и гибридные модели механики разрушения. СПб.: изд-во СПбГУ, 1995. 157 с.

24. Новоселов B.C. Аналитическая механика систем с переменными массами. Л.: ЛГУ, 1969.

25. Партон В.З., Борисковский В.Г. Динамика хрупкого разрушения. М.: Машиностроение, 1988. 239 с.

26. Поляхов H.H., Зежда С.А., Юшков М.П. Теоретическая механика. М.:Высш. шк. 2000. 592 с.

27. Прочность. Устойчивость. Колебания. Справочник в трех томах. Т. 3. Под редакцией Биргера И.А. Пановко Я.Г. 1963. 736 с.

28. Синильщикова Г.А. Ударное расклинивание тонкого бруса// Международная конференции по механике «Четвертые Поляховские чтения». 710 февраля 2006 г., Санкт-Петербург: Тезисы докладов. СПб.: Изд-во «ВВМ», 2006.-С.215.

29. Слепян Л.И. Механика трещин. Л.Судостроение, 1990. 296 с.

30. Сысик В.П. Использование уравнений Лагранжа в теории удара и в динамике развития трещин. Автореферат диссертации на соискание ученой степени кандидата физико-математических наук. СПб. 2002.

31. Сиратори М., Миеси Т., Мацусита X. Вычислительная механика разрушения. М.: Мир, 1986. 334 с.

32. Тимошенко С.П. К вопросу о действии удара на балку // Изв. С.-Петербургского политех, ин-та. Т. 17. Вып.2. 1912.

33. Тимошенко С.П. Колебания в инженерном деле. М., 1967.

34. Филиппов А.П. Колебания деформируемых систем. М. Машинострое- * • ние. 1975,170 с.

35. Хеллан К. Введение в механику разрушения. М.: Мир, 1988. 364 с.

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

37. Хечумов P.A., Кеплер X., Прокопьев В.И. Применение метода конечных элементов к расчету конструкций. М.: Изд-во АСВ, 1994. 351с.

38. Шардин X. Исследование скорости разрушения // Атомарный механизм разрушения. М.: Металлургиздат, 1963. С. 297-330.

39. Achenback J.D., Bazant Z.R. Electrodynamics near-tip stress and* displaceiment fields for rapidly propagating cracks in anisotropic materials // Trans. ASME. Ser. E. Appl. Mech. 1975. V. 42. №1. P. 183-189.

40. Arakawa K., Takahashi K. Interlaminar fracture analysis of composite DCB specimens. // International Journal of Fracture 74, 1995.

41. Broberg K.B. The propagation of brittle crack //Ark. Fys. 1960. V. 18. №2. P.159-192.

42. Freund L.B., Clifton R. On the uniqueness of plane electrodynamics solution for running cracks // J. Elast. 1974. V.4. №4. P. 293-299.

43. Hertz H. Uber die berührung fester elastischer Korper//Zs. f. Math. (Crelle). 1881. Bd. 92.

44. Juschkov M.P. Anwendung der Lagrangeschen Gleichungen 1. An zur Untersuchung der nichtlinearen Querschwingungen von Balken mit unverschieblichen Lagern. //Technische Mechanik, Bd 18, HI, (1998), P.79-84.

45. Yoffe E.H. The moving Griffith's crack. // Phil. Mag. Ser. 7. V. 42, No. 330, pp.739-750,* 1951.

46. Viggo Tvergaard. Cleavage crack growth resistance due to plastic flow around a near-tip dislocation-free region. // J. Mech. Phys. Solids, Vol. 45, No. 6, pp. 1007-1023, 1997.

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