Пример выполнения задания

Постановка задачи и получение расчетных формул

На рис.С5а представлена схема составной конструкции, два элемента которой (АС и СВ) соединены между собой цилиндрическим шарниромC. На конструкцию в целом наложены внешние связи в виде шарнирно-неподвижных опор A и B и внутренняя связь в шарнире С. Размеры элементов указаны в метрах.

Дано: F₁ = 10 кH, F₂ = 7,0 кН, М = 11 кН×м, q = 1 кН/м, b = 60⁰.

Определить, при каких значениях угла a опора А будет иметь наименьший вес и потенциальная энергия деформации всех опор минимальной.

Решение. Заменим равномерно распределенную нагрузку интенсивностью q равнодействующей , модуль которой Q = q×3 = 3 кН, а точка приложения находится на середине вертикального участка элемента СВ.

Расчленим заданную конструкцию на два элемента АС и СВ, представленные на рис. С5б и рис. С5в соответственно.

Рис. С5а

Рис. С5б Рис. С5в

К каждому элементу прикладываем заданные активные силы и реакции внешних и внутренней связей. Примем для обоих элементов единую систему координат О xy. В точке А элемента АС прикладываем составляющие (компоненты) и реакции внешней связи, а в точке С – составляющие и реакции внутренней связи. В точке В элемента ВС прикладываем составляющие и реакции внешней связи, а в точке С – составляющие и реакции внутренней связи ( = - , = - ).

Составим уравнения, выражающие равновесие элемента АС:

(32)

(33)

(34)

Составим уравнения равновесия элемента ВС:

(35)

(36)

(37)

С учетом, что , уравнения (32)–(37) перепишем в виде системы шести уравнений, перенеся в их правые части известные члены, считая, что значение угла a также задано:

(38)

Система уравнений (38) позволяет определить реакции связей как функции угла a ().

С целью упрощения уравнений (38) подставим числовые значения заданных постоянных величин:

(39)

Матричная форма уравнений равновесия

Правые части уравнений (39) имеют вид

(40)

Представим систему уравнений (39) в матричной форме:

(41)

где X – матрица-столбец неизвестных,

(42)

A – матрица коэффициентов при неизвестных,

(43)

матрица-столбец правых частей системы уравнений (39)

(44)

Для любого варианта задания каждый элемент этого столбца может иметь вид , где коэффициенты a, b и с зависят от заданных сил и их плеч.

В матрицах (43) и (44) целые числа являются точными, дробное число должно быть записано с относительной погрешностью, соответствующей пяти верным значащим цифрам. Такая же запись промежуточных результатов должна быть сохранена во всех вариантах задания. Она вытекает из принятой в инженерной практике точности расчетов в 5 %, что соответствует записи окончательных результатов с тремя верными значащими цифрами.

Решением матричного уравнения (41) будет

где А^-1 – матрица, обратная матрице А.

Алгоритм вычислений

Для решения матричного уравнения (41) в Mathcad имеется функция lsolve (A, B), которая возвращает вектор Х для системы уравнений А×Х = В при заданной матрице коэффициентов А и векторе свободных членов В. В нашей задаче матрица В является переменной, зависящей от угла a. Выбираем шаг изменения угла Da = p/12.

Модули реакций в точках определяются через составляющие

Из графика R _i (a) находим значение угла, при котором данная реакция имеет минимум.

Результаты вычислений

Значения модулей составляющих реакций всех связей полной искомой реакции должны быть сведены в табл. С5.1.

Таблица С5.1

Обозначения, принятые в расчете	RA x	RA y	RC x	RC y	RB x	RB y	RA
Их иденти-фикаторы	X(1)	X(2)	X(3)	X(4)	X(5)	X(6)	F	П
a = 0
a = p/12
a = p/6
.	.	.	.	.	.	.
2p

По данным табл. С5.1 вначале строим график R_A в функции от угла a (рис. C5г). Для этого наносим на координатную плоскость R_A – a точки, соответствующие вычисленным значениям реакции. Полученные точки соединяем плавной кривой. На этой кривой отмечаем точку, соответствующую значению глобального минимума реакции R_Amin и по ней находим значение a_min, т. е. угла a, соответствующего значению R_Amin.

Рис. С5г

Рис. С5д

Примечание. В тех вариантах задания, где в точке A расположена шарнирно-подвижная опора или невесомый стержень, кривая зависимости R_A = R_A(a) может пересекать ось a и часть ее располагаться ниже этой оси (рис. С5д). Это означает, что для этой части кривой значения модулей оказались отрицательными. Эта часть графика показана пунктиром.

Поскольку модуль не может быть отрицательным, нижнюю часть графика следует зеркально перенести вверх относительно оси a. Точки пересечения кривой (сплошной и пунктирной) с осью a дают значения углов (a_min)₁, (a_min)₂, соответствующих значению реакции R_A, равное нулю (глобальный минимум).

Для оценки оптимальности конструкции по второму критерию строим график изменения сумм квадратов модулей составляющих реакции в зависимости от угла a (рис. С5е).

Рис. С5е

На графике отмечаем точку A, соответствующую глобальному минимуму целевой функции f (a) и определяем угол a_min.

Содержание отчета

Отчет по работе должен содержать:

1) схему конструкции, исходные данные и постановку задачи;

2) изложение расчета с расчетными схемами и краткими пояснениями (уравнениями равновесия, решение уравнений в общем виде, запись целевых функций, методику нахождения оптимального угла a);

3) схему алгоритма, программу;

4) распечатку числовых результатов, графики, анализ результатов.

Вычисляем потенциальную энергию деформации:

П (a) = (Х(a)₁)² + (Х(a)₂)² + (Х(a)₃)² + (Х(a)₄)² + (Х(a)₅)² + (Х(a)₆)².

Таблица C5.2

Рис. С5ж

Рис. С5з

Табл. С5.2 содержит результаты вычислений, а на рис. С5ж и С5з представлены графики зависимостей R_A = R_A(a) и П = П(a), построенные по этим результатам.

Поиск по сайту: