МЕХАНИКА ДЕФОРМИРУЕМОГО ТВЕРДОГО ТЕЛА
УДК 539.374; 539.389.2
ПРИМЕНЕНИЕ ТЕОРИИ СОБСТВЕННЫХ НАПРЯЖЕНИЙ К ОПИСАНИЮ НЕЛИНЕЙНОГО ДЕФОРМИРОВАНИЯ РАЗНОСОПРОТИВЛЯЮЩИХСЯ МАТЕРИАЛОВ
Б.М. Пахомов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Предложена модель нелинейного деформирования разносопротивляющихся материалов, основанная на схеме разделения обобщенной жесткости на жесткость, определяющую наличие связей между различными направлениями деформирования, и собственную жесткость. Получено новое условие начала пластического течения для материалов, имеющих разные значения пределов текучести при растяжении и сжатии. В пространстве главных напряжений это условие определяет эллипсоид вращения со смещенным относительно начала координат центром. Представленная модель позволяет описать нелинейное поведение таких материалов, которые в упругой области работают как обычные, но имеют разные пределы текучести и разные диаграммы деформирования в неупругой зоне при растяжении и сжатии.
Ключевые слова: схема разделения обобщенной жесткости, условие текучести, деформационная модель, разносопротивляемость.
APPLICATION OF INHERENT STRESS THEORY TO ANALYSIS OF NON-LINEAR DEFORMATION OF MATERIALS WITH MULTIMODULUS BEHAVIOUR
B.M. Pakhomov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]
A model of non-linear deformation of materials is suggested. The model is based on the decomposition of generalized stiffness on the stiffness defining connections between different deformation directions, and self-stiffness. New condition of plastic flow beginning is obtained for materials with different yield points under tension and compression loads. This condition defines in the space of principal stresses an ellipsoid of rotation with rotational center shifted from the beginning of coordinates. The suggested model allows to describe nonlinear behavior of materials working in elastic domain as conventional materials but having different yield stresses and deformation diagrams in non-elastic domain under tensile and compressing loads.
Keywords: decomposition of generalized stiffness, plastic flow condition, deformation model, multimodulus behaviour.
Материалы, по-разному реагирующие на растягивающие и сжимающие нагрузки, условно можно разделить на две группы. Одни материалы можно называть разномодульными. Диаграмма растяжение-сжатие у них возле нулевой точки обычно имеет резко нелинейный
характер. Эти диаграммы хорошо аппроксимируются двумя прямыми, исходящими из нуля под разными углами к координатным осям. Примером материалов с такими свойствами могут служить многие композиционные материалы типа углепластиков, углерод-углеродные композиты и материалы, имеющие пористую структуру [1, 2].
Ко второй группе можно отнести такие материалы, которые до определенных уровней напряжений ведут себя как обычные изотропные линейно-упругие среды, одинаково реагирующие на растягивающие и сжимающие нагрузки, но при этом имеют разные значения пределов текучести и прочности при растяжении и сжатии и разные диаграммы деформирования в неупругой области. К ним можно отнести некоторые металлы и полимеры. В работах [3-5] приведены экспериментальные данные и модели, описывающие поведение таких материалов при сложном напряженном состоянии. В настоящей работе речь будет идти только о материалах второй группы. Будем называть их разносопротивляющимися.
Предлагаемая в настоящей работе модель основана на схеме разделения обобщенной жесткости материала на жесткость, определяющую связи между различными направлениями деформирования, и так называемую собственную жесткость, т.е. обобщенную жесткость материала, лишенного этих связей. Основные положения этой схемы изложены в работах [6, 7].
Под направлениями деформирования в данном случае понимаются растяжение, сжатие по осям х, у, г и сдвиги в плоскостях ху, уг и гх. Математическая формулировка такого разделения может быть представлена на примере записи закона Гука для изотропного, линейно-упругого тела в следующем виде:
= Г0ец + Ье^ц , (1)
где оц — компоненты тензора напряжений; ец — компоненты тензора деформаций; е = ехх + еуу + егг; 5ц — символ Кронекера; г, ] = х, у, г. Параметры в0 и Ь, называемые собственным модулем и модулем связи, в дальнейшем считаются основными характеристиками среды. С техническими характеристиками материала эти параметры связаны очевидными соотношениями
Г = Е Т = Еи
Го — -, Ь — —-:—;-т ,
0 А+У (1 + V)(1 - 2v), Е = Оо(Го + зь) v = _ Ь (2)
во + 2Ь , во + 2Ь'
Здесь Е — модуль упругости; V — коэффициент Пуассона.
Также вводятся некоторые новые понятия: собственные напряжения рц и напряжения связи Ц, определяемые по следующим формулам:
Pij — Go£ij ; (3)
Q — Le,
Рг = \/ 1 Ргз Ргз (4)
— интенсивность собственных напряжений;
и.р = — ргз £ц;
2 (5)
^ = 2 ^
— энергия собственных напряжений или просто собственная энергия ир и энергия связи
Под интенсивностью деформаций в данном случае будем понимать величину €г, определяемую по формуле
€ = \/1 €гз . (6)
Компоненты тензора напряжений при этом равны сумме собственных напряжений и напряжений связи
&гз = Ргз + Я^гз , (7)
а удельная потенциальная энергия равна сумме собственной энергии и энергии связи
и = и + ид. (8)
Легко показать, что условие аддитивности в данном случае выполняется.
Для большей ясности можно провести аналогию между предложенной схемой разделения обобщенной жесткости и классической схемой разделения на жесткость формоизменения и жесткость изменения объема [8, 9]. Собственные напряжения ргз являются аналогами компонент девиатора напряжений; напряжения связи Я — компонент шарового тензора; коэффициент С0 — аналог удвоенного модуля сдвига; Ь — модуля изменения объема; интенсивность рг — интенсивности обычных напряжений; собственная энергия ир — энергии формоизменения.
Понятия собственных модулей, собственных напряженных состояний используются в работах Н.И. Остросаблина, Ю.И. Димитриенко и др. [10-12].
В этих работах, используя чисто математический подход, исследуются свойства тензоров напряжений и деформаций, тензоров коэффициентов упругости и податливостей. Здесь же понятия с определением
"собственный" имеют принципиально другой физический смысл. Так, например, под собственным модулем понимается физическая величина, характеризующая жесткость материала только в данном направлении, т.е. как бы принадлежащая только этому направлению.
В классическом подходе, разработанном для металлов и сплавов, считается, что объем изменяется упруго, а нелинейным считается механизм формоизменения. В предложенной в работах [6, 7] схеме основным постулатом принимается предположение о том, что упруго работают связи между направлениями деформирования, а все нелинейные эффекты заключены в механизме собственной жесткости материала. Отсюда следует вывод о том, что пластическое течение в материале наступает тогда, когда критического значения достигает собственная энергия ир или когда критического значения достигает интенсивность рг, что легко показать, используя соотношения (3)-(5). Очевидно, что это критическое значение для изотропного материала должно равняться пределу текучести при сдвиге т*. Используя формулы (2)-(4), получаем условие текучести через параметры во и Ь
(Г + 4ГоЬ + 6Ь2) (о2х + Оу2у + о2*) -
— 2Ь (2Го + 3Ь) (оххоуу+оуу°гг+оггохх) +
+ 2(Го + 3Ь)2 (оху + о2г + о2х) = 2(Го + 3Ь)2т*2 (9) и через модуль упругости Е и коэффициент Пуассона V
(1 + 2*2) (охх + о2у + о2г) - 2v (2 - V) (
оххоуу + оуу огг + °22 о хх
)+
+ 2(1 + V)2 (оху + о2г + о2х) = 2(1 + V)2т*2. (9')
Это условие инвариантно по отношению к напряженному состоянию и его можно записать через первый и второй /2, инварианты тензора напряжений
(1 + 2^) + 2(1 + V)2/2„ = 2(1 + V)2т*2. (10)
На основе предложенной схемы разделения обобщенной жесткости материала была построена нелинейная модель деформационного типа [7]. При этом были приняты следующие допущения.
1. При рг < т* материал работает упруго с начальным модулем во, а при рг > т* значение собственного модуля в монотонно убывает и является функцией только интенсивности собственных напряжений рг. Это означает, что для всех видов напряженно-деформированного состояния существует единая диаграмма деформирования рг - ег.
Уравнение (1) перепишем в новом виде:
оц = вец + Ье£ц, (11)
где _
д [ С при рг < П, \ С (рг) при Рг > п.
Напряжения р^ будут определяться по формуле
2. Значение модуля связи Ь остается постоянным.
3. Разгрузка происходит упруго с начальным модулем С0.
(12)
Обратив соотношения (11), получим, что полные ец, упругие е(е и
пластические едеформации определяются через конечные значения напряжений по следующим формулам:
1
L
Je)
G(G + 3L) L
aSj;
Go
Go (Go + 3L)
a8j;
e(- = ,___
1
G
i
Go
o,jj — L
G(G + 3L) Go (Go + 3L)
а5.
U
(13)
где а = а гг.
Особое отличие предложенного подхода от классического заключается в том, что в условии текучести (10) присутствует первый инвариант тензора напряжений 11а и при неупругом деформировании происходит пластическое изменение объема Де(р), пропорциональное а,
Де(-) = -0-, K-
(14)
где K- =
(Go + 3L) (G + 3L)
будем называть модулем пластического
Со — С изменения объема.
Из (14) следует, что знаки Де(р) и а всегда совпадают. В процессе деформации тело всегда накапливает упругую энергию, за счет которой оно возвращается в естественное состояние после того, как силы, вызвавшие деформацию, исчезают. Поэтому удельная потенциальная энергия и в общем случае должна быть положительной.
Исходя из (3), (5) и (8), и выражается формулой
1
U = 1 Ge 2
i3eij + 2 Le
т.е. и является положительно определенной квадратичной формой относительно компонент тензора деформаций. Применяя критерий Сильвестра, легко установить, что параметры й0 и Ь должны удовлетворять условиям С0 > 0 и С0 + 3Ь > 0 одновременно.
1
1
1
Если считать, что модуль упругости Е всегда должен быть больше нуля, то, используя выражения (2), получаем границы изменения коэффициента Пуассона для изотропного линейно-упругого материала
1
-1 < V < -.
< < 2
У некоторых разносопротивляющихся материалов (например, у серых чугунов) наблюдается разрыхление, т.е. пластическое увеличение объема при сжимающих нагрузках. Если не отказываться от условия положительности и, оставаясь в рамках описанной ранее модели деформационного типа, то описать "разбухание" материала при сжатии невозможно. Развитие данного подхода применительно к таким материалам возможно, но не является темой настоящей работы.
При разгрузке напряжения оц, Q и рц будут равны:
Я = + Ь Де(р) = -^о + --^ о , Д е(р);
+ 3Ь 1 + V (1 + V)(1 - 2v) '
Go + oL
v Ei и
Pij = aj - гт"— a^ - T;—:—Г71-^ Ae(p).
(15)
1 + V ц (1 + V) (1 - 2v)
Подставив выражения для собственных напряжений в условие начала текучести = т*и, где т*и — значение предела текучести на сдвиг, достигнутое к началу разгрузки, получим условие текучести для повторного нагружения
(Г + 4ГоЬ + 6Ь2) (о2х + оу2у + о2г) -
- 2Ь (2Г о+
3Ь) (оххоуу+оууогг+оггохх
)+
+ 2(Го + 3Ь)2 (оху + + ог2х) - 2Го (Го + 3Ь) ЬоДе(р) +
+ 3Ь2(Го + 3Ь)2(Де(р))2 = 2(Го + 3Ь)2т*2и (16) или через модуль упругости Е и коэффициент Пуассона V
(1 + 2v2) (аХж + а^ + а^) - 2v (2-v) (ажжата+атаa^ + a^ажж) +
3E 2v 2
(1—2V)2
q p2 2
+ 2(1 + v)2 (aXy + a^ + aL) - 2£vaAe(p) + ^(p))2 =
= 2(1 + V )2т*2и. (17)
В пространстве главных напряжений о1, о2, о3 уравнения (16) и (17) определяют эллипсоид вращения, у которого наибольший из главных диаметров равнонаклонен к координатным осям, а координаты
центра выражаются через остаточное изменение объема
Ь(д0р + 3Ь)Де (р);
ЕР
а1 s = а2 s = a3s = as = -^г
(18)
Go
Os =
-Де(-).
(1 — )2
Полудлины главных диаметров ё1, ё2, ё3 определяются формулами
¿1 = ¿2 = \[2т*и] = У2 (1+ V) ¿3 = (1 — 2v)
Легко показать, что условиями существования этого эллипсоида, так же как и эллипсоида, который определяется условиями (9) и (9'), являются полученные ранее ограничения, накладываемые на параметры С0, Ь и V.
Из (14) следует, что модуль пластического изменения объема стремится при С ^ 0 к величине
кр = ^Е
т1П (1 — 2v )2'
Тогда с учетом (14) и (18) максимальное значение перемещения центра эллипсоида не может быть больше (1/3)а.
Таким образом, согласно предложенной деформационной модели, если происходит пластическое изменение объема, то поверхность пластичности изотропно расширяется и одновременно перемещается вдоль своего главного диаметра, при этом материал приобретает анизотропию свойств.
Используя (4), (6) и (13) можно показать, что при полной разгрузке возникнут остаточные "внутренние" напряжения ргхх, ргуу, ртгг и . Значения этих остаточных напряжений будут определяться формулой
рХх = ргуу = рг22 = —СТ = —ЬДе(р). (19)
Здесь надо сказать, что любая теория, допускающая неупругое изменение объема и пользующаяся понятием поверхности текучести, обязана описывать этот эффект без введения дополнительных предположений о существовании ориентированных микронапряжений или еще каких-либо предположений. И хотя правильное количественное описание анизотропного упрочнения реальных материалов невозможно исключительно только с помощью предлагаемой модели, нельзя не учитывать вклад в перемещение поверхности текучести в пространстве напряжений, который может внести сколько-нибудь заметное пластическое изменение объема.
Далее рассмотрим случай, когда материал сначала подвергался одноосному растяжению до растягивающего напряжения а^ причем в
материале появлялись пластические деформации, а затем нагруже-ние проводилось в противоположном направлении. Согласно условиям (16) и (17) предел текучести при сжатии ос должен выражаться через о формулами
(55 + б5оЬ + 6Ь2) 5 + Ь(52 + 12ёсЬ + 18Ь2)
ос =--—-—=2-=-о;
(в + 3Ь) (50 + 4воЬ + 6Ь2)
_ (1 - 2^) (1 + 2^ - 2^2) ф + V(1 + 8^ - 2^2)
0с = (ГГ^Т^Т—2^ГфГзй
Параметр ф = характеризует степень развития пластических
деформаций.
Здесь можно сделать предположение о том, что разносопроти-вляющиеся материалы второй группы — это материалы, анизотропия свойств которых появляется в процессе изготовления образца, когда при пластическом деформировании появляется пластическое изменение объема и возникают внутренние напряжения, определяющие эту анизотропию, т.е. причиной различия свойств у некоторых материалов при растягивающих и сжимающих нагрузках могут являться не поры и микротрещины, как у большинства разномодульных материалов, а остаточные, внутренние напряжения, вызванные структурными изменениями при неупругом деформировании.
Применив описанную ранее схему, построим модель деформирования разносопротивляющегося материала, включающую в себя предысторию его создания.
Представим, что сначала материал пластически деформировался, например, под действием гидростатического давления для материалов, у которых предел текучести при сжатии ос больше предела текучести при растяжении о, или под действием всестороннего растяжения, если ос < о. Пусть после этого нагрузка была снята полностью и при этом в материале образовалось остаточное изменение объема ДеЦ^, которое определяется по формуле (20), где р и в — среднее напряжение и собственный модуль, достигнутые к моменту начала разгрузки:
« = 3Р (Со- 5) (20)
Ц (5о + 3Ь) (5 + 3Ь)' (20)
Интенсивность собственных напряжений в момент начала разгрузки обозначим т*Ц. Согласно (3) и (5) она может быть вычислена по формуле
т* и — А т и-
2 в + 3Ь
При этом в материале должны возникнуть остаточные собственные
напряжения р^ж, р^, р!^ и остаточные напряжения связи .
или через E и v
(21)
Условиями пластичности для случая повторной нагрузки будут условия (16) и (17).
Последовательно применив эти условия для случаев одноосного растяжения и одноосного сжатия, получим выражения для неизвестных параметров Де!р) и тт через известные из экспериментов характеристики разносопротивляющегося материала а г и ас:
Деи.) = С2 + 4С0Ь + 6Ь (а, — ас)"
и 2С0Ь(С0 + 3Ь) ^ ' '
Де? = ^ <а. — ^ 2
2 С2 + 4СЬ + 6Ь2 3 (С + 4СЬ + 6Ь2)' >2
т*и = 2(С + 3Ь)2 агас + 8 С2(С + 3Ь)2 (а — ас) ;
2 1 + 2v2 3 (1 + 2v2)2 2
т*и = агас + 8 2^2(а — ас).
(22)
Подставив в (16) и (17) соответственно выражения (21) и (22), получим условия текучести разносопротивляющегося материала через параметры аг и ас
2,2,2 2L(2Go + 3L) , s
axx + ayy + — =2-=-- (oxx0yy+ayyazz + axx) +
Go + 4GoL + 6L2
2 (Go + 3L)2
{°ly + ^z + aL) — (at — ac) (oxx + ayy+azz) = °t°c
1 -=2 „ \ x,y 1 uyz 1 uzxJ t v cj \ "xx
Go + 4GoL + 6L2 y y ; и через E и v
2,2,2 2v (2 — V) , , , X ,
ахх + ауу + ^--1 + 2^ (аххауу + ауу аzz + ^ ахх) +
2(1 + V)2 ( 2
1 + 2V2 \аху + ^ + uzx) — (а — °с) (ахх
(23)
Из (23) следует, что предел текучести разносопротивляющегося материала при чистом сдвиге т^с должен выражаться через аг и ас по формуле
/(ТГ^2)
т"с = агас. (24)
2(1 + V)- 222 + Т . г,7<2 [°xy + 0yz + °zx) — (at — ac) (oxx + ayy + azz) = at°c.
Проведем сравнение полученного условия с известными критериями предельного состояния для разносопротивляющихся материалов [13,14].
Критерий Шлейхера определяет предельное состояние по значению полной удельной потенциальной энергии деформации и включает в себя зависимость от среднего напряжения
аЖж + + а^ — + аота^ + а^ажж) +
+ 2(1 + V) (аЖу + а^ + а^) — — ас) (ажж + аот + а^) = ага^
(25)
Критическое значение касательного напряжения при сдвиге равно
= I а^ас У2(1 + V).
Особенностью данного критерия является то, что в пространстве главных напряжений условию (25) соответствует эллипсоид вращения со смещенным центром относительно координат с главным диаметром, одинакового наклоненным к координатным осям, т.е. это выпуклая и замкнутая поверхности. Достоинством можно считать и то, что условие (25) включает в себя только две независимые константы а4 и ас, определяемые из простых опытов на одноосные растяжение и сжатие. Также представляется логичным тот факт, что условие (25), как и условие (23), содержит коэффициент Пуассона V, так как в упругой области деформирование среды определяется упругими соотношениями. Предлагаемый в настоящей работе критерий (23) получится, если из полной удельной потенциальной энергии деформирования выделить собственную энергию.
Критерий Мизеса-Шлейхера [15] является модифицированным критерием Хубера-Мизеса
аг + ва = к, (26)
где а» — интенсивность напряжений; в — коэффициент внутреннего трения, к — сцепление, хорошо подходит для описания текучести хрупких металлов и горных пород. В этом критерии три независимых величины а4, ас и т*4с, получаемые из эксперимента.
Критерий П.П. Баландина, так же как и критерий Мизеса - Шлейхера, в качестве основы использует удельную потенциальную энергию формоизменения, связанную со средним напряжением параметров:
аЖж + + а22 — (аЖЖауу + ауу а22 + а22°ЖЖ) +
+ 3 (аЖу + а^ + а^ж) — (о4 — ас) (ажж + аот + а^) = а4ас. (27)
Величина т^с через иг и ис определяется по формуле
т* tc
OtOc
3
При и = ис условие (27) так же, как и условие (25), переходит в условие Мизеса. По форме это условие совпадает с условием (23), только здесь вместо собственной энергии используется энергия изменения формы. Критерий П.П.Баландина получается из условия (23) при стремлении коэффициента Пуассона к 0,5.
В пространстве главных напряжений критерию П.П. Баландина соответствует разомкнутая поверхность — параболоид вращения, что сильно уменьшает его универсальность.
Критерий И.Н. Миролюбова включает в себя по сравнению с критерием П.П. Баландина еще и квадратичную зависимость от среднего напряжения:
К + ис)2
8atac
[(oxx - avv)2 + (avv - ozz)2 + (ozz - vxxf +
+ 6 (Oxy + Oyz + ^x) -
(Vt - ^c)2 4ütOc
(Oxx+Oyy + Ozz )
- (vt - oc)(oxx + Oyy + O-z-z)= OtOc. (28)
Здесь для т* 1с имеем
tc
OtOc
а/3 V и г + и с
Предельная поверхность, описываемая условием (28), — однополосный гиперболоид, т.е. разомкнутая поверхность с отрицательной гауссовой кривизной, что является большим недостатком. Более универсальным является критерий Ю.И. Ягна:
OtOc t 2
27^ (°2
2Т* tc
(Olx + Oly + OL) + 2 (o2xy + °2z + OL) +
+ - иИ^) (ихх+иуу)2- (иг - ис) (ихх+иуу+иzz) =игис.
Этот критерий содержит три независимых параметра т^с, иг и ис. Меняя соотношения между ними, из него можно получить и критерий удельной энергии формоизменения, и критерий П.П. Баландина, и критерий И.Н. Миролюбова, но нельзя получить предлагаемый нами критерий (23).
В работе [13] представлен критерий
3/2 = [OtOc- (Ot - Oc) Ii] X
x < 1 - I 1 - ^
OtOc
X
1 - !>( /2
2 V 3
-3/2
, (29)
где Д — первый инвариант тензора напряжений; /2 и /3 — второй и третий инварианты девиатора напряжений, а И.И. Гольденблат и В.А. Копнов предложили критерий
1 ( 1 1 \ , , ч \ ^х+^да + + ^2 + , ----(^хж + ^уу + ^22) + \ -—2--+
2 °с) У 2^
+
i(± + М 2_ 1
4 V ^ ' °с) 2r2ic
1/2
(Ожж + ^yy + )2 > =1. (30)
Критерии (29) и (30) имеют три независимых параметра r*ic, at и ас, определяемых из экспериментов на одноосные растяжение и сжатие и чистый сдвиг, и при выполнении условий at = ас и r*ic = at/\/3 приводят к условию предельного состояния по теории формоизменения.
В дальнейшем многие авторы предлагали различные соотношения в качестве критерия предельного состояния материала, но все они использовали в качестве основы в этих соотношениях выражения для энергии формоизменения, так что при применении к обычному изотропному материалу, одинаково сопротивляющемуся растяжению и сжатию, из этих критериев получается условие Хубера-Мизеса. Нам не удалось найти работы, где содержались бы результаты, приводящие к условию (23). Обширный обзор литературы по критериям предельного состояния представлен в работе [13].
После того, как мы установили вид условия начала пластического течения, можно перейти к рассмотрению вопроса о построении модели нелинейного деформирования изотропного разносопротивля-ющегося материала. За основу берем описанную ранее модель деформационного типа, т.е. модель, в которой связь устанавливается между конечными значениями напряжений и деформаций.
Процесс деформирования представим проходящим в два этапа. Второй этап соответствует реальному нагружению разносопротивля-ющегося материала (все параметры, относящиеся к реальному материалу, имеют индекс 2). Первый этап соответствует неупругому деформированию гипотетического (изначально одинаково сопротивляющегося растяжению и сжатию) изотропного материала под действием всестороннего растяжения, если предел текучести на растяжение at разносопротивляющегося материала больше предела текучести на сжатие ас, либо под действием всестороннего сжатия, если at < ас. Гипотетическому материалу соответствуют параметры с индексом 1.
Упругие характеристики материала на первом и втором этапах считаются одинаковыми.
Для построения модели необходимо знать значение параметров at и ас, а также иметь диаграмму деформирования разносопротивляюще-
1
гося материала при чистом сдвиге Ttc--Ytc, которая для изотропного
2
материала не зависит от направления кручения. Если диаграмма известна, то известно и значение предела текучести при сдвиге т^е. На этом этапе необходимо провести проверку условия (24). Если условие (24) выполняется с достаточной точностью, то можно предположить, что применение данного подхода окажется корректным. При большой погрешности применение данной модели будет проблематичным.
Далее необходимо построить диаграмму зависимости интенсивности собственных напряжений р от интенсивности деформаций е для гипотетического материала. Это нетрудно сделать для каждой точки
диаграммы. Для чистого сдвига можно записать выражения для гипо-
- (1) -(1)
тетического материала для интенсивностей р\ и е\ , которые реализовались бы, если бы действительно существовала описанная ранее предыстория пластического деформирования
р(1) = </2 Ь2(Де%)2 + т-
г—--(31)
е = ^М'«)2 + 4 ^.
Деформации гипотетического материала при повторном упругом деформировании определяются соотношением
е3 = ^ иц — =—=-и8ц +—Де(?\8ц,
3 во 3 Со(Оо + 3Ь) 3 3 (1) 3
а при неупругом —
еЦ = 1 из -„^-и8г,, + -Де{р1Аз. (32)
3 в 3 в (в + 3Ь) 3 3 !(1) 13 v 7
Деформации реального разносопротивляющегося материала определяются по формуле
j = Oj -o813 . (33)
j G G(G + 3L) j V 7
Неизвестные параметры: деформации и собственный модуль в находятся из решения уравнений (31)-(33), совместно с условием, свя-
(1) (1)
зывающим параметры р\ и . Это решение находится методом последовательных приближений.
С помощью предложенной модели была проведена обработка некоторых экспериментальных данных, приведенных в работе [16]. На рисунке приведены диаграммы деформирования среднеуглеродистой стали, полученные при пропорциональном нагружении трубчатых образцов осевой силой и внутренним давлением. Сплошные линии соответствуют эксперименту, штриховые — расчету. Кривые 1 соответствуют одноосному растяжению, кривая 2 — одноосному сжатию,
кривая 3 — чистому сдвигу, кривые 4, 5 — двухосному растяжению при a/ai = 0,577; 0,667. Здесь ai — интенсивность напряжений. Диаграммы представлены в координатах — обычные интенсивности напряжений ai и деформаций ei. За опорную кривую, по которой проводились описанные ранее расчеты, бралась кривая для чистого сдвига.
Проведенные расчеты для аналогичных экспериментальных данных показали, что предложенная модель может быть использована для описания нелинейного поведения некоторых разносопротивляющихся материалов. Очевидно, что возможность использования данного подхода должна анализироваться в каждом конкретном случае.
Автор благодарит профессора МГТУ им. Н.Э. Баумана В.С. Зарубина за полезные замечания.
ЛИТЕРАТУРА
1. Мясников В.П., Олейников А.И. Основы механики гетерогенно-сопротивляющихся сред. Владивосток: Дальнаука, 2007. 172 с.
2. Мэтьюз Ф., Ролингс Р. Композитные материалы. М.: Техносфера, 2004. 408 с.
3. Максимов Р.Д., Плуме Э.З., Янсонс Ю.О. Сравнительное исследование механических свойств термореактивного полимера при растяжении и сжатии // Механика композитных материалов. 2005. Т. 41. № 5. С. 633-650.
4. Бессонов Д.Е., Зезин Ю.П., Ломакин Е.В. Разно сопротивляемость зернистых композитов на основе ненасыщенных полиэфиров // Изв. Саратовского ун-та. Сер. Математика. Механика. Информатика. 2009. Т. 9. Вып. 4. Ч. 2. С. 9-13.
5. Бессонов Д.Е., Ершова А.Ю., Зезин Ю.П., Мартиросов М.И., Рыбинский Л.Н. Экспериментальное исследование деформирования и разрушения зернистых композитов на основе полиэфирных смол // Механика композиционных материалов и конструкций. 2008. Т. 14. № 1. С. 111-125.
6. Пахомов Б.М. Условие пластического течения, включающее коэффициент Пуассона // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2014. № 2. С. 15-27.
7. Пахомов Б.М. Применение теории собственных напряжений к описанию нелинейного деформирования металлов и сплавов // Инженерный журнал: наука и инновации. 2013. Вып. 7. URL: http://engjournal.ru/catalog/machin/rocket/854.html
8. Хажинский Г.М. Деформирование. Разрушение. Надежность. Задачи деформирования и разрушения стали. Методы оценки прочности энергетического оборудования и трубопроводов / М.: Эдиториал УРСС, 2014. 544 с.
9. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 512 с.
10. Остросаблин Н.И. Анизотропия и общие решения уравнений линейной теории упругости. Дисс. ... д-ра техн. наук. Тула, 2000. 215 с.
о,-, МПа
Диаграмма деформирования сред-неуглеродистой стали
11. Димитриенко Ю.И. Нелинейная механика сплошной среды. М.: ФИЗМАТЛИТ. 2009. 624 с.
12. Димитриенко Ю.И.Тензорный анализ. Т. 1. Механика сплошной среды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 463 с.
13. Божанов П.В. Задачи пластического деформирования тонких пластинок из ди-латирующих разносопротивляющихся материалов. Дисс. ... канд. техн. наук. Тула, 2004. 233 с.
14. Потапова Л.Б., Ярцев В.П. Механика материалов при сложном напряженном состоянии. М.: Машиностроение, 2005. 244 с.
15. Коврижных А.М. Уравнения плоского напряженного состояния при условии пластичности Мизеса-Шлейхера // Прикладная механика и техническая физика. 2004. Т. 45. № 6. С. 144-153.
16. Трещев А.А. Вариант подхода к построению определяющих соотношений разносопротивляющихся материалов и использование его при расчете элементов конструкций. Дисс. ... д-ра техн. наук. Тула, 1995. 501 с.
REFERENCES
[1] Myasnikov V.P., Oleynikov A.I. Osnovy mekhaniki geterogenno-soprotiv-lyayushchikhsya sred [Grounds of mechanics of heterogeneously resisting media]. Vladivostok, Dal'nauka Publ., 2007. 172 p.
[2] Met'yuz F., Rolings R. Russ. ed.: Kompozitnye materialy [Compositematerials]. Moscow, Tekhnosfera Publ., 2004. 408 p.
[3] Maksimov R.D., Plume E.Z., Yansons Yu.O. Comparative research of thermosetting polymer mechanical properties at tension and compression.Mekh. Kompozitn. Mater. [Mechanics of Composite Materials], 2005, vol. 41, no. 5, pp. 633-650 (in Russ.).
[4] Bessonov D.E., Zezin Yu.P., Lomakin E.V. Multimodal resistance ofgranular composites based on unsaturated polyesters.Izv. Saratov. Univ. Mat. Mekh. Inform., 2009, vol. 9, iss. 4, part. 2, pp. 9-13 (in Russ.).
[5] Bessonov D.E., Ershova A.Yu., Zezin Yu.P., Martirosov M.I., Rybinskiy L.N. Experimental research of deforming and destruction of granular composites based on polyester resins. Mekhanika kompozitsionnykh materialov i konstruktsiy [Journal on Composite Mechanics and Design], 2008, vol. 14, no. 1, pp. 111-125 (in Russ.).
[6] Pakhomov B.M. Plastic yielding condition including Poisson's ratio. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2014, no. 2, pp. 15-27 (in Russ.).
[7] Pakhomov B.M. Application of the inherent stresses theory to nonlinear deforming of metals and alloys. Jelektr. nauchno-tehn. Izd. "Inzhenernyj zhurnal: nauka i innovacii" MGTU im. Baumana [El. Sc.-Techn. Publ. "Eng. J.: Science and Innovation" of Bauman MSTU], 2013, no. 7(19). Available at: URL: http://engjournal.ru/catalog/machin/rocket/854.html (in Russ.).
[8] Khazhinskiy G.M. Deformirovanie. Razrushenie. Nadezhnost'. Zadachi deformirovaniya i razrusheniya stali. Metody otsenki prochnosti energeticheskogo oborudovaniya i truboprovodov [Deformation. Fracture. Reliability. Problems of Deformation and Fracture of Steel. Methods for Assessing the Strength of Power Equipment and Pipelines.]. Moscow, Editorial URSS Publ., 2014. 544 p.
[9] Zarubin V.S., Kuvyrkin G.N. Matematicheskie modeli mekhaniki i elektrodinamiki sploshnoy sredy [Mathematical Models of Continuum Mechanics and Electrodynamics]. Moscow MGTU im. N.E. Baumana Publ., 2008. 512 p.
[10] Ostrosablin N.I. Anizotropiya i obshchie resheniya uravneniy lineynoy teorii uprugosti [Anisotropy and general solutions of linear elasticity theory]. Diss. dokt. tekhn. nauk [Dr. tehn. sci. diss.]. Tula, 2000. 215 p.
[11] Dimitrienko Yu.I. Nelineynaya mekhanika sploshnoy sredy [Nonlinear mechanics of continuum]. Moscow, Fizmatlit Publ., 2009. 624 p.
[12] Dimitrienko Yu.I. Tenzornyy analiz. T. 1. Mekhanika sploshnoy sredy [Tensor analysis. Vol. 1. Continuum mechanics. Moscow, MGTU im. N.E. Baumana Publ., 2011. 463 p.
[13] Bozhanov P.V. Zadachi plasticheskogo deformirovaniya tonkikh plastinok iz dilatiruyushchikh raznosoprotivlyayushchikhsya materialov [Problems of plastic deformation of thin plates of dilatant materials with multimodal resistance]. Diss. kand. tekhn. nauk [Cand. tehn. sci. diss.]. Tula, 2004. 233 p.
[14] Potapova L.B., Yartsev V.P. Mekhanika materialov pri slozhnom napryazhennom sostoyanii [Mechanics of materials in the complex stressed state]. Moscow, Mashinostroenie Publ., 2005. 244 p.
[15] Kovrizhnykh A.M. Plane stress equations for the Von Mises-Schleicher yield criterion. Prikl. Mat. iMekh. [J. Appl. Math. Mech.], 2004, vol. 45, no. 6, pp. 144-153 (in Russ.).
[16] Treshchev A.A. Variant podkhoda k postroeniyu opredelyayushchikh sootnosheniy raznosoprotivlyayushchikhsya materialov i ispol'zovanie ego pri raschete elementov konstruktsiy [Version of approach to obtain the defining relationships for materials with multimodal resistance, and application of this approach for designing the construction elements]. Diss. dokt. tekhn. nauk [Dr. tehn. sci. diss.]. Tula, 1995. 501 p.
Статья поступила в редакцию 16.05.2014
Пахомов Борис Максимович — канд. техн. наук, доцент кафедры "Космические аппараты и ракеты-носители" МГТУ им. Н.Э. Баумана. Автор более 10 научных работ в области механики нелинейного деформирования материалов. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Pakhomov B.M. — Cand. Sci. (Eng.), assoc. professor of "Spacecrafts and Space Vehicles" department of the Bauman Moscow State Technical University. Author of more than 10 publications in the field of mechanics of nonlinear deforming of materials. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.