★有限要素法(Finite Element Method)★ (392ﾚｽ)
★有限要素法(Finite Element Method)★ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/
上
下
前
次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
191: 名無しさん＠５周年 [] 2012/08/09(木) 01:23:40.38 固有振動数レベルだとANSYS、NASTRAN、ABAQUSどれも結果変わらないし。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/191
192: 名無しさん＠５周年 [sage] 2012/08/09(木) 01:41:20.54 >>189 六面体の2次要素はまだ実装していないのです >>190,191 すみません、レスを混同して変なレスをしています http://kizuna.5ch.net/test/read.cgi/sim/1284635327/192
193: 名無しさん＠５周年 [] 2012/08/10(金) 12:36:54.72 テトラは2次でないと全然駄目だけど ６面体は一次で十分です 固有値グチャグチャは気になるけど 固定が駄目なのかなー http://kizuna.5ch.net/test/read.cgi/sim/1284635327/193
194: 名無しさん＠５周年 [] 2012/08/11(土) 07:48:08.08 今時二次がなんていつの時代のソフトだよｗ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/194
195: 183 [sage] 2012/08/11(土) 13:13:16.14 >>193 固定か〜固定している所の剛体行列と質量行列をどのように扱ったら良いか分からないから固定条件無しでやっています まず剛体モードが綺麗に直行で出てくれないし何か変だ〜 今は固有ベクトルの計算は逆反復法らしきものでやっています (上手く収束しないんだよねと言うかほとんど収束しない、振動してたり反復の上限一万回に引っかかってたりするよ) >>194 手作りだからね・・・まずは一次要素からかなと思って作ってます (ちょこちょこ弄っていて大分汚くなってきたよw) http://kizuna.5ch.net/test/read.cgi/sim/1284635327/195
196: 名無しさん＠５周年 [] 2012/08/11(土) 13:18:46.78 仕事のためじゃなくて勉強か http://kizuna.5ch.net/test/read.cgi/sim/1284635327/196
197: 名無しさん＠５周年 [] 2012/08/15(水) 10:58:16.91 >>193 いやいやいやいや、 ６面体だろうが、なんだろうが、２次要素と１次要素じゃ精度が全然違う。 普通の構造解析で積極的に１次要素を使う理由は、全く無い。 ↓ http://www.42ch.net/UploaderSmall/source/1344995786.png http://kizuna.5ch.net/test/read.cgi/sim/1284635327/197
198: 名無しさん＠５周年 [] 2012/08/15(水) 15:57:16.77 ２次なんて言ってる時点で終わり。分からないの？ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/198
199: 名無しさん＠５周年 [] 2012/08/15(水) 18:18:06.88 >>198 わからない。 学問板なので、理論的に説明してくれ。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/199
200: 名無しさん＠５周年 [] 2012/08/16(木) 00:49:01.29 >>197 その差は、特に長手方向に粗いのが原因っぽい メッシュ細かくしてきゃ、どっちも値大きくなって、双方一緒っぽい値になる筈 テトラ一次は細かくして行っても、2次と全然違う値になったりするけど http://kizuna.5ch.net/test/read.cgi/sim/1284635327/200
201: 名無しさん＠５周年 [] 2012/08/16(木) 09:53:53.82 >>200 まぁそうなんだけどねｗ つまり、二次要素は、メッシュ切りがテキトーでもそれなりの精度が出る。 一次はメッシュの細かさが、モロ、精度に出るので、メッシュ密度とか形状に非常に神経を使わなきゃならんので 実務的には向いていないという考えです。 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/201
202: 名無しさん＠５周年 [] 2012/08/16(木) 12:39:52.17 >>201 テトラの2次と1次は 3割4割違ってくるので 比べるとマシだけど ちょっと差があるのは気になるね ウーム http://kizuna.5ch.net/test/read.cgi/sim/1284635327/202
203: 名無しさん＠５周年 [] 2012/08/16(木) 20:22:52.36 >>197 あ、書き忘れていた、これの理論値は -0.75mmです。 つまり 理論値 -0.75 1次６面体（8節点）要素 -0.72 ２次６面体（（20節点）要素 -0.63 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/203
204: 名無しさん＠５周年 [] 2012/08/16(木) 20:24:58.37 >>203 う、逆だった 理論値 -0.75 ２次６面体（（20節点）要素 -0.72 1次６面体（8節点）要素 -0.63 http://kizuna.5ch.net/test/read.cgi/sim/1284635327/204
205: 名無しさん＠５周年 [] 2012/08/18(土) 06:46:14.80 倍メッシュ切ればいいじゃない http://kizuna.5ch.net/test/read.cgi/sim/1284635327/205
206: 名無しさん＠５周年 [] 2012/08/18(土) 11:24:42.70 >>205 その「倍」っていうのはどこから来る数値なのか？ たとえばこういう単純ハリの場合、一次要素で解析するには、 縦横高さで何等分、アスペクト比をいくつにすればいい、 っていう目安ってありますか？ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/206
207: 名無しさん＠５周年 [sage] 2012/08/18(土) 11:32:30.25 「倍」というのが辺の長さを指しているのであれば、 3次元だから2の3乗で8倍のメッシュを必要とするという理解で良い？ メッシュ数が「倍」であれば、辺の長さを2の立方根で割れば良い？ http://kizuna.5ch.net/test/read.cgi/sim/1284635327/207
208: 名無しさん＠５周年 [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: 名無しさん＠５周年 [] 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
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 182 ﾚｽあります
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.862s*