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

1167

.pdf
Скачиваний:
62
Добавлен:
22.03.2015
Размер:
1.44 Mб
Скачать

 

1

 

 

 

∂ϕ

 

 

1

 

2ϕ

 

 

σr =

 

 

 

 

 

+

 

 

 

 

;

 

r

 

r

r2

∂θ2

 

 

 

 

 

 

 

 

 

σθ =

2

ϕ;

 

 

 

 

 

 

(159)

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

∂ϕ

 

 

τrθ = −

 

 

 

 

 

 

∂θ

.

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

 

Подставляя эти соотношения в уравнения равновесия (143), убеждаемся, что при отсутствии объемных сил уравнения равновесия обращаются в тождества. Чтобы преобразовать уравнение сплошности (144), сложим почленно формулы для нормальных напряжений (159)

σr θ =

1

 

∂ϕ

+

1

 

2ϕ

+

2ϕ

.

r

r

r2

∂θ2

r2

 

 

 

 

 

 

В правой части этой суммы согласно формуле (145) стоит оператор Лапласа над функцией ϕ(r,θ) . Следовательно,

σr θ = 2ϕ.

Подставляя это соотношение в уравнение сплошности (144), получаем:

2 2ϕ = 0

или

4ϕ = 0.

(160)

В развернутом виде это уравнение представляется следующим образом:

 

2

1

 

 

1 2

 

2ϕ 1

 

ϕ

 

1 2ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0 .

 

 

2 +

 

 

 

+

 

 

2

 

2

 

2

+

 

 

 

+

 

2

 

(161)

 

r

r

r

 

r

∂θ

 

r

r

r

r

∂θ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

функция напряжений ϕ(r,θ)

для плоской за-

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

6.8. Осесимметричные задачи. Решение в перемещениях

Остановимся на плоских задачах, в которых напряжения не зависят от полярного угла θ. Как следует из зависимостей (159) в

141

этих задачах функцию напряжений можно считать не зависящей от θ. Бигармоническое уравнение (161) примет более простой вид

 

d 2

1 d

d 2

ϕ

1 dϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0

 

2 +

 

 

 

 

2 +

 

 

 

 

r

 

 

dr

r

 

dr

 

 

 

dr

 

 

 

dr

 

или после дифференцирования

d 4ϕ

+

2

 

d 3ϕ

1

 

d 2ϕ

+

1

 

dϕ

= 0.

(162)

dr4

r

dr3

r2

dr2

r3

dr

 

 

 

 

 

 

 

 

Также упростятся выражения напряжений через функцию

ϕ(r) из (159)

 

 

 

 

 

 

 

 

 

 

 

σr =

1

 

 

dϕ

;

 

 

σθ =

 

d 2ϕ

; τrθ = 0.

(163)

 

 

 

 

 

 

 

 

 

 

 

r

 

dr

 

 

 

dr2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнений равновесия (143) при отсутствии объемных сил

останется только одно

 

 

 

 

 

 

dσr

 

 

 

 

σr

−σθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

= 0.

 

 

(164)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формулы Коши (147) упростятся, так как составляющая пе-

ремещения uθ =0 в силу симметрии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εθ =

ur

;

 

εr =

dur

;

γrθ = 0.

 

 

 

(165)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из формул (148) закона Гука останутся лишь две

 

 

 

 

 

εθ =

 

1

(σθ −µσr ); εr =

1

 

 

(σr −µσθ).

 

 

(166)

 

 

 

 

E

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Осесимметричную задачу в перемещениях можно решить в

общем виде. Для этого из формул (166) находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σr =

 

 

 

 

 

 

 

 

 

 

(εr +µεθ);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1−µ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σθ =

1−µ

2

(εθ +µεr ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью формул Коши (165) исключаем из уравнений (а)

составляющие относительных деформаций εθ и εr

 

 

 

 

 

 

 

E

 

du

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

u

 

 

 

 

du

 

 

σ

r

=

 

 

 

 

 

 

 

r

 

 

r

;

 

σ

θ

=

 

 

 

 

 

 

 

 

 

 

 

r

 

r .

(б)

1

−µ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

r

 

 

 

 

 

 

 

1−µ2

r

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

142

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя (б) в уравнение равновесия (164), получаем следующее дифференциальное уравнение относительно составляющей перемещения ur:

d 2u

r

+

1

 

du

r

u

r

= 0.

(167)

 

 

r

 

 

dr2

 

dr

r2

 

 

Это дифференциальное уравнение с переменными коэффициентами. Для решения приведем его к уравнению с постоянными коэффициентами с помощью подстановки

t = ln r

 

 

(в)

или

 

 

 

r = et .

 

 

(г)

Дифференцируя соотношение (в) по переменной r, получаем:

 

dt

=

1 .

(д)

 

dr

 

 

r

 

Найдем связь между производными от функции ur по старой и новой переменным

dudrr = dudtr drdt .

С учетом соотношения (д) получаем

dur

=

1

 

dur

.

(е)

dr

r

 

 

 

dt

 

Вторая производная

d 2u

r

 

 

 

 

d

du

 

 

d

 

 

 

 

=

 

 

 

 

 

 

 

 

 

r

=

 

 

dr2

 

dr

 

 

 

dr

 

 

 

 

 

 

dr

 

 

1

 

 

 

 

d

 

du

 

 

 

 

1

 

+

 

 

 

 

 

 

 

 

 

 

r = −

 

 

r

dr

 

 

r2

 

 

 

 

 

 

dt

 

 

 

 

1

 

 

 

2

ur

 

 

 

 

dur

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

r

2

 

 

 

 

2

 

 

dt

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

1

 

du

r

 

= −

 

1

 

 

du

r

+

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

r

 

 

dt

 

 

 

 

 

dt

 

 

du

r

+

1

 

d 2u

r

 

 

dt

=

(ж)

 

 

 

r

dt2

 

 

dr

 

 

dt

 

 

 

 

 

 

 

 

Подставляя производные (е) и (ж) в уравнение (167), нахо-

дим:

d 2ur ur = 0. dt2

143

Решение этого уравнения имеет вид ur = Aet + Bet .

Переходя к старой переменной r согласно зависимостям (в) и

(г), получаем следующее решение уравнения (167):

 

 

 

 

ur = Ar +

B

.

(168)

 

 

 

 

 

 

 

 

 

r

 

Так как

dur

= A

B

, то после подстановки в уравнение (165)

dr

 

 

 

r2

 

значения ur и ее производной по r, получим

εr = A rB2 ;

εθ = A + B2 , r

а из формул (б) – составляющие напряжений

 

 

 

E

 

 

 

 

)A (1−µ)

B

 

 

 

σr =

 

 

 

 

 

(1

 

 

 

;

1

−µ

2

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

B

 

 

 

 

 

 

 

)A +(1−µ)

 

σθ =

 

 

 

 

 

 

(1

 

 

 

.

 

 

 

 

2

 

 

2

 

 

1−µ

 

 

 

 

 

r

 

 

 

 

Постоянные A и B определяются из граничных условий.

(169)

(170)

Контрольные вопросы

1.Что называется простым радиальным напряженным состоянием?

2.Как определяются постоянные интегрирования после решения уравнений равновесия и сплошности для клина, нагруженного продольной силой?

3.В чемсмысл задачи Буссинеска?

4.Чтотакое круг Буссинеска и каковы его свойства?

5.Как записываются функции напряжений для плоской задачи теории упругости в полярныхкоординатах?

144

ЗАКЛЮЧЕНИЕ

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

Принимая во внимание, что студенты только что закончили изучение курса «Механика материалов и конструкций» («Сопротивление материалов»), в котором особенно подробно рассматривалась проблема внутренних сил, в приложении дается краткое объяснение метода сил в применении к замкнутым контурам. Это требуется для решения методом конечных разностей балки-стенки.

Выражаю свою искреннюю признательность официальным рецензентам, а также коллегам по работе, оказавшим неоценимую помощь при подготовке рукописи.

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

600000, Владимир, ул Горького, 87, кафедра «Тепловые двигатели и энергетические установки».

E-mail: hotz@mail.ru; hotz@vlsu.ru

145

ПРИЛОЖЕНИЕ

1. РАСЧЕТ ЗАМКНУТЫХ КОНТУРОВ16

1.1. Замкнутые контуры с одноосной симметрией

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

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

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

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

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

16Расчет замкнутых контуров методом сил рассматривается в курсе «Механика материалов и конструкций».

17Предполагается, что студенты уже прослушали курс «Механика материалов и конструкций» и знакомы с расчетом статически неопределимых систем методом сил и умеют составлять канонические уравнения (см., например, [6]).

146

необходимо
Рис. 1. Схема расчета замкнутого контура с одноосной симметрией

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

щим моментом X2.

Система канонических уравнений для этих двух неизвестных будет иметь вид [6]

δ11 X1

δ21 X1

+δ12 X 2

+δ22 X 2

+1,P

+2,P

= 0,

(а)

 

= 0.

 

 

 

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

(в направлении моментов X2). Другими словами, при действии неизвестных пока сил X1 и моментов X2 основная система (с разрезом по нижнему стержню, см. рис. 1, б) будет эквивалентна заданной (см. рис. 1, а), т. е. сечения в разрезе не будут иметь ни линейных, ни угловых перемещений.

На рис. 1, б, в и г построены эпюры моментов M p , M1 и M2. Напомним, что после построения эпюр M p , M1 и M2

определить перемещения от действия каждой неизвестной силы или, другими словами, для определения δ11 в (а) перемножить эпюру M1 саму на себя, для δ12 = δ21 M1 – на M2, а далее – для определения 1, p и 2, p необходимо M p перемножить соответственно

на M1 и M2, используя, например, правило Верещагина. Оно гла-

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

147

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

Вычисляя перемещения по этому правилу, определяем коэффициенты системы канонических уравнений (а):

EIδ =

1

 

3a 3a

2

3a 2 +3a 2a 3a = 36a2 ;

 

 

 

 

 

11

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EIδ22 = 2 3a 1+ 2 2a 1 =10a;

 

 

 

 

EIδ =

EIδ

 

= 2

1

3a 3a 1+3a 2a 1 =15a2 ;

21

 

 

 

12

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

26

 

 

EI1,P = −

 

 

1

 

Pa a

 

8

a 2 Pa 2a 3a = −

Pa3

;

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

EI2,P = −

 

1

Pa a 1 2 Pa 2a 1 = −3Pa2 .

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система канонических уравнений (а) примет вид

 

 

 

36a2 X1 +15a2 X 2 =

26

 

Pa3 ,

 

 

(б)

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

15a

2

X1 +10aX 2 = 3Pa

2

 

.

 

 

 

 

 

 

 

 

 

 

 

Решение системы (б) дает величины

 

 

 

 

 

X 1=

25

P = 0,309P и

 

X 2= −

22

 

Pa = −0,163Pa.

 

 

 

 

 

 

81

 

 

 

 

 

 

 

 

 

 

 

135

 

 

 

 

 

На рис. 1, д и е построены эпюры M1 X1, M 2 X 2 , а на рис. 1, ж

суммарная эпюра для заданной статически неопределимой системы, которая строится суммированием ординат от заданной нагрузки, а также найденных нагрузок X1 и X 2 : M = M P + M1 X1 + M 2 X 2 . Проверка правильности построения эпюры M проводится путем

перемножения этой эпюры с любой единичной. Если суммарная эпюра сложной формы, то перемножать можно эпюры M P , M1 X1, M 2 X 2 , например, на эпюру M2 (ординаты ее постоянны). В

нашем случае угловое перемещение в сечении разреза (см. рис. 1) равно:

EIθ = −2 12 Pa a 1 Pa 2a 1 + 2 12 0,926 Pa 2a 1 + 0,926 Pa 2a 1 − − 2 0,163 Pa 2a 1 2 0,163 Pa 3a 1 = −4,63 + 4,63 0.

148

Таким образом, решение проведено правильно.

Наибольший изгибающий момент M max = 0,454Pa получается в

сечениях, совпадающих с точками приложения сосредоточенных сил P. Напомним, что по суммарной эпюре моментов можно представить себе форму стержней после деформации (эпюра моментов строится на сжатой стороне стержня, а нулевая точка на эпюре моментов соответствует месту перегиба упругой линии). На рис. 1, ж форма упругой линии показана пунктиром.

1.2. Замкнутые контуры с двухосной симметрией

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

 

 

 

по конструкции, так и по составу

 

 

 

 

 

 

нагрузки, для решения целесооб-

 

 

 

разно применять сквозной разрез

 

 

 

по линии симметрии, разделяя

 

 

 

контур на две отдельные части.

 

 

 

Например, системаввидеквадратно-

 

 

 

го замкнутого контура (рис. 2, а),

 

 

 

а пусть стержни, составляющие

 

 

 

контур – постоянной жесткости.

 

 

 

Рама

загружена

самоуравнове-

 

 

 

шенной сплошной

нагрузкой,

 

 

 

равномерно

распределенной по

 

 

 

боковым стержням, и решается с

 

 

 

помощью

основной системы,

 

 

 

показанной на рис. 2, б. По усло-

 

 

 

виям

симметрии,

неизвестные

 

Рис. 2. Метод решения замкнутого

 

силовые факторы в этом случае

контура с двухосной симметрией

получатся в верхнем и нижнем

сечении одинаковыми, причем осевые силы

N = qa

2

определяют-

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

На рис. 2, в после условного разреза заданная система показана в деформированном состоянии. Согласно предыдущему, замеча-

149

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

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

в плоскости разреза не может быть поперечных сил вследствие симметрии;

осевые силы определяются непосредственно по условиям равновесия основной системы под заданной нагрузкой (в данном случае они равны нулю);

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

δ11 X1 + ∆1,P = 0.

(а)

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

На рис. 2, г, д построены эпюры моментов МР , М1 для левой половины основной системы, показанной на рис. 2, в (эпюра Мр построена по правилам расслоения – относительно среднего сечения определяются моменты отдельно от продольной силы и распределенной нагрузки).

Согласно эпюрам, по правилу Верещагина получим

EIδ11

= 4 1 a

1 = 2a;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI1,P = 2

1

 

qa2

 

a

1

2

1

 

qa2

 

a

1 =

qa3

,

2

4

2

3

8

2

12

а потому согласно (а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qa3

 

 

= − qa2 .

 

 

 

X1 = −

 

1,P

= −

 

 

 

 

 

 

 

 

12 2a

 

 

 

 

 

 

 

 

δ11

 

 

 

 

 

 

24

 

 

Вследствие симметрии вычисления проводили только на по-

ловину системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 2, е построена эпюра M1 X1

и на рис. 2, ж – эпюра M ,

ординаты которой получены на прежних основаниях по формуле

M = M P + M1 X1.

150