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

Новакович В.И. Математическое моделирование систем и процессов. Учеб. пособ. 2019

.pdf
Скачиваний:
0
Добавлен:
15.03.2026
Размер:
3.23 Mб
Скачать

Наибольшие напряжения

 

Mz

 

M

,

(2.3)

J

W

 

 

 

 

где W – момент сопротивления сечения стержня, J Wz .

Одной из типичных для элементов верхнего строения пути является математическая модель для расчета стержня с заделанным концом (консоль), изображенная на рис. 2.2.

Рис. 2.2. Математическая модель стержня с заделанным концом

Изгибающий момент в любом сечении консоли будет равен:

M P l x

Соответственно напряжения в кромках

.

(2.4)

 

консоли можно определить по

(2.3).

В расчетах может быть использовано понятие «балки равного сопротивления», когда в любом сечении:

σM W

const

.

(2.5)

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

W bh62 , где b – ширина, h – высота сечения (рис. 2.3).

Рис. 2.3. Прямоугольное сечение стержня

В консоли с условием (2.5) прогиб на ее конце можно определить по следующей формуле:

21

 

2Pl

3

 

f

,

 

3EJ

 

 

0

 

 

 

 

где J0 – момент инерции в плоскости закрепления консоли,

J

Кстати, по (2.6) можно определить, что прогиб на конце два раза большим, чем для стержня с одинаковым сечением a

(2.6)

 

bh

3

 

 

 

.

12

 

 

консоли будет в b .

2.2 Типичные математические модели, применяемые для расчета рельсов и шпал

Рельс – это упругий стержень на упругом основании (рис. 2.4).

Рис. 2.4. Упругий стержень на упругом основании

Уравнение изогнутой оси стержня будет:

M EJ

d

2

y

 

dx

2

 

.

(2.7)

Если при распределенной вдоль стержня нагрузке p, то в соседнем сечении, отстоящем от перерезывающей силы на dx:

Q dQ Q pdx , тогда

dQ

dx

 

Дифференцируя два раза уравнение изогнутой оси

p .

(2.7) с учетом

(2.8)

 

 

2

 

 

 

dM

Q и

d

M

dQ p ,

 

 

 

2

dx

dx

dx

 

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

EJ

d 4 y

p Uy ,

(2.9)

dx4

 

 

 

где U – модуль упругости основания под стержнем. При обозначении

22

 

 

K

U

.

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4EJ

 

 

 

 

Общий интеграл дифференциального уравнения (2.9) будет:

y C e

kx

coskx C e

kx

sin kx

C e

kx

coskx C e

kx

sin kx .

 

 

 

 

1

 

2

 

 

 

3

 

4

 

 

Если, как обычно, допускается, что стержень очень длинный, то очень большим числом, тогда при y x 0 .

Этому условию удовлетворяют в общем интеграле

e

kx

 

(2.10)

(2.11)

будет

С1 С2

0

,

в этом случае уравнение изогнутой оси стержня запишется так: y e kx C3 cos kx C4 sin kx .

(2.12)

(2.13)

Для определения постоянных С3 и С4 используем условия в точке приложения груза. Здесь в силу симметрии

dy

0

,

(2.14)

dx

x 0

 

 

 

 

 

а перерезывающая сила в этой же точке равна:

 

 

 

 

 

 

 

 

d

3

y

 

 

 

Р

 

 

 

 

 

 

 

 

 

EJ

 

 

 

 

.

 

 

 

 

 

 

 

 

dx

3

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kx

C3

cos kx sin kx C4 cos kx sin kx

,

Так как y e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то из условия (2.14) следует:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С3

С4

С

,

 

 

 

Тогда y Сe

kx

cos kx sin kx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее дифференцированием получаем:

 

 

 

 

 

 

 

 

 

y 2Сke

kx

cos kx sin kx ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

2Сk

2

e

kx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin kx cos kx , .

 

 

 

 

 

 

y 4Сk

 

e

kx

cos kx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.15)

(2.16)

(2.17)

(2.18)

Пользуясь выражением для третьей производной, можем на основании

условия (2.15) записать

 

 

 

 

 

C

P

 

1

.

(2.19)

 

 

 

8EJ k 3

 

Для упругой линии изогнутого стержня имеем:

y

P 1

e kx cos kx sin kx

P

η.

(2.20)

 

 

 

 

8EJ k 3

8EJk 3

 

 

 

 

23

Тогда изгибающий момент в любом сечении стержня:

y

P

e

kx

sin kx cos kx

P

μ .

 

 

4k

 

4k

 

 

 

 

 

 

 

 

Наибольший прогиб в сечении x = 0 будет:

 

 

 

 

 

 

f

P

 

.

 

 

 

 

 

 

8EJk

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для ординат

 

 

 

 

 

 

 

 

 

 

η cos kx sin kx ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ sin kx cos kx .

 

 

 

 

 

 

 

 

 

 

 

 

(2.21)

(2.22)

(2.23)

называемых линиями влияния (или инфлюентами), для конкретных условий составлены таблицы их значений, зависящих от коэффициента относительной жесткости k [13].

Линии влияния позволяют получить результат решения по (2.20) или (2.21) сразу от нескольких внешних сил.

2.3 Математическая модель изогнутых коротких стержней на упругом основании

Таким стержнем (балкой) можно считать железнодорожную шпалу (рис. 2.5).

Рис. 2.5. Математическая модель изогнутых коротких стержней на упругом основании

Здесь необходимо рассматривать два участка упругой линии прогибов от x l2 до x l2 a . На каждом участке имеет место общий интеграл (2.11), но

произвольные постоянные для каждого участка должны быть определены отдельно. Нужно составить восемь уравнений для нахождения восьми произвольных постоянных. На концах шпалы изгибающие моменты и перерезывающие силы равны нулю, следовательно:

24

y 0 и

y 0 при

x

l

2

.

По середине шпалы в силу симметрии нагрузки получим еще два уравне-

ния:

y

 

0 и

y

 

0 при

x 0 .

 

 

Остальные четыре уравнения могут быть составлены на основании условий в месте сопряжения двух участков. Очевидно, что оба участка упругого изгиба имеют: 1) одинаковый прогиб y; 2) общую касательную, следовательно, одинаковые значения y ; 3) одинаковые значения изгибающего момента, т. е. y

Кроме того, перерезывающая сила при переходе от одного участка к другому

сразу меняется на Р, т.е.

/ y

y /

при x l 2 a

 

11

1

 

P EJ

.

2.4Математические модели состояния упругих стержней

вупругой среде

Механические модели под действием продольных сжимающих сил изогнутой оси упругих стержней могут быть представлены в следующем виде (рис. 2.6).

Рис. 2.6. Математические модели состояния упругих стержней в упругой среде

Эти модели известны, как «задачи Эйлера» – задачи на устойчивость продольно сжатых стержней. В этих задачах все факторы, влияющие на напряженнодеформированное состояние, являются детерминированными, а именно – с достаточной точностью известны их длина l и изгибная жесткость в плоскости изгиба EJ. Продольная устойчивость упругих стержней под действием критической силы Fкр по Эйлеру определяется по следующей формуле

F

μ

EJ

,

 

2

кр

 

l

 

 

 

 

 

где μ – коэффициент, зависящий от конфигурации изгиба,

n – число полуволн искривления;

(2.24)

μ π2 2n 1 2 ;

4

25

Fкр – наименьшее значение продольной сжимающей силы, при которой становится возможным искривление стержня.

Для стержня (рис. 2.6, а) n = 0

 

 

2

EJ

 

 

 

 

F

 

π

·

 

(2.25)

 

 

2

 

кр

 

4l

 

 

 

 

 

 

 

 

 

 

 

Для стержня (рис. 2.6, б) n = 1

 

 

 

 

 

 

 

 

F

 

2 EJ

·

(2.26)

 

l2

 

 

кр

 

 

 

 

 

 

Для стержня (рис. 2.6, в) n = 1,5

 

 

2

EJ

F

 

 

 

2

кр

 

l

 

 

 

 

Для стержня (рис. 2.6, г) n = 2,35

 

 

 

 

F

 

2 EJ

кр

 

l2

·

·

(2.27)

(2.28)

2.5 Математическая модель упруго изогнутого стержня под действием продольной силы в упругой среде

Расчетная схема его представлена на рис. 2.7.

Рис. 2.7. Математическая модель упруго изогнутого стержня

Его дифференциальное уравнением будем:

EJ

4 y

F

2 y

Uy 0

 

x4

x2

 

 

 

.

(2.29)

Решение (2.29) быстрее и проще получить энергетическим методом, основанном на принципе Лагранжа: «Сумма работ всех внешних и внутренних сил на возможных бесконечно малых перемещениях равна нулю», т. е.

dAi 0 .

(2.30)

Для математической модели для схемы рис. 2.6 элементарная работа продольной силы

26

 

 

F

l

 

 

2

 

 

 

 

 

dA1

 

 

 

dx ,

2

y

 

 

 

0

 

 

 

 

элементарная работа внутренних сил изгиба стержня

(2.31)

 

 

EJ

l

 

y

 

2

 

,

 

 

 

 

2

 

 

 

dx

dA

2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

элементарная работа реактивных погонных сил сопротивления

dA3 U l y2dx .

2 0

При подстановке (2.31), (2.23) и (2.33) в (2.30) и учитывая, что

(2.32)

(2.33)

y am sin

mπx

,

l

 

 

где m – число полуволн; l – длина полуволны,

(2.34)

получим

F

2 4

2

 

 

EJ am

 

 

 

 

U am

 

 

 

 

 

 

 

l

 

 

 

 

 

.

(2.35)

2

2

 

 

am

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

Критическая продольная сила, при которой произойдет искривление упругого стержня, равна:

 

 

EJ mπ

4

Ul

4

 

 

 

 

F

 

 

 

 

2

 

2

 

2

 

кр

 

π

l

 

 

 

m

 

 

 

.

(2.36)

2.6Математические модели состояния упругих стержней

впластической среде

Вжелезнодорожной инфраструктуре в качестве упругого стержня чаще всего выступает рельс, а средой – его подрельсовое основание, состоящее из скреплений шпал и балласта. Считается, что скрепления со шпалами соединены абсолютно жестко. Если допустить, что возможны продольные перемещения рельсов относительно шпал – угон, то это значит допустить аварийное состояние верхнего строения железнодорожного пути. Таким образом, в расчетах следует исходить из того, что возможны продольные перемещения рельса только вместе со шпалами, которым препятствует сопротивление балласта. Это сопротивление балласта сдвигу в нем шпал, прикрепленных к рельсу, считаем пластическим (или фрикционным), если воздействие поездов не учитывается.

Достаточно много примеров, когда воздействие поездов не возникает. Например, когда производятся путевые работы в период прекращения движения поездов, т.е. в так называемое «оконное» время. В «оконное» время выполняются

27

путевые работы с применением выправочно-подбивочно-рихтовочных (ВПР) машин, при выполнении сварочных работ с применением передвижных рельсосварочных машин (ПРСМ) или с применением оборудования для алюминотермитной сварки используют гидравлические натяжные устройства (ГМУ) или нагревательные приборы. Во всех перечисленных случаях в рельсах возникают продольные силы, под действием которых происходят продольные, а иногда и поперечные перемещения отдельных рельсов или всей рельсошпальной решетки.

Обычно в таких случаях для вывода расчетных формул рационально воспользоваться энергетическим методом. В этом методе задаются конфигурацией кривой изогнутой оси упругого стержня. Чаще всего этот изгиб соответствует «одностороннему изгибу», который можно отобразить смещенной синусоидой

y

f 2

 

1

cos

x

 

l

 

 

 

 

.

(2.37)

Из математического анализа известно, что разница λ для кривой (2.37) и стягивающей ее хордой может быть найдена из выражения:

l /2

2

2

(y') dx

0

 

,

(2.38)

После подстановки

y

πfx

sin

x

l

l

 

в (2.38) и интегрирования в пределах

стягивающей кривую хорды получим:

 

π

2

f

2

λ

 

.

 

4l

 

 

 

(2.39)

Если, например, необходимо изогнуть упругий стержень силой Р на стрелу f, на длине l в пластической среде (рис. 2.8), в которой есть погонные сопротивления сдвигу стержня поперек оси q и вдоль оси r, один конец стержня заделан, а другой имеет длину L, то для того, чтобы воспользоваться принципом Лагранжа (2.30), нужно определить величину элементарных работ dAi .

Рис. 2.8. Упругий стержень силой в пластической среде

28

Стержень считаем нерастяжимым (такое условие принимается, когда под действием продольных сил сжатием или растяжением можно пренебречь).

Работа внешней силы Р на бесконечно малом перемещении df:

dA

Рdf

1

 

.

(2.40)

Работа распределенных сил сопротивления поперечным перемещениям q на длине изогнутого участка l:

dA

ql

df

 

2

2

 

 

 

.

(2.41)

Работа сил погонного сопротивления продольным перемещениям r участка стержня на длине L:

dA3 rLdλ

Работа внутренних сил изгиба стержня:

.

(2.42)

 

 

4

EJf

dA

 

 

3

3

 

l

 

 

 

df

.

(2.43)

Для того чтобы установить связь между бесконечно малыми перемещениями вдоль (dλ) и поперек (dr), продифференцируем (2.39):

π2 f

2

df .

(2.44)

2l

 

 

 

 

 

Тогда при подстановке (2.40)–(2.43) в (2.30) с учетом (2.44) получим:

 

4

EJf

 

2

rLf

P

 

π

l

3

2l

 

 

 

 

 

ql 2

.

(2.45)

Здесь одна сила P – активная, остальные реактивные.

Рассмотрим другой вариант, например, под действием продольной силы выпрямляется предварительно изогнутый участок упругого (на изгиб), но нерастяжимого продольной силой стержня (рис. 2.9).

Рис. 2.9. Расчетная схема для определения силы F

 

Работа внешней силы Р на бесконечно малом перемещении dλ:

 

dA1 Рdλ .

(2.44)

29

Поскольку начальная стрела f0, то выпрямление происходит от стрелы f0 по переменной φ f0 f . Связь между продольными и поперечными перемещени-

ями определяется следующей зависимостью:

 

π

2

2 f

 

φ φ

2

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

0

 

 

 

 

 

 

 

 

4l

 

 

 

 

 

 

 

 

 

 

тогда

π

 

 

f

 

φ dφ

dλ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2l

 

 

,

.

(2.45)

(2.46)

Работа распределенных сил сопротивления поперечным перемещениям q на длине изогнутого участка l:

dA

ql

.

(2.47)

 

2

2

 

 

 

 

 

Работа внутренних сил изгиба рельсов будет иметь тот же знак «‒» как и в (2.44), поскольку упруго изогнутый стержень стремится выпрямиться под действием потенциальной энергии, возникающей в нем при изгибе

 

4

EJ f

 

dA

0

 

 

 

 

 

 

 

3

 

l

3

 

 

 

 

φ dφ

.

(2.48)

Работа сил погонного сопротивления участке L определяется по:

dA4 rLdλ

продольным перемещениям r на

.

(2.49)

Тогда с учетом (2.46) при подстановке в (2.30), (2.44), (2.47), (2.48) и (2.49)

получим:

 

ql

2

 

2

EJ

P

 

 

 

πf

l

2

 

 

 

 

 

 

rL

.

(2.50)

Окончательное выпрямление предварительно изогнутого участка со стрелой f0 силой F невозможно, так как при f 0, F . Тогда для выпрямления его

необходимо приложить поперечную силу G по схеме (рис. 2.10).

Рис. 2.10. Расчетная схема для определения силы G

30