Fundamental Study on Nonlinear Finite Element Method to Achieve High Accuracy

实现高精度非线性有限元方法的基础研究

• 批准号: 02650062
• 负责人: WATANABE Osamu
• 金额: $ 1.15万
• 依托单位: University of Tsukuba
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for General Scientific Research (C)
• 财政年份: 1990
• 资助国家: 日本
• 起止时间: 1990 至 1992
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-02650062/
• 关键词: Finite Element Method; Constitutive Equation; Variational Principle; Numerical Analysis; Hypoelasticity; Hyperelasticity; 有限変形理論; 混合法; 有限変形論; 大変形シミュレ-ション

The finite element method is widely used in the structural analysis of engineering products for the reasonable design and developments. The present study is concerned with the constitutive equation and variational principles and their effects on numerical analysis in order to achieve high accuracy in the non-linear finite element analysis during three years from 1990 to 1992. Hypoelasticity, defined as the rate type constitutive equation, is used to describe plastic behaviors of metal alloys, on the other hand, hyperelasticity, defined as the potential type constitutive equation, is used for rubber materials. Hypoelasticity widely used with Jaumann stress derivative gives oscillation in stress-deformation curve if linear kinematic hardening law is adopted in simple shear problem. Hyperelasticity, based on Green-Lagrange strain measured from the undeformed configuration, cannot predict accurately the experimental fact that shear stress respond behind shear deformation in simple shear problem. For the analysis of hyperelasticity, the employed mixed variational principle does not include the initial stress term, which significantly affects bifurcation analysis. In order to overcome such shortcomings in constitutive equation and mixed variational principles, the following points are studied consecutively.(1) By using rotationless strain regarded as the generalized logarithmic strain, hyperelastic constitutive equation is proposed.(2) Theoretical aspects of hypoelasticity is studied from the viewpoint of corotational equation and non-coaxiality relating to intermediate configurations.(3) Constitutive equation for large plastic strain is formulated based on internal time theory and compared to the experiments.(4) Mixed variational principles taking account of the initial stress term and stability of solutions are studied theoretically.(5) Numerical analysis of hypo and hyperelasticity is carried out by using rotationless strain.

有限元方法被广泛应用于工程产品的结构分析中，以便进行合理的设计和开发。为了提高非线性有限元分析的精度，本文对1990~1992年三年间的本构方程和变分原理及其对数值分析的影响进行了研究。亚弹性被定义为速率型本构方程，用于描述金属合金的塑性行为，而超弹性被定义为势型本构方程，用于橡胶材料。在简单剪切问题中，如果采用线性运动硬化定律，亚弹性中广泛使用的Jaumann应力导数会使应力-变形曲线出现振荡。基于未变形构型的格林-拉格朗日应变的超弹性不能准确地预测简单剪切问题中剪切变形后的剪切应力响应这一实验事实。对于超弹性分析，所采用的混合变分原理没有考虑初始应力项，这对分叉分析有很大影响。为了克服本构方程和混合变分原理中的这些不足，本文对以下几点进行了连续的研究：(1)将无旋转应变作为广义对数应变，提出了超弹性本构方程。(2)从共转方程和与中间构型相关的非共轴性的角度研究了亚弹性的理论问题。(3)基于内时理论建立了塑性大应变的本构方程，并与实验结果进行了比较。(4)从理论上研究了考虑初应力项和解的稳定性的混合变分原理。(5)利用无旋转应变进行了亚弹性和超弹性的数值分析。

项目成果

木下 伸・渡部 修: "非回転ひずみを用いた超弾性体の有限要素解析(2次元平面ブロックのせん断変形)" 日本機械学論文集 A編. 58-556. 2328-2335 (1992)

Shin Kinoshita 和 Osamu Watanabe：“使用非旋转应变的超弹性体有限元分析（二维平面块的剪切变形）”日本机械工程师学会，第 A. 58-556 卷（1992 年）。

郭 智宏,渡部 修: "亜弾性構成式と塑性硬化則が数値シミュレ-ションに及ぼす影響(2次元平面ブロックのせん断変形)" 日本機械学会論文集.

Tomohiro郭，OsamuWatanabe：“亚弹性本构方程和塑性硬化定律对数值模拟的影响（二维平面块的剪切变形）”日本机械工程学会会议录。

渡部 修,原田 隆: "Green型構成式を用いた平面ブロックの弾性大変形シミュレ-ション" 日本機械学会論文集Vo.91ー0142.

Osamu Watanabe、Takashi Harada：“使用格林型本构方程模拟平面块的大弹性变形”日本机械工程师学会会议录，卷 91-0142。

木下 伸,渡部 修: "非回転ひずみを用いた超弾性体の有限要素解析(2次元平面ブロックのせん断変形)" 日本機械学会論文集,A編. 58. 2328-2335 (1992)

Shin Kinoshita、Osamu Watanabe：“使用非旋转应变的超弹性体有限元分析（二维平面块的剪切变形）”日本机械工程师学会会刊，A 版 58. 2328-2335（1992 年）。

Watanabe, O.: "Hyperelasticity Theory based on Rotationless Strain (Constitutive Equation for Incompressible Materials)" JSME International Journal, Series I. Vol.35, No.1. 53-61 (1992)

Watanabe, O.：“基于无旋转应变的超弹性理论（不可压缩材料的本构方程）” JSME 国际期刊，系列 I. 第 35 卷，第 1 期。