Nonstationary 1D Dynamics Problems for Heteromodular Elasticity with Piecewise-Linear Approximation of Boundary Conditions

Abstract


The paper provides the investigation of a heteromodular elastic medium under dynamic loading. The heteromodularity (when the stress - strain relation depends on the deformation direction) is a distinctive feature of many natural and structural materials: rocks, porous and cohesive bulk media, fibrous and granular composites, some metal alloys, etc. The fact that the listed materials show the heteromodular property at the stage of elastic deformation should be especially taken into account when solving problems of their shock dynamics. To describe the heteromodular behavior of an elastic medium in terms of small strains we use the physically nonlinear model of V.P. Myasnikov. The accepted assumption about the one-dimensional straining reduces the nonlinear relationship of stresses and small strains to piecewise linear equations. In the case of dynamic shock deformation, the initial nonlinearity of the model is concentrated in the equations which define the velocity of the shock wave abruptly transforming the heteromodular medium from a stretched to a compressed state. In this paper we investigate the processes of generation, motion, and possible interactions of plane one-dimensional deformation waves (including shock ones) in a heteromodular elastic half-space. The points of the half-space boundary undergo one-dimensional motions according to a given non-linear law corresponding to the “stretching-compression” mode. We suggest replacing the nonstationary boundary condition of the problem by its piecewise linear approximation and constructing a connected sequence of analytical solutions with a linear boundary condition at each local time interval. The proposed approach is the basis of the numerical solving algorithm for a boundary value problem with a given nonlinear condition. It is shown that the general solution behind the shock wave consists of several local layers, which number is related to the quantity of nodes in the piecewise linear decomposition of the boundary condition. In these layers, the compression deformation is defined by the relevant part of the boundary condition and simultaneously “stores” information on the preliminary tension, which should be considered an important feature of the heteromodular medium dynamics.

Full Text

Введение Известно, что большинство природных и конструкционных материалов обладают существенной нелинейностью связи между напряжениями и деформациями. Разномодульные среды занимают особое место в этом ряду, поскольку их нелинейность имеет практически сингулярный характер и проявляется уже на стадии упругого деформирования. В первую очередь это касается горных пород [1-5], всевозможных пористых и связных сыпучих сред [6-8], волокнистых и зернистых композитов [9-12], а также некоторых современных сплавов [13-16]. Распространенность и промышленная ценность таких материалов являются причиной активного экспериментального и теоретического изучения их физико-механических свойств. К настоящему времени разработан целый ряд нелинейных математических моделей разномодульных сред (например, [17-23]), поскольку классических представлений линейной теории здесь оказывается недостаточно. Учет физической нелинейности разномодульного материала достигается в модельных соотношениях самыми разными способами: зависимостью упругих модулей от типа напряженно-деформированного состояния [21, 22, 24], постулированием упругого потенциала особой формы [17-20, 23], построением различных реологических схем [25-29] и т.д. В [30] на основе компиляции подходов [17, 18, 20, 24] модельные соотношения разномодульной упругой среды в простейшем случае ее одномерного движения записаны в форме кусочно-линейных функций напряжения от деформации. Замена гладкой нелинейной связи на ее кусочно-линейную аппроксимацию существенно упрощает процедуру решения ряда нестационарных краевых задач одномерной динамики разномодульной среды. Так, в [31] представлены обобщенные решения с одиночными плоскими одномерными волнами деформаций, возникающими в разномодульной среде [18] при различных гладких нелинейных краевых условиях; в [32] для этой же модели получено решение одномерной задачи о гармоническом нагружении границы полупространства с чередованием ударных волн и жестких слоев; в [33-35] рассмотрено движение одномерных плоских и сферических волн в кусочно-линейной разномодульной среде при кусочно-гладком изменении граничной нагрузки с растяжения на сжатие или наоборот. Однако до сих пор для разномодульных материалов не затрагивался вопрос о нестационарном взаимодействии нелинейных волновых фронтов друг с другом, поскольку даже кусочно-линейные модельные соотношения не позволяют в этом случае обойтись без привлечения специальных приближенных методов [36]. Данную проблему можно в некоторых частных случаях снять, если в дополнение к линеаризации модельных соотношений разномодульной среды выполнить кусочно-линейную аппроксимацию нестационарных краевых условий задачи. Продемонстрируем это в настоящей статье. 1. Определяющие соотношения Запишем систему модельных уравнений динамики деформирования разномодульной упругой среды в прямоугольных декартовых координатах, считая, что процесс адиабатический, деформации малые и массовые силы отсутствуют: (1) Здесь - компоненты тензора напряжений Коши; - компоненты тензора малых деформаций; и - начальная и текущая плотность среды; и - компоненты вектора перемещений и вектора скорости перемещений точек среды; латинским индексом после запятой обозначена частная производная функции по независимой пространственной переменной , точкой - частная производная по времени . Различное сопротивление материала растяжению и сжатию в системе (1) постулируем упругим потенциалом среды [18, 20], где , - параметры Ламе линейной среды, дополнительное слагаемое с константой отвечает за проявление эффекта разномодульности при смене типа деформированного состояния. Многоточием обозначены невыписанные слагаемые неаналитичного при разложения [20] упругого потенциала в ряд по сферическим функциям, для которого в [20] дополнительно принято и показано существование пределов , . В случае одномерного движения точек среды ( , , ) упругий потенциал в форме (1) приводит к кусочно-линейной зависимости : (2) где при . Аналогично обобщенной модели [20] второе равенство (2) обеспечивает сохранение непрерывности связи в нуле. Используя (2), из системы (1) получаем уравнение движения (3) В (3) характеристическая скорость может принимать различные значения в зависимости от типа деформации: в областях сжатия среды, в областях растяжения ( при ). Решение уравнения (3) в форме Даламбера (4) с неизвестными функциями , определяется для каждой конкретной задачи с учетом заданных краевых и начальных условий. Если в обобщенном решении уравнения движения (3) возникает одномерный фронт сильного разрыва , то множество краевых условий задачи необходимо дополнить условиями совместности разрывов [37] и требованием непрерывности перемещений на поверхности : (5) В (5) и далее разрывы функций на фронте обозначены как , где , . Из условий (5) для модели (2) можно вычислить скорость сильного разрыва: (6) которая оказывается зависимой от деформированного состояния ( , ) в малых окрестностях с обеих сторон от фронта . Как будет показано далее, такая модельная особенность приводит к существенным отличиям динамики разномодульной среды от известных результатов линейной упругости даже в простейшем случае одномерных движений. Для классификации возможных сильных разрывов в обобщенном решении уравнения (3) используем терминологию, аналогичную принятой в [38]. Если и , то в среде распространяется быстрый или медленный простой разрыв ; если при , то такая волна - быстрый или медленный полусигнотон . Из (5), (6) следует, что простые разрывы и полусигнотоны в среде (2) движутся с характеристической скоростью , равной для быстрых фронтов , и для медленных , [31, 38] (знаками «+» и «-» у верхних индексов ( ) обозначим противоположные направления движения быстрых (медленных) фронтов; далее знак «+» будем опускать). Ударной волной со скоростью будем называть такой сильный разрыв, у которого . Используя (5), (6), для модели (2), (3) можно показать, что [31, 38]. Согласно приведенной классификации, простой разрыв и ударная волна распространяются только по ненулевому предварительному полю деформации . Ударная волна скачком меняет тип деформации с растяжения ( , ) на сжатие ( , ). Следует отметить, что ударный переход от сжатия к растяжению невозможен в силу неэволюционности такого разрыва [39]. Простой разрыв, в отличие от ударной волны, не изменяет тип деформированного состояния среды. Полусигнотон несет ударное граничное возмущение (растяжение или сжатие) в предварительно недеформированную область ( ). Кроме того, в разномодульной среде (2) быстрый и медленный полусигнотоны могут двигаться парой, образуя между собой жесткий слой при переходе предварительно сжатого материала в растянутое состояние [31, 32, 38]. Перечисленные свойства сильных разрывов в обобщенном решении уравнения движения (3) далее используем для анализа возможных волновых картин в задаче нестационарной одномерной динамики разномодульной среды. 2. Движение и взаимодействие плоских одномерных волн деформации в разномодульном упругом полупространстве Рассмотрим задачу об одномерном динамическом деформировании разномодульного упругого полупространства под действием нестационарного граничного нагружения в режиме «растяжение-сжатие». Считаем, что изначально среда находилась в недеформированном состоянии: , , . Пусть в момент на границу полупространства начинает действовать нагрузка, вызывающая движение граничных точек по гладкому нелинейному закону (пунктирная линия на рис. 1): Такое одномерное движение точек границы приводит к растяжению среды при и к сжатию при В [34, 35] для описанного режима граничного воздействия показано, что в момент в разномодульной среде (2) возникает ударная волна , которая движется со скоростью за первичным медленным фронтом и скачком переводит среду из растянутого состояния в сжатое. Однако временной интервал представленных в [34, 35] решений не включает момент, когда ударная волна догонит идущий впереди нее медленный фронт растяжения. Здесь мы не будем накладывать таких ограничений, так как изучение взаимодействия волновых фронтов различного типа в разномодульной среде (2) является одной из основных целей нашего исследования. Рис. 1. Заданное перемещение граничных точек полупространства Fig. 1. Displacement specified for boundary points of a half-space Для заданной гладкой функции граничных перемещений можно построить ее непрерывное кусочно-линейное приближение: (7) В (7) количество и расположение узловых точек функции выбирается в зависимости от желаемой точности аппроксимации. Примем в (7) и построим функцию в форме трех линейных участков (сплошная линия на рис. 1): (8) Константы , в (8) вычисляются согласно (7). Три узловые точки функции фиксируем при постановке задачи: и соответствуют нулям функции , , ; - момент изменения граничной нагрузки с растяжения на сжатие ( ). Положение точки может быть любым в интервале Таким образом, нестационарную начально-краевую задачу с гладкой функцией граничных перемещений разбиваем на несколько последовательных временных этапов, выбирая для каждого из них соответствующую часть кусочно-линейного краевого условия (8). Прогноз волновой картины на текущем этапе строим исходя из результатов предыдущего этапа. В дополнение к (8) в систему краевых условий каждого локального этапа включаем условия (5) на подвижных границах - фронтах сильных разрывов. Действуя по описанному алгоритму, получаем последовательность автомодельных решений, объединение которых составляет кусочно-гладкую аппроксимацию полей и во всей области деформирования на интервале . Момент , позже которого волновая картина уже не меняется при заданном формулой (8) движении граничной плоскости, является здесь точкой остановки вычислительной процедуры. Этап 1: нестационарное растяжение разномодульного полупространства при . Согласно (8) движение точек границы на интервалах времени и определяется функциями и соответственно. Результатом такого граничного воздействия являются два медленных волновых фронта (рис. 2): полусигнотон и простой разрыв . а б Рис. 2. Диаграммы сильных разрывов и областей гладкого решения: a - сценарий AB ( , ); б - сценарий BA ( , ) Fig. 2. Diagrams for strong discontinuities and regions of a smooth solution: a - AB scenario ( , ); b - BA scenario ( , ) Область 0 соответствует начальному (недеформированному) состоянию материала, в областях I и II среда подвергается растяжению ( ). Полусигнотон (волна первичного растяжения) возникает в начальный момент . Появление простого разрыва обусловловлено исключительно наличием излома функции в точке . Здесь стоит отметить, что увеличение количества фронтов на этапе 1 за счет увеличения числа узловых точек повышает точность аппроксимации всего решения, однако для наших целей достаточно одного простого разрыва . С учетом краевого условия (8) на этапе 1 при решение в областях I и II (см. рис. 2) принимает вид (9) Этап 2: нестационарное сжатие разномодульного полупространства с момента и до первого взаимодействия разрывов различного типа. В момент режим граничной нагрузки меняется с растягивающих усилий на сжимающие. Согласно (8) движение граничных точек теперь задано функцией . В результате в момент на границе возникает ударная волна (см. рис. 2), которая движется со скоростью ( ) в предварительно растянутую область II, создавая за собой новую область сжатия III ( при ). Такая волновая картина существует в полупространстве до момента , когда ударная волна догонит медленный простой разрыв (см. рис. 2). Как и перемещение , скорость входит в число неизвестных функций задачи. Eсли условие на границе задано в форме гладкого нелинейного закона , то в (5), (6) будет зависеть, кроме , еще и от неизвестного скачка . В таком случае решение даже простой одномерной краевой задачи можно получить только численно или с привлечением специальных приближенных аналитических методов [36]. Замена нелинейной функции ее кусочно-линейной аппроксимацией (8) позволяет существенно упростить процедуру решения, так как перемещение граничных точек по линейному закону инициирует ударную волну , которая движется по предварительному полю с постоянной скоростью: (10) Перемещение в форме (4), вычисляемое совместно с из системы краевых условий (5), (8), с учетом решения предыдущего этапа (9) принимает вид (11) Этап 3: эволюция ударной волны в результате ее взаимодействия с медленным простым разрывом . Область растяжения II перестает существовать при , когда ударная волна догоняет простой разрыв (см. рис. 2). Положение, которое в момент занимают волны и определяется координатой . В результате такого попутного столкновения и транформируются в два новых фронта, расходящихся от плоскости (см. рис. 2): быстрый простой разрыв и ударную волну со скоростью . Решение в новой области сжатия IV ( ) имеет вид (12) В (12) скорость ударной волны вычисляется из (5) с учетом (8), (9), (11) по формуле (13) Этап 4: эволюция ударной волны в процессе нарастающего сжатия на границе полупространства. Система разрывов существует до момента наступления одного из событий: A = { } - попутное столкновение ударной волны с медленным полусигнотоном в момент ( ); B = { } - отражение быстрого простого разрыва от границы в момент . Таким образом, дальнейшее решение может развиваться по различным сценариям: AB ( ) или BA ( ). Случай не рассматриваем ввиду малой вероятности одновременного наступления событий A и B. Опуская подробности дальнейших вычислений, представим решения для каждого из указанных вариантов. На рис. 2, a показана диаграмма сильных разрывов для случая AB ( ), на рис. 2, б - для случая BA ( ). Области, где волновые картины различаются, отмечены штриховкой. При (см. рис. 2, a) волновой пакет состоит только из быстрых разрывов , разделенных областями сжатия III, IV и V ( ), причем полусигнотон и простой разрыв есть результат события A. В случае (рис. 2b) волновая картина включает попутные разрывы , движущиеся с различными скоростями, и области с разными типами деформации: в IV и VI - сжатие ( ), в I - растяжение ( ). Таким образом, ударная волна может как закончить свое существание в самом начале этапа 4 в момент вследствие ее трансформации в быстрые разрывы и (вариант AB), так и продолжать участвовать в формировании поля деформации в областях сжатия IV, VI на всем протяжении этапа 4 вплоть до момента (вариант ВА). Это приводит к двум различным наборам краевых условий (5) на сильных разрывах для случаев AB и BA. С учетом данного факта перемещение в областях V и VI вычисляется в виде (14) Константы , в (14) известны из решения (12), (13) предыдущего этапа 3. Этап 5: продолжение сжатия на границе полупространства при или в зависимости от сценария предыдущего этапа (AB или BA). На данном этапе, начало которого совпадает с моментом (вариант AB, рис. 2, a) или (вариант BA, рис. 2, б), движение точек граничной плоскости по-прежнему определяется согласно (8) функцией . В любой момент этапа 5 волновая картина состоит из трех быстрых разрывов (попутных или разнонаправленных, рис. 2), тип деформации во всех актуальных областях гладкого решения соответствует сжатию. Передним фронтом сжатия является быстрый полусигнотон который движется в недеформированную область 0. Помимо , в формировании полей и в областях V-VIII участвуют быстрые простые разрывы , , . Неизменность типа деформации на протяжении всего этапа 5 приводит к постоянству характеристической скорости всюду между передним фронтом и границей ( ), а также к сохранению характера сильных разрывов при их взаимодействии друг с другом (при ) и с границей полупространства (при ). Это позволяет легко получить решение задачи в областях VII, VIII: (15) Как показано на рис. 2, в обоих случаях AB и BA при возникает волновой пакет из попутных быстрых разрывов , который уже не изменяется с течением времени, если граничные точки продолжают двигаться по закону . Таким образом, процедуру решения описанной задачи можно остановить в любой момент , например, руководствуясь оценкой масштаба реального граничного перемещения , достигнутого к этому моменту времени. 3. Анализ решения На рис. 3 схематически изображена последовательность моментальных кусочно-линейных диаграмм перемещения и деформации , соответствующих полученным решениям (9) - (15) при возникновении случая AB ( ). На каждой диаграмме момент принадлежит одному из описанных выше этапов деформирования, условная пространственная координата имеет свой нормирующий множитель, характерный для выбранного интервала времени: (a) , ; (б) , ; (в) , ; (г) , ; (д) , ; (е) , . a б в г д е Рис. 3. Результаты поэтапного решения краевой задачи (сценарий AB) Fig. 3. Results of the stepwise solution for the boundary value problem (AB scenario) Несомненно, рис. 3 не отражает всех возможностей развития системы разрывов и полей , в рассмотренной задаче. Так, на рис. 2 показано, что кусочно-линейная геометрия ударного фронта которая определяется значениями скоростей и может изменяться в зависимости от положения узловой точки функции в (8). Устремляя данную точку к началу процесса растяжения (моменту ), мы получаем вогнутую форму фронта , так как (см. рис. 2, a); если точка расположена ближе к моменту , то разрыв принимает выпуклую форму при (см. рис. 2, б). Аналогичного эффекта можно добиться при фиксированном положении узла путем изменения значения , т.е. фактически варьируя кривизну функции на интервале растяжения (рис. 4). Полученные с учетом (8), (10), (12)-(14) зависимости скоростей , и моментов , от величины схематически показаны на рис. 5, где , ; в области при , в при ; при (сценарий АВ), при (сценарий BA). Отметим, что значения и в общем случае не совпадают, т.е. существует узкий диапазон м, где при волновая картина на этапе 4 (нарастающее сжатие) может развиваться по сценарию BA при (заштрихованные области на рис. 5). Рис. 4. Функция граничного перемещения с вариацией кривизны на стадии растяжения Fig. 4. Boundary displacement function with curvature variation at the stretching stage Рис. 5. Зависимости параметров решения от значения Fig. 5. Dependencies of solution parameters on the value Заключение Итак, в статье показано, что кусочно-линейная аппроксимация нелинейных краевых условий динамической задачи позволяет описать процесс возникновения, движения и взаимодействия плоских одномерных волн деформаций в разномодульной упругой среде, не обращаясь к специальным методам типа метода возмущений и др. Полученное с помощью такого подхода кусочно-линейное решение моделирует нелинейный эффект, характерный для динамики деформирования нелинейно-упругих материалов и не описываемый в рамках линейной упругости - возникновение ударной волны сжатия со скоростью, зависящей от динамического поля деформаций по обе стороны от волнового фронта. Данный подход можно использовать для решения нестационарных краевых задач упругой динамики нелинейных и разнопрочных сред с более сложными граничными условиями, чем представленное здесь (например, режимы циклической нагрузки-разгрузки или вибрации), при этом точность решения может быть улучшена путем добавления дополнительных узловых точек аппроксимирующей кусочно-линейной граничной функции. Также результаты работы могут послужить теоретическим дополнением для экспериментальных исследований нелинейного поведения реальных природных и конструкционных материалов под действием высокоскоростных динамических нагрузок.

About the authors

O V Dudko

Institute of Automation and Control Processes

A A Lapteva

Institute of Automation and Control Processes

V E Ragozina

Institute of Automation and Control Processes

References

  1. Капустянский С.М. Анизотропия геоматериалов // Итоги науки и техники. Механика деформируемого твердого тела. - 1987. - Т. 18. - С. 53-113.
  2. Ставрогин А.Н., Протосеня А.Г. Прочность горных пород и устойчивость выработок на больших глубинах. - М.: Недра, 1985. - 271 с.
  3. Experimental investigation of the anisotropic mechanical behavior of phyllite under triaxial compression / Guowen Xu, Chuan He, Ang Su, Ziquan Chen // International Journal of Rock Mechanics and Mining Sciences. - 2018. - Vol. 104. - Р. 100-112. doi: 10.1016/j.ijrmms.2018.02.017
  4. Динамическое деформирование мягких грунтовых сред: экспериментальные исследования и математическое моделирование / В.В. Баландин [и др.] // МТТ. - 2015. - № 3. - С. 69-77.
  5. Баклашов И.В., Картозия Б.А. Механика горных пород. - М.: Недра, 1975. - 271 с.
  6. Mechanical properties and non-homogeneous deformation of open-cell nickel foams: application of the mechanics of cellular solids and of porous materials / Badiche X. [et al.] // Materials Science and Engineering. - 2000. - Vol. 289. - Р. 276-288. doi: 10.1016/S0921-5093(00)00898-4
  7. Краснощеков П.И., Федотов А.Ф. Упругие модули изотропных порошковых и пористых материалов // Вестн. Самар. гос. техн. ун-та. Сер. Физ.-мат. науки. - 2006. - Вып. 43. - С. 81-87. doi: 10.14498/vsgtu455
  8. Федотов А.Ф Прогнозирование эффективных модулей упругости пористых композиционных материалов // Изв. вузов. Порошковая металлургия и функциональные покрытия. - 2015. - № 1. - С. 32-37. doi: 10.17073/1997-308X-2015-1-32-37
  9. Капитонов А.М., Редькин В.Е. Физико-механические свойства композиционных материалов. Упругие свойства. - Красноярск: Изд-во Сиб. федер. ун-та, 2013. - 532 с.
  10. Березин А.В. Влияние повреждений на деформационные и прочностные характеристики твердых тел. - М.: Наука, 1990. - 135 с.
  11. Бессонов Д.Е., Зезин Ю.П., Ломакин Е.В. Разносопротивляемость зернистых композитов на основе ненасыщенных полиэфиров // Изв. Сарат. ун-та. Сер. Математика. Механика. Информатика. - 2009. - Т. 9, вып. 4, ч. 2. - С. 9-13.
  12. Lomakin E.V., Fedulov B.N. Nonlinear anisotropic elasticity for laminate composites // Meccanica. - 2015. -Vol. 50. - Iss. 6. - P. 1527-1535. doi: 10.1007/s11012-015-0104-5
  13. Ambroziak A., Kłosowski P., Pyrzowski L. Elasto-viscoplastic properties of aa2017 aluminium alloy // Task Quarterly. - 2015. - No. 1. - P. 5-20.
  14. Моделирование поведения титанового сплава вт20 при ударном взаимодействии / Е.В. Ломакин [и др.] // Вестн. Нижегород. ун-та им. Н.И. Лобачевского. - 2011. - № 1. - С. 129-132.
  15. Баев Л.В. Распространение продольных и поперечных волн в упругой разномодульной среде // Прикладная механика и техническая физика. - 2009. - Т. 50, № 4 (296). - С. 176-182.
  16. Механические свойства конструкционных материалов при сложном напряженном состоянии / А.А. Лебедев [и др.]. - Киев: Ин Юре, 2003. - 540 с.
  17. Ломакин Е.В., Работнов Ю.Н. Соотношения теории упругости для изотропного разномодульного тела // Изв. АН СССР. Механика твердого тела. - 1978. - № 6. - С. 29-34.
  18. Мясников В.П., Топалэ В.И. Моделирование сейсмической анизотропии в литосфере как разномодульном упругом теле // Изв. АН СССР. Физика Земли. - 1987. - № 5. - С. 22-30.
  19. Ломакин Е.В. Определяющие соотношения деформационной теории для дилатирующих сред // Изв. АН СССР. Механика твердого тела. - 1991. - № 6. - С. 66-75.
  20. Мясников В.П., Олейников А.И. Основы механики гетерогенно-сопротивляющихся сред. - Владивосток: Дальнаука, 2007. - 172 с.
  21. Ломакин Е.В. Механика сред с зависящими от вида напряженного состояния свойствами // Физическая мезомеханика. - 2007. - Т. 10, № 5. - С. 41-52.
  22. Пахомов Б.М. Вариант модели изотропного разномодульного материала // Вестн. Моск. гос. техн. ун-та им. Н.Э. Баумана. Сер.: Машиностроение. - 2017. - № 6 (117). - С. 35-48. doi: 10.18698/0236-3941-2017-6-35-48
  23. Пирогов С.А. Модель динамического деформирования и разрушения разносопротивляющихся материалов при интенсивных ударных воздействиях // Проблемы прочности и пластичности. - 2012. - Вып. 74. - С. 40-48.
  24. Амбарцумян С.А., Хачатрян А.А. Основные уравнения теории упругости для материалов, разносопротивляющихся растяжению и сжатию // Изв. АН СССР. МТТ. - 1966. - № 2. - С. 44-53.
  25. Sadovskaya O., Sadovskii V. Mathematical modeling in mechanics of granular materials. - Germany: Heidelberg, Springer, 2012. - 390 p. doi: 10.1007/978-3-642-29053-4
  26. Садовский В.М., Садовская О.В. Анализ деформации пористой среды с учетом схлопывания пор // Прикладная механика и техническая физика. - 2016. - Т. 57, № 5. - С. 53-65. doi: 10.15372/PMTF20160507
  27. Кузоватова О.И., Садовский В.М. Численное исследование задачи о продавливании связной сыпучей среды // Вестн. СибГАУ. - 2009. - № 4. - С. 21-25.
  28. Садовская О.В., Садовский В.М. Реологические модели одноосного деформирования пористых сред // Вестн. КрасГУ. Физико-математические науки. - 2006. - № 9. - C. 202-206.
  29. Sadovskii V.M., Sadovskaya O.V., Luk'yanov A.A. Radial expansion of a cylindrical or spherical cavity in an infinite porous medium // Journal of Applied Mechanics and Technical Physics. - 2014. - Vol. 55. - Iss. 4. - P. 689-700. doi: 10.1134/S0021894414040154
  30. Дудко О.В., Лаптева А.А., Чигарев А.В. К построению математической модели разномодульной изотропно-упругой среды // Вестник ЧГПУ им. И.Я. Яковлева. Сер. Механика предельного состояния. - 2013. - № 2 (16). - С. 61-67.
  31. Дудко О.В., Лаптева А.А., Семенов К.Т. О распространении плоских одномерных волн и их взаимодействии с преградами в среде, по-разному сопротивляющейся растяжению и сжатию // Дальневосточ. мат. журн. - 2005. - Т. 6, № 1-2. - С. 94-105.
  32. Gavrilov S.N., Herman G.C. Wave propagation in a semi-infinite heteromodular elastic bar subjected to a harmonic loading // Journal of Sound and Vibration. - 2012. - Vol. 331 (20). - P. 4464-4480. doi: 10.1016/j.jsv.2012.05.022
  33. Дудко О.В., Лаптева А.А., Рагозина В.Е. О возникновении одномерных сферических центрированных волн в разномодульной упругой среде при ее динамическом деформировании // Вестн. ЧГПУ им. И.Я. Яковлева. Сер.: Механика предельного состояния. - 2016. - № 1 (27). - С. 64-69.
  34. Рагозина В.Е., Дудко О.В. Движение сходящихся сферических волн деформаций в разномодульной упругой среде // Прикладная механика и техническая физика. - 2016. - Т. 57, № 4. - С. 149-157. doi: 10.15372/PMTF20160414
  35. Dudko O.V., Ragozina V.E. On the motion of shock waves at a constant speed in multimodulus elastic media // Mechanics of Solids. - 2018. - Vol. 53. - Iss. 1. - P. 111-119. doi: 10.3103/S0025654418010132
  36. Найфэ А.Х. Методы возмущений: пер. с англ. - М.: Мир, 1976. - 456 с.
  37. Бленд Д. Нелинейная динамическая теория упругости. - М.: Мир, 1972. - 183 с.
  38. Маслов В.П., Мосолов П.П. Общая теория уравнений движения разномодульной упругой среды // ПММ. - 1985. - Т. 49, вып. 3. - С. 419-437.
  39. Куликовский А.Г., Свешникова Е.И. Нелинейные волны в упругих средах. - М.: Московский лицей, 1998. - 412 с.

Statistics

Views

Abstract - 429

PDF (Russian) - 1236

Cited-By


PlumX


Copyright (c) 2019 Dudko O.V., Lapteva A.A., Ragozina V.E.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies