![]() |
Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | |
|
С целью снижения материалоемкости круглые пластины часто изготавливают с переменной по радиусу толщиной. Для расчета таких пластин используют численный метод начальных параметров. При решении задачи будем определять в каждой точке пластины три основные величины - прогиб , угол поворота нормали
, и произведение момента
на текущий радиус. Составив матрицу- столбец из этих величин
)T, назовем её вектором состояния. Действительно, три компоненты этого вектора определяют не только перемещения
,
, но и изгибающие моменты, а, следовательно, и напряжения. В самом деле, подставив во второе уравнение (16.6)
из первого, вы разим
через
и
:
Равенства (16.5), (16.6), (16.10) позволяют записать три дифференциальных уравнения первого порядка относительно компонент вектора состояния:
![]() | ![]() | ![]() |
![]() | ||
![]() |
Первое из уравнений (16.11) совпадает с уравнением (16.10), второе есть разрешенное относительно первое из равенств (16.6), а третье - уравнение равновесия (16.5), в которое подставлено значение
по формуле (16.11). Поскольку левые части уравнений (16.12) - производные компонент вектора со стояния, а правые части линейно зависят от этих компонент, систему уравнений (16.12) можно записать в виде одного матричного уравнения:
где матрица переменных коэффициентов F и вектор g выражаются формулами:
Решение системы уравнений (16.13) должно удовлетворять граничным условиям. Так как эта система третьего порядка, то число граничных условий равно трём. При этом на каждом из контуров пластины (внутреннем и внешнем) должно быть задано ограничение либо на угол поворота q, либо на соответствующий этому углу изгибающий момент , либо, при упругой заделке, граничное условие ставится на линейную комбинацию угла q и момента
. Третьим граничным условием является условие
на окружности, где расположена опора. В табл.16.1 приведены возможные варианты граничных условий (ось z направлена вниз).
Таблица 16.1
Для пластины без центрального отверстия граничное условие в центре
() имеет вид:
Однако при численном интегрировании это условие не используется, так как при некоторые коэффициенты системы (16.13) обращаются в бесконечность. В связи с этим начальное значение независимой переменной r принимают равным малому числу а (а - "разгонный участок" обычно равный шагу интегрирования), а граничное условие
заменяют следующим:
где D - жесткость на изгиб в центре пластины.
Приведем краткий вывод условия (16.14). Вблизи центра пластины угол q пропорционально радиусу (последующие члены разложения q в степе нной ряд можно при малых r не учитывать). Поэтому при малых r имеем
,
и из первого уравнения (16.6) следует условие (16.14).
Таким образом, необходимо найти численное решение системы дифференциальных уравнений (16.13), удовлетворяющее граничным условиям, два из которых наложенные на q или , или на их комбинацию, формулируются на внутреннем и наружном контурах пластины, и одно, наложенное на прогиб - в месте ограничения прогиба.
Для решения этой краевой задачи используем метод начальных параметров [4], который позволяет наиболее простым способом свести краевую задачу к задаче Коши. Система уравнений (16.13) представляет собой систему линейных неоднородных (из-за наличия слагаемого с Q) дифференциальных уравнений третьего порядка. Известно, что решение такой системы складывается, вообще говоря, из трех линейно независимых (
) решений однородной системы и одного
решения неоднородной системы. Неоднородное решен
ие удовлетворяет уравнению (16.13), а все однородные - уравнению
Если вычислить все решения по отдельности и, подчинив выражение
граничным условиям, определить постоянные , то задача будет решена.
Метод начальных параметров позволяет сократить количество вычисляемых решений однородной задачи. Действительно, так как, по крайней мере, одно граничное условие задано в начале интервала интегрирования (при или
- для пластинки без отверстия), то из трех постоянных
независимыми оказываются лишь две, и решение, удовлетворяющее указанному граничному условию, можно представить в виде:
При этом начальные значения должны быть заданы так, чтобы граничное условие при
выполнялось при любых
и
.
Дополнительное упрощение связано с тем, что одно из решений однородной задачи изгиба пластины очевидно - оно соответствует перемещению пластины как жесткой () при
,
Таким образом, окончательно
Рассмотрим теперь вопрос о выборе значений векторов ,
, обеспечивающих тождественное выполнение наложенных на q или
граничных условий при
:
1. Поворот (контур
жестко заделан или пластина связана с жестким центром радиуса
). Условие выполняется, если принять
2. Момент (контур пластины
свободен или шарнирно оперт), в этом случае
3. Внутренний контур () нагружен распределенным моментом
:
4. Пластина без отверстия. На границе "разгонного участка" () принимают:
где - жесткость пластины в центре.
Постоянную интегрирования находят, подчиняя вектор решения на внешнем контуре (
)
граничному условию, наложенному на угол поворота q или радиальный момент . Например, в случае свободного внешнего контура
постоянная
определяется из следующего алгебраического уравнения:
(здесь первый из индексов указывает номер компоненты, второй-номер вектора).
Вычислив , определяют
из граничного условия, наложенного на прогиб. Пусть на окружности радиуса
прогиб равен нулю. Тогда справедливо следующее тождество:
или
где ,
-первые компоненты векторов решений
и
при
.
Если закреплен внутренний контур пластины , то всегда
.
После определения постоянных интегрирования ,
вектор состояния в произвольном сечении r вычисляют по формуле (16.15). Решения
и
можно находить либо последовательно - интегрированием уравнения (16.13), причем при определении
следует положить
, либо на ходить эти два решения параллельно интегрированием следующего дифференциального уравнения:
где - матрица частных решений;
- матрица, первый столбец которой нулевой, а вторым столбцом является вектор g.
Заметим, что вектор частного решения - это второй столбец матрицы Y. Поэтому в программных обозначениях компоненты этого вектора имеют вторым индексом цифру 2 (т.е.
). В приведенной программе интегрируется матричное уравнение (16.16).
Описание программы. Для численного расчета пластины методом начальных параметров разработана программа, позволяющая рассчитывать пластину с одной кольцевой опорой при осесимметричной поперечной нагрузке.
При расчете пластина разбивается на участки. Участки отличаются друг от друга либо различным законом изменения толщины, либо различным законом изменения поперечной силы, либо и тем и другим.
Дата публикования: 2015-01-23; Прочитано: 644 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!