Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике Во многих важных для практики задачах функции, описывающие возмущающие силы, не удается выразить в аналитическом виде, поэтому их представляют либо в виде набора точек на диаграмме, либо в виде таблиц. В подобных случаях иногда можно аппроксимировать исходные данные с помощью формул, применяемых в методах построения кривых по точкам, и затем подставлять полученные зависимости в интеграл Дюамеля. Однако более общий подход для определения динамического поведения систем состоит в использовании некоторых простых интерполяционных функций в периодически повторяющихся сериях вычислений. Последний упомянутый метод численного исследования и будет обсуждаться в данном параграфе применительно к нескольким типам интерполирующих функций. Предположим, что на систему с одной степенью свободы с демпфированием (рис. 1.54) действует сила $Q$, изменяющаяся во времени некоторым произвольным образом, аналогичным показанному на рис. 1.55. Эту непрерывную функцию, описывающую возмущающую силу, можно приближенно представить в виде набора ступенчатых функций с различными значениями в различные моменты времени, как показано на рис. 1.55. Первое значение ступенчатой функции равно $\Delta Q_{0}$ в момент времени $t=0$, второе $-\Delta Q$ в момент времени $t=t_{1}$ и т. д. В произвольный интервал времени $t_{i-1} \leqslant t<t_{i}$ динамическое поведение системы при воздействии возмущающей силы, представляемой указанной ступенчатой функцией, можно описать выражением (вм. пример 3 из п. 1.12) В момент времени $t_{i}$ это перемещение что в случае отсутствия демпфирования дает В данном методе величина $\Delta Q_{i}$ некоторой типичной степени может быть либо положительной, либо отрицательной, что зависит от угла наклона касательной к рассматриваемой кривой. Для получения достаточной точности в рассматриваемом методе следует выбирать достаточно малые шаги и самокомпенсирующиеся погрешности площади области, лежащей под графиком функции, описывающей возмущающую силу. Речь идет о том, чтобы заштрихованные площади на рис. 1.55, лежащие выше кривой, были примерно равны незаштрихованным площадям, лежащим ниже кривой. Использование такого приема означает, что каждый (после первого) шаг начинается в момент времени, когда ордината кривой равна средней высоте ординат на заданном шаге. Это можно видеть на рисунке. Разумеется, если функция, описывающая силу, представляет импульс, действительно ограниченный горизонтальной и вертикальными линиями, то метод приведет к точному результату. Другой метод состоит в использовании линий, параллельных осям координат (рис. 1.56). В этом случае кривая аппроксимируется рядом импульсов прямоугольной формы, различной величины и длительности. Для получения достаточной точности величина $Q_{i}$ типичного импульса должна быть выбрана такой, чтобы она равнялась ординате кривой в середине временного интервала $\Delta t_{i}$, как показано на рис. 1.56. В произвольный интервал времени $t_{i-1} \leqslant t \leqslant t_{i}$ реакцию системы с одной степенью свободы с демпфированием можно вычислить, представив ее как сумму влияний начальных условий в момент времени $t_{i-1}$ и влияния импульса, действующего на интервале времени $\Delta t_{i}$, что дает В конце интервала это выражение принимает вид что в случае отсутствия демпфирования дает Кроме того, можно найти скорость $\dot{x}_{i}$ в конце интервала времени, которая, будучи поделенной на частоту $p$, имеет вид Выражения (1.76в) и (1.76г) представляют собой рекуррентные формулы для определения динамического перемещения системы при отсутствии демпфирования в конце $i$-го шага и тем самым начальных условий в начале шага ( $i+1)$. Последовательно используя эти формулы, можно проследить, как изменяются во времени перемещения и скорость системы с одной степенью свободы, но больший интерес представляет перемещение *. Полученные выше рекуррентные формулы (1.76в) и (1.76г) позволяют находить перемещения в конце $i$-го интервала времени путем последовательных вычислений. Другой подход заключается в определении перемещений (в момент времени $t_{i}$ ) при действии всех предыдущих импульсов прямоугольной формы. Для случая колебаний без демпфирования данный подход основывается на использовании следующей формулы: Поскольку в последнем слагаемом формулы (1.76д) необходимо производить суммирование членов ряда, при использовании этой формулы требуется выполнять большее число арифметических операций, чем по формулам (1.76в) и (1.76г). Поэтому с целью облегчения проведения расчетов при определении окончательного вида зависимости перемещений от времени предпочтительнее использовать выражения (1.76в) и (1.76г). Однако, если требуется определить только перемещения в конкретное время, то лучше выбрать формулу (1.76д). При использовании интерполяции кусочно-постолнного типа, описанной выше, не всегда удобно делать равными погрешности площади областей, лежащих над графиком функции возмущающей силы и под ним. Более грубым подходом является выбор ординат кривой, относящихся к началу (или концу) интервала времени, в качестве значения импульса прямоугольной формы (или ступенчатой функции). При этом для сохранения заданной точности решения может потребоваться большее число шагов по времени, и при вычислении может стать значительной ошибка округленил. Для того чтобы избежать указанных трудностей, можно воспользоваться интерполирующими функциями более высокого порядка. На рис. 1.57 показан логически вытекающий из сказанного способ представления импульсного возмущения с помощью наклонных линий и вертикальных полос. Для этой интерполяции кусочно-линейного типа переме- где $\Delta Q=Q_{i}-Q_{i-1}$. Последнее слагаемое выражения (1.77a) переписано из решения для задачи о колебаниях с демпфированием, когда возмущающая сила описывается линейной функцией (см. задачу 1.12.9). В конце $i$-го интервала времени выражение (1.77a) будет иметь форму Если демпфированием пренебречь, из этого выражения получим откуда находим выражение для скорости Выражения (1.77в) и (1.77г) представляют собой рекуррентные формулы, аналогичные формулам (1.76в) и (1.76г) для импульса прямоугольной формы. Для того чтобы определить перемещение только в момент времени $t_{i}$ для случая отсутствия демпфирования при кусочно-линейной интерполяции, можно взять следующий вариант формулы: где последняя группа слагаемых взята из задачи 1.13.4. этих рисунках применяются постоянные шаги по времени $\Delta t=t_{1} / 10$. Первые три случая относятся к импульсам прямоугольной формы, тогда как последний к импульсу трапецеидальной формы (с использованием кусочно-линейной интерполяции). В подходах, к которым относятся рис. $1.58, a-\varepsilon$, величины импульсов определялись значениями ординат кривой соответственно в начале, конце и середине шага. Для удобства сравнения указанных случаев была выбрана описываемая уравнением кривая, для которой известно точное решение при отсутствии демпфирования (см. задачу 1.13.7): где В качестве специального случая положим $t_{1}=\tau / 2$, и тогда выражение (б) примет вид Вычислить и представить графически изменения перемещений во времени для системы при колебаниях без демпфирования на интервале времени $0 \leqslant t<t_{1}$ для четырех способов представления кривых, показанных на рис. 1.58, a-2. Предполагается, что начальные условия при $t=0 x_{0}=0$ и $\dot{x}_{0}=0$ и что величины $Q_{1}$ и $k$ равны единице. Решение. Для первых трех случаев воспользуемся выражениями (1.76в) и (1.76г); в свою очередь, выражения (1.77в) и (1.77г) применимы для случая 4. Согласно исходным параметрам в этом примере имеем $\Delta t_{i}=\tau / 20, p \Delta t_{i}=\pi / 10=$ $=18, \cos p \Delta t_{i}=0,951, \sin p \Delta t_{i}=0,309,1-\cos p \Delta t_{i}=0,0489, p \Delta t-\sin p \Delta t_{i}=$ $=0,00514$. Вначале можно заполнить первые шесть столбцов в таблице, поместив начальные значения $x_{0}$ и $\dot{x}_{0} / p$ (в данном примере равные нулю) на первой строке столбцов 1.1. Решение для случая 1 из примера 1 с номерами соответственно 11 и 12. Тогда величины для второй строки столбцов 7-10 подсчитываем по следующим формулам: 1.2. Перемещения для примера 1 В этой точке величины $x_{1}$ и $\dot{x}_{1} / p$ можно подсчитать по выражениям (1.76в) и (1.76r), после чего подставить соответственно в столбцы 11 и 12; при этом имеем Все последующие строки таблицы записываем по тем же правилам, что и для первой строки. Аналогичным образом можно рассмотреть остальные три случая из данного примера, полученные таким образом для всех четырех случаев перемещения. Их точные значения приведены в табл. 1.2. Графики, построенные по полученным результатам, приведены на рис. 1.59, где сплошной линией представлено точное решение. Қак и следовало ожидать, кривая, относяцаяся к случаю 1 (см. рис. $1.58, a$ ), лежит выше кривой для точного решения, тогда как в случае 2 (см. рис. 1.58, б) – кривая лежит ниже точного решения; в случаях 3 и 4 (см. рис. 1.58 , в и а) эти кривые практически совпадают с точным решением, но из сопоставления числовых данных из табл. 1.2 следует, что оба способа интерполяции дают значения, несколько меньшие, чем точные. Пример 2. Пусть система с одной степенью свободы без демпфирования имеет жесткость пружины $k=1,787 \cdot 10^{3} \mathrm{H} /$ м и период собственных колебаний $\tau=1,2 \mathrm{c}$. Возмущающая сила $Q$ изменяется во времени так, как показано на рис. 1.60, $a$. Числовые данные для силы $Q$ приведены в столбце 3 табл.1.3. Как видно из рис. $1.60, a$, в качестве интерполирующей выбрана кусочно-постоянная функция с постоянным шагом по времени $\Delta t_{i}=0,1$ с. В качестве значения каждого импульса прямоугольной формы выбрано значение ординаты кривой в середине шага. Кроме того, предполагается, что начальные условия в данной задаче: $x_{0}=0$ и $\dot{x}_{0}=0$ при $t=0$. Вычислить и построить график перемещения системы на интервале времени $0 \leqslant t \leqslant 1,6$ с. Решение. В рассмотренном примере имеем $\Delta t_{i} / \tau=1 / 12, p \Delta t_{i}=\pi / 6=30^{\circ}$, $\cos p \Delta t_{i}=0,866, \sin p \Delta t_{i}=0,5,1-\cos p \Delta t_{i}=0,134$. В табл. 1.3 приведены все этапы вычислений, а на рис. $1,60,6$ построен график зависимости перемещения от времени. 1.3. Решения для примера 2 ЗАДАЧИ* 1.15.1. Определить выражение для скорости $\dot{x}_{i}$, поделенной на частоту $p$, для случая, когда реакция системы при колебаниях без демпфирования определяется выражением (1.76д).
|
1 |
Оглавление
|