CAE掲示板 | 挑戦する製造業のために/ＮＣネットワーク

2006-12-19 10:31 投稿者： fuyo リンク：

削除キー (？)

> > 単純な梁に、振幅aの単振動を与えたときの周波数特性を知りたい時
> >
> > ・周波数応答解析の基礎励起を使用すれば良いんでしょうか？
>
> 「基礎励起」もしくは「加速度荷重」にします。
>
> > ・周波数依存の関数は、a*sin(2*pi*f)で良いんでしょうか？
>
> 違います。
> 横軸を「周波数」、縦軸を「加速度」とした時の入力の周波数特性
> ですから、全周波数体に関して一定の加速度をとりたいのであれば、
> a*(2*PI*f)^2と言う式になります。
>
> > ・周波寸依存のデフォルトの一様てなんでしょう？
>
> 全周波数帯に渡って単位加速度入力と言うことになります。
>
> > ・方向は、何を入れれば良いんでしょう？？
>
> 拘束部に加速度をかける方向です。
> ベクトルですので、合成されて所定の方向を向くようになればいい
> はずです。
>
> また、「次のアイテムに相対的」は「サポート」にしておいた方が
> いいと思います。

ありがとうございます。

No.69484 Re: 9349 バッチ処理

2006-12-18 23:14 投稿者：梅里リンク：

削除キー (？)

> ABAQUS/StandardVer6.5-6では
> nlge abq656 job=ジョブ名(もしくはn1ge?)でバッチ処理が可能らしい
> のですが、コマンドを入力しても「内部コマンドまたは外部コマンド、
> 操作可能なプログラムまたはバッチファイルとして認識されていません」
> と表示されます
ABAQUSは、全く知らないですが、上記のメッセージは、Windowsのメッセージです。Pathがとおせば良いだけのように思います。
インストールは、正常に終了しているのでしょうか。インストールが正常終了しているのでしたら、もう一度インストールマニュアルを読み返したほうが良いでしょう。

No.69483 Re: 9360 メカニカで周波数応答解析

2006-12-18 22:21 投稿者： imada リンク：

削除キー (？)

> 単純な梁に、振幅aの単振動を与えたときの周波数特性を知りたい時
>
> ・周波数応答解析の基礎励起を使用すれば良いんでしょうか？

「基礎励起」もしくは「加速度荷重」にします。

> ・周波数依存の関数は、a*sin(2*pi*f)で良いんでしょうか？

違います。
横軸を「周波数」、縦軸を「加速度」とした時の入力の周波数特性
ですから、全周波数体に関して一定の加速度をとりたいのであれば、
a*(2*PI*f)^2と言う式になります。

> ・周波寸依存のデフォルトの一様てなんでしょう？

全周波数帯に渡って単位加速度入力と言うことになります。

> ・方向は、何を入れれば良いんでしょう？？

拘束部に加速度をかける方向です。
ベクトルですので、合成されて所定の方向を向くようになればいい
はずです。

また、「次のアイテムに相対的」は「サポート」にしておいた方が
いいと思います。

No.69481 最近の学会・これからの学会

2006-12-18 21:46 投稿者：ハッピーリンク：

削除キー (？)

私は機械学会に入っていますが、最近の講演会は企業からの参加者が少ないという話を
聞きました。（私自身も２、３年ほどご無沙汰していますが）
CAE関連では計算工学会や設計工学部門とかもありますが、どんな状況なのでしょうね。
一方で、ベンダーさんのユーザー会やCAE懇話会など企業よりのところは賑やかのようです。
とは言え、エンジニアの卵を育てる場であり、また基盤的な研究開発を行う場である大学・学会
がもっと盛り上がって欲しいですよね。ちょっと堅い話でスミマセン。

No.69475 メカニカで周波数応答解析

2006-12-18 18:30 投稿者： fuyo リンク：

削除キー (？)

単純な梁に、振幅aの単振動を与えたときの周波数特性を知りたい時

・周波数応答解析の基礎励起を使用すれば良いんでしょうか？
・周波数依存の関数は、a*sin(2*pi*f)で良いんでしょうか？
・周波寸依存のデフォルトの一様てなんでしょう？
・方向は、何を入れれば良いんでしょう？？

No.69441 Re: 9354 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-17 14:23 投稿者：ハッピーリンク：

削除キー (？)

> YiやδXの値そのままを判定するのが絶対判定で、一方、変化の程度で判定するのが相対判定。
> つまり、Yi+1/Yi＜判定値、δXi+1/δXi＜判定値
また自己レスです。

相対判定は、上の記述は誤り。変化の程度が小さくなったことを判定するので、
(Yi－Yi+1)/Yi＜判定値、(δXi－δXi+1)/δXi＜判定値
でした。すみません。
何を判定基準に使うか（残差荷重、変位の変化）、絶対判定か相対判定か、そして基準値をいくらにするか
で、打ち切るXi、即ち解が変わってきます。

No.69437 Re: 9354 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-17 10:35 投稿者：ハッピーリンク：

削除キー (？)

> 紙にY＝X^2のグラフを、X＝０から１の範囲で描いて、X軸上にX0、グラフ上にY0点をおいて
> X0からY0に直線を引く、これがスタートラインです。
> ここから（Xi、Yi）を動かしていき、Yi≒０となるXiを探すわけです。

自己レスです。
Y＝X^2のグラフは薄く書いて下さい。最初からグラフがあってその上を動かすのではなく、
あくまで、Xiに対するYiを、Xi^2の式で逐次求め、結果として（Xi、Yi）の点群が得られます。
応力解析で言えば、結果として荷重変位曲線が得られます。

No.69435 Re: 9355 CAE技術者認定試験

2006-12-17 09:14 投稿者：ハッピーリンク：

削除キー (？)

> テキストだけの勉強では受かりそうもありません。
> みなさんどうクリアされましたか？
by むうううさん

右上に小さな「検索」欄があります。「認定試験」で検索すると諸先輩の
ご経験談がいくつか見つかると思います。受験勉強の方法については少ないようですが。

No.69431 CAE技術者認定試験

2006-12-17 01:18 投稿者：むうううリンク：

削除キー (？)

固体力学分野２級を受けました。
はっきりいって（泣）...です。
テキストだけの勉強では受かりそうもありません。
みなさんどうクリアされましたか？
熱流体もテキストだけでは難しいものなんでしょうか。

No.69429 Re: 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-17 00:47 投稿者：ハッピーリンク：

削除キー (？)

> Newton-Raphson法を文字で分かりやすく説明するのは、かなりの表現力が必要かも．．．
> 式なら→http://kiso.dj.kit.ac.jp/~nakamori/class/Numerical/html/chap2/Newton.html　　

無理があるかもしれませんが．．．
一番単純な非線形方程式式として、X^2(2乗)＝０をニュートンラフソン法で解いてみると。
答えはいうまでもなく、X＝０です。
収束計算を行うには、まず初期値が必要です。仮にX0＝１とします。
紙にY＝X^2のグラフを、X＝０から１の範囲で描いて、X軸上にX0、グラフ上にY0点をおいて
X0からY0に直線を引く、これがスタートラインです。
ここから（Xi、Yi）を動かしていき、Yi≒０となるXiを探すわけです。
準備として、微分して接線の傾きの式＝2Xを導いておきます。

Step0：（X0、Y0）＝（１、１）　接線の傾き：２
　点（１、１）から傾き２の接線を引いてX軸との交点を求めると、X＝1/2となり、これが初めの第一歩。
Step1：（X1、Y1）＝（1/2、1/4）　　接線の傾き：１
　点（1/2、1/4）から傾き１の接線を引いてX軸との交点を求めると、X＝1/4となり、これが２歩目。
Step2：（X2、Y2）＝（1/4、1/16）　　接線の傾き：1/2
　点（1/4、1/16）から傾き1/2の接線を引いてX軸との交点を求めると、X＝1/8となり、これが３歩目。
このように、のこぎりのように描いていくと、Yｉが１、1/4、1/16、1/64どんどん小さくなっていきます。
一方、Xについて見れば、変化の幅δXが、
　　　X0→X1：　δX＝1/2
　　　X1→X2：　δX＝1/4
　　　X2→X3：　δX＝1/8　
とやはり小さくなっていきます。
数値解析ですので、δX、Yiがぴったりゼロになることはなく、どこかで「収束した」と打ち切る判定をする
必要があります。そこで、「Yiが十分ゼロに近づいた」で判定するのが残差荷重判定、
「δXが十分ゼロに近づいた」で判定するのが変位判定。
YiやδXの値そのままを判定するのが絶対判定で、一方、変化の程度で判定するのが相対判定。
つまり、Yi+1/Yi＜判定値、δXi+1/δXi＜判定値
応力解析に置き換えて言えば、Xは変位、Yが荷重。傾きは剛性マトリックスKです。
また、接線を引いて交点を求めるのは、荷重に対する変位増分を求めていること、
また、X＝？に対する、Y＝X^2を求めているのは、等価節点力と荷重との差＝残差。
と言うか、いま描いたグラフを180度回転させると非線形の荷重変位曲線に見えると思います。
どうでしょうか。

No.69409 Re: 9351 9338 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-15 12:37 投稿者：よし☆彡リンク：

削除キー (？)

> 世の中そこまで来てますか。うちにはマダマダ抵抗勢力がいたりします。
> 井の中の蛙にならぬよう気合を入れないと．．．
by ハッピーさん

状況はハッピーさんと一緒だと思いますよ,ある領域は確立して,
新規領域は抵抗という構図です.逆に言えばここ２，３年新規領域の
開発を行ってないということで,誇れる状態ではないんです.(泣)

でも,そんな状態でも皆さんの話題がわかるのはやっぱり,基礎を
勉強したからかもしれません.理論的に想像できますから,,,(笑)
初心者の方には特に,基礎をやっておくことをお勧めしますね.

No.69396 Re: 9338 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-14 22:13 投稿者：ハッピーリンク：

削除キー (？)

> プレゼに載せる実測データ収集に困ったり,逆に実験データと同時に
> 解析データをプレゼンで使うと,”なんで,解析したのに実験するの?”
> なんて質問になったるするんですよねぇ～
> 解析が珍しいなんてのは,過去の時代だし,ベンダーの開発スピードも
> 早いので,”普通じゃん”と言われないための努力も大変です.ホント!
by よし☆彡さん

世の中そこまで来てますか。うちにはマダマダ抵抗勢力がいたりします。
井の中の蛙にならぬよう気合を入れないと．．．

No.69395 Re: 9348 9337 熱応力解析の収束判定について

2006-12-14 22:07 投稿者：ハッピーリンク：

削除キー (？)

> >収束設定のルール化進めている次第です。
> >そこで、熱応力時は変位判定が経験上良かったので、それをルール化する
> >上で、理論的な理屈付けがほしかった次第です。
> by temmaさん
> 　収束判定のルール化は確かに重要で難しいですね．
by Peroさん

今年の品質工学研究発表大会で、シェル構造の溶接部のモデル化について解析精度を安定化する
モデル化（制御因子）の発表があって注目を浴びていました。自動車では、そこら中に溶接が
あり、いちいち細かな吟味ができないので、ルール化をしたいというのが動機だったと思います。
収束条件のような経験を要する問題にも適用できると面白いですね。
パラメーター設計で、設計の要求に応じてQCDを如何様にでも調整できるかも。

No.69394 バッチ処理

2006-12-14 20:33 投稿者：くるりリンク：

削除キー (？)

ABAQUS/Standardでバッチ処理を行いたいのですが、デフォルトでは
使用できないのでしょうか。

ABAQUS/StandardVer6.5-6では
nlge abq656 job=ジョブ名(もしくはn1ge?)でバッチ処理が可能らしい
のですが、コマンドを入力しても「内部コマンドまたは外部コマンド、
操作可能なプログラムまたはバッチファイルとして認識されていません」
と表示されます

バッチ処理を可能にするには何らかの作業を
サブルーチン？に足さないと行えないんでしょうか

No.69369 Re: 9337 9336 9332 9331 9329 熱応力解析の収束判定について

2006-12-14 01:29 投稿者： Pero リンク：

削除キー (？)

>収束設定のルール化進めている次第です。
>そこで、熱応力時は変位判定が経験上良かったので、それをルール化する
>上で、理論的な理屈付けがほしかった次第です。
by temmaさん

　収束判定のルール化は確かに重要で難しいですね．
ルール化に際し，思ったことがあるので，ちょっと書き込みしたいと思います．

例として簡単な塑性加工問題を取り上げます．

まず，少々乱暴な表現ですが，この場合，変位判定については，

Incrementにおける最大節点変位＝剛体の変位＝工具の1Incあたりの移動量・・・A
Iterationにおける修正変位＝収束計算中に修正された変位の最大値・・・ B

となります．

　考えてる問題は，板厚10mmで，板厚方向分割が10とします．
すると，板厚方向の1要素の長さは1mmですね．これをCとします．

さて次に，1Incあたりの工具の変位量を考えます．

　案①として，Aを10mm，相対しきい値A/B＝0.1とした場合，Bは10×0.1=1mm
となりますね．
今度は案②として，Aを0.1mm，相対しきい値A/B＝0.1とした場合，Bは0.1×0.1=0.01mm
となります．

　このときのBが感覚的に許容できるかどうかです．

つまり，板厚方向の1要素の長さCが1mmなのですから，
案①では，Bが1mmですので，BとCが同等なので，精度的には，どうかな？と判断できるかと
思います．案②の方が，妥当ですね．

以上，うまく言えないのですが，要は，相対しきい値A/Bも重要なのですが，Aの値とBの値の
絶対値に対しても，妥当なのかどうか把握する必要があるのではないかと思います．

その上で，時間がかかるような問題なのであれば，それは仕方がないことであり，
設計部隊の方に，理解してもらう努力も必要かと思います．

　ご参考まで．

No.69366 Re: 9342 9341 9339 メカニカで固有値

2006-12-14 00:24 投稿者：ハッピーリンク：

削除キー (？)

> > > 　メカニカで水中の物体の固有値解析て、可能なんでしょうか？
> 液中の薄板状の固有値を知りたいのですが、
by fuyoさん

昔、同僚がタンクのスロッシングの解析を行っていましたが、付加質量では難しいでしょうか？
「付加質量　振動解析」あたりで検索してはどうでしょう。
最初にヒットしたのが→　http://www.east.jasnaoe.or.jp/snaj/tm/0103/ronbun31.html

No.69364 Re: 9337 解析のQCD

2006-12-14 00:17 投稿者：ハッピーリンク：

削除キー (？)

> 私の仕事には解析Quality 向上活動というのが関係します。
> 私の会社では、設計者向けに受託解析を行っていますが、
> 最近の短納期には、本当に悩まされます。
> もっとデリバリーを！！設計の人の耳にはあまり届きません。。。。
by temmaさん

私も前は受託解析を商売にしていましたので少しはお気持分かります。
どこかで読みましたが、昔は計算機コストが非常に高く、実際、費用を見積もる時は
計算機使用料が人件費と肩を並べていました。解析技術者は高級な装置を操るある種の
ステータスがあったかも知れません。ところが計算機のコストパフォーマンスの飛躍的
向上で人件費が際立つようになると共に解析の普及もあって解析技術者の立場は苦しく
なってきた．．．と。
　私は解析に関しては、QCDでなくQDCと思っています。やはり開発はTime Is Moneyですから
デリバリーの短縮は避けがたい（temmaさんゴメン）。
　じゃあ計算機を数段ランクアップするのは如何でしょう？当然、その費用は設計に請求する。
「要求されるQとDを満足するにはCアップが必要」と、設計者に解析コストの意識を持たせる
のも意味があるように思います。
　（ろくに考えずに解析に回す、なんて不埒なケースを回避できる効果もあるかも）
　あと設計者を巻き込む。3DCADは当たり前になっているかも知れませんが、もっと早い時点で
解析チームと連携する「コンカレント」なんていうのも。

No.69358 Re: 9342 9341 9339 メカニカで固有値

2006-12-13 21:57 投稿者： imada リンク：

削除キー (？)

> 液中の薄板状の固有値を知りたいのですが、
> 周辺液体のダンピング効果て問題になるのでしょうか？
>
> 　そもそも、静解析では意味がないのでしょうか？？

そんな気がしますね。
水の流体としての挙動を考慮して解析するとしたら
メカニカの線形解析では無理だと思います。
水中での固有値がどの程度変化するのものなかは良く知らないのですが・・・。

No.69354 Re: 9341 9339 メカニカで固有値

2006-12-13 20:17 投稿者： fuyo リンク：

削除キー (？)

> > 　メカニカで水中の物体の固有値解析て、可能なんでしょうか？
>
> 液体との連成解析と言うことですか？
> 具体的にどんな問題なのでしょう？？？
>

液中の薄板状の固有値を知りたいのですが、
周辺液体のダンピング効果て問題になるのでしょうか？

　そもそも、静解析では意味がないのでしょうか？？

No.69347 Re: 9339 メカニカで固有値

2006-12-13 18:35 投稿者： imada リンク：

削除キー (？)

> 　メカニカで水中の物体の固有値解析て、可能なんでしょうか？

液体との連成解析と言うことですか？
具体的にどんな問題なのでしょう？？？