小川の修論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*2021 時間割 [#kd898b2c]
**前期 [#j74bf09b]
,曜日,1-2,3-4,5-6,7-8,9-10
,月,,,,ゼミ,
,火,,都市システム計画特論,土質工学特論,,
,水,,外国語文献講読,,,
,木,,材料設計学特論,,,科学技術者倫理特論(後半)
,金,,,,地域防災学特論,
**後期 [#b3222793]
,曜日,1-2,3-4,5-6,7-8,9-10
,月,Introduction to Systems Design Engineering(後半),,,ゼ...
,火,地域エネルギー特論,,,,
,水,,,,地域産業アントレプレナー論(前半),
,木,,,,,
,金,,構造力学特論,水理学特論,,
*やること [#v4ce1bc4]
[[このサイト:http://opencae.gifu-nct.ac.jp/pukiwiki/index...
に書いてある荷重ステップごとのデータを一枚のシートにまと...
*めり込み応力参考文献 [#s6b04714]
https://clta.jp/wp-content/uploads/2013/01/e7b5704203d0fb...
*仮想材料② [#y2e32ea2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Ka2...
引張側鞘管と支柱，支柱とプレート，鋼材の丸棒と鞘管，丸棒...
初期剛性はかなり近くなったが，2段階の降伏は再現できていな...
*支柱基部に仮想材料を入れた場合 [#kefae55c]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/kas...
Cecro3が仮想材料を入れる前，kasouが仮想材料を入れた後
仮想材料を入れても変位はほとんど変わらなかった．
*等値面積置換の種類 [#z4db4868]
http://korankyokai.jp/doc/unified_view_2019.pdf
*変位制御のforc_nodaの合計が荷重を表しているか検討 [#u865...
解析モデルは10*10*100の直方体片持ち梁
**弾性解析 [#cf90e6a5]
,,変位,荷重,相対誤差
,荷重制御,6.44286530701755,100,
,変位制御,6.5,100.8885526,0.89%
,,変位,荷重,相対誤差
,荷重制御,6.44286530701755,100,
,変位制御,6.44286530701755,100.001754,0.001754%
**弾塑性解析 [#w62bde28]
,,変位,荷重,相対誤差
,荷重制御,7.31816723684211,800,
,変位制御,7.3,840.204427,5.03%
**弾塑性解析(防護柵の円柱部分) [#k1c0cd9e]
,,変位,荷重,相対誤差
,荷重制御,225.047847058824,37720,
,変位制御,225,44938,19.14%
**円柱の高さ600mmの部分に仮想材料（ヤング率は円柱の10倍）...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
仮想材料を設けた部分の拡大図（左），荷重制御でのモデル全...
,,変位,荷重,相対誤差
,荷重制御,110.60235,40000,
,変位制御,110.60235,76013.865,90.03%
*変位制御の弾塑性解析 [#h3be51eb]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
載荷線に300mmの強制変位をかけて計算．線全体が300mm変位し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
天端プレートを強制変位させてみた．きれいに曲がりはしたが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
支柱の高さ600mmの位置に円を作り，それを300mm強制変位させ...
*接触解析 [#e17d2286]
[[ボルト締結部のモデリング方法:https://www.fem-vandv.net/...
[[Salome-Mecaを用いた普通ボルト接合の接触解析:http://open...
*ボルトの強度計算 [#r897a809]
[[アンウィンの安全率:https://d-monoweb.com/blog/howto-det...
*LSB接合 [#y6979f15]
実験ではLSB接合とボルト接合
→LSB接合のほうが取り替えが効率的
→佐々木先生の実験では引き抜き耐力が設計荷重を下回る．
→文章でしか書いていないのでこっちでもsalomeで引き抜き力(...
*Cタイプやり直し [#z799495c]
,,設計図,実験
,寸法,橋軸方向960mm 幅員方向960mm,橋軸方向955mm 幅員方向9...
,等級,Mx60($E_{x}=780$MPa $E_{y}=4748.8$MPa $E_{z}=130$MP...
上の表のように設計図のCLTと実際の実験で使ったCLTの寸法，...
**解析結果 [#w86ed6fe]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/C.png
Cecro3が設計図，Cが実験で使った方の寸法に合わせたモデル．
実験で使ったモデルの方が変位が小さくなっている．
等級の違いによりモデル全体の剛性が大きくなったことが原因...
*Cタイプ [#ba1c2e8a]
**接触面 [#v46009de]
接触面のボルト固定はなし。写真のような感じで支柱とチャン...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/IMG...
接触面の幅を[[このサイト:https://gijyutsu-keisan.com/weba...
**解析結果 [#efb3a78a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/IMG...
実験では鞘管部分の2つのボルトが一本ずつ破壊したため、この...
解析ではその部分を一体化させているため、2段階の降伏を再現...
**ECRO_LINE [#y6d6a2fa]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c4....
座標で降伏後の挙動を決めるのではなく、降伏点とその後のグ...
降伏前の挙動はあまり変わらないが、降伏後の変位は大きくな...
降伏後の傾きをECRO_LINEに合わせて解析してみた。グラフはピ...
先週のやり方でも降伏するまでの挙動は一緒のはずなので降伏...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c5....
どちらを使うか
→座標の方を使いたい
→多直線近似など応用が効くから
*弾塑性解析の精度を上げる [#s87a8c92]
東北支部で解析と実験があまり合っていないのではないかとい...
変更点：ひずみを測るグループを長方形のものから線形のもの...
理由：長方形だと破断していない部分も含んでしまい、そこが...
結果：下の図の通り。灰色の線が支部のときで、紫色の線が今...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Bb....
*新しい接合部の提案 [#r8906915]
**鋼製地覆(ブリッジプラスアルファ) [#x29850af]
鋼製地覆の上に防護柵を取り付ければ、車両が衝突しても地覆...
有効幅員が増加するなど他にも利点はあるが、コストや施工方...
https://www.terrearmee.com/others/%E3%80%90%E3%83%96%E3%8...
https://www.sekisuijushi.co.jp/products/koran/product/oth...
*CLTのめり込みについて [#yb9bb658]
https://clta.jp/wp-content/uploads/2013/01/e7b5704203d0fb...
https://clta.jp/wp-content/uploads/2014/10/d71ccd51294e80...
日本CLT協会によると、5層5プライのCLT3つに対してめり込み試...
この3つの平均をとった5.25MPaをめり込み応力として使おうと...
*仮想材料でのめり込み再現 [#x8204799]
-計算はできたがめり込みは確認できない
-原因：仮想材料の材料定数、ボルトと木材の接合部分が邪魔を...
*木材のめり込みについて [#r73b44c9]
https://www.pu-kumamoto.ac.jp/kokenkenkyu/h12gaiyo/12ohas...
http://library.jsce.or.jp/jsce/open/00074/2012/56-01-0039...
Hillの降伏条件を用いて木材のめり込みを再現してる
*CLTのめり込み応力 [#ha557aaa]
https://www.hro.or.jp/list/forest/research/fpri/dayori/20...
ここによると、スギ材のめり込み応力は6MPa
木構造計算基準(8ページ目の表1.7)によると、スギ材の材端の...
短期許容応力度は長期の2倍なので、32kgf/cm2
これを単位換算して3.1MPa
*弾塑性解析(CLTに降伏点を入れない) [#eeecb6b0]
CLTのDefine a materialは直交異方性の設定
STAT_NON_LINE:COMPORTEMENTのGroup of elementに弾塑性にし...
RESI_CPLAN_RELA, ITER_CPLAN_MAXI, RESI_INTE_RELA, ITER_IN...
これでRunすると鋼材と仮想材料は弾塑性で、CLTは弾性で計算...
**計算結果 [#yceec017]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/clt...
r:直交異方性の両端
k:直交異方性のスパン中央
Bs:等方性の両端
Bt:等方性のスパン中央
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/bou...
h:直交異方性のときの防護柵のひずみ
b:等方性のときの防護柵のひずみ
載荷部の変位
等方性：140mm程度
異方性：145mm程度
*CLTのポアソン比 [#a45bf966]
弱軸方向x, 強軸方向y, 板厚方向zとすると、$E_{x}=0.78$GPa,...
\[\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
\end{array} \right)\]
に上記のヤング率と$\nu_{xy}=\nu_{yz}=\nu_{xz}=0.016$を代...
\[\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
\left[ \begin{array}{ccc}
\frac{1}{0.78} & -\frac{0.016}{0.78} & -\frac{0.016}{0.78...
-\frac{\nu_{yx}}{4.75} & \frac{1}{4.75} & -\frac{0.016}{4...
-\frac{\nu_{zx}}{0.13} & -\frac{\nu_{zy}}{0.13} & \frac{1...
\left( \begin{array}{c} \sigma_{x}\\ \sigma_{y}\\ \sigma_...
\]
よって、
$-\frac{0.016}{0.78}=-\frac{\nu_{yx}}{4.75}$より、$\nu_{y...
$-\frac{0.016}{0.78}=-\frac{\nu_{zx}}{0.13}$より、$\nu_{z...
$-\frac{0.016}{4.75}=-\frac{\nu_{zy}}{0.13}$より、$\nu_{z...
したがって、salomeで与えるポアソン比は全て0.016で良いと考...
*CLTのヤング率(Shear Analogy Method) [#oaabe9f5]
*CLTのせん断弾性係数 [#je30522a]
$G_{強}=\frac{E_{弱}}{15}$
$G_{弱}=\frac{E_{強}}{15}$
$G_{厚}=\frac{G_{強}+G_{弱}}{2}$
*補剛材なし [#zce21b6d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/hog...
**geometry [#p64aca8a]
-cut:自己交差の検出にチェック
-fuse:自己交差の検出、不要なエッジを削除にチェック
-partition:下側のタイプの形状を保持、サブ図形の交差部はあ...
**解析結果 [#b5b5f57a]
防護柵基部のひずみを、Bタイプの解析結果と共にグラフに示し...
hizumiが補剛材なし、dansoseiがBタイプ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/hiz...
CLT床版の変位を、Bタイプの床版の結果と共にグラフに示した。
ryoutan, tyuuouが補剛材なし、Bs, BtがBタイプ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/syo...
プレート縁部付近の最小主応力は15.1MPaだった。Bタイプの場...
**やり直し [#h492f358]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/b.png
bがBタイプ、hizumiが補剛材なし。
誤差の原因としては要素数の違いが考えられる。
2万くらいの差があった。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bty...
やり直し前：Aタイプとの比較のためにAタイプと同じ位置にも...
やり直し後：補剛材なしのときとグループの分け方を一緒にし...
*vi [#g5d0c523]
ファイル全体の置換
:%s/置換対象文字列/置換後文字列/g
*gnuplot [#p513dda8]
グラフの色を変えるコマンド
lc '色の名前'
https://yutarine.blogspot.com/2018/12/gnuplot-colorname...
色の名前一覧
http://ayapin-film.sakura.ne.jp/Gnuplot/Primer/Misc/col...
線の太さを変える
lw(linewidth) 5（何倍の太さにするかの数値）
http://graph.pc-physics.com/linethickness.html
目盛りの文字の大きさを変える
set tics font "Arial,20"("フォント,サイズ")
https://mio.yokohama/?p=103
*実験 [#pe23a59f]
http://library.jsce.or.jp/jsce/open/00035/2000/55-01-a/55...
sm400の引張試験．形状は違うが，8.5%のひずみで破壊するもの...
https://www.jstage.jst.go.jp/article/jscejseee/68/2/68_41...
https://www.jstage.jst.go.jp/article/jscejseee/71/2/71_18...
防護柵の静荷重試験が行われていた．防護柵の形状や載荷方法...
https://www.jstage.jst.go.jp/article/jscejseee/70/3/70_38...
選定グラフ，ここから引用
*CLT床版に降伏点を設ける [#g3daafbf]
課題：今までの弾塑性解析ではCLTを等方性として解析していた...
CLTのめり込み応力3.1MPaを降伏点として解析
→Marcではこのやり方だったのでSalomeでも適用できないかと考...
下図はCLT床版の荷重とひずみの関係を表した図である．
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/kou...
降伏点を下げたことにより，実験値よりも変位が大きくなって...
最小主応力は14.6MPa. 等方性のときは17.3MPaだったので少し...
防護柵のひずみを比較した．
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/sak...
防護柵の応力ひずみ曲線にはほとんど変化がない．
*弾塑性解析 [#p7a9087b]
**解析結果 [#t8345eff]
Aタイプ
[[Atype.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/sal...
Bタイプ
[[Btype.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/sal...
hdfファイルを開くにはそれぞれAtype_FilesとAtype.med，Btyp...
gFTPのhttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/oga...
**結果の比較(Aタイプ) [#t067714d]
実験値との比較
実験では降伏後すぐに載荷器具が破壊してしまったため，ひず...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Ahi...
**結果の比較(Bタイプ) [#y5d787ac]
実験値と比較した．比較的実験値と近い曲線を描いているのが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bhi...
**材料定数(応力ひずみ曲線)の設定 [#x871830f]
https://www.jstage.jst.go.jp/article/structcivil/56A/0/56...
http://library.jsce.or.jp/jsce/open/00061/2014/41-01-0057...
http://library.jsce.or.jp/jsce/open/00127/1987/33-0141.pdf
一個目の論文の式(5)と表2のパラメータを用いて下図のような...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bko...
*修論日誌 [#ne83ea39]
**salome20 [#nbce8658]
hdfファイルを開こうとするとエラー
OpenGl_Window::CreateWindow: glXCreateContext failed.
原因は分からない．salome19では起こらない．
**salome19 [#i1cc73ae]
メッシュの切り方
https://qiita.com/Jun_Tatsuno/items/d7246abcebaff3e99b72
終了行:
#contents
*2021 時間割 [#kd898b2c]
**前期 [#j74bf09b]
,曜日,1-2,3-4,5-6,7-8,9-10
,月,,,,ゼミ,
,火,,都市システム計画特論,土質工学特論,,
,水,,外国語文献講読,,,
,木,,材料設計学特論,,,科学技術者倫理特論(後半)
,金,,,,地域防災学特論,
**後期 [#b3222793]
,曜日,1-2,3-4,5-6,7-8,9-10
,月,Introduction to Systems Design Engineering(後半),,,ゼ...
,火,地域エネルギー特論,,,,
,水,,,,地域産業アントレプレナー論(前半),
,木,,,,,
,金,,構造力学特論,水理学特論,,
*やること [#v4ce1bc4]
[[このサイト:http://opencae.gifu-nct.ac.jp/pukiwiki/index...
に書いてある荷重ステップごとのデータを一枚のシートにまと...
*めり込み応力参考文献 [#s6b04714]
https://clta.jp/wp-content/uploads/2013/01/e7b5704203d0fb...
*仮想材料② [#y2e32ea2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Ka2...
引張側鞘管と支柱，支柱とプレート，鋼材の丸棒と鞘管，丸棒...
初期剛性はかなり近くなったが，2段階の降伏は再現できていな...
*支柱基部に仮想材料を入れた場合 [#kefae55c]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/kas...
Cecro3が仮想材料を入れる前，kasouが仮想材料を入れた後
仮想材料を入れても変位はほとんど変わらなかった．
*等値面積置換の種類 [#z4db4868]
http://korankyokai.jp/doc/unified_view_2019.pdf
*変位制御のforc_nodaの合計が荷重を表しているか検討 [#u865...
解析モデルは10*10*100の直方体片持ち梁
**弾性解析 [#cf90e6a5]
,,変位,荷重,相対誤差
,荷重制御,6.44286530701755,100,
,変位制御,6.5,100.8885526,0.89%
,,変位,荷重,相対誤差
,荷重制御,6.44286530701755,100,
,変位制御,6.44286530701755,100.001754,0.001754%
**弾塑性解析 [#w62bde28]
,,変位,荷重,相対誤差
,荷重制御,7.31816723684211,800,
,変位制御,7.3,840.204427,5.03%
**弾塑性解析(防護柵の円柱部分) [#k1c0cd9e]
,,変位,荷重,相対誤差
,荷重制御,225.047847058824,37720,
,変位制御,225,44938,19.14%
**円柱の高さ600mmの部分に仮想材料（ヤング率は円柱の10倍）...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Stu...
仮想材料を設けた部分の拡大図（左），荷重制御でのモデル全...
,,変位,荷重,相対誤差
,荷重制御,110.60235,40000,
,変位制御,110.60235,76013.865,90.03%
*変位制御の弾塑性解析 [#h3be51eb]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
載荷線に300mmの強制変位をかけて計算．線全体が300mm変位し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
天端プレートを強制変位させてみた．きれいに曲がりはしたが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/hen...
支柱の高さ600mmの位置に円を作り，それを300mm強制変位させ...
*接触解析 [#e17d2286]
[[ボルト締結部のモデリング方法:https://www.fem-vandv.net/...
[[Salome-Mecaを用いた普通ボルト接合の接触解析:http://open...
*ボルトの強度計算 [#r897a809]
[[アンウィンの安全率:https://d-monoweb.com/blog/howto-det...
*LSB接合 [#y6979f15]
実験ではLSB接合とボルト接合
→LSB接合のほうが取り替えが効率的
→佐々木先生の実験では引き抜き耐力が設計荷重を下回る．
→文章でしか書いていないのでこっちでもsalomeで引き抜き力(...
*Cタイプやり直し [#z799495c]
,,設計図,実験
,寸法,橋軸方向960mm 幅員方向960mm,橋軸方向955mm 幅員方向9...
,等級,Mx60($E_{x}=780$MPa $E_{y}=4748.8$MPa $E_{z}=130$MP...
上の表のように設計図のCLTと実際の実験で使ったCLTの寸法，...
**解析結果 [#w86ed6fe]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/C.png
Cecro3が設計図，Cが実験で使った方の寸法に合わせたモデル．
実験で使ったモデルの方が変位が小さくなっている．
等級の違いによりモデル全体の剛性が大きくなったことが原因...
*Cタイプ [#ba1c2e8a]
**接触面 [#v46009de]
接触面のボルト固定はなし。写真のような感じで支柱とチャン...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/IMG...
接触面の幅を[[このサイト:https://gijyutsu-keisan.com/weba...
**解析結果 [#efb3a78a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/IMG...
実験では鞘管部分の2つのボルトが一本ずつ破壊したため、この...
解析ではその部分を一体化させているため、2段階の降伏を再現...
**ECRO_LINE [#y6d6a2fa]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c4....
座標で降伏後の挙動を決めるのではなく、降伏点とその後のグ...
降伏前の挙動はあまり変わらないが、降伏後の変位は大きくな...
降伏後の傾きをECRO_LINEに合わせて解析してみた。グラフはピ...
先週のやり方でも降伏するまでの挙動は一緒のはずなので降伏...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/c5....
どちらを使うか
→座標の方を使いたい
→多直線近似など応用が効くから
*弾塑性解析の精度を上げる [#s87a8c92]
東北支部で解析と実験があまり合っていないのではないかとい...
変更点：ひずみを測るグループを長方形のものから線形のもの...
理由：長方形だと破断していない部分も含んでしまい、そこが...
結果：下の図の通り。灰色の線が支部のときで、紫色の線が今...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/ogawa/Bb....
*新しい接合部の提案 [#r8906915]
**鋼製地覆(ブリッジプラスアルファ) [#x29850af]
鋼製地覆の上に防護柵を取り付ければ、車両が衝突しても地覆...
有効幅員が増加するなど他にも利点はあるが、コストや施工方...
https://www.terrearmee.com/others/%E3%80%90%E3%83%96%E3%8...
https://www.sekisuijushi.co.jp/products/koran/product/oth...
*CLTのめり込みについて [#yb9bb658]
https://clta.jp/wp-content/uploads/2013/01/e7b5704203d0fb...
https://clta.jp/wp-content/uploads/2014/10/d71ccd51294e80...
日本CLT協会によると、5層5プライのCLT3つに対してめり込み試...
この3つの平均をとった5.25MPaをめり込み応力として使おうと...
*仮想材料でのめり込み再現 [#x8204799]
-計算はできたがめり込みは確認できない
-原因：仮想材料の材料定数、ボルトと木材の接合部分が邪魔を...
*木材のめり込みについて [#r73b44c9]
https://www.pu-kumamoto.ac.jp/kokenkenkyu/h12gaiyo/12ohas...
http://library.jsce.or.jp/jsce/open/00074/2012/56-01-0039...
Hillの降伏条件を用いて木材のめり込みを再現してる
*CLTのめり込み応力 [#ha557aaa]
https://www.hro.or.jp/list/forest/research/fpri/dayori/20...
ここによると、スギ材のめり込み応力は6MPa
木構造計算基準(8ページ目の表1.7)によると、スギ材の材端の...
短期許容応力度は長期の2倍なので、32kgf/cm2
これを単位換算して3.1MPa
*弾塑性解析(CLTに降伏点を入れない) [#eeecb6b0]
CLTのDefine a materialは直交異方性の設定
STAT_NON_LINE:COMPORTEMENTのGroup of elementに弾塑性にし...
RESI_CPLAN_RELA, ITER_CPLAN_MAXI, RESI_INTE_RELA, ITER_IN...
これでRunすると鋼材と仮想材料は弾塑性で、CLTは弾性で計算...
**計算結果 [#yceec017]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/clt...
r:直交異方性の両端
k:直交異方性のスパン中央
Bs:等方性の両端
Bt:等方性のスパン中央
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/bou...
h:直交異方性のときの防護柵のひずみ
b:等方性のときの防護柵のひずみ
載荷部の変位
等方性：140mm程度
異方性：145mm程度
*CLTのポアソン比 [#a45bf966]
弱軸方向x, 強軸方向y, 板厚方向zとすると、$E_{x}=0.78$GPa,...
\[\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
\end{array} \right)\]
に上記のヤング率と$\nu_{xy}=\nu_{yz}=\nu_{xz}=0.016$を代...
\[\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{y...
\left[ \begin{array}{ccc}
\frac{1}{0.78} & -\frac{0.016}{0.78} & -\frac{0.016}{0.78...
-\frac{\nu_{yx}}{4.75} & \frac{1}{4.75} & -\frac{0.016}{4...
-\frac{\nu_{zx}}{0.13} & -\frac{\nu_{zy}}{0.13} & \frac{1...
\left( \begin{array}{c} \sigma_{x}\\ \sigma_{y}\\ \sigma_...
\]
よって、
$-\frac{0.016}{0.78}=-\frac{\nu_{yx}}{4.75}$より、$\nu_{y...
$-\frac{0.016}{0.78}=-\frac{\nu_{zx}}{0.13}$より、$\nu_{z...
$-\frac{0.016}{4.75}=-\frac{\nu_{zy}}{0.13}$より、$\nu_{z...
したがって、salomeで与えるポアソン比は全て0.016で良いと考...
*CLTのヤング率(Shear Analogy Method) [#oaabe9f5]
*CLTのせん断弾性係数 [#je30522a]
$G_{強}=\frac{E_{弱}}{15}$
$G_{弱}=\frac{E_{強}}{15}$
$G_{厚}=\frac{G_{強}+G_{弱}}{2}$
*補剛材なし [#zce21b6d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/hog...
**geometry [#p64aca8a]
-cut:自己交差の検出にチェック
-fuse:自己交差の検出、不要なエッジを削除にチェック
-partition:下側のタイプの形状を保持、サブ図形の交差部はあ...
**解析結果 [#b5b5f57a]
防護柵基部のひずみを、Bタイプの解析結果と共にグラフに示し...
hizumiが補剛材なし、dansoseiがBタイプ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/hiz...
CLT床版の変位を、Bタイプの床版の結果と共にグラフに示した。
ryoutan, tyuuouが補剛材なし、Bs, BtがBタイプ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/syo...
プレート縁部付近の最小主応力は15.1MPaだった。Bタイプの場...
**やり直し [#h492f358]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/b.png
bがBタイプ、hizumiが補剛材なし。
誤差の原因としては要素数の違いが考えられる。
2万くらいの差があった。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bty...
やり直し前：Aタイプとの比較のためにAタイプと同じ位置にも...
やり直し後：補剛材なしのときとグループの分け方を一緒にし...
*vi [#g5d0c523]
ファイル全体の置換
:%s/置換対象文字列/置換後文字列/g
*gnuplot [#p513dda8]
グラフの色を変えるコマンド
lc '色の名前'
https://yutarine.blogspot.com/2018/12/gnuplot-colorname...
色の名前一覧
http://ayapin-film.sakura.ne.jp/Gnuplot/Primer/Misc/col...
線の太さを変える
lw(linewidth) 5（何倍の太さにするかの数値）
http://graph.pc-physics.com/linethickness.html
目盛りの文字の大きさを変える
set tics font "Arial,20"("フォント,サイズ")
https://mio.yokohama/?p=103
*実験 [#pe23a59f]
http://library.jsce.or.jp/jsce/open/00035/2000/55-01-a/55...
sm400の引張試験．形状は違うが，8.5%のひずみで破壊するもの...
https://www.jstage.jst.go.jp/article/jscejseee/68/2/68_41...
https://www.jstage.jst.go.jp/article/jscejseee/71/2/71_18...
防護柵の静荷重試験が行われていた．防護柵の形状や載荷方法...
https://www.jstage.jst.go.jp/article/jscejseee/70/3/70_38...
選定グラフ，ここから引用
*CLT床版に降伏点を設ける [#g3daafbf]
課題：今までの弾塑性解析ではCLTを等方性として解析していた...
CLTのめり込み応力3.1MPaを降伏点として解析
→Marcではこのやり方だったのでSalomeでも適用できないかと考...
下図はCLT床版の荷重とひずみの関係を表した図である．
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/kou...
降伏点を下げたことにより，実験値よりも変位が大きくなって...
最小主応力は14.6MPa. 等方性のときは17.3MPaだったので少し...
防護柵のひずみを比較した．
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/sak...
防護柵の応力ひずみ曲線にはほとんど変化がない．
*弾塑性解析 [#p7a9087b]
**解析結果 [#t8345eff]
Aタイプ
[[Atype.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/sal...
Bタイプ
[[Btype.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/sal...
hdfファイルを開くにはそれぞれAtype_FilesとAtype.med，Btyp...
gFTPのhttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/oga...
**結果の比較(Aタイプ) [#t067714d]
実験値との比較
実験では降伏後すぐに載荷器具が破壊してしまったため，ひず...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Ahi...
**結果の比較(Bタイプ) [#y5d787ac]
実験値と比較した．比較的実験値と近い曲線を描いているのが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bhi...
**材料定数(応力ひずみ曲線)の設定 [#x871830f]
https://www.jstage.jst.go.jp/article/structcivil/56A/0/56...
http://library.jsce.or.jp/jsce/open/00061/2014/41-01-0057...
http://library.jsce.or.jp/jsce/open/00127/1987/33-0141.pdf
一個目の論文の式(5)と表2のパラメータを用いて下図のような...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/ogawa/Bko...
*修論日誌 [#ne83ea39]
**salome20 [#nbce8658]
hdfファイルを開こうとするとエラー
OpenGl_Window::CreateWindow: glXCreateContext failed.
原因は分からない．salome19では起こらない．
**salome19 [#i1cc73ae]
メッシュの切り方
https://qiita.com/Jun_Tatsuno/items/d7246abcebaff3e99b72
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
