有限要素法
微分方程式 |
---|
![]() |
分類 |
解 |
計算物理学 |
---|
![]() |
データ解析 · 可視化 |
有限要素法(ゆうげんようそほう、英語: Finite Element Method, FEM)は数値解析手法の一つ。 解析的に解くことが難しい微分方程式の近似解を数値的に得る方法の一つであり[1][2][3][4]、 Turner-Clough-Martin-Toppによって導入された[5]。 方程式が定義された領域を小領域(要素)に分割し、各小領域における方程式を比較的単純で共通な補間関数で近似する[3]。 構造力学分野で発達し[6]、他の分野でも広く使われている手法である。その背景となる理論は、関数解析(リースの表現定理、ラックス=ミルグラムの定理など)と結びついて、数学的に整然としている[1][2][注釈 1][7]。
FEMを用いて現象を研究・分析することを「有限要素解析(FEA)」と呼ぶことがある。
特徴
[編集]- 各小領域内を一次関数で補間(近似空間が元の解空間の部分空間になる場合はある種の射影を求めることになる)した場合、全領域では適切なノルムに対して最良近似であることが示される[注釈 2]。
- 線形問題[8]・非線形問題[9][10]・動的解析[11][12][13]など、さまざまな問題に対応できる。これは、近似方程式の作り方や領域形状について、自由度が高いことに起因する[1][2][3]。
- FEMでは、変動微分法を用いて、誤差関数を最小化することで解を近似する。
メッシュ
[編集]- 領域全体をより単純な部分に細分化することには、以下のような利点がある。
- 複雑な形状を正確に表現できる。異種材料の特性を与えることができる。トータルソリューションの表現が容易。局所的な効果を把握できる。
- 計算のため小領域に分割することを「メッシュを切る」という言い方をする。メッシュの適切さは解析結果の正確さに大きく影響する。
- FEMは工学的な解析を行うための計算ツールとしても広く用いられる。FEMを用いたCAEソフトウエアでは、複雑な幾何形状を微小要素に分割するメッシュを自動生成する機能をもつものがある。ただしそうした場合も形状の特異点の扱いに気をつける必要がある。(構造解析における隅部など)
- 自動車や石油パイプラインのような複雑な解析や、境界が移動する固体反応のように領域が変化する場合や、必要な精度が領域全体で変化する場合や、解が滑らかでない場合などにはFEMは特に適している。FEMではメッシュの粗さを調整することで解析の計算コストを抑えることができる。(変化が大きく重要な部分はメッシュを細かくし予測精度を高め、変化の小さい部分は予測精度を上げる必要がないのでメッシュを粗くし、計算量を減らす。)例えば、自動車の正面衝突シミュレーションでは、車の前部など重要な領域はメッシュを細かくし、後部のメッシュを粗くする。また数値気象予報では、比較的穏やかな場所よりも、高度に非線形な現象が発生している場所(大気中の熱帯低気圧や海洋の渦など)を正確に予測することが重要になる。
アルゴリズム
[編集]- 解析対象領域内で成り立つ方式(ポアソン方程式など)に対してある重み関数の積を施し、それを領域内で積分した弱形式を形成する。
- 解析領域内部を小さな有限範囲の要素に分割する。一般的に、要素はその境界に節点が配置され、要素内部の物理量は各節点に対応する形状関数と節点の値の積の和として表現される[注釈 3]。
- 有限要素法では多くの種類の要素が定式化されていて、問題に依って使い分けられるようになっている。要素の種類の違いは、要素の形状、要素内での解の近似に用いる多項式の次数や、隣り合う要素の間の境界での近似解の連続性などによる。
- 解析領域全体の弱形式は積分で表されるので、それぞれの要素内の積分の総和として表すことができる。つまり、各要素の節点における未知数に対してこの積分を適用することによって、各要素の係数行列(連立一次方程式の左辺行列)を作成する(未知数は変位、速度、圧力など。右辺ベクトルも同時に形成される)。この係数行列は要素剛性行列と呼ばれる。
- 実際の複雑な問題では要素領域内に対する積分の値を解析的な式計算で求めるのは難しいので、領域の補間関数の次数に応じてガウス・ルジャンドル法などの数値積分を用いて近似することが多い[3]。
- 各要素における係数行列(要素係数行列)の総和を取って領域全体の係数行列(全体剛性行列[14]と呼ばれる)を作成し、解を求めることができる。
多くの場合に有限要素法では、近似解を求めることが連立一次方程式を解くことに帰着される (つまり最終的には数値線形代数の知識が必要になる)[3]。得られる全体の係数行列は一般に疎行列となる。使用記憶領域の削減と計算速度向上のため、行列のデータ構造には様々な形式が用いられ、その格納形式に対応して効率よく解くソルバーが存在する。たとえば、直接法で解く場合のスカイライン法などが知られている[15]。
形状関数
[編集]形状関数とは、節点における物理量(変位など)から要素内の物理量を内挿するために用いられる関数である。たとえば四面体一次要素の場合、4つの頂点に節点i = 1, ... , 4 がとられ、節点i に対する形状関数Ni とそれぞれの点における物理量ui を用いて、要素内の任意の点 p における物理量up は形状関数の線形結合として
と表される。
形状関数Ni には、
- 節点i の位置においてNi = 1
- それ以外の節点位置においてNi = 0
という性質がある。
複雑な構造物を小さな要素の集合体として、(静的解析の場合)一次方程式に各節点の変位量の境界条件(ディリクレ境界条件やノイマン境界条件等)を代入して解く。
対象の構造に外力が加わって変形する場合などを解析する際に、構造解析には大きく分けて、変位を未知数にとる変位法と応力を未知数にとる応力法があり、有限要素構造解析では変位法が主流である。その理由は、応力法に比べてアルゴリズムが機械的に実行でき、プログラミングに適しているからである。機械設計分野ではCADモデルを用いた解析が浸透している[16]。一例を挙げれば、日産自動車は、厚木にあるNTC(日産テクニカルセンター)で、スーパーコンピューターを用いて自社の製品を解析している。
その他の分野への応用
[編集]構造解析では使用している式に意味づけをしているが、その他の分野では手法として使用することが多い。電子状態計算(→実空間法[17])・電気工学[18]・電磁場解析[19]・流体解析[20][21][22]など、微分方程式で記述されるあらゆる場の問題に適用可能であって、近年ではそれらの連成解析[23](流体構造連成、電磁場構造解析など)も盛んに研究されている。 また、従来は取扱いが難しかったクラックや大変形問題に対して、格子を用いないメッシュフリー法の研究も行われている[24][25][26]。
関連項目
[編集]手法
[編集]ソフトウェア
[編集]脚注
[編集]注釈
[編集]出典
[編集]- ^ a b c Brenner, S., & Scott, R. (2007). The mathematical theory of finite element methods (Vol. 15). Springer Science & Business Media.
- ^ a b c Oden, J. T., & Reddy, J. N. (2012). An introduction to the mathematical theory of finite elements. Courier Corporation.
- ^ a b c d e 山本哲朗『数値解析入門』(増訂版)サイエンス社〈サイエンスライブラリ 現代数学への入門 14〉、2003年6月。ISBN 4-7819-1038-6。
- ^ Ciarlet, P. G. (2002). The finite element method for elliptic problems (Vol. 40). SIAM.
- ^ Clough, R. W., Martin, H. C., Topp, L. J., & Turner, M. J. (1956). Stiffness and deflection analysis of complex structures. Journal of the Aeronautical Sciences, 23(9).
- ^ a b Zienkiewicz, O. C., & Taylor, R. L. (2005). The finite element method for solid and structural mechanics. Elsevier.
- ^ 桂田祐史、Poisson方程式に対する有限要素法の解析超特急
- ^ Johnson, C., Navert, U., & Pitkaranta, J. (1984). Finite element methods for linear hyperbolic problems. CMAME, 45, 285-312.
- ^ Wriggers, P. (2008). Nonlinear finite element methods. Springer Science & Business Media.
- ^ Akyuz, F. A., & Merwin, J. E. (1968). Solution of nonlinear problems of elastoplasticity by finite element method. AIAA Journal, 6(10), 1825-1831.
- ^ Bathe, K. J., Ramm, E., & Wilson, E. L. (1975). Finite element formulations for large deformation dynamic analysis. International journal for numerical methods in engineering, 9(2), 353-386.
- ^ Jonker, B. (1990). A finite element dynamic analysis of flexible manipulators. The International Journal of Robotics Research, 9(4), 59-74.
- ^ Hughes, T. J. (2012). The finite element method: linear static and dynamic finite element analysis. Courier Corporation.
- ^ 土倉泰. (2005). 全体剛性行列を用いた粒状要素解析. 応用力学論文集, 8, 549-556.
- ^ Elwi, A. E., & Murray, D. W. (1985). Skyline algorithms for multilevel substructure analysis. International journal for numerical methods in engineering, 21(3), 465-479.
- ^ Foucault, G., Cuillière, J. C., François, V., Léon, J. C., & Maranzana, R. (2008). Adaptation of CAD model topology for finite element analysis. Computer-Aided Design, 40(2), 176-196.
- ^ Garcıa-Cervera, C. J. (2007). An efficient real space method for orbital-free density-functional theory. Commun. Comput. Phys, 2(2), 334-357.
- ^ 中田高義, & 高橋則雄. (1980). 電気工学のための有限要素法. 電氣學會雜誌, 100(1), 45-48.
- ^ Salon, S. J. (1995). Finite element analysis of electrical machines (Vol. 101). Boston: Kluwer Academic Publishers.
- ^ Reddy, J. N., & Gartling, D. K. (2010). The finite element method in heat transfer and fluid dynamics. CRC Press.
- ^ Chung, T. J. (1978). Finite element analysis in fluid dynamics. NASA STI/Recon Technical Report A, 78.
- ^ Sackinger, P. A., Brown, R. A., & Derby, J. J. (1989). A finite element method for analysis of fluid flow, heat transfer and free interfaces in Czochralski crystal growth. International journal for numerical methods in fluids, 9(4), 453-492.
- ^ 張群, & 久田俊明. (2001). 流体・構造連成有限要素解析における連成手法に関する検討. 日本機械学会論文集 A 編, 67(662), 1555-1562.
- ^ Liu, G. R. (2009). Meshfree methods: moving beyond the finite element method. Taylor & Francis.
- ^ Yagawa, G., & Yamada, T. (1996). Free mesh method: a new meshless finite element method. Computational Mechanics, 18(5), 383-386.
- ^ Puso, M. A., Chen, J. S., Zywicz, E., & Elmer, W. (2008). Meshfree and finite element nodal integration methods. International Journal for Numerical Methods in Engineering, 74(3), 416-446.
- ^ Y. Renard, K. Poulios GetFEM: Automated FE modeling of multiphysics problems based on a generic weak form language. Preprint, https://hal.archives-ouvertes.fr/hal-02532422/document .
- ^ 大塚厚二, & 高石武史. 日本応用数理学会監修, 有限要素法で学ぶ現象と数理―FreeFem++ 数理思考プログラミング―,(2014).
- ^ 大塚厚二. (2001). 数理指向有限要素計算プログラミング言語 FreeFem++. 応用数理, 11(3), 249-252.
- ^ Dupont, T., Hoffman, J., Johnson, C., Kirby, R. C., Larson, M. G., Logg, A., & Scott, L. R. (2003). The fenics project. Chalmers Finite Element Centre, Chalmers University of Technology.
- ^ Logg, A., Mardal, K. A., & Wells, G. (Eds.). (2012). Automated solution of differential equations by the finite element method: The FEniCS book (Vol. 84). Springer Science & Business Media.
- ^ Kwon, Y. W., & Bang, H. (2018). The finite element method using MATLAB. CRC Press.
- ^ Pepper, D. W., & Heinrich, J. C. (2017). The finite element method: basic concepts and applications with MATLAB, MAPLE, and COMSOL. CRC Press.
- ^ Alberty, J., Carstensen, C., Funken, S. A., & Klose, R. (2002). Matlab implementation of the finite element method in elasticity. Computing, 69(3), 239-263.
- ^ Dabrowski, M., Krotkiewski, M., & Schmid, D. W. (2008). MILAMIN: MATLAB‐based finite element method solver for large problems. Geochemistry, Geophysics, Geosystems, 9(4).
- ^ Chen, L. (2008). iFEM: an innovative finite element methods package in MATLAB. Preprint, University of Maryland.
- ^ Alberty, J., Carstensen, C., & Funken, S. A. (1999). Remarks around 50 lines of Matlab: short finite element implementation. Numerical algorithms, 20(2-3), 117-137.
- ^ Ferreira, A. J. (2008). MATLAB codes for finite element analysis: solids and structures (Vol. 157). Springer Science & Business Media.
関連文献
[編集]和書
[編集]- 大地羊三:「有限要素法とその応用」、森北出版(数学ライブラリー 38)(1975年2月20日)。
- G. ストラング、G.J.フィックス、三好哲彦・藤井宏(共訳):「有限要素法の理論」、培風館(1976年4月30日)。
- 加川幸雄:「電気・電子のための 有限要素法入門」、オーム社(1977年5月30日)。
- A.R.Mitchell、R.Wait、森正武(訳):「偏微分方程式の有限要素法」、科学技術出版(1978年5月10日)。
- K.J.Bathe、E.L.Wilson、菊地文雄(訳):「有限要素法の数値計算」、科学技術出版社(1979年9月15日)。
- J.T. オーデン、山田嘉昭(訳):「非線形連続体の有限要素法1」、培風館(1979年11月10日)。
- J.T. オーデン、山田嘉昭(訳):「非線形連続体の有限要素法2」、培風館(1980年1月20日)。
- 鷲津久一郎、宮本博、山田嘉昭、山本善之、川井忠彦(共編)「有限要素法ハンドブック I基礎編」、培風館、ISBN 4-563-03169-0(1981年9月25日)。
- 加川幸雄:「電気・電子のための 有限要素法の実際」、オーム社、ISBN 4-274-02923-9(1982年5月30日)。
- 中田高義、高橋則雄:「電気工学の有限要素法」、森北出版(1982年7月15日)。
- 水木久夫、原平八郎:「有限要素法 理論編」、森北出版、ISBN 4-627-07180-9 (1983年1月25日)。
- 鷲津久一郎、宮本博、山田嘉昭、山本善之、川井忠彦(共編)「有限要素法ハンドブック II応用編」、培風館、ISBN 4-563-03180-1(1983年1月25日)。
- 加川幸雄:「閉領域問題のための 有限/境界要素法」、サイエンス社(1983年7月10日)。
- 森正武:「有限要素法とその応用」、岩波書店(応用数学叢書)(1983年9月9日)。
- C.C.ツィエンキーヴィッツ、𠮷識雅夫(訳)、山田嘉昭(訳):「マトリックス有限要素法(三訂版)」、培風館、ISBN 4-563-03168-2(1984年9月30日)。
- スハス V. パタンカー、水谷幸夫・香月正司(共訳):「コンピュータによる熱移動と流れの数値解析」、森北出版、ISBN 4-627-91190-4(1985年2月27日)。
- 川原睦人:「有限要素法流体解析」、日科技連、ISBN 4-8171-6011-X(1985年3月8日)。
- 鷲津久一郎、池川昌弘:「有限要素法」、岩波書店、ISBN: 4-00-006065-1 (1987年6月24日)。
- 小柴正則:「光・波動のための有限要素法の基礎」、ISBN 4-627-91350-8(1990年10月31日)。
- 久田俊明:「非線形有限要素法のためのテンソル解析の基礎」、丸善、ISBN 978-4621045817(1992年11月1日)。
- 矢川元基、半谷裕彦:「有限要素法の基礎」、朝倉書店、ISBN 4-254-23079-6(1994年6月15日)。
- 日本塑性加工学会 (編):「非線形有限要素法-線形弾性解析から塑性加工解析まで」 、コロナ社、ISBN 978-4339043181(1994年12月1日)。
- 久田俊明、野口裕久:「非線形有限要素法の基礎と応用」、丸善、ISBN 978-4621041253(1996年1月1日)。
- 棚橋隆彦:「流れの有限要素法解析 I」、朝倉書店、ISBN 4-254-11406-0(1997年9月10日)。
- 棚橋隆彦:「流れの有限要素法解析 II」、朝倉書店、ISBN 4-254-11407-9(1997年9月10日)。
- 三好俊郎、白鳥正樹、坂田信二:「有限要素法解析」、朝倉書店、ISBN 4-254-11405-2(1998年9月1日)。
- 矢川元基、塩谷隆二:「超並列有限要素解析」、朝倉書店、ISBN 4-254-23662-X(1998年10月20日)。
- 久田俊明:「非線形有限要素法のためのテンソル解析の基礎」、丸善、ISBN 978-4621045817(1999年9月1日)。
- 矢川元基、青山裕司: 「有限要素法固有値解析」、森北出版、ISBN 4-627-91761-9(2001年9月28日)。
- 日本計算工学会(編)、手塚明、土田英二:「アダプティブ有限要素法」、丸善(2003年8月30日)。
- 矢川元基、吉村忍:「計算固体力学」、岩波書店(シリーズ現代工学入門)、ISBN 4-00-006942-X(2005年7月7日)。
- 福森栄次:「よくわかる有限要素法」、オーム社、ISBN 4-274-06628-2(2005年11月15日)。
- 林正:「ハイアラーキ有限要素法―大型要素による高精度解析法」、技報堂出版、ISBN 978-4765517119 (2006年12月)。
- 川面恵司、渡邉隆之、岡本紀明:「有限要素法のモデル化技術と応用解析」、養賢堂、ISBN 978-4-8425-0422-3(2007年6月27日)。
- 邵長城:「基本からわかる有限要素法」、森北出版、ISBN 978-4-627-91991-4(2008年10月1日)。
- Jacob Fish, Ted Belytschko: 「有限要素法」、丸善、ISBN 978-4-621-07996-6(2008年12月25日)。
- 浦川肇:「ラプラシアンの幾何と有限要素法」、朝倉書店(朝倉数学大系3)、ISBN 978-4-254-11823-0(2009年10月25日)。
- 板根政男:「例題で学ぶ Marc有限要素法解析入門」、丸善、ISBN 978-4-621-08362-8(2011年3月30日)。
- EA de Souza Neto、D. Peric, DRJ Owen:「非線形有限要素法-弾塑性解析の理論と実践」、森北出版、ISBN 978-4-627-92021-7 (2012年6月30日)。
- 竹内則雄、樫山和男、寺田賢二郎:「計算力学(第2版)-有限要素法の基礎」、森北出版、ISBN 978-4-627-91802-3(2012年12月5日)。
- 藤井文夫、田中真人、佐藤繊美:「Fortran90/95による有限要素法プログラミング」、丸善出版、ISBN 978-4-621-08784-8(2014年1月20日)。
- 大塚厚二、高石武史:「有限要素法で学ぶ現象と数理:FreeFEM++ 数理思考プログラミング」、共立出版、ISBN 978-4-320-01953-9 (2014年2月15日)。
- 森正武:「有限要素法とその応用」、岩波書店、ISBN 978-4007303685(2016年2月10日)。
- Javier Bonet, Richard D. Wood: 「非線形有限要素法のための連続体力学(第2版)」、森北出版、ISBN 978-4-627-67512-4 (2017年6月21日)。
- 日本計算工学会(編):「第3版 有限要素法による流れのシミュレーション」、丸善出版、ISBN 978-4-621-30183-8(2017年7月30日)。
- 久田俊明、野口裕久:「復刊 非線形有限要素法の基礎と応用」、丸善、ISBN 978-4621305676(2020年10月26日)。
洋書
[編集]- Susanne C. Brenner , L. Ridgway Scott: "The Mathematical Theory of Finite Element Methods", Springer-Verlag New York, ISBN 978-0-387-75933-3 (2008).
- Claes Johnson、Mathematics:"Numerical Solution of Partial Differential Equations by the Finite Element Method"、Dover Publications、ISBN 978-0486469003 (2009年1月15日)。
- Mats G. Larson, Fredrik Bengzon: "The Finite Element Method: Theory, Implementation, and Applications"、Springer、ISBN 978-3642332869 (2013年1月12日).
- Douglas N. Arnold: "Finite Element Exterior Calculus", SIAM (2018), ISBN 978-1-61197-553-6, doi:10.1137/1.9781611975543.
- Leszek F. Demkowicz: "Mathematical Theory of Finite Elements", SIAM, ISBN 978-1-61197-772-1 (2024).
外部リンク
[編集]- Flux2D ソフトウェア
- 日本大百科全書(ニッポニカ)『有限要素法』 - コトバンク
- Analysis3D ソフトウェア
- FEM と幾何学の機械学習
- 非線形有限要素法特論(UTokyo OpenCourseWare、2005年度開講
- 永井学志・橋本一輝:「有限要素計算における全体剛性行列の作成法−疎行列データ構造の視点から−」 その1 その2 その3
実装
[編集]- FEniCS - GitHub
- FEniCS Project (@fenicsproject) - X(旧Twitter)
- femhub - GitHub
- MFEM