Пространственные двумерные решения

Автор: Сенашов С.И., Савостьянова И.Л.

Журнал: Сибирский аэрокосмический журнал @vestnik-sibsau

Рубрика: Информатика, вычислительная техника и управление

Статья в выпуске: 3 т.22, 2021 года.

Бесплатный доступ

В работе рассматриваются стационарные пространственные уравнения идеальной пластичности с условием текучести Мизеса. Материал предполагается несжимаемым. Подробно изучен случай, когда все три компоненты вектора скорости и гидростатическое давление зависят только от двух координат x, y. Для этого случая введено новое название - пространственная двумерная система уравнений, чтобы отличить ее от общепринятых двумерных систем уравнений, когда от нуля отличны только две компоненты вектора скорости и гидростатическое давление. Доказано, что система допускает, в смысле С. Ли, алгебру Ли размерности 10. Показано, что пространственное двумерное деформированное состояние - это есть суперпозиция плоского напряженного состояния и пластического кручения вокруг оси z. Построены два инвариантных решения уравнений, описывающих пространственное двумерное деформированное состояние. Первое решения можно использовать для описания пластических течений между двумя жесткими плитами, которые сближаются с разными скоростями. Второе решение служит для описания напряженно-деформированного состояния материала внутри плоского канала, образованного сходящимися плитами.

Еще

Пространственные решения уравнений идеальной пластичности, точечные симметрии, инвариантные решения

Короткий адрес: https://sciup.org/148323911

IDR: 148323911   |   DOI: 10.31772/2712-8970-2021-22-3-452-456

Текст научной статьи Пространственные двумерные решения

В название статьи вынесено понятие «пространственное двумерное решение». Механикам известно плоское деформированное состояние – это случай, когда в двумерной декартовой системе координат две компоненты вектора скорости деформации и гидростатическое давление зависят от x , y . Известно также плоское напряженное состояние – это когда компоненты тензора напряжений с z , т xz , т yz равны нулю, а компоненты с x , с y , т xy не зависят от z .

В нашем случае все компоненты тензора напряжений не зависят от z , именно такой случай мы и назвали пространственным двумерным состоянием.

Система пространственных уравнений пластичности в декартовой системе координат x 1 = x , x 2 = y , x3 = z в стационарном случае имеет вид

d js = д iP, sij = X(d iuj 1 d A-)/2, d iuj = 0, sysy = 2 kL i, j =1,2,3                (1)

Здесь с j , s ij - компоненты тензора и девиатора тензора напряжений; и 1 = u , и 2 = v , и 3 = w -компоненты вектора скорости деформаций; X - неотрицательная функция; ks - постоянная пластичности; p – гидростатическое давление, по повторяющимся индексам проводится суммирование.

Исключая из системы уравнений (1) компоненты девиатора тензора напряжений, получаем следующую нелинейную систему уравнений

2 k 22 k 22

d i P = "-J5 jj u i -      e ij e mn d mj u n , e ij e ij = A , d i u i = 0.                        (2)

Известно, что система уравнений (2) имеет эллиптический тип. Опишем известные решения этой системы.

Решения этой системы построены Р. Хиллом в 1948 г. [1], В. Прагером в 1954 г. [2], Д. Д. Ивлевым в 1960 г. [3; 4], М. А. Задояном в 1964 г. [5–8], а также авторами этой статьи [9–13]. Отметим также ряд точных решений, построенных Б. Д. Анниным [14], для уравнений пластичности в пространственном случае с условием текучести Треска.

Симметрии системы (2)

Группа точечных симметрий системы уравнений (2) порождается следующими операторами

X, = d ., Y = 6„ ., N = x dx ., M = u au ., i = 1,2,3. xi             ui                 xi                 ui

T 1 = x 2 d и 3 - x3 d и 2 , T 2 = x 3 d u 1 - Х1 д и 3 , T 3 = x1 d и 2 - x 2 d u 1 ,

Z 1 = x2 d x 3 - x 3 d x 2 + и2 д и 3 - и 3 d и 2 , Z 2 = x 3 d x 1 - x 1 d x 3 + и 3 d и 1 - и1 д и 3 , Z 3 = x 1 d x 2 - x 2 d x 1 + и1 д и 2 - и2 д и 1 , S = d p .

Операторам X i , Y i , N , M соответствуют следующие непрерывные преобразования x ' i = x i + a i , и ' i = и 4 + b i , x ' i = x i exp a , и ' i = и^ exp b , i = 1,2,3.

Это переносы по координатам и компонентам вектора скоростей деформации, а также растяжения.

Преобразования Ti показывают, что система уравнений (2) не меняется при жестких перемещениях

u 2 = u 2

- C 1 x 3,

u 3 — u 3

+ C 1 x 2 ,

u 1 u 1

+ C 2 x 3

u 3 — u 3

- с 2 x 1 ,

u 2 u 2

+ C 3 x 1 ,

u 1 — u 1

- C 3 x - 2

Группы, порождаемые операторами Zi, – это вращения вокруг координатных осей x '2 — x2 cos ф1 + x3 sin ф1, x '3 — -x2 sin ф1 + x3 cos ф1, x '3 — x3 cos ф2 + x1 sin ф2, x '1 = -x3 sin ф2 + x1 cos ф2, x '1 = x1 cos ф3 + x2 sin ф3, x '2 = -x1 sin ф3 + x2 cos ф3.

Последнее преобразование описывает инвариантность гидростатического давления относительно сдвигов

p ‘ — p + d.

Во всех этих формулах a i , b i , c i , ф i , d - групповые параметры. Обычно предполагается, что они непрерывно меняются в окрестности нуля.

Проведенные исследования показали, что все построенные решения Р. Хиллом, В. Прагером, Д. Д. Ивлевым и М. А. Задояном суть инвариантные, решения относительно некоторых одномерных подгрупп точечных преобразований, порождаемых операторами (3). Инвариантность здесь означает, что решения не меняются при некоторых преобразованиях, порождаемых группой симметрий (3). Так, решение Р. Хилла инвариантно относительно подалгебры, порожденной оператором 2 C 0 S + X 1 + a Y 1 + в T 1 , решение Д. Д. Ивлева инвариантно относительно подалгебры 2 C 0 S + X 1 + a Y , решение Прагера инвариантно относительно подалгебры aS + X 1 + T 1 + a T 2 , решение Задояна инвариантно относительно этой же подалгебры. О чем говорит этот факт? Он говорит о том, что фактически все эти решения «двумерны», т. е. в подходящей системе координат их можно записать как функции только от двух независимых переменных. Это же можно сказать и о решениях, построенных авторами этой работы. Тогда возникает вопрос: что надо понимать под пространственным решением? Исходя из приведенных здесь решений, ответ получается такой: пространственные решения – это такие решения, которые имеют три компоненты вектора скорости, давление, фактически зависящие от двух переменных в подходящей системе координат. Эти решения являются инвариантными решениями ранга 2. В этом случае задача поиска пространственных решений может быть формально поставлена так: построение новых инвариантных решений ранга 2 для пространственных уравнений идеальной пластичности. Вид таких решений можно без труда перечислить, если перечислить все различные, с точностью до внутренних автоморфизмов, одномерные подалгебры алгебры (3).

Подалгебр, на которых можно построить инвариантные решения ранга 2, существует несколько. Перечислим их.

X 3 + у S , X 1 + Z 1 + у S , a M + N + у S , N + Y + у S ,

Z i + a N + Y + уS, Z i + aN + в M + уS,

X 1 + a Z 1 + M + у S ,

X 1 + Y + a T 1 + у S , X 1 + Z 1 + T 1 + у S , M + N + T 1 + y S , Z i + Y + a T + y S , X i + a X 2 + T 2 + в T 3 + y S , X 1 + Z 1 + Y + a y + y s .

Здесь α , β , γ – произвольные постоянные, разным значениям этих постоянных соответствуют не подобные подалгебры.

В этой работе рассмотрим только решения инвариантные относительно подалгебры X 3 + γ S .

Замечание. Другие решения, построенные на подалгебрах (4), а также вид всех инвариантных решений, которые могут быть построены для системы уравнений (1), можно найти в [11].

Решения. инвариантные относительно этой подалгебры, следует искать в виде u = u(x, y), v = v(x, y), w = w(x, y), p=p(x,y)+γz.                        (5)

Подставляя соотношения (5) в систему уравнений (1) получим

x σ x +∂ y τ xy = 0, x τ xy +∂ y σ y = 0, x τ xz +∂ y τ yz , x u +∂ y v = 0,

( σ x y ) 2 + ( σ x z ) 2 + ( σ y z ) 2 + 6( τ 2 xy 2 xz 2 yz ) = 6 k s 2 ,

(6) σ x - p =λ∂ x u , σ y - p =λ∂ y v , σ z - p = 0, 2 τ xy = λ ( x v + ∂ y u ),

2 τ xz = λ∂ x w , 2 τ yz =λ∂ y w .

Сделаем замену переменных в (6) по следующим формулам

σ x = k ( у/ 3 cos ω + sin ω cos2 ϕ ), σ y = k ( 3 cos ω - sin ω cos2 ϕ ), τ xy = k sin ω sin2 ϕ .          (7)

Здесь k k s , 0 <δ< 1 – некоторая постоянная.

Подставляя эти соотношения в (6) получаем

( - cos ω+ 3sin ω cos2 ϕ ) x ω+ 3sin ω sin2 ϕ∂ y ω - 2sin ω∂ y ϕ = 0,

(cos ω+ 3sin ω cos2 ϕ ) y ω+ 'J 3sin ω sin2 ϕ∂ x ω + 2sin ω∂ x ϕ = 0,

∂ τ+∂τ=γ , u +∂ v = 0, τ 2 2 = k 2 - k 2 = K 2 ,                 (8)

xxz y yz xy       xz yz s

σ x - p =λ∂ x u , σ y - p =λ∂ y v , σ z - p = 0, 2 τ xy = λ ( x v + ∂ y u ), 2 τ xz = λ∂ x w , 2 τ yz = λ∂ y w .

Из (8) видим, что исходная система распалась на две подсистемы: два первых уравнения по сути совпадают с уравнениями, описывающими плоское напряженное состояние

( - cos ω+ 3sin ω cos2 ϕ ) x ω + 3sin ω sin2 ϕ∂ y ω - 2sin ω∂ y ϕ= 0,              (9)

(cosω+3 sinωcos2ϕ)∂yω+ 3sinωsin2ϕ∂xω +2sinω∂xϕ = 0, и уравнениями, напоминающими уравнения, описывающие пластическое кручение стержня, при γ=0 и другом пределе текучести

∂τ +∂τ =γ , τ 2 2 = k 2 - k 2 = K 2 .                        (10)

xxz y yz xz yz s

Решая уравнения (9) и (10) можно найти компоненты тензора напряжений, при этом

σ z = p = 1/2( σ x y ). Единственная проблема – это определение постоянной k .

Для определения компонент вектора скорости получаем следующие уравнения

x u   =    y v   = x v +∂ y v ,   x (    x w   ) +∂ y ( y w  ) .

2 σ x y 2 σ y x     6 τ xy           w x 2 + w y 2          w x 2 + w y 2

Отметим, что последнее уравнение при γ ≠ 0 до сих пор недостаточно исследовано, оно даже не вошло в справочник [15].

В итоге получилось, что решение системы уравнений (6) – это фактически суперпозиция плоского напряженного состояния и пластического кручения вокруг оси z .

Приведем некоторые другие решения уравнений (6). Для этого найдем группу точечных симметрий системы уравнений (6).

Эта группа порождается операторами

X 1 =∂ x , X 2 =∂ y , Y i =∂ u i , N = x x + y y , M = u i ui , i = 1,2,3.

T = x 1 u 2 - x 2 u 1 , Z 3 = x 1 x 2 - x 2 x 1 + u 1 u 2 - u 2 u 1 , S =∂ p .

Ищем решение, инвариантное относительно подалгебры, порождаемой операторами X 1 T + Y 1 Y 3 S .

Решение следует искать в виде

u = x xy + U ( y ), v =-α /2 x 2 + V ( y ), w x + W ( y ), p x + P ( y ).

Из уравнения несжимаемости получаем v=-α/2(x+y)-y+C1.

Подставляем (12) в (6) и получаем систему обыкновенных дифференциальных уравнений d (λU')=γ, d (λW')=0,6ks2λ-2=2(1+αy)2+6((U')2+β2+(W')2).          (13)

dy          dy

Из (13) имеем

λ U ' y + C 2, λ W ' = C 3.

Здесь C 1, C 2, C 3 – произвольные постоянные.

Получаем

U '/ W ' = ( γ y + C 2 )/ C 3 ,

6 k s 2 λ- 2 = 2(1 y ) 2 + 6( β 2 + (1 + ( C 4 y ) 2 )( W ') 2 ), ν=γ / C 3, C 4 = C 2/ C 3.

Систему уравнений (13) удается свести к квадратурам, которые выражаются через эллипти- ческие интегралы.

Построенное решение можно использовать для описания пластического течения слоя, сжи- маемого жесткими плитами ортогональными оси oz.

Запишем систему (6) в цилиндрической системе координат r , θ , z .

∂rσr+r-1∂θσrθ+(σr-σθ)/r=0,∂rσrθ+r-1∂θσθ+2σrθ/r=0,

r σ rz + r - 1 θσθ z + σ rz / r = 0, σ r - p =λ∂ r u , σθ- p ( u / r +∂θ v / r ), σ z - p = 0, 2 σ r θ = λ ( θ u / r + r r ( v / r )), 2 σ z θ r - 1 θ w , 2 σ zr = λ∂ r w ,

( σ r θ ) 2 + ( σ r z ) 2 + ( σθ z ) 2 + 6( σ r 2 θ 2 θ z r 2 z ) = 6 k s 2 , r u + u / r +∂θ v / r = 0.

Ищем инвариантное решение на подалгебре Z 3 + M , оно имеет вид u = U ( θ ), v = V ( θ ), w = W ( θ ), p = P ( θ ) ln r .                     (15)

Подставляя (15) в (14) получаем систему обыкновенных дифференциальных уравнений d θσ r + ( σ r θ ) , d θσθ + 2 σ r θ = 0, d θσ r θ + σ rz = 0, σ r - p = 0, σθ - p ( U / r + V '/ r ), σ z - p = 0, 2 σ r θ = λ ( U '/ r + r r ( V / r )),

2 σ z θ r - 1 W ', 2 σ zr / r ,                                 (16)

( σ r θ ) 2 + ( σ r z ) 2 + ( σθ z ) 2 + 6( σ r 2 θθ 2 z r 2 z ) = 6 k s 2 ,

U + V '=0.

Здесь штрих означает производную по θ .

Отсюда получаем систему обыкновенных дифференциальных уравнений

( X ( U ' - V ))' = a , X ( U ' - V ) = P , ( X W ')' = 0, 6 k ^ X"1 = ( U '- V ) 2 + ( W ') 2 + 1.

Система (17) решается полностью аналогично системе (13).

Найденное решение можно использовать для описания пластического течения в сходящемся плоском канале с жесткими и шероховатыми стенками.

Другие решения системы уравнений (1) можно найти в [11].

Заключение

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

Список литературы Пространственные двумерные решения

  • Хилл Р. Математическая теория пластичности. М. : ГИТТЛ, 1956. 408 с.
  • Прагер В. Трехмерное пластическое течение при однородном напряженном состоянии. Механика // Сб. переводов и обзоров иностр. лит-ры. 1958. № 3. С. 23–27.
  • Предельное состояние деформируемых тел и горных пород / Д. Д. Ивлев, Л. А. Максимова, Р. И. Непершин и др. М. : Физматлит, 2008. 832 с.
  • Ивлев Д. Д. Теория идеальной пластичности. М. : Наука, 1966. 232 с.
  • Ольшак В., Мруз З., Пежина П. Современное состояние теории пластичности. М. : Мир, 1964. 243 с.
  • Задоян М. А. Частное решение уравнений идеальной пластичности // Докл. АН СССР СССР. 1964. Т. 156, № 1. С. 38–39.
  • Задоян М. А. Частное решение уравнений идеальной пластичности в цилиндрических координатах // Докл. АН СССР СССР. 1964. Т. 157. № 1. С. 73–75.
  • Задоян М. А. Пространственные задачи теории пластичности. М. : Наука, 1992. 382 с.
  • Сенашов С. И., Савостьянова И. Л. Новые трехмерные пластические течения, соответствующие однородному напряженному состоянию // Сиб. журн. индуст. матем. 2019. Т. 22, № 3. С. 114–117.
  • Сенашов С. И. Пластические течения среды Мизеса со спирально-винтовой симметрией // Прикладная матем. и механика. 2004. Т. 68, № 1. С. 150–154.
  • Аннин Б. Д., Бытев В. О., Сенашов С. И. Групповые свойства уравнений упругости и пластичности. Новосибирск : Наука, 1983. 140 с.
  • Сенашов С. И. Решение уравнений пластичности в случае спирально-винтовой симметрии // Докл. АН СССР, 1991. Т. 317, № 1. С. 57–59.
  • Сенашов С. И. Решение уравнений пластичности в случае спирально-винтовой симметрии // Известия РАН. Механика твердого тела. 1991. № 5. С. 167–171.
  • Аннин Б. Д. Новые точные решения пространственных уравнений пластичности Треска // Доклады Академии наук. 2007. Т. 415. № 4. С. 482–485.
  • Polyanin A. D., Zaitsev V. F. Handbook of nonlinear partial differential equations. 2nd edition, 2012, Taylor&Francis Group. 1875 p.
Еще
Статья научная