Перейти к публикации

Полезная информация для всех расчетчиков


Рекомендованные сообщения

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

Ссылка на сообщение
Поделиться на других сайтах


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

Рассмотрим вынужденные колебания механической системы, вызванные внешними силами.

Скрытый текст

2.1.4 Вынужденные колебания

 

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

image004.gif                        (2.1.4.1)

Пусть нам известны собственные частоты image006.gif и собственные формы image008.gif первых image010.gif колебательных мод image012.gif. Так  как  во  многих задачах число степеней свободы механической системы image014.gif, то определяемое на практике число первых колебательных мод image016.gif. Первые моды соответствуют низшим частотам механической системы и поэтому во многих задачах являются определяющими. Поэтому искомое решение "уравне­ний движения" (2.1.4.1) image018.gif можно приближенно представить в виде суперпозиции первых image010.gif колебательных мод [6]:

image020.gif                               (2.1.4.2)

Собственные формы колебаний image008.gif удовлетворяют следующим условиям (см. п. 2.1.3):

image023.gif                                                                      (2.1.4.3)

Подставив разложение image025.gif по колебательным модам (2.1.4.2) в (2.1.4.1), получим:

image027.gif

откуда после домножения слева на image029.gif с учетом (2.1.4.3) получим систему image031.gif независимых дифференциальных уравнений относительно image033.gif:

image035.gif       (2.1.4.4)

где: image037.gif        (2.1.4.5)

Считая, что при image039.gif механическая система покоится в положении равновесия, получим следующие начальные условия для уравнений (2.1.4.4):

image041.gif                   (2.1.4.6)

Заметим, что уравнения (2.1.4.1) и (2.1.4.4) справедливы только для консервативной системы. В общем случае надо  учитывать  присутствующее в механической системе трение. Для силы трения, пропорциональной скорости, уравнения (2.1.4.4) следует записать в виде:

image043.gif       (2.1.4.7)

Здесь введены коэффициенты трения image045.gif. Рассмотрим общий метод аналитического решения задачи Коши:

image047.gif                    (2.1.4.8)

Пусть сила image049.gif действует в течение бесконечно малого промежутка времени image051.gif и передает единичный "импульс" image053.gif:

image055.gif                             (2.1.4.9)

Произведя интегрирование (2.1.4.8) на промежутке времени image051.gif с правой частью image058.gif, найдем, что "скорость" image060.gif испытывает при image062.gif скачек:

image064.gif               (2.1.4.10)

Так как при image066.gif правая часть равна нулю, то для image060.gif получаем следующую задачу Коши:

image068.gif

(2.1.4.11)  (2.1.4.12)

Решение однородного уравнения (2.1.4.11) ищем в виде:

image070.gif                                       (2.1.4.13)

Подставив (2.1.4.13) в (2.1.4.11), получим:

image072.gif                          (2.1.4.14)

Уравнение (2.1.4.14) будет справедливо при значениях image074.gif:

image076.gif                                                                    (2.1.4.15)

Здесь предполагается, что image078.gif (малое трение).

Поэтому общее решение (2.1.4.11) имеет вид:

image080.gif                       (2.1.4.16)

Используя начальные условия (2.1.4.12), найдем image082.gif и image084.gif:

image086.gif                        (2.1.4.17)

С помощью формул (2.1.4.15) - (2.1.4.17) получим следующий вид решения image088.gif задачи Коши (2.1.4.11) - (2.1.4.12):

image090.gif (2.1.4.18)

Решение image092.gif для импульсной силы image094.gif будем называть функцией Грина дифференциального уравнения  (2.1.4.11). Для линейных уравнений справедлив принцип суперпозиции решений. Так как для произвольной правой части image096.gif верно интегральное представление:

image098.gif                     (2.1.4.19)

то из принципа суперпозиции следует, что решение (2.1.4.11) с начальными условиями (2.1.4.12) и правой частью image058.gif равно:

image101.gif                  (2.1.4.20)

Здесь: image103.gif - функция Грина (2.1.4.18);

image105.gif - правая часть (2.1.4.11).

Таким образом, мы нашли решение задачи Коши (2.1.4.11) - (2.1.4.12):

image107.gif                                                                    (2.1.4.22)

где: image109.gif

Рассмотрим несколько примеров. Пусть правая часть image096.gif имеет ступенчатый характер:

image112.gif                           (2.1.4.23)

Подставив правую часть (2.1.4.23) в (2.1.4.22), получим:

image114.gif

(2.1.4.24)

Для случая гармонической силы image116.gif, в отсутствии трения image118.gif, найдем:

image120.gif                                                                    (2.1.4.25)

В случае возбуждения, меняющегося по гармоническому закону image122.gif, решение задачи Коши (2.1.4.8) можно получить более просто. Общее решение неоднородного линейного дифференциального уравнения (2.1.4.8) равно сумме общего решения однородного уравнения (уравнения (2.1.4.8) без правой части) и частного решения неоднородного уравнения. Для нахождения искомого решения удобно перейти к следующему дифференциальному уравнению относительно комплексной перемен­ной image124.gif:

image126.gif               (2.1.4.26)

где: image128.gif

image130.gif - комплексная сила image132.gif.

Если взять мнимую часть  уравнения (2.1.4.26), то получим уравнение (2.1.4.8). Найдем сначала общее решение однородного уравнения:

image134.gif                        (2.1.4.27)

Ищем частное решение линейного однородного уравнения (2.1.4.27) в виде:

image136.gif                                     (2.1.4.28)

После подстановки image138.gif(2.1.4.28) в уравнение (2.1.4.27) в силу линейности уравнения получим следующее квадратное  уравнение для определения частот image140.gif:

image142.gif                       (2.1.4.29)

откуда найдем два корня image144.gif:

image146.gif                                                                    (2.1.4.30)

Здесь введена частота:

image148.gif                                (2.1.4.31)

Используя формулу (2.1.4.28), получим общее решение однородного уравнения (2.1.4.27) как суперпозицию двух  линейно  независимых частных решений с произвольными комплексными множителями image150.gif:

image152.gif          (2.1.4.32)

Тогда общее решение соответствующего однородного уравнения относительно действительной переменной  может  быть представлено в виде image154.gif:

image156.gif                                                                    (2.1.4.33)

где: image158.gif и image160.gif - произвольные действительные числа.

Найдем теперь частное решение неоднородного уравнения (2.1.4.26). Будем искать решение в виде:

image162.gif                                     (2.1.4.34)

 Для амплитуды image164.gif получим с помощью (2.1.4.26):

image166.gif                           (2.1.4.35)

 

Поэтому частное решение (2.1.4.34) будет равно:

image168.gif               (2.1.4.35)

 

где: image170.gif - сдвиг фазы.

Сдвиг фазы image170.gif может быть определен из следующей формулы:

image172.gif                                                                    (2.1.4.36)

Если image174.gif (демпфирование отсутствует), то формула (2.1.4.36) значительно упрощается:

image176.gif                     (2.1.4.37)

Таким образом, при малых значениях возбуждающей частоты image178.gif изме­нение силы image058.gif и перемещения image060.gif происходят в одной фазе image182.gif, а при больших значениях image184.gif в противофазе image186.gif.

Поэтому общее решение уравнения (2.1.4.8) будет равно:

image188.gif

(2.1.4.38)

Значения констант image158.gif и image160.gif определяются по начальным условиям задачи Коши (2.1.4.8):

image192.gif                              (2.1.4.39)

Присутствие затухающего множителя image194.gif приводит к тому,  что при  достаточно  больших  значениях  времени image196.gif решение (2.1.4.38) равно:

image198.gif               (2.1.4.40)

то есть не зависит от начальных условий! Следовательно, после окончания переход­ного процесса, определяющегося начальными условиями, колебания устанавливаются и становятся гармоническими со сдвигом фазы image170.gif относительно возбуждающей силы image058.gif. Амплитуда установившихся колебаний (2.1.4.40) равна:

image201.gif                  (2.1.4.41)

 

Так как обычно image203.gif то амплитуда как функция возбуждающей частоты image074.gif имеет ярко выраженный максимум при image206.gif, где: image140.gif - собственная частота колебаний.

Заметим, что в случае произвольной периодической силы:

image209.gif, где: image211.gif - период  гармонических  колебаний) ее можно разложить в ряд Фурье:

image213.gif                                                                    (2.1.4.42)

где: image215.gif - ряд возбуждающих частот.

Для image217.gif предельное решение, отвечающее окончанию переходного процес­са, равно (см. (2.1.4.24)):

image219.gif                                           (2.1.4.43)

Следовательно, можно сделать вывод, что в случае произвольной периодической возбуждающей силы image058.gif сначала имеет место переходный процесс, определяющийся начальными условиями. После завершения переходного процесса устанавливается периодическое решение image222.gif, не зависящее от начальных условий. Это решение может быть представлено следующим рядом Фурье:

image224.gif                                                                    (2.1.4.44)

где: image226.gif - значения функции image228.gif (2.1.4.36) при image230.gif.

 

 

 

Если объяснять на пальцах, то при переходе в систему координат собственных векторов, матрицы [М] и [К] будут диагональные. И для каждого уравнения можно найти точное решение для некоторых видов нагрузок. Далее если сделать обратный переход – получаем решение в обычных координатах.

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

На этом пальцы закончились.:beee:

Ссылка на сообщение
Поделиться на других сайтах

Вроде это штатный метод решения систем ДУ через подстановку решения в виде экспоненты от матриц и описан с куче букварей. Лень даже в старые тетрадки с конспектами заглядывать :)

Ссылка на сообщение
Поделиться на других сайтах
15 часов назад, Fedor сказал:

Лень даже в старые тетрадки с конспектами заглядывать

В таких вопросах не нужно лениться.:)

Это будет полезная информация для всех расчетчиков.

Изменено пользователем ДОБРЯК
Ссылка на сообщение
Поделиться на других сайтах

Бесполезная если они не собираются в программисты подаваться которые пишут программы про математику .   :) 

Достаточно знать что есть две проблемы алгебры - решение систем и нахождение корней полиномом. Собственные числа - это вторая. Смысл собственных векторов в том, что существует направление движения деформируемого тела в котором  упругие силы и инерционные совпадают по направлению  F = CU    и Fi= MU  а коэффициенты связей этих векторов и есть собственные числа и чтобы их найти надо найти корни полинома то есть решать вторую алгебраическую проблему. В силу ее нелинейности - с мужеством и оптимизмом итерировать :) 

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах

Если матрицу масс [М] заменить на матрицу [С] (геометрическую матрицу) и решить туже задачу  с мужеством и оптимизмом, то найдем точки бифуркации и формы потери устойчивости.

То есть решим задачу начальной потери устойчивости. 

 

Изменено пользователем ДОБРЯК
Ссылка на сообщение
Поделиться на других сайтах
Цитата

 [С] (геометрическую матрицу)

Вообще-то ее называют матрицей жесткости :)

Если заменим то получим (1-лябда) CU=0   . Ясно что вектор совпадает с самим собой по направлению  CU= лямбда CU :)

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах
44 минуты назад, Fedor сказал:

Вообще-то ее называют матрицей жесткости

Матрица жесткости [K], матрица масс [M]. 

5 часов назад, ДОБРЯК сказал:

Если матрицу масс [М] заменить на матрицу [С] (геометрическую матрицу) ...

 

Ссылка на сообщение
Поделиться на других сайтах

Можно и так смотреть , мелкий вопрос. У Галлагера все расписано насколько помню  :)

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах

Информация о собственных колебаниях в параллельной ветке частично потеряна.

Обновил здесь.

Скрытый текст

2.1.3 Собственные колебания

 

Собственные колебания механической системы совершаются при отсутствии внешних сил. Полагая в "уравнениях движения" (2.1.1.18) правую часть image002.gif равной нулю, получим уравнения собственных колебаний [6]:

image004.gif                             (2.1.3.1)

Решение уравнений (2.1.3.1) можно искать в виде гармоник с частотой колебания image006.gif:

image008.gif                                                                      (2.1.3.2)

Так как справедливы равенства:

image010.gif

то после подстановки image012.gif или image014.gif в уравнения (2.1.3.1) получим:

image016.gif                         (2.1.3.3)

В уравнении (2.1.3.3) члены, зависящие от времени, опущены. Система уравнений (2.1.3.3) имеет фундаментальное значение в теории колебаний. Ее решение позволяет найти собственные частоты image018.gif и собственные формы колебаний image020.gif механической  системы, image022.gif, где: image024.gif - число степеней свободы.

Общее решение системы уравнений (2.1.3.1) имеет вид:

image026.gif

(2.1.3.4)

Заметим, что в силу положительной определенности матриц image028.gif и image030.gif из (2.1.3.3) следует:

image032.gif                              (2.1.3.5)

т.е. собственные частоты image034.gif, (image022.gif) являются действительными (положи­тельными) числами, как и должно быть для консервативной системы. За исключением особых случаев, например, наличия определенной симметрии механической системы, все собственные частоты различны.

Покажем, что если

image037.gif           (2.1.3.6)

то справедливы следующие формулы ортогональности собственных форм колеба­ний:

image039.gif     (2.1.3.7)

Действительно, домножая векторные равенства

image041.gif

слева на image043.gif соответственно, получим:

image045.gif       (2.1.3.8)

Так  как  матрицы image028.gif и image030.gif симметричные, то из (2.1.3.8) следует:

image049.gif              (2.1.3.9)

Но image051.gif (см. (2.1.3.6)). Поэтому: image053.gif т.е. доказана справед­ливость второй формулы из (2.1.3.7). Отсюда с помощью (2.1.3.3) можно доказать справедливость первой формулы из (2.1.3.7). Для однозначного определения векторов собственных форм image055.gif необходимо ввести их нормировку. Введем сле­дующую нормировку собственных форм:

image057.gif           (2.1.3.10)

В этом случае справедливы равенства:

image059.gif         (2.1.3.11)

Таким образом, задача нахождения собственных колебаний механической системы свелась к задаче определения собственных частот и собственных форм для системы алгебраических уравнений (2.1.3.3):

image061.gif                         (2.1.3.12)

Методы решения задачи (2.1.3.12) рассмотрены в главе 8.

 

Ссылка на сообщение
Поделиться на других сайтах

Зачем учебники копировать на форум, замусоривать ?  Можно же просто ссылку дать :)

Ссылка на сообщение
Поделиться на других сайтах
3 минуты назад, Fedor сказал:

Зачем учебники копировать на форум, замусоривать ?  Можно же просто ссылку дать :)

Нет ссылки. Уже объяснял.

 

Ссылка на сообщение
Поделиться на других сайтах

 

1 час назад, ДОБРЯК сказал:

Таким образом, задача нахождения собственных колебаний механической системы свелась к задаче определения собственных частот и собственных форм для системы алгебраических уравнений (2.1.3.3):

Ничего у вас не свелось. Вы просто нашли собственные частоты и формы. Вам еще нужно найти "a" и "b" в выражении 2.1.3.4.

 

Ссылка на сообщение
Поделиться на других сайтах
9 минут назад, Fedor сказал:

Учебники все есть в сети. Или похожие. Не смешите с утра :) 

Не все учебники полезны для расчетчиков.

http://old.exponenta.ru/soft/Mathemat/pinega/main.asp

Но полезны для разработчиков. :)

 

Ссылка на сообщение
Поделиться на других сайтах

PS. И правда, не смешите математиков. Это же элементарные вещи курса Линейной_Алгебры.. и многих других.

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, Борман сказал:

PS. И правда, не смешите математиков. Это же элементарные вещи курса Линейной_Алгебры.. и многих других.

Так вы математик.  Теперь понятно почему вы пишите про одну нулевую собственную частоту незакрепленной конструкции в параллельной ветке. :beee:

 

Ссылка на сообщение
Поделиться на других сайтах
12 минуты назад, ДОБРЯК сказал:

Теперь понятно почему вы пишите про одну нулевую собственную частоту незакрепленной конструкции в параллельной ветке. :beee:

Опять буквы начали путать ?

Ссылка на сообщение
Поделиться на других сайтах
13 минуты назад, Борман сказал:

Опять буквы начали путать ?

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

И нигде вы не поправили Соклакова, что нулевых частот шесть. Значит вы с ним согласны.:5a33a36a07342_3DSmiles(142):

А иначе зачем такие глупости выносить на самое видное место. И это не единственная глупость. 

Изменено пользователем ДОБРЯК
Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




×
×
  • Создать...