Механика деформируемого твердого тела
.pdfУравнения теории старения
Пренебрегая упругой деформацией, уравнение теории ползучести можно
постулировать в виде |
σ = |
ϕ(ε) |
|
где функция ϕ и коэффициент β |
1+ βt |
|
|
, вообще говоря, зависят от температуры. |
Содерберг использовал в расчетах на ползучесть следующий вариант теории старения, учитывающий упругую деформацию среды:
ε=σE + Ω(t)σ m.
При больших временах первым слагаемым (упругой деформацией) можно пренебречь по сравнению со вторым (деформацией ползучести). В этом случае диаграмма ползучести описывается уравнением
ε= Ω(t)σ0m
Поэтому график функции получается путем сжатия диаграммы ползучести по оси
ординат. Уравнение кривой релаксации имеет в рамках этой теории вид
σ= (ε0 Ω(t))1m.
Беляев и Малинин применили к решению ряда прикладных задач типа задачи о ползучести вращающегося диска более общее определяющее уравнение теории старения, которое можно рассматривать также с позиции теории течения
& |
1 |
|
& |
& |
ε |
& |
n |
|
|
= |
|
|
|
|
, |
||||
ε |
+ψ σ +ψ σ = |
σ |
σ + Bσ |
|
|||||
|
E |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
Теория течения
В теории течения скорость деформации ползучести представляет собой функцию, зависящую от напряжения, температуры и времени. Если опустить параметр температуры, то определяющее уравнение примет вид
ε&c = f (σ,t) = B(t)σ n .
Здесь n >1 – показатель ползучести, который обычно полагается равным 2 . Впервые вариант этой теории с учетом упругой деформации материала
рассмотрел Дейвенпорт. В его модели постулировалось уравнение
ε& = ε&e +ε&c =σ&E + B(t)σ n ,
Если функция B зависит только от температуры и не зависит от t , то модель допускает преобразование сдвига по времени.
Диаграмма ползучести задается в теории течения уравнением
& |
n |
|
n |
t |
n |
|
: |
∫B(ϑ)dϑ =σ0 |
|||||
ε = Bσ0 |
ε =σ0 E +σ0 |
E + Bσ0 t. |
0
Процесс релаксации напряжений подчиняется обыкновенному дифференциальному уравнению
σ& = −BE σ n
с начальным условием σ0 = Eε0 .
Следовательно, ) −
σ(t) = (σ01−m + (m −1)BEt 1 (1 m) → 0 (t → ∞).
Теория упрочнения
Теория упрочнения при ползучести учитывает зависимость свойств материала от достигнутой деформации. Определяющее уравнение этой теории имеет вид
ε&c = f (σ)g(εc ).
Один из вариантов задания функций f и g , характеризующих свойства материала, состоит в следующем: f (σ) = K eσ A, g(ε) = εα .
При таком задании, функция f (σ) приσ → 0 стремится к отличному от нуля конечному
пределу. Следовательно, скорость деформации ползучести не равна нулю даже в ненапряженном состоянии образца.
В рамках теории упрочнения диаграмма ползучести материала описывается обыкновенным дифференциальным уравнением (ε c ≈ ε):
εαε& = K eσ0 A,
решение которого, удовлетворяющее начальному условию ε(0) = 0 , имеет вид
ε(t) = ((α +1)K eσ0 At)1(α+1).
Кривые релаксации определяются в результате интегрирования дифференциального
уравнения |
|
|
σ α σ& |
|
|
|
ε0 |
− |
= K eσ A |
||
− |
|
E |
|||
|
|
|
E |
|
Интегрируя уравнение с учетом начального условияσ0 = Eε0 , имеем
& |
α+1 σ A |
|
|
1 |
σ |
α |
−ς A |
|
|
|
∫ |
|
|||||
−(σ0 −σ)σ = KE |
e |
, |
t = − |
|
(σ0 −ς) e |
|
dς. |
|
KEα+1 |
|
|||||||
|
|
|
|
|
σ0 |
|
|
|
Задача обратной ползучести
В рамках модели вязкоупругой среды Поинтинга–Томсона определить зависимость деформации от времени при ступенчатом нагружении образца, считая, что на интервале времени (0,2τ) , гдеτ =η E – время релаксации, действовало напряжениеσ0 , которое в момент времени t = 2τ мгновенно упало до нуля.
Реологическая схема модели Поинтинга–Томсона
|
представляет собой последовательное соединение упругого |
|||||
|
элемента с модулем Юнга E0 и схемы Кельвина–Фойхта.На |
|||||
|
этапе нагружения t (0,2τ) , имеем |
|||||
|
|
ε(t) = |
|
σ0 + σ0 (1−e−t τ ). |
||
|
|
|
|
E0 |
E |
|
|
К началу разгрузки t → 2τ −0 деформация образца, |
|||||
|
вычисленная по этой формуле, равна |
|||||
|
|
ε = σ0 + |
σ0 |
(1−e−2 ). |
||
|
|
E0 |
|
E |
|
|
|
В процессе разгрузки деформация упругого элемента с |
|||||
|
модулем мгновенно падает до нуля, поэтому к моменту |
|||||
|
t → 2τ +0 : |
ε =σ0 (1− e−2 ) E . |
||||
Дальнейшее деформирование описывается уравнением |
||||||
|
|
ε(t) = |
σ |
(e2 −1) e−t τ . |
||
|
|
|
0 |
|
||
|
|
|
E |
|||
|
|
|
|
|
||
|
График зависимости ε(t) изображен на рисунке внизу. |
Реологическая схема Поинтинга–Томсона
Диаграмма обратной ползучести
Задача ползучести стержневой системы
Три одинаковых стержня длиной l и сечением S , образующие два правильных треугольника, нагружены силойP в общей вершине. Необходимо определить, как со временем изменяются напряжения и деформации стержней.
Решение будем строить в рамках теории течения при ползучести,
|
учитывая упругие деформации. Напряжение и деформацию в |
|
|
||||||||||||
|
стержнях отметим соответствующими индексами. Уравнение |
|
|
||||||||||||
|
равновесия сил в точке приложения имеет вид |
|
|
|
|
||||||||||
|
1S |
+ 2 2S cosπ |
P : |
1 |
+ |
|
2 |
P . |
|
|
|
Схема нагружения |
|||
|
|
3 |
|
|
|
|
|
S |
|
|
|
||||
В начальный момент времени стержни находятся в упругом состоянии, отсюда |
|||||||||||||||
|
|
σ1 = Eε1, |
σ2 = Eε2 : |
ε1 +ε2 = |
P |
. |
|
|
|||||||
|
|
|
|
|
|||||||||||
|
Из геометрических соображений ε1 = 2ε2 , |
|
|
|
|
ES |
|
|
|||||||
|
При t = 0 |
σ10 = |
2P |
, |
σ20 = |
P |
|
, |
ε10 = |
2P |
, |
ε20 = |
P |
. |
|
|
|
3S |
3S |
|
|
||||||||||
|
|
|
|
|
|
|
3ES |
|
|
3ES |
При t > 0 стержни деформируются в условиях ползучести. По теории течения
|
Кроме того, σ2 = P S −σ1 |
ε&j =σ& j |
E + Bσ 2j . |
( j =1,2) |
|
|
|
|||||||
|
1 |
|
|
2 |
|
2 |
|
|
P |
|
|
2 |
||
|
Справедливо соотношение |
|
|
|
|
|
|
|
||||||
|
|
|
σ&1 |
+ Bσ1 |
= − |
|
|
σ&1 |
+ 2B |
−σ1 |
. |
|||
|
E |
E |
||||||||||||
|
|
|
|
|
|
|
|
S |
|
|
|
Ползучесть приводит к перераспределению напряжений в стержневой системе: напряжение
в центральном стержне снижается, а в боковых – растет.
Характерное время релаксации напряжения |
τ = |
1 |
= |
3S |
обратно пропорционально |
величине условного напряжения P S . |
|
k |
|
2 2BEP |
|
|
|
|
|
|
|
Задача разрушения стержня при ползучести
Определить время, за которое стержень, имеющий начальную длину l0 и поперечное сечение S0 , растянется за счет деформации ползучести до бесконечных размеров под действием постоянной осевой силы P .
В этой задаче упругую деформацию стержня можно не учитывать.Определяющее
уравнение теории течения запишем в виде
ε& = Bσ n .
Для простоты будем считать материал несжимаемым. В этом случае объем стержня не меняется: S l = S0l0, зато меняется его длина и поперечное сечение. Напряжение в стержне также является величиной переменной:
σ = |
P |
= |
Pl |
=σ0 |
l |
. |
S |
S0l0 |
|
||||
|
|
|
l0 |
Деформация за бесконечно малое время |
t равна относительному удлинению |
|||||||
стержня, поэтому |
& |
l |
|
& |
l& |
|
d |
|
|
ε t = |
|
: |
ε = |
|
= |
|
ln l. |
|
|
|
|
ll dt
Искомая зависимость l(t) принимает следующий вид: l = l0 n 1−nBσ0n.
При некоторой достаточно большой длине стержня происходит его разрушение.
Считая длину в момент разрушения бесконечной, по этой формуле можно
определить долговечность и предел длительной прочности – начальное напряжение, под действием которого стержень сохраняет ресурс заданное время :
|
|
|
|
|
|
|
|
|
|
t0 = |
1 |
, |
σ0 |
= n 1 . |
|
|
|
|
|
nBσ0n |
|
|
|
nBt |
|
|
Задача ползучести вращающегося диска
Требуется определить напряженное и деформированное состояние диска радиуса a переменной толщины h(r), вращающегося длительное время с постоянной угловой скоростьюω . Упругой деформацией диска можно пренебречь.
Решение задачи
Запишем уравнение движения диска в проекции на радиальную ось r с учетом центробежной силы инерции,
играющей в данном случае роль массовой силы:
d
dr (rhσr )−hσϕ + ρhω2r2 = 0.
Задача о
вращающемся диске
Здесь σr и σϕ – отличные от нуля радиальное и окружное напряжения.
Уравнения совместности имеют вид
ddrε&ϕ + ε&ϕ −r ε&r = 0.
Для замыкания к системе уравнений необходимо добавить определяющие уравнения
ε&r = Bσin−1 σr −σϕ 2 , ε&ϕ = Bσin−1 σϕ −σr 2 ,
причемσi = σr2 −σrσϕ +σϕ2
Полученная замкнутая система четырех уравнений с четырьмя неизвестными функциями решается с учетом граничного условия σr = 0 на внешнем контуре r = a и условия
ограниченности радиального напряжения в центре диска r = 0 . После этого еще одна
отличная от нуля компонента тензора скоростей деформации ε&z = h& h по оси может быть найдена какε&z = −ε&r −ε&ϕ .
При произвольной заданной зависимости h(r) решение задачи строится численно с привлечением методов решения краевых задач для систем обыкновенных дифференциальных уравнений.
Задача течения вязкопластической среды в трубе
Требуется определить распределение скорости при
собственного веса в вертикальной трубе радиуса a . Заданы плотность среды ρ , предел текучести τs и вязкостьη . На стенках трубы выполняется условие полного прилипания.
Пусть r и z – радиальная и осевая координаты точки. В этой |
|
|
задаче отлично от нуля касательное напряжение τrz =τ |
, остальные |
|
|
||
компоненты тензора напряжений в цилиндрической системе |
|
|
координат тождественно равны нулю. Выполняется |
|
|
|
|
|
дифференциальное уравнение равновесия в проекции на ось z , в |
Течение в трубе |
|
которое входит ускорение свободного падения g : |
|
1 d(rτ)
rdr − ρg = 0.
Отсюда, с учетом ограниченности касательного напряжения при r = 0 следует, что оно линейно зависит от радиуса:τ = ρgr2.
Если a ≤ c ≡ 2τs (ρg) , то напряжение всюду не превосходит предела текучести, и так как на
выполняется условие прилипания, среда не движется. При a >вблизиc оси
трубы возникает застойная зона, движущаяся как жесткое целое. Радиус этой зоны равен |
||||||||
. |
c |
|
|
|
|
|
|
|
Вертикальная скорость vz = v |
в области течения находится из определяющего уравнения |
|||||||
|
& |
|
2 |
2 |
|
|
τs (r −a) |
|
|
|
ρg(r −a |
|
) |
|
|
||
для скорости сдвигаγrz = dv dr: |
v = |
|
− |
. |
||||
|
|
4η |
|
|
η |
|||
|
|
|
|
|
|
|
Такое распределение скорости реализуется только в области течения r > c. В застойной
зоне скорость постоянна и равна v(c) = ρg(c2 −a2 ) −τs (c −a) .
4η η
Задача течения в вискозиметре Шведова
Требуется описать вращательное движение невесомой вязкопластической среды между двумя коаксиальными цилиндрами радиуса a и радиуса b > a , считая, что внутренний цилиндр
неподвижен, а внешний вращается вокруг своей оси с постоянной
угловой скоростью ω , если на стенках цилиндров выполняется условие прилипания.
При вращательном движении отлично от нуля касательное напряжениеτrϕ =τ , которое удовлетворяет дифференциальному
цилиндрической системы координат: dτ |
+ |
2τ |
= 0. |
Вращательное движение |
|
|
dr |
r |
Вблизи внешнего цилиндра может образоваться застойная зона, в которой напряжение не достигает предела текучести, и среда не деформируется.
Напряжение записывается в видеτ =τ0 a 2 r2 . Здесь τ0 - касательное
напряжение на внутреннем цилиндре |
|
|
||||||||
Решение задачи |
|
v(b) |
|
τ0a2 1 |
|
1 |
τs |
b |
||
|
ω = |
b |
= |
2η |
|
− |
|
− |
η ln a . |
|
|
a2 |
b2 |
||||||||
построено полностью. Оно получено в параметрической форме. Задавая величину τ0 в |
||||||||||
пределах от τs до |
τs b2 |
a2, получаем распределение касательного напряжения и скорости в |
||||||||
области a ≤ r ≤ c . |
|
|
|
|
|
|
|
|
|
|
5. Введение в механику разрушений