• 締切済み

FEM（有限要素法）の問題の解答を教えてください

2次元空間直交座標系ポテンシャルAについて ((∂^2)A)/(∂x^2)+((∂^2)A)/(∂y^2)=0 が成立する 0≦x≦1,0≦x≦1 A(0,0)=1,A(x=1)=A(y=1)=0 x=0,y=0の境界で (∂A)/(∂h)=0 この時 X(1/2,1/2),Y(3/4,1/4) におけるAをFEMを用いて求めよ

みんなの回答

  • hashioogi
  • ベストアンサー率25% (102/404)
回答No.1

私の認識が間違っていたら申し訳ないのですが、ＦＥＭというのは微分方程式を解析的に解けない場合に計算機プログラムを使用して微分方程式の近似解を求める手法ではないかと思うのです。 ですからプログラムを作成するか、既成のプログラムを準備する必要があるように思いますし、プログラムを作成するにしてもここで披露できるほど簡単ではないだろうし。 今回の問題は境界条件がそれほど複雑ではないようですからＦＥＭなんていう牛刀を持ち出さなくても差分方程式だとか別の小刀的な方法が使用できるのでは？ もしかしたらそういう手法を含めてＦＥＭと言っているのかな？

  1. カテゴリ
  2. 学問・教育
  3. 数学・算数

関連するQ&A

  • FEM有限要素法について

    物理シミュレーション手法の１つのFEMについてお尋ねします。 ※いつもは差分法を用いて計算をしている人間からみた漠然とした感覚での質問です。 １．FEMの基本的な考え方は内挿関数（いろいろ名称がありますが）を使って節点上で定義された未知数を領域内部に連続的に展開し、方程式（重み関数をかけて）を領域内部で積分するという展開になると思います（がラーキン法）。このとき部分積分を用いて弱形式に変換したりするわけですが、この部分積分は紙と鉛筆を使って数学理論的に実施するものなのでしょうか。もともと解く方程式の複雑さによってはなかなか難しいのではないかと思いますが。どのような場合でもできるのでしょうか。 （差分だと３回微分でも何とかやることができそうです。もちろん困難もありますが見通しはつきやすい。） ２．自然境界条件（境界での微分値が指定されている）ところの説明でも部分積分が出てきます。自然境界条件の設定の考え方のコツを教えて頂きたいのですが。本を読むとポテンシャル理論を使う実例が示されていたりしており、ポテンシャル問題でない場合はどうするのかなと思います。 基本境界条件（境界での値が指定されている）についてはマトリックスの一部に000010000など使うと問題なく指定できます。これはすごくわかりやすいのですが。 ３．FEMはラプラス方程式、ポアソン方程式、などよく知られた割りと形式的に簡単に見える方程式の説明が多いのですが、微分の回数が増えるとか自分で式を作ったような場合実行するのが難しいものでしょうか。前述しましたが、３回微分の場合、１回積分して２回微分になるから内挿関数が直線的だと２回微分がゼロとなりますから使えないということになりそうです。差分とはだいぶ違うな思います。解きやすい問題は行列にしたら差分と同じということになるようですが。 以上、よろしくお願いします。漠然とした質問と言いながら長くなってしまいました。

  • 三次元ユークリッド空間上の直線の方程式は？

    三次元ユークリッド空間上で，直交座標を x, y, z とする時， 任意の平面は，a, b, c, d を実数として（abc ≠ 0）， ax ＋ by ＋ cz ＋ d ＝ 0 で表されます． では，三次元ユークリッド空間上の任意の或る一つの直線の方程式は， 直交座標を x, y, z とする時，一般的に，どの様に表現されるのでしょうか？ どなたか，教えて下さい．

  • 直交変換と回転は同じものなの？

    座標の回転が直交変換なのは任意の回転がｘ軸の回転とｙ軸の回転とｚ軸の回転の組み合わせであることから理解できます しかしその逆が分からないのです つまり直交変換は座標の回転なのかどうかです ３次元直交座標Ａと３次元直交座標Ｂがある 空間に点Ｐがある ＰのＡによる座標を（ｘ，ｙ，ｚ）＝ａ＾Ｔとし ＰのＢによる座標を（Ｘ，Ｙ，Ｚ）＝ｂ＾Ｔとする そこで質問します 「Ｕ＾Ｔ・Ｕ＝Ｅかつ｜Ｕ｜＝１である３次正方実行列Ｕがあり任意のＰについてａ＝Ｕ・ｂならばＢはＡを原点を中心に回転したものである」 は正しいのですか？ 正しければどうしてなのですか？ 正しくなければどうしてなのですか？ よろしくお願いします

  • ラプラス方程式の解析解

    電磁気学を勉強しているのですが，分からないことがあるので質問させてもらいます． 静電場内にある電荷が作る電位分布を示す方程式としてラプラス方程式(∇^2*V=0)があると思います． ラプラス方程式とポアソン方程式の違いまでは理解できていると思います． 2次元のラプラス方程式は以下の式を変数分離法を用いて解くことで，直交座標系や球面座標系として考えることで，解析解が得られると理解しています． (ここまではたどり着くことが出来ました) (∂^2/∂^2x)V(x,y)+(∂^2/∂^2y)V(x,y)=0 分からないのは，ここから実際の電位分布を求める方法です． 具体的には，xy平面上の原点にポテンシャルV0がある場合，このV0による電位分布を求めることが出来ません． 直交座標系で考えると一般解は，A,B,C,D,kを定数として，次のようになると思います． V(x,y)=(A*exp(kx)+B*exp(-kx))*(Csinky+Dcosky) 境界条件から未知定数を求めたいのですが，うまくいきません・・・． 原点にポテンシャルがあるので，x→∞でV→0，y→∞でV→0，x=0,y=0でV=V0が境界条件になると思ったのですが，y→∞で(Csinky+Dcosky)は0に収束しません． 境界条件の設定が間違っているのでしょうか？ 数値解では原点にポテンシャルを設定している解説は見つけられたのですが，解析解では資料がなく，どうすればいいか困っています． すみませんが，教えてください．

  • 4次元空間問題

    4次元ベクトル空間(変数はxyzu) x+y+z+u=1 において、 この式を満たす空間上にあり、この空間と直交し、お互いに直交する3つのベクトル空間を求めて下さい。

  • 三次元空間の３点のなす角度

    三次元空間の３点のなす角度を知る公式が知り無たいのです。直交座標で、３点の x, y, z 座標値はわかっているものとします。 自分でいろいろ考えたのですが、かなりややこしくなってしまいこれは公式を見つけないとだめだなと思いました。よろしくお願いします。

  • 3次元空間上の2点を結ぶ線分の中点を知りたい

    3次元空間上の点A(x1, y1, z1)と点B(x2, y2 z2)を結んで出来る線分の中点を知りたいのですが、 完全な文系出身であまり数学に詳しくないため、公式の見方がよくわかりません。 Wikipediaの中点のページにあるn次元ユークリッド空間上の中点の公式がそれのようですが、 「n 次元ユークリッド空間上の2点 A, B を直交座標系であらわし、それぞれを (a0, ..., an－1), (b0, ..., bn－1) とすると」 の時点ですでに理解できないので、単純な公式で教えて下さると助かります。

  • 線形代数の問題です

    3次元空間R^3の点をxyz座標を用いて縦ベクトル (ベクトルx)=(x,y,z) ∈R^3によって表示する 以下ではR^3内の二つの平面、α：z=0、 β：x+y+z=0への直交射影を考える。 以下の手順で平面βへの直交射影が次式で与えられることを示せ (1)点(ベクトルx)を通る直線：(ベクトルx)+t(ベクトルu) (t∈R) が平面βと直交したとする。このような(ベクトルu)を求めよ (2)(ベクトルx)から平面βへの垂線の足(ベクトルb)を求めよ この２問がどうしてもわかりません。 どなたかご回答よろしくお願いします

  • 直交行列とノルムについて

    n次元Euclid空間のベクトルxと任意の直交行列Aに対して， y≠Ax ⇒ |y|≠|Ax|　（|x|はEuclidノルム） は直接示せますか？ 一般的に「x≠y⇒|x|≠|y|」は言えませんが，ここではAが直交行列で大きさを変えない変換ということから|y|≠|Ax|と言ってしまったのですがこれは言えないのでしょうか？ どんな直交行列を持ってきてもy≠Axとなるということは大きさが異なるということではないのですか？ よろしくお願いします．

  • 直交補空間に関する問題です。

    直交補空間に関する問題です。 当方あまり数学を厳密に理解をしておりません。 以下の問題が進まずに困っております。 R4内で a1=(4 3 2 1), a2=(1 3 5 7), a3=(1 2 3 4)によって生成される部分空間をVとして、通常のユークリッド内積に関する<a1>の直交補空間をWとするとき、 x*a1 + y*a2 ∈ W となるための実数x, yに対する条件を求めよ。 -------------------------- 直交補空間自体あまり理解できていないのでそこら辺からやさしめに教えていただけると助かります。 よろしくお願いします。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する