有限要素法のフリーのソフトウエア を試す

CAEで使われるフリーの有限要素法(FEM)ソフトウエアであるCalculixやOpenfoamを使ってみようとチャレンジしてたその足跡を残す。。。ついでに他のフリーソフトや商用ソフトの無償版にも手を出してみる。

calculix熱応力解析機能の確認-1要素編(その2)

昨日のデータではx=0の面を拘束しましたが、さらにx=1の面も拘束します。

これにより、x方向に伸びることができず、熱応力が発生するはずです。

*NODE, NSET=Nall
1,1.000000000000e+000,1.000000000000e+000,0.000000000000e+000
2,1.000000000000e+000,0.000000000000e+000,0.000000000000e+000
3,1.000000000000e+000,0.000000000000e+000,1.000000000000e+000
4,1.000000000000e+000,1.000000000000e+000,1.000000000000e+000
5,0.000000000000e+000,1.000000000000e+000,0.000000000000e+000
6,0.000000000000e+000,0.000000000000e+000,0.000000000000e+000
7,0.000000000000e+000,0.000000000000e+000,1.000000000000e+000
8,0.000000000000e+000,1.000000000000e+000,1.000000000000e+000
9,1.000000000000e+000,5.000000000000e-001,0.000000000000e+000
10,1.000000000000e+000,0.000000000000e+000,5.000000000000e-001
11,1.000000000000e+000,5.000000000000e-001,1.000000000000e+000
12,1.000000000000e+000,1.000000000000e+000,5.000000000000e-001
13,5.000000000000e-001,1.000000000000e+000,0.000000000000e+000
14,5.000000000000e-001,0.000000000000e+000,0.000000000000e+000
15,5.000000000000e-001,0.000000000000e+000,1.000000000000e+000
16,5.000000000000e-001,1.000000000000e+000,1.000000000000e+000
17,0.000000000000e+000,5.000000000000e-001,0.000000000000e+000
18,0.000000000000e+000,0.000000000000e+000,5.000000000000e-001
19,0.000000000000e+000,5.000000000000e-001,1.000000000000e+000
20,0.000000000000e+000,1.000000000000e+000,5.000000000000e-001
*ELEMENT, TYPE=C3D20, ELSET=Eall
1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
11, 12, 17, 18, 19, 20, 13, 14, 15, 16
*MATERIAL, Name=mat1
*ELASTIC
1000, 0.0001
*EXPANSION,ZERO=273.
1.E-6
*SOLID SECTION, Elset=EALL, Material=mat1
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,273.
*STEP
*STATIC
*TEMPERATURE
NALL,373
*BOUNDARY
1,1,3,0
2,1,3,0
3,1,3,0
4,1,3,0
5,1,3,0
6,1,3,0
7,1,3,0
8,1,3,0
9,1,3,0
10,1,3,0
11,1,3,0
12,1,3,0
17,1,3,0
18,1,3,0
19,1,3,0
20,1,3,0
*NODE FILE
U
*EL FILE
S, E
*END STEP

結果はこうなりました。

1要素熱応力
クリックすると大きくなります。

綺麗なコンター模様が出ていますが、レンジを見るとすべて-0.1になっています。

手計算してみると、

拘束なしの状態でののびは昨日の計算から、1.0e-4

1辺長さ1なので、ひずみも1.0e-4

これが拘束された時に起きる応力は、ヤング率を1000、圧縮方向の応力になることを考慮して

-1.0e-4x1000=-0.1

となり、シミュレーションと一致します。



***これはフリー有限要素法ソフトウエアCalculixについての説明です。ABAQUSコマンドの説明ではありません***

にほんブログ村 科学ブログ 技術・工学へ

技術・工学ブログランキング





スポンサーサイト

テーマ:宇宙・科学・技術 - ジャンル:学問・文化・芸術

  1. 2008/08/27(水) 21:01:32|
  2. Calculix熱応力解析、伝熱解析
  3. | トラックバック:0
  4. | コメント:0
<<熱応力解析の入力説明 | ホーム | calculix熱応力解析機能の確認-1要素編>>

コメント

コメントの投稿


管理者にだけ表示を許可する

トラックバック

トラックバックURLはこちら
http://freecaetester.blog62.fc2.com/tb.php/68-ef06a10c
この記事にトラックバックする(FC2ブログユーザー)

プロフィール

tsunodako

Author:tsunodako
tsunodakoのプロフィール
無料でCAE環境を構築に挑戦中。
が、現在少し休止中。そのうち復活します。
ご意見ご感想をブログコメント大歓迎です!!!
メールの場合は下のほうにあるフォームに入力していただくか

tsunodako[at]yahoo.co.jpまで([at]=@に変換してください)。
このブログにあるソフトやデータを使用して生じた結果は責任負いませんのであしからず。
静解析CAE技術者のための動解析講座は http://vibrationcae.blog.fc2.com/ に移動いたしましたのでよろしく。

ブログ内検索

カレンダー(月別)

04 ≪│2017/05│≫ 06
- 1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30 31 - - -

最近の記事

カテゴリー

tsunodakoへのメールはここからでもOK

名前:
メール:
件名:
本文:

最近のコメント

最近のトラックバック

フリーエリア

にほんブログ村 科学ブログ 技術・工学へ
にほんブログ村

RSSフィード

リンク

このブログをリンクに追加する