コンピュータによる構造解析結果を手計算により実証したいです。
問題は単純に筒にかかる応力の問題です。
添付図のような筒に解析ソフト(ANSYS workbench)によりランダム振動を与えると図で赤に色を付けた部分に最大応力が発生しました。そこでその最大応力値の妥当性を手計算(曲げ応力計算)により実証したいのです。その場合以下の方法であっているのでしょうか?
固定部分以外の重さ:W=0.120[kg]
固定部以外の重心距離:X=0.020[m] (図参照)
加速度:G=10[Grms]
すると固定面のそばにかかるモーメントMは
M=W*X*G
=0.120*0.020*(10*9.8)
=0.2352[N・m]
また筒の断面係数Zは
Z=π(D^4-d^4)/32D (図参照)
=0.000000231[m3]
したがって相当応力σは
σ=M/Z
=0.2352/0.000000231
=1018260[Pa]
≒1.0[MPa]
ここまで、筒にかかる応力の計算として一般的にあっていますでしょうか?
またANSYSによる計算結果はもっと大きい値を示しているので単純にランダム振動の二乗平均の10Grmsを使用するのではなく、違ったファクターを加えるのではないかと思っています。
以上、どなたかお詳しい方がおりましたらどうぞ宜しくお願いします。
No.1ベストアンサー
- 回答日時:
回答するための情報が不足しています。
下記の点について、明確にしてください。
ランダム振動を与えたとありますが、
a.本当にランダム波形を与えて、振動解析を行ったのですか?
b.それとも、単に、10G相当の体積力を円筒軸に垂直な方向に与えて、静解析を行っただけなのですか?
ANSYSでは、どのようにして解析したのですか?
a.マップトメッシュを使って、きれいに分割したのですか?
b.それとも、Workbenchに勝手にメッシュ切りさせたのですか?
この質問は、次の質問に関係するだけであって、あまり重要ではありません。
最大応力発生位置での円筒の肉厚は、半径方向に、何分割ぐらい、されていますか?
a.1分割?(=半径方向には分割されていない)
b.数分割?
c.2分割程度?
上記の質問に対する回答は、すべてb.だと予想しています。
もしそれが当たっているなら、あなたの手計算に比べ、FEM解が大きくなったのは、応力集中のせいで、極めて正常です。
ある領域を完全固定すると、固定端の所には、無限大の応力が発生します。(応力が無限大になるような点を、特異点と呼びます。特殊なケースでは、変位も無限大になって、これも特異点になります。)
の質問に対する分割数が大きければ大きいほど、FEM解は、この特異点でも応力を、何とか表示しようとして、大きい値の応力値を計算して出してきます。
FEMの計算結果が、手計算と一致するのは、(2)の質問で、マップトメッシュを使い、半径方向を1分割した場合のみです。
この場合に限り、FEMも特異点の応力(に近い値)を計算することはできません。
一方で、あなたが利用した梁の計算式では、特異点での応力が発生することへの配慮はされていません。
梁の応力の公式では、断面に対し、線形変化する応力しか考慮されていないからです。
梁の応力をFEMで計算し、梁の式で得られる解と比較しようとする人の多くは、この特異点の現象が理解できていないために、つまずくのです。
要するに、両者を比較することには、あまり意味がないのです。
なお、あなたの手計算の式は、合っていると思います。
この回答への補足
ご回答ありがとうございます。質問にも足りないところがありましてすみません。
回答者様からのご質問への回答ですが、実は(たぶん)全てa.です。
>ランダム振動を与えたとありますが、
>a.本当にランダム波形を与えて、振動解析を行ったのですか?
>b.それとも、単に、10G相当の体積力を円筒軸に垂直な方向に与えて、静解析を行っただけなのですか?
a.です。質問では省略してしまったのですが、周波数に応じて加速度が変わるランダム振動の条件で解析を行っています。Gの条件はパワースペクトル密度PSD[G^2/Hz]で表わされています。(例:0-50Hzは3G、50-100Hzは6G・・・といった感じで2乗平均?が10Grmsのようです。)
>ANSYSでは、どのようにして解析したのですか?
>a.マップトメッシュを使って、きれいに分割したのですか?
>b.それとも、Workbenchに勝手にメッシュ切りさせたのですか?
>この質問は、次の質問に関係するだけであって、あまり重要ではありません。
どちらかというとa.だと思います。トライアングルメッシュ?を使用し、全体を1.0mmのメッシュで切って解析しています。(ただし今考えてみると筒の肉厚が0.5mmなのでメッシュサイズが大きいですか?)
>最大応力発生位置での円筒の肉厚は、半径方向に、何分割ぐらい、されていますか?
>a.1分割?(=半径方向には分割されていない)
>b.数分割?
>c.2分割程度?
たぶんa.です。すみません。まだ解析を学び始めたばかりで分からないことが多いのですが、構造解析では肉厚(メッシュで?)も分割して解析を行うのが通常ですか?肉厚の分割に関しては何も設定した覚えがありません。
ちなみに材料定義?は線形です。線形解析ですかね?その場合必ずa.ですか?
また、特異点に関してはどこかで聞いた覚えもあります。いくらメッシュを細かくしていっても無限に応力値が大きくなるような例ですよね?
静解析ならまだしもランダム振動をかけた解析を手計算と比較するのはやはり厳しいのでしょうか。
宜しくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# C言語 3 2022/10/04 15:07
- 物理学 風車が回転する理由 5 2023/05/08 08:03
- 建築学 梁の計算で手こずってます。建築系の事に詳しい人、誰か助けてください 僕は工学系の学生なんですが、梁に 1 2022/11/15 12:53
- 化学 [例題①]プロパンのOHとの反応速度定数を k=1.0×10⁻¹²cm³s⁻¹、 [OH]=1.0× 0 2023/05/28 13:12
- その他(Microsoft Office) 従業員増減対応で当番種類の増減対応な当番表 21 2022/07/19 07:30
- 統計学 回帰分析検証不一致 5 2022/05/04 20:48
- 宇宙科学・天文学・天気 AIが答えた方程式 1 2023/02/20 00:12
- 生物学 【生命科学】ヒトが1日に消費するATP量?(精度を変えて再計算) 3 2022/10/07 18:48
- 物理学 流体力学の勉強方法 2 2023/04/25 23:18
- 発達障害・ダウン症・自閉症 中学の時にIQ82の境界知能と診断されました。 今の私も、やはり境界知能でしょうか? そしてこれは、 3 2023/02/19 00:37
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
降伏点が明確でない場合の降伏...
-
FEMのミーゼス応力と主応力の違い
-
応力と凸凹
-
真空容器の強度計算方法を教え...
-
応力、垂直抗力について
-
JIS K 6301ゴムの物性計算の...
-
ゴムの静的せん断弾性率とは?
-
銅の応力‐ひずみ線図
-
物理の問題です<(_ _)>助けてく...
-
法線応力差についてわかる方い...
-
薄肉円筒の軸応力について
-
0.2%耐力とは?(弾性係数,応力...
-
片持ち梁の破壊荷重について
-
ミーゼス相当ひずみについて
-
引っ張り変形におけるくびれの...
-
弾性力学で
-
降伏応力、弾性率とガラス転移...
-
斜面の垂直応力とせん断応力に...
-
3点集中荷重の最大曲げ応力の計...
-
等方性材料の構成式について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ゴムの静的せん断弾性率とは?
-
FEMのミーゼス応力と主応力の違い
-
応力と凸凹
-
真空容器の強度計算方法を教え...
-
銅の応力‐ひずみ線図
-
引張試験の負荷速度
-
3点集中荷重の最大曲げ応力の計...
-
吊り具の強度計算について教え...
-
材料のネッキングは何故おこる?
-
0.2%耐力とは?(弾性係数,応力...
-
『構造粘性』とは何でしょうか...
-
ポアソン比と張力の関係!?
-
木材の曲げ弾性率と曲げ強度に...
-
段付き棒(梁)の曲げについて
-
応力勾配とはどう言う意味ですか?
-
法線応力差についてわかる方い...
-
樹脂材料の曲げ弾性率について
-
ゴム材にかかる応力計算
-
曲げ弾性率から引張弾性率を求...
-
薄肉円筒の耐圧計算に関して
おすすめ情報