- ベストアンサー
はりのモード解析
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
差分法は微分方程式を近似して解く方法です。微分方程式に解析解があるとは限らないので、通常差分法や有限要素法で近似して解きます。
関連するQ&A
- 有限要素法 FEMについて
有限要素の変位法についての質問です。 有限要素の変位法で応力を求める場合 どのような順序で計算されるのでしょうか？ このような順序であってますか？よろしくおねがいします。 (1)各メッシュの剛性を求める (2)各メッシュの剛性をモデル全体と重なり合わせて、モデル全体の剛性を求める (3)拘束条件及び荷重条件を入力して、モデル全体の変位を求める (4)モデル全体の変位量を各要素に展開して、各節点の変位を求める (5)節点の変位量を使い、要素の歪量、応力値を求める。
- 締切済み
- 物理学
- FEM薄板曲げ要素の節点断面力から理論解への換算
有限要素法で長方形の板曲げ要素を自作し、周辺固定板の計算をして理論解と比較したところ、節点変位は数％の誤差で合いますが、節点の力（断面力）が合いません。そこで気がついたのですが、FEMの節点の断面力の定義は理論解のそれと違うのではないかということです。FEMの断面力は次のように計算しました。「FEM断面力=要素剛性×要素座標系の節点変位」 一方、理論解の方は「辺の単位長さあたりの力」のようで、また記号Mx等の意味も違うようです。 そこで質問ですが、FEMの節点の断面力から理論解の断面力に換算するにはどうすればよいのでしょうか。
- ベストアンサー
- 物理学
- FEM有限要素法について
物理シミュレーション手法の１つのFEMについてお尋ねします。 ※いつもは差分法を用いて計算をしている人間からみた漠然とした感覚での質問です。 １．FEMの基本的な考え方は内挿関数（いろいろ名称がありますが）を使って節点上で定義された未知数を領域内部に連続的に展開し、方程式（重み関数をかけて）を領域内部で積分するという展開になると思います（がラーキン法）。このとき部分積分を用いて弱形式に変換したりするわけですが、この部分積分は紙と鉛筆を使って数学理論的に実施するものなのでしょうか。もともと解く方程式の複雑さによってはなかなか難しいのではないかと思いますが。どのような場合でもできるのでしょうか。 （差分だと３回微分でも何とかやることができそうです。もちろん困難もありますが見通しはつきやすい。） ２．自然境界条件（境界での微分値が指定されている）ところの説明でも部分積分が出てきます。自然境界条件の設定の考え方のコツを教えて頂きたいのですが。本を読むとポテンシャル理論を使う実例が示されていたりしており、ポテンシャル問題でない場合はどうするのかなと思います。 基本境界条件（境界での値が指定されている）についてはマトリックスの一部に000010000など使うと問題なく指定できます。これはすごくわかりやすいのですが。 ３．FEMはラプラス方程式、ポアソン方程式、などよく知られた割りと形式的に簡単に見える方程式の説明が多いのですが、微分の回数が増えるとか自分で式を作ったような場合実行するのが難しいものでしょうか。前述しましたが、３回微分の場合、１回積分して２回微分になるから内挿関数が直線的だと２回微分がゼロとなりますから使えないということになりそうです。差分とはだいぶ違うな思います。解きやすい問題は行列にしたら差分と同じということになるようですが。 以上、よろしくお願いします。漠然とした質問と言いながら長くなってしまいました。
- 締切済み
- 情報工学
- 数値解析ソフトを用いた実験の授業のあと…
数値解析ソフトを用いた実験の授業のあと、次のような課題が出ました。どんな本を調べて良いのか分からないままレポート提出が近づいて来てヤバいです…(T_T)どうかお助けくださいませ… (1)計算手法として多くの流れ解析に用いられている有限要素法と有限差分法について，特徴を述べよ。 (2)流れのシュミレーションに関し、実験結果と比較する方法としてどのような方法で行われているか述べよ。 どうかよろしくお願いします。
- 締切済み
- その他（学問・教育）
- 曲がりはりの変形をたわみの基礎式で求められるか
タイトルのとおりですが、曲がりはりの変形は通常エネルギー法を使用した方が便利と習いましたが たわみの基礎式でもたわみを求めることはできるのでしょうか 例えば下記URL先のような1／4円の先端に鉛直方向の荷重が作用する場合ではいかがでしょうか https://bloodystream.hatenadiary.jp/entry/2021/04/25/011358 モーメントがxの関数ではないので、そもそも適用できないのでしょうか。
- 締切済み
- 機械設計
- 微分法・積分法は知ってるけど「差分法」って一体何ですか？
差分法というのを聞いたんですけど、一体差分法って何ですか？ 微分積分系（？）ですか？それとも全く関係の無い数学ですか？ 差分法とは何で、どういう計算法（？）か軽く教えてください。
- ベストアンサー
- 数学・算数
- 数値解析の手法（差分法）について
現在、とある2元の1階偏微分方程式（解はu,vでそれぞれ右と左に進行する波）を数値解析によって解こうと考えています。 数値解析の手段として、差分法がよく用いられると思いますが、 現在、私は、場所に関してはuは後退差分、vは前進差分を使い、 時間に関しては前進差分を使って解いています。 ネットでは場所に関しては中心差分、時間に関してはルンゲ=クッタやリープフロッグなどが 使われていることが多く、私もこの２つを用いて解いてみました。 偏微分方程式には線形項が含まれていたため、 線形問題に対する制約であるΔt/Δx << 1は最低満たすように刻み幅をいろいろ取り、 計算時間も辞さず計算機を動かしてみましたが、 ノイズが消えず、解析解に限りなく近づくには至りませんでした。。 Δt/Δx=0.0001なども試したのですが・・・ そこで、諦めて違う差分法を試し、 場所に関して、uは後退差分、vは前進差分を使い、 時間に関しては前進差分を使って見たところ、 Δt/Δx=0.01程度で解析解に近い、なかなか精度の良い数値解を得ることが出来ました。 2次の差分では上手くいかず、1次の差分だとわりかし上手くいく・・・ 精度的には中心差分やルンゲ=クッタなどの方がいいと思うのですが・・ 正直不思議でなりませんでした。。。 最初に試した差分法のコードミスかと思い、何回もコードを確認し直しましたが、 やはり解析解に近づくには至りませんでした。 こんなことってあるのでしょうか？？ 差分法でも場合によって使い分ける必要があるということでしょうか・・？ その場合分けするときの指標など、知っておられる方、教えて頂けると助かります。 問題によって
- 締切済み
- 数学・算数
- 節点応力の求め方（有限要素法）
今、有限要素法をアイソパラメトリック要素（４節点）でやっています。 で、要素ごとの応力・ひずみは、同じようにヤコビアンなどを計算して、できるんですが よく考えたら節点ごとの応力・ひずみを求める方法がわかりません。 どのように求めるのでしょうか。 参考サイトでも構いませんのでよろしくお願いします。
- ベストアンサー
- 物理学
- 有限要素法の手順とその原理に関する初歩的な質問
有限要素法という数値計算の方法は、まず方程式(A=B)があって、それをA-B=0とし、任意の重み関数ωを乗じて方程式の成立する範囲で積分（∫ω(A-B)dA=0）し、さらにそれを部分積分し、面積分（領域全体の）と線積分（境界上の）とします（弱形式）。そして、領域全体を細かく分割した有限要素内部での未知変数（重み関数も）を内挿関数とそれらの節点値（三角形だったら基本３個とか）を使って内挿して表現し、求めた式に代入してωiでくくって代数方程式を作り、それを全体マトリックスとして組み上げて行列計算を行う、というものだと思います（ガラーキン法: 汎関数と変分法を使うリッツ法と等価になる？）。式を使わず言葉で長々と書きましたが。そういう風に理解しています。（間違っているかもしれませんが）。この一連の流れの中でどうしても１つ腑に落ちないところがあります。それは部分積分のところです。この部分積分は必須なのでしょうか。部分積分をすることによって境界積分が生じるので境界条件を課するのにちょうどいいからということなのでしょうか。部分積分をしないと先に進めないというところがやや理解しずらいのですが。よろしくお願いします。
- 締切済み
- 数学・算数
お礼
回答ありがとごさいます！ 助かりました^^