
книги из ГПНТБ / Биргер И.А. Резьбовые соединения
.pdfУравнение (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 |
|
ш |
||
|
дѵі |
~ |
Zà дѵі |
' |
|
||
|
|
|
|
і = |
і |
|
|
|
дП |
і |
= |
у |
|
|
(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