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

книги из ГПНТБ / Биргер И.А. Резьбовые соединения

.pdf
Скачиваний:
203
Добавлен:
24.10.2023
Размер:
10.55 Mб
Скачать

Уравнение (358) и условие (359) могут быть записаны в виде неоднородного интегрального уравнения

 

q = K(q)

+ q0,

где

 

 

К ( < 7 ) = ; - Ъ jj ß ( 2 i ) jj < 7 ( z 2 ) d 2 2 d 2 l -

"Y (г)

о

о

z,

Ü5ß(2i) Ü9(г2)dZädzidz

ïf e

О О

 

H

 

0

<?o = -

_Q

я

 

v ( z ) $ y ! )

(360)

(361)

(362)

Уравнение (360) решается методом последовательных приближений по схеме

</„Чі> =

^(<7<г>) + <7о

( ( = 0 , 1 , 2 , 3 , . . . ) ,

 

(363)

где<7<(-|И <7(,-+1, — исходное

и последующее

приближения. Нулевое

приближение

можно принимать в виде

 

 

 

 

 

 

 

 

 

 

 

 

9(0)

'H^

= const.

 

 

 

(364)

В первом приближении

(при заданном усилии Q) материал считается упругим

и распределение усилий может быть

найдено также

по уравнению (235). Далее

 

по известным напряжениям находим секущие мо­

 

дули. Например, если в сечении z в болте имеется

 

напряжение о( 1 ,(г) (рис.

124),

то

секущий

мо­

 

дуль

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£* (г) =

^ .

 

 

 

 

Значения \\* и Ы в формуле (357) можно прини­

 

мать такими же, как в упругой области [см. фор­

 

мулы (216) и (217)]. Значения £*(z) и £*(z),£*e (z) и

 

Effi(z)

следует выбирать в соответствии

с напряжен­

 

ным

состоянием

в теле

болта

и гайки, а также

 

в витках

резьбы.

 

 

 

 

 

 

 

Для приближенного решения можно принять в

 

качестве такого характерного напряжения средне­

 

квадратичное

напряжение

в основании витка

 

 

 

 

 

 

 

с (z)

=«1,9,9 (г).

 

 

Рис. 124. Кривая дефор­

Более точное решение требует

использования

за­

мирования

дачи

об

изгибе

клина

в

упруго-пластической

об­

ласти.

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

НО

(см. рис. 125, где цифрами /, 2, 3 обозначены номера приближений). Сходимость улучшается, если в качестве последующего приближения (і > : 3) принять полу­

сумму двух предыдущих приближений:

1

<7<г+1> = у (<7<і> + <7<і-і>)-

Рис. 125.

Схема,

иллюстри­

Рис. 126. Схема распределения нагрузки по

рующая

процесс

последова­

виткам с учетом пластических деформаций

тельных

приближений при

 

расчете в упруго-пластиче­ ской области

На рис. 126 дано распределение нагрузки по виткам резьбы М24Х 1,5 с учетом пластических деформаций в резьбе.

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

Рис. 127. Распределение нагрузки по виткам резьбы с учетом пластических деформаций при различной высоте гаек из сплава АЛ5 (резьба М24 X 1,5; Q = 27 ООО кгс)

пластического участка практически не зависит от

длины свинчивания и

яв­

ляется

ограниченной. Этот вывод согласуется с

результатами

эксперимен­

тов, в

которых выявлена ограниченность несущей

способности

резьбовых

со­

единений.

 

 

 

Ш

Рис. 128. Модели резьбовых соединений для исследований оптическим методом:

/

— гайка с опорой по н а р у ж н о м у

контуру; 2

гайка со сферической шайбой; 3

гайка

с

конической

резьбой; 4

гайка

с

коническим

участком;

5,

6 и 7 — обычные

гайки;

8

— гайка с

д в у х з а х о д н о й

резьбой

(абсолютные

размеры

в

д ю й м а х , резьба

Витворта

 

 

 

восемь ниток

на

1")

 

 

 

 

 

 

V

 

 

2

 

 

 

S

 

 

Рис.

129.

Распределение

нагрузки

по

виткам

резьбы:

 

/ — обычная гайка;

2

гайка с опорой

по

н а р у ж н о м у

контуру;

3 — гайка

со сфери­

ческой шайбой; 4 — гайка

с д в у х з а х о д н о й

резьбой;

5

— гайка

с

конической

резьбой;

 

 

6

— гайка с коническим

участком

 

 

 

112

11.СРАВНЕНИЕ ТЕОРЕТИЧЕСКИХ ДАННЫХ

СЭКСПЕРИМЕНТАЛЬНЫМИ ИССЛЕДОВАНИЯМИ

Из формулы (235) следует, что в резьбовом соединении типа болт—гайка наи­ более нагруженным является нижний виток. Этот результат подтверждается дан­ ными практики и испытаний болтов.

Качественное подтверждение возрастания нагрузки на нижние витки было получено ранее на резиновых моделях. Тензометрирование соединения гайки с по­ лым болтом также выявило увеличение нагрузки на нижних витках.

Более обстоятельные данные были получены оптическим методом. Опыты про­ водились на плоских и пространственных моделях. Пространственные модели ис­ следовались с замораживанием образцов [4]. Размеры образцов показаны на рис. 128, а результаты исследования приведены на рис. 129. Из сравнения различ­ ных гаек видно, что существенные преимущества дают лишь коническая и сжаторастянутая гайки.

 

 

 

 

 

 

 

////?

 

Y

 

І

#

 

 

 

 

 

 

 

Ш

 

 

 

 

 

 

 

 

 

 

 

 

 

il

 

'

_

Г

 

 

 

Рис.

130.

Сравнение

теоре­

Рис. 131.

Сравнение теорети­

 

 

тического

(сплошная

линия)

ческого

(сплошная

линия)

и

 

 

и экспериментального (точки

экспериментального

 

(точки)

 

 

на

кривых)

распределения

распределения

нагрузки

по

 

 

нагрузки

по

виткам

обычной

виткам

гайки

с коническим

 

 

 

 

гайки

 

 

 

участком

 

 

 

 

Кривые показывают значения максимальных напряжений (во впадинах на­

 

резки), которые пропорциональны нагрузке на виток. Буквой k обозначено отно­

 

шение этих напряжений к напряжению в поперечном сечении болта. Две кривые

 

на рисунках соответствуют впадинам левой и правой сторон.

 

 

 

 

На рис. 130 и 131 показаны экспериментальные точки для обычной и сжато-

 

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

 

построении кривых опытные и теоретические значения q (H) для первого работаю­

 

щего витка

принимались равными.

 

 

 

 

 

 

 

 

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

Г)

первый

ее участок

считался

цилиндрическим с диаметром, равным среднему диа­

 

метру

конуса.

 

 

 

 

 

 

 

 

 

 

 

Общий вид графиков показывает, что теоретическое решение правильно охва­

 

тывает

особенности

задачи. Этот вывод можно сделать и на основании

ряда дру-

I

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

 

 

 

 

 

 

 

 

Значения коэффициентов концентрации к, приведенные на рис. 129,

являются

 

приближенными,

что объясняется малыми

размерами

моделей.

 

 

 

Г л а в а VI.

НАПРЯЖЕННОЕ СОСТОЯНИЕ И КОНЦЕНТРАЦИЯ НАПРЯЖЕНИЙ В РЕЗЬБОВЫХ СОЕДИНЕНИЯХ

1. НАПРЯЖЕННОЕ СОСТОЯНИЕ В РЕЗЬБЕ

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

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

нений.

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

сложностью геометрической формы тел болта и гайки, а также граничных условий.

В

работе

Р. Б. Хейвуда [31]

для

определения максимальных напряжений

в резьбе

предложена формула

 

°мп

 

 

 

 

 

о =

а„ +

,

(365)

 

 

 

 

 

1 . 1 / .

°о

 

где а 0

=

k'aH

— максимальное напряжение

в

середине впадины резьбы

только

от осевой нагрузки на болт; ан — номинальное напряжение в стержне болта по резьбовой части; k' — теоретический коэффициент концентрации напряжений в растягиваемом стержне с кольцевыми выточками, приближенные значения кото­ рого можно получить по формуле Нейбера; ам — максимальное напряжение только от местной нагрузки на виток—клин.

114

с = ІУбО--?-\/44

— коэффициент,

g — половина угла

профиля резьбы в гра­

дусах.

 

ом

 

 

 

Для определения

напряжения

дана

эмпирическая

зависимость

 

Р

/1,5а ,

0,45 , sin ß

(366)

 

 

 

 

 

Здесь Р — сосредоточенная нагрузка, приложенная к витку (клину) (рис. 133);

 

 

/ е

\0J

— коэффициент,

 

'

 

и еобозначены

на рис. 133;/- —

£ " = 1 + 0 , 2 6 —

)

 

размеры а,Ь

радиус закругления во впадине резьбы; S — шаг резьбы; ß — угол, определяющий

направление нагрузки относительно

касательной

к

по­

 

 

 

верхности клина

в рассматриваемой

точке.

 

 

 

 

 

 

90°\

 

Формула (366) основана на результатах

исследований

 

 

 

 

 

 

на фотоупругих моделях зубьев шестерен.

 

 

 

 

 

 

 

 

Более

точные

результаты могут

быть

получены

при

 

 

 

применении теории

упругости.

 

 

 

 

 

 

 

 

 

 

Рассмотрим плоский аналог задачи для резьбового

 

 

 

соединения типа болт — гайка .(рис. 134). Такой подход

 

 

 

вполне оправдан потому, что резьбу можно

рассматривать

 

 

 

как мелкие выточки, а кольцевыми

напряжениями

можно

 

 

 

пренебречь в сравнении с осевыми и радиальными

напря­

 

 

 

жениями.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

решения

задачи

применим

вариационный

прин­

 

 

 

цип — вариация полной потенциальной энергии тела рав­

 

 

 

на

нулю:

 

6П = Ь(Пдс)

= 0,

 

 

(367)

 

 

 

 

 

 

 

 

 

Рис. 133. Схема на-

где

Пд — потенциальная

энергия

деформации

(энергия

гружения

клина

внутренних сил); Пс

— потенциальная

энергия

внешних

 

 

 

сил.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет по уравнению

(367) может быть реализован

на ЭЦВМ

путем

решения

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

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

Для составления разностных уравнений используем уравнение для потенциаль­

ной

энергии деформации

тела

 

 

 

 

 

 

 

 

 

 

fdvY

, (ди>\Ц

,

[дѵ dw\

 

 

 

 

 

ду

 

дг J

 

г[дудг)

+

 

 

, (дода>\,

1 Пд+ѵу./дшу

 

»

 

(368)

где

G — модуль сдвига;

а,

 

 

2(і

— безразмерные

коэффи-

 

 

 

 

циенты; ц — коэффициент

 

Пуассона;

ѵ н w — компоненты

смещений соответст­

венно по осям у и г.

 

 

 

 

 

 

 

 

 

 

Потенциальная энергия

внешних

сил

может

быть определена по

формуле

 

 

 

 

 

+

Zw)

dL,

 

(369)

где Y я Z — проекции на оси ( / и г

компонентов, напряжений на границе области

с контуром L . Подставляя

соотношения (368) и (369) в равенство (367),

получим

115

вариацию полной потенциальной энергии системы

(370)

Действительными перемещениями ѵ к w будут те перемещения, которые мини­ мизируют полную потенциальную энергию при любых виртуальных перемещениях 6Ч> и ow.

6,

а)

Рис. 134. Сеточная разметка (а) и напряженное состояние в резьбе (б)

Из уравнения (370) можно получить уравнения равновесия и граничные усло­ вия для плоской задачи. Следовательно, вариационная задача (370) может заме­ няться краевой задачей.

Для получения приближенного выражения уравнения (370) контур области аппроксимируется конечным числом прямолинейных отрезков. Область покрыва­ ется согласованной нерегулярной (регулярной) сеткой, состоящей из линий у = = const и г = const (см. рис. 134, а и рис. 135) так, что точки пересечения прямо­ угольной сетки совпадают с точками пересечения отрезков, аппроксимирующих контур. Кроме того, на область наносится дополнительная сетка, линии которой

116

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

Вариационную задачу (370), записанную для непрерывных переменных V (у, г) и w (у, г) заменим вариационной задачей для дискретных значений vt

"—г—*

 

 

1

 

1

 

 

 

 

 

 

 

 

 

-

t

~

 

1

,

 

 

 

 

 

 

 

 

 

 

1

 

1

 

6

i

\2

 

5"

 

'

 

!

< 1

1

t

 

1

1

4 T .

 

 

 

 

1

 

1

r

 

f

4

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

1

 

 

1

 

 

 

- t -

 

•—(—<

? v\

 

 

 

 

1

'

1

' '

 

1

 

/

 

 

 

1

I

 

1'

 

 

 

 

1

 

1 \

1

 

 

 

 

 

 

 

 

 

 

 

- — h -

 

 

 

 

 

 

 

-

I- 4 ——

 

 

 

 

1

1

1

1

1 1

1

1

 

 

 

 

 

-

t

-

1

1

« h

 

1

1

 

 

1 1

 

 

 

 

 

 

'

 

t

'

i

1 t

1

 

 

 

1

 

1

,

1

.1

„ | _

i,

1

 

1

 

1 • 1

6)

Рис. 135. Сеточная область

и Wi в 1-х узловых точках. Последнее достигается аппроксимацией полной потен­ циальной энергии для каждой ячейки параболической функцией от ч,- и а»,-, взятых в трех узловых точках, связанных с каждой ячейкой (см. рис. 135).

Рассмотрим типичную прямоугольную внутреннюю ячейку (заштрихована на рис. 135, а). Производные для смещений в такой ячейке заменим разностными выражениями «вперед»:

дѵ _ V! —- ѵ 0

_

dv _ ti2 — fо

 

ду ~~

V

'

дг

Лог . '

(371)

dw w1 — ai0

dw w2 — w0

 

dz -

ft01

 

' 'd'y ~~

Лог

 

где hoi и Л0 2 — шаг сетки между узлами 0 и 1, 0 и

2.

 

Предполагая, что функции и их производные постоянны в ячейке и, подстав­ ляя соотношения (371) в уравнение (370), запишем выражение для потенциальной энергии деформации элементарной ячейки 201 (см. рис. 135)

" * - - e h ( a ? ) , + - ( a i ? } , + - ( ! i c s a e ) +

+ І

^

^

)

+ т ( Ѵ ) , +

т ( й

Е а ) '

Так как аппроксимации

для каждой ячейки имеют вид

(372), то,

суммируя

по всем ячейкам (см. рис. 135), получим уравнение для полной потенциальной де­

формации тела

 

 

 

ПЪ-ЪП*

 

(373)

 

 

 

 

 

117

С учетом равенства

(373) можно

записать

 

 

 

д

П

 

р

 

 

 

 

д

^ у

dnàjik

 

ш

 

дѵі

~

дѵі

'

 

 

 

 

 

і =

і

 

 

 

дП

і

=

у

 

 

(374)

 

_

ÔWi

'

 

дШі

mm.

 

 

 

 

 

1=1

 

 

 

где i— номер узла, t =

1, 2, 3,

4,

p.

 

 

 

Записав соотношения типа (372) для каждой из 12 ячеек, прилежащих к цент­ ральному (рис. 135, а) узлу 0, продифференцировав их отдельно по ѵ0 и w„ и про­ суммировав полученные выражения согласно равенствам (374), получим

 

Пл _ 9

 

 

 

W

0

W;

 

 

 

 

 

 

д

 

4«і (ftoi + Лоз)

(• =

 

І

+

2 ( Л О 2

+ Л О 4 )

2

 

w - ^ +

dw0

~~ 4

 

 

2,4

 

 

 

 

 

 

i=i,3

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

( а 2

+ 1 ) ( і \ , - и 6

+

 

я 8 - і / 7 ) ] ;

 

 

(375)

 

 

G Г

4ai(A„, + A M )

2

\

 

p

+

2(Aoi + Ao»)

J

 

+

 

 

4

 

 

( = 1 , 3

 

 

 

 

 

 

( = 2,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

(aa + 1 ) (oi„ — w b

+

w a

w7)

 

 

 

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

дѵ VI VQ _

дѵ

 

и т. д.

(376)

 

dz''

h,,

 

 

 

Однако процесс алгоритмизации существенно облегчается путем введения фик­ тивных и законтурных узлов. Типичный законтурный узел 5 показан на рис. 135,6.

В таком случае соотношения

(375) для внутреннего узла 0 сохраняют

свой вид,

а фактическая площадь ячейки учитывается введением коэффициентов q> (<р = 1

для прямоугольной ячейки,

(р^=-^и 0 соответственно для треугольной

и закон­

турной ячеек).

С учетом указанных выше положений выражение для потенциальной энергии

внешних

сил можно

представить

в

виде

 

 

 

 

 

Пс=

S

(YpM +

Zfttjblj),

(377)

 

 

 

 

/ =

і

 

 

 

 

где Yj и Zj — компоненты

внешних

 

сил, относящиеся к внешнему

граничному

узлу (/ =

1, 2, 3,

/); \

l j

— длина элемента

границы.

 

Производные от П с по

ѵ 0 и

w 0

 

 

 

 

 

 

 

д П

с

Yn\L

 

 

 

 

 

 

 

дѵ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dwr,

 

Zn&L.

 

(378)

 

 

 

 

 

 

 

 

 

118

Принимая во внимание соотношения (370), (372) и (377), из условия стационар­ ности получим разностные уравнения, связывающие смещения в некотором узле і с неизвестными смещениями в близлежащих узлах:

| ^ = 0

( і = 1 , 2 , 3

р);

'

 

(379)

где р и q — соответственно количество узлов со смещениями ѵ и w, отличными от нуля.

Уравнения для типичного внутреннего узла 0 (см. рис. 135, а), полученные из условий (379), имеют вид

 

 

1 = 2 ,4

 

1

=

1,3

'

+

( а 2

+ 1) {ve — Vs +

Vg —

v,) ••

Z0\l0;

 

4а,(Л„2 +

/г 4 )

2 -jf +

2 ( А о і + Л и )

 

2VjTÏr +

 

 

1=1, 3

 

1 =

2,4

 

+

( а 2

+ 1 ) ( ш в — a ) 5 + t t ) s

a;7 )j =

F 0 A/ 0 .

(380)

Нумеруя последовательно узлы и записывая уравнения (380), получим р + q линейных разностных уравнений, которые можно записать в матричной форме

 

 

Ах^В,

(381)

где А — симметричная

положительно определенная матрица

коэффициентов,

ранг которой равен +

<7+

<?); ^ — вектор-столбец с неизвестными компо­

нентами и,- и Wi, В — известный

вектор-столбец, характеризующий внешние

нагрузки.

 

 

 

Следует отметить, что А — сингулярная матрица, и решение не будет единст­ венным, если не известны три компонента перемещений: у; в двух точках и ОУ;В од­

ной точке

или наоборот.

(381), получим значения для перемещений Vi и Ш;

Решая

систему уравнений

в узлах области.

 

Напряжения вычисляются

через перемещения по формулам:

ѵг0 , w1—w0\

r n o г _ 2Gp,

 

— 2ц"

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

(см.

разметку области на рис. 134, а) может превысить тысячу.

 

Решение системы уравнений с таким числом неизвестных затруднительно

на

ЭЦВМ.

Настоящая задача решалась * на ЭЦВМ БЭСМ-4 с использованием в качестве подпрограммы программы для решения систем линейных алгебраических урав­ нений по методу квадратного корня. Подпрограмма ограничивала количество внут­ ренних узлов сеточной области і < 130. Поэтому при расчете напряженного сооо-

Расчеты выполнены И. В . Рокитянской .

119