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

книги / Строительная механика и металлоконструкции строительных и дорожных машин

..pdf
Скачиваний:
8
Добавлен:
19.11.2023
Размер:
27.58 Mб
Скачать

t2 =-го°с

+

h - l/W

h

 

h

 

. ( t, = *10‘C

, {

~

С

 

 

Vw77.

В

V

 

 

a)

 

 

Рис, 1,34. Определение переме­ щения под действием темпера­ туры

Рис. 1.35. Определение переме­ щения от дислокаций

лить от ее действия реакции /у во всех i-x связях. При перемножении этих реакций с дислокациями берется положительный знак, если реакция и перемещение имеет одно направление.

Определим, например, вертикальное перемещение Д ^ точки к рамы (рис. 1.35, а) , основание которой одновременно переместилось вертикально на ёд

и повернулось на угол \рд. Для этого найдем

в заделке А (в связи, получившей

дислокации)

реакции от силы

= 1, приложенной по направлению искомого пе­

ремещения

(рис. 1.35, б). В

соответствии с

формулой (1.64) получим Д ^ =

= -

~ ^ 4 ) = dA +

ltpA ‘

 

1.4.4.ОСНОВЫ ВАРИАЦИОННЫХ ПРИНЦИПОВ

ИМЕТОДОВ РАСЧЕТА

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

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

Полная энергия Э (потенциал всех сил) равна работе, которую со­ вершают как внутренние, так и внешние силы при переводе (а не при разгрузке!) системы из ее деформированного состояния в начальное, недеформированное:

Э = W + Т7

(1.65)

где W — энергия упругой деформации (упругий потенциал) (W > 0); Г - энергия (потенциал) внешних сил — величина отрицательная и равная произведению обоб­ щенной силы на обобщенное перемещение без множителя 0,5 (в силу того, что при переводе системы из деформированного состояния' в недефор миров энное внеш­ ние силы остаются постоянными).

Энергию деформации W так же, как и потенциал внешних сил Г, можно представить через функции перемещения точек системы. Для этого внутренние силы в выражении (1.33) надо выразить (в соответст­ вии с формулами сопротивления материалов) через соответствующие перемещения. Например, при изгибе балки Af(z) = —EJv"(z) и, следо­ вательно,

W = 0,5 E /£ 7 (v " )J dz.

(1.66)

/

 

Таким образом, полную энергию Э системы можно представить как функцию перемещений v (z ):

Э = 3(v).

(1.67)

Отметим, что величина, зависящая не от z, а от v(z), в математике называется функционалом.

В ряде случаев энергию деформации W отдельного стержня удобно выразить не через функцию v(z) (при изгибе), а через дискретные пе­ ремещения одного или нескольких его сечений. Например, для изогну­ той балки постоянной жесткости (рис. 1.36) с учетом того, что прогиб fB = PI31(3EJ) , можно записать:

W = 0,5PfB = 0,5(3 E J/l3) fB .

Следовательно, энергию деформации W и полную энергию системы Э можно представить через функции перемещений v(z) или через не­ сколько дискретных параметров, определяющих состояние системы. Например, для рассмотренной балки этим параметром является пере­ мещение fB в точке В. Согласно (1.65)

1

*EJ

_

 

 

Э -

*2

-

pfB -

= ------------ П

В

2

 

В

Как видим, здесь полная энергия системы определяется через одно перемещение fB.

В общем же случае величина Э будет зависеть от выбора функций v(z) или от дискретных перемещений, обозначаемых в общем случае через Z,-. Перемещения Z/ называются обобщенными координатами. В стержневых системах функцию перемещений, например v(z), можно точно выразить через обобщенные координаты Z/. Потенциал всех сил системы Э = Э (Z х, Z 2 ...) в этом случае вычисляется точно.

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

Рис. 1.36. Схема для определения полной энергии

 

 

системы

 

р и

EJ

В'

1

 

 

Рис. 1.37. Схемы, поясняющие вариационные прин-

t

 

ципы

 

 

Принцип вариации перемещений (принцип Лагранжа) читается так: из всех мыслимых перемещений упругой системы истинные перемеще­ ния сообщают энергии минимальное значение.

Поясним его на примере. Выразим потенциал сил для балки, изобра­ женной на рис. 1.37, а через функцию прогибов v (z ). С учетом формулы ( 1 .66) получим

Э(v) = W(w) + T(v) = 0,5 JEJ(v")2 dz - Pf. I

Исследуем, как изменяется потенциал сил Э в зависимости от произ­ вольного изменения (вариации) изогнутой оси балки. Увеличим, напри­ мер, все ординаты v(z) в к раз (рис. 1.37, а), так что \ к = ку и fk = kf. Тогда функционал (1.67)

Эк(v) =

к 2 0 ,5 /E J(y")2 dz

- kPf.

 

 

 

/

 

 

 

Так как

потенциальная энергия упругой

деформации

W численно

равна действительной работе внешней силы А ,

 

 

эк (v) =

к 2 0,5P f - kP f =

Pf{0,5k2 -

к).

( 1 .68)

Как видим, функционал Э представляет собой функцию второй степени от коэффициента к = f k/f. На рис. 1.37, б изображена зависи­

мость Э (к) , из которой

следует, что в действительном состоянии (при

к = 1) полная

энергия

имеет отрицательное минимальное

значение:

Э (v) = min. Здесь приращение энергии равно нулю:

 

ЬЭ = 0.

 

 

(1.69)

Сказанное (для устойчивых систем, которые мы рассматриваем)

можно выразить и в такой форме:

 

d 3 /d fk = 0;

d 23 / d f 2 > 0.

(1.70)

Вариационное уравнение (1.69) используется для приближенного определения деформированного состояния сложных систем (пластин, оболочек и др.). Применение его для дискретных систем заключается

в следующем. Функцию перемещений, например v (z ), можно предста­ вить в виде совокупности п базисных функций f.( z ) :

v(z) = e i/,(z )+ a2fi(z)+

+ a„/„(z),

(1.71)

где каждая из базисных функций f^ z ) соответствует

z-му состоянию

системы и не противоречит граничным условиям, а неизвестные коэф­

фициенты

(параметры) д/ представляют собой обобщенные координаты

Z. В этом

случае функция Э будет зависеть от параметров д,- и от на­

грузки Р:

 

 

 

Э = 3 (a l f a2, . . . a n,P ).

(1.72)

Условие равенства нулю полной вариации 5 Э:

 

ЪЭ

ЪЭ

дЭ

8Э = ----- da\ +

------da2 +

+ ----- dan = О

 

Ъах

Ъаг

Ъап

при произвольных и независимых вариациях 6д/ может быть выполнено только в случае соблюдения равенств:

ЪЭ

ъм

ът

—=

ип

Ъах

Ъах

-----

Ъах

 

 

ЪЭ

ЪW

ЪТ

=

о,

----- =

+ —

 

Ъап

Ъап

систему п алгебраических уравнений, каждое

образующих линейную

из которых выражает условие равновесия, соответствующее /-му состоя­ нию системы. Решая эту систему уравнений, находят неизвестные пара­ метры Д/.

Аналогично формулируется и доказывается принцип вариации внутренних усилий (принцип Кастильяно) : энергия системы в действи­ тельном состоянии имеет минимальное значение. Для линейно-дефор- мируемых систем при использовании принципа Кастильяно сохраняются те же уравнения (1.73), только в них роль обобщенных координат Z выполняют не перемещения, а неизвестные силы - ’’силовые параметры”

Если в выражении (1.72) нагрузку Р надо определить, например, в задачах устойчивости, то из нетривиального решения системы урав­ нений (1.73) находят искомый параметр нагрузки, а параметры д/ ос­ таются неопределенными. В этом заключается применение метода Ритца для отыскания так называемых собственных значений (критических сил или частот свободных колебаний) (см. пп. 1.8.2, 2.2.3).

1.5. РАСЧЕТ СТАТИЧЕСКИ НЕОПРЕДЕЛИМЫХ

СТЕРЖНЕВЫХ СИСТЕМ

1.5.1.СТАТИЧЕСКИ НЕОПРЕДЕЛИМЫЕ СИСТЕМЫ

ИОБЩИЕ ПРЕДПОСЫЛКИ МЕТОДОВ ИХ РАСЧЕТА

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

псистемы. Одновременное разрезание (удаление) лишних связей превра­ щает систему в статически определимую./

Если лишними являются только внешние связи, то систему назы­ вают внешне статически неопределимой. Если лишние связи имеются только во внутреннем образовании свободной (отделенной от опор) системы, то она называется внутренне статически неопределимой. Такое деление условно, поскольку оно зависит от того, какие связи прини­ маются за лишние. Системы могут быть одновременно внешне и внут­ ренне статически неопределимые.

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

( Расчет статически неопределимой системы начинается с установления ее степени статической неопределимости п.\ Для плоских рам, состоящих из замкнутых контуров, удобно использовать формулу

п = Ък - ш,

 

(1.74)

где к число замкнутых контуров;

ш -

число врезанных шарниров с учетом

их кратности.

 

 

При

разрезании (’’размыкании”)

плоского замкнутого стержня

(контура)

удаляются три связи

(рис. 1.38, а) : две связи, препятствую­

щие взаимным линейным перемещениям, и одна связь, препятствующая взаимному угловому перемещению сечений контура в месте разреза. Этим обобщенным перемещениям соответствуют обобщенные внутрен­ ние силы М, QnN.

На рис. 1.38, б показана рама с двумя контурами (земля считается замыкающей связью). Шарниры А и С — однократные, а шарнир D — двукратный. Поэтому согласно выражению (1.74) п = 3 2 —4 = 2 .

Рис. 1.38. Основные системы метода сил

Рама, изображенная на рис. 1.38, в, 7 раз статически неопределима (п = = 3 3 —2 = 7). На рис. 1.38, г, д показаны системы, полученные из за­ данных после удаления лишних связей. Эти системы статически опреде­ лимы.

Отметим, что при удалении шарнира удаляются две внутренние (линейные) связи; при разрезании затяжки CD (рис. 1.38, б, г) (это делается обязательно) удаляется одна (линейная) связь; также одна связь (угловая) удаляется при врезании в контур шарнира.

Общие предпосылки методов расчета статически неопределимых систем. Расчету статически неопределимой системы подлежит не задан­ ная, а преобразованная система, которая должна быть эквивалентна ей в статическом й кинематическом отношениях —все усилия и перемещения в эквивалентной системе должны быть такими же, как и в заданной.

Условия эквивалентности заданной и преобразованной систем вы­ полняются при соблюдении принципа заменяемости связей, согласно которому:

любая /-я связь может быть разрезана (удалена), если вместо нее приложены обобщенные силы Я/ и принято, что взаимное перемещение,

соответствующее этим силам, равно нулю,т.е.

 

А/ = 0;

(1.75)

любая /-я связь может быть введена, если при этом введенная связь получит соответствующее перемещение Z/ и будет принято, что реакция (усилие) в этой связи равна нулю, т.е.

Ri = 0.

(1.76)

Система, полученная из заданной путем удаления связей или введе­ ния их, называется основной системой, а соответствующие ей неизвест­ ные Х{ nZj - основными неизвестными.

Условие кинематической эквивалентности представляет собой равенство нулю обобщенного перемещения А/ от действия на преобра­ зованную систему как внешних воздействий, так и основных неизвест-

ных X j или Z j . Обозначая для общности рассуждений основные неизвест­ ные через У/, представим это условие в виде

Д,- = ДtYi + AiYi +

+ д.

+ Д/У. + Д/m = о,

 

 

(1.77)

где п - число основных неизвестных У/, a m - фактор внешнего воздействия.

Любой член этого уравнения, за исключением последнего, можно представить в виде произведения единичного перемещения 8.к на соот­ ветствующее неизвестное Yk : А /У,. _ 5ik Yk- Следовательно, выражение

(1.77)

можно записать так:

 

 

 

А/

— 5/1 У'1 + 6/2 У2 +

^ik^k +

+

Ул + Ai m — О,

 

 

 

 

(1.78)

где Дim - обобщенное перемещение Д/ от внешнего воздействия т на основную систему.

Исходя' из аналогичных рассуждений, запишем условие статической эквивалентности (1.76):

 

“* ri\ Y \ "** ri2 У2

г{к У к

*

*

rin Уд +

0 »

 

 

 

 

 

 

(1.79)

где

- реакция в /-й введенной связи от Уд. =

1;

-

реакция в той же связи

от внешнего воздействия.

 

 

 

 

 

Число уравнений типа (1.78) и (1.79) равно числу п основных неизвестных. В совокупности они образуют систему канонических уравнений:

С1 1 У1 +

с 1 2 У2 +

+ c\k Yk + . . .

+

С\п Yn +

С\т —0;

 

с2 1 Yi +

с 22 У2 +

+ С2кУк +

+

ст Yn +

Сгт = 0 ;

(1 .8 0 )

У1 +

СП2 у 2 +

. ■ +cnkYk + * ♦ • +

спп Ул +

Спт = 0 »

 

где cik - единичные коэффициенты (б.к или г -к) , не зависящие от внешнего воз­

действия; С(т - свободные члены (Д/ш или Rim)* зависящие от внешнего воз­ действия.

Размерность единичного коэффициента с.к определяется отношением размерности свободного члена С,*т к размерности соответствующего неизвестного Yk : [с.к] = [Qm]/[Yk ]. Отметим, что в соответствии с теоремами о взаимности [см. формулы (1.43) и (1.46)] Cik = СкГ Вычислив коэффициенты cik и свободные члены Qm из решения системы п уравнений с п неизвестными (1.80), находят величины У/. Затем в соответствии с принципом суперпозиции определяют усилия

S. в любой /-й связи или перемещения Д;.любого/-го узла:

Yy

+ s-

У

+ S- у

(1.81)

 

п

]ТП7

 

У. + «/ 2 r 1 + . . . + S/„yn + д / т ,

(1.82)

Где Sjm - усилие в /-й связи;

Д^т - перемещение /-го узла ог внешнего воздейст­

вия в основной системе.

 

 

 

 

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

Если основная система получена из заданной только путем удаления связей и основными неизвестными являются силы Л/, то метод расчета называется методом сил\ если основная система получена из заданной только путем наложения связей и основными неизвестными являются перемещения Z/, то метод расчета называется методом перемещений; если же основная система получена путем одновременного и удаления, и

наложения связей и основными неизвестными являются силы

и

перемещенияZ/, то метод расчета называется смешанным.

 

Универсальный контроль решения по любому из этих методов заключается в том, что перемещение Ду по направлению любой отбро­ шенной /-й связи, вычисленное по формуле (1.82), должно быть равно нулю; также должна быть равна нулю реакция R., вычисленная по формуле (1.81) в любой введенной /-й связи. I

1.5.2._МСТОДсил

Особенности метода. Основной системой метода сил является неиз­ меняемая и статически определимая (хотя это и не обязательно) систе­ ма, полученная из заданной путем устранения всех лишних связей. От­ метим, что основных систем в методе сил множество, так как разрезать замкнутый стержень или врезать шарнир можно в различных точках. Каждой основной системе соответствуют свои основные неизвестные - силы Х{.

Число уравнений кинематической эквивалентности всегда равно числу удаленных связей. Группируя п уравнений (1.77), получают систему п алгебраических уравнений с п неизвестными Л/. Система ка­ нонических уравнений метода сил при внешнем воздействии в виде нагрузки = р) имеет вид:

б11 <^1 + 6 1 2 ^2 + • • • + \блХ п + Д рi = 0 ;

6 2 1 ^ 1 + 622-^2 +

+ 62/2 Хп + А2р = 0;

(1.83)

бл 1 Х\ + дп2Х2 + . . . + Ьпп Хп + ДПр —0.

 

Единичные коэффициенты

и свободные члены Д/р представляют

собой перемещения и определяются по формуле Мора, согласно которой (для плоских рам):

6*7 2 /

Mi2.dz

*ik

2 /

М М к dz

EJ

~EJ

/

 

/

 

М. М°

dz

 

 

2

1 Р

 

 

(1.84)

/

 

 

/ EJ

Для рам, имеющих затяжку —ненагруженный стержень с шарнир­ ными концами (см. рис. 1.38, б), собственное перемещение

M'dz

*зат

5,7

2 /

+

(1.85)

 

/

EJ

EF

 

 

 

где /зат - длина стержня-затяжки.

Если система имеет упругоподатливые связи, то они так же, как и затяжки, разрезаются, и б,*/ определяется по формуле, аналогичной (1.82):

M*dz

 

5// = 6i7 „ + Ф = 2 / - J — + ф,

( 1 .86)

!EJ

где б//и - перемещение, зависящее только от изгибных деформаций системы;

ф- коэффициент податливости связи.

Вматричной форме система канонических уравнений (1.83) запи­ сывается так:

D X + Dp = 0 ,

( 1.87 )

—►

—►

где X - матрица-столбец основных неизвестных; Dp - матрица-столбец свободных членов; D - матрица единичных коэффициентов:

 

6 ц

6 12

8 . л

 

D =

62 1

622

62п

(1.88)

 

5 л 1

6л2

блл

 

Элементы б//, расположенные на главной диагонали, называются главными (собственными) единичными коэффициентами; они согласно выражению (1.84) всегда положительны. Остальные элементы б .к(к Ф/) называются побочными, они могут иметь любой знак и быть равными нулю.

Основную систему следует выбирать так, чтобы как можно больше побочных коэффициентов обратилось в нуль. Этого всегда можно до­ биться при расчете симметричных систем (см. ниже).

Определив по формулам (1.84) коэффициенты 5jk и члены А/р,из решения системы (1.83) находят силы Я",-. После этого любую внутрен­ нюю силу S. можно найти в соответствии с принципом суперпозиции [см. формулу (1,81)]:

S. = 5 ? + 5. X, +

/ 2 Л 2

+

+ sjkx k +

( 1.8 9 )

)

JP

11 1

 

 

щ е SJp -

внутренняя сила в /-м сечении основной системы от нагрузки;

Sjk - то

же, но от соответствующей силы

= 1.

 

 

В соответствии с выражением (1.89) строят эпюры внутренних сил в заданной системе.

Обратим внимание на то, что на каждом этапе расчета можно вы­ полнять контроль вычислений. Например, правильность определения коэффициентов д-к и свободны* членов Д/р можно проверить путем умножения так называемой суммарно-единичной эпюры моментов М5 = 2 Л/,- саму на себя и, отдельно, на грузовую эпюру моментов Л/°. В результате должны выполняться равенства:

М г dz

 

 

M s M ° d z

 

 

~ i r =

S 5 «*

и

* { ~ j —

= БА<р-

(1>90)

Для рам с затяжками

(или с податливыми опорами) в левой части

первого равенства

надо

учесть деформацию

затяжки, как

это сделано

в формуле (1.85).

Подчеркнем, что контроль по выражениям (1.90) является лишь этапным и не отражает правильности построения единичных и грузовой эпюры.

Универсальный контроль расчета системы методом сил заключается в кинематической проверке. Если все вычисления правильны, то пере­ мещение любой точки к по направлению соответствующей удаленной связи должно быть равно нулю, т.е. должно выполняться условие (для плоских рам):

 

М .

Л /

dz

 

 

К

р

= о,

(1.91)

 

А*Р - 3 /

EJ

 

 

 

 

щ е

момент в любой

(желательно новой) основной системе от силы

= 1,

приложенной по направлению удаленной связи; Мр —момент в заданной системе от нагрузки.

Для рамы с затяжками формула универсального контроля имеет

вид

М. Л / dz

X

 

I

 

 

зат

 

 

к р

 

зат

 

Д*р = * /

EJ

+ 2

 

= о,

(1.92)

/

EF

 

 

 

 

 

 

 

где Я'зат, F3aT - усилия в затяжках и площади их поперечных сечений.

Порядок и примеры расчета. Порядок расчета поясним на примере расчета плоской рамы, изображенной на рис. 1.39, а.

1. Определяем степень статической неопределимости системы. Здесь две лишние связи. Следовательно, п = 2.

2. Выбираем основную систему с основными неизвестными (рис. 1.39,6).

3.Строим грузовую и единичную эпюру моментов, соответствующие грузовому и единичным состояниям системы (рис. 1.39, в ) .

4.Вычисляем единичные коэффициенты и свободные члены системы

канонических

уравнений

по формуле Мора (1.51): EJbn = 8/3;

Е Л 12 = —/3/3 ; Е Л 22 =

/3/ 18; EJAlp= 5,5Pl3j EJA2p = -0,25/У3

5. Строим суммарно-единичную эпюру Ms =Mi + М2 (рис. 1.39, ё)

и по формулам

(1.90) проверяем правильность вычисленных коэффи­

циентов и свободных членов.