千代岡の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
2024年4月からは[[千代岡の修論日誌:https://www.str.ce.akit...
3月までに書いたことはそのままこのページに残しておきますが...
研究内容：様々な外力が作用する斜張橋の動的応答解析
*ちょっとした春課題 [#u272b919]
model300のケーブル7本モデルのC1破断時の主桁中央の鉛直変位...
減衰率の設定法は[[秋山さんの修論日誌:https://www.str.ce.a...
・減衰0%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰0.5%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰1%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰2%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰5%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰10%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
基本的に減衰率の値以外の設定条件は先輩方が設定したものと...
減衰率5,10%だと主桁中央変位の最もたわむ値が変わっているが...
ステップ条件を複数条件→荷重時間に対する初期比率を0.01にや...
先生に結果をお見せしたところ,問題なさそうとおっしゃていた...
*データ類の置き場所 [#uc29fb98]
R3年度バックアップの根本さんのお名前が書いてある外付けHDD...
ファイル名のDは死荷重,Cは集中荷重を示しています.また,どの...
*課題 [#e9b3591e]
自分の研究成果の課題点（院でやること）をまとめておく。
・実験的検討を行う。解析結果とSalome_Mecaで出したデータを...
・橋モデルをSalome_Mecaで作成して,今回得られた結果とどう...
*卒論内容 [#a313a74b]
研究室で編み出した連鎖崩壊のメソッドの妥当性を確認する。...
しかし、その妥当性の評価を行なっていなかったので、その評...
Salome-mecaで行う理由としては...~
①ビーム要素を用いた動的応答解析の手法ができていないのでそ...
②2種類以上の解析ソフトで解析できるようにした上で、同じよ...
卒論としての目標は、簡易モデルでのDAFを用いた連鎖崩壊に関...
DAFとは？~
動的増幅率を表す指標。式等は[[こちら:https://www.str.ce.a...
*今後の予定 [#b61dcef6]
3/22 卒業式~
3/23 一旦実家に帰省~
4月以降については修論日誌に書く。
*やることリスト [#naf093be]
大学卒業までにやること
まずは東北支部の発表に備える！~
・本論作成（3/15まで）~
・引継ぎ関連（先輩方にお願いしないと...）~
・以前、Salomeで作成した橋モデル作成の続き~
・wiki整理~
・謝恩会準備（研究とは関係ないが...）~
卒論まで残り1週間 残りの期間でやること
・2回目練習の動画を見て,反省会（メモしきれなかったところ...
・スライドの修正~
・概要は卒論用に手直ししたので、一度指導教員にお見せする~
・ひたすら練習.個人練習を含め,15回くらいはやりたい~
2024年 1,2月にやること~
・0.01sから短めにtimestepを切ると途中で張力が落ちる問題を...
・東北支部、卒論に向けての概要とスライド作成~
・(研究とは関係ないけど...)謝恩会の準備
ラストスパート！！自分のできる限りのことは尽くす！
2023年 12月にやること~
・簡易モデルで得られる張力の値がMarcと大きく異なる謎を解...
・橋モデル:model300をsalomeで作成する。(橋はすでにgeometr...
・PRE_EPSIを使ったモーダル解析(山本君と共同研究)~
・12/22に向けての中間発表の準備~
今年もいよいよ今月で終わり、これまでの成果を自信を持って...
2023年 11月にやること~
・動的解析(線形&弾塑性を含んだver.)でプレストレスを考慮し...
・model300をsalomeで作成する。~
・PRE_EPSIの解析手法の検討(山本君と共同研究)~
・ビーム張力を見れるようにする。~
2023年 10月にやること~
・簡易モデルの弾塑性周辺の設定見直し→ケーブルには何も設定...
・ピン接合の設定確認(ばねが本当に一軸方向のみに作用してい...
・簡易モデルにプレストレス入れてみる~
・張力だけを落として集中荷重のみが載荷し続ける挙動再現(片...
・Salome-mecaの簡易モデルを用いた動的応答解析のやり方を学...
・Salome知識(ソリッド絡みの知識)復習 10月に来る3年にSalom...
2023年夏休みにやること
8月は実家に戻るまではピン接合問題(Salome,calculixの両方の...
Salomeを使えるノートパソコンを借りることができたので、実...
何か進捗があれば、先生に報告することを忘れないように。~
秋田に戻ったら、中間発表に向けて今まで以上に気合を入れて...
2023年7月にやること
(1)Salome-Mecaでビーム要素の弾塑性解析をできるようにし、...
⇒ケーブルには引張しか存在しないので、それを再現するため。
(2)Salome-Mecaでピン接合をできるように考え、Marcの簡易モ...
今は上記2つを重点的に行う。今のままだと、10月の中間発表で...
(3)英語ゼミの中間発表に向けてスライド作り~
⇒これは研究背景と今取り組んでいることとこれからどうしてい...
後は、自分にとって必要だと思う力(コミュニケーション能力と...
焦ってもしょうがないが、自分のやっていることなどを人に分...
2023年6月にやること
(1)Salome-Mecaでビーム要素の弾塑性解析をできるようにする。
(2)Salome-Mecaでソリッド要素の弾塑性解析のグラフ改善（Mar...
(3)Salome-Mecaでピン接合をできるように考える。
(4)振動解析の勉強（動画のモデルの実際に作る）
後は、来月に大学院入試があるので、それに向けて準備する。
2023年5月にやること
(1)先月同様、salomeでの振動解析について学ぶため、code_ast...
(2)Marcで弾塑性解析（モデルは10×10×100mm）を行い、235MPa...
Salome-MecaとMarcの結果比較をするために、Salomeで梁要素,M...
(3)Salome-Mecaでピン接合をできるように考える。
後は、時間があるときに院試に向けて勉強や本などを読んで知...
2023年4月にやること
(1)code_asterの動画[[https://www.youtube.com/@DigvijayPat...
（字幕は英語のほうがよい。...が、話している人の声が小さい...
目標はGW明けに終わらせる！（どんなに時間がかかっても5月中...
参考資料
-[[Salome-Meca説明書:https://sites.google.com/site/codeas...
-[[code-asterの説明:https://code-aster.org/V2/doc/default...
(2)時間があるときにMarcの使い方に慣れておく。
*進捗状況 [#be3fe568]
**ゼミ進捗報告 [#l52104c3]
12/12 温度考慮した動的解析ができるようになったこととケー...
12/6 現状の問題点と打開策を話す。~
11/29 軸力絡みの話~
11/22 ビーム軸力の話~
11/15 特になし~
11/8 PRE_EPSIを使ったプレストレス再現のお話~
11/1 線形での動的解析の設定の仕方は分かった。(FORC,REACの...
10/25 線形解析ではプレストレスが導入できたこと、動的解析...
10/18 プレストレス導入の話~
10/11 プレストレスとピン接合絡みの話~
10/4 後期1回目のゼミ 報告することは特になし~
**ビーム要素の簡易モデルを用いた連鎖崩壊のメソッドの妥当...
現時点で分かっている問題点と分かったこと
-曲げモーメントの出力できない．応力もできない?軸力は可能...
→SIPO_NOEU,ELNOで曲げモーメントを見ることはできるみたい。...
ただ、REAC_NODAのDRX,DRY,DRZでモーメントを見ることはでき...
[[こちらの資料:https://code-aster.org/V2/doc/default/en/m...
-現時点でケーブルはトラス要素を用いている．おそらく圧縮力...
-非線形解析と動的解析を組み合わせた解析がまだ１回も回った...
-動的の精度の問題
（以上、[[秋山さんの修論日誌:https://www.str.ce.akita-u.a...
***簡易モデル(Marc編) [#ybe00db2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑簡易モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑桁断面
簡易モデルを作成して、梁の先端変位を簡易モデルのケーブル...
ケーブルはピン接合、梁は原点を完全固定、ケーブルとの節点...
桁の長さは1000mm,完全固定から図の上のピン接合までの距離は...
材料特性は以下に示す。~
,,ヤング率(GPa),ポアソン比,密度(ton/mm^3)
,梁,200,0.3,7.9e-9
,ケーブル,195,0.3,7.9e-9
ケーブルに掛けたプレストレスは-100℃,熱膨張係数は1.25e-5
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ケーブル破断前の先端変位グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ケーブル破断後の先端変位グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル破断後の先端変位グラフ拡大図
先端変位の理論値(片持ち梁のたわみ式から)を導出すると、15....
破断後の最大変位が理論値の1.5〜2倍程度になればいいみたい...
ちなみに、ケーブル破断後に振動したグラフができたのは、梁...
***簡易モデル(Salome-meca編) [#ga6e819b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ビーム要素で作成し、モデルの向きも合わせた。 また、メッ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑梁断面(Marcと同じ) 写真には載せてないが、ケーブル断面も...
モデルの作り方は桁とケーブルのモデルを作成し、メッシュを...
***コマンドファイル [#uc67854a]
・[[ケーブル張力を得る前の設定(DYNA_LINE.ver):http://www....
***Salome-Mecaでできること [#xc3c201b]
・動的応答解析は？~
簡易モデルではDYNA_LINE,DYNA_NON_LINEを使って解析する。(...
・弾塑性解析は?~
再現できた。しかし、降伏応力より大きめの値を設定すると収...
・プレストレスは?~
線形静的解析の場合、温度からでも導入できる。しかし、弾塑...
・途中で張力を落とす(0にする)やり方~
ケーブル破断した際、張力は0になるので、それを再現するため...
***ある時間に力(荷重)を0にする方法 [#k7a3caa6]
ケーブル破断後の挙動を再現するためにある時間でケーブルの...
本来、まず簡易モデルに荷重をかけて、そのときのケーブル張...
そこで、試しに片持ち梁(10✕10✕1000mm)を作成してやってみた。~
梁の先端に0.1mmのバネ要素の線をつけて、バネの先端(1000,0....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
やり方として、DEFI_FONCTION(関数の定義)でNOM_PARA[INST]に...
表の左側にNOM_PARAで定義した数値(今回はINST)を入れ、表の...
図は張力の設定だが、INST0のとき、0倍の荷重が載荷、INST50...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
次にDEFI_LIST_REEL(データ定義)での設定を行う。~
DEBUTは0にする。これは解析を始めるスタート地点(初期値)で...
INTERVALLEはuntil(JUSQU_Aと同じ意味。横軸の最終値)に100,P...
ちなみに、分割数は(until - DEBUT)/PAS で求められる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
これらの設定を踏まえた上で以下のような結果が得られた。~
10/25 動的解析に直し、集中荷重を100N,張力を50Nに直した結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
70s後に梁の先端が揺れる挙動ができた！~
ただ、変位だけは出せるが、何故か力を導出する設定をすると...
プログラムミスということは、どこか設定がおかしいとは思う...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
10/23 静的解析でやってしまった...~
始めの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
70sの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
100sの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
集中荷重を載荷した場所でのy方向の力の様子~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
張力を掛けた場所でのy方向の力の様子~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
梁の先端変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
集中荷重と張力の状態を見る限り、載荷の仕方は求めたいもの...
しかし、張力が0になると、その衝撃で梁が上に持ち上がり、上...
原因を探ってみたところ、載荷再現は問題なかったが、動的解...
静的解析だと、時間を考慮しない解析になり、時間に依存する...
また、挙動を見る限り、荷重を掛けすぎていると思うので、荷...
***ピン接合のやり方(Salome-meca) [#a7708b00]
ピン接合とは？~
→節点（部材と支点の節点）が自由に回転するように接合する方...
バネ要素を用いた線で再現できる。(AFFE_MODELEのDIS_TR,AFFE...
部材と部材の距離はできる限り短くすること。(下の写真のモデ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
回転が自由にできるピン接合を再現するなら、回転に対する剛...
REPEREはGLOBALでよい。
バネ要素を用いたピン接合は、本来節点を共有して接合するの...
・ばねが一軸方向に作用しているかの確認~
以下の3枚の画像には桁の先端に荷重Fy=-500N掛けている。白色...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑x方向変位に対する剛性のみを緩めた場合~
y軸負の向きに力を掛けることで、ケーブルとバネ要素の線の共...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(1000,-14.5803,0.000175109)~
ケーブル側の共有節点(1000,0.1,0)→(985.425,-14.4802,0.0001...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑y方向変位に対する剛性のみを緩めた場合~
y軸負の向きに力を掛けることによって、バネ要素の線が伸びた...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(1000,-14.5802,0.000175108)~
ケーブル側の共有節点(1000,0.1,0)→(1000,0.0953226,0.000175...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑z方向変位に対する剛性のみを緩めた場合~
z軸には影響がほとんど見られなかったため、ほぼ変形せず。(...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(999.989,-0.16365,1.965e-6)~
ケーブル側の共有節点(1000,0.1,0)→(999.984,-0.0587042,1.96...
y方向変位に対する剛性を緩めた画像を見れば分かりやすいと思...
画像を見るかぎりは一軸方向に作用していることが分かる（バ...
***プレストレス [#ld0cda14]
・PRE_EPSIを使ったプレストレス再現(線形静的解析 MECA_STAT...
AFFE_CHAR_MECA → PRE_EPSI を使ってプレストレスを再現して...
PRE_EPSIは変形の負荷を与えるコマンドで、説明書を見るとεと...
[[PRE_EPSIの詳しい説明はこちら:https://code-aster.org/V2/...
今回、ケーブルにプレストレスを適用したいので、EPX（梁のロ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
上記の設定に変えて、解析を行なったら以下の結果が得られた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑scalefactor1倍だと変化が見にくかったので、30倍で見てみた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑桁の先端変位
Marcで得た先端変位：2.49607mm~
Salomeで得た先端変位：2.49608mm~
Marcだと節点に温度応力-100℃与えてプレストレス再現していた...
変形図と先端変位のグラフを見る限り、再現性は高いと思うの...
注意点として、何回か解析を回していると、計算自体は回るが...
・温度で再現する際の問題点~
12/12追記~
CREA_CHAMPとCREA_RESUを使うことで動的解析でも温度を考慮す...
CREA_CHAMPやAFFE_VARCの設定は一切変えずにSTAT_NON_LINEを...
線形静力学を計算してくれる(MECA_STATICS)は、材料の機械的...
説明書やSalome-Mecaの掲示板なども見て、試せるものは全てや...
(現在の設定)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・やり方~
CREA_CHAMPとAFFE_VARCを用いて再現可能。~
過去の先輩方が行なっていた設定を参考にしながらやってみた。~
線形解析で桁に死荷重Fy=-0.0183456N/mmを掛けているが、桁が...
(白線が変形前、色付きが変形後)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Marcでも線形解析でのプレストレスを導入した簡易モデルのデ...
,,Marc,Salome
,桁の先端変位(mm),2.4954,2.49607
SalomeでもMarcとほぼ同じ結果が得られたのでプレストレスの...
後、ピン接合の箇所を確認したところ、変形は見られなかった。~
***軸力の見方 [#df51b2bf]
軸力はEFGE_ELNOの[N]で見る事ができることが分かった。Nは法...
見る方法だが、medファイルをParaViSで閲覧し、フィルター→me...
[[参考資料:https://code-aster.org/V2/doc/v15/en/man_u/u2/...
↑Salomeで出力される力やひずみなどについて詳しく書かれてい...
<模索中のメモ>~
節点力から求める軸力はあまりにも値が大きすぎることから、...
しかし、得られた軸力は4Nほどでケーブル破断後のモデルにそ...
[[PRE_EPSIを使って力を見てみる人が質問していた:https://co...
ただ、回答していた人が[[あるテストケース:https://code-ast...
12/12追記：残念ながら自分の使っているSalomeのversionでは...
***あるtimestepにすると張力が落ちる問題 [#ke57d6cd]
timestepを0.01より小さな値で刻むと安定状態に入る直前に張...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑初期モデルから得た張力(timestepを0.01にしている) 要素分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑99sと100sの間の様子 一瞬だけ張力が落ちているのが分かる。
原因としては温度を用いたプレストレスが問題であることが分...
ただ、PRE_EPSIは設定したひずみになるような外力を与えてい...
そこで、ケーブルを10要素に分割して、1要素（(0,1000,0)側）...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑白色のところにPRE_EPSIを適用した。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑変位も張力も温度で解析していたときとほぼ一致しつつ、張力...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ただ、個人的に気になる点として、0.1Nほどの振れ幅が見られ...
(追記)無視していいそうです。
***10mの簡易モデルについて [#h6683bd4]
1m(1000mm)の桁長で解析を行なったが、思っていた以上の衝撃...
そこで、桁長を10m(10000mm)に変えて解析を行なってみた。(断...
・問題点
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑1m簡易モデルの初期モデルは上記のような挙動が得られていた。
しかし、10m簡易モデルでは下にたわんでしまった。(120sの様子)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル1本モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル2本モデル
この状態でケーブル1本モデルの解析を行なってみたところ、以...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ただ、ケーブル2本モデルで解析を行うと,初期モデルのC1ケー...
<対策>~
桁がたわむのは問題ないらしい。まず、1本モデルの結果をまと...
2本モデルに関しては、C2のプレストレスを落として再解析。
**解析結果 [#h7cbc853]
***ケーブル要素なしモデルの説明 [#p1b7be48]
このモデルは初期モデルで得た初期張力を桁の先端に元々ケー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
2024年のデータ~
12/20~
120s後、10s,5s,1s,0.1s,0.05sの時間を掛けて張力を0に落とし...
変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Salomeの結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcの結果
10s〜0.1sまではよさそうな結果が得られているが...0.05sを見...
・問題点~
本来であれば、0.01sとも比較するつもりだったが、初期モデル...
そのため、そのデータをケーブル要素なしモデルに入れても安...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑初期モデルから得た張力(timestepを0.01にしている) 要素分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑99sと100sの間の様子 一瞬だけ張力が落ちているのが分かる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑得た張力をケーブル要素なしモデルにあてはめた変位
ちなみに、0.005sの場合~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
50sで下がって、100sで上がっている。これ...何が起こってい...
また、0.05s以降から衝撃力が出過ぎな気がするけど...これは...
上フランジ応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
<模索中に得られた結果>~
12/13~
PRE_EPSI、もしくは温度を用いた場合でもケーブル内の応力は...
試しに、ケーブル破断モデルに入れてみたが、桁が下にたわむ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑50sの様子 (scalefactorは100倍)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の拡大図
11/29~
Marcの設定を参考にしながら、Salome-Mecaでプログラムを組ん...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑張力を安定させた直後(50s)のモデルの様子 上がりすぎだろ.....
力を安定させた直後(50s)の桁の先端変位の比較~
,,Marc,Salome
,先端変位(mm),2.508,4321
あまりにもMarcとの差がありすぎるので、試しにMarcで得たビ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Salomeで得た張力データを入れた際の先端変位結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcで得た張力データを入れた際の先端変位結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcの結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の振動の拡大図 Salomeに入力した張力データの結果では...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の振動の拡大図 Marcで作成したグラフの70s後の振動は...
<ケーブル破断後モデルについて、現時点で考えること>~
・当初、Marcの出力結果がkNだと思っていたが、得られた結果...
また、SalomeもMarcも張力が[N]で出ていることは分かったので...
FORC_NODAでビーム軸力を得られるなら、何故値がずれるのか考...
仮にSalomeで正しいビーム軸力を得られなかった場合、Marcの...
⇒張力の見方を改めて考える必要がある。案としては、ケーブル...
それか、MULTIFIBERを用いてケーブルモデルの応力を求め、そ...
・力と時間のグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:FORCの合力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:応力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:応力から力に変換
12/6追記~
色々探ってみたが、バネと桁の影響を無視してケーブルのみのP...
一つ思いついた策として、ケーブルをcable要素に変えてやって...
ただ、問題点として、salomeのケーブル要素はコンクリートに...
また、cable要素は温度を考慮した解析をしてくれることは分か...
・Marcの軸力データをsalomeに入力したとしても、50sまでに桁...
***モデルBについて [#k829fe9a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ケーブルを2本にして桁を支えたモデルで外側のケーブルを破断...
手順としては初期モデルで外側のケーブル(C1)の張力を取って...
解析時間は150sで100sまでプレストレスと死荷重を掛け、その...
まず、桁の先端変位とC1の張力を取ってみたが、張力はマイナ...
先生に意見を頂いたところ、2本のケーブルのプレストレスが同...
,,変位,張力
,Salome,2.50mm,-32.28N
,Marc,2.51mm,-32.08N
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑先端変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑張力
***解析結果 卒論用 [#z22a073e]
グラフの縦軸は桁の先端変位(mm),横軸は時間(s) 119s〜126s間...
また、死荷重はFy=-0.017248N/mm,集中荷重は桁の先端にFy=-10...
・桁1mモデルA,Bの先端変位の動的応答値とDAF~
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,1.95,1.47
,0.01s,1.93,1.27
,1s,1.02,1.00
,5s,1.01,1.00
・桁10mモデルAの先端変位の動的応答値とDAF（モデルBは健全...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA
,0.005s,2.00
,0.01s,2.00
,1s,1.96
,5s,1.57
ただし、120sまでの健全モデルとケーブル破断モデルで桁の先...
・桁の先端に集中荷重を載荷した桁1mモデルA,Bの先端変位の動...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,1.97,1.47
,0.01s,1.95,1.27
,1s,1.03,1.00
,5s,1.00,1.00
・桁の先端に集中荷重を載荷した桁10mモデルA,Bの先端変位の...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,2.00,1.95
,0.01s,2.00,1.95
,1s,1.96,1.26
,5s,1.57,1.04
ただし、モデルAについては120sまでの健全モデルとケーブル破...
○考察~
DAFの値としては,1を超えているため大きな衝撃力が発生してる...
しかし、桁長10mの場合、健全モデルで得た桁の先端変位と120s...
そもそも、健全モデルで得た張力を元々ケーブルがあった箇所...
値が違うということは,何かしら原因があると思うが...。
**Marc Mentatで行う弾塑性解析とSalome_Mecaとの結果比較 [#...
***Salome_Mecaのビーム要素での弾塑性解析の手法 [#afd3f525]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ビーム要素のみのモデル（片持ち梁のモデルは10✕10✕100mmで...
分かったこと~
・1D要素(ビーム要素)で非線形解析を行う際は、Mulifiber(POU...
[[このページ:https://code-aster.org/forum2/viewtopic.php?...
・断面モデルもGeometryで作成してMeshで切る必要がある。(Mu...
・結果は_ELGAで見ること。NOEU,ELNOでは結果を見ることがで...
・反復回数や分割数(step数)は細かいほうがよい。(載荷値を増...
7/27~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑グラフが折れ曲がる付近の拡大図
[[応力とひずみのデータはこちら:http://www.str.ce.akita-u....
載荷側の断面上の真ん中の球体の応力ひずみデータを取った。
[[ビーム要素の弾塑性解析の設定方法:http://www.str.ce.akit...
DEFI_MATERIAU,DEFI_MATERIAU_GC、どちらでも可。
step数を細かくすると、235MPaで降伏していることが確認でき...
235MPaを超えると、ひずみが大きく増加するので、降伏してい...
応力がずっと同じ値で増える理由だが、恐らく周辺の要素(断面...
ちなみに、降伏後のひずみだが、降伏後の応力ひずみ曲線の傾...
***解析結果 [#bf4e606c]
10✕10✕1000mmの片持ち梁で曲げによる弾塑性解析を行い、235MP...
ビーム要素~
完全固定で設定した節点から応力ひずみ曲線を求めた。
salome-meca~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Marc~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
ソリッド要素~
Salome-meca~
7/24~
time0.5sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
time0.7sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
time1.0sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
応力(ミーゼス応力)-ひずみグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑SIEQ_NOEUとEPSI_NOEUを用いた解析結果~
[[グラフの数値はこちら:http://www.str.ce.akita-u.ac.jp/~g...
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
応力(ミーゼス応力)-ひずみ(相当ひずみ)グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑SIEQ_NOEUとEPEQ_NOEUを用いた解析結果~
[[グラフの数値はこちら:http://www.str.ce.akita-u.ac.jp/~g...
見たところは断面真ん中上で(②)、step数は50分割しています。~
SIEQ_NOEUのVMISでミーゼス応力を見れることを知ったので、試...
ミーゼス応力は相当応力とも呼ぶらしく、ひずみも相当ひずみ...
これらが正しいグラフなのか分からないが(そもそも、何故これ...
ちなみに、SIEQ_NOEUではミーゼス応力の他にもトレスカ応力や...
***ミーゼス応力 [#r131c6c9]
後藤メモ(23/9/8)：Salome-Mecaのミーゼス応力は、相当応力$\...
$\overline{\sigma}$は偏差応力のノルムみたいなもので、1軸...
せん断降伏応力を$\tau_{Y}$と書くと、1軸引張の引張降伏応力...
$\overline{\sigma}$の方を相当応力という場合もあり（更には...
この辺は注意が必要である。
-[[構造力学特論テキストにこの辺の話を書きました:https://w...
-[[岩熊先生・小山先生のテキスト参照:http://mechanics.civi...
***相当塑性ひずみ(EPEQ_NOEU) [#d51a592f]
このコマンドはCALC_CHAMP→CRETERIAから設定できる.~
意味は,相当塑性ひずみ（累積塑性ひずみ）で,EPSIから計算さ...
ParavisではINVA_2（ミーゼス相当塑性ひずみ）を使えば,相当...
ちなみに...相当塑性ひずみとは,各成分の塑性ひずみを用いて...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑片持ち梁の断面図
メッシュ長さ0.5 要素数394121でメッシュを切った。下の図は...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
timestepは50分割しており、引張側の真ん中の点が降伏したの...
,,降伏応力(MPa),ひずみ
,①,264.193,0.001099
,②,238.937,0.000947
,③,261.468,0.001107
,⑤,-259.924,-0.001051
,⑥,-248.081,-0.000970
,⑦,-267.69,-0.001079
timestep100の場合(春課題と同じ)~
引張側の真ん中の点が降伏したのがstep40
,,降伏応力(MPa),ひずみ
,①,264.201,0.001099
,②,238.937,0.000947
,③,261.412,0.001107
,⑤,-259.925,-0.001051
,⑥,-248.081,-0.000970
,⑦,-267.692,-0.001051
春課題のデータ
,,降伏応力(MPa),ひずみ
,①,249.267,0.001042
,②,234.025,0.000936
,③,236.751,0.000984
,⑤,-252.657,-0.001037
,⑥,-233.743,-0.000935
,⑦,-252.314,-0.001067
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
Marc
5/22(月)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑断面4隅の点を取った場合の応力-ひずみ曲線。4隅どの1点を取...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑②（引張側）の点を取って、応力-ひずみ曲線を求めた。
引張側の真ん中の点が降伏したincrement14で断面の4隅の状況...
,,降伏応力(MPa),ひずみ
,①③,240.8,0.001451
,②,236.8,0.001186
,⑤⑦,-240.8,-0.001451
,⑥,-236.8,-0.001186
グラフに関してはどうしても上降伏点と下降伏点が現れるグラ...
一応，設定している降伏値に近い値であればよいらしいので，...
しかし，1つ問題点?があり，固定面からずれた位置のほうが大...
②(引張側の真ん中の点)
,位置,降伏応力(MPa),ひずみ
,0mm,236.8,0.001186
,1mm,248.8,0.001116
,2mm,239,0.001091
,3mm,236.7,0.00114
4隅（数値は絶対値で表記）
,位置,降伏応力(MPa),ひずみ
,0mm,240.8,0.001451
,1mm,237.9,0.001304
,2mm,236.6,0.001117
,3mm,238.3,0.001098
5/17(水)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑②（引張側）の点を取って、応力-ひずみ曲線を求めた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑①,②,③の応力-ひずみ曲線をまとめたグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑④,⑧,⑨の応力-ひずみ曲線をまとめたグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑⑤,⑥,⑦の応力-ひずみ曲線をまとめたグラフ
・上降伏点と下降伏点が現れるグラフができた。~
グラフとしては間違っていないが、設定した降伏応力を超えな...
***結果比較 [#sa835e43]
ビーム要素~
ソリッド要素~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
,解析ツール,降伏応力(MPa),ひずみ
,Salome,238.937,0.000947
,Marc,236.8,0.001186
4隅のデータ~
Salome-meca
,,降伏応力(MPa),ひずみ
,①,264.193,0.001099
,③,261.468,0.001107
,⑤,-259.924,-0.001051
,⑦,-267.69,-0.001079
Marc(数値は絶対値)
,降伏応力(MPa),ひずみ
,240.8,0.001451
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
Marcのビーム要素でできるグラフ（指定した降伏応力で曲がる...
どちらも降伏応力はそれなりに近い値が出るが、グラフ形状が...
Salomeは指定した降伏応力を設定しても応力がずっと増えてい...
一方、Marcは上降伏点と下降伏点が現れるグラフができるが、...
が、固定面の応力より固定面から離れた位置のほうが少し大き...
また、4隅のデータについてだが、Salomeは4隅のデータそれぞ...
現段階で言えることは、ソリッド要素の弾塑性解析はSalome-me...
(追記：一旦、ビーム要素の弾塑性解析を重視で行っているので...
***ビーム要素での片持ち梁モデルの作り方 [#wedc87f7]
片持ち梁曲げの解析方法は以下の通りに行なった。（一部、[[M...
-長さ単位 (基本操作)~
長さを選ぶ(私はmm)~
-形状とメッシュ (基本操作)~
節点を2点追加 (0 0 0 と100 0 0)(x y z の順，スペース開け...
要素を線(2)に変えて追加 → 節点nodeを入力 → 右クリック確定...
-再分割 (オペレーション)~
分割 → 上からx,y,zの順に値を入力 → 要素 → 分割したいとこ...
-スイープ (オペレーション)~
節点 → 全て → OK (必要のない節点を消す？みたいな)~
-再番号付け (基本操作)~
全ての形状/メッシュ → OK (節点の番号をつけ直す？みたいな)~
形状設定
-新規(構造)▼ (形状特性)~
3ーD ▶ → 中実断面梁(中が詰まっている断面) → 形状選ぶ(長方...
X,Y,Zの値は変えたり変えなかったり(1が橋軸直角方向になるよ...
追加したいところ左クリック(全体選択でも可) → 右クリック確...
テーブルと座標系
-新規(1次独立)~
独立変数V1のタイプをtime(荷重を増分荷重にする)に → データ...
F(縦軸)の数値を入力してenter → データポイントの追加で先ほ...
自分は(0,0),(1,1)と設定した。
ちなみに、バイリニアだとここで応力ひずみの値を入れる必要...
なお、最初の点の塑性ひずみは0でその時の応力が降伏応力であ...
材料設定
-新規(構造)▼ (材料特性)~
有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/...
塑性にチェックを入れて降伏応力を入力(テーブルに応力-ひず...
追加したいところ左クリック(全体選択でも可) → 右クリック確...
境界条件設定
-(固定)~
新規(構造)▼ (境界条件)~
変位指定 → 変位X,Y,Zと回転X,Y,Zにチェックを入れて0にする(...
固定端にする方の節点を左クリック → 右クリック → 矢印出て...
-(載荷)~
新規(構造)▼ (境界条件)~
集中荷重 → 力Yに値を入力し、テーブルにtimeのtableを適用 →...
解析
-新規▼ (解析ジョブ)~
構造 → チェック → 初期荷重をクリックし、境界条件の載荷の...
解析結果をクリック → 要素テンソルのstress,Total Strain,El...
要素スカラーのEquivalent Von Mises Stressにチェック →~
OK → Job check complete になったらOK → 実行 → Marc実行(1)...
結果見る
-ポストファイルを開く(モデルプロット結果処理メニュー)~
→ 履歴プロット → 位置設定で見る点を囲んで右クリック → 全...
X軸データキャリアタイプを位置に → 位置を見たい節点 → 変数...
Y軸データキャリアタイプを位置に → 位置を見たい節点 → 変数...
カーブの追加 → フィット → OK
***ソリッド要素での片持梁モデルの作り方[#xc3ecfbb]
形状設定より先の設定は梁要素と同じなので、それより前の設...
-長さ単位(基本操作)~
基本はmmで行う。
-形状とメッシュ (基本操作)~
節点を2点追加 (0 0 0 と100 0 0)(x y z の順，スペース開け...
要素を線(2)に変えて追加 → 節点nodeを入力 → 右クリック確定...
-拡張 (オペレーション)~
線を拡張してソリッドをつくる。~
「拡張」→回転角（0 0 0）、並進（0 10 0）、繰り返し(1)と入...
回転角（0 0 0）、並進（0 0 10）、繰り返し(1)と入力→「要素...
これで100mm✕10mm✕10mmのソリッド要素が完成
-再分割 (オペレーション)~
ここでメッシュを作成する。~
分割 → 上からx,y,zの順に値を入力（好きな値でいいが、私は1...
→ 要素 → 分割したいところを左クリックで囲って右クリック →...
-スイープ (オペレーション)~
節点 → 全て → OK
形状設定
-新規(構造)▼ (形状特性)~
3ーD ▶ → ソリッド→ 要素を選択 → 追加したいところを左クリ...
ここからの設定は梁要素と同じです。注意点として、境界条件...
（例:点が100個あって、荷重を800N掛けるなら、入力値は8）
**Salome-Mecaの振動解析動画 [#m7a07b3d]
***記録用メモ [#sd5b7728]
5/19(金)遅くなりましたが、Tutorial9,10視聴完了。Tutorial9...
5/4(木)Tutorial8視聴完了。内容は高調波ベース加振を受けるD...
4/30(日)Tutorial7視聴完了。内容は非減衰MDOFシステムの周波...
4/27(木)Tutorial4で使用したモデルを作成し、動画と同じ挙動...
4/26(水)Tutorial6視聴完了。内容はMDOF法（多自由度振動系）...
4/24(月)Tutorial5,5A視聴完了。内容はSDOF法（1自由度系の振...
4/21(金)Tutorial4視聴完了。片持ち梁のモード解析の内容だっ...
4/13(木)Tutorial2B,3視聴完了。データの保存方法(.resuファ...
4/11(火)Tutorial2で使用したモデルを作成し、動画と同じ挙動...
4/8(土)Tutorial2 一通り視聴完了。点と線で片持ち梁モデルを...
*Salome-meca 知識 [#m3ab994c]
ここでは自分が学んだSalome-mecaの知識についてまとめる。ビ...
code-aster名称[[http://opencae.gifu-nct.ac.jp/pukiwiki/in...
**解析手法等の設定まとめ [#bb74bb5a]
***ビーム断面の適用方法 [#g75f9b28]
1D要素で弾塑性解析を始めとした材料非線形の解析を行う場合...
(1)DEFI_GEOM_FIBER~
まず、このコマンドを使う前にビーム要素の断面として使う断...
SECTIONは断面定義、FIBREは断面の穴あき部分を定義する?(例...
ビーム断面が厚肉断面もしくは空白部分に何も入れる必要がな...
TOUT_SECTはビーム要素に作成した平面モデルを全て適用したい...
GROUP_FIBERは好きな名前をつけてよい。~
MAILAGE＿SECTはLIRE_MAILAGEで適用した平面モデルを選択。~
COOR_AXE_POUTREとANGLEは0でよい。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
(2)AFFE_CARA_ELEM~
GEOM_FIBREが適用されていることを確認して、MULTIFIBREの選...
(ただし、AFFE_MODELEで_*Mと書いてあるキーワードを設定する...
(3)DEFI_COMPOR~
ここでは断面の材料特性と状態(弾塑性を入れるとか)を定義す...
このコマンドは画面上の{a}+の横にある{}+(Show All Commands...
スクロールして探すなら、Otherの欄にあるので、それを選択し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
適用すると、下図のような画面が現れるが、GEOM_FIBREはDEFI_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
MULTIFIBREを選択したら、GROUP_FIBREは(1)と同じ名前,mater...
最後に、DEFI_COMPORをAFFE_MATERIAUのAFFE_COMPORで適用する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
***動的応答解析のやり方 [#t4e8361f]
・注意点~
(1)DEFI_FONCTIONとAFFE_CHAR_MECAは、Pre Analysis・Analysi...
(2)DISCRETに質量マトリクスを用意しておくこと。DISCRETには...
この設定をしないと、ASSEMBLAGEでバネ要素に質量マトリクス...
[[詳しくはこちらを参照:https://code-aster.org/forum2/view...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)DYNA_VIBRAとCALC_CHAMPの名前は同じにすること。(同じ名...
(4)CALC_CHAMPのEXCITにDYNA_LINEで設定したCHARGEの分だけ入...
設定方法~
[1]Pre analysisのASSEMBLAGEを選択~
・ここでやること~
①NUME_DDLは{}押し、nddlと入力。~
②MATR_ASSEで質量マトリクスと剛性マトリクスを定義する。(質...
③CHARGEはAFFE_CHAR_MECAで適用した境界条件分だけ追加する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2-1]AnalysisのDYNA_LINEを選択~
・ここでやること~
①BASE_CALCULをPHYS,TYPE_CALCULをTRANに設定する。~
②EXCITに載荷条件(PRE_EPSIも同様に)の設定を適用。DEFI_FONC...
③SCHEMA TEMPSでNEWMARK,BETAを0.25と設定する。~
④INCREMENTはDEFI_FONCTIONで作ったテーブルの時間に合わせる...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2-2]DYNA_NON_LINEの場合~
ここで弾塑性を入れた動的解析ができると思われるが、そのあ...
・最低限やること
①SCHEMA TEMPSでNEWMARK,FORMULATIONをDEPLACEMENT(変位),BET...
ちなみに,[[説明書:https://code-aster.org/V2/doc/v13/en/ma...
更に詳しく調べると,これは直接積分法（正確には直接時間積分...
（Marcもニューマークβ法を用いて連鎖崩壊に関する解析をして...
[[参考サイト:http://jikosoftcom.blog25.fc2.com/blog-entry...
②INCREMENTはDEFI_LIST_INSTで設定したものを選ぶ。~
③EXCITにAFFE_CHAR_MECAで作成した数だけ設定する。載荷条件...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[3]Post ProcessingのCALC_CHAMPを選択
・ここでやること~
①CALC_CHAMPの名前をDYNA_LINEと同じ名前にすること!これをし...
②EXCITにAFFE_CHAR_MECAで適用した境界条件分だけ追加する。...
***PRE_EPSIについて [#c71fb860]
[[PRE_EPSIの説明書:https://code-aster.org/V2/doc/default/...
モデルをひずみを用いて変形させることができるコマンドであ...
1Dモデル(Multifiber要素も含む)ではEPX（梁のローカル軸に応...
(ky,kzはたわみの微分方程式から得られるdθ/dx = -M/EI の値...
X,Y,Zの軸の向きはグローバル座標(salome画面の色がついてる...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
単純梁モデル(10✕10✕120mm)にEPXX：0.01を与えると、以下の挙...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
線に色がついてるところが変形前のモデルの端であり、写真を...
・簡易モデルに用いた際に分かったこと~
現在、ケーブルにPRE_EPSIを入れており、縮む挙動(ケーブルが...
その理由について考えてみたが、要素を構成している節点にモ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑PRE_EPSIの意味 日本語に訳すと要素に変形の負荷を与えるみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑簡単な例 PRE_EPSIで与えたひずみになるような外力をかけて...
**Geometry [#i477a78c]
(1)マウス操作について~
マウスホイールを動かす or ctrlキー+マウス左クリック押し続...
ctrlキー+マウスホイールを押し続ける：移動~
ctrlキー+マウス右クリック押し続ける：回転~
(2）つくったモデルを消さずに長さ等を調整する方法~
作成したモデルの長さなどを変えたいとき、わざわざ作り直さ...
Geometryで左上のファイル→スクリプトの保存→ファイルの種類...
すると、下図のような画面が現れるので、点の位置や線の長さ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
**Mesh [#j28db1c4]
(1)メッシュのcompoundについて~
Meshにもcompound機能があり、別々に切ったメッシュを1つにく...
これを使えば、Marcの簡易モデルと同じメッシュの切り方がで...
**Aster_study [#qf331ace]
***Model Definition [#i67b38f6]
・MACR_CARA_POUTREとDEFI_GEOM_FIBERの違い~
どちらのコマンドもビーム要素の断面特性を定義できるが、MAC...
・AFFE_MODELEとAFFE_CARA_ELEMの対応コマンド~
AFFE_MODELEで設定したmodelisationはAFFE_CARA_ELEMのコマン...
対応しているものを以下にまとめたので、参考にしてください。~
POUTRE：POU_D_E, POU_D_T, POU_D_TG, POU_D_T_GD, FLUI_STRU...
BARRE：BAR, CABLE_GAINE~
CABLE：CABLE, CABLE_POULIE~
COQUE：COQUE_AXIS, COQUE_C_PLAN, COQUE_D_PLAN, DKT, DST, ...
DISCRETE：DIS_T, DIS_TR~
DISCRETE_2D：2D_DIS_T, 2D_DIS_TR~
MULTI_FIBRE：POU_D_EM, POU_D_TGM~
・AFFE_MODELE ビーム要素でのMODELISATIONの選び方~
梁の種類によって、選択するものが異なる。~
・真っ直ぐな梁：POU_D_E(オイラー梁), POU_D_T(ティモシェン...
・マルチファイバーの梁(鉄筋コンクリートや非線形材料の解析...
・棒(鉄筋とか断面が円のもの？) 2D_BARRE, BARRE 引張と圧...
・ケーブル CABLE 引張荷重 （引張だけ出せるが、斜張橋など...
・AFFE_CARA_ELEM DISCRETの使い方
DISCRETは離散という意味でバネ要素を定義する場所である。こ...
アルファベットの意味~
A:減衰~
K:剛性~
M:質量~
T：並進運動~
TR:並進運動&回転~
D:対角行列~
L:線~
N:節点~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑K_TR_D_Lにした場合の設定画面~
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
GLOBALはグローバル座標,LOCALはローカル座標を示す。
***BC & Load [#pcf9b829]
(1)AFFE_CHAR_MECA DDL_IMPO 完全固定の別のやり方
これまで完全固定をする場合、DX=DY=DZ=0と設定していたが、L...
(2)AFFE_CHAR_MECA FORCE_の意味
FORCE_NODAL ノードあたりの力
FORCE_FACE 1面あたりの力
FORCE_ARETE エッジ1本あたりの力(1D要素では使えなかった。)
FORCE_CONTOUR 輪郭の長さあたりの力
FORCE_POUTRE 梁の長さあたりの力
***Analysis [#eb94d3a9]
***Post Processing [#v995178b]
ELGA：ガウス点のelementごとに応力または力またはモーメント
ELNO：node点におけるelementごとの応力または力またはモーメ...
EFGE_ELNO:要素ローカル軸の節点における要素ごとの力または...
・CALC_CHAMPとCREA_CHAMPの違い
CALC_CHAMP を使用して、既存のデータからデータ構造を計算し...
CREA_CHAMP を使用して新しいフィールドを初期化し、後で分析...
CREA_CHAMPでは、構造物にどれだけ温度を与えるか設定できる...
**エラーメッセージ一覧 [#eb4955ab]
・このエラーメッセージは構造的に不安定なときに出る。~
これが出た場合は基本的にAsterstudyの設定を間違えているこ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_CHAR_MECAの設定をモデルに適用できない場合に出る。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_MODELEとAFFE_CARA_ELEMの条件が一致しない場合に出る...
AFFE_MODELEで設定したMODELISATIONは、それぞれAFFE_CARA_EL...
(例) AFFE_MODELEでBARREを選択した場合→AFFE_CARA_ELEMではB...
AFFE_MODELEでPOU_D_Eを選択した場合→AFFE_CARA_ELEMで...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Paravis [#zdf3c7e9]
(1)評価項目(magnitudeとか選択するところ)について~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
magnitude：全項目のΣ2乗の平方根(ただし、DEPLだけは変位量...
DEPLのDX：X方向変位量~
SIGM_NOEUのSIXX：X方向垂直応力~
SIGM_NOEUのSIXY：XYせん断応力~
SIEQ_NOEUのVMIS：ミーゼス応力(ちなみに、ミーゼス応力とは...
今載せたのは一例だが、他にも[[このページ:https://qiita.co...
(2).csvファイルの数字とParavisの評価項目~
.csvファイルの1行目の()の数字(例えば,DEPL(0)とか)はPavari...
例を挙げると、SIGM_NOEUの場合、SIXXはSIGM_NOEU(0),SIYYはS...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*TeX [#wcbe1bb5]
久し振りに使うとコマンド等を色々忘れているので、ここにメ...
(1)余白の埋め方
入れる場所はタイトルや名前を書くところより上に(こんな感じ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
{0pt}にすると余白がなくなる。図などを貼り付けた際にできる...
\setlength\textfloatsep{0pt} %本文と図の間の余白
\setlength\floatsep{0pt} %図と図の間の余白
\setlength\intextsep{0pt} %本文中の図の余白
\setlength\abovecaptionskip{0pt} %図とキャプション（図や...
(2)画像の貼り方
ファイルの語尾が.pdf_texの場合(inkscapeで編集したグラフな...
.png(スクリーンショットした画像など)の場合、\includegraph...
(3)数式について
数式で表すには$で挟む必要がある。
(例)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
*参考資料等 [#a2f72b97]
・研究活動で使う参考資料
[[code_asterの動画:https://www.youtube.com/@DigvijayPatan...
[[code-asterの説明:http://dugi6514.odns.fr/doc/v15/en/ind...
[[beam要素の弾塑性解析マニュアル:https://www.str.ce.akita...
[[CAE用語辞典:https://www.cybernet.co.jp/ansys/learning/g...
[[Marc Mentat2018メモ:https://www.str.ce.akita-u.ac.jp/cg...
[[例題で学ぶMarc有限要素法解析入門:https://www.lib.akita-...
[[Salome-Mecaのチュートリアル:https://code-aster.org/V2/s...
[[Pavaris wiki:https://www-paraview-org.translate.goog/Wi...
[[Salome-meca beginning with code_aster:https://archives....
・UNIXコマンド一覧 リンク集
[[https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html...
[[https://proengineer.internous.co.jp/content/columnfeatu...
[[https://amorphous.tf.chiba-u.jp/lecture.files/seminar-I...
・弾塑性解析
春課題で行なった弾塑性解析については以下のURLにまとめてあ...
[[https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E6%...
*日記 [#e5e63424]
3/31(日)~
3年生から更新してきた卒論日誌も本日で最終更新です。明日か...
これからの2年間、色々なことが待ち受けていると思いますが、...
そして、これまで鋼橋の研究をされていた先輩方という自分に...
3/22(金)~
卒業式が無事に終了しました。今日で3年生の後期からずっと一...
けど、またいつかどこかで会う機会はあるはずです。そのとき...
3/8(金)~
無事に謝恩式が終了しました。皆さんが協力していただいたお...
挨拶のとき、皆さんの前に立ったとき、緊張もありましたが、...
秋田での生活はまだまだこれからも続きますが、まずは卒論本...
来年度から鋼橋の研究する先輩が誰もいなくなるから、自分が...
2/21（水）~
久々の更新です。卒論発表が無事に終わり、先生方からの評価...
そしたらなんと、学部生の中では1番良い評価をいただいたので...
ただ、自分の研究成果はまだまだ途中段階ではあるし、まだこ...
12/31（日）~
2023年も残りあとわずかです。~
今年も色々ありましたが、楽しく過ごすことができました。~
特に今年は自分のことを色々見直す機会が多く、もっと人とし...
来年は1月から3月にかけて、卒論に東北支部、それに謝恩会の...
2024年はどんなことが待っているのだろうか…~
来年もよろしくお願いしますm(_ _)m
11/28(火)~
11月ももうすぐ終わり...早いな(・・;)~
来月でいよいよ2023年も終わり、中間発表に卒論が近づいてき...
卒業できないということはないと思うけど、このままだと面白...
残りの時間でどこまでやれるか分かりませんが、先生方や先輩...
11/17(金)~
今日、創造工房実習でSalome-Mecaを教えました。~
皆さんの前で立って教えましたが、人にわかりやすく教えるの...
あの説明で、ある程度理解することはできたのだろうか...。~
次は前で教えることはありませんが、反省点もいくつか見つか...
11/5(日)~
10月下旬あたりからずっと考えていたプレストレス問題だが、...
何で線形静的解析でしかできない(対応していない)のか不思議...
...卒論まで残り3ヶ月だし、いい研究成果が発表できるように...
10/24(火)~
家の近くに熊が出没するという中々恐ろしいことがありました...
思えば10月も来週で終わり、時が経つのが早いなと感じます。~
来年の4月から就職する同級生や先輩と過ごすのもあと半年もあ...
10/20(金)~
ようやく家のネット環境が復帰して、早速自宅で日記を書いて...
…さて、今日は構造研に配属された3年生と顔合わせをしました。~
3年生達の姿を見て、去年の自分達もこんな感じだったのかなと...
自分は大学院に進むので、少なくとも1年半は一緒に過ごすこと...
よろしくお願いしますm(_ _)m
10/18(水)~
本日、3年生の研究室配属先が決まったそうです。~
すでに何人かは顔合わせをしましたが、全員での顔合わせは今...
同時に、自分もしっかりしなきゃいけないとも思います。~
卒業された先輩方や今の院生のようにちゃんとサポートできる...
10/11(水)~
学部生での研究活動も後半戦、他の研究室でもそろそろ卒論に...
いよいよ今週末に3年生も研究室に入ってくるので、悪い手本に...
ただ...最近、細かいミスをしたり、人の話を1回で正しく理解...
あまりこういう状況が続くのはさすがにまずいと思う(場合によ...
9/29(金)~
中間発表が無事に終わり、そして夏休みも終わりました。~
発表に対しても初めて発表をした2月の頃と比べると成長したな...
...少しは変われたかな？~
来週から後期が始まりますので、引き続き研究に取り組みたい...
9/12(火)~
今週は昼過ぎまで集中講義を受講し、夕方から研究をするとい...
calculixでプログラムを作っている最中で、どうしても分から...
それでも簡易モデルの再現ができなければ、どうするべきか相...
9/8(金)~
昨日は帰ってきて何もせずに寝てしまったので、部屋の掃除や...
ひとまず、実家でやっていたデータを研究室のパソコンに移し...
中間発表までは、平日は解析をメインで、休日は概要とスライ...
9/7(木)~
今日、秋田に戻ってきます。~
実家での夏休みは父の実家に行ったり、インターンシップに参...
秋田に戻ったら、再び研究一筋の日々になりますが、卒論に向...
8/6(日)~
秋田は竿燈祭りで盛り上がっています(^^)~
自分は平日に見に行ったので、昨日今日はのんびり過ごしつつ...
後は...台風がどうなるか心配です。無事に帰省できるといいけ...
8/2(水)~
前期のゼミは本日で終わりました。...早いな(^_^;)~
Salomeを使えるノートパソコンを借りることができ、これで実...
ちなみに、明日から竿燈祭りがあるらしいです。人はたくさん...
7/26(水)~
久々の更新です。（4月〜7月の日誌はwikiの整理の際に消して...
今日、外国語文献の講義での発表が終了しました。質問の際、...
けど、4月のときより聞き取る力は上がっていると思うので、英...
気が向いたら、こんな風に更新していこうと思います。
*創造工房実習 課題まとめ [#r8b139fc]
**創造工房実習 発表準備進捗状況記録用メモ（提出済み） [#i...
2/3（金）
課題は完成しました。後はスライド作りをするのみです。（現...
1/31(火)
一通り確認し、後はサンドウィッチ梁の考察について思うとこ...
1/27(金)
概要は提出できる状況までできましたが、提出期限まであと1週...
今日からスライド作成に入りたいと思います。（2/10までには...
1/26(木)
まとめまで一通り書き終えました。後は図の配置や適切な大き...
1/24(火)
異方性2次単純梁の解析完了。サンドウィッチ梁について書く材...
1/23(月)
異方性2次単純梁の解析を行っている途中である。（サンドウィ...
後、LaTeXに貼り付けた図の位置がいまいちなので、調整するこ...
1/20(金)
解析結果と考察を1つにまとめることにした。
残りの作業は解析結果と考察に関してはサンドウィッチ梁を書...
1/17(火)
論文の構成を変更（はじめに→材料性質と解析方法→解析結果→考...
はじめに、材料性質と解析方法を一通り書いた。（図の作成はi...
1/16(月)
創造工房実習の課題の解析結果を一部編集しました。（単純梁...
解析結果の編集に入りました。
1/13(金)
課題で使う全てのグラフをinkscapeで修正しました。
とりあえず、課題の論文の構成は（はじめに→解析結果→考察→ま...
また、論文のタイトルと名前の修正を行い、「はじめに」の部...
1/11(水)
wikiに貼り付けたグラフをsvgに変換させました。
メモ：inkscapeで凡例のみを消す方法はグラフをノードツール...
**12/2(金)課題 サンドウィッチ梁（木材を鋼材で挟む）[#wb04...
異方性の単純梁(1次)と2次要素のサンドウィッチ梁（鋼材：等...
2次要素のサンドウィッチ梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,583289,-,-,千代岡
,0.6,214850,-,-,高井
,0.7,155266,0.0861,-13.0,関合
,0.8,138453,0.083487,-15.7,岡田
,0.9,82766,0.083312,-15.8,松田
,1.2,32279,0.083574,-15.6,青野
,1.3,28343,0.083668,-15.49,山口
,1.4,23667,0.083680,-15.48,山本
,1.5,19958,0.083516,-15.6,進藤
,1.6,19451,0.086037,-13.1,河合
,1.8,10933,0.084022,-15.13,山口
,2,10764,0.083324,-15.8,進藤
,3,3618,0.083497,-15.66,山本
,4,1623,0.0852,-13.9,関合
,5,1007,0.083104,-16.1,千代岡
,6,842,0.0821,-17.1,高井
,7,554,0.080750,-18.4,青野
,8,289,0.079715,-19.5,岡田
,9,261,0.078427,-20.78,松田
,10,232,0.082495,-16.67,河合
異方性1次単純梁と2次要素サンドウィッチ梁の比較グラフ （ih...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：鋼材を挟むのと挟まない場合では理論値が約5倍異なり、...
***おまけ1 異方性単純梁を2次要素にして比較してみる [#qc21...
2次要素の異方性単純梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%)
,0.5,604167,-,-
,0.6,202805,-,-
,0.7,145019,0.516,5.1
,0.8,141896,0.515,4.8
,0.9,91974,0.514,4.6
,1.2,24520,0.512,4.1
,1.3,23810,0.512,4.1
,1.4,17725,0.511,4.0
,1.5,15433,0.511,3.9
,1.6,15746,0.510,3.7
,1.8,11647,0.509,3.5
,2,10460,0.509,3.5
,3,2298,0.504,2.5
,4,1482,0.500,1.6
,5,431,0.496,0.8
,6,356,0.497,1.1
,7,196,0.487,-0.8
,8,104,0.491,-0.14
,9,81,0.485,-1.3
,10,78,0.488,-0.75
（ihou2：11/25の異方性1次解析のデータ ihou2ji:2次要素の異...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
***おまけ2 1次要素のサンドウィッチ梁の場合は？[#gd69b506]
1次要素のサンドウィッチ梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%)
,0.5,583289,0.08164,-17.5
,0.6,214850,0.07724,-22
,0.7,155276,0.07681,-22.4
,0.8,139020,0.07642,-22.8
,0.9,83327,0.07678,-22.4
,1.2,32182,0.07602,-23.2
,1.3,28507,0.07652,-22.7
,1.4,24105,0.07563,-23.6
,1.5,19934,0.07602,-23.2
,1.6,19653,0.07603,-23.2
,1.8,12468,0.07620,-23
,2,10647,0.07373,-25.5
,3,3599,0.06616,-33.2
,4,1665,0.05568,-43.8
,5,1067,0.05320,-46.2
,6,834,0.03034,-69.4
,7,554,0.01788,-81.9
,8,289,0.01806,-81.8
,9,261,0.01413,-85.7
,10,232,0.03341,-66.3
異方性1次単純梁と1次要素&2次要素サンドウィッチ梁の比較グ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：サンドウィッチ梁を1次要素に変えたらどうなるか気にな...
2次要素のサンドウィッチ梁の解析結果を比較すると、理論値よ...
11/25(金)の課題でも述べたが、2次要素の解析が1次要素より正...
**11/25(金)課題 木材の単純梁の2次要素と異方性1次[#xdf73a23]
＜等方性1次のデータは11/18（金）課題を参照＞
***一次と二次要素の等方性比較 [#j95e0189]
二次要素の等方性のデータ(メッシュ長さ0.5は要素数が多すぎ...
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,-,-,千代岡
,0.6,203209,0.42383,0.98,高井
,0.7,145234,0.43301,3.22,関合
,0.8,140987,0.43006,3.2,岡田
,0.9,91974,0.429913,3.18,松田
,1.2,24800,0.429777,3.14,青野
,1.3,23132,0.42989,3.16,山口
,1.4,17617,0.429745,3.13,山本
,1.5,15433,0.429844,3.2,進藤
,1.6,15900,0.429754,3.13,河合
,1.8,11677,0.42962,3.1,山口
,2,10460,0.429605,3.1,進藤
,3,2486,0.429217,3.0,山本
,4,1453,0.4293,3.02,関合
,5,431,0.427885,2.69,千代岡
,6,360,0.4282,2.78,高井
,7,196,0.42606,2.25,青野
,8,104,0.42631,2.3,岡田
,9,81,0.42513,2.03,松田
,10,78,0.424466,1.8,河合
一次と二次要素の等方性比較のグラフ（niji：二次要素のデー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：[1次要素は頂点の節点だけで要素の変形を表現する、2次...
（引用元：日経XTECH 第2回精度はメッシュで決まる https://x...
↑2次要素の解析での相対誤差が1次要素より小さいのは恐らくこ...
グラフを見ると、1次要素よりも理論値に近い解析ができており...
実際、メッシュ長さ0.5の解析はできなかった。（他のパソコン...
***等方性1次と異方性1次の比較 [#l937a0df]
異方性1次のデータ(メッシュ長さ1.0と1.1追加)
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.50919,3.56,千代岡
,0.6,203209,0.50472,2.6,高井
,0.7,145234,0.5036,2.42,関合
,0.8,140987,0.50283,2.3,岡田
,0.9,91974,0.50053,1.8,松田
,1.0,77582,0.4999,1.6,千代岡
,1.1,27533,0.4886,-0.63,千代岡
,1.2,24800,0.48739,-0.9,青野
,1.3,23132,0.48841,-0.67,山口
,1.4,17617,0.48403,-1.56,山本
,1.5,15433,0.48202,-2.0,進藤
,1.6,15900,0.48329,-1.7,河合
,1.8,11677,0.47855,-2.67,山口
,2,10460,0.47906,-2.6,進藤
,3,2436,0.42787,-12.98,山本
,4,1453,0.42772,-13.02,関合
,5,431,0.27364,-44.3,千代岡
,6,360,0.33927,-31.0,高井
,7,196,0.21363,-58.5,青野
,8,104,0.22574,-54.1,岡田
,9,81,0.22750,-53.7,松田
,10,78,0.20327,-58.7,河合
等方性1次と異方性1次の比較のグラフ（ihou：異方性1次のデー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：解析結果から異方性1次単純梁の変位は約0.50mmに落ち着...
後、等方性1次のグラフと比較すると、お互いグラフの形が似て...
**11/18(金)課題 木材の単純梁 (メッシュ長さ1.0と1.1追加)[#...
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.428982,2.94,千代岡
,0.6,361584,0.421233,1.09,高井
,0.7,145234,0.4225,1.4,関合
,0.8,140987,0.422627385,1.4,岡田
,0.9,91857,0.420351606,0.88,松田
,1.0,77582,0.4192,0.6,千代岡
,1.1,27533,0.4057,-2.6,千代岡
,1.2,24520,0.404744325,-2.87,青野
,1.3,23132,0.4045,-2.93,山口
,1.4,17530,0.3986,-4.34,山本
,1.5,15433,0.396317756757,-4.9,進藤
,1.6,15900,0.399049,-4.24,河合
,1.8,11677,0.404457,-2.94,山口
,2,10460,0.394818715517,-5.3,進藤
,3,2344,0.32447,-22.13,山本
,4,1453,0.3329,-20.1,関合
,5,431,0.136240,-67.3,千代岡
,6,360,0.2130486,-48.9,高井
,7,196,0.1019892,-75.5,青野
,8,104,0.1158624,-72.2,岡田
,9,81,0.1247076,-70.1,松田
,10,78,0.07733,-81.4,河合
単純梁のグラフ：縦軸 変位(mm) 横軸 要素数 （rironti1118...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：解析結果から、単純梁のグラフの傾きが0に近づくのは約...
**11/11(金) 課題 木材の片持梁[#g9ad393e]
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
,メッシュ長さ,要素数,先端変位(mm),相対誤差(%),計算者
,0.5,59504,6.56,-1.5,千代岡
,0.6,45512,6.48774,-2.69,高井
,0.7,39075,6.54133,-2.0,関合
,0.8,13397,6.43695,-3.5,岡田
,0.9,9903,6.36315,-4.6,松田
,1.2,6256,6.3043375,-5.4,青野
,1.3,5767,6.29784,-5.6,山口
,1.4,5146,6.286015,-5.76,山本
,1.5,3935,6.24807,-6.3,進藤
,1.6,3400,6.20446,-6.98,河合
,1.8,2952,6.17161,-7.5,山口
,2,1632,5.64585,-15.3,進藤
,3,667,5.4053975,-18.96,山本
,4,264,3.6161,-45.8,関合
,5,191,3.86,-42,千代岡
,6,190,2.5077325,-62.4,高井
,7,75,1.41225,-78.8,青野
,8,56,1.2887175,-80.7,岡田
,9,49,1.28799,-80.9,松田
,10,44,1.226075,-81.6,河合
片持ち梁のグラフ：縦軸 変位(mm) 横軸 要素数 （katamochi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ:解析結果から、片持ち梁のグラフの傾きが0に近づくのは...
**11/4（金）グラフ作成 [#v961224d]
・課題で作成したグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
・創造工房実習時に作成したグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
*3年次の日記(2022/10〜2023/2) [#g01e77dc]
2月10日(金)
今日で発表含め創造工房実習の課題が全て終わりました。ただ...
次は春休み課題（弾塑性解析）を行う必要があるので、それに...
1/31(火)
今日は中村先生の講義に参加しました。
特に印象に残ったのは、ケーブルに乾燥した空気を送りこんで...
ケーブルに関する知識は鋼橋の研究をする際に必要になる知識...
1/11(水)
今日の創造工房実習でLaTeXについて一通り学びました。
後は課題提出に向けて作業するのみなので、頑張っていきたい...
1/9(月)
遅くなりましたが、2023年が始まりました。今年もよろしくお...
明後日、創造工房実習があるので、LaTeXについて昨年学んだこ...
また、創造工房実習の課題の締切まであと約1ヶ月後なので、計...
12/16(金)
LaTeXについて学びました。色々覚えることはありそうですが、...
12/14(水)
これまで行なってきた解析データやグラフを少し整理しました...
12/9(金)
本日は中村先生の特別講義を聴講しました。東京湾アクアライ...
明石海峡大橋のお話は別の機会でしてくださるそうなので、次...
12/8(木)
12/2の課題を終わらせました。明日は特別講義ですが、どんな...
12/5(月)
salome-mecaの単純梁の復習とサンドウィッチ梁の1次要素解析...
2次要素の結果はまだ出ていませんが、結果がでた後に1次とど...
12/3(土）
salome-mecaで行う片持梁解析の復習を行いました。恐らく、来...
12/2(金)
本日で創造工房で習うsalome-mecaの演習が終わりました。今ま...
12/1(木)
11/25の課題を終わらせました。明日はサンドウィッチ梁の解析...
今年も残り1か月です。(月初めから大雪でしたが...(-_-;))
2022年の残りの期間も楽しく過ごしつつ、自分の知識を高めて...
2022年11月の記録 [#o8b95629]
11/28(月)
2次要素メッシュ長さ0.5の解析を違うパソコンを使って行いま...
11/25(金)
今日も先週に引き続きsalome解析を学びました。ただ、2次要素...
また、今週は研究室で色々な方からお話を聴きました。中には...
11/18(金)
salomeで単純梁の解析について学びました。やることは単純梁...
11/11(金)
salomeについて学びました。次回もsalomeを使った解析を行う...
11/9(水)
viiとgnuplotの使い方を軽く復習しました。
11/4（金）
初めて昼休憩での英語会話に参加しましたが、皆さんの会話を...
理解して自分の考えを英語で話すことができるように英語力を...
また、本日の創造工房でgnuplotとviについて学びました。
今後、salomeなど様々なことを学びますが、これらの機能を使...
11/2（水）
UNIXコマンドを実際に使ってみました。昨日調べたものを含め...
これから様々な場面でUNIXコマンドを使う機会があるはずなの...
（ちなみに...ctrl+zで実行しているコマンドを終了させること...
11/4追記：ctrl+zはどうやらコマンドの強制終了らしいので、...
11/1（火)
UNIXコマンドについて調べました。使えそうなコマンドをいく...
＜取り組むこと＞
・UNIXコマンドについて学ぶ。
・構造・材料系の本を読む
2022年10月の記録 [#z5e925c3]
10/28（金）
UNIXコマンドについて学びました。資料に載っているコマンド...
10/27(木)
タッチタイプの一通りの流れをストップウォッチで計測したら...
しかし、キーボードを見ながらの記録なので、できる限り、下...
（追記）木材についての勉強を始めました。（本：プロでも意...
10/24（月）
タッチタイプの練習を行いました。大文字と小文字を交互に打...
10/21(金）
タッチタイプの練習を行いました。まずは、3分以内に打つこと...
＜取り組むこと＞
・タッチタイプの練習
・構造・材料系の本を読む
終了行:
#contents
2024年4月からは[[千代岡の修論日誌:https://www.str.ce.akit...
3月までに書いたことはそのままこのページに残しておきますが...
研究内容：様々な外力が作用する斜張橋の動的応答解析
*ちょっとした春課題 [#u272b919]
model300のケーブル7本モデルのC1破断時の主桁中央の鉛直変位...
減衰率の設定法は[[秋山さんの修論日誌:https://www.str.ce.a...
・減衰0%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰0.5%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰1%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰2%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰5%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・減衰10%~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
基本的に減衰率の値以外の設定条件は先輩方が設定したものと...
減衰率5,10%だと主桁中央変位の最もたわむ値が変わっているが...
ステップ条件を複数条件→荷重時間に対する初期比率を0.01にや...
先生に結果をお見せしたところ,問題なさそうとおっしゃていた...
*データ類の置き場所 [#uc29fb98]
R3年度バックアップの根本さんのお名前が書いてある外付けHDD...
ファイル名のDは死荷重,Cは集中荷重を示しています.また,どの...
*課題 [#e9b3591e]
自分の研究成果の課題点（院でやること）をまとめておく。
・実験的検討を行う。解析結果とSalome_Mecaで出したデータを...
・橋モデルをSalome_Mecaで作成して,今回得られた結果とどう...
*卒論内容 [#a313a74b]
研究室で編み出した連鎖崩壊のメソッドの妥当性を確認する。...
しかし、その妥当性の評価を行なっていなかったので、その評...
Salome-mecaで行う理由としては...~
①ビーム要素を用いた動的応答解析の手法ができていないのでそ...
②2種類以上の解析ソフトで解析できるようにした上で、同じよ...
卒論としての目標は、簡易モデルでのDAFを用いた連鎖崩壊に関...
DAFとは？~
動的増幅率を表す指標。式等は[[こちら:https://www.str.ce.a...
*今後の予定 [#b61dcef6]
3/22 卒業式~
3/23 一旦実家に帰省~
4月以降については修論日誌に書く。
*やることリスト [#naf093be]
大学卒業までにやること
まずは東北支部の発表に備える！~
・本論作成（3/15まで）~
・引継ぎ関連（先輩方にお願いしないと...）~
・以前、Salomeで作成した橋モデル作成の続き~
・wiki整理~
・謝恩会準備（研究とは関係ないが...）~
卒論まで残り1週間 残りの期間でやること
・2回目練習の動画を見て,反省会（メモしきれなかったところ...
・スライドの修正~
・概要は卒論用に手直ししたので、一度指導教員にお見せする~
・ひたすら練習.個人練習を含め,15回くらいはやりたい~
2024年 1,2月にやること~
・0.01sから短めにtimestepを切ると途中で張力が落ちる問題を...
・東北支部、卒論に向けての概要とスライド作成~
・(研究とは関係ないけど...)謝恩会の準備
ラストスパート！！自分のできる限りのことは尽くす！
2023年 12月にやること~
・簡易モデルで得られる張力の値がMarcと大きく異なる謎を解...
・橋モデル:model300をsalomeで作成する。(橋はすでにgeometr...
・PRE_EPSIを使ったモーダル解析(山本君と共同研究)~
・12/22に向けての中間発表の準備~
今年もいよいよ今月で終わり、これまでの成果を自信を持って...
2023年 11月にやること~
・動的解析(線形&弾塑性を含んだver.)でプレストレスを考慮し...
・model300をsalomeで作成する。~
・PRE_EPSIの解析手法の検討(山本君と共同研究)~
・ビーム張力を見れるようにする。~
2023年 10月にやること~
・簡易モデルの弾塑性周辺の設定見直し→ケーブルには何も設定...
・ピン接合の設定確認(ばねが本当に一軸方向のみに作用してい...
・簡易モデルにプレストレス入れてみる~
・張力だけを落として集中荷重のみが載荷し続ける挙動再現(片...
・Salome-mecaの簡易モデルを用いた動的応答解析のやり方を学...
・Salome知識(ソリッド絡みの知識)復習 10月に来る3年にSalom...
2023年夏休みにやること
8月は実家に戻るまではピン接合問題(Salome,calculixの両方の...
Salomeを使えるノートパソコンを借りることができたので、実...
何か進捗があれば、先生に報告することを忘れないように。~
秋田に戻ったら、中間発表に向けて今まで以上に気合を入れて...
2023年7月にやること
(1)Salome-Mecaでビーム要素の弾塑性解析をできるようにし、...
⇒ケーブルには引張しか存在しないので、それを再現するため。
(2)Salome-Mecaでピン接合をできるように考え、Marcの簡易モ...
今は上記2つを重点的に行う。今のままだと、10月の中間発表で...
(3)英語ゼミの中間発表に向けてスライド作り~
⇒これは研究背景と今取り組んでいることとこれからどうしてい...
後は、自分にとって必要だと思う力(コミュニケーション能力と...
焦ってもしょうがないが、自分のやっていることなどを人に分...
2023年6月にやること
(1)Salome-Mecaでビーム要素の弾塑性解析をできるようにする。
(2)Salome-Mecaでソリッド要素の弾塑性解析のグラフ改善（Mar...
(3)Salome-Mecaでピン接合をできるように考える。
(4)振動解析の勉強（動画のモデルの実際に作る）
後は、来月に大学院入試があるので、それに向けて準備する。
2023年5月にやること
(1)先月同様、salomeでの振動解析について学ぶため、code_ast...
(2)Marcで弾塑性解析（モデルは10×10×100mm）を行い、235MPa...
Salome-MecaとMarcの結果比較をするために、Salomeで梁要素,M...
(3)Salome-Mecaでピン接合をできるように考える。
後は、時間があるときに院試に向けて勉強や本などを読んで知...
2023年4月にやること
(1)code_asterの動画[[https://www.youtube.com/@DigvijayPat...
（字幕は英語のほうがよい。...が、話している人の声が小さい...
目標はGW明けに終わらせる！（どんなに時間がかかっても5月中...
参考資料
-[[Salome-Meca説明書:https://sites.google.com/site/codeas...
-[[code-asterの説明:https://code-aster.org/V2/doc/default...
(2)時間があるときにMarcの使い方に慣れておく。
*進捗状況 [#be3fe568]
**ゼミ進捗報告 [#l52104c3]
12/12 温度考慮した動的解析ができるようになったこととケー...
12/6 現状の問題点と打開策を話す。~
11/29 軸力絡みの話~
11/22 ビーム軸力の話~
11/15 特になし~
11/8 PRE_EPSIを使ったプレストレス再現のお話~
11/1 線形での動的解析の設定の仕方は分かった。(FORC,REACの...
10/25 線形解析ではプレストレスが導入できたこと、動的解析...
10/18 プレストレス導入の話~
10/11 プレストレスとピン接合絡みの話~
10/4 後期1回目のゼミ 報告することは特になし~
**ビーム要素の簡易モデルを用いた連鎖崩壊のメソッドの妥当...
現時点で分かっている問題点と分かったこと
-曲げモーメントの出力できない．応力もできない?軸力は可能...
→SIPO_NOEU,ELNOで曲げモーメントを見ることはできるみたい。...
ただ、REAC_NODAのDRX,DRY,DRZでモーメントを見ることはでき...
[[こちらの資料:https://code-aster.org/V2/doc/default/en/m...
-現時点でケーブルはトラス要素を用いている．おそらく圧縮力...
-非線形解析と動的解析を組み合わせた解析がまだ１回も回った...
-動的の精度の問題
（以上、[[秋山さんの修論日誌:https://www.str.ce.akita-u.a...
***簡易モデル(Marc編) [#ybe00db2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑簡易モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑桁断面
簡易モデルを作成して、梁の先端変位を簡易モデルのケーブル...
ケーブルはピン接合、梁は原点を完全固定、ケーブルとの節点...
桁の長さは1000mm,完全固定から図の上のピン接合までの距離は...
材料特性は以下に示す。~
,,ヤング率(GPa),ポアソン比,密度(ton/mm^3)
,梁,200,0.3,7.9e-9
,ケーブル,195,0.3,7.9e-9
ケーブルに掛けたプレストレスは-100℃,熱膨張係数は1.25e-5
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ケーブル破断前の先端変位グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ケーブル破断後の先端変位グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル破断後の先端変位グラフ拡大図
先端変位の理論値(片持ち梁のたわみ式から)を導出すると、15....
破断後の最大変位が理論値の1.5〜2倍程度になればいいみたい...
ちなみに、ケーブル破断後に振動したグラフができたのは、梁...
***簡易モデル(Salome-meca編) [#ga6e819b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ビーム要素で作成し、モデルの向きも合わせた。 また、メッ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑梁断面(Marcと同じ) 写真には載せてないが、ケーブル断面も...
モデルの作り方は桁とケーブルのモデルを作成し、メッシュを...
***コマンドファイル [#uc67854a]
・[[ケーブル張力を得る前の設定(DYNA_LINE.ver):http://www....
***Salome-Mecaでできること [#xc3c201b]
・動的応答解析は？~
簡易モデルではDYNA_LINE,DYNA_NON_LINEを使って解析する。(...
・弾塑性解析は?~
再現できた。しかし、降伏応力より大きめの値を設定すると収...
・プレストレスは?~
線形静的解析の場合、温度からでも導入できる。しかし、弾塑...
・途中で張力を落とす(0にする)やり方~
ケーブル破断した際、張力は0になるので、それを再現するため...
***ある時間に力(荷重)を0にする方法 [#k7a3caa6]
ケーブル破断後の挙動を再現するためにある時間でケーブルの...
本来、まず簡易モデルに荷重をかけて、そのときのケーブル張...
そこで、試しに片持ち梁(10✕10✕1000mm)を作成してやってみた。~
梁の先端に0.1mmのバネ要素の線をつけて、バネの先端(1000,0....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
やり方として、DEFI_FONCTION(関数の定義)でNOM_PARA[INST]に...
表の左側にNOM_PARAで定義した数値(今回はINST)を入れ、表の...
図は張力の設定だが、INST0のとき、0倍の荷重が載荷、INST50...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
次にDEFI_LIST_REEL(データ定義)での設定を行う。~
DEBUTは0にする。これは解析を始めるスタート地点(初期値)で...
INTERVALLEはuntil(JUSQU_Aと同じ意味。横軸の最終値)に100,P...
ちなみに、分割数は(until - DEBUT)/PAS で求められる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
これらの設定を踏まえた上で以下のような結果が得られた。~
10/25 動的解析に直し、集中荷重を100N,張力を50Nに直した結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
70s後に梁の先端が揺れる挙動ができた！~
ただ、変位だけは出せるが、何故か力を導出する設定をすると...
プログラムミスということは、どこか設定がおかしいとは思う...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
10/23 静的解析でやってしまった...~
始めの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
70sの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
100sの状態~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
集中荷重を載荷した場所でのy方向の力の様子~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
張力を掛けた場所でのy方向の力の様子~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
梁の先端変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
集中荷重と張力の状態を見る限り、載荷の仕方は求めたいもの...
しかし、張力が0になると、その衝撃で梁が上に持ち上がり、上...
原因を探ってみたところ、載荷再現は問題なかったが、動的解...
静的解析だと、時間を考慮しない解析になり、時間に依存する...
また、挙動を見る限り、荷重を掛けすぎていると思うので、荷...
***ピン接合のやり方(Salome-meca) [#a7708b00]
ピン接合とは？~
→節点（部材と支点の節点）が自由に回転するように接合する方...
バネ要素を用いた線で再現できる。(AFFE_MODELEのDIS_TR,AFFE...
部材と部材の距離はできる限り短くすること。(下の写真のモデ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
回転が自由にできるピン接合を再現するなら、回転に対する剛...
REPEREはGLOBALでよい。
バネ要素を用いたピン接合は、本来節点を共有して接合するの...
・ばねが一軸方向に作用しているかの確認~
以下の3枚の画像には桁の先端に荷重Fy=-500N掛けている。白色...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑x方向変位に対する剛性のみを緩めた場合~
y軸負の向きに力を掛けることで、ケーブルとバネ要素の線の共...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(1000,-14.5803,0.000175109)~
ケーブル側の共有節点(1000,0.1,0)→(985.425,-14.4802,0.0001...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑y方向変位に対する剛性のみを緩めた場合~
y軸負の向きに力を掛けることによって、バネ要素の線が伸びた...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(1000,-14.5802,0.000175108)~
ケーブル側の共有節点(1000,0.1,0)→(1000,0.0953226,0.000175...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑z方向変位に対する剛性のみを緩めた場合~
z軸には影響がほとんど見られなかったため、ほぼ変形せず。(...
変化前と後の座標~
桁側の共有節点(1000,0,0)→(999.989,-0.16365,1.965e-6)~
ケーブル側の共有節点(1000,0.1,0)→(999.984,-0.0587042,1.96...
y方向変位に対する剛性を緩めた画像を見れば分かりやすいと思...
画像を見るかぎりは一軸方向に作用していることが分かる（バ...
***プレストレス [#ld0cda14]
・PRE_EPSIを使ったプレストレス再現(線形静的解析 MECA_STAT...
AFFE_CHAR_MECA → PRE_EPSI を使ってプレストレスを再現して...
PRE_EPSIは変形の負荷を与えるコマンドで、説明書を見るとεと...
[[PRE_EPSIの詳しい説明はこちら:https://code-aster.org/V2/...
今回、ケーブルにプレストレスを適用したいので、EPX（梁のロ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
上記の設定に変えて、解析を行なったら以下の結果が得られた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑scalefactor1倍だと変化が見にくかったので、30倍で見てみた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑桁の先端変位
Marcで得た先端変位：2.49607mm~
Salomeで得た先端変位：2.49608mm~
Marcだと節点に温度応力-100℃与えてプレストレス再現していた...
変形図と先端変位のグラフを見る限り、再現性は高いと思うの...
注意点として、何回か解析を回していると、計算自体は回るが...
・温度で再現する際の問題点~
12/12追記~
CREA_CHAMPとCREA_RESUを使うことで動的解析でも温度を考慮す...
CREA_CHAMPやAFFE_VARCの設定は一切変えずにSTAT_NON_LINEを...
線形静力学を計算してくれる(MECA_STATICS)は、材料の機械的...
説明書やSalome-Mecaの掲示板なども見て、試せるものは全てや...
(現在の設定)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・やり方~
CREA_CHAMPとAFFE_VARCを用いて再現可能。~
過去の先輩方が行なっていた設定を参考にしながらやってみた。~
線形解析で桁に死荷重Fy=-0.0183456N/mmを掛けているが、桁が...
(白線が変形前、色付きが変形後)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Marcでも線形解析でのプレストレスを導入した簡易モデルのデ...
,,Marc,Salome
,桁の先端変位(mm),2.4954,2.49607
SalomeでもMarcとほぼ同じ結果が得られたのでプレストレスの...
後、ピン接合の箇所を確認したところ、変形は見られなかった。~
***軸力の見方 [#df51b2bf]
軸力はEFGE_ELNOの[N]で見る事ができることが分かった。Nは法...
見る方法だが、medファイルをParaViSで閲覧し、フィルター→me...
[[参考資料:https://code-aster.org/V2/doc/v15/en/man_u/u2/...
↑Salomeで出力される力やひずみなどについて詳しく書かれてい...
<模索中のメモ>~
節点力から求める軸力はあまりにも値が大きすぎることから、...
しかし、得られた軸力は4Nほどでケーブル破断後のモデルにそ...
[[PRE_EPSIを使って力を見てみる人が質問していた:https://co...
ただ、回答していた人が[[あるテストケース:https://code-ast...
12/12追記：残念ながら自分の使っているSalomeのversionでは...
***あるtimestepにすると張力が落ちる問題 [#ke57d6cd]
timestepを0.01より小さな値で刻むと安定状態に入る直前に張...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑初期モデルから得た張力(timestepを0.01にしている) 要素分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑99sと100sの間の様子 一瞬だけ張力が落ちているのが分かる。
原因としては温度を用いたプレストレスが問題であることが分...
ただ、PRE_EPSIは設定したひずみになるような外力を与えてい...
そこで、ケーブルを10要素に分割して、1要素（(0,1000,0)側）...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑白色のところにPRE_EPSIを適用した。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑変位も張力も温度で解析していたときとほぼ一致しつつ、張力...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ただ、個人的に気になる点として、0.1Nほどの振れ幅が見られ...
(追記)無視していいそうです。
***10mの簡易モデルについて [#h6683bd4]
1m(1000mm)の桁長で解析を行なったが、思っていた以上の衝撃...
そこで、桁長を10m(10000mm)に変えて解析を行なってみた。(断...
・問題点
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑1m簡易モデルの初期モデルは上記のような挙動が得られていた。
しかし、10m簡易モデルでは下にたわんでしまった。(120sの様子)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル1本モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑ケーブル2本モデル
この状態でケーブル1本モデルの解析を行なってみたところ、以...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ただ、ケーブル2本モデルで解析を行うと,初期モデルのC1ケー...
<対策>~
桁がたわむのは問題ないらしい。まず、1本モデルの結果をまと...
2本モデルに関しては、C2のプレストレスを落として再解析。
**解析結果 [#h7cbc853]
***ケーブル要素なしモデルの説明 [#p1b7be48]
このモデルは初期モデルで得た初期張力を桁の先端に元々ケー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
2024年のデータ~
12/20~
120s後、10s,5s,1s,0.1s,0.05sの時間を掛けて張力を0に落とし...
変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Salomeの結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcの結果
10s〜0.1sまではよさそうな結果が得られているが...0.05sを見...
・問題点~
本来であれば、0.01sとも比較するつもりだったが、初期モデル...
そのため、そのデータをケーブル要素なしモデルに入れても安...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑初期モデルから得た張力(timestepを0.01にしている) 要素分...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑99sと100sの間の様子 一瞬だけ張力が落ちているのが分かる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑得た張力をケーブル要素なしモデルにあてはめた変位
ちなみに、0.005sの場合~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
50sで下がって、100sで上がっている。これ...何が起こってい...
また、0.05s以降から衝撃力が出過ぎな気がするけど...これは...
上フランジ応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
<模索中に得られた結果>~
12/13~
PRE_EPSI、もしくは温度を用いた場合でもケーブル内の応力は...
試しに、ケーブル破断モデルに入れてみたが、桁が下にたわむ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑50sの様子 (scalefactorは100倍)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の拡大図
11/29~
Marcの設定を参考にしながら、Salome-Mecaでプログラムを組ん...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑張力を安定させた直後(50s)のモデルの様子 上がりすぎだろ.....
力を安定させた直後(50s)の桁の先端変位の比較~
,,Marc,Salome
,先端変位(mm),2.508,4321
あまりにもMarcとの差がありすぎるので、試しにMarcで得たビ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Salomeで得た張力データを入れた際の先端変位結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcで得た張力データを入れた際の先端変位結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑Marcの結果
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の振動の拡大図 Salomeに入力した張力データの結果では...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑70s後の振動の拡大図 Marcで作成したグラフの70s後の振動は...
<ケーブル破断後モデルについて、現時点で考えること>~
・当初、Marcの出力結果がkNだと思っていたが、得られた結果...
また、SalomeもMarcも張力が[N]で出ていることは分かったので...
FORC_NODAでビーム軸力を得られるなら、何故値がずれるのか考...
仮にSalomeで正しいビーム軸力を得られなかった場合、Marcの...
⇒張力の見方を改めて考える必要がある。案としては、ケーブル...
それか、MULTIFIBERを用いてケーブルモデルの応力を求め、そ...
・力と時間のグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:FORCの合力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:応力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
横軸:時間 縦軸:応力から力に変換
12/6追記~
色々探ってみたが、バネと桁の影響を無視してケーブルのみのP...
一つ思いついた策として、ケーブルをcable要素に変えてやって...
ただ、問題点として、salomeのケーブル要素はコンクリートに...
また、cable要素は温度を考慮した解析をしてくれることは分か...
・Marcの軸力データをsalomeに入力したとしても、50sまでに桁...
***モデルBについて [#k829fe9a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
ケーブルを2本にして桁を支えたモデルで外側のケーブルを破断...
手順としては初期モデルで外側のケーブル(C1)の張力を取って...
解析時間は150sで100sまでプレストレスと死荷重を掛け、その...
まず、桁の先端変位とC1の張力を取ってみたが、張力はマイナ...
先生に意見を頂いたところ、2本のケーブルのプレストレスが同...
,,変位,張力
,Salome,2.50mm,-32.28N
,Marc,2.51mm,-32.08N
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑先端変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑張力
***解析結果 卒論用 [#z22a073e]
グラフの縦軸は桁の先端変位(mm),横軸は時間(s) 119s〜126s間...
また、死荷重はFy=-0.017248N/mm,集中荷重は桁の先端にFy=-10...
・桁1mモデルA,Bの先端変位の動的応答値とDAF~
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,1.95,1.47
,0.01s,1.93,1.27
,1s,1.02,1.00
,5s,1.01,1.00
・桁10mモデルAの先端変位の動的応答値とDAF（モデルBは健全...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA
,0.005s,2.00
,0.01s,2.00
,1s,1.96
,5s,1.57
ただし、120sまでの健全モデルとケーブル破断モデルで桁の先...
・桁の先端に集中荷重を載荷した桁1mモデルA,Bの先端変位の動...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,1.97,1.47
,0.01s,1.95,1.27
,1s,1.03,1.00
,5s,1.00,1.00
・桁の先端に集中荷重を載荷した桁10mモデルA,Bの先端変位の...
モデルA~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
モデルB~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
DAF~
,時間,モデルA,モデルB
,0.005s,2.00,1.95
,0.01s,2.00,1.95
,1s,1.96,1.26
,5s,1.57,1.04
ただし、モデルAについては120sまでの健全モデルとケーブル破...
○考察~
DAFの値としては,1を超えているため大きな衝撃力が発生してる...
しかし、桁長10mの場合、健全モデルで得た桁の先端変位と120s...
そもそも、健全モデルで得た張力を元々ケーブルがあった箇所...
値が違うということは,何かしら原因があると思うが...。
**Marc Mentatで行う弾塑性解析とSalome_Mecaとの結果比較 [#...
***Salome_Mecaのビーム要素での弾塑性解析の手法 [#afd3f525]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑ビーム要素のみのモデル（片持ち梁のモデルは10✕10✕100mmで...
分かったこと~
・1D要素(ビーム要素)で非線形解析を行う際は、Mulifiber(POU...
[[このページ:https://code-aster.org/forum2/viewtopic.php?...
・断面モデルもGeometryで作成してMeshで切る必要がある。(Mu...
・結果は_ELGAで見ること。NOEU,ELNOでは結果を見ることがで...
・反復回数や分割数(step数)は細かいほうがよい。(載荷値を増...
7/27~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑グラフが折れ曲がる付近の拡大図
[[応力とひずみのデータはこちら:http://www.str.ce.akita-u....
載荷側の断面上の真ん中の球体の応力ひずみデータを取った。
[[ビーム要素の弾塑性解析の設定方法:http://www.str.ce.akit...
DEFI_MATERIAU,DEFI_MATERIAU_GC、どちらでも可。
step数を細かくすると、235MPaで降伏していることが確認でき...
235MPaを超えると、ひずみが大きく増加するので、降伏してい...
応力がずっと同じ値で増える理由だが、恐らく周辺の要素(断面...
ちなみに、降伏後のひずみだが、降伏後の応力ひずみ曲線の傾...
***解析結果 [#bf4e606c]
10✕10✕1000mmの片持ち梁で曲げによる弾塑性解析を行い、235MP...
ビーム要素~
完全固定で設定した節点から応力ひずみ曲線を求めた。
salome-meca~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
Marc~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
ソリッド要素~
Salome-meca~
7/24~
time0.5sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
time0.7sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
time1.0sのとき~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
応力(ミーゼス応力)-ひずみグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑SIEQ_NOEUとEPSI_NOEUを用いた解析結果~
[[グラフの数値はこちら:http://www.str.ce.akita-u.ac.jp/~g...
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
応力(ミーゼス応力)-ひずみ(相当ひずみ)グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑SIEQ_NOEUとEPEQ_NOEUを用いた解析結果~
[[グラフの数値はこちら:http://www.str.ce.akita-u.ac.jp/~g...
見たところは断面真ん中上で(②)、step数は50分割しています。~
SIEQ_NOEUのVMISでミーゼス応力を見れることを知ったので、試...
ミーゼス応力は相当応力とも呼ぶらしく、ひずみも相当ひずみ...
これらが正しいグラフなのか分からないが(そもそも、何故これ...
ちなみに、SIEQ_NOEUではミーゼス応力の他にもトレスカ応力や...
***ミーゼス応力 [#r131c6c9]
後藤メモ(23/9/8)：Salome-Mecaのミーゼス応力は、相当応力$\...
$\overline{\sigma}$は偏差応力のノルムみたいなもので、1軸...
せん断降伏応力を$\tau_{Y}$と書くと、1軸引張の引張降伏応力...
$\overline{\sigma}$の方を相当応力という場合もあり（更には...
この辺は注意が必要である。
-[[構造力学特論テキストにこの辺の話を書きました:https://w...
-[[岩熊先生・小山先生のテキスト参照:http://mechanics.civi...
***相当塑性ひずみ(EPEQ_NOEU) [#d51a592f]
このコマンドはCALC_CHAMP→CRETERIAから設定できる.~
意味は,相当塑性ひずみ（累積塑性ひずみ）で,EPSIから計算さ...
ParavisではINVA_2（ミーゼス相当塑性ひずみ）を使えば,相当...
ちなみに...相当塑性ひずみとは,各成分の塑性ひずみを用いて...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑片持ち梁の断面図
メッシュ長さ0.5 要素数394121でメッシュを切った。下の図は...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
timestepは50分割しており、引張側の真ん中の点が降伏したの...
,,降伏応力(MPa),ひずみ
,①,264.193,0.001099
,②,238.937,0.000947
,③,261.468,0.001107
,⑤,-259.924,-0.001051
,⑥,-248.081,-0.000970
,⑦,-267.69,-0.001079
timestep100の場合(春課題と同じ)~
引張側の真ん中の点が降伏したのがstep40
,,降伏応力(MPa),ひずみ
,①,264.201,0.001099
,②,238.937,0.000947
,③,261.412,0.001107
,⑤,-259.925,-0.001051
,⑥,-248.081,-0.000970
,⑦,-267.692,-0.001051
春課題のデータ
,,降伏応力(MPa),ひずみ
,①,249.267,0.001042
,②,234.025,0.000936
,③,236.751,0.000984
,⑤,-252.657,-0.001037
,⑥,-233.743,-0.000935
,⑦,-252.314,-0.001067
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
Marc
5/22(月)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑断面4隅の点を取った場合の応力-ひずみ曲線。4隅どの1点を取...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑②（引張側）の点を取って、応力-ひずみ曲線を求めた。
引張側の真ん中の点が降伏したincrement14で断面の4隅の状況...
,,降伏応力(MPa),ひずみ
,①③,240.8,0.001451
,②,236.8,0.001186
,⑤⑦,-240.8,-0.001451
,⑥,-236.8,-0.001186
グラフに関してはどうしても上降伏点と下降伏点が現れるグラ...
一応，設定している降伏値に近い値であればよいらしいので，...
しかし，1つ問題点?があり，固定面からずれた位置のほうが大...
②(引張側の真ん中の点)
,位置,降伏応力(MPa),ひずみ
,0mm,236.8,0.001186
,1mm,248.8,0.001116
,2mm,239,0.001091
,3mm,236.7,0.00114
4隅（数値は絶対値で表記）
,位置,降伏応力(MPa),ひずみ
,0mm,240.8,0.001451
,1mm,237.9,0.001304
,2mm,236.6,0.001117
,3mm,238.3,0.001098
5/17(水)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑②（引張側）の点を取って、応力-ひずみ曲線を求めた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑①,②,③の応力-ひずみ曲線をまとめたグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑④,⑧,⑨の応力-ひずみ曲線をまとめたグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
↑⑤,⑥,⑦の応力-ひずみ曲線をまとめたグラフ
・上降伏点と下降伏点が現れるグラフができた。~
グラフとしては間違っていないが、設定した降伏応力を超えな...
***結果比較 [#sa835e43]
ビーム要素~
ソリッド要素~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
,解析ツール,降伏応力(MPa),ひずみ
,Salome,238.937,0.000947
,Marc,236.8,0.001186
4隅のデータ~
Salome-meca
,,降伏応力(MPa),ひずみ
,①,264.193,0.001099
,③,261.468,0.001107
,⑤,-259.924,-0.001051
,⑦,-267.69,-0.001079
Marc(数値は絶対値)
,降伏応力(MPa),ひずみ
,240.8,0.001451
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
Marcのビーム要素でできるグラフ（指定した降伏応力で曲がる...
どちらも降伏応力はそれなりに近い値が出るが、グラフ形状が...
Salomeは指定した降伏応力を設定しても応力がずっと増えてい...
一方、Marcは上降伏点と下降伏点が現れるグラフができるが、...
が、固定面の応力より固定面から離れた位置のほうが少し大き...
また、4隅のデータについてだが、Salomeは4隅のデータそれぞ...
現段階で言えることは、ソリッド要素の弾塑性解析はSalome-me...
(追記：一旦、ビーム要素の弾塑性解析を重視で行っているので...
***ビーム要素での片持ち梁モデルの作り方 [#wedc87f7]
片持ち梁曲げの解析方法は以下の通りに行なった。（一部、[[M...
-長さ単位 (基本操作)~
長さを選ぶ(私はmm)~
-形状とメッシュ (基本操作)~
節点を2点追加 (0 0 0 と100 0 0)(x y z の順，スペース開け...
要素を線(2)に変えて追加 → 節点nodeを入力 → 右クリック確定...
-再分割 (オペレーション)~
分割 → 上からx,y,zの順に値を入力 → 要素 → 分割したいとこ...
-スイープ (オペレーション)~
節点 → 全て → OK (必要のない節点を消す？みたいな)~
-再番号付け (基本操作)~
全ての形状/メッシュ → OK (節点の番号をつけ直す？みたいな)~
形状設定
-新規(構造)▼ (形状特性)~
3ーD ▶ → 中実断面梁(中が詰まっている断面) → 形状選ぶ(長方...
X,Y,Zの値は変えたり変えなかったり(1が橋軸直角方向になるよ...
追加したいところ左クリック(全体選択でも可) → 右クリック確...
テーブルと座標系
-新規(1次独立)~
独立変数V1のタイプをtime(荷重を増分荷重にする)に → データ...
F(縦軸)の数値を入力してenter → データポイントの追加で先ほ...
自分は(0,0),(1,1)と設定した。
ちなみに、バイリニアだとここで応力ひずみの値を入れる必要...
なお、最初の点の塑性ひずみは0でその時の応力が降伏応力であ...
材料設定
-新規(構造)▼ (材料特性)~
有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/...
塑性にチェックを入れて降伏応力を入力(テーブルに応力-ひず...
追加したいところ左クリック(全体選択でも可) → 右クリック確...
境界条件設定
-(固定)~
新規(構造)▼ (境界条件)~
変位指定 → 変位X,Y,Zと回転X,Y,Zにチェックを入れて0にする(...
固定端にする方の節点を左クリック → 右クリック → 矢印出て...
-(載荷)~
新規(構造)▼ (境界条件)~
集中荷重 → 力Yに値を入力し、テーブルにtimeのtableを適用 →...
解析
-新規▼ (解析ジョブ)~
構造 → チェック → 初期荷重をクリックし、境界条件の載荷の...
解析結果をクリック → 要素テンソルのstress,Total Strain,El...
要素スカラーのEquivalent Von Mises Stressにチェック →~
OK → Job check complete になったらOK → 実行 → Marc実行(1)...
結果見る
-ポストファイルを開く(モデルプロット結果処理メニュー)~
→ 履歴プロット → 位置設定で見る点を囲んで右クリック → 全...
X軸データキャリアタイプを位置に → 位置を見たい節点 → 変数...
Y軸データキャリアタイプを位置に → 位置を見たい節点 → 変数...
カーブの追加 → フィット → OK
***ソリッド要素での片持梁モデルの作り方[#xc3ecfbb]
形状設定より先の設定は梁要素と同じなので、それより前の設...
-長さ単位(基本操作)~
基本はmmで行う。
-形状とメッシュ (基本操作)~
節点を2点追加 (0 0 0 と100 0 0)(x y z の順，スペース開け...
要素を線(2)に変えて追加 → 節点nodeを入力 → 右クリック確定...
-拡張 (オペレーション)~
線を拡張してソリッドをつくる。~
「拡張」→回転角（0 0 0）、並進（0 10 0）、繰り返し(1)と入...
回転角（0 0 0）、並進（0 0 10）、繰り返し(1)と入力→「要素...
これで100mm✕10mm✕10mmのソリッド要素が完成
-再分割 (オペレーション)~
ここでメッシュを作成する。~
分割 → 上からx,y,zの順に値を入力（好きな値でいいが、私は1...
→ 要素 → 分割したいところを左クリックで囲って右クリック →...
-スイープ (オペレーション)~
節点 → 全て → OK
形状設定
-新規(構造)▼ (形状特性)~
3ーD ▶ → ソリッド→ 要素を選択 → 追加したいところを左クリ...
ここからの設定は梁要素と同じです。注意点として、境界条件...
（例:点が100個あって、荷重を800N掛けるなら、入力値は8）
**Salome-Mecaの振動解析動画 [#m7a07b3d]
***記録用メモ [#sd5b7728]
5/19(金)遅くなりましたが、Tutorial9,10視聴完了。Tutorial9...
5/4(木)Tutorial8視聴完了。内容は高調波ベース加振を受けるD...
4/30(日)Tutorial7視聴完了。内容は非減衰MDOFシステムの周波...
4/27(木)Tutorial4で使用したモデルを作成し、動画と同じ挙動...
4/26(水)Tutorial6視聴完了。内容はMDOF法（多自由度振動系）...
4/24(月)Tutorial5,5A視聴完了。内容はSDOF法（1自由度系の振...
4/21(金)Tutorial4視聴完了。片持ち梁のモード解析の内容だっ...
4/13(木)Tutorial2B,3視聴完了。データの保存方法(.resuファ...
4/11(火)Tutorial2で使用したモデルを作成し、動画と同じ挙動...
4/8(土)Tutorial2 一通り視聴完了。点と線で片持ち梁モデルを...
*Salome-meca 知識 [#m3ab994c]
ここでは自分が学んだSalome-mecaの知識についてまとめる。ビ...
code-aster名称[[http://opencae.gifu-nct.ac.jp/pukiwiki/in...
**解析手法等の設定まとめ [#bb74bb5a]
***ビーム断面の適用方法 [#g75f9b28]
1D要素で弾塑性解析を始めとした材料非線形の解析を行う場合...
(1)DEFI_GEOM_FIBER~
まず、このコマンドを使う前にビーム要素の断面として使う断...
SECTIONは断面定義、FIBREは断面の穴あき部分を定義する?(例...
ビーム断面が厚肉断面もしくは空白部分に何も入れる必要がな...
TOUT_SECTはビーム要素に作成した平面モデルを全て適用したい...
GROUP_FIBERは好きな名前をつけてよい。~
MAILAGE＿SECTはLIRE_MAILAGEで適用した平面モデルを選択。~
COOR_AXE_POUTREとANGLEは0でよい。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
(2)AFFE_CARA_ELEM~
GEOM_FIBREが適用されていることを確認して、MULTIFIBREの選...
(ただし、AFFE_MODELEで_*Mと書いてあるキーワードを設定する...
(3)DEFI_COMPOR~
ここでは断面の材料特性と状態(弾塑性を入れるとか)を定義す...
このコマンドは画面上の{a}+の横にある{}+(Show All Commands...
スクロールして探すなら、Otherの欄にあるので、それを選択し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
適用すると、下図のような画面が現れるが、GEOM_FIBREはDEFI_...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
MULTIFIBREを選択したら、GROUP_FIBREは(1)と同じ名前,mater...
最後に、DEFI_COMPORをAFFE_MATERIAUのAFFE_COMPORで適用する...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
***動的応答解析のやり方 [#t4e8361f]
・注意点~
(1)DEFI_FONCTIONとAFFE_CHAR_MECAは、Pre Analysis・Analysi...
(2)DISCRETに質量マトリクスを用意しておくこと。DISCRETには...
この設定をしないと、ASSEMBLAGEでバネ要素に質量マトリクス...
[[詳しくはこちらを参照:https://code-aster.org/forum2/view...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
(3)DYNA_VIBRAとCALC_CHAMPの名前は同じにすること。(同じ名...
(4)CALC_CHAMPのEXCITにDYNA_LINEで設定したCHARGEの分だけ入...
設定方法~
[1]Pre analysisのASSEMBLAGEを選択~
・ここでやること~
①NUME_DDLは{}押し、nddlと入力。~
②MATR_ASSEで質量マトリクスと剛性マトリクスを定義する。(質...
③CHARGEはAFFE_CHAR_MECAで適用した境界条件分だけ追加する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2-1]AnalysisのDYNA_LINEを選択~
・ここでやること~
①BASE_CALCULをPHYS,TYPE_CALCULをTRANに設定する。~
②EXCITに載荷条件(PRE_EPSIも同様に)の設定を適用。DEFI_FONC...
③SCHEMA TEMPSでNEWMARK,BETAを0.25と設定する。~
④INCREMENTはDEFI_FONCTIONで作ったテーブルの時間に合わせる...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[2-2]DYNA_NON_LINEの場合~
ここで弾塑性を入れた動的解析ができると思われるが、そのあ...
・最低限やること
①SCHEMA TEMPSでNEWMARK,FORMULATIONをDEPLACEMENT(変位),BET...
ちなみに,[[説明書:https://code-aster.org/V2/doc/v13/en/ma...
更に詳しく調べると,これは直接積分法（正確には直接時間積分...
（Marcもニューマークβ法を用いて連鎖崩壊に関する解析をして...
[[参考サイト:http://jikosoftcom.blog25.fc2.com/blog-entry...
②INCREMENTはDEFI_LIST_INSTで設定したものを選ぶ。~
③EXCITにAFFE_CHAR_MECAで作成した数だけ設定する。載荷条件...
ビーム要素なら、CARA_ELEMにチェックを入れること。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
[3]Post ProcessingのCALC_CHAMPを選択
・ここでやること~
①CALC_CHAMPの名前をDYNA_LINEと同じ名前にすること!これをし...
②EXCITにAFFE_CHAR_MECAで適用した境界条件分だけ追加する。...
***PRE_EPSIについて [#c71fb860]
[[PRE_EPSIの説明書:https://code-aster.org/V2/doc/default/...
モデルをひずみを用いて変形させることができるコマンドであ...
1Dモデル(Multifiber要素も含む)ではEPX（梁のローカル軸に応...
(ky,kzはたわみの微分方程式から得られるdθ/dx = -M/EI の値...
X,Y,Zの軸の向きはグローバル座標(salome画面の色がついてる...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
単純梁モデル(10✕10✕120mm)にEPXX：0.01を与えると、以下の挙...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
線に色がついてるところが変形前のモデルの端であり、写真を...
・簡易モデルに用いた際に分かったこと~
現在、ケーブルにPRE_EPSIを入れており、縮む挙動(ケーブルが...
その理由について考えてみたが、要素を構成している節点にモ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑PRE_EPSIの意味 日本語に訳すと要素に変形の負荷を与えるみ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑簡単な例 PRE_EPSIで与えたひずみになるような外力をかけて...
**Geometry [#i477a78c]
(1)マウス操作について~
マウスホイールを動かす or ctrlキー+マウス左クリック押し続...
ctrlキー+マウスホイールを押し続ける：移動~
ctrlキー+マウス右クリック押し続ける：回転~
(2）つくったモデルを消さずに長さ等を調整する方法~
作成したモデルの長さなどを変えたいとき、わざわざ作り直さ...
Geometryで左上のファイル→スクリプトの保存→ファイルの種類...
すると、下図のような画面が現れるので、点の位置や線の長さ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
**Mesh [#j28db1c4]
(1)メッシュのcompoundについて~
Meshにもcompound機能があり、別々に切ったメッシュを1つにく...
これを使えば、Marcの簡易モデルと同じメッシュの切り方がで...
**Aster_study [#qf331ace]
***Model Definition [#i67b38f6]
・MACR_CARA_POUTREとDEFI_GEOM_FIBERの違い~
どちらのコマンドもビーム要素の断面特性を定義できるが、MAC...
・AFFE_MODELEとAFFE_CARA_ELEMの対応コマンド~
AFFE_MODELEで設定したmodelisationはAFFE_CARA_ELEMのコマン...
対応しているものを以下にまとめたので、参考にしてください。~
POUTRE：POU_D_E, POU_D_T, POU_D_TG, POU_D_T_GD, FLUI_STRU...
BARRE：BAR, CABLE_GAINE~
CABLE：CABLE, CABLE_POULIE~
COQUE：COQUE_AXIS, COQUE_C_PLAN, COQUE_D_PLAN, DKT, DST, ...
DISCRETE：DIS_T, DIS_TR~
DISCRETE_2D：2D_DIS_T, 2D_DIS_TR~
MULTI_FIBRE：POU_D_EM, POU_D_TGM~
・AFFE_MODELE ビーム要素でのMODELISATIONの選び方~
梁の種類によって、選択するものが異なる。~
・真っ直ぐな梁：POU_D_E(オイラー梁), POU_D_T(ティモシェン...
・マルチファイバーの梁(鉄筋コンクリートや非線形材料の解析...
・棒(鉄筋とか断面が円のもの？) 2D_BARRE, BARRE 引張と圧...
・ケーブル CABLE 引張荷重 （引張だけ出せるが、斜張橋など...
・AFFE_CARA_ELEM DISCRETの使い方
DISCRETは離散という意味でバネ要素を定義する場所である。こ...
アルファベットの意味~
A:減衰~
K:剛性~
M:質量~
T：並進運動~
TR:並進運動&回転~
D:対角行列~
L:線~
N:節点~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
↑K_TR_D_Lにした場合の設定画面~
VALE[0],[1],[2]は変位に対する剛性(0,1,2の順にx,y,z方向を...
[3],[4],[5]は回転に対する剛性(3,4,5の順にDRX,DRY,DRZを示...
GLOBALはグローバル座標,LOCALはローカル座標を示す。
***BC & Load [#pcf9b829]
(1)AFFE_CHAR_MECA DDL_IMPO 完全固定の別のやり方
これまで完全固定をする場合、DX=DY=DZ=0と設定していたが、L...
(2)AFFE_CHAR_MECA FORCE_の意味
FORCE_NODAL ノードあたりの力
FORCE_FACE 1面あたりの力
FORCE_ARETE エッジ1本あたりの力(1D要素では使えなかった。)
FORCE_CONTOUR 輪郭の長さあたりの力
FORCE_POUTRE 梁の長さあたりの力
***Analysis [#eb94d3a9]
***Post Processing [#v995178b]
ELGA：ガウス点のelementごとに応力または力またはモーメント
ELNO：node点におけるelementごとの応力または力またはモーメ...
EFGE_ELNO:要素ローカル軸の節点における要素ごとの力または...
・CALC_CHAMPとCREA_CHAMPの違い
CALC_CHAMP を使用して、既存のデータからデータ構造を計算し...
CREA_CHAMP を使用して新しいフィールドを初期化し、後で分析...
CREA_CHAMPでは、構造物にどれだけ温度を与えるか設定できる...
**エラーメッセージ一覧 [#eb4955ab]
・このエラーメッセージは構造的に不安定なときに出る。~
これが出た場合は基本的にAsterstudyの設定を間違えているこ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_CHAR_MECAの設定をモデルに適用できない場合に出る。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
・AFFE_MODELEとAFFE_CARA_ELEMの条件が一致しない場合に出る...
AFFE_MODELEで設定したMODELISATIONは、それぞれAFFE_CARA_EL...
(例) AFFE_MODELEでBARREを選択した場合→AFFE_CARA_ELEMではB...
AFFE_MODELEでPOU_D_Eを選択した場合→AFFE_CARA_ELEMで...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
**Paravis [#zdf3c7e9]
(1)評価項目(magnitudeとか選択するところ)について~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
magnitude：全項目のΣ2乗の平方根(ただし、DEPLだけは変位量...
DEPLのDX：X方向変位量~
SIGM_NOEUのSIXX：X方向垂直応力~
SIGM_NOEUのSIXY：XYせん断応力~
SIEQ_NOEUのVMIS：ミーゼス応力(ちなみに、ミーゼス応力とは...
今載せたのは一例だが、他にも[[このページ:https://qiita.co...
(2).csvファイルの数字とParavisの評価項目~
.csvファイルの1行目の()の数字(例えば,DEPL(0)とか)はPavari...
例を挙げると、SIGM_NOEUの場合、SIXXはSIGM_NOEU(0),SIYYはS...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/...
*TeX [#wcbe1bb5]
久し振りに使うとコマンド等を色々忘れているので、ここにメ...
(1)余白の埋め方
入れる場所はタイトルや名前を書くところより上に(こんな感じ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
{0pt}にすると余白がなくなる。図などを貼り付けた際にできる...
\setlength\textfloatsep{0pt} %本文と図の間の余白
\setlength\floatsep{0pt} %図と図の間の余白
\setlength\intextsep{0pt} %本文中の図の余白
\setlength\abovecaptionskip{0pt} %図とキャプション（図や...
(2)画像の貼り方
ファイルの語尾が.pdf_texの場合(inkscapeで編集したグラフな...
.png(スクリーンショットした画像など)の場合、\includegraph...
(3)数式について
数式で表すには$で挟む必要がある。
(例)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
*参考資料等 [#a2f72b97]
・研究活動で使う参考資料
[[code_asterの動画:https://www.youtube.com/@DigvijayPatan...
[[code-asterの説明:http://dugi6514.odns.fr/doc/v15/en/ind...
[[beam要素の弾塑性解析マニュアル:https://www.str.ce.akita...
[[CAE用語辞典:https://www.cybernet.co.jp/ansys/learning/g...
[[Marc Mentat2018メモ:https://www.str.ce.akita-u.ac.jp/cg...
[[例題で学ぶMarc有限要素法解析入門:https://www.lib.akita-...
[[Salome-Mecaのチュートリアル:https://code-aster.org/V2/s...
[[Pavaris wiki:https://www-paraview-org.translate.goog/Wi...
[[Salome-meca beginning with code_aster:https://archives....
・UNIXコマンド一覧 リンク集
[[https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html...
[[https://proengineer.internous.co.jp/content/columnfeatu...
[[https://amorphous.tf.chiba-u.jp/lecture.files/seminar-I...
・弾塑性解析
春課題で行なった弾塑性解析については以下のURLにまとめてあ...
[[https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E6%...
*日記 [#e5e63424]
3/31(日)~
3年生から更新してきた卒論日誌も本日で最終更新です。明日か...
これからの2年間、色々なことが待ち受けていると思いますが、...
そして、これまで鋼橋の研究をされていた先輩方という自分に...
3/22(金)~
卒業式が無事に終了しました。今日で3年生の後期からずっと一...
けど、またいつかどこかで会う機会はあるはずです。そのとき...
3/8(金)~
無事に謝恩式が終了しました。皆さんが協力していただいたお...
挨拶のとき、皆さんの前に立ったとき、緊張もありましたが、...
秋田での生活はまだまだこれからも続きますが、まずは卒論本...
来年度から鋼橋の研究する先輩が誰もいなくなるから、自分が...
2/21（水）~
久々の更新です。卒論発表が無事に終わり、先生方からの評価...
そしたらなんと、学部生の中では1番良い評価をいただいたので...
ただ、自分の研究成果はまだまだ途中段階ではあるし、まだこ...
12/31（日）~
2023年も残りあとわずかです。~
今年も色々ありましたが、楽しく過ごすことができました。~
特に今年は自分のことを色々見直す機会が多く、もっと人とし...
来年は1月から3月にかけて、卒論に東北支部、それに謝恩会の...
2024年はどんなことが待っているのだろうか…~
来年もよろしくお願いしますm(_ _)m
11/28(火)~
11月ももうすぐ終わり...早いな(・・;)~
来月でいよいよ2023年も終わり、中間発表に卒論が近づいてき...
卒業できないということはないと思うけど、このままだと面白...
残りの時間でどこまでやれるか分かりませんが、先生方や先輩...
11/17(金)~
今日、創造工房実習でSalome-Mecaを教えました。~
皆さんの前で立って教えましたが、人にわかりやすく教えるの...
あの説明で、ある程度理解することはできたのだろうか...。~
次は前で教えることはありませんが、反省点もいくつか見つか...
11/5(日)~
10月下旬あたりからずっと考えていたプレストレス問題だが、...
何で線形静的解析でしかできない(対応していない)のか不思議...
...卒論まで残り3ヶ月だし、いい研究成果が発表できるように...
10/24(火)~
家の近くに熊が出没するという中々恐ろしいことがありました...
思えば10月も来週で終わり、時が経つのが早いなと感じます。~
来年の4月から就職する同級生や先輩と過ごすのもあと半年もあ...
10/20(金)~
ようやく家のネット環境が復帰して、早速自宅で日記を書いて...
…さて、今日は構造研に配属された3年生と顔合わせをしました。~
3年生達の姿を見て、去年の自分達もこんな感じだったのかなと...
自分は大学院に進むので、少なくとも1年半は一緒に過ごすこと...
よろしくお願いしますm(_ _)m
10/18(水)~
本日、3年生の研究室配属先が決まったそうです。~
すでに何人かは顔合わせをしましたが、全員での顔合わせは今...
同時に、自分もしっかりしなきゃいけないとも思います。~
卒業された先輩方や今の院生のようにちゃんとサポートできる...
10/11(水)~
学部生での研究活動も後半戦、他の研究室でもそろそろ卒論に...
いよいよ今週末に3年生も研究室に入ってくるので、悪い手本に...
ただ...最近、細かいミスをしたり、人の話を1回で正しく理解...
あまりこういう状況が続くのはさすがにまずいと思う(場合によ...
9/29(金)~
中間発表が無事に終わり、そして夏休みも終わりました。~
発表に対しても初めて発表をした2月の頃と比べると成長したな...
...少しは変われたかな？~
来週から後期が始まりますので、引き続き研究に取り組みたい...
9/12(火)~
今週は昼過ぎまで集中講義を受講し、夕方から研究をするとい...
calculixでプログラムを作っている最中で、どうしても分から...
それでも簡易モデルの再現ができなければ、どうするべきか相...
9/8(金)~
昨日は帰ってきて何もせずに寝てしまったので、部屋の掃除や...
ひとまず、実家でやっていたデータを研究室のパソコンに移し...
中間発表までは、平日は解析をメインで、休日は概要とスライ...
9/7(木)~
今日、秋田に戻ってきます。~
実家での夏休みは父の実家に行ったり、インターンシップに参...
秋田に戻ったら、再び研究一筋の日々になりますが、卒論に向...
8/6(日)~
秋田は竿燈祭りで盛り上がっています(^^)~
自分は平日に見に行ったので、昨日今日はのんびり過ごしつつ...
後は...台風がどうなるか心配です。無事に帰省できるといいけ...
8/2(水)~
前期のゼミは本日で終わりました。...早いな(^_^;)~
Salomeを使えるノートパソコンを借りることができ、これで実...
ちなみに、明日から竿燈祭りがあるらしいです。人はたくさん...
7/26(水)~
久々の更新です。（4月〜7月の日誌はwikiの整理の際に消して...
今日、外国語文献の講義での発表が終了しました。質問の際、...
けど、4月のときより聞き取る力は上がっていると思うので、英...
気が向いたら、こんな風に更新していこうと思います。
*創造工房実習 課題まとめ [#r8b139fc]
**創造工房実習 発表準備進捗状況記録用メモ（提出済み） [#i...
2/3（金）
課題は完成しました。後はスライド作りをするのみです。（現...
1/31(火)
一通り確認し、後はサンドウィッチ梁の考察について思うとこ...
1/27(金)
概要は提出できる状況までできましたが、提出期限まであと1週...
今日からスライド作成に入りたいと思います。（2/10までには...
1/26(木)
まとめまで一通り書き終えました。後は図の配置や適切な大き...
1/24(火)
異方性2次単純梁の解析完了。サンドウィッチ梁について書く材...
1/23(月)
異方性2次単純梁の解析を行っている途中である。（サンドウィ...
後、LaTeXに貼り付けた図の位置がいまいちなので、調整するこ...
1/20(金)
解析結果と考察を1つにまとめることにした。
残りの作業は解析結果と考察に関してはサンドウィッチ梁を書...
1/17(火)
論文の構成を変更（はじめに→材料性質と解析方法→解析結果→考...
はじめに、材料性質と解析方法を一通り書いた。（図の作成はi...
1/16(月)
創造工房実習の課題の解析結果を一部編集しました。（単純梁...
解析結果の編集に入りました。
1/13(金)
課題で使う全てのグラフをinkscapeで修正しました。
とりあえず、課題の論文の構成は（はじめに→解析結果→考察→ま...
また、論文のタイトルと名前の修正を行い、「はじめに」の部...
1/11(水)
wikiに貼り付けたグラフをsvgに変換させました。
メモ：inkscapeで凡例のみを消す方法はグラフをノードツール...
**12/2(金)課題 サンドウィッチ梁（木材を鋼材で挟む）[#wb04...
異方性の単純梁(1次)と2次要素のサンドウィッチ梁（鋼材：等...
2次要素のサンドウィッチ梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,583289,-,-,千代岡
,0.6,214850,-,-,高井
,0.7,155266,0.0861,-13.0,関合
,0.8,138453,0.083487,-15.7,岡田
,0.9,82766,0.083312,-15.8,松田
,1.2,32279,0.083574,-15.6,青野
,1.3,28343,0.083668,-15.49,山口
,1.4,23667,0.083680,-15.48,山本
,1.5,19958,0.083516,-15.6,進藤
,1.6,19451,0.086037,-13.1,河合
,1.8,10933,0.084022,-15.13,山口
,2,10764,0.083324,-15.8,進藤
,3,3618,0.083497,-15.66,山本
,4,1623,0.0852,-13.9,関合
,5,1007,0.083104,-16.1,千代岡
,6,842,0.0821,-17.1,高井
,7,554,0.080750,-18.4,青野
,8,289,0.079715,-19.5,岡田
,9,261,0.078427,-20.78,松田
,10,232,0.082495,-16.67,河合
異方性1次単純梁と2次要素サンドウィッチ梁の比較グラフ （ih...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：鋼材を挟むのと挟まない場合では理論値が約5倍異なり、...
***おまけ1 異方性単純梁を2次要素にして比較してみる [#qc21...
2次要素の異方性単純梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%)
,0.5,604167,-,-
,0.6,202805,-,-
,0.7,145019,0.516,5.1
,0.8,141896,0.515,4.8
,0.9,91974,0.514,4.6
,1.2,24520,0.512,4.1
,1.3,23810,0.512,4.1
,1.4,17725,0.511,4.0
,1.5,15433,0.511,3.9
,1.6,15746,0.510,3.7
,1.8,11647,0.509,3.5
,2,10460,0.509,3.5
,3,2298,0.504,2.5
,4,1482,0.500,1.6
,5,431,0.496,0.8
,6,356,0.497,1.1
,7,196,0.487,-0.8
,8,104,0.491,-0.14
,9,81,0.485,-1.3
,10,78,0.488,-0.75
（ihou2：11/25の異方性1次解析のデータ ihou2ji:2次要素の異...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
***おまけ2 1次要素のサンドウィッチ梁の場合は？[#gd69b506]
1次要素のサンドウィッチ梁のデータ
,メッシュ長さ,要素数,変位(mm),相対誤差(%)
,0.5,583289,0.08164,-17.5
,0.6,214850,0.07724,-22
,0.7,155276,0.07681,-22.4
,0.8,139020,0.07642,-22.8
,0.9,83327,0.07678,-22.4
,1.2,32182,0.07602,-23.2
,1.3,28507,0.07652,-22.7
,1.4,24105,0.07563,-23.6
,1.5,19934,0.07602,-23.2
,1.6,19653,0.07603,-23.2
,1.8,12468,0.07620,-23
,2,10647,0.07373,-25.5
,3,3599,0.06616,-33.2
,4,1665,0.05568,-43.8
,5,1067,0.05320,-46.2
,6,834,0.03034,-69.4
,7,554,0.01788,-81.9
,8,289,0.01806,-81.8
,9,261,0.01413,-85.7
,10,232,0.03341,-66.3
異方性1次単純梁と1次要素&2次要素サンドウィッチ梁の比較グ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：サンドウィッチ梁を1次要素に変えたらどうなるか気にな...
2次要素のサンドウィッチ梁の解析結果を比較すると、理論値よ...
11/25(金)の課題でも述べたが、2次要素の解析が1次要素より正...
**11/25(金)課題 木材の単純梁の2次要素と異方性1次[#xdf73a23]
＜等方性1次のデータは11/18（金）課題を参照＞
***一次と二次要素の等方性比較 [#j95e0189]
二次要素の等方性のデータ(メッシュ長さ0.5は要素数が多すぎ...
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,-,-,千代岡
,0.6,203209,0.42383,0.98,高井
,0.7,145234,0.43301,3.22,関合
,0.8,140987,0.43006,3.2,岡田
,0.9,91974,0.429913,3.18,松田
,1.2,24800,0.429777,3.14,青野
,1.3,23132,0.42989,3.16,山口
,1.4,17617,0.429745,3.13,山本
,1.5,15433,0.429844,3.2,進藤
,1.6,15900,0.429754,3.13,河合
,1.8,11677,0.42962,3.1,山口
,2,10460,0.429605,3.1,進藤
,3,2486,0.429217,3.0,山本
,4,1453,0.4293,3.02,関合
,5,431,0.427885,2.69,千代岡
,6,360,0.4282,2.78,高井
,7,196,0.42606,2.25,青野
,8,104,0.42631,2.3,岡田
,9,81,0.42513,2.03,松田
,10,78,0.424466,1.8,河合
一次と二次要素の等方性比較のグラフ（niji：二次要素のデー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：[1次要素は頂点の節点だけで要素の変形を表現する、2次...
（引用元：日経XTECH 第2回精度はメッシュで決まる https://x...
↑2次要素の解析での相対誤差が1次要素より小さいのは恐らくこ...
グラフを見ると、1次要素よりも理論値に近い解析ができており...
実際、メッシュ長さ0.5の解析はできなかった。（他のパソコン...
***等方性1次と異方性1次の比較 [#l937a0df]
異方性1次のデータ(メッシュ長さ1.0と1.1追加)
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.50919,3.56,千代岡
,0.6,203209,0.50472,2.6,高井
,0.7,145234,0.5036,2.42,関合
,0.8,140987,0.50283,2.3,岡田
,0.9,91974,0.50053,1.8,松田
,1.0,77582,0.4999,1.6,千代岡
,1.1,27533,0.4886,-0.63,千代岡
,1.2,24800,0.48739,-0.9,青野
,1.3,23132,0.48841,-0.67,山口
,1.4,17617,0.48403,-1.56,山本
,1.5,15433,0.48202,-2.0,進藤
,1.6,15900,0.48329,-1.7,河合
,1.8,11677,0.47855,-2.67,山口
,2,10460,0.47906,-2.6,進藤
,3,2436,0.42787,-12.98,山本
,4,1453,0.42772,-13.02,関合
,5,431,0.27364,-44.3,千代岡
,6,360,0.33927,-31.0,高井
,7,196,0.21363,-58.5,青野
,8,104,0.22574,-54.1,岡田
,9,81,0.22750,-53.7,松田
,10,78,0.20327,-58.7,河合
等方性1次と異方性1次の比較のグラフ（ihou：異方性1次のデー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：解析結果から異方性1次単純梁の変位は約0.50mmに落ち着...
後、等方性1次のグラフと比較すると、お互いグラフの形が似て...
**11/18(金)課題 木材の単純梁 (メッシュ長さ1.0と1.1追加)[#...
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.428982,2.94,千代岡
,0.6,361584,0.421233,1.09,高井
,0.7,145234,0.4225,1.4,関合
,0.8,140987,0.422627385,1.4,岡田
,0.9,91857,0.420351606,0.88,松田
,1.0,77582,0.4192,0.6,千代岡
,1.1,27533,0.4057,-2.6,千代岡
,1.2,24520,0.404744325,-2.87,青野
,1.3,23132,0.4045,-2.93,山口
,1.4,17530,0.3986,-4.34,山本
,1.5,15433,0.396317756757,-4.9,進藤
,1.6,15900,0.399049,-4.24,河合
,1.8,11677,0.404457,-2.94,山口
,2,10460,0.394818715517,-5.3,進藤
,3,2344,0.32447,-22.13,山本
,4,1453,0.3329,-20.1,関合
,5,431,0.136240,-67.3,千代岡
,6,360,0.2130486,-48.9,高井
,7,196,0.1019892,-75.5,青野
,8,104,0.1158624,-72.2,岡田
,9,81,0.1247076,-70.1,松田
,10,78,0.07733,-81.4,河合
単純梁のグラフ：縦軸 変位(mm) 横軸 要素数 （rironti1118...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ：解析結果から、単純梁のグラフの傾きが0に近づくのは約...
**11/11(金) 課題 木材の片持梁[#g9ad393e]
[[コマンドファイルはこちら:http://www.str.ce.akita-u.ac.j...
,メッシュ長さ,要素数,先端変位(mm),相対誤差(%),計算者
,0.5,59504,6.56,-1.5,千代岡
,0.6,45512,6.48774,-2.69,高井
,0.7,39075,6.54133,-2.0,関合
,0.8,13397,6.43695,-3.5,岡田
,0.9,9903,6.36315,-4.6,松田
,1.2,6256,6.3043375,-5.4,青野
,1.3,5767,6.29784,-5.6,山口
,1.4,5146,6.286015,-5.76,山本
,1.5,3935,6.24807,-6.3,進藤
,1.6,3400,6.20446,-6.98,河合
,1.8,2952,6.17161,-7.5,山口
,2,1632,5.64585,-15.3,進藤
,3,667,5.4053975,-18.96,山本
,4,264,3.6161,-45.8,関合
,5,191,3.86,-42,千代岡
,6,190,2.5077325,-62.4,高井
,7,75,1.41225,-78.8,青野
,8,56,1.2887175,-80.7,岡田
,9,49,1.28799,-80.9,松田
,10,44,1.226075,-81.6,河合
片持ち梁のグラフ：縦軸 変位(mm) 横軸 要素数 （katamochi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
メモ:解析結果から、片持ち梁のグラフの傾きが0に近づくのは...
**11/4（金）グラフ作成 [#v961224d]
・課題で作成したグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
・創造工房実習時に作成したグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/...
*3年次の日記(2022/10〜2023/2) [#g01e77dc]
2月10日(金)
今日で発表含め創造工房実習の課題が全て終わりました。ただ...
次は春休み課題（弾塑性解析）を行う必要があるので、それに...
1/31(火)
今日は中村先生の講義に参加しました。
特に印象に残ったのは、ケーブルに乾燥した空気を送りこんで...
ケーブルに関する知識は鋼橋の研究をする際に必要になる知識...
1/11(水)
今日の創造工房実習でLaTeXについて一通り学びました。
後は課題提出に向けて作業するのみなので、頑張っていきたい...
1/9(月)
遅くなりましたが、2023年が始まりました。今年もよろしくお...
明後日、創造工房実習があるので、LaTeXについて昨年学んだこ...
また、創造工房実習の課題の締切まであと約1ヶ月後なので、計...
12/16(金)
LaTeXについて学びました。色々覚えることはありそうですが、...
12/14(水)
これまで行なってきた解析データやグラフを少し整理しました...
12/9(金)
本日は中村先生の特別講義を聴講しました。東京湾アクアライ...
明石海峡大橋のお話は別の機会でしてくださるそうなので、次...
12/8(木)
12/2の課題を終わらせました。明日は特別講義ですが、どんな...
12/5(月)
salome-mecaの単純梁の復習とサンドウィッチ梁の1次要素解析...
2次要素の結果はまだ出ていませんが、結果がでた後に1次とど...
12/3(土）
salome-mecaで行う片持梁解析の復習を行いました。恐らく、来...
12/2(金)
本日で創造工房で習うsalome-mecaの演習が終わりました。今ま...
12/1(木)
11/25の課題を終わらせました。明日はサンドウィッチ梁の解析...
今年も残り1か月です。(月初めから大雪でしたが...(-_-;))
2022年の残りの期間も楽しく過ごしつつ、自分の知識を高めて...
2022年11月の記録 [#o8b95629]
11/28(月)
2次要素メッシュ長さ0.5の解析を違うパソコンを使って行いま...
11/25(金)
今日も先週に引き続きsalome解析を学びました。ただ、2次要素...
また、今週は研究室で色々な方からお話を聴きました。中には...
11/18(金)
salomeで単純梁の解析について学びました。やることは単純梁...
11/11(金)
salomeについて学びました。次回もsalomeを使った解析を行う...
11/9(水)
viiとgnuplotの使い方を軽く復習しました。
11/4（金）
初めて昼休憩での英語会話に参加しましたが、皆さんの会話を...
理解して自分の考えを英語で話すことができるように英語力を...
また、本日の創造工房でgnuplotとviについて学びました。
今後、salomeなど様々なことを学びますが、これらの機能を使...
11/2（水）
UNIXコマンドを実際に使ってみました。昨日調べたものを含め...
これから様々な場面でUNIXコマンドを使う機会があるはずなの...
（ちなみに...ctrl+zで実行しているコマンドを終了させること...
11/4追記：ctrl+zはどうやらコマンドの強制終了らしいので、...
11/1（火)
UNIXコマンドについて調べました。使えそうなコマンドをいく...
＜取り組むこと＞
・UNIXコマンドについて学ぶ。
・構造・材料系の本を読む
2022年10月の記録 [#z5e925c3]
10/28（金）
UNIXコマンドについて学びました。資料に載っているコマンド...
10/27(木)
タッチタイプの一通りの流れをストップウォッチで計測したら...
しかし、キーボードを見ながらの記録なので、できる限り、下...
（追記）木材についての勉強を始めました。（本：プロでも意...
10/24（月）
タッチタイプの練習を行いました。大文字と小文字を交互に打...
10/21(金）
タッチタイプの練習を行いました。まずは、3分以内に打つこと...
＜取り組むこと＞
・タッチタイプの練習
・構造・材料系の本を読む
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
