Численный алгоритм решения задачи о больших упругопластических деформациях МКЭ

Автор: Давыдов Руслан Лаврентьевич, Султанов Ленар Усманович

Журнал: Вестник Пермского национального исследовательского политехнического университета. Механика @vestnik-pnrpu-mechanics

Статья в выпуске: 1, 2013 года.

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

Работа посвящена методике исследования конечных упругопластических деформаций. В качестве тензоров, описывающих деформацию и скорость деформации, используются левый тензор Коши–Грина, тензор пространственного градиента скорости и тензор деформации скорости. Вводится удельная потенциальная энергия деформации, которая зависит от левого тензора Коши–Грина. Рассмотрен изотропный материал. Напряженное состояние описывается тензором истинных напряжений Коши–Эйлера, который определяется в актуальном состоянии. Получены линеаризованные определяющие соотношения упругого деформирования в виде зависимости производной Трузделла тензора напряжении Коши–Эйлера от деформации скорости. В рамках теории течения используются аддитивное представление для полной деформации скорости. Предполагается справедливость ассоциированного закона течения. Критерием упругого деформирования является условие Мизеса–Губера. Алгоритм исследования основан на методе последовательных нагружений. В качестве базового уравнения принимается уравнение мощностей в актуальном состоянии. После линеаризации получена разрешающая система линейных уравнений, где неизвестным является приращение перемещений в текущем временном слое. При моделировании пластических деформаций применяется метод проецирования напряжений на поверхность текучести с итерационным уточнением текущего напряженно-деформированного состояния, основанным на введение в разрешающие уравнение мощности дополнительных напряжений. В качестве примера рассмотрено построение алгоритма решения для материала второго порядка. Выбрано соответствующее выражение потенциала упругих деформаций, критерием пластического течения служит условие Губера–Мизеса с изотропным упрочнением. Получены линеаризированные определяющие соотношения. Численная реализация основана на методе конечных элементов. Используется восьмиузловой конечный элемент. Созданный алгоритм исследования больших упругопластических деформаций опробован на решении тестовой задачи о растяжении круглого стержня с образованием шейки. Приводятся результаты решения и сравнение с результатами, полученными другими авторами. Также исследовалось деформирование квадратной плиты под действием внутреннего давления.

Еще

Большие деформации, нелинейная упругость, пластичность, метод конечных деформаций

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

IDR: 146211468   |   УДК: 539.3

Numerical algorithm of solving the problem of large elastic-plastic deformation by fem

A numerical algorithm of the investigation of stress-strain state of the elastic–plastic solids with large deformations is described. The left Cauchy-Green tensor and velocity gradient tensor are used as the tensors describing the deformation and deformation rate. The potential of the elastic strain the specific strain energy, which depends on the left Cauchy-Green tensor is introduced. An isotropic material is considered. The state of stress is described by the Cauchy stress tensor. The linearized constitutive equations of elastic deformation are obtained as a function of the derivative of the Truesdell stress rate on the strain rate. The theory of flow and an additive representation of the total deformation rate are used. The von Mises yield criterion is applied. The research algorithm is based on the incremental method. The principle of virtual work in terms of the virtual velocity is used. After linearization the system of linear equations is obtained, where the unknown is the increment of displacement in the current state. The radial return method with an iterative refinement of the current mode of deformation, based on the introduction to the governing equations of additional power voltages is applied. As an example the potential of elastic deformation is considered. The von Mises yield criterion with isotropic hardening is used. The linearized constitutive relations is obtained. The numerical implementation is based on the finite element method. An 8-node brick element is used. Developed algorithm of investigation of large elastoplastic deformations is tested on the solution of the necking of circular bar problem. The results of solutions and comparison with results obtained by other authors is reduced. Also the deformation of a square plate under internal pressure is investigated.

Еще