1994 Fiscal Year Annual Research Report
| Project/Area Number |
06650073
|
|---|
| Research Category |
Grant-in-Aid for General Scientific Research (C)
|
| Research Institution | The University of Tokyo |
|---|
Principal Investigator |
速水 謙 東京大学, 工学部, 助教授 (20251358)
|
|---|
| Keywords | 数値解析 / 境界要素法 / 偏微分方程式 / 数値積分 / 連立一時方程式 / 反復解法 / 静弾性問題 |
|---|
| Research Abstract |
境界要素法は、境界のみを離散化すればよい、開放領域が扱いやすい等の利点から、近年工学で有力な偏微分方程式の数値シミュレーション手法として定着してきた。しかしながら、解決すべき数値解析上の問題も多い。その中で今年度は、1.数値成分の高精度化、2.連立一次方程式解法の高速化、について検討した。
1.数値積分の高精度化
境界要素法で薄い構造や、境界付近での場を計算する際に生じる近特異積分(有界ではあるが、急激に変化する関数の積分)に対して、筆者らは近特異性を適切な変数変換によって緩和してから積分する手法を開発して来た。
しかし、従来の手法はGauss-Legendre法に基づいていたため、自動積分に向いていなかった。そこで今年度は、この変数変換手法に精度保証つきの自動積分の昨日を持たせるため、(動径方向の積分に対して)台形公式に基づいた新たな変数変換(logL_2-DE変換)を開発し、その有効性を数値実験及び理論誤差解析により検証した。
2.連立一次方程式解法の高速化
境界要素法では、密行列を扱うので直接法を用いるのが常套手段であるが、要素数nが大きくなると計算量(n^3)・記憶量(n^2)が莫大になり、実用化において問題になっている。そこで、従来の直接法の代わりに反復法を用いる事に着目し、更に高速化するため、遠い要素間の影響に近似的に取り扱うクラスタリング法を検討した。この手法はポテンシャル問題に対しては知られているが、これを2次元静弾性問題に適用する方法の開発に成功した。これにより、計算量・記憶量は共にn log nで済み、従来不可能だったワークステーション上での3000要素の計算に成功した。
|
Research Products
(6 results)
-
[Publications] K.Hayami: "A Numerial Quadrature for Nearly Singular Boundary Element Integrals" Engineering Analysis with Boundary Elements. Vol.13,No.2. 143-154 (1994)
-
[Publications] K.Hayami: "Improvement of Quadrature for Nearly Singular Integrals in 3D-BEM" Boundary Element Method XVI. 201-210 (1994)
-
[Publications] 国廣昇: "境界要素法の変数変換型の自動数値積分法とその誤差解析" 日本応用数理学会論文誌. 第5巻,第1号. (1995)
-
[Publications] T.Washio: "Parallel Block Preconditioning Based on SSOR and MILV" Numerical Linear Algebra with Applications. Vol.1(6). 533-553 (1994)
-
[Publications] T.Washio: "Overlapped Multicolor MILV Preconditioning" SIAM Journal of Scientific Computing. Vol.16,No.3. (1995)
-
[Publications] 山田賢志: "クラスタリング法の2次元静弾性問題への適用" 境界要素法研究会、境界要素法論文集. 第17巻. 31-36 (1994)
