☆有限要素法(Finite ..
[2ch|▼Menu]
2:名無しさん@5周年
10/09/19 11:46:57 .net
これ一から作ることになって死にそうなんだが。

3:名無しさん@5周年
10/09/19 15:46:56 .net
入出力部をいちから作るのはたいへんだが、計算部だけならそうでもないでしょ

4:名無しさん@5周年
10/09/19 20:42:28 .net
MATLABとかを使うと、行列計算はしやすいよ

5:名無しさん@5周年
10/09/20 00:25:11 .net
>>4
なんでそういう商用アプリを使うかな。
FreeMATなりOctaveを使えよ!

6:名無しさん@5周年
10/09/20 17:18:44 .net
オーソドックスにC言語でやるのは、かなりきついよね?

7:名無しさん@5周年
10/09/20 17:52:52 .net
メッシュ作成のサンプルどっかに転がってないかな

8:名無しさん@5周年
10/09/21 10:44:27 .net
オープンソースのADVENTUREとかは?


9:名無しさん@5周年
10/09/21 13:39:56 .net
不連続ガラーキン法って流行ってんの?

10:けいた
10/09/21 20:45:13 .net
 なんで、そういう発想になるのかな?
 物理が、分からないと、何やっても、ダメだって!!!
 差分法!!!!!


11:名無しさん@5周年
10/09/23 13:48:11 .net
>>5
それって、windows上で動きますか?

12:名無しさん@5周年
10/09/23 21:26:52 .net
Octaveは動く


13:名無しさん@5周年
10/09/24 21:07:24 .net
>>11
URLリンク(freemat.sourceforge.net)

Linux, MacOS Windows で動くらしいですね。

14:名無しさん@5周年
10/09/24 21:13:10 .net
普通 scicoslab

15:名無しさん@5周年
10/09/25 09:02:28 .net
abaqus student がインストールできないんですが
誰かわかりませんか?涙

16:DNS厨
10/09/29 17:59:13 .net
Spectral Element Method
とか
FEM
なんかGPLでソース配布してるプロジェクトが
少なからずある

17:名無しさん@5周年
10/10/03 01:52:06 .net
FEMの良書を天ぷらしてくれないかな

18:名無しさん@5周年
10/10/08 10:31:05 .net
>>15

普通にインストールできたけど…?
当方の環境はWindows XP + Abaqus Student Editon 6.9(書籍の付録)


19:名無しさん@5周年
10/10/08 14:23:36 .net
Abaqus Studentってのはフリーなの?有料ならサポがあるよな

20:こうちゃん
10/10/08 20:47:04 .net

 クソガキの、巣窟!!!!!


21:名無しさん@5周年
10/10/09 10:35:23 .net
>>19
フリーっていうか、1万円の本の付録扱い

22:名無しさん@5周年
10/10/10 18:25:23 .net
>>21
本に付録してるだけなのか?
広める気あるのかないのか・・・・w

23:名無しさん@5周年
10/10/10 22:07:31 .net
>>22
個人でAbacus使おうなんていうヒトはいないでしょwww

ちなみにAbacusっぽいやつなら CALCULIX

URLリンク(www.calculix.de)
URLリンク(www.bconverged.com) (FreeのWindows版あり)

へどうぞ。

24:名無しさん@5周年
10/10/19 23:50:06 .net
>>23
もっと詳しく!

25:名無しさん@5周年
10/10/24 01:03:08 .net
>>15
1万円でAbaqusが使えると聞いて、私も購入しました。
ですが私もインストールできません。

windows7だから?

26:名無しさん@5周年
10/11/13 19:41:54 .net
URLリンク(www.google.co.jp)

27:名無しさん@5周年
10/11/22 22:23:53 .net
よほど大規模解析でもしないのならFlexPDEなるソフトはどう?
評価版(要素数や支配方程式の数の制限あり)ならフリーでwindows上で動くよ。

28:名無しさん@5周年
10/11/25 21:14:20 .net
.

29:名無しさん@5周年
11/01/05 13:07:42 .net


30:名無しさん@5周年
11/01/12 23:01:12 .net


31:名無しさん@5周年
11/01/26 18:46:34 .net


32:名無しさん@5周年
11/02/05 19:30:42 .net
.

33:名無しさん@5周年
11/02/17 21:58:46 .net


34:名無しさん@5周年
11/02/21 15:36:32.25 .net
easy-σとかいうの大学で使った
原理はよくわからんけどとりあえず結果は出た

35:名無しさん@5周年
11/02/21 20:42:22.47 .net
使う側は原理なんて考えなくていいよ。
イレギュラーな結果を見定める時に多少知ってればいい程度だ。

36:名無しさん@5周年
11/03/06 17:48:37.29 .net
>>34
東大乙

37:名無しさん@5周年
11/04/15 01:37:38.95 .net
2次元構造解析で異方性材料使えるソフトってある?

38:名無しさん@5周年
11/04/15 18:14:07.95 .net
>>37
Ansysとか


39:名無しさん@5周年
11/04/16 01:39:09.66 .net
3次元図化したいのですが、フリーのソフトでオススメありますか?


40:名無しさん@5周年
11/04/16 21:10:57.23 .net
>>39
gnuplotとか


41:名無しさん@5周年
11/04/17 12:52:33.93 .net
gnuplot でfEMの結果を表示でキルンでしょうか?

単純なグラフならよさそうだけど。


42:名無しさん@5周年
11/04/18 20:05:30.33 .net
>>41
質問を小出しにするやつにFEMをやる資格はないと思え


43:名無しさん@5周年
11/04/18 20:11:44.15 .net
普通、VTKフォーマットにして paraview


44:名無しさん@5周年
11/04/18 23:31:19.08 .net
>>42
FEMのスレなんだから

45: 忍法帖【Lv=11,xxxPT】
11/04/19 04:43:15.50 .net
図のインパクトで考えるなら
やっぱParaViewだわな


46:名無しさん@5周年
11/04/21 15:54:22.71 .net
商用ソフト無駄機能多すぎて使いづらいよ…

47:名無しさん@5周年
11/04/22 09:50:55.59 .net
>>42
 品のないやつにエンジニアの資格はない。


48:名無しさん@5周年
11/05/09 19:44:50.11 .net
100万要素ぐらいでもサクサク動いて
くれるソフトありますか?

49:名無しさん@5周年
11/05/09 19:53:08.19 .net
そこまで行ったらソフトの良し悪しじゃないもんなぁ
粒子法でGPUバリバリって方があってるかもねw

50:名無しさん@5周年
11/05/10 19:58:40.80 .net
フリー、有料どちらでもいいんですが、粒子法のおすすめ
ソフトありますか?


51:名無しさん@5周年
11/06/06 22:06:54.20 .net
!ninja

52:名無しさん@5周年
11/07/31 19:50:46.55 .net
C言語とかで、地道にFEMを計算している人っています?

53:名無しさん@5周年
11/08/01 07:02:00.50 .net
います


54:スミレ
11/08/07 08:16:07.15 .net
層化の矢倉の婆が朝から異常行動繰り返して迷惑しているんですが
日に日にエスカレートしてまったくやめる気配がまったくありません



55:名無しさん@5周年
11/08/07 09:33:33.33 .net
FEMの場合は、どう考えてもFortranが多いでしょう


56:名無しさん@5周年
11/08/07 14:40:04.81 .net
てす

57:名無しさん@5周年
11/08/08 19:41:00.61 .net
境界要素法はどうでしょうか?


58:名無しさん@5周年
11/09/22 20:50:41.75 .net
ここを見ている人はいるのかな・・・

59:名無しさん@5周年
11/09/22 20:53:14.79 .net
なに

60:名無しさん@5周年
11/09/22 20:54:18.38 .net
>>57みたいに何を解くのか
なにを求めるのかも書かずに
アホのような質問されても
誰も答えないっしょ。

61:名無しさん@5周年
11/09/22 22:14:14.50 .net
今、有限要素法の勉強をちまちまやっているのだけど、難しくてさっぱり分からんw
Octaveでプログラムを書いているのだけど、境界条件を入れた全体剛性行列のdetがゼロになって解けん・・・
正六面体の一次要素一つの形状を作って片面を全固定(123自由度)、反対面に荷重をかけているのだけど
正六面体だと解けない理由ってのが有るのかな?まあプログラムの何処かが悪いのだろうけど
何から調べていけば良いのかアドバイスお願いします

62:名無しさん@5周年
11/09/22 22:21:28.48 .net
OpenFoamのコードを読め

63:名無しさん@5周年
11/09/22 22:25:09.16 .net
detがゼロって・・・・

64:名無しさん@5周年
11/09/22 22:34:11.21 .net
>>62
そんな大それたもの読めません(本気なのかー!)

65:名無しさん@5周年
11/09/22 22:39:49.80 .net
自分で組んでんなら
他人のコードはおもしろいよ。

66:名無しさん@5周年
11/09/22 23:38:05.51 .net
>>61
2次元問題だと解けるの?


67:名無しさん@5周年
11/09/23 00:10:15.05 .net
>>66
2次元の要素をまだ実装していないのです。申し訳有りません。
実装するにしてもまだ時間がかかりそうです
使える要素は四面体と六面体の一次要素です。
四面体も変なのですがまだ細かくつめていません
(時間はかかると思いますが、2次元問題の方も作ってみます)

68:名無しさん@5周年
11/09/24 00:28:00.68 .net
Microsoft Office Excelでゆうげんようそほうはかのうですか?

69:名無しさん@5周年
11/09/24 00:49:45.97 .net
>>68
かのうです

70:名無しさん@5周年
11/09/24 00:51:16.44 .net
ありがと
頑張ってみるわノシ

71:名無しさん@5周年
11/09/24 01:46:13.50 .net
手計算でも可能。

72:名無しさん@5周年
11/09/30 20:11:42.74 .net
CalculiXの遣い方をおしえてけろ

73:名無しさん@5周年
11/09/30 21:32:53.61 .net
URLリンク(www.str.ce.akita-u.ac.jp)
URLリンク(freecaetester.blog62.fc2.com)

74:名無しさん@5周年
11/10/01 09:36:05.83 .net
なんだGoogle検索ですぐ見つかるやつか(笑)

75:名無しさん@5周年
11/10/01 10:31:07.63 .net
それすらもできずに聞くバカがいるから仕方なくな

76:名無しさん@5周年
11/10/01 10:31:43.57 .net
それも理解できずにレスする更にバカのためにもな

77:名無しさん@5周年
11/10/01 11:01:47.00 .net
馬鹿ばっかりだな(^^;

78:名無しさん@5周年
11/10/01 12:27:41.04 .net
>>73
URLリンク(freecaetester.blog62.fc2.com)
ここから読んでいけば、ミトコンドリアでもわかりそうだ。
ブログ作者に感謝だな。


79:名無しさん@5周年
11/10/02 00:10:42.40 .net
コマンド一覧でもあればな

80:名無しさん@5周年
11/10/02 10:07:31.63 .net
>>79
コマンド一覧は
URLリンク(www.usennet.ne.jp)


81:名無しさん@5周年
11/10/02 10:14:16.53 .net
つうかその程度自分で見つけられない人間って・・・。

82:名無しさん@5周年
11/10/02 10:24:45.93 .net
まぁまぁ
Caliculixに関する情報は結構イッパイあるから、
全く初めての人にとっては、どこから読めばいいのか
さっぱりわからん場合もあるので。
 自分で見つけたれば、ググればいいし、
教えてあげたい人は教えればいいでしょ。


83:名無しさん@5周年
11/10/02 11:08:39.48 .net
dualcorepcでも動いた!感動した!

84:名無しさん@5周年
11/10/02 16:31:03.61 .net
誰も答えないっしょ。

85:名無しさん@5周年
11/10/05 01:36:50.84 .net
今日から訓練してる

86:名無しさん@5周年
11/10/06 22:47:29.92 .net
プログラムいじってるんだけど、
質問ある?

87:名無しさん@5周年
11/10/08 21:20:42.21 .net
並列FEMのフリープログラムって
どこかでdownloadできますか?


88:名無しさん@5周年
11/10/08 21:34:24.37 .net
はい

89:名無しさん@5周年
11/10/11 23:42:36.10 .net
>>87
さしあげるよ
PHPコードだけど

90:名無しさん@5周年
11/10/12 02:53:09.96 .net
どこでdownloadできますか?

91:名無しさん@5周年
11/10/12 17:09:53.34 .net
>>90
インターネッツ

92:名無しさん@5周年
11/10/16 20:48:19.73 .net
>>79
2ch病なのかなあ
コマンド一覧

コマンドー覧 と誤読したorz

93:名無しさん@5周年
11/10/22 20:41:14.40 .net
えばっばばっばばっばあっば。

BEMってどうんな

94:名無しさん@5周年
11/10/24 13:26:50.19 .net


95:名無しさん@5周年
11/10/25 19:06:55.40 .net
山梨ってスタッドレスいらなくね

96:名無しさん@5周年
11/11/03 03:14:36.04 .net
隠れるようにないけどもないでのねと

97:名無しさん@5周年
11/11/03 06:43:35.54 .net
URLリンク(www.caelinux.com)

98:名無しさん@5周年
11/11/19 23:27:03.63 .net
スタッドレスいらないの?ねぇ

99:名無しさん@5周年
11/11/20 02:02:55.10 .net
いらないんじゃねえの?

100:名無しさん@5周年
11/11/20 08:56:30.99 .net
Gmshから読み込んで3次元のtet4のコードを書いているんだが、paraview可視化の部分がうまくいかない。
手動で座標とconnectivityを書たテストは、うまく出ているので書き方は間違ってないと思うのだが?

でいま積んでる。知ってるヤツいる?


101:名無しさん@5周年
11/11/21 07:22:20.99 .net
電波テロ装置の戦争(始)
魂は幾何学、コピー出来る公安はサリンオウム信者の子供を40歳まで社会から隔離している
オウム信者が地方で現在も潜伏している
それは新興宗教を配下としている公安の仕事だ
発案で盗聴器を開発したら霊魂が寄って呼ぶ来た
<電波憑依>
スピリチャル全否定なら江原三輪氏、高橋佳子大川隆法氏は、幻聴で強制入院矛盾する日本宗教と精神科
<コードレス盗聴>
2004既に国民20%被害250〜700台数中国工作員3〜7000万円2005ソウルコピー2010ソウルイン医者アカギ絡む<盗聴証拠>
今年5月に日本の警視庁防課は被害者SDカード15分を保持した有る国民に出せ!!<創価幹部>
キタオカ1962年東北生は二十代で2人の女性をレイプ殺害して入信した創価本尊はこれだけで潰せる<<<韓国工作員鸛<<<創価公明党 <テロ装置>>東芝部品)>>ヤクザ<宗教<同和<<公安<<魂複<<官憲>日本終Googl検索

102:名無しさん@5周年
11/11/21 07:52:33.87 .net
>>100
使ったこともないわ〜。

103:名無しさん@5周年
11/11/23 23:22:47.95 .net

947 名前:彼氏いない歴774年[sage] 投稿日:2010/07/05(月) 15:45:01 ID:7xkp916Q

女子大生風二人。電車にて。
A「最近まとめサイト見てんだけどさー」
B「なんのまとめだろう」
A「コピペのまとめサイトだよ。わかる?」
B「あーー…わかる」
A「わかってなさそうだなー。ま、いいけど。モジョの人たちが聞いた会話のコピペが多くて」
B「モジョ…?」
A「超面白いよ。一回モジョ板も見てみ」
B「の…奇妙な冒険?」
A「は?まあとにかくさ、どっかでモジョの人が私たちの会話聞いたりしてんだよ」
B「ネアンデルタールの次あたり…?」
A「何言ってんの?つーかさ、私たちの会話がコピペのネタになってるとか超素敵じゃん!皆に教えたいじゃん!あ私ここで降りるね!また明日ー!」
B「(手を振りながら)モジョとは…なんじょ」
お嬢ちゃん、後ろ後ろ。

104:名無しさん@5周年
11/11/25 19:11:33.68 .net
>> 102

自己解決すた。vtkのconnectivityのノード番号って0から始めなきゃいけないのね

105:名無しさん@5周年
11/11/25 19:14:28.22 .net
日本語ユーザーマニュアルでも作ってみれば?

106:名無しさん@5周年
11/12/05 22:49:37.15 .net
strand7 ってソフト使わなきゃいけないんだけど、
使い方を簡単にまとめたサイトとか誰か知りませんか?


107:名無しさん@5周年
11/12/06 05:02:10.48 .net
delFEMってソフト使ってみようと思ったけど全然コンパイル通んないですね・・
これ推奨の環境ってMACなの?

108:名無しさん@5周年
11/12/06 06:39:40.45 .net
ページも読めない人には無理
商用使いな。

109:名無しさん@5周年
11/12/06 13:34:06.00 .net
んー一応VSとかminGWとか使ってみてるんだけど
エラー出まくる・・これはソース直して使うのがデフォ?

110:名無しさん@5周年
11/12/06 19:30:30.62 .net
エラーの内容が判らんようなやつには無理

111:名無しさん@5周年
11/12/06 22:42:20.41 .net
直せばコンパイル通すくらいはできると思うが、、
コンパイル通らない状態で公開されてるって認識で良いのだろうかってことが知りたい。

112:名無しさん@5周年
11/12/07 02:26:25.13 .net
メアド公開してんだから本人に聞けよ

113:名無しさん@5周年
11/12/07 19:28:32.98 .net

亀田も真っ青!前人未到の犯罪5階級制覇や!!

違法駐車王 塩ホ容疑者

器物損壊王 仕置容疑者

飲酒運転ほう助王 死沖容疑者

無法バーベキュー王 Jおき容疑者

児童虐待王 SALTWOOD容疑者



114:名無しさん@5周年
11/12/07 20:45:40.51 .net
>>112
ここできいたほうがはやいよ

115:名無しさん@5周年
11/12/08 06:21:58.20 .net
未だに回答ないじゃねえか。
アホか

116:名無しさん@5周年
11/12/10 04:46:29.74 .net
1.2.0でやったら動きました〜♪
お世話になりましたノシ

衝突解析でっきるっかな。

117:電脳プリオン 【大吉】 【392円】
12/01/01 21:56:45.57 ?2BP(1960).net
そういうのもあるのか

118:名無しさん@5周年
12/01/05 12:49:46.22 .net
URLリンク(peacelandia.blogspot.com)
2011年1月24日
田原総一郎の講演料
去年の話になりますが、10月19日に弘前市で行われた「原子燃料サ
イクル意見交換会」に関して、安藤晴美県議が情報開示請求し、その結
果、わかったことです。
企画は、資源エネルギー庁と青森県の共催で、(株)RABサービスに
委託されて行われており、委託料は129万6750円。その内訳は、
ジャーナリストの田原総一郎氏の講演料が110万円。交通費8万円。司会料3万円。管理手数料(食費含む)2万5千円。


119:名無しさん@5周年
12/01/06 11:27:33.20 .net
前に話題になってる、 Abaqus Student Editon 6.9(本のおまけのやつ)ってmacでも動きますか?

対応OSがわかれば教えて下さい。

120:名無しさん@5周年
12/01/06 19:48:58.12 .net
URLリンク(www.simulia.com)

121:名無しさん@5周年
12/01/07 01:37:12.58 .net
>>17
これとか URLリンク(www.amazon.com)

122:名無しさん@5周年
12/01/13 00:54:16.16 .net
Abaqusの固有振動数抽出がまったくうまくいかん
単位もN mm t 系でやってんのに
メッシュをめちゃくちゃ細かくしたら少しだけ理論値に近づくが・・・

123:名無しさん@5周年
12/01/13 01:02:09.70 .net
Abaqusだろうがどれでも
基本的に数値は大差ない。

変ならモデルに妙な接合部があるとか
材料データ間違ってる。

124:名無しさん@5周年
12/01/13 01:16:41.11 .net
単純な片持ち梁なんだけど、どこかおかしいのかな。また明日みてみるか
スケッチで50*50*1000の梁を作って、密度7.8e-9、剛性206e3、ポアソン比0.3
を入れてるだけなんだけど・・・

125:名無しさん@5周年
12/01/13 06:49:13.48 .net
その数値にそれぞれ単位を入れてみて

126:名無しさん@5周年
12/01/13 06:53:04.34 .net
あと断面形状がなにかと
現状でてる結果の値も

127:名無しさん@5周年
12/01/13 23:27:09.84 .net
>>124
URLリンク(www.42ch.net)

poor man's abaqus であるところの calculixでやってみたけど、
41.57Hzじゃない?


128:名無しさん@5周年
12/01/13 23:47:58.62 .net
メッシュ図をつけたもの
URLリンク(www.42ch.net)

ちなみに理論値は

(1.875/1000)^2*sqrt((206000*50*50^3/12)/(7.8e-9*50*50))/(2*pi)
=41.504 Hz

だから十分あっていると思う。

あ、念のためにだけど、要素は二次要素(*ELEMENT, TYPE=C3D20, ELSET=Eall )
にしたほうがいいよん






129:名無しさん@5周年
12/01/14 00:19:18.17 .net
固有値解析のポスト処理について質問です。
諸々の事情によりParaViewを使いたいと考えています。
変形表示はWarpByVectorでできそうなのですが、
これをさらにクネクネアニメーションさせたいです。
どんなやり方がありますか?

130:名無しさん@5周年
12/01/14 00:49:36.40 .net
>>125-128
レスありがとうございます!
今日(きのう)もう一度、最初からモデルを作って解析をすると、

Mode1 41.439
Mode2 41.448
Mode3 256.72
Mode4 256.78
Mode5 706.19
Mode6 706.34

断面:50*50[mm]正方形
長さ:1000[mm]
密度:7.8e-9 [t/mm3]
剛性:206e3 [N/mm]
Poisson:0.3

となって、理論値に非常に近くなりました。どうやらメッシュの切り方等に問題があったようです・・・
単位系ばかりに気を取られていて

しかしAbaqus、メッシュを少し荒目にするとMode1が30Hz位になって誤差が大きくなりますね。
静的な応力解析の時は、メッシュの大きさが大きくても、そこそこの精度は出るのに。
詳しく説明をしていただき、本当にありがとうございました。
ソフトの使い方だけでなく有限要素法についても、さらに勉強しようと思います。

131:名無しさん@5周年
12/01/14 03:13:22.84 .net
おつかれ!

132:名無しさん@5周年
12/01/14 03:13:57.65 .net
メッシュが一次要素になってない?
二次じゃないとおかしな値になるよ。

133:名無しさん@5周年
12/01/14 03:17:41.18 .net
>>129
URLリンク(d.hatena.ne.jp)

134:名無しさん@5周年
12/01/14 09:07:16.31 .net
>>130
うーん。全体を4要素という極めて少ない要素数にしても
URLリンク(www.42ch.net)

43Hzくらいで、まぁまぁの値になるんですが。
URLリンク(www.42ch.net)


135:129
12/01/14 09:43:05.09 .net
>>133
単なる非定常計算結果のアニメーションは
できるんですが、固有値解析のように
モードシェイプに対応する変位データが
ひとつしかない時に、変形前と変形後の
ノード座標値を(マイナス変位側まで)
ParaView上で補完して、物体が振動して
いるようなアニメーションがしたいのです。
FEMAPなどで固有値解析結果を表示するのと
同じことを、ParaViewでできないかなぁと。

変位データは自分で作成しているので、
連番ファイルで補完データを書き出すことも
不可能ではないのですが、数十個ある
モードシェイプ情報に対して全て補完データを
出力すると連番ファイルが1000個近くなるので、
何かうまいやり方があれば助かります。

136:名無しさん@5周年
12/01/14 10:07:08.76 .net
>>132
1次要素って、2次要素より「曲がらない」=「硬い」ので、周波数としては上がる傾向に あるみたいです。

4要素だと、2次要素なら43Hzですが、1次要素(C3D8)にすると77Hzになっちゃいます。


137:名無しさん@5周年
12/01/14 11:31:58.53 .net
そんなに高くなる?

138:名無しさん@5周年
12/01/14 11:33:42.41 .net
あ、4要素の場合か。

139:名無しさん@5周年
12/01/14 13:59:38.61 .net
神道教育が起源の、創価学会から、はじまった公明党には、票を入れられません。なぜなら、日本は、国教国家ではなく、世俗国家
なのです。これは、ほぼ違憲であり、戦後の、反省の立場から、公明党が政権をとることを想定すると、神道国家に回帰するため、
アジアや旧連合諸国にたいする国際的なプレゼンスや国際関係を著しく悪化するものと思われます。また、別の観点から、創価大学を持っているといえども、「学会」と言う名称も、国際的には、本来の学会の意味でないため、非常に不適切な団体です。
創価学会との無関係を訴える公明党ですが、フジテレビの選挙報道時は、「創価学会のみなさん」と不自然なくらい必ず、放送します。
ブッシュドクトリンに協調した小泉自民・公明は、「テロとの戦い」を表明するも「オウム真理教」に破防法を適用できない失敗をし、野党に追いやられました。
つまり、宗教団体は、特に、世俗国家においては、政党を持つべきではありません。
小泉元総理は、公明党と公約を守って、「自民党をぶっこわす」だけでなく、原発も
ぶっこわしたのです。これを小泉の奇跡と言いましょう。世襲制にも反対しましょう。

140:名無しさん@5周年
12/01/19 03:22:23.57 .net
質問させていただきます。
Abaqusで、粘着挙動を選択しますとデフォルトでペナルティー法による粘着となっています。
このペナルティー法というのは、どういったものなのでしょうか。
調べてみますと、FEM解析における接触の基礎としてのペナルティー法や、計算が速い代わりに食い込みを許すということ等、おおまかなことは分かったのですが、
このペナルティー法でA部材とB部材を粘着させると、完全に一体として接着されたものとみるのでしょうか?

具体的には、A部材をStep1で伸ばしておき、そこにStep2でB部材を張り付けるということをしたいのですが、張り付けに、ペナルティー法の粘着を使用するとどうなるのかがわかりません。

141:名無しさん@5周年
12/01/19 06:22:30.22 .net
ヘルプぐらい嫁。
購入先のサポートに問い合わせろ。


142:名無しさん@5周年
12/01/26 18:06:20.34 .net
C言語で自分で作ろうよ

143:名無しさん@5周年
12/02/06 19:24:54.48 .net
すいません、この中でFemap with NXNastranで接触解析をしたことがある人はいますか!?

SOL601(NX Adv Nonlin)で二次元の弾性体半球同士の接触解析を強制変位を使ってしようとしているのですが、Connectionsが上手く働いていないのさ、両物体が思い切り貫入するような結果になってしまいます。収束も最後までしません。

Slide Lineを使った二次元解析ではなく、ConnectとSOL601を使って接触解析をするなら、プロパティやリージョンをどのように設定すれば良いのでしょうか!?アドバイス下さい!

144:名無しさん@5周年
12/02/09 01:19:33.60 .net
板要素のDKMQ要素の理論について知っている人いますか?

離散kirchhoff板理論(DKT,DKQ)とひずみ仮定(MITC8)を
拡張された理論であるのは分かったのですが、
ここから先はさっぱり分かりません。

なにか参考になるソースはありますか?


今見ている文献は
A new discrete Kirchhoff-Mindlin element based on Mindlin-Reissner
plate theory and assumed shear strain fields. PartU

145:名無しさん@5周年
12/03/20 09:45:59.80 .net
有限要素法の理論的な部分に詳しい著書を教えてください

146:名無しさん@5周年
12/03/20 11:00:38.11 .net
読み込むならこのへんかな
URLリンク(books.google.co.jp)





147:名無しさん@5周年
12/03/20 16:42:07.01 .net
URLリンク(www.brandany.com)<)
URLリンク(www.lv-market.net)<)

148:名無しさん@5周年
12/03/20 20:15:37.62 .net
>>146
ありがとうございます
その本が図書館にあるか調べてみます
もし、アマゾンなどで簡単に購入できる本があれば教えてください

149:名無しさん@5周年
12/04/11 18:02:08.05 .net
C言語とかで組み立てている人います?

150:名無しさん@5周年
12/04/12 00:00:49.11 .net
FEMはFORTRAN。 例外は認めない!
なんつったってAnsysだってFORTRANだ。


151:名無しさん@5周年
12/04/15 03:58:31.06 .net
この外人、日本語うますぎワロタwww
URLリンク(www.youtube.com)

152:名無しさん@5周年
12/04/21 22:04:52.88 .net
有限要素法ってナニに使うの?

153:名無しさん@5周年
12/04/23 21:08:01.92 .net
有限要素法が必要な計算に使うよ

154:名無しさん@5周年
12/05/03 11:03:26.71 .net
有限要素法はね、
個体の荷重ー変位曲線を割り出すのに使います

155:名無しさん@5周年
12/05/04 19:32:44.52 .net
<Sale>中国語版の書籍<IT|3D|NN|理系|上古音などなど>(格安) 280円より
http◎//lang-8○com/194279/journals/1457398/

156: 【38.1m】 電脳プリオン
12/05/08 22:54:02.60 ?PLT(12079).net
無限要素法ってある?

157:名無しさん@5周年
12/05/09 13:42:39.90 .net
#EANF#

158:名無しさん@5周年
12/05/13 00:59:26.81 .net
致違いで申し訳ございません。
どこにもこれから質問するものに関連するものがなかったのでここで質問させていただきます。

今現在研究で(大学生4年です)解析をしています。
解析手法としては線形加速度法という陰解法を使用していますが、
この手法の欠点としては、ステップ間隔を短くしなければなりません。

実際研究で使う地震波は、100秒もあり、扱う対象のものは
固有円振動数だけで330個、その中の最小のものより1/60程度にしなければ
なりません。そうしなければ、解は拡散してしまうからです。

そこで、0.0005秒間隔で100秒までとなると20万ステップしなければなりません。
構造解析のソフトウェアではそれにどれだけいったいかかるのか
図ったことはありませんが、最適化ですのでこれを少なくとも100回は回していかなければ
なりません。線形加速度法で扱う行列は、330行列。。。
この前0.0005秒間隔で20000ステップやるのですら、
1時間かかていました。。 20万ステップですから単純計算で
一回終えるのに10時間。これが最低でも100回回さなければならないため、
1000時間かかってしまいます。

そこで、解析手法を変えようと思っているのですが、
どの方法が適切でしょうか?何卒ご協力お願い申し上げます

159:名無しさん@5周年
12/05/13 15:23:32.07 .net
>>158
>最適化ですので
何をどう最適化したいのか、エスパーしようにも
どうしたものやら。
地盤?それとも構造物?

真面目に変位量の時間履歴を取りたいのなら
愚直に解くしかないんじゃないの?
それをどうやったら他の手法で置き換えられるのか、
それはむしろ君自身が研究テーマとして取り組むべき
問題なんじゃないかと思うのだが。

いま京の研究公募してるから、それに応募してみたら?
今の手法でも結構頑張れちゃうかもよ?

160:名無しさん@5周年
12/05/13 18:26:15.98 .net
>>158
時間が限られてるなら最適化の項目を減らしなさい。

161:名無しさん@5周年
12/05/13 18:26:38.22 .net
どうしてもやる必要があるなら数年掛かろうがやりなさい。

162:名無しさん@5周年
12/05/19 20:06:02.61 .net
>>158
そのまえに0.005秒のデータってあんの?
俺は15万要素くらいで60秒解析してたけど、実際0.05秒とかでも、大して結果かわんなかったぞ

163:名無しさん@5周年
12/05/19 21:29:22.63 .net
そもそも地震関係やっていて
そんな周波数を解析するって時点で無能君。

164:名無しさん@5周年
12/05/19 21:44:28.77 .net
いや、こいつのテーマは
様々な地震動に対して最適な構造を見つけるってことじゃねえのか
最適化ってそういうことだと思ったんだが
つーか書いてて思ったが1000時間って大したことないよね
パソコンがやるんだし。連続して解析できるんだろ?
2台使えば500時間だ。一ヶ月もかからねえぞ

165:名無しさん@5周年
12/05/19 21:46:50.18 .net
そういえば、数年前にノートPCで解析してるペーパーがあった。
防衛大だったと思う。飛行機とか衝突してたから
ただ、終わりのほうに、ノートPCでも3ヶ月かければ解析は可能だって書いてて、
こいつは何を伝える気なんだと思ってワロタ


166:名無しさん@5周年
12/05/24 23:46:58.24 .net
>>95
要らないよ。
チェーンで十分だ。

167:名無しさん@5周年
12/05/25 05:26:40.14 .net
>>166
良かった!
心配でずっと眠れなかったよ。

168:名無しさん@5周年
12/05/27 09:07:06.21 .net
シーズンすぎてるだろ…

169:名無しさん@5周年
12/07/04 00:11:21.07 .net
OpenFOAMのStressAnalysisに詳しい方は居ませんか?
集中荷重の指定方法とか何故繰り返し計算を強いられるのかとか教えてください
分かっていない事が多すぎるのでソースを読めはご勘弁を

170:名無しさん@5周年
12/07/05 19:35:22.36 .net
>>169
OpenFOAMでStressAnalysisをする時点で、アレ。
calculix とか、そーいうの使うほうがいいかと。


171:169
12/07/06 20:03:37.71 .net
>>170
StressAnalysisはオマケみたいなものなのですかね・・・
Calculixも平行して為してみているのですがプリ部分でまだ試行錯誤していますw
FreeCADでSTLファイル作成→netgenでメッシュ作成→cgxで荷重、境界条件作成→ccxで計算→cgxで結果を確認
と言う感じにやってますが、単一部材のケースはわりとすんなり行くのですが複数部材のケースがどうしたら良いのか
チラ裏でし

172:名無しさん@5周年
12/07/07 02:56:03.46 .net
>>171
アセンブリは素直にcgxでモデリングする方がてっとりばやいかと。

pythonが使えれば
URLリンク(dip28p.web.fc2.com)
っていうのが使えるらしい。


173: ◆uJTwg2na76
12/07/07 08:40:55.47 .net
どなたか塩谷豊という人物をご存じありませんか?
どんな情報でも結構です
どうかよろしくお願いいたします

174:169
12/07/07 18:11:54.59 .net
>>172
cgxか〜複雑なモデリングには向いてなさそうだけど使ってみるよ

netgenまでは複数部材は考慮されているみたいで、サンプルで提供されている
複数部材を定義されたgeoファイルを読み込み、abaqus形式で吐き出すと
要素の集合名が分けられて出てきます

CAD系で複雑なモデリングをしていくのが理想に思えるけどここらが自分の限界なのかな
STEP形式ファイルは複数の部材を考慮した出力って出来るのでしょうか?
あ、FreeCADで複数部材のモデリング方法が分からないからダメか・・・チラ裏でした

175:名無しさん@5周年
12/07/08 01:30:51.56 .net
>>174
Ass’yの難しさとして、

□■
A B

部材Aと部材Bが一つの面でくっついているとすると、
CGXみたくモデラーで作るときは共通平面ということで作っていけばいいけど
CADからもっていくとAとBの共通平面は別々の平面になっているはず。
そこでAとBをくっつけなきゃならないんだが、たいてい、その別々の平面を
接触(固着)させて解析させることになる。
 ここがちょっと難しいところですね。


176:名無しさん@5周年
12/07/08 01:39:06.96 .net
モデリングでは境界条件が難しいよね。
平面と平面が面で完全に接合している物体なんて現実的ではないし。
割り切るしかないんだろうね…。

177:169
12/07/11 23:42:31.43 .net
複数部材でできたモデル(片持ち梁)をうまく作ることが出来たからメモ代りにカキコ
1) FreeCADで梁を2つ(接触した形)作成し、2つ選択した状態でSTEP214形式でexport
2) gmshでSTEPファイルをオープンする。一旦geoファイルに保存しそれで開き直す。
3) Geometry:Elementary entities→Coherenceで隣接したジオメトリの節点を結合する。
4) Geometry:Physical groups→Add→Volumeで2つPhysical groupを作成する。
5) Mesh:3Dでメッシュ作成、メニューのFile→Save asでAbaqus INPを選択して適当な名前でメッシュファイルを作る
6) cgx -c メッシュファイル名 で開き境界条件、荷重を作る。適当にccx用のインプットファイルを作成して解析、結果確認

分かっている人には当たり前の事かもしれなかったけど、これを分かるまで大分時間がかかったよ
複雑な形状だと節点が重ならなかったり上手く行かないんだろうな・・・

178:名無しさん@5周年
12/07/12 22:48:50.21 .net
>>177
>3) Geometry:Elementary entities→Coherenceで隣接したジオメトリの節点を結合する。

お、これは知らなかった。一回やってみよ。さんくすです。


179:名無しさん@5周年
12/07/15 23:51:43.05 .net
四面体一次要素の形状関数を教えてください。とても複雑になる様な気がしてます・・・
要素の体積を A として,要素内の点が与えられたとき,対応する節点の対面と
その点が作る四面体の体積をA1,A2,A3,A4 として
N1=A1/A, N2=A2/A, N3=A3/A, N4=A4/A
となるらしいのです
節点をP1,P2,P3,P4としてmaxia風に書くと
A: 1/6*determinant( transpose(matrix((P2-P1), (P3-P1), (P4-P1))) );
要素内に与えられる点をP0とすると
A1: 1/6*determinant( transpose(matrix((P2-P0), (P3-P0), (P4-P0))) );
A2: 1/6*determinant( transpose(matrix((P1-P0), (P3-P0), (P4-P0))) );
A3: 1/6*determinant( transpose(matrix((P1-P0), (P2-P0), (P4-P0))) );
A4: 1/6*determinant( transpose(matrix((P1-P0), (P2-P0), (P3-P0))) );
となって複雑に・・・何かが間違えている

180:名無しさん@5周年
12/07/17 21:26:05.61 .net
 > 169
OpenFOAMってVOF法みたいので構造も有限体積で計算していたはず。
弱連成もできると聞いたことあるので、有限要素のソルバーもあるかも知れんが。

詳しくは開発元が書いた論文よめ
URLリンク(onlinelibrary.wiley.com)


181:名無しさん@5周年
12/07/22 13:21:24.64 .net
URLリンク(jp.alpha-lab.info)
ニューラルネット超入門

182:名無しさん@5周年
12/08/04 21:34:49.30 .net
機械部品の強度解析するのに、FEMのソフトを導入しようと思うのですが、
ANSYS、NASTRANなどを検討していますが、何かお勧めはありますか?

メインは線形で、非線形も少しやろうかと思っています。

183:名無しさん@5周年
12/08/04 22:07:19.48 .net
運用実績が無いなら、まわりで使っているソルバー又は聞ける人がいるソルバー
なんちって

固有値計算プログラムが上手く動かん、固有ベクトルがぐちゃぐちゃだ
固有値はハウスホルダー+QR法で何とかoctaveのeig関数の結果と一部一致したけど
(ゼロに持っていく値が振動して収束しない時があるんだよね…)
三重対角行列の固有ベクトルを求めて(怪しい実装)それをハウスホルダーの副産物を使って逆変換(とても怪しい実装)
見事に撃沈、誰かエスパー求む

184:名無しさん@5周年
12/08/07 16:21:53.17 .net
>固有ベクトルがぐちゃぐちゃだ

メッシュがグチャグチャなのが原因と予想だな

185:名無しさん@5周年
12/08/07 20:44:37.89 .net
>>184
片持ち梁、片方を3成分固定、片方を集中荷重をあたえた時にそれらしい結果が出たデータを使ってます
要素は8節点6面体の一次要素です。メッシュがぐちゃぐちゃとは要素座標とかそこらの考え方がぐちゃぐちゃしているのか・・・
剛性マトリクスを作る時についでで質量マトリクスを作っているからそこらが怪しいと思っているんだけど
積分点の位置で形状関数を使って質量を計算してはいけないのかな?
仕事とか学生とかではないから焦ってはいないんだけど、資料が無くて分からない事が多いんですよね
暇な時に少し知恵を貸してくれると幸いです

186:名無しさん@5周年
12/08/07 23:04:53.86 .net
>>182
まずCalculixみたいなフリーのやつ入れて、勉強したほうがいいと思う。
Salome-Mecaが英語だったら、こっちがオススメなんだけど、フランス語なんで日本人には辛い。
3DCADは何かつかってますか?
SolidWorksだったら、SolidWorksアドオンのCAE(COSMOS?)を載せたほうがいいと思うんだけどさ。


187:名無しさん@5周年
12/08/08 01:44:58.48 .net
>>186
レス有難う
calculixは一応入れてソースを見てるんだけど・・・コンパイルが悪いのか固有振動数計算が解けないんですよね
(*ERROR in arpack: info=-3 # of converged eigenvalues=70 とか出ます)
3DCADはFreeCADを使ってますが、形状がそんな複雑でないのとメッシャーが四面体要素中心なのでテストデータは手で作ってます
使える要素で静解析の結果が一番良かったのが六面体一次要素なのでそれを中心に進めています
使っているプラットフォームがdebian(squeeze 32bit)なのと貧乏なので商用の物は使えないのですが
SolidWorksですかちょっと調べてみます

188:名無しさん@5周年
12/08/08 18:35:14.23 .net
>>182
済まない、レスを混同してしまいました>>187は無かったことに

189:名無しさん@5周年
12/08/08 23:06:50.32 .net
>>187
構造解析なら、メッシュ数おとしてでも、2次要素を使ったほうが絶対いいです。
絶対っていうのは言いすぎかw
だけど、自分の経験からいうと、メッシュ数がテキトーでも二次要素だとかなり精度よく解けます。


190:名無しさん@5周年
12/08/09 01:21:56.74 .net
>>187
二次要素とか気にしなくて良いANSYS DesignSpaceにしとけば?
非線形は必要になったらライセンスアップグレードすれば良い。

191:名無しさん@5周年
12/08/09 01:23:40.38 .net
固有振動数レベルだとANSYS、NASTRAN、ABAQUSどれも結果変わらないし。

192:名無しさん@5周年
12/08/09 01:41:20.54 .net
>>189
六面体の2次要素はまだ実装していないのです

>>190,191
すみません、レスを混同して変なレスをしています

193:名無しさん@5周年
12/08/10 12:36:54.72 .net
テトラは2次でないと全然駄目だけど 6面体は一次で十分です

固有値グチャグチャは気になるけど 固定が駄目なのかなー

194:名無しさん@5周年
12/08/11 07:48:08.08 .net
今時二次がなんていつの時代のソフトだよw

195:183
12/08/11 13:13:16.14 .net
>>193
固定か〜固定している所の剛体行列と質量行列をどのように扱ったら良いか分からないから固定条件無しでやっています
まず剛体モードが綺麗に直行で出てくれないし何か変だ〜
今は固有ベクトルの計算は逆反復法らしきものでやっています
(上手く収束しないんだよねと言うかほとんど収束しない、振動してたり反復の上限一万回に引っかかってたりするよ)

>>194
手作りだからね・・・まずは一次要素からかなと思って作ってます
(ちょこちょこ弄っていて大分汚くなってきたよw)

196:名無しさん@5周年
12/08/11 13:18:46.78 .net
仕事のためじゃなくて勉強か

197:名無しさん@5周年
12/08/15 10:58:16.91 .net
>>193
いやいやいやいや、
6面体だろうが、なんだろうが、2次要素と1次要素じゃ精度が全然違う。
普通の構造解析で積極的に1次要素を使う理由は、全く無い。

URLリンク(www.42ch.net)


198:名無しさん@5周年
12/08/15 15:57:16.77 .net
2次なんて言ってる時点で終わり。分からないの?

199:名無しさん@5周年
12/08/15 18:18:06.88 .net
>>198
わからない。
学問板なので、理論的に説明してくれ。


200:名無しさん@5周年
12/08/16 00:49:01.29 .net
>>197
その差は、特に長手方向に粗いのが原因っぽい
メッシュ細かくしてきゃ、どっちも値大きくなって、双方一緒っぽい値になる筈  

テトラ一次は細かくして行っても、2次と全然違う値になったりするけど

201:名無しさん@5周年
12/08/16 09:53:53.82 .net
>>200
まぁそうなんだけどねw
つまり、二次要素は、メッシュ切りがテキトーでもそれなりの精度が出る。
一次はメッシュの細かさが、モロ、精度に出るので、メッシュ密度とか形状に非常に神経を使わなきゃならんので
実務的には向いていないという考えです。


202:名無しさん@5周年
12/08/16 12:39:52.17 .net
>>201
テトラの2次と1次は 3割4割違ってくるので 比べるとマシだけど
ちょっと差があるのは気になるね ウーム

203:名無しさん@5周年
12/08/16 20:22:52.36 .net
>>197
あ、書き忘れていた、これの理論値は -0.75mmです。
つまり

理論値               -0.75
1次6面体(8節点)要素      -0.72
2次6面体((20節点)要素    -0.63


204:名無しさん@5周年
12/08/16 20:24:58.37 .net
>>203
う、逆だった


理論値               -0.75
2次6面体((20節点)要素    -0.72
1次6面体(8節点)要素      -0.63


205:名無しさん@5周年
12/08/18 06:46:14.80 .net
倍メッシュ切ればいいじゃない

206:名無しさん@5周年
12/08/18 11:24:42.70 .net
>>205
その「倍」っていうのはどこから来る数値なのか?
たとえばこういう単純ハリの場合、一次要素で解析するには、
縦横高さで何等分、アスペクト比をいくつにすればいい、
っていう目安ってありますか?


207:名無しさん@5周年
12/08/18 11:32:30.25 .net
「倍」というのが辺の長さを指しているのであれば、
3次元だから2の3乗で8倍のメッシュを必要とするという理解で良い?
メッシュ数が「倍」であれば、辺の長さを2の立方根で割れば良い?

208:名無しさん@5周年
12/08/18 13:16:14.14 .net
正六面体に近い形でメッシュを切れば良いと思うの

209:183
12/08/18 17:37:04.75 .net
皆遊んでいる様なので質問を・・・
固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか?
ググっているのですがそれらしい文献を見つけられずに困っています。

詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください
固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを
逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで
プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに
したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード
を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると
サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると
ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー+QR法+逆反復法で
出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します
ここらもエスパーレスしてくれると有難いです

テストデータは次の感じです
・方持ち梁
・六面体一次要素(要素数1 と5これが主 と 625これは時間が足りなくて計算できない)
・拘束条件無し(なんかこれがいけないっぽい)

210:名無しさん@5周年
12/08/18 18:22:03.30 .net
>>209
>・拘束条件無し(なんかこれがいけないっぽい)

その通り。

固定条件の入れ方って、単に、固定する節点の変位を0にするだけじゃないの?

URLリンク(www.jsce.or.jp)

211:名無しさん@5周年
12/08/18 18:25:46.19 .net
>>208
その通りだと思います。
ただ、ちょっと複雑な形状をきると、アスペクト比1を目指してメッシュを切ると
要素数が非常に多くなることがあるので、1次要素を使ったほうが計算時間が
かかってしまう事はあるようです。


212:名無しさん@5周年
12/08/18 18:48:00.38 .net
>>210
スライドできる状態=摩擦ゼロ状態いうのは、厳密に起こり得ないので、
条件に対し解を求めるいう点では合ってるけど

スライドできるいう条件設定が、果たして現実的か?  そこはかなり難しい

例えば土木や建築の分野だと 地盤をモデル化したり、剛性体を挟んだり
色々しないと、単品+拘束やすべり条件設定では、合わない

213:183
12/08/18 19:01:45.09 .net
>>210
変形の計算(静計算って言うのかな?)はその通りだと思うのだけど、固有値、固有モード計算と同じなんですか?
変形計算の時って [K]{U}={F}で固有値計算の時は[K]{U}=λ[M]{U}ですよね。これをどう節点変位ゼロと
絡めたら良いのか分からないです。頭が悪くて申し訳ない

214:名無しさん@5周年
12/08/18 19:49:21.46 .net
>>213
固有振動の場合は固定点がなくても解が出るはず。
収束計算させるんじゃなくて、ガウス吐き出し法とかで解けば?


215:183
12/08/18 21:13:36.78 .net
>>214
境界条件がある固有値、固有ベクトル計算も知りたいので宜しくお願いします。

216:名無しさん@5周年
12/08/19 21:51:17.23 .net
>215
固有値=フリーでの計算ばかり

固定など条件設定を伴ったものは見た事がない気もします  できないのかも知れません



217:名無しさん@5周年
12/08/19 22:09:42.27 .net
>>216
出来ないってこたーないだろ。
AnsysだろうがNastranだろうが、出来ないCAEを見たことねぇぞ



次ページ
最新レス表示
スレッドの検索
類似スレ一覧
話題のニュース
おまかせリスト
▼オプションを表示
暇つぶし2ch

386日前に更新/101 KB
担当:undef