> 途中までの解析結果の変位対時間曲線を調べてみましたら
> 計算終了直前あたりから変位（荷重ステップ）が
> 塑性領域にはいります。つまり塑性領域にちょっと入ってから
> エラーが出るんですよね～。
> もしかしたらこれが原因じゃないかなと思っているのですが。。

材料定義は適切でしょうか？
解析中に材料が完全塑性(剛性ゼロ)になっていたり、硬化係数が応力値より
小さくなって、つい合いが得られなくなった可能性もあるかもしれません。
もしそうなら、応力-歪み曲線を拡張するなどで解決(計算)できるかと。。。
見当違いだったらゴメンなさい。

> よし☆彡さん
> > たぶん,対象となる製品などがあると思いますから,その事を踏まえて
> > 何をもってMARC＞NASTRANとするかが大切かもね!
> > 失敗すると場合によっては逆かも,,(笑)
>
> 御教授ありがとうございます。
> うちはいわゆる外注でして、現時点では対象となる製品がないため、
> 安易に熱解析でもMARC＞NASTRANが言えればと考えていました。
> 不勉強で恥ずかしい限りです。（汗）

蛇足ですが、私は勉強不足ですが、
今現在では、接触境界(移動境界）は、M***に分があると考えてます。
但し、得られる解の精度ということではないですので、、、(笑）
いろんなケーススタディがありそうで、困りますね。
久々にカードやプロパティといった議論があったので新鮮でした。(涼）

よし☆彡さん,ハッピーさんご教授ありがとうございます。
盆明けにさっそくやってみたいと思います。
でわー。

ハッピーさん
> FEMAPだと、
> 境界条件の設定は
> ・model/load/elementalまたはon surfaceで、
> 　Convectionを選んでCoefficient（熱伝達係数）とTemperature（雰囲気温度）を設定
> 書き出すときに
> ・file/export/analysis modelで、nastranのsteady state heat transferを選べばOK。

FEMAPでも設定できたんですね。（汗）
ご丁寧にご指導頂きありがとうございました。
休み明けにさっそくやってみたいと思います。
重ね重ねありがとうございました！

> そしてその伝達条件(熱伝達係数)をCHBDY等の面要素に対する物性としてを与える！
> との解釈でよろしいのでしょうか？
by 木下さん

良いと思いますよ。

FEMAPだと、
境界条件の設定は
・model/load/elementalまたはon surfaceで、
　Convectionを選んでCoefficient（熱伝達係数）とTemperature（雰囲気温度）を設定
書き出すときに
・file/export/analysis modelで、nastranのsteady state heat transferを選べばOK。
すると、
CHBDYG 2 AREA4 +CB 2
+CB 2 2 1 3 4
PCONV 2 2 0 0.
MAT4 2 1000.
SPOINT 100001
SPC 1 100001 1 200.
CONV 2 2 100001
と、一つの条件を指定するだけでもたくさんのデータが出てきます。
SPOINTがダミーのスカラーポイントで、この温度をSPCで与えたり、Tempで与えたり。

ハッピーさん
> Natranは、境界条件を表すための「面要素(CHBDY、CHBDYG、QHBDY)」を
> ソリッド要素面に張り付け、この面要素に対する、物性として熱伝達係数や、
> 節点拘束として雰囲気温度を与えます。
> Nastranでは、データのヒエラルキー
> 　Cxxx（要素）-->Pxxx（Property）-->MATxxx（材料特性）
> 　　　　　　　-->Grid（節点）-->SPC、TEMP（拘束条件）
> を守り、何とCHBDY要素に対応するMATカードに熱伝達係数が入ります。
> 熱伝達条件では、さらにPxxxで参照されるCONVカードが加わって、
> これで雰囲気温度を指定するダミーの節点を指定し、この節点に対し、
> SPCで温度拘束（雰囲気温度）を設定します。

御教授ありがとうございます。
ほんとに何のこっちゃわからなかったものが、(汗）
おおよそ理解できました。
「ダミーの節点」が発熱源みたいなもので、
それを熱源側の任意の位置に作って雰囲気温度を与える。
そしてその伝達条件(熱伝達係数)をCHBDY等の面要素に対する物性としてを与える！
との解釈でよろしいのでしょうか？


> これは、手作業で組めるものではなく、Patranなどのプリプロセッサーで
> 作成して、出来上がったデータとマニュアルを見比べて驚くことになると思います。
> よし☆彡さん、おっしゃるように添付の例題集もじっくり分析しましょう。

FEMAP(NASTRAN for Windows)でデータを組もうとしたんですが、
設定方法がわからなかったのでExcelでつくろうとしたものの
カードが理解できず途方に暮れてました。(汗）
patranは少ししか触ったことがないのですが、
例題集含め、盆明けに格闘してみます。


よし☆彡さん
> たぶん,対象となる製品などがあると思いますから,その事を踏まえて
> 何をもってMARC＞NASTRANとするかが大切かもね!
> 失敗すると場合によっては逆かも,,(笑)

御教授ありがとうございます。
うちはいわゆる外注でして、現時点では対象となる製品がないため、
安易に熱解析でもMARC＞NASTRANが言えればと考えていました。
不勉強で恥ずかしい限りです。（汗）

> 一様な板圧（温度依存物性）の両側に異なる(非)定常温度境界条件(雰囲気温度、熱伝達係数)
> を与える程度の簡単解析でその差（MARC＞NASTRAN）が示せればと思うのですが、

たぶん,対象となる製品などがあると思いますから,その事を踏まえて
何をもってMARC＞NASTRANとするかが大切かもね!　
失敗すると場合によっては逆かも,,(笑)
方法はNASTRANの入門例題集によりどりみどり載ってると思います.

> （温度条件が急変する際にはNASTRANは精度が落ちると聞いていますが、
> 実際のところどうなのでしょうか？）

私は初耳です.(笑)

>そして、条件として増分制御方法は荷重増分法、
>繰り返し計算方法はニュートンラプソン法で行っています。

接線剛性から算出して仮に予測した結果と実際のモデルの結果を
比較して許容値以下であれば,仮に予測した結果が等しいと
いうことになります.
　このとき比較される対象が幾つかありますが,その比較対象の
何かが満足しなかった場合,ニュートン法の収束圏外と判断し,
最小増分を下げてきます.ですから,メッセージにどの節点,要素
が原因で下げたか残っていると思いますよ.
　さらに,ニュートン法の場合計算毎に2乗の割合で近似していく
手法ですから,メッセージに出た原因の節点がどのように近似して
いくかによってその対応方法も違ってきます.

難しければ,以下の方法を試してみたら良いかも!

1.　間違った要素の選定をしてないか調べる.(幾何学的非線形に対応した要素か?)
2.　修正ニュートン法から弧長法に変更するか,または荷重制御から変位制御に変更.
3.　塑性された付近のメッシュを粗くするか,不規則なメッシュにする.

> そのまま下げても,1×10＾-19のように
> 不必要に下げてる場合は時間ばかりロスし
> 無理だと思います.

そうなんですよね。1×10＾-19という極端に
さげる方法はやたら時間がたってしまいます。
こういったことは知っていたのですが
自分の今の知識ではこうするしか方法がなかったんですよね。

> 　その場合は何が原因で増分が下がったか,
> その原因を調べ,対策します.

途中までの解析結果の変位対時間曲線を調べてみましたら
計算終了直前あたりから変位（荷重ステップ）が
塑性領域にはいります。つまり塑性領域にちょっと入ってから
エラーが出るんですよね～。
もしかしたらこれが原因じゃないかなと思っているのですが。。

> しかし,その対策手法が統一されてるわけでない
> ので,実はノウハウということになると思います.

やっぱりこういった事はノウハウの領域に入ってしまうんですね。
やっぱしノウハウというのはその人が苦労して手に入れるわけですから
試行錯誤の繰り返しだと思います。今の自分はどこの部分を試行錯誤していけば
分からなかったものですから困っていました。
そして、条件として増分制御方法は荷重増分法、
繰り返し計算方法はニュートンラプソン法で行っています。

>各ステップでの最小増分量を一番小さな値を設定して
>うまくまわらなかった場合、解析は無理なんでしょうか？
>なんとか計算させたいのですが。
そのまま下げても,1×10＾-19のように
不必要に下げてる場合は時間ばかりロスし
無理だと思います.
　その場合は何が原因で増分が下がったか,
その原因を調べ,対策します.
しかし,その対策手法が統一されてるわけでない
ので,実はノウハウということになると思います.
とりあえずモデルがどうであるかは必要ないので,
計算がどのような原因(残差力が規定以下で
カットバックしたなど)とその場所を調べて
みて下さい.

> メジャーにはグループとして選ばれた要素やメジャーポイント
> （１点）が選択可能と思います。しかし、界面の応力を目的
> 関数にする場合等、いくつかのメジャーポイント（もしくは
> グループ化されている幾つかのポイント）を目的関数にしたい
> 場合があると思いますが、MECHANICAでは、これは無理なんで
> しょうか？

by promechanica初心者 さん

私のわかる範囲でお答えします。
MECHANICAでは指定した範囲での最大値か、指定したある一点の値しか
目的関数にできません。
ご質問は、複数の点の中の最大値を目的関数にしたいとのことだと
思いますが、複数の点を指定することはできません。
この場合はご存知のように、対象とする要素を選択して指定するか
中心点と半径を入れて、その範囲内での最大値をとるようにするしか
できません。
ちなみに要素を指定してその中の最大値をとる方法は「独立モード」でしか
できません。
上記２つのどちらかを使って代替させることができるといいのですが・・・。

>iSIGHTでは、探索的手法なら焼きなまし法(SA:Simulated Annealing)
>を勧めているようです。

そうですか,やっぱり1回の計算時間がネックになってる感じがしますね.
詳しく書いていただいてありがとうございます.

余談ですが,一時期GA(遺伝的アルゴリズム)がものすごいブームで
様々な問題に応用されてきましたが,GAが流行るかはやらないかのころ,
生物の進化は,そんな単純なものでないと,激怒された先生がおられました.
もちろん,その方は生物の進化に思い入れがあったのだと思いますが,,,
GAも結構パラメータに依存してしまいますからねぇ
本当に真解だけにやっかい!　　
　すみません,つまらない余談でした.