- ベストアンサー
有限要素法解析について
有限要素法では、線膨張係数をどのような計算によって考慮しているのですか？ 線膨張係数に温度変化をかけて、歪みが求められると思うのですが、それから ε＝〔B〕u という関係から変位をもとめて、仮想仕事の原理より求めた剛性マトリックス〔K〕をつかって F＝〔K〕u の関係が全体を満たすように計算しているのでしょうか？ 間違いのご指摘、正答の提示、なんでもかまいません。 どなたかご返答をお願い致します。
- wandeful12
- お礼率100% (2/2)
- 物理学
- 回答数1
- ありがとう数1
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
「初期ひずみ」という概念を入れているらしいですね． 私も使ったことは無いので，本をもう少しよまないと 確かなことはわかりません． O.C.ツィンキービッツの本です． 温度による膨張は，単一材料の場合，拘束されていなければ，それによって，応力は生じません． それをどのように表現しているのかが興味のあるところですね． 単なる強制変位ではないと言うことです． それでは．
関連するQ&A
- 有限要素法 FEMについて
有限要素の変位法についての質問です。 有限要素の変位法で応力を求める場合 どのような順序で計算されるのでしょうか？ このような順序であってますか？よろしくおねがいします。 (1)各メッシュの剛性を求める (2)各メッシュの剛性をモデル全体と重なり合わせて、モデル全体の剛性を求める (3)拘束条件及び荷重条件を入力して、モデル全体の変位を求める (4)モデル全体の変位量を各要素に展開して、各節点の変位を求める (5)節点の変位量を使い、要素の歪量、応力値を求める。
- 締切済み
- 物理学
- 有限要素法の多点拘束の要素段階での組み込み
FEMプログラムの中に多点拘束の機能を組み込みたいのですが、全体マトリクスの段階で組み込むと、従属関係を表すマトリクスが大きくなるので、要素マトリクスの段階で組み込む方法があれば、教えてくれないでしょうか。なお私の理解している全体マトリクスの場合は次です。 (1)元の剛性方程式：[K]{d}={F} (2)従属関係：{d}=[C]{d0}、ここに[C]は従属関係を表すマトリクス、 {d0}は主（master)変位 ここで(2)を(1)に代入し、両辺に前から[C]の転値を乗じると (3) [C]T[K][C]{d0}=[C]T{F} （３）式を解いて主変位を求め、さらに（２）から従変位を求める。
- ベストアンサー
- 物理学
- 有限要素法の要素コネクティビティの影響
平面ひずみと仮定して三角形定ひずみ要素を用いて片持ち梁の解析を行いました。 解析条件は2要素4節点です。図としては下のURLの節点1と2を完全に固定して 節点4に下向きの力をかけるというものです。 端的に言いますと正方形を三角形二つにわけた状態です。 http://kentiku-kouzou.jp/yugenyoso-gouseimatorikusu.html そこで質問なのですが、要素の切り方を図のように 要素1(1,2,3),要素2(1,3,4) としたときと 要素1(1,2,4),要素2(2,3,4) としたときでは、でてくる変位の解が異なりました。 手計算でも行ってみたのでプログラム自体の誤りではないと思います。 同じ要素数、自由度数なのに何故解が異なってしまうのでしょうか？ 力のかかる部分を二つの要素で共有してるかしてないかの問題なんでしょうか？ それともそれぞれの全体剛性マトリクスが異なることが原因なのでしょうか？ よろしくお願いします。
- ベストアンサー
- 物理学
- 有限要素法において応力分布図の作成
２次元有限要素法で平面応力の解析をしています。４節点アイソパラメトリック要素を用いています。接点変位は計算できました。接点応力、接点ひずみも計算しました。これらからどのような方法で応力分布やひずみ分布を計算し図を描けばよいのでしょうか？市販のソフトは使いたくないのですが。
- 締切済み
- 物理学
- 骨組の解析について
フレームモデルを魚骨形や多質点系モデルへと置換して骨組を解析する際、わからないことがあります。 骨組を設計し、各部材に部材の応力ーひずみ関係を設定して、部材レベルで応答解析をする場合はこのまま解析を行えば良いのですが、より簡易なモデルへと置換した際は各層で集約された柱、梁に応力ーひずみ関係を設定する事になります。 このような集約された骨組を解析するとき、集約柱や集約梁に与える応力ひずみ関係は部材レベルで応答解析をする際に設定したものをそのまま使ってよいのでしょうか？ 元の骨組の部材剛性、耐力を元に集約柱や集約梁の剛性、耐力を算出し、それに見合うような断面を持つ仮想部材をつくり、それに応力ーひずみ関係を与えると思うのですが、この応力ーひずみ関係は層を代表する応力ーひずみ関係になるはずです。 層の応力ーひずみ関係と言われてもいまいちピンと来ないもので、剛性や耐力は集約するために計算を行うのに、応力ーひずみ関係は何も触らずにそのまま使用して良い者なのか、どうもよくわかりません。
- 締切済み
- 建築・土木・環境工学
- 有限要素法による解析～熱変形と硬化収縮～
光で硬化する材料を扱っています。 硬化するときに材料が収縮し、変形してしまいます。 そのため硬化時の収縮による変形を 有限要素法で解析したいと思っております。 収縮だけなら収縮率を熱膨張係数に置き換えて、 解析するという策もあるようなのですが （何とか合わせ込んでるといった感じ） 今回やろうと考えているのが、 加熱環境下での硬化を視野に入れた計算したいと考えています。 加熱下でも硬化は進行するとともに、熱による変形も同時に発生します。 また、反応性の材料ですので反応熱による温度変化もあります。 熱パラメータと硬化パラメータを同時に 解析できないかと考えあぐねております。 ちなみに現在はＡＢＡＱＵＳを使っております。 （最近使い始めたばかりの超初心者ですが・・・・） ちゃんと意味を理解していただけるか分かりませんが これらが実現できそうなソルバーなどの情報がありましたら お教えいただきたいと思います。 （ＡＢＡＱＵＳでこういった解析ってできるんですか？）
- ベストアンサー
- 物理学
- 質量行列と剛性行列について（ラーメン構造）
失礼します。 今、図のような条件のもとでプログラムを作っています。 これは、ラーメン構造ですので、非常に説明が長くなっていますし、 専門的なのでもしよければご回答お願いします。非常に困っています。 まず、私のプログラムでは、非減衰のものを作ります（本当はレーリー減衰で決めるけれど、まずは剛性、質量マトリックスを作ってから） 自分は下記の式を用いて、振動を解析しています。 Mu''+Ku=-mg(地震動） M：質量マトリックス Ｋ：剛性マトリックス u:変位 u'':加速度（変位２回微分） -mg(これはどこからか資料をとってきてこの中にいれます）つまり、gは加速度です。 まずは、要素マトリックスですが、剛性マトリックスは以下の通りにしています。 各節点には、3つの力があります（軸力：Ｎ、垂直力：Ｖ、曲げモーメント：Ｍ） よって、6行6列の行列になります。 k=[EA/L 0 0 -EA/L 0 0;0 12EI/L^3 6EI/L^3 0 -12EI/L^3 6EI/L^2; .....] ここに E:ヤング係数、A：部材断面積 L:長さ これは調べていただければわかると思います。 問題は、質量マトリックスのほうです。 調べたところ、以下の式が出てきました。 m=ρAL/420*[210 0 0 0;0 L^2 0 0; 0 0 210 0; 0 0 0 L^2] と4行4列の式が出てきました。（対角行列） ρ；密度 Ａ、Ｌ：上記と同じ ここで最初の質問です。 (1)上式は、おそらく私の考えでは、軸力あるいは垂直力、と曲げモーメントに対しての マトリックスであるから、4行4列だと思います。 もしここに、欠けているであろう、軸力、あるいは垂直力をつけるとするならば下の式であっていますか？ m=ρAL/420*[210 0 0 0 0 0;0 210 0 0 0 0; 0 0 L^2 0 0 0; 0 0 0 210 0 0;0 0 0 0 210 0; 0 0 0 0 0 L^2] の6行6列の対角行列の式であっていますか？ また、この次に、変換行列を使って、この図でいえば、 (5)、(6)の要素を垂直に戻すか、 (1)、(2)、(3)、(4)を水平に戻すかです。 このときに、mは変換しなくてよいと思うので、残りはｋだけです。 ここで２つ目の質問です。 (2)これは、普通に静的なときにやる座標変換でも大丈夫でしょうか？（今は動的です） なお、そのあとは、座標変換したものを、要素ごとに重ね合わせの原理を用いて、全体剛性マトリックスを作っていきます。もちろん、全体質量マトリックスも作っていきます。 この二つさえできれば、減衰定数は作れます（レーリー減衰）が、今はこの段階です。 専門的になりすぎてなにとぞご不便をおかけしますがよろしくお願いします。 非常に困っています。 失礼します。
- 締切済み
- 物理学
- 質量マトリックスの導出
今、弦の振動を有限要素法で解く勉強をしているのですが、 途中で運動エネルギーから質量マトリックスを導くところがあるのですが。。。 Ｔ=∫ρ・ｕ(t)＾Ｔ・Ｌ(x)・Ｌ(x)＾Ｔ・ｕ(t) ｄx =ｕ(t)＾Ｔ・∫ρ・Ｌ・Ｌ＾Ｔ ｄｔ・ｕ(t) =ｕ(t)＾Ｔ・[M]・ｕ(t) [M]＝∫ρ・Ｌ・Ｌ＾Ｔ ｄｔ T；運動エネルギー L(x)；形状関数 u(t)；変位ベクトル ○^T；○の転置行列 [M]；質量マトリックス という部分なのですが、式の変形自体はわかるのですが、 なぜ運動エネルギーから質量マトリックスが導出されるのか、 いまいちわかりません。 運動方程式から、仮想仕事の原理やハミルトンの原理を用いて、 表されるというのを聞いたのですが、これらと運動エネルギーの式を どうつなげるのか。。。どなたか教えてください。 よろしくお願いします！！！ わかりにくい表現が多くてすみません・・・
- ベストアンサー
- 物理学
- 自由端の境界条件（二次元有限要素法）
自由端の境界条件（二次元有限要素法）を教えてください。 二次元有限要素法を使って、薄いテープ状の線材の断面（厚みー長さ平面）の・熱応力・を計算したいのですがうまくいきません。 この線材は厚み方向に３層構造になっているんですが、 知りたいのは 中心の層と上下の層の界面に働くせん断応力・ひずみ 中心の層に働く長手方向の応力・ひずみ です。 最初は剛体壁に両端を固定された状態で線材全体を一様に温度上昇させるモデルで 解析していたのですが、両端を固定（両端の接点変位＝０）しているためか 長手方向応力はありえない程大きい値になるし、 せん断ひずみや応力はほとんど０になってしまいます。 そこで境界条件を自由端にしたいのですが、 自由端の境界条件を教えていただけないでしょうか？ よろしくお願いいたします。 また、このようなテープ状の材料の強度を評価するためには何を調べればいいんでしょうか？ 上ではせん断応力～と書きましたが自分で半信半疑です。曲げモーメント（？）とかも計算するべきでしょうか？ 私、専門は電気工学で突然降ってきた仕事にとても困っています。 後、このようなテープ状の材料を解析するのに適した要素をしっておられましたらアドバイスお願い致します。 ちなみにサイズが厚さが数100umで、幅が数mmオーダー、長さが10cm程度です。
- ベストアンサー
- 物理学
- 仮想仕事の原理で出てくるδとは
エネルギー原理や仮想仕事の原理で出てくるδ（デルタ）は、どのように理解すれば良いですか？ いままでは、微小な仮想変位という理解で δU とあってもこれは演算子ではなく、 二つでひとつの記号のように理解していたのですが 例えば http://www.md.ams.eng.osaka-u.ac.jp/~nakatani/Lectures/Fundamentals_of_Solid_Mechanics/TEXT/HTML/node58.html の式（331）などで、私にとっては急に、 くっついたり離れたりし始めてなにこれ？となってしまいました。 有限要素法を理解するための通過点で、深入りする必要はないのですが 漠然とで良いので、薄っすら納得したいです。 よろしくお願いします。
- ベストアンサー
- 数学・算数
お礼
返事が遅くなってしまってすみません。返答ありがとうございます。温度変化から温度荷重（外力）として接点にかけることで、仮想仕事の原理から有限要素法を用いているのだということまでわかりました。短い説明でわかりにくいかと思いますが・・。また、この見解が間違いということがありましたら、よろしくお願いします。何分独学なので。