Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

133992

.pdf
Скачиваний:
2
Добавлен:
15.11.2022
Размер:
2.93 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

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

X 2 ,…,
X 2 ,…,

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]