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

2768.Несущая способность и расчёт деталей машин на прочность

..pdf
Скачиваний:
42
Добавлен:
15.11.2022
Размер:
28.78 Mб
Скачать

Напряженное состояние при упруго-пластическом деформировании

61

V еъ + еЬ+ еч>ев»

Рис. 31. Значения функций Фи; Фр и Фир при ц = 0,5 и % = 1,2

ряемое вдоль меридиана оболочки от произвольной точки в направлении интегрирования, D — цилиндрическая жесткость, коэффициент Пуассона ц = = 0,5; интегральные функции пластич­ ности

где

осевая

деформация

е ^ = ~ —

w

dd

 

 

 

ds

 

 

 

 

5

- ±

у -г- и кольцевая деформация

A l

wS

 

 

 

 

£0=

и

.

ш .

,

±

и

— cosip---- - sm

 

у — cosoj).

Эту систему уравнений решают ме­ тодом последовательных приближений на ЭЦВМ, используя в частности метод Рунге-Кутта с последующей ортогонализацией решений [31. В последователь­ ных приближениях .определяют значе­ ния интегральных функций пластично­ сти для последовательно вычисляемых значений деформаций.

Для цилиндрических оболочек урав­ нения (1.200) существенно упрощаются. Уравнение равновесия при отсутствии осевой силы имеет вид:

 

 

 

d / _ d2w

\

 

 

 

 

d^K ^

Фи)+ ИМ)Р + /,= 0' (1 202)

ф и= - |

§ Ф01-

( 1.201)

Eh3

 

 

Дпо)2^ ;

где D — -д—(цилиндрическая жесткость

1

 

при (л =

п п

,

Eh

 

\

 

 

 

0,5),

R=

- ^ , р — нормальное

Фир = 2

} фТ)ФВ

 

давление.

 

 

 

смещение срединной поверхности обо­

 

 

 

 

лочки в

сторону

внешней нормали

 

 

 

 

 

1

 

 

 

 

 

При приближенных расчетах в силу малости отношения Фир/Фр можно принять Дг]0 = 0 и определить функ­

ции Фи, Фр по рис. 27, гл. 11. Интен­

 

 

сивность деформаций определяется по

Рис. 32.

Схема размеров, перемещений

формуле

и усилий

оболочки вращения

62 Расчет на прочность при статическом нагружении

В этом уравнении интегральные получим функции

1

или

 

 

(1.206)

 

 

 

* <

= ■ / 4 + 4

,)2-

(1-207)

1

 

 

Относительные усилия

через дефор­

мации

 

 

Ту = — ё рФ р

И Мх= -----g ^ и ш а х ® ц ,

 

 

(1.204)

где

 

 

w

_d2w

h

еР = ае~ e* max = 'dx2

2е~*

Максимальная интенсивность дефор­ маций (при Т] = 1)

e i m a x = ] / r ^р + ^ н т а х у *

<1208)

Интенсивность деформаций на нейт­ ральной поверхности

^i0= ^р-

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

г»

_______

ei шах

J ё; У Ц -

 

 

{

Si dei

1

J

J i /

j T ? F +

 

y

i j + e i

(e i max e i0)

Интегральные функции пластично­ сти удобно выразить через интенсив­ ность деформаций при изгибе цилинд­ рической оболочки:

= |7 з ~ ^ ё 'хх ё ххё уу ё уу'

Подставив в это уравнение выраже­ ния для компонентов деформаций

ехх =- - ё х -

f d 2w

h

л;

\ dx~

'* т

 

(1.205)

ё уу = ■ ё2 =

W

1

а

ет

 

 

 

^ фё,- У ё \ - ё \ й dё2i|

(1.209)

С^ i d S j _

J. V n - * h

При ё,-о> 1 первый интеграл следует

считать равным нулю, а интегрирова­ ние второго вести в пределах от ё;0Д°

е.гшах*

После интегрирования запишем вы­ ражения для функций пластичности при изгибе оболочки для полигональ­ ной аппроксимации:

при г 1о < 1

Фи =

,

.

v .

 

{еСш

■*?.)

 

+

(е;

3-

е п

1 м + 1 ’Л ) - ,

 

Vе I шах

c i0)

■“

€ zz —

(е Х Х ~\~е уу)>

(1.210)

Напряженное состояние при упруго-пластическом деформировании

63

ф р=

eI0

® « = --------5—

т 1 (1 -

 

У ^ т а х

{^1 шах

^(а)

 

 

 

 

+ VH 1

2 м + ьЛ );

V el max

W„

при el0 5= l

 

лз

и _

 

 

з х

 

 

(^i max

*1.)

 

X 2

K ^ ' +

V

;;)I

(1 .2 П)

Фр =

У е?тах -

е?0 X

 

Функции У„ и Ур зависят лишь от е,- и ё( , поэтому, вычислив их значения

для различных деформаций и экспери­ ментально определив ап и Ьп легко найти Ф„ и Фр. Значения функций У„ и У для различных ё;тау и ё(.= 1; 1,25;

1,5; 2; 3; 4; 5, а так же значения всех других величин, входящих в функции Фи и Фр, приведенные в работе [22].

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

диаграммы ап = 1 — GT и bn = GT. За начальное значение ё. принимают ё(.

(или 1), за конечное значение — et.max.

Используя данные работы (21) удоб­ нее записать выражения для функций

пластичности в следующем виде:

Фи = - О -* ? .)

 

Т +

(е ?шах

2

ei0)

3

£

 

(е ? шах

52 Ч2

е Г0)

х [ ( 1 - 0 т ) 2 ^ + 0 т 2 ^ 1 [<L212>

фр=

r^Fmax

+ T / g2

_ ^ Х

К Wmax

£0

х [(1 - о, ) 2 ^р + О т1 ; ^ ] ;

- о т) 2 ^ + 5

2 Я .

 

 

• (1.212а)

ф р -

_ 1 _ = [ ( 1 -

V

е i шах

е 1а

— GT) ^ y p + GT ^ У р ],

Графики функций Фи и Фр для мо­ дуля упрочнения Gx = 0 показаны на рис. 27, гл. 11.

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

В результате получим

av

+ \ р ы * х + тх[оф" Ш = < ,;-

при этом надо иметь в виду, что

s [ D ® " S ] = - Q’ r i e ° - nepepe-

зывающая сила в сечении оболочки. Изгибающий момент в сечении обо­

лочки

d-w

— М = ДФ„ dx1

XXt

=— ^ ^ &Фрш dx dx

av av

-f ^

p (A:) dx dx — J Q (av) dx-f

av a*v

aV

 

(1.213)

64 Расчет на прочность при статическом нагружении

отсюда

X Xi

аv аv

X X

+oWS Sp(x)dxdx~

аv аv

-c k iQMdx- W : (,-24)

X i x г

=5 * ® » » * * * +

S

flv av

XX i x%

+$ 5 ^ - $ [ p W i x d x t o -

°V

a\

a\

 

 

X

X i

 

— Q(ev» = 5экг S ‘ted*—

 

- M <‘ v

) $

^ + ( £ )

.<>.2 Ш,

Интегральное уравнение для прогиба

Система уравнений (1.213) — (1.216) дает решение задачи об упруго-пласти­ ческом деформировании оболочки. В этой системе геометрические пара­ метры оболочки характеризуются ве-

Е№ Eh

личинами D = - - и k = — . Уравне­

ние (1.213) можно интегрировать при плавно меняющейся толщине оболоч­ ки h, при этом D и К будут функ­ циями длины оболочки; входящая в уравнение нагрузка р (х) может иметь произвольный закон изменения с тем ограничением, что функция нагрузки должна принимать на интервале интег­ рирования конечное значение и иметь конечное число разрывов первого рода. При известных значениях ш и d'2w по

формулам (1.205) легко определить значения напряжений ахх и ауу в лю­ бом сечении оболочки.

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

X Xi XЯ х3

" ““ И в к

И к®*1»

 

 

 

а.,

а.,

 

а., а,,

 

 

 

v

v

 

v

v

 

 

X

X i

 

Х г

Хз

 

X

X i

+ 5

 

 

 

 

W d x d x d x i x - Q ^

 

- ^ t o d t d x -

а

а

 

а

а

а

 

а

у “у

 

у

у

 

у иу

- М^ П Ш +[яЪ-а.5 dX+'"

(1.216)

 

Интегральное уравнение для про­ гиба w можно решить методом последо­ вательных приближений, если заданы условия на границах интервала интег­ рирования. В этом уравнении av — (значение границы интервала) равно а или Ь.

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

Напряженное Состояние при упруго-пластическом деформировании

65

линейность уравнения, отличаются от 1

методов. Для инженерных расчетов не­

больше, чем в других возможных слу­

обходимо найти достаточно точный и

чаях упрочнения.

 

оказывается

простой способ вычисления коэффициен­

Сходимость

процесса

тов концешрации напряжений и дефор­

весьма медленной, и лишь четвертое

маций, позволяющий оценить уровень

приближение

может

быть признано

напряжений в зонах

концентрации.

удовлетворительным,

так

как переме­

Наибольшее распространение для та­

щения w в третьем и в четвертом при­

ких целей получила формула Нейбера

ближении различаются

примерно на

КаКе= а*,'

 

(1.217)

3%. Для улучшения сходимости про­

 

выражающая коэффициенты концентра­

цесса вместо простой

итерации можно

воспользоваться подобной. В этом слу­

ции напряжений Ка

идеформаций К,е

чае в л-м приближении можно запи­

через значение коэффициента концент­

сать

 

 

 

рации в упругой области аст.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф,

 

где

Ш(£т) =

ы) (Е ) Здесь wTt0) — про-

Формула

Нейбера

была

получена

 

 

 

w t

( 0 )

 

 

 

 

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

 

 

 

 

 

 

 

[24]. Распространение ее на все другие

гиб в нулевом сечении при дости­

жении в нем предела текучести.

 

случаи

концентрации

носит

прибли­

Кроме того, для улучшения сходимо­

женный характер.

 

 

сти

последовательных

приближений

Изучение этой формулы, проведенное

можно улучшить

процесс отысканий

Н. А. Махутовым, показало, что соот­

параметров упругости,

определив зна­

ношение (1.217) можно уточнить, если

чение

_

/d % \

Фи о,,

 

учесть зависимость коэффициентов кон­

(d2w\

 

центрации от уровня номинальных на­

\ d ^

)п +1

 

\ d ^

)п Фц (л ц> ’

 

пряжений

и

степени

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

 

 

ского упрочнения [12]

 

а по нему и по ранее определенному

 

значению

прогиба

wn получить

окон­

}S ^± = F[cc0o j { o nen)\,

(1.218)

чательное значение Фи для п-го при­

 

 

 

 

 

 

ближения.

 

з а д а ч

о к о н ц е н ­

где оп, ёп —соответственно номиналь­

Р е ш е н и е

т р а ц и и

н а п р я ж е н и й

при

ные напряжения и деформации, отне­

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

деформировании

сенные к величинам, соответствующим

связано с существенными трудностями,

пределу текучести.

 

 

поэтому

получили

распространение

На основе обработки ряда теоретиче­

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

методы,

эксперимен­

ских решений и большого экспери­

тальные методы (с помощью оптически

ментального

материала в работе [ 12]

активных покрытий, метода муара и

показано, что функция F для степен­

малобазных тензорезисторов) и методы

ного упрочнения tfmax= ^max имеетв|,Д

решения

краевых задач

(вариационно­

 

 

 

 

 

 

разностный различных модификаций и

 

-

v O .S d - m J J l - ^ - 1^ ) ] '

конечных элементов) с помощью ЭЦВМ

 

14,11, 14]. Дальнейший прогресс будет,

 

 

 

 

 

( 1.2

по-видимому, достигнут на пути сочета­

[

°

п )

 

 

а для линейного упрочнения

 

ния

расчетных

и

экспериментальных

 

 

 

F =

 

 

 

1

 

 

 

 

 

 

(1 .2 2 0 )

 

 

(

~ '»0'5

С - о т) [ | - ( 5 „ -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\а0ап)

 

 

 

 

 

 

 

 

 

3 Серенсен и др,

66

Расчет на прочность при статическом нагружении

Зависимости функции F для тепло­ устойчивой стали = 0,08) в широком диапазоне значений коэффициента кон­ центрации аа и максимальных напря­ жений приведена на рис. 33. Точками обозначены результаты экспериментов, штриховыми линиями — значения F, вычисленные по уравнению (1.219) при

оп = 0,5 и 0,55. Можно отметить, что при значениях максимальных упругих

напряжений а 0а„ в диапазоне до 5, что соответствует обычно встречающимся в практикезначениям, функция F соста­ вляет около 0,7 и выше, при значе­

ниях

оп = 1

функция

F

не превы­

шает 0.8.

линейного

упрочнения

Для

случая

между

коэффициентами

концентрации

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

К„ = =!- +

& - < К А .- 1).

(1 -221)

0/1

Оп

 

Используя уравнение (1.221), можно записать значения коэффициентов кон­ центрации в упруго-пластической об­ ласти в следующем виде:

Зависимость между максимальными деформациями на контуре отверстия и номинальными напряжениями, полу­ ченная расчетом и из эксперимента для полосы с отверстием при растяжении показана на рис. 34. Соответствие ре­ зультатов расчета по точному методу, по приближенному и эксперимента ока­ залось достаточно хорошим. На этом рисунке верхняя шкала относится к двум верхним кривым, построенным по расчету по приближенным формулам для теплостойкой стали двух различ­ ных термообработок. Точки соответст­ вуют экспериментальным данным, по­ лученным методами муара (А) и се­

ток (Л).

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

Хорошее соответствие результатов наблюдается и при расчете оболочек с неукрепленными отверстиями и пат­ рубками по приближенным формулам и точным методам с помощью ЭВМ (рис. 35) [12].

КР=

 

 

Л - с т V 1/2

_ 1 — GT

 

 

 

2GjOn

 

A ( < v J 0,‘ ( l - a T ) [ l " ( i * " ,/“e)1

K2Glial

 

 

при

оп ^

1;

 

 

 

 

a -о.

1/2

Ке =

 

i - а ,

1 —От

 

°.5 (1 - GT) [! - (о„ - 1/а0)] + V 2GTё п

2G1en

 

Gren (aoOri)

 

 

при

on

1 ;

 

( 1. 222)

 

 

а Ъ ° Т

 

Ka =

 

1 - С т у 1/2

1 - GT

 

 

/

2оп

 

 

 

при

оп ^

1;

 

 

К0 =

а аепс2 т

1 - С ТХ2 1/2

1 —от

 

+

 

2оп

при On ^ 1.

Напряженное состояние при упруго-пластическом деформировании

67

 

 

 

 

 

 

 

 

ке =1

 

 

 

\

 

 

 

f.77iп •

 

 

 

 

 

 

N

 

 

 

 

\

 

 

 

 

5^а

 

 

 

 

 

 

 

 

 

'' РЬ";--

 

 

 

 

\

 

 

 

 

 

 

 

 

7771

 

 

\

 

 

 

 

 

 

 

щ

 

 

\

\

 

 

 

 

 

 

 

 

 

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ K G=CK6

 

 

0,1 --------------------------

 

 

-----

1

I

I

I

N\\-------------

------- ----------

1-

1 М М

/

2

3

О 5

6

7

8

9 1 0

1 0 3 0

0-0 5 0

60 70 <Х6 ,0П

Рис. 33.

Функции F

для

теплоустойчивой стали

 

 

 

3*

Напряженное состояние при упруго-пластическом деформировании

69

известно это соотношение для некото­ рого модуля.

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

например М = Ф^шах* N= Фрер при изгибе пластинок

Л// /. = Фц-д- ( 2 ^ г т а х + ^в m a x);

— Фи (2^о тах ^^гтах );

Mi = Ф(?,-тах и т. д. Рассмотрим свой­

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

Функция пластичности в точке при линейном упрочнении с модулем GTl

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

Так как для Gt2 и Gt1 при деформациях ёу = ёг усилия в сечении Q2 — Ф2?х и Qi — Ф^х, то из уравнения (1.228) можно получить выражение (1.225), которое, следовательно, справедливо не /только для предельного случая, но и для любого значения деформаций.

Соотношения (1.225) и (1.226) можно записать в форме комбинации реше­ ния упругой и упруго-пластической задачи:

Q _Q

^Tl

I

и/2

— Vynp —

-=

г

 

1

— GTl

 

+

1 - G T I

П.229)

 

 

При одном и том же значении дефор­ мации функция пластичности ф2 при

модуле <?та определяется через функ­ цию фх:

1

(' —Ч>0-

(1-227)

1 —GT1

Имея в виду, что интегральные функ­ ции пластичности Ф связаны с функ­ цией ф выражениями типа

Ф = К ^ фГ)2 dr) или Ф = К ^ фг) dr),

л л

можно записать уравнение, аналогич­

ное уравнению

(1.227):

 

Фв= 1 —-J— ^

(1 —Фх);

(1.228)

 

1

GTl

 

причем Ф =

1

при ф =

1 .

Уравнение

(1.228) справедливо при

условии,

что

границы

пластической

области

не зависят от модуля линей­

W - g-

GTl

(1.230)

GT2

1 — GT

 

в предположении, что тело деформи­ руется упруго усилием и С?упр = ёупр.

Следует еще раз подчеркнуть, что соотношения (1.225) и (1.229) справед­ ливы, когда при переходе к другому модулю упрочнения сохраняются по­ стоянными все составляющие деформа­ ций, а соотношения (1.226) и (1.230) — когда сохраняются постоянными все составляющие усилий. Например, при совместном действии изгиба и растяже­ ния выражение (1.229) можно исполь­ зовать для одинаковых ё[1ГПах (или ёр)

@Р

их — ,а выражение (1.230)—для

е и max

_

-

N

одинаковых М

(или N) и А, =

— ; при

 

 

М

изгибе пластинок выражение (1.229) — для одинаковых е.тах (или ёгтак и

ёд max) и выражение (1.230)— для оди­

наковых Mi = Vr^ r ~ МгМ$-1- Мв(или

МГ И 2Ид) И Т. Д.

70 Расчет на прочность при статическом нагружении

Пользуясь уравнениями (1.229) и (1.230) можно выразить остаточные де­ формации в следующем виде:

в предположении равенства харак­ терных деформаций

(1-231)

1 — GTl

и в предположении равенства усилий

г ш в = т ^ й - ^ г М1-

(1-232)

1 ит1 иТ2

 

Для приближенного определения на­ пряжений можно использовать уравне­ ние (1.230), подставив его в выражение

о2 = (1 — Gt2) + GT2е2 (при линейном

упрочнении Gx2). Напряжение при усилии Q, соответствующее характер­

ной деформации при модуле Gt2, опре­

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

с модулем GTl и напряжение аупр при упругом деформировании:

G-г2

GTi

<*2

ступр"Ь

1 - G n

+

(1.233)

1 — GTt

 

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

Для статически определимых задач, когда нагрузки и усилия в сечении про­ порциональны, формула (1.226) для пе­ ремещений принимает вид

Д2 = Q I - 12

Gt1

7Г" +

 

CTI

Gt2

 

Gt2

(1.234)

1 - G , i

GT1

 

По этой формуле получают завышен­ ные значения нагрузок и заниженные

значения

перемещений

по сравнению

с точным решением.

 

Используя выражения (1.225) для оп­

ределения

зависимости

нагрузок от

перемещений, можно записать

Q2= Q 1

1 GJ2

А, ° « - Сч , (1.235)

 

1-G,x

1 — <5„

Рис. 37. Сопоставление результатов при­ ближенных решений и точных при изгибе балки

В этом случае, как показывают вы­ числения, по сравнению с точным ре­ шением значения нагрузок оказы­ ваются заниженными, а перемещения (при одинаковых нагрузках) — завы­ шенными, при этом погрешность может составлять до 10%. На рис. 37 приве­ дена зависимость, полученная точным решением для изгиба балки на двух опорах силой, приложенной посредине

при Gt2 = 0,2 (сплошная линия). Эта зависимость сопоставлена с другими, полученными по приближенным фор­ мулам (1.234) и (1.235) (для того же зна­

чения Gt2). Приближенные решения ограничивают точное сверху и снизу.

Для статически неопределимых задач упруго-пластического деформирования внешние нагрузки и усилия в сечении не пропорциональны, поэтому погреш­ ность формулы (1.235) может оказаться больше, если пластическая деформация достаточно развита. На рис. 38 сопоста­ влены зависимости нагрузок в диске с отверстием (пропорциональные квад­ рату угловой скорости) от перемеще­ ний на внутреннем контуре, получен­ ные по формулам (1.234) и (1.235), с за­ висимостями, полученными в резуль­ тате решения интегральных уравнении диска. Погрешность небольшая, при­ чем формула (1.234) дает завышенные, а формула (1.235) заниженные значе-