Hyperelastic material model
неогуковское твердое тело [1] [2] — это модель гиперупругого материала , подобная закону Гука , которую можно использовать для прогнозирования нелинейного поведения напряжения и деформации материалов, подвергающихся большим деформациям . Модель была предложена Рональдом Ривлином в 1948 году с использованием инвариантов, хотя Муни уже описал версию в растянутой форме в 1940 году, а Уолл отметил эквивалентность сдвига модели Гука в 1942 году.
В отличие от линейных упругих материалов кривая растяжения неогуковского материала не является линейной . Вместо этого взаимосвязь между приложенным напряжением и деформацией изначально линейна, но в определенной точке кривая напряжения-деформации выходит на плато. Модель нео-Гука не учитывает диссипативное выделение энергии в виде тепла при деформации материала, и предполагается идеальная эластичность на всех стадиях деформации. Помимо использования для моделирования физических материалов, стабильность и сильно нелинейное поведение при сжатии сделали неогуковские материалы популярным выбором для вымышленных подходов к средам, таких как третий метод контакта со средой .
Модель нео-Гука основана на статистической термодинамике сшитых полимерных цепей и пригодна для пластмасс и резиноподобных веществ. Сшитые полимеры будут действовать по принципу неогука, поскольку первоначально полимерные цепи могут перемещаться относительно друг друга при приложении напряжения. Однако в определенный момент полимерные цепи будут растянуты до максимальной точки, которую позволят ковалентные поперечные связи, и это приведет к резкому увеличению модуля упругости материала. Модель материала Нео-Гука не предсказывает такого увеличения модуля при больших деформациях и обычно точна только для деформаций менее 20%. [3] Модель также неадекватна для двухосных напряженных состояний и была заменена моделью Муни-Ривлина .
Функция плотности энергии деформации несжимаемого вид неогуковского материала в трехмерном описании имеет
W = C 1 ( I 1 − 3 ) {\displaystyle W=C_{1}(I_{1}-3)} где C 1 {\displaystyle C_{1}} является материальной константой, а I 1 {\displaystyle I_{1}} — первый инвариант ( след ) правого тензора деформации Коши-Грина , т. е.
I 1 = λ 1 2 + λ 2 2 + λ 3 2 {\displaystyle I_{1}=\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}} где λ i {\displaystyle \lambda _{i}} являются основными участками . [2]
Для сжимаемого неогуковского материала функция плотности энергии деформации определяется выражением
W = C 1 ( I 1 − 3 − 2 ln J ) + D 1 ( J − 1 ) 2 ; J = det ( F ) = λ 1 λ 2 λ 3 {\displaystyle W=C_{1}~(I_{1}-3-2\ln J)+D_{1}~(J-1)^{2}~;~~J=\det({\boldsymbol {F}})=\lambda _{1}\lambda _{2}\lambda _{3}} где D 1 {\displaystyle D_{1}} является материальной константой и F {\displaystyle {\boldsymbol {F}}} – градиент деформации . Можно показать, что в 2D функция плотности энергии деформации равна
W = C 1 ( I 1 − 2 − 2 ln J ) + D 1 ( J − 1 ) 2 {\displaystyle W=C_{1}~(I_{1}-2-2\ln J)+D_{1}~(J-1)^{2}} Существует несколько альтернативных формулировок для сжимаемых неогуковских материалов, например
W = C 1 ( I ¯ 1 − 3 ) + ( C 1 6 + D 1 4 ) ( J 2 + 1 J 2 − 2 ) {\displaystyle W=C_{1}~({\bar {I}}_{1}-3)+\left({\frac {C_{1}}{6}}+{\frac {D_{1}}{4}}\right)\!\left(J^{2}+{\frac {1}{J^{2}}}-2\right)} где I ¯ 1 = J − 2 / 3 I 1 {\displaystyle {\bar {I}}_{1}=J^{-2/3}I_{1}} является инвариантом изохорной части первым C ¯ = ( det C ) − 1 / 3 C = J − 2 / 3 C {\displaystyle {\bar {\boldsymbol {C}}}=(\det {\boldsymbol {C}})^{-1/3}{\boldsymbol {C}}=J^{-2/3}{\boldsymbol {C}}} правого тензора деформации Коши–Грина .
Для обеспечения соответствия линейной эластичности,
C 1 = μ 2 ; D 1 = λ L 2 {\displaystyle C_{1}={\frac {\mu }{2}}~;~~D_{1}={\frac {{\lambda }_{L}}{2}}} где λ L {\displaystyle {\lambda }_{L}} – первый параметр Ламе и μ {\displaystyle \mu } – модуль сдвига или второй параметр Ламе . [4] Альтернативные определения C 1 {\displaystyle C_{1}} и D 1 {\displaystyle D_{1}} иногда используются, особенно в коммерческом программном обеспечении для анализа методом конечных элементов, таком как Abaqus . [5]
Для сжимаемого Огдена неогуковского материала напряжение Коши определяется выражением
σ = J − 1 P F T = J − 1 ∂ W ∂ F F T = J − 1 ( 2 C 1 ( F − F − T ) + 2 D 1 ( J − 1 ) J F − T ) F T {\displaystyle {\boldsymbol {\sigma }}=J^{-1}{\boldsymbol {P}}{\boldsymbol {F}}^{T}=J^{-1}{\frac {\partial W}{\partial {\boldsymbol {F}}}}{\boldsymbol {F}}^{T}=J^{-1}\left(2C_{1}({\boldsymbol {F}}-{\boldsymbol {F}}^{-T})+2D_{1}(J-1)J{\boldsymbol {F}}^{-T}\right){\boldsymbol {F}}^{T}} где P {\displaystyle {\boldsymbol {P}}} – первое напряжение Пиолы-Кирхгофа. Упрощая правую часть, приходим к
σ = 2 C 1 J − 1 ( F F T − I ) + 2 D 1 ( J − 1 ) I = 2 C 1 J − 1 ( B − I ) + 2 D 1 ( J − 1 ) I {\displaystyle {\boldsymbol {\sigma }}=2C_{1}J^{-1}\left({\boldsymbol {F}}{\boldsymbol {F}}^{T}-{\boldsymbol {I}}\right)+2D_{1}(J-1){\boldsymbol {I}}=2C_{1}J^{-1}\left({\boldsymbol {B}}-{\boldsymbol {I}}\right)+2D_{1}(J-1){\boldsymbol {I}}} которое для бесконечно малых деформаций равно
≈ 4 C 1 ε + 2 D 1 tr ( ε ) I {\displaystyle \approx 4C_{1}{\boldsymbol {\varepsilon }}+2D_{1}\operatorname {tr} ({\boldsymbol {\varepsilon }}){\boldsymbol {I}}} Сравнение с законом Гука показывает, что C 1 = μ 2 {\displaystyle C_{1}={\tfrac {\mu }{2}}} и D 1 = λ L 2 {\displaystyle D_{1}={\tfrac {\lambda _{L}}{2}}} .
Для сжимаемого Ривлина неогуковского материала напряжение Коши определяется выражением
J σ = − p I + 2 C 1 dev ( B ¯ ) = − p I + 2 C 1 J 2 / 3 dev ( B ) {\displaystyle J~{\boldsymbol {\sigma }}=-p~{\boldsymbol {I}}+2C_{1}\operatorname {dev} ({\bar {\boldsymbol {B}}})=-p~{\boldsymbol {I}}+{\frac {2C_{1}}{J^{2/3}}}\operatorname {dev} ({\boldsymbol {B}})} где B {\displaystyle {\boldsymbol {B}}} – левый тензор деформации Коши–Грина,
p := − 2 D 1 J ( J − 1 ) ; dev ( B ¯ ) = B ¯ − 1 3 I ¯ 1 I ; B ¯ = J − 2 / 3 B . {\displaystyle p:=-2D_{1}~J(J-1)~;~\operatorname {dev} ({\bar {\boldsymbol {B}}})={\bar {\boldsymbol {B}}}-{\tfrac {1}{3}}{\bar {I}}_{1}{\boldsymbol {I}}~;~~{\bar {\boldsymbol {B}}}=J^{-2/3}{\boldsymbol {B}}~.} Для бесконечно малых деформаций ( ε {\displaystyle {\boldsymbol {\varepsilon }}} )
J ≈ 1 + tr ( ε ) ; B ≈ I + 2 ε {\displaystyle J\approx 1+\operatorname {tr} ({\boldsymbol {\varepsilon }})~;~~{\boldsymbol {B}}\approx {\boldsymbol {I}}+2{\boldsymbol {\varepsilon }}} а напряжение Коши можно выразить как
σ ≈ 4 C 1 ( ε − 1 3 tr ( ε ) I ) + 2 D 1 tr ( ε ) I {\displaystyle {\boldsymbol {\sigma }}\approx 4C_{1}\left({\boldsymbol {\varepsilon }}-{\tfrac {1}{3}}\operatorname {tr} ({\boldsymbol {\varepsilon }}){\boldsymbol {I}}\right)+2D_{1}\operatorname {tr} ({\boldsymbol {\varepsilon }}){\boldsymbol {I}}} Сравнение с законом Гука показывает, что μ = 2 C 1 {\displaystyle \mu =2C_{1}} и κ = 2 D 1 {\displaystyle \kappa =2D_{1}} .
показывать Доказательство: The Cauchy stress in a compressible hyperelastic material is given by
σ = 2 J [ 1 J 2 / 3 ( ∂ W ∂ I ¯ 1 + I ¯ 1 ∂ W ∂ I ¯ 2 ) B − 1 J 4 / 3 ∂ W ∂ I ¯ 2 B ⋅ B ] + [ ∂ W ∂ J − 2 3 J ( I ¯ 1 ∂ W ∂ I ¯ 1 + 2 I ¯ 2 ∂ W ∂ I ¯ 2 ) ] I {\displaystyle {\boldsymbol {\sigma }}={\cfrac {2}{J}}\left[{\cfrac {1}{J^{2/3}}}\left({\cfrac {\partial {W}}{\partial {\bar {I}}_{1}}}+{\bar {I}}_{1}~{\cfrac {\partial {W}}{\partial {\bar {I}}_{2}}}\right){\boldsymbol {B}}-{\cfrac {1}{J^{4/3}}}~{\cfrac {\partial {W}}{\partial {\bar {I}}_{2}}}~{\boldsymbol {B}}\cdot {\boldsymbol {B}}\right]+\left[{\cfrac {\partial {W}}{\partial J}}-{\cfrac {2}{3J}}\left({\bar {I}}_{1}~{\cfrac {\partial {W}}{\partial {\bar {I}}_{1}}}+2~{\bar {I}}_{2}~{\cfrac {\partial {W}}{\partial {\bar {I}}_{2}}}\right)\right]~{\boldsymbol {I}}} For a compressible Rivlin neo-Hookean material,
∂ W ∂ I ¯ 1 = C 1 ; ∂ W ∂ I ¯ 2 = 0 ; ∂ W ∂ J = 2 D 1 ( J − 1 ) {\displaystyle {\cfrac {\partial {W}}{\partial {\bar {I}}_{1}}}=C_{1}~;~~{\cfrac {\partial {W}}{\partial {\bar {I}}_{2}}}=0~;~~{\cfrac {\partial {W}}{\partial J}}=2D_{1}(J-1)} while, for a compressible Ogden neo-Hookean material,
∂ W ∂ I ¯ 1 = C 1 ; ∂ W ∂ I ¯ 2 = 0 ; ∂ W ∂ J = 2 D 1 ( J − 1 ) − 2 C 1 J {\displaystyle {\cfrac {\partial {W}}{\partial {\bar {I}}_{1}}}=C_{1}~;~~{\cfrac {\partial {W}}{\partial {\bar {I}}_{2}}}=0~;~~{\cfrac {\partial {W}}{\partial J}}=2D_{1}(J-1)-{\cfrac {2C_{1}}{J}}} Therefore, the Cauchy stress in a compressible Rivlin neo-Hookean material is given by
σ = 2 J [ 1 J 2 / 3 C 1 B ] + [ 2 D 1 ( J − 1 ) − 2 3 J C 1 I ¯ 1 ] I {\displaystyle {\boldsymbol {\sigma }}={\cfrac {2}{J}}\left[{\cfrac {1}{J^{2/3}}}~C_{1}~{\boldsymbol {B}}\right]+\left[2D_{1}(J-1)-{\cfrac {2}{3J}}~C_{1}{\bar {I}}_{1}\right]{\boldsymbol {I}}} while that for the corresponding Ogden material is
σ = 2 J [ 1 J 2 / 3 C 1 B ] + [ 2 D 1 ( J − 1 ) − 2 C 1 J − 2 3 J C 1 I ¯ 1 ] I {\displaystyle {\boldsymbol {\sigma }}={\cfrac {2}{J}}\left[{\cfrac {1}{J^{2/3}}}~C_{1}~{\boldsymbol {B}}\right]+\left[2D_{1}(J-1)-{\cfrac {2C_{1}}{J}}-{\cfrac {2}{3J}}~C_{1}{\bar {I}}_{1}\right]{\boldsymbol {I}}} If the isochoric part of the left Cauchy-Green deformation tensor is defined as B ¯ = J − 2 / 3 B {\displaystyle {\bar {\boldsymbol {B}}}=J^{-2/3}{\boldsymbol {B}}} , then we can write the Rivlin neo-Heooken stress as
σ = 2 C 1 J [ B ¯ − 1 3 I ¯ 1 I ] + 2 D 1 ( J − 1 ) I = 2 C 1 J dev ( B ¯ ) + 2 D 1 ( J − 1 ) I {\displaystyle {\boldsymbol {\sigma }}={\cfrac {2C_{1}}{J}}\left[{\bar {\boldsymbol {B}}}-{\tfrac {1}{3}}{\bar {I}}_{1}{\boldsymbol {I}}\right]+2D_{1}(J-1){\boldsymbol {I}}={\cfrac {2C_{1}}{J}}\operatorname {dev} ({\bar {\boldsymbol {B}}})+2D_{1}(J-1){\boldsymbol {I}}} and the Ogden neo-Hookean stress as
σ = 2 C 1 J [ B ¯ − 1 3 I ¯ 1 I − I ] + 2 D 1 ( J − 1 ) I = 2 C 1 J [ dev ( B ¯ ) − I ] + 2 D 1 ( J − 1 ) I {\displaystyle {\boldsymbol {\sigma }}={\cfrac {2C_{1}}{J}}\left[{\bar {\boldsymbol {B}}}-{\tfrac {1}{3}}{\bar {I}}_{1}{\boldsymbol {I}}-{\boldsymbol {I}}\right]+2D_{1}(J-1){\boldsymbol {I}}={\cfrac {2C_{1}}{J}}\left[\operatorname {dev} ({\bar {\boldsymbol {B}}})-{\boldsymbol {I}}\right]+2D_{1}(J-1){\boldsymbol {I}}} The quantities
p := − 2 D 1 J ( J − 1 ) ; p ∗ = − 2 D 1 J ( J − 1 ) + 2 C 1 {\displaystyle p:=-2D_{1}~J(J-1)~;~~p^{*}=-2D_{1}~J(J-1)+2C_{1}} have the form of pressures and are usually treated as such. The Rivlin neo-Hookean stress can then be expressed in the form
τ = J σ = − p I + 2 C 1 dev ( B ¯ ) {\displaystyle {\boldsymbol {\tau }}=J~{\boldsymbol {\sigma }}=-p{\boldsymbol {I}}+2C_{1}\operatorname {dev} ({\bar {\boldsymbol {B}}})} while the Ogden neo-Hookean stress has the form
τ = − p ∗ I + 2 C 1 dev ( B ¯ ) {\displaystyle {\boldsymbol {\tau }}=-p^{*}{\boldsymbol {I}}+2C_{1}\operatorname {dev} ({\bar {\boldsymbol {B}}})}
Для несжимаемого неогуковского материала с J = 1 {\displaystyle J=1}
σ = − p I + 2 C 1 B {\displaystyle {\boldsymbol {\sigma }}=-p~{\boldsymbol {I}}+2C_{1}{\boldsymbol {B}}} где p {\displaystyle p} это неопределенное давление.
Напряжение Коши с точки зрения главных растяжений [ редактировать ] Для сжимаемого неогуковского гиперупругого материала главные компоненты напряжения Коши определяются выражением
σ i = 2 C 1 J − 5 / 3 [ λ i 2 − I 1 3 ] + 2 D 1 ( J − 1 ) ; i = 1 , 2 , 3 {\displaystyle \sigma _{i}=2C_{1}J^{-5/3}\left[\lambda _{i}^{2}-{\cfrac {I_{1}}{3}}\right]+2D_{1}(J-1)~;~~i=1,2,3} Следовательно, разности главных напряжений равны
σ 11 − σ 33 = 2 C 1 J 5 / 3 ( λ 1 2 − λ 3 2 ) ; σ 22 − σ 33 = 2 C 1 J 5 / 3 ( λ 2 2 − λ 3 2 ) {\displaystyle \sigma _{11}-\sigma _{33}={\cfrac {2C_{1}}{J^{5/3}}}(\lambda _{1}^{2}-\lambda _{3}^{2})~;~~\sigma _{22}-\sigma _{33}={\cfrac {2C_{1}}{J^{5/3}}}(\lambda _{2}^{2}-\lambda _{3}^{2})} показывать Доказательство: For a compressible hyperelastic material , the principal components of the Cauchy stress are given by
σ i = λ i λ 1 λ 2 λ 3 ∂ W ∂ λ i ; i = 1 , 2 , 3 {\displaystyle \sigma _{i}={\cfrac {\lambda _{i}}{\lambda _{1}\lambda _{2}\lambda _{3}}}~{\frac {\partial W}{\partial \lambda _{i}}}~;~~i=1,2,3} The strain energy density function for a compressible neo Hookean material is
W = C 1 ( I ¯ 1 − 3 ) + D 1 ( J − 1 ) 2 = C 1 [ J − 2 / 3 ( λ 1 2 + λ 2 2 + λ 3 2 ) − 3 ] + D 1 ( J − 1 ) 2 {\displaystyle W=C_{1}({\bar {I}}_{1}-3)+D_{1}(J-1)^{2}=C_{1}\left[J^{-2/3}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2})-3\right]+D_{1}(J-1)^{2}} Therefore,
λ i ∂ W ∂ λ i = C 1 [ − 2 3 J − 5 / 3 λ i ∂ J ∂ λ i ( λ 1 2 + λ 2 2 + λ 3 2 ) + 2 J − 2 / 3 λ i 2 ] + 2 D 1 ( J − 1 ) λ i ∂ J ∂ λ i {\displaystyle \lambda _{i}{\frac {\partial W}{\partial \lambda _{i}}}=C_{1}\left[-{\frac {2}{3}}J^{-5/3}\lambda _{i}{\frac {\partial J}{\partial \lambda _{i}}}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2})+2J^{-2/3}\lambda _{i}^{2}\right]+2D_{1}(J-1)\lambda _{i}{\frac {\partial J}{\partial \lambda _{i}}}} Since J = λ 1 λ 2 λ 3 {\displaystyle J=\lambda _{1}\lambda _{2}\lambda _{3}} we have
λ i ∂ J ∂ λ i = λ 1 λ 2 λ 3 = J {\displaystyle \lambda _{i}{\frac {\partial J}{\partial \lambda _{i}}}=\lambda _{1}\lambda _{2}\lambda _{3}=J} Hence,
λ i ∂ W ∂ λ i = C 1 [ − 2 3 J − 2 / 3 ( λ 1 2 + λ 2 2 + λ 3 2 ) + 2 J − 2 / 3 λ i 2 ] + 2 D 1 J ( J − 1 ) = 2 C 1 J − 2 / 3 [ − 1 3 ( λ 1 2 + λ 2 2 + λ 3 2 ) + λ i 2 ] + 2 D 1 J ( J − 1 ) {\displaystyle {\begin{aligned}\lambda _{i}{\frac {\partial W}{\partial \lambda _{i}}}&=C_{1}\left[-{\frac {2}{3}}J^{-2/3}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2})+2J^{-2/3}\lambda _{i}^{2}\right]+2D_{1}J(J-1)\\&=2C_{1}J^{-2/3}\left[-{\frac {1}{3}}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2})+\lambda _{i}^{2}\right]+2D_{1}J(J-1)\end{aligned}}} The principal Cauchy stresses are therefore given by
σ i = 2 C 1 J − 5 / 3 [ λ i 2 − I 1 3 ] + 2 D 1 ( J − 1 ) {\displaystyle \sigma _{i}=2C_{1}J^{-5/3}\left[\lambda _{i}^{2}-{\cfrac {I_{1}}{3}}\right]+2D_{1}(J-1)}
С точки зрения главных растяжений разность напряжений Коши для несжимаемого гиперупругого материала определяется выражением
σ 11 − σ 33 = λ 1 ∂ W ∂ λ 1 − λ 3 ∂ W ∂ λ 3 ; σ 22 − σ 33 = λ 2 ∂ W ∂ λ 2 − λ 3 ∂ W ∂ λ 3 {\displaystyle \sigma _{11}-\sigma _{33}=\lambda _{1}~{\cfrac {\partial {W}}{\partial \lambda _{1}}}-\lambda _{3}~{\cfrac {\partial {W}}{\partial \lambda _{3}}}~;~~\sigma _{22}-\sigma _{33}=\lambda _{2}~{\cfrac {\partial {W}}{\partial \lambda _{2}}}-\lambda _{3}~{\cfrac {\partial {W}}{\partial \lambda _{3}}}} Для несжимаемого неогуковского материала
W = C 1 ( λ 1 2 + λ 2 2 + λ 3 2 − 3 ) ; λ 1 λ 2 λ 3 = 1 {\displaystyle W=C_{1}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}-3)~;~~\lambda _{1}\lambda _{2}\lambda _{3}=1} Поэтому,
∂ W ∂ λ 1 = 2 C 1 λ 1 ; ∂ W ∂ λ 2 = 2 C 1 λ 2 ; ∂ W ∂ λ 3 = 2 C 1 λ 3 {\displaystyle {\cfrac {\partial {W}}{\partial \lambda _{1}}}=2C_{1}\lambda _{1}~;~~{\cfrac {\partial {W}}{\partial \lambda _{2}}}=2C_{1}\lambda _{2}~;~~{\cfrac {\partial {W}}{\partial \lambda _{3}}}=2C_{1}\lambda _{3}} что дает
σ 11 − σ 33 = 2 ( λ 1 2 − λ 3 2 ) C 1 ; σ 22 − σ 33 = 2 ( λ 2 2 − λ 3 2 ) C 1 {\displaystyle \sigma _{11}-\sigma _{33}=2(\lambda _{1}^{2}-\lambda _{3}^{2})C_{1}~;~~\sigma _{22}-\sigma _{33}=2(\lambda _{2}^{2}-\lambda _{3}^{2})C_{1}} Истинное напряжение как функция одноосного растяжения, предсказанное сжимаемым неогуковским материалом для различных значений C 1 , D 1 {\displaystyle C_{1},D_{1}} . Свойства материала соответствуют натуральному каучуку . Для сжимаемого материала, подвергающегося одноосному растяжению, основные растяжения равны
λ 1 = λ ; λ 2 = λ 3 = J λ ; I 1 = λ 2 + 2 J λ {\displaystyle \lambda _{1}=\lambda ~;~~\lambda _{2}=\lambda _{3}={\sqrt {\tfrac {J}{\lambda }}}~;~~I_{1}=\lambda ^{2}+{\tfrac {2J}{\lambda }}} Следовательно, истинные напряжения (Коши) для сжимаемого неогуковского материала определяются выражением
σ 11 = 4 C 1 3 J 5 / 3 ( λ 2 − J λ ) + 2 D 1 ( J − 1 ) σ 22 = σ 33 = 2 C 1 3 J 5 / 3 ( J λ − λ 2 ) + 2 D 1 ( J − 1 ) {\displaystyle {\begin{aligned}\sigma _{11}&={\cfrac {4C_{1}}{3J^{5/3}}}\left(\lambda ^{2}-{\tfrac {J}{\lambda }}\right)+2D_{1}(J-1)\\\sigma _{22}&=\sigma _{33}={\cfrac {2C_{1}}{3J^{5/3}}}\left({\tfrac {J}{\lambda }}-\lambda ^{2}\right)+2D_{1}(J-1)\end{aligned}}} Различия в напряжениях определяются выражением
σ 11 − σ 33 = 2 C 1 J 5 / 3 ( λ 2 − J λ ) ; σ 22 − σ 33 = 0 {\displaystyle \sigma _{11}-\sigma _{33}={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\tfrac {J}{\lambda }}\right)~;~~\sigma _{22}-\sigma _{33}=0} Если материал неограничен, мы имеем σ 22 = σ 33 = 0 {\displaystyle \sigma _{22}=\sigma _{33}=0} . Затем
σ 11 = 2 C 1 J 5 / 3 ( λ 2 − J λ ) {\displaystyle \sigma _{11}={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\tfrac {J}{\lambda }}\right)} Приравнивая два выражения для σ 11 {\displaystyle \sigma _{11}} дает отношение к J {\displaystyle J} как функция λ {\displaystyle \lambda } , то есть,
4 C 1 3 J 5 / 3 ( λ 2 − J λ ) + 2 D 1 ( J − 1 ) = 2 C 1 J 5 / 3 ( λ 2 − J λ ) {\displaystyle {\cfrac {4C_{1}}{3J^{5/3}}}\left(\lambda ^{2}-{\tfrac {J}{\lambda }}\right)+2D_{1}(J-1)={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\tfrac {J}{\lambda }}\right)} или
D 1 J 8 / 3 − D 1 J 5 / 3 + C 1 3 λ J − C 1 λ 2 3 = 0 {\displaystyle D_{1}J^{8/3}-D_{1}J^{5/3}+{\tfrac {C_{1}}{3\lambda }}J-{\tfrac {C_{1}\lambda ^{2}}{3}}=0} Приведенное выше уравнение можно решить численно с использованием Ньютона – Рафсона итеративной процедуры поиска корня .
Сравнение экспериментальных результатов (точки) и прогнозов для закона Гука (1), твердого тела нео-Гука (2) и твердотельных моделей Муни-Ривлина (3) При одноосном растяжении λ 1 = λ {\displaystyle \lambda _{1}=\lambda \,} и λ 2 = λ 3 = 1 / λ {\displaystyle \lambda _{2}=\lambda _{3}=1/{\sqrt {\lambda }}} . Поэтому,
σ 11 − σ 33 = 2 C 1 ( λ 2 − 1 λ ) ; σ 22 − σ 33 = 0 {\displaystyle \sigma _{11}-\sigma _{33}=2C_{1}\left(\lambda ^{2}-{\cfrac {1}{\lambda }}\right)~;~~\sigma _{22}-\sigma _{33}=0} Предполагая отсутствие сцепления по бокам, σ 22 = σ 33 = 0 {\displaystyle \sigma _{22}=\sigma _{33}=0} , поэтому мы можем написать
σ 11 = 2 C 1 ( λ 2 − 1 λ ) = 2 C 1 ( 3 ε 11 + 3 ε 11 2 + ε 11 3 1 + ε 11 ) {\displaystyle \sigma _{11}=2C_{1}\left(\lambda ^{2}-{\cfrac {1}{\lambda }}\right)=2C_{1}\left({\frac {3\varepsilon _{11}+3\varepsilon _{11}^{2}+\varepsilon _{11}^{3}}{1+\varepsilon _{11}}}\right)} где ε 11 = λ − 1 {\displaystyle \varepsilon _{11}=\lambda -1} это инженерное напряжение . Это уравнение часто записывают в альтернативных обозначениях как
T 11 = 2 C 1 ( α 2 − 1 α ) {\displaystyle T_{11}=2C_{1}\left(\alpha ^{2}-{\cfrac {1}{\alpha }}\right)} Приведенное выше уравнение предназначено для истинного напряжения (отношения силы растяжения к деформированному поперечному сечению). Для инженерного напряжения уравнение имеет вид:
σ 11 e n g = 2 C 1 ( λ − 1 λ 2 ) {\displaystyle \sigma _{11}^{\mathrm {eng} }=2C_{1}\left(\lambda -{\cfrac {1}{\lambda ^{2}}}\right)} При небольших деформациях ε ≪ 1 {\displaystyle \varepsilon \ll 1} у нас будет:
σ 11 = 6 C 1 ε = 3 μ ε {\displaystyle \sigma _{11}=6C_{1}\varepsilon =3\mu \varepsilon } Таким образом, эквивалентный модуль Юнга неогуковского твердого тела в одноосном растяжении равен 3 μ {\displaystyle 3\mu } , что соответствует линейной эластичности ( E = 2 μ ( 1 + ν ) {\displaystyle E=2\mu (1+\nu )} с ν = 0.5 {\displaystyle \nu =0.5} несжимаемость).
Истинное напряжение как функция двухосного растяжения, предсказанное сжимаемым неогуковским материалом для различных значений C 1 , D 1 {\displaystyle C_{1},D_{1}} . Свойства материала соответствуют натуральному каучуку . В случае равноосного растяжения
λ 1 = λ 2 = λ ; λ 3 = J λ 2 ; I 1 = 2 λ 2 + J 2 λ 4 {\displaystyle \lambda _{1}=\lambda _{2}=\lambda ~;~~\lambda _{3}={\tfrac {J}{\lambda ^{2}}}~;~~I_{1}=2\lambda ^{2}+{\tfrac {J^{2}}{\lambda ^{4}}}} Поэтому,
σ 11 = 2 C 1 [ λ 2 J 5 / 3 − 1 3 J ( 2 λ 2 + J 2 λ 4 ) ] + 2 D 1 ( J − 1 ) = σ 22 σ 33 = 2 C 1 [ J 1 / 3 λ 4 − 1 3 J ( 2 λ 2 + J 2 λ 4 ) ] + 2 D 1 ( J − 1 ) {\displaystyle {\begin{aligned}\sigma _{11}&=2C_{1}\left[{\cfrac {\lambda ^{2}}{J^{5/3}}}-{\cfrac {1}{3J}}\left(2\lambda ^{2}+{\cfrac {J^{2}}{\lambda ^{4}}}\right)\right]+2D_{1}(J-1)\\&=\sigma _{22}\\\sigma _{33}&=2C_{1}\left[{\cfrac {J^{1/3}}{\lambda ^{4}}}-{\cfrac {1}{3J}}\left(2\lambda ^{2}+{\cfrac {J^{2}}{\lambda ^{4}}}\right)\right]+2D_{1}(J-1)\end{aligned}}} Различия в стрессах
σ 11 − σ 22 = 0 ; σ 11 − σ 33 = 2 C 1 J 5 / 3 ( λ 2 − J 2 λ 4 ) {\displaystyle \sigma _{11}-\sigma _{22}=0~;~~\sigma _{11}-\sigma _{33}={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\cfrac {J^{2}}{\lambda ^{4}}}\right)} Если материал находится в состоянии плоского напряжения, то σ 33 = 0 {\displaystyle \sigma _{33}=0} и у нас есть
σ 11 = σ 22 = 2 C 1 J 5 / 3 ( λ 2 − J 2 λ 4 ) {\displaystyle \sigma _{11}=\sigma _{22}={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\cfrac {J^{2}}{\lambda ^{4}}}\right)} У нас также есть отношения между J {\displaystyle J} и λ {\displaystyle \lambda } :
2 C 1 [ λ 2 J 5 / 3 − 1 3 J ( 2 λ 2 + J 2 λ 4 ) ] + 2 D 1 ( J − 1 ) = 2 C 1 J 5 / 3 ( λ 2 − J 2 λ 4 ) {\displaystyle 2C_{1}\left[{\cfrac {\lambda ^{2}}{J^{5/3}}}-{\cfrac {1}{3J}}\left(2\lambda ^{2}+{\cfrac {J^{2}}{\lambda ^{4}}}\right)\right]+2D_{1}(J-1)={\cfrac {2C_{1}}{J^{5/3}}}\left(\lambda ^{2}-{\cfrac {J^{2}}{\lambda ^{4}}}\right)} или,
( 2 D 1 − C 1 λ 4 ) J 2 + 3 C 1 λ 4 J 4 / 3 − 3 D 1 J − 2 C 1 λ 2 = 0 {\displaystyle \left(2D_{1}-{\cfrac {C_{1}}{\lambda ^{4}}}\right)J^{2}+{\cfrac {3C_{1}}{\lambda ^{4}}}J^{4/3}-3D_{1}J-2C_{1}\lambda ^{2}=0} Это уравнение можно решить для J {\displaystyle J} используя метод Ньютона.
Для несжимаемого материала J = 1 {\displaystyle J=1} а разности главных напряжений Коши принимают вид
σ 11 − σ 22 = 0 ; σ 11 − σ 33 = 2 C 1 ( λ 2 − 1 λ 4 ) {\displaystyle \sigma _{11}-\sigma _{22}=0~;~~\sigma _{11}-\sigma _{33}=2C_{1}\left(\lambda ^{2}-{\cfrac {1}{\lambda ^{4}}}\right)} В условиях плоского напряжения имеем
σ 11 = 2 C 1 ( λ 2 − 1 λ 4 ) {\displaystyle \sigma _{11}=2C_{1}\left(\lambda ^{2}-{\cfrac {1}{\lambda ^{4}}}\right)} В случае чистой дилатации
λ 1 = λ 2 = λ 3 = λ : J = λ 3 ; I 1 = 3 λ 2 {\displaystyle \lambda _{1}=\lambda _{2}=\lambda _{3}=\lambda ~:~~J=\lambda ^{3}~;~~I_{1}=3\lambda ^{2}} Следовательно, главные напряжения Коши для сжимаемого неогуковского материала определяются выражением
σ i = 2 C 1 ( 1 λ 3 − 1 λ ) + 2 D 1 ( λ 3 − 1 ) {\displaystyle \sigma _{i}=2C_{1}\left({\cfrac {1}{\lambda ^{3}}}-{\cfrac {1}{\lambda }}\right)+2D_{1}(\lambda ^{3}-1)} Если материал несжимаем, то λ 3 = 1 {\displaystyle \lambda ^{3}=1} а главные напряжения могут быть произвольными.
На рисунках ниже показано, что для достижения большого трехосного растяжения или сжатия необходимы чрезвычайно высокие напряжения. Аналогичным образом, относительно небольшие состояния трехосного растяжения могут вызвать развитие очень высоких напряжений в резиноподобном материале. Величина напряжения весьма чувствительна к объемному модулю, но не к модулю сдвига.
Истинное напряжение как функция равнотрехосного растяжения, предсказанное сжимаемым неогуковским материалом для различных значений C 1 , D 1 {\displaystyle C_{1},D_{1}} . Свойства материала соответствуют натуральному каучуку . Истинное напряжение как функция J, предсказанное с помощью сжимаемого неогуковского материала для различных значений C 1 , D 1 {\displaystyle C_{1},D_{1}} . Свойства материала соответствуют натуральному каучуку .
Для случая простого сдвига градиент деформации в компонентах относительно эталонного базиса имеет вид [2]
F = [ 1 γ 0 0 1 0 0 0 1 ] {\displaystyle {\boldsymbol {F}}={\begin{bmatrix}1&\gamma &0\\0&1&0\\0&0&1\end{bmatrix}}} где γ {\displaystyle \gamma } – сдвиговая деформация. Следовательно, левый тензор деформации Коши-Грина равен
B = F ⋅ F T = [ 1 + γ 2 γ 0 γ 1 0 0 0 1 ] {\displaystyle {\boldsymbol {B}}={\boldsymbol {F}}\cdot {\boldsymbol {F}}^{T}={\begin{bmatrix}1+\gamma ^{2}&\gamma &0\\\gamma &1&0\\0&0&1\end{bmatrix}}} В этом случае J = det ( F ) = 1 {\displaystyle J=\det({\boldsymbol {F}})=1} . Следовательно, σ = 2 C 1 dev ( B ) {\displaystyle {\boldsymbol {\sigma }}=2C_{1}\operatorname {dev} ({\boldsymbol {B}})} . Сейчас,
dev ( B ) = B − 1 3 tr ( B ) I = B − 1 3 ( 3 + γ 2 ) I = [ 2 3 γ 2 γ 0 γ − 1 3 γ 2 0 0 0 − 1 3 γ 2 ] {\displaystyle \operatorname {dev} ({\boldsymbol {B}})={\boldsymbol {B}}-{\tfrac {1}{3}}\operatorname {tr} ({\boldsymbol {B}}){\boldsymbol {I}}={\boldsymbol {B}}-{\tfrac {1}{3}}(3+\gamma ^{2}){\boldsymbol {I}}={\begin{bmatrix}{\tfrac {2}{3}}\gamma ^{2}&\gamma &0\\\gamma &-{\tfrac {1}{3}}\gamma ^{2}&0\\0&0&-{\tfrac {1}{3}}\gamma ^{2}\end{bmatrix}}} Следовательно, напряжение Коши определяется выражением
σ = [ 4 C 1 3 γ 2 2 C 1 γ 0 2 C 1 γ − 2 C 1 3 γ 2 0 0 0 − 2 C 1 3 γ 2 ] {\displaystyle {\boldsymbol {\sigma }}={\begin{bmatrix}{\tfrac {4C_{1}}{3}}\gamma ^{2}&2C_{1}\gamma &0\\2C_{1}\gamma &-{\tfrac {2C_{1}}{3}}\gamma ^{2}&0\\0&0&-{\tfrac {2C_{1}}{3}}\gamma ^{2}\end{bmatrix}}} Используя соотношение для напряжения Коши для несжимаемого неогуковского материала, получаем
σ = − p I + 2 C 1 B = [ 2 C 1 ( 1 + γ 2 ) − p 2 C 1 γ 0 2 C 1 γ 2 C 1 − p 0 0 0 2 C 1 − p ] {\displaystyle {\boldsymbol {\sigma }}=-p~{\boldsymbol {I}}+2C_{1}{\boldsymbol {B}}={\begin{bmatrix}2C_{1}(1+\gamma ^{2})-p&2C_{1}\gamma &0\\2C_{1}\gamma &2C_{1}-p&0\\0&0&2C_{1}-p\end{bmatrix}}} Таким образом, неогуковское тело демонстрирует линейную зависимость касательных напряжений от сдвиговой деформации и квадратичную зависимость разности нормальных напряжений от сдвиговой деформации. Выражения для напряжения Коши для сжимаемого и несжимаемого неогуковского материала при простом сдвиге представляют одну и ту же величину и позволяют определить неизвестное давление. p {\displaystyle p} .
^ Трелоар, LRG (1943). «Эластичность сети длинноцепочечных молекул — II» . Труды Фарадеевского общества . 39 : 241–246. ^ Перейти обратно: а б с Огден, RW (26 апреля 2013 г.). Нелинейные упругие деформации . Курьерская корпорация. ISBN 978-0-486-31871-4 . ^ Гент, АН, изд., 2001, Инженерное дело с резиной , Карл Хансер Верлаг, Мюнхен. ^ Пенс, Ти Джей, и Гоу, К. (2015). О сжимаемых вариантах несжимаемого неогуковского материала. Математика и механика твердого тела , 20(2), 157–182. [1] ^ Abaqus (Версия 6.8) Теоретическое руководство