Материал: 1973

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Решение систем (2.4), (2.5), (2.6), (2.7) определяет внутренние силы, перемещения и угол поворота сечения во всем диапазоне текущей координаты x при любых заданных законах изменения внешних нагрузок и размеров сечений при условии приложения нагрузок в изгибе с учетом больших перемещений представляется в конечноразностном виде

Ni+1 = Ni – qxi x,

Qzi+1 = Qzi – qzi x,

M i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

(N

 

sin

 

Q

 

cos

) (1

Ni cos i Qzi sin i

) x,

i

i

i

zi

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

EAi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

V (

Ni

cos i

Qzi sin i

) sin

 

x,

 

 

 

 

 

 

 

i

i 1

 

i

 

 

 

 

 

 

 

EAi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M i

x,

 

 

 

 

 

 

 

 

 

 

 

i 1

i

EI i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

W ((1

Ni cos i Qyi sin i

) cos

 

1) x.

 

i

 

i 1

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EAi

Для случая изгиба при малых перемещениях в виде

Ni+1 = Ni – qxi x,

Qzi+1 = Qzi – qzi x,

Мi+1 = Мi + (Qzi + Ni i) x, Vi+1 = Vi - i x,

i+1 = i + (Mi / EIi) x,

Wi+1 = Wi +((Ni- i ·Qzi)/(EAi)) x.

Если рассматривать случай поперечного изгиба, то система (2.9) упрощается

Qzi+1 = Qzi – qzi x,

Мi+1 = Мi + Qzi x,

Vi+1 = Vi - i x,

i+1 = i + (Mi /EIi) x.

Если же рассматривать случай растяжения сжатия, то система (2.9) упрощается и принимает вид

Ni+1 = Ni – qxi x,

Wi+1 = Wi +((Ni)/(EAi)) x.

(2.8)

(2.9)

(2.10)

(2.11)

26

Такая форма записи дифференциальных уравнений позволяет в электронной среде MS Excel пошагово интегрировать представленные системы и наглядно иллюстрировать результаты. Использование стандартной процедуры Excel «Поиск решения» позволяет с заданной точностью определять неизвестные начальные параметры No, Qo, Mo, Wo, Vo, o, соответствующие стандартному набору граничных условий, сформулированных в виде целевой функции. Поиск производится методом Ньютона или сопряженных градиентов. Процедура «Поиск решения» адаптируется для решения различного рода проектных и оптимизационных задач, например определение напряженнодеформированного состояния стержневых элементов, поиск формы распределения сечения по длине стержня, доставляющей равные напряжения по длине, оптимизация стержня по параметру минимального веса стержня.

Пример 1. Расчет центрально растянутого стержня методом начальных параметров

Для стального стержня сечением bxh = 5x5 см пролетом Lx=3 м построить эпюры продольных усилий N и горизонтальных перемещений W. Стержень жестко защемлен в начале, конец свободный. Стержень растянут равномерной погонной нагрузкой qx = 200 кгс/м.

Для получения результатов, с достаточной для инженерных расчетов точностью, разобьем стержень на отрезки с шагом x =Lx/100.

Рис.2.3. Расчетная схема стержня

Решение

Запишем все исходные данные в ячейках листа MS Excel, как показано в табл. 2.1.

27

Таблица 2.1

Исходные данные

L

300

см

ΔL

3

 

см

qx

2

 

кгс/см

b

5

 

см

h

5

 

см

A

25

см2

E

2 100

000

кгс/см2

Ry

2 450

кгс/см2

Для решения задачи на растяжение-сжатие воспользуемся системой уравнений (2.25). Для удобной записи данной математической модели потребуется 6 столбцов на листе Excel (табл. 2.2).

Таблица 2.2

Наименование столбцов таблицы

№т п/п L, см qx, кгс/см P, кгс N, кгс W, см

Где №т п/п – номер точки по порядку, поскольку стержень разбиваем на сто отрезков, то получаем сто одну точку на стержне;

L, см – координата соответствующей точки;

qx, кгс/см – распределенная погонная нагрузка на стержень; Р, кгс – сосредоточенная нагрузка в точке стержня;

N, кгс – продольное усилие в стержне;

W, см – горизонтальное смещение точек стержня.

Для удобства расчетов начало стержня помещаем в начало координат. Таким образом у точки №1, которая располагается в начале координат, будет координата L1 = 0. Координаты последующих точек будут записаны следующим образом:

Li+1 = Li+ ΔL,

то есть, координата последующей точки получается прибавлением длины отрезка ΔL к координате предыдущей точки. При такой записи при изменении длины координаты точек будут рассчитываться автоматически.

Поскольку распределенная погонная нагрузка приложена по всей длине стержня (рис. 2.3), то значение qx будет записано в каждой ячейке для каждой точки. Поскольку направление нагрузки совпадает

28

с направлением оси, то нагрузка записывается со знаком “+”. Если нагрузка была бы направлена против оси, то запись шла со знаком “-”.

Поскольку к стержню не приложено никакой сосредоточенной нагрузки, то во всем столбце будут стоять значения “0”. Правило знаков для сосредоточенных нагрузок аналогично распределенным.

Продольные усилия записываются в следующем порядке:

-в первой ячейке продольных усилий ставится “0”;

-во второй ячейке продольных усилий вводится формула из системы уравнений (2.11) с учетом сосредоточенных сил. То есть продольное усилие в каждой последующей точке равно значению усилия в предыдущей точке плюс изменение продольного усилия на отрезке ΔL

(рис.2.4);

-формула из второй ячейки растягивается на весь столбик.

Рис.2.4. Запись расчетной формулы для продольных усилий

Горизонтальные перемещения записываются в следующем порядке:

-в первой ячейке горизонтальных перемещений ставится “0”;

-во второй ячейке горизонтальных перемещений вводится формула из системы уравнений (2.11). То есть горизонтальные перемещения в каждой последующей точке равны значению перемещения в предыдущей точке плюс деформации отрезка ΔL (рис. 2.5);

-формула из второй ячейки растягивается на весь столбик.

29

Рис.2.5. Запись расчетной формулы для горизонтальных перемещений

Врезультате получаем взаимосвязанную систему дифференциальных уравнений.

Дальнейшее решение происходит при помощи процедуры «Поиск решения» (рис. 2.6), которая добавляется в рабочее меню через «Параметры Excel», графа «Надстройки».

Воткрывшемся окне «Поиск реше-

ния» (рис. 2.7) требуется указать изменяемую и целевую ячейки. Изменяемой ячейкой является неизвестный параметр в начале координат. Целевой ячейкой является

известное значение в конце координат. Из- Рис.2.6. Поиск решения меняемые и целевые ячейки назначаются в соответствии с граничными условиями.

Поскольку в начале координат располагается жесткая заделка, то горизонтальное перемещение на опоре равно нулю W1 = 0. Данное значение напрямую вводится в первую ячейку горизонтальных перемещений. Продольное усилие является неизвестным параметром. Таким образом, изменяемой ячейкой будет продольное усилие в начале координат N1.

Поскольку на конце стержня свободный конец, то N101 = 0, это значение и будет целевой ячейкой.

30