133992
.pdf
|
|
|
|
|
|
|
|
|
|
|
m Li |
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
Ω |
|
y |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
M |
ij |
iP |
dx |
|
|
m |
iP |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
jP = ∑ ∫ |
|
|
|
|
|
|
|
|
|
= ∑ |
|
|
|
|
|
|
|
Cij |
(11) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EJ i |
|
|
|
|
|
|
|
|
EJ i |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
i =1 0 |
|
|
|
|
|
|
|
|
|
|
i =1 |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
где Ω iP - площадь эпюры M P на |
|
|
i - м участке; yCij - орди- |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ната эпюры M j |
|
на i - м участке, взятая под центром тяжести |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
грузовой эпюры M P . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
При вычислении единичных и грузовых перемещений |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
можно воспользоваться так же матричным способом, основан- |
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ном на формуле Симпсона. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
Рис.11 |
|
|
|
Преобразуем формулу (9) к матричному виду |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
m Li |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
M |
ij |
M |
ik |
dx |
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
δ jk = ∑ ∫ |
|
|
|
|
|
|
|
|
|
∑ M ijT Li M ik , |
12) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
EJ i |
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
i =1 0 |
|
|
|
|
|
|
|
|
|
|
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Л |
1 0 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л |
|
||||||||||||||
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ij |
|
|
|
|
|
|
|
|
|
|
Li |
|
|
|
|
|
|
|
|
|
|
ik |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
M ij = M ijC ; |
Li = 0 4 0 |
|
|
|
; M ik |
= M ikC |
, |
||||||||||||||||||||||||||||||||
|
|
|
|
|
EJ i |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
П |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
П |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
M ij |
0 0 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M ik |
|||||||||||||||||||||
|
|
T = {M Л , M C , M П }- транспонированная матрица |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
M |
M |
ij |
. |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
ij |
ij |
|
|
ij |
|
ij |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
Дальнейшее упрощение формулы (12) сводится к объединению векторов и мат- |
||||||||||||||||||||||||||||||||||||||||
риц по всем участкам упругой системы, |
так что |
|
окончательно эта формула принимает |
||||||||||||||||||||||||||||||||||||||||
вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
m Li |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
M |
|
M |
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
δ jk = ∑ ∫ |
|
|
|
ij |
|
|
|
ik |
|
|
|
M Tj |
L M k . |
|
|
(13) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EJ i |
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
i =1 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
Канонические уравнения метода сил (5) удобно решать в системе MathCAD, |
||||||||||||||||||||||||||||||||||||||||
предварительно записав их в матричной форме |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δ X = |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(14) |
|||||||||
где δ |
-матрица податливости (матрица единичных перемещений); |
|
X - вектор неиз- |
вестных реакций связей; - вектор грузовых перемещений взятых с обратным знаком
|
δ |
11 |
δ |
12 |
K δ |
1n |
|
|
X |
1 |
|
− |
1 P |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
δ 21 |
δ 22 |
K δ 2 n |
|
|
|
||||||||||
|
X |
2 |
|
− |
2 P |
||||||||||
δ = |
|
|
|
|
|
|
|
|
; |
X = |
|
; |
= |
|
. |
|
K K K K |
|
|
K |
|
|
K |
|
|||||||
|
δ n1 |
δ n2 |
K |
|
|
|
|
|
|
|
|
|
|
||
|
|
δ nn |
|
X n |
− |
nP |
Решением системы уравнений является вектор |
|
X = δ −1 , |
(15) |
где δ −1 - обратная матрица для матрицы податливости δ .
11
3.СОДЕРЖАНИЕ ЗАДАНИЯ
«РАСЧЕТ СТАТИЧЕСКИ НЕОПРЕДЕЛИМОЙ РАМЫ МЕТОДОМ СИЛ»3
o Требуется:
1.Раскрыть статическую неопределимость заданной рамы, учитывая деформа- цию участков только от изгибающего момента.
2.Построить эпюру изгибающих моментов.
3.Определить перемещения указанной в задании точки рамы.
Примечание
Вычисления выполнить в системах SCAD и MathCAD. Результаты реше- ния, полученные в системе SCAD, использовать как контрольные для проверки правильности решения задачи в системе MathCAD.
o |
Заданы: |
1. |
Схема плоской статически неопределимой рамы, размеры ее элементов, зна- |
чение и расположение внешних сил. |
|
2. |
Модуль продольной упругости материала E и осевые моменты инерции J |
площади поперечных сечений элементов рамы. |
|
o |
Варианты заданий |
Необходимые для расчета данные приведены в табл.1. Во всех заданных рамах модуль продольной упругости материала стержней E , осевой момент инерции пло- щади поперечного сечения стержней J , а также величины, указанные в табл.1, счи- тать известными: a = 2 м, q = 10 кН/м, q1 = 5 кН/м. Искомые перемещения точек системы отмечены на схемах (табл.2) вопросительным знаком.
|
|
|
|
|
Таблица1 |
|
№ |
Исходные данные |
№ |
Исходные данные |
№ |
Исходные данные |
|
п/п |
п/п |
п/п |
||||
|
|
|
1 |
q1 , P = q1a , L = a 9 |
q1 , L = 2a |
17 |
q , P1 = qa , L = a |
||||
2 |
q1 , L = 2a |
10 |
q1 , P = q1a , L = a 18 |
q , P = qa , L = 2a |
||||
3 |
q , P = qa , L = a |
11 |
q , L = 2a |
19 |
M , a |
|
||
|
|
|
|
|
|
|
|
|
4 |
q , P = qa , a |
12 |
q , |
M = q a 2 , |
20 |
q , M = qa |
2 |
, L = a |
1 |
1 |
|
||||||
|
|
|
|
|
|
|
||
|
|
|
|
L = a |
|
|
|
|
5 |
q , L = a |
13 |
q1 , P = q1a , L = 2a |
21 |
q , P = qa , a |
|||
6 |
q , L = a |
14 |
q , P1 = qa , L = a |
22 |
q , P , L = a |
|||
7 |
q , P = qa , L = 2a |
15 |
q1 , P = q1a , L = 2a |
23 |
q , P , a , ϕ = 30 o |
|||
8 |
q , P = qa , L = a |
16 |
q , M = qa 2 , L = a |
24 |
q , P = qa , a |
|||
25 |
q , P = qa , L = 2a 26 |
q , P = qa , L = 2a |
|
|
|
|
3 Расчетные и курсовые работы по сопротивлению материалов /Ф.З. Алмаметов, С.И. Арсеньев, С.А. Енгалычев и др. – М.: Высш. шк., 1992.- 320 с.
12
Таблица 2
13
14
15
4.МЕТОДИЧЕСКИЕ УКАЗАНИЯ
4.1.Последовательность расчета статически неопределимых систем по методу сил
4.1.1. Выбор основной системы. Основная система должна быть статически определимой и геометрически неизменяемой, т.е. связей должно быть отброшено не больше, чем статическая неопределимость n . При выборе рабочего варианта основной системы руководствуются также соображениями простоты расчетов, включая вычис- ление единичных и грузовых перемещений, решение системы канонических уравнений метода сил и т.д.
За лишние неизвестные при этом могут быть приняты как внутренние силовые факторы в некоторых сечениях, так и реакции опор.
4.1.2.Построение единичных и грузовых эпюр. Для выбранной основной сис-
темы для каждого участка записываются аналитические выражения для внутренних силовых факторов и строятся эпюры этих внутренних силовых факторов, которые ис- пользуются при определении единичных и грузовых перемещений.
4.1.3.Вычисление коэффициентов канонических уравнений метода сил.
Вычисляются интегралы в выражениях (6) или (7). В системах, состоящих из прямоли- нейных стержней, при перемножении эпюр по методу Верещагина, площадь вычисля- ется для криволинейной эпюры, а ордината под её центром тяжести берется с линей- ной эпюры.
16
4.1.4. Решение системы канонических уравнений метода сил. Для любых правильно выбранных основных систем уравнения всегда имеют единственное реше-
ние X 1 , X n . Если в результате решения системы канонических уравнений метода сил какое-то из неизвестных X i имеет знак минус, это означает, что истинное
направление этой реакции противоположно первоначально выбранному.
4.1.5. Построение эпюр внутренних силовых факторов. После того как вы-
числены значения лишних неизвестных X 1 , X n могут быть определены внут- ренние силовые факторы в любых сечениях стержней системы и построены их эпюры. Например, эпюру изгибающих моментов M y в статически неопределимой системе
можно построить следующим образом. Построим эпюры изгибающих моментов от ка- ждого из найденных усилий X 1 , X 2 ,…, X n . Для этого можно использовать постро-
енные ранее единичные эпюры M i . Просуммировав по характерным точкам ординаты эпюр от действия всех сил X i с ординатами грузовой эпюры M P , т.е.
n
M y = M P + ∑X i M i ,
i =1
получим окончательную эпюру M y для заданной статически неопределимой системы. Другой способ построения эпюры M y заключается просто в применении обыч- ного метода сечений к основной системе, загруженной внешними силами и силами
X1 , X 2 ,…, X n .
4.1.6.Деформационная проверка. После раскрытия статической неопредели-
мости перемещения k какого-либо сечения вычисляются с помощью интеграла Мак- свелла-Мора (6). При этом единичное воздействие прикладывается к основной систе- ме. Строится соответствующая эпюра M y и, вычислив интегралы в формуле (6), оп-
ределяют k . Это возможно в силу того, что, как уже отмечалось при формулировке содержания метода сил, основная система под действием внешней нагрузки и сил X 1 , X 2 ,…, X n деформируется так же, как исходная статически неопределимая система.
Таким способом можно определить любые перемещения в основной системе, в том числе и перемещения в направлении лишних неизвестных. Поскольку эти перемеще- ния должны равняться нулю в силу наложенных связей, то деформационная проверка и
заключается в проверке равенства нулю перемещений j в основной системе, т.е.
m Li |
M |
ij |
M |
iP |
dx |
|
||
jP = ∑ ∫ |
|
|
|
= 0 , (j = 1,2,3Kn). |
(16) |
|||
|
|
EJ i |
|
|
||||
i=1 0 |
|
|
|
|
|
Деформационная проверка контролирует правильность расчета статически не- определимой системы от этапа определения единичных и грузовых перемещений до
построения окончательной эпюры изгибающего момента M y . Она проверяет также
правильность решения системы канонических уравнений метода сил, но не дает указа- ний о корректности выбора основной системы и правильности построения для этой системы единичных и грузовых эпюр.
17
4.2. Пример расчета
Дано (рис. 12, а):
oгеометрические размеры рамы: R = 2 м; h = 3 м ;
oжесткость поперечного сечения стержней рамы: EJ = const = 1 ;
oнагрузка: P = 50 кН; q = 10 кН/м.
Рис.12
Требуется:
1.Раскрыть статическую неопределимость упругой системы методом сил, учиты- вая деформацию участков только от изгибающего момента, и проверить пра- вильность определения дополнительных неизвестных.
2.Построить эпюры нормальных и поперечных сил, изгибающих моментов.
3. Определить линейное перемещение Dx точки D .
4.2.1 Расчет статически неопределимой рамы в системе SCAD
Система SCAD позволяет сразу рассчитать любую раму (статически определи- мую или статически неопределимую), загруженную заданной нагрузкой, но мы на эту систему возложили функцию преподавателю – проверку результатов поэтапного стан- дартного расчета статически неопределимых систем методом сил. Поэтому приводи- мая ниже методика расчета не соответствует традиционной методике прочностных расчетов в системе SCAD, а соответствует процедуре создания «шпаргалки» для про- верки результатов расчета в системе MathCAD.
Прежде чем мы начнем осваивать процедуру расчета упругих статически неоп- ределимых систем в программном продукте SCAD, давайте договоримся о методике создания конечноэлементной модели рамы.
Заданная для расчета рама имеет криволинейный участок, положение которого удобно описать в полярной системе координат. Так как в SCAD применяются только правые системы координат, то положительные дуговые координаты отсчитываются в направлении против хода часовой стрелки. А так как в МКЭ необходимо соблюдать определенный порядок обхода узлов, так чтобы местные системы координат по отно- шению к элементам имели одинаковую ориентацию, ибо усилия и напряжения выво- дятся в местной системе координат, то целесообразно ввод узлов и элементов начинать
18
от точки B по направлению к точке A (рис.12). Кроме того, задание распределенной нагрузки на криволинейных участках рамы легче реализуется в местной системе коор- динат.
Стержневой конечный элемент универсального вида изображен на рис.13. При- нято следующее обозначение концевых сечений: номером 1 обозначается сечение у начала стержня, номером 2 - у его конца.
Рис.13
В общем случае предполагается, что стержень может быть произвольным обра- зом расположен по отношению к принятой общей (глобальной) системе координат XYZ . В местной (локальной) системе координат X 1Y1 Z1 , относительно которой вы-
числяются усилия, оси Y1 и Z1 , являются главными осями инерции поперечного сече-
ния.
Еще до ввода информации о форме поперечного сечения система SCAD строит оси локальной системы координат X 10Y10 Z10 («предварительный проект») с использова- ние следующих соображений:
•ось X 10 направлена вдоль геометрической оси упругой части сечения от начала стержня к его концу;
•ось Z10 ориентирована в верхнее полупространство (в сторону возрастания коорди-
наты Z общей системы координат XYZ , в которой представлена расчетная схема;
• ось Y10 дополняет локальные координаты таким образом, что система образует пра- вую тройку X 10Y10 Z10 .
Переходим к описанию алгоритма создания конечноэлементной модели рамы. Приведенная на рис.12, а рама три раза статически неопределима внутренним
образом. За «лишние» неизвестные выбираем внутренние силовые факторы: N − X 1 , Q − X 2 , M − X 3 . Эквивалентная система приведена на рис.12, б.
Сначала создадим конечно-элементную модель основной системы.
1.Загружаем комплекс Structure CAD.
2.Создаём новый проект.
3.Задаём единицы измерения: линейные размеры – м; размеры сечения – мм; силы
- кН.
4.Активизируем раздел Расчетная схема и щелкаем по кнопкам Узлы и Ввод узлов .
5.Задаём координаты узлов, соответствующих характерным точкам основной системы, начиная с узла B (рис.12): (4 ,0 ,0 ) - узел 1 , (4 ,0 ,3) - узел 2 ,
19
(2.05 ,0 ,3) - узел 3 , (2 ,0 ,3) - узел 4 , (1.95,0 ,3)- узел 5 , (0 ,0 ,3) - узел 6 ,
(0 ,0 ,0 ) - узел 7 . Включаем фильтры Узлы и Номера узлов . Окно отображения графической информации будет выглядеть так, как показано на
рис.14.
6.Отжимаем кнопку Узлы и нажимаем кнопку Элементы . Щелкая по первому, а затем по второму узлу вводим первый стержневой элемент.
7.Активизируем кнопку Добавление стержней .
Рис.14
8.Нажимаем на кнопку Разбивка стержня . Открывается одноимённая па- нель, в поле ввода которой заносим цифру 15 и щелкаем по клавише ОК.
9.В графическом окне выделяем стержень 1 (он изменит черный цвет на крас- ный) и в главном меню, нажимая кнопку , подтверждаем наш выбор. В ре- зультате стержень будет разбит на пятнадцать стержневых элементов.
10.Активизируем кнопку Ввод элементов по дуге . Появится панель Ввод эле- ментов по дуге окружности. Вид панели с заполненными полями ввода приве- ден на рис.15. После ввода информации щёлкаем по кнопке ОК.
11.В поле отображения графической информации выделяем 4 узел и, нажимая кнопку , подтверждаем наш выбор. На дисплее появится изображение ко- нечноэлементной модели, двух участков рамы (рис.16, а). Если включить фильтр Местные оси элементов , отключить Номера элементов и
включить фильтр Номера узлов , то можно видеть, что местные оси элемен- тов, как на стойке, так и на криволинейной части рамы направлены одинаково – ось Z1 направлена наружу рамы (рис.16, б и в).
20