★有限要素法(Finite Element Method)★ (392レス)
★有限要素法(Finite Element Method)★ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/
上
下
前
次
1-
新
通常表示
512バイト分割
レス栞
200: 名無しさん@5周年 [] 2012/08/16(木) 00:49:01.29 >>197 その差は、特に長手方向に粗いのが原因っぽい メッシュ細かくしてきゃ、どっちも値大きくなって、双方一緒っぽい値になる筈 テトラ一次は細かくして行っても、2次と全然違う値になったりするけど http://kizuna.5ch.net/test/read.cgi/sim/1284635327/200
201: 名無しさん@5周年 [] 2012/08/16(木) 09:53:53.82 >>200 まぁそうなんだけどねw つまり、二次要素は、メッシュ切りがテキトーでもそれなりの精度が出る。 一次はメッシュの細かさが、モロ、精度に出るので、メッシュ密度とか形状に非常に神経を使わなきゃならんので 実務的には向いていないという考えです。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/201
202: 名無しさん@5周年 [] 2012/08/16(木) 12:39:52.17 >>201 テトラの2次と1次は 3割4割違ってくるので 比べるとマシだけど ちょっと差があるのは気になるね ウーム http://kizuna.5ch.net/test/read.cgi/sim/1284635327/202
203: 名無しさん@5周年 [] 2012/08/16(木) 20:22:52.36 >>197 あ、書き忘れていた、これの理論値は -0.75mmです。 つまり 理論値 -0.75 1次6面体(8節点)要素 -0.72 2次6面体((20節点)要素 -0.63 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/203
204: 名無しさん@5周年 [] 2012/08/16(木) 20:24:58.37 >>203 う、逆だった 理論値 -0.75 2次6面体((20節点)要素 -0.72 1次6面体(8節点)要素 -0.63 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/204
205: 名無しさん@5周年 [] 2012/08/18(土) 06:46:14.80 倍メッシュ切ればいいじゃない http://kizuna.5ch.net/test/read.cgi/sim/1284635327/205
206: 名無しさん@5周年 [] 2012/08/18(土) 11:24:42.70 >>205 その「倍」っていうのはどこから来る数値なのか? たとえばこういう単純ハリの場合、一次要素で解析するには、 縦横高さで何等分、アスペクト比をいくつにすればいい、 っていう目安ってありますか? http://kizuna.5ch.net/test/read.cgi/sim/1284635327/206
207: 名無しさん@5周年 [sage] 2012/08/18(土) 11:32:30.25 「倍」というのが辺の長さを指しているのであれば、 3次元だから2の3乗で8倍のメッシュを必要とするという理解で良い? メッシュ数が「倍」であれば、辺の長さを2の立方根で割れば良い? http://kizuna.5ch.net/test/read.cgi/sim/1284635327/207
208: 名無しさん@5周年 [sage] 2012/08/18(土) 13:16:14.14 正六面体に近い形でメッシュを切れば良いと思うの http://kizuna.5ch.net/test/read.cgi/sim/1284635327/208
209: 183 [sage] 2012/08/18(土) 17:37:04.75 皆遊んでいる様なので質問を・・・ 固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか? ググっているのですがそれらしい文献を見つけられずに困っています。 詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください 固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを 逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー+QR法+逆反復法で 出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します ここらもエスパーレスしてくれると有難いです テストデータは次の感じです ・方持ち梁 ・六面体一次要素(要素数1 と5これが主 と 625これは時間が足りなくて計算できない) ・拘束条件無し(なんかこれがいけないっぽい) http://kizuna.5ch.net/test/read.cgi/sim/1284635327/209
210: 名無しさん@5周年 [] 2012/08/18(土) 18:22:03.30 >>209 >・拘束条件無し(なんかこれがいけないっぽい) その通り。 固定条件の入れ方って、単に、固定する節点の変位を0にするだけじゃないの? http://www.jsce.or.jp/committee/struct/civil_edu_home/prog/anketo/H19_femdoc02.pdf http://kizuna.5ch.net/test/read.cgi/sim/1284635327/210
211: 名無しさん@5周年 [sage] 2012/08/18(土) 18:25:46.19 >>208 その通りだと思います。 ただ、ちょっと複雑な形状をきると、アスペクト比1を目指してメッシュを切ると 要素数が非常に多くなることがあるので、1次要素を使ったほうが計算時間が かかってしまう事はあるようです。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/211
212: 名無しさん@5周年 [] 2012/08/18(土) 18:48:00.38 >>210 スライドできる状態=摩擦ゼロ状態いうのは、厳密に起こり得ないので、 条件に対し解を求めるいう点では合ってるけど スライドできるいう条件設定が、果たして現実的か? そこはかなり難しい 例えば土木や建築の分野だと 地盤をモデル化したり、剛性体を挟んだり 色々しないと、単品+拘束やすべり条件設定では、合わない http://kizuna.5ch.net/test/read.cgi/sim/1284635327/212
213: 183 [sage] 2012/08/18(土) 19:01:45.09 >>210 変形の計算(静計算って言うのかな?)はその通りだと思うのだけど、固有値、固有モード計算と同じなんですか? 変形計算の時って [K]{U}={F}で固有値計算の時は[K]{U}=λ[M]{U}ですよね。これをどう節点変位ゼロと 絡めたら良いのか分からないです。頭が悪くて申し訳ない http://kizuna.5ch.net/test/read.cgi/sim/1284635327/213
214: 名無しさん@5周年 [] 2012/08/18(土) 19:49:21.46 >>213 固有振動の場合は固定点がなくても解が出るはず。 収束計算させるんじゃなくて、ガウス吐き出し法とかで解けば? http://kizuna.5ch.net/test/read.cgi/sim/1284635327/214
215: 183 [sage] 2012/08/18(土) 21:13:36.78 >>214 境界条件がある固有値、固有ベクトル計算も知りたいので宜しくお願いします。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/215
216: 名無しさん@5周年 [] 2012/08/19(日) 21:51:17.23 >215 固有値=フリーでの計算ばかり 固定など条件設定を伴ったものは見た事がない気もします できないのかも知れません http://kizuna.5ch.net/test/read.cgi/sim/1284635327/216
217: 名無しさん@5周年 [] 2012/08/19(日) 22:09:42.27 >>216 出来ないってこたーないだろ。 AnsysだろうがNastranだろうが、出来ないCAEを見たことねぇぞ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/217
218: 名無しさん@5周年 [] 2012/08/19(日) 22:44:51.38 そもそも何の研究をしようとしてるのか わけわからん。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/218
219: 名無しさん@5周年 [sage] 2012/08/19(日) 22:46:37.18 ロケットサイエンス http://kizuna.5ch.net/test/read.cgi/sim/1284635327/219
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 173 レスあります
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.944s*