Критерий прочности Друкера — Прагера

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

Назван по именам Даниэля Друкера и Прагера, разработавшим эту модель в 1952 году.

Формулировка

Критерий описывается следующей формулой:

$\sqrt{J_2}=A+B{I}_1$

где $I_1$ — первый инвариант тензора напряжений, а $J_2$ — второй инвариант девиатора тензора напряжений. Константы $A,B$ определяются экспериментально.

В терминах эквивалентных напряжений (или напряжений по Мизесу) и гидростатических напряжений, критерий Друкера — Прагера может быть записан как:

${\sigma }_e=a+b{\sigma }_m$

где ${\sigma }_e$ — эквивалентное напряжение, ${\sigma }_m$ — гидростатическое напряжение, и $a,b$ константы материала. Критерий Друкера — Прагера, выраженный в координатах Хейга — Вестергаарда следующим образом:

$\frac{1}{\sqrt{2}}\rho -\sqrt{3}B\xi =A$

Поверхность текучести Друкера — Прагера есть сглаженная версия поверхности текучести Мора — Кулона.

Выражения для A и B

Модель Друкера — Прагера может быть записана в терминах главных напряжений:

$\sqrt{\frac{1}{6}[({\sigma }_1-{\sigma }_2{)}^2+({\sigma }_2-{\sigma }_3{)}^2+({\sigma }_3-{\sigma }_1{)}^2]}=A+B({\sigma }_1+{\sigma }_2+{\sigma }_3).$

Если ${\sigma }_t$ — предел прочности при одноосном растяжении, критерий Друкера — Прагера означает:

$\frac{1}{\sqrt{3}}{\sigma }_t=A+B{\sigma }_t.$

Если ${\sigma }_c$ предел прочности при одноосном сжатии, критерий Друкера — Прагера означает:

$\frac{1}{\sqrt{3}}{\sigma }_c=A-B{\sigma }_c.$

Решая эти 2 уравнения, получаем

$A=\frac{2}{\sqrt{3}}\left(\frac{{\sigma }_c{\sigma }_t}{{\sigma }_c+{\sigma }_t}\right);B=\frac{1}{\sqrt{3}}\left(\frac{{\sigma }_t-{\sigma }_c}{{\sigma }_c+{\sigma }_t}\right).$

Одноосный асимметричный коэффициент

Различные одноосные критерии прочности на растяжение и сжатие были предсказаны с помощью модели Друкера — Прагера. Одноосный асимметричный коэффициент для модели Друкера — Прагера:

$\beta =\frac{{\sigma }_{\mathrm{c}}}{{\sigma }_{\mathrm{t}}}=\frac{1-\sqrt{3}B}{1+\sqrt{3}B}.$

Выражение в терминах угла трения и когезии

Поскольку поверхность текучести Друкера — Прагера — сглаженная версия поверхности текучести Мора — Кулона, то он часто выражается в терминах когезии ($c$) и угла внутреннего трения ($\varphi$), которые используются в теории Мора — Кулона. Если допустить, что поверхность текучести Друкера — Прагера описана около поверхности текучести Мора — Кулона, тогда выражения для $A$ и $B$ следующие:

$A=\frac{6c\cos \varphi }{\sqrt{3}(3+\sin \varphi )};B=\frac{2\sin \varphi }{\sqrt{3}(3+\sin \varphi )}$

Если поверхность текучести Друкера-Прагера вписана в поверхность текучести Мора-Кулона, то

$A=\frac{6c\cos \varphi }{\sqrt{3}(3-\sin \varphi )};B=\frac{2\sin \varphi }{\sqrt{3}(3-\sin \varphi )}$

Модель Друкера — Прагера для полимеров

Модель Друкера — Прагера используется для моделирования таких полимеров, как полиформальдегид и полипропилен. Для полиформальдегида критерий прочности есть линейная функция от нагрузки. Однако, для полипропилена наблюдается квадратичная зависимость от нагрузки.

Модель Друкера-Прагера для пен

Для пен модель GAZT использует:

$A=\pm \frac{{\sigma }_y}{\sqrt{3}};B=\mp \frac{1}{\sqrt{3}}\left(\frac{\rho }{5{\rho }_s}\right)$

где ${\sigma }_y$ — критическое напряжение для разрушения при растяжении или сжатии, $\rho$ — плотность пены, и ${\rho }_s$ — плотность базового материала(из которого получена пена).

Выражения для изотропной модели Друкера — Прагера

Критерий Друкера — Прагера также может быть использован в альтернативной формулировке:

$J_2=(A+B{I}_1{)}^2=a+b{I}_1+c{I}_1^2.$

Критерий прочности Дешпанде — Флека

Критерий прочности Дешпанде — Флека для пен имеет форму приведенного выше уравнения. Параметры $a,b,c$ для критерии Дешпанда-Флека

$a=(1+{\beta }^2){\sigma }_y^2,b=0,c=-\frac{{\beta }^2}{3}$

где $\beta$ -параметр, определяющий форму поверхности текучести, а ${\sigma }_y$ предел прочности на растяжение или сжатие.

Анизотропный критерий прочности Друкера — Прагера

Анизотропная форма критерия прочности Друкера — Прагера совпадает с критерием прочности Лю — Хуана — Стаута. Этот критерий прочности выражен в обобщенном критерии текучести Хилла:

$\begin{array}{cc}f:=& \sqrt{F({\sigma }_{22}-{\sigma }_{33}{)}^2+G({\sigma }_{33}-{\sigma }_{11}{)}^2+H({\sigma }_{11}-{\sigma }_{22}{)}^2+2L{\sigma }_{23}^2+2M{\sigma }_{31}^2+2N{\sigma }_{12}^2}\\ & +I{\sigma }_{11}+J{\sigma }_{22}+K{\sigma }_{33}-1\le 0\end{array}$

Коэффициенты $F,G,H,L,M,N,I,J,K$ есть:

$\begin{array}{cc}F=& \frac{1}{2}[{\mathrm{\Sigma }}_2^2+{\mathrm{\Sigma }}_3^2-{\mathrm{\Sigma }}_1^2];G=\frac{1}{2}[{\mathrm{\Sigma }}_3^2+{\mathrm{\Sigma }}_1^2-{\mathrm{\Sigma }}_2^2];H=\frac{1}{2}[{\mathrm{\Sigma }}_1^2+{\mathrm{\Sigma }}_2^2-{\mathrm{\Sigma }}_3^2]\\ L=& \frac{1}{2({\sigma }_{23}^y{)}^2};M=\frac{1}{2({\sigma }_{31}^y{)}^2};N=\frac{1}{2({\sigma }_{12}^y{)}^2}\\ I=& \frac{{\sigma }_{1c}-{\sigma }_{1t}}{2{\sigma }_{1c}{\sigma }_{1t}};J=\frac{{\sigma }_{2c}-{\sigma }_{2t}}{2{\sigma }_{2c}{\sigma }_{2t}};K=\frac{{\sigma }_{3c}-{\sigma }_{3t}}{2{\sigma }_{3c}{\sigma }_{3t}}\end{array}$

где

${\mathrm{\Sigma }}_1:=\frac{{\sigma }_{1c}+{\sigma }_{1t}}{2{\sigma }_{1c}{\sigma }_{1t}};{\mathrm{\Sigma }}_2:=\frac{{\sigma }_{2c}+{\sigma }_{2t}}{2{\sigma }_{2c}{\sigma }_{2t}};{\mathrm{\Sigma }}_3:=\frac{{\sigma }_{3c}+{\sigma }_{3t}}{2{\sigma }_{3c}{\sigma }_{3t}}$

и ${\sigma }_{ic},i=1,2,3$ пределы прочности при одноосном сжатии по трем главным направлениям анизотропии, ${\sigma }_{it},i=1,2,3$ пределы прочности при одноосном растяжении, и ${\sigma }_{23}^y,{\sigma }_{31}^y,{\sigma }_{12}^y$ пределы прочности при чистом сдвиге. Выше было допущено, что значения ${\sigma }_{1c},{\sigma }_{2c},{\sigma }_{3c}$ положительные, а ${\sigma }_{1t},{\sigma }_{2t},{\sigma }_{3t}$ — отрицательные.

Критерий текучести Друкера

Критерий Друкера — Прагера не должен вступать в противоречие с более ранним критерием Друкера, который независим от нагрузок ($I_1$). Критерий Друкера имеет запись

$f:=J_2^3-\alpha J_3^2-k^2\le 0$

где $J_2$ — второй инвариант девиатора тензора напряжения, $J_3$ — третий инвариант девиатора тензора напряжения, $\alpha$ — константа, находящаяся между −27/8 и 9/4 (чтобы поверхность текучести была выпуклой), $k$ — константа, меняющаяся в зависимости от $\alpha$. Для $\alpha =0$, $k^2=\frac{{\sigma }_y^6}{27}$, где ${\sigma }_y$ критерий прочности при одноосном растяжении.

Анизотропный критерий Друкера

Анизотропная версия критерия текучести Друкера — критерий текучести Казаку — Барлата, который имеет вид

$f:=(J_2^0{)}^3-\alpha (J_3^0{)}^2-k^2\le 0$

где $J_2^0,J_3^0$ — обобщенные формы девиатора тензора напряжения, определенные как:

$\begin{array}{cc}{J}_2^0:=& \frac{1}{6}[a_1({\sigma }_{22}-{\sigma }_{33}{)}^2+a_2({\sigma }_{33}-{\sigma }_{11}{)}^2+a_3({\sigma }_{11}-{\sigma }_{22}{)}^2]+a_4{\sigma }_{23}^2+a_5{\sigma }_{31}^2+a_6{\sigma }_{12}^2\\ J_3^0:=& \frac{1}{27}[(b_1+b_2){\sigma }_{11}^3+(b_3+b_4){\sigma }_{22}^3+\{2(b_1+b_4)-(b_2+b_3)\}{\sigma }_{33}^3]\\ & -\frac{1}{9}[(b_1{\sigma }_{22}+b_2{\sigma }_{33}){\sigma }_{11}^2+(b_3{\sigma }_{33}+b_4{\sigma }_{11}){\sigma }_{22}^2+\{(b_1-b_2+b_4){\sigma }_{11}+(b_1-b_3+b_4){\sigma }_{22}\}{\sigma }_{33}^2]\\ & +\frac{2}{9}(b_1+b_4){\sigma }_{11}{\sigma }_{22}{\sigma }_{33}+2b_{11}{\sigma }_{12}{\sigma }_{23}{\sigma }_{31}\\ & -\frac{1}{3}[\{2b_9{\sigma }_{22}-b_8{\sigma }_{33}-(2b_9-b_8){\sigma }_{11}\}{\sigma }_{31}^2+\{2b_{10}{\sigma }_{33}-b_5{\sigma }_{22}-(2b_{10}-b_5){\sigma }_{11}\}{\sigma }_{12}^2\\ & \{(b_6+b_7){\sigma }_{11}-b_6{\sigma }_{22}-b_7{\sigma }_{33}\}{\sigma }_{23}^2]\end{array}$

Критерий текучести Казаку — Барлата для плоского напряженного состояния

Для тонких металлических пластин напряжения могут быть рассмотрены как в случае плоского напряженного состояния. В этом случае критерий текучести Казаку-Барлата сокращается до своей двумерной версии:

$\begin{array}{cc}{J}_2^0=& \frac{1}{6}[(a_2+a_3){\sigma }_{11}^2+(a_1+a_3){\sigma }_{22}^2-2a_3{\sigma }_1{\sigma }_2]+a_6{\sigma }_{12}^2\\ J_3^0=& \frac{1}{27}[(b_1+b_2){\sigma }_{11}^3+(b_3+b_4){\sigma }_{22}^3]-\frac{1}{9}[b_1{\sigma }_{11}+b_4{\sigma }_{22}]{\sigma }_{11}{\sigma }_{22}+\frac{1}{3}[b_5{\sigma }_{22}+(2b_{10}-b_5){\sigma }_{11}]{\sigma }_{12}^2\end{array}$

Для тонких пластин из металла и сплавов параметры критерия текучести Казаку — Барлата могут быть найдены в соответствующих таблицах

Таблица 1. Параметры критерия текучести Казаку-Барлата для металлов и сплавов

Материал	$a_1$	$a_2$	$a_3$	$a_6$	$b_1$	$b_2$	$b_3$	$b_4$	$b_5$	$b_{10}$	$\alpha$
6016-T4 сплав алюминия	0.815	0.815	0.334	0.42	0.04	-1.205	-0.958	0.306	0.153	-0.02	1.4
2090-T3 сплав алюминия	1.05	0.823	0.586	0.96	1.44	0.061	-1.302	-0.281	-0.375	0.445	1.285