Гранично-элементное решение трехмерных динамических задач анизотропной вязкоупругости и изотропной поровязкоупругости
Автор: Марков Иван Петрович, Ипатов Александр Александрович, Белов Александр Александрович, Литвинчук Светлана Юрьевна
Журнал: Вычислительная механика сплошных сред @journal-icmm
Статья в выпуске: 4 т.9, 2016 года.
Бесплатный доступ
Обсуждается динамическое поведение анизотропных вязкоупругих и изотропных поровязкоупругих тел. Поровязкоупругая постановка опирается на полную модель насыщенной пороупругой среды Био. Теория Био является расширением классической теории упругости на случай двухфазной среды, состоящей из упругого скелета с порами и наполнителя. Применяется принцип соответствия упругой и вязкоупругой реакций. Для описания вязкоупругих свойств скелета пористого материала используется модель стандартного вязкоупругого тела. Приводится система дифференциальных уравнений для полной модели Био в преобразованиях Лапласа. Решение исходной задачи строится в пространстве преобразований Лапласа с последующим обращением интегрального преобразования с помощью численного алгоритма. Для отыскания решения в изображениях по Лапласу записывается система граничных интегральных уравнений прямого подхода. Рассматриваются регуляризованные граничные интегральные уравнения, и производится согласованное гранично-элементное разбиение для получения дискретных аналогов. Коллокационные точки решения граничного интегрального уравнения совпадают с узлами интерполяции неизвестных граничных функций. Анизотропные фундаментальные решения представляются как сумма статической и динамической частей, которые записываются в виде интегралов по единичной окружности и единичной полусфере соответственно. Численное обращение преобразования Лапласа реализуется на основе шагового по времени метода на узлах схемы Рунге-Кутты. На решении, найденном методом граничных элементов, продемонстрировано влияние вязкоупругих свойств поровязкоупругого и анизотропного вязкоупругого материалов на амплитуды и формы откликов при переходе с мгновенных модулей на длительные. Приведены численные решения задач определения волновых полей в Г-образном анизотропном вязкоупругом теле при действии силы на его торец и в поровязкоупругом кубе, содержащем сферическую полость, подверженную равномерно распределенному внутреннему давлению.
Трехмерные краевые задачи, метод граничных элементов, поровязкоупругость, анизотропная вязкоупругость, обращение преобразования лапласа, шаговый метод, схема рунге-кутты
Короткий адрес: https://sciup.org/14320826
IDR: 14320826 | УДК: 539.3 | DOI: 10.7242/1999-6691/2016.9.4.39
Boundary element formulation for 3d dynamic problems of anisotropic viscoelasticity and isotropic poroviscoelasticity
The dynamic behavior of anisotropic viscoelastic and poroviscoelastic solids is considered. A poroviscoelastic formulation is based on the Biot model of fully saturated poroelastic media. The elastic-viscoelastic correspondence principle is applied to describe the viscoelastic properties of a porous material skeleton. The system of differential equations of the full Biot model in Laplace transforms and formulas for elastic modules are given. A solution to the original problem is constructed using Laplace transforms, and numerical inversion yields the solution in the time domain. The system of direct boundary integral equations is introduced, and the system of regularized boundary integral equations is considered. Discrete analogues are obtained through mixed boundary element discretizations. The collocation points of the boundary integral equation coincide with the interpolation nodes of unknown boundary functions. Anisotropic fundamental solutions are represented as a sum of static and dynamic parts expressed in terms of integrals over the unit circle and the half of the unit sphere, respectively. Numerical inversion of Laplace transform is realized at the nodes of the Runge-Kutta scheme using the time-step method. The boundary element analysis is conducted to demonstrate the influence of the viscoelastic properties of poroviscoelastic and anisotropic viscoelastic materials on the amplitude and shape of responses during the transition from instantaneous to equilibrium moduli. A numerical solution to the problem of force acting on a prismatic anisotropic viscoelastic solid is given. The problem of growing pressure in a spherical cavity inside a poroviscoelastic solid is also numerically solved.
Список литературы Гранично-элементное решение трехмерных динамических задач анизотропной вязкоупругости и изотропной поровязкоупругости
- Frenkel J. On the theory of seismic and seismoelectric phenomena in a moist soil//J. Eng. Mech. -2005. -Vol. 131, no. 9. -P. 879-887.
- Biot M.A. General theory of three-dimensional consolidation//J. Appl. Phys. -1941. -Vol. 12, no. 2. -P. 155-164.
- Schanz M. Wave propagation in viscoelastic and poroelastic continua. -Berlin: Springer, 2001. -170 p.
- De Boer R. Highlights in the historical development of the porous media theory: toward a consistent macroscopic theory//Appl. Mech. Rev. -1996. -Vol. 49, no. 4. -P. 201-262.
- Nikolaevskiy V.N. Biot-Frenkel poromechanics in Russia (Review)//J. Eng. Mech. -2005 -Vol. 131, no. 9. -P. 888-897.
- Garg S.K., Nayfeh A.H., Good A.J. Compressional waves in fluid-saturated elastic porous media//J. Appl. Phys. -1974. -Vol. 45, no. 5. -P. 1968.
- Beskos D.E. Boundary element methods in dynamic analysis: Part II (1986-1996)//Appl. Mech. Rev. -1997. -Vol. 50, no. 3. -P. 149-197.
- Carini A., Gioda G. A boundary integral equation technique for visco-elastic stress analysis//Int. J. Numer. Anal. Meth. Geomech. -1986. -Vol. 10, no. 6. -P. 585-608.
- Hwu C., Chen Y.C. Analysis of defects in viscoelastic solids by a transformed boundary element method//Procedia Engineering. -2011. -Vol. 10. -P. 3038-3043.
- Sim W.J., Kwak B.M. Linear viscoelastic analysis in time domain by boundary element method//Comput. Struct. -1988. -Vol. 29, no. 4. -P. 531-539.
- Баженов В.Г., Игумнов Л.А. Методы граничных интегральных уравнений и граничных элементов в решении задач трехмерной динамической теории упругости с сопряженными полями. -М.: Физматлит, 2008. -352 с.
- Игумнов Л.А., Ратаушко Я.Ю. Шаговый метод численного обращения преобразования Лапласа на узлах схемы Рунге-Кутты с использованием переменного шага интегрирования//Проблемы прочности и пластичности. -2013. -Т. 75, № 4. -С. 280-287.
- Wang C.-Y., Achenbach J.D. Elastodynamic fundamental solutions for anisotropic solids//Geophys. J. Int. -1994. -Vol. 118, no. 2. -P. 384-392.
- Wang C.-Y., Achenbach J.D. Three-dimensional time-harmonic elastodynamic Green’s functions for anisotropic solids//Proc. R. Soc. A. -1995. -Vol. 449, no. 1937. -P. 441-458.
- Banjai L., Schanz M. Wave propagation problems treated with convolution quadrature and BEM//Lecture Notes in Applied and Computational Mechanics. -2012. -Vol. 63. -P. 145-184.