#author("2022-05-20T12:28:55+09:00","default:kouzouken","kouzouken")
#author("2022-06-23T15:15:37+09:00","default:kouzouken","kouzouken")
#contents
*やること [#ua236cae]
**異種材料(鋼材と木材)の接合モデルについて [#ccb9afdb]
-異種材料(鋼材と木材)の接合モデルについて、有限要素法を使って解析するための手法を検討する
-応力や挙動について、仮想材を用いて実際の実験に近づける方法
**木製治山ダムの解析[#cf825553]
***一面せん断モデル作成 [#m10bf313]
-打ち込み深さ100mm
-  〃   150mm
-  〃   200mm
-上記の結果と実験結果との比較、妥当性の検討
***引き抜き試験モデル作成 [#w3b32ee1]
-打ち込み深さ100mm
-  〃   150mm
-  〃   200mm
-ラグスクリュー、異型棒鋼D16
-上記の結果と実験結果との比較、妥当性の検討
***接合具についての検討 [#adbc776c]
-打ち込み深さ100〜200mm
-ラグスクリュー又は異型棒鋼
-異型棒鋼の場合の径(D16、D14、D12・・・)
***タイトル [#x49ed128]
-鋼棒の木部へのめり込みを考慮した木質構造接合部の数値モデル化(仮)
**木材サイズ、鋼棒径を変えたときの挙動 [#j83481f3]
-□120×120および□105×105
-D16とD13
*改良モデル [#rb29be0a]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokudam2022/kairyou.png
-応力分布

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokudam2022/kairyou2.png

-解析

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokudam2022/kairyou3.png

青が解析 緑と紫は実験結果
*□120×120 D16のときの挙動 [#p85861d8]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokudam2022/kyodou1.png
-よくわからない荷重変位曲線になった。
-ある程度載荷すると、そこからさらに載荷しても木材が動かなくなる(どこか引っかかっている?)
-http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokudam2022/kyodou2.png
-ひずみのコンター図
*めり込みの再現 [#f74cfbae]
-めり込み用の降伏点を仮想材に設定し、ある程度の応力がかかったときに仮想材が降伏→めり込みの再現 と仮定する
-最小主応力(paravis上はPRIN_1)が降伏点に達したかどうかでその部材が降伏したか判断する
-せん断のときの動きを確認するため、ミーゼス応力も確認する
**めり込みの状況 [#n6469e10]
降伏応力2MPaのときの降伏点のひずみは、ヤング率E=7000のとき0.0003になる
-降伏ひずみ0.0003のとき
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/danmen2.png
仮想材の厚さは2mmなので、1mm程度めり込んでいる
-降伏ひずみ0.005のとき(比較用)
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/danmen.png
2mmめりこんでいる
-降伏ひずみを大きくすると、めり込みが大きくなる傾向にある→salome上で降伏時のひずみを調整することで、めり込みの再現につながるのではないか?
-実際のめり込みを確認できれば、それに合わせためり込みを再現するための降伏ひずみを求めることができる?
**仮想材2mmのとき 降伏応力を変えたときの応力ひずみグラフ [#rdbcc373]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/ouryokuhizumi.png
縦軸が木材上面の応力(kN)、横軸がひずみ、紫が降伏応力1N/mm2 緑が2N/mm2 青が5N/mm2 黄が10N/mm2
**最小主応力とミーゼス応力 [#ae48bc9f]
***降伏応力2N/mm2 [#c2e6f046]
-グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/PRIN_1_VMIS.png
-図(最小主応力)

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/prin1.png
-図(ミーゼス応力)

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/vmis.png
-めり込み状況断面図

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/danmen2.png
** [#i426e1cc]
-仮想材料を鋼棒周りに二重に作成(両方0.5mm)し、モデルA、Bを作成
-A・・・やわらかいゴムの仮想材料(内側)、内側に降伏点を入れた木材(等方)の仮想材料(外側)
-B・・・     〃      (外側)        〃             (内側)
**モデルA [#d2970eca]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/merikomi1.png
**モデルB [#h3e0782a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/merikomi2.png

**メモ [#see9716c]
-どうやって「めり込みを再現した」と判断するか
-モデルAはめり込んでいる感じに見えるけどモデルBは鋼棒がすり抜けている?ように見える
-どこの数字を見ればよいか?変位?応力?
*気づいたこととか [#nfe2e2b4]
**Salome-Meca関連メモ [#o7bad890]
***7/2 異方性を持つ材料の弾塑性解析(Salome-Meca) [#q2b2350d]
-やっぱり異方性材料の弾塑性解析はSalome-Meca2019でもサポートしていないっぽいです。(マニュアルをかなり読んだがどこにも書いていなかった)
-圧縮と引張の特性を変えれるDEFI_MATERIAUの設定は見つけることができたけど、何回やってもABNORMAL_ABOUTのエラーが出て解析がうまくいかない。
--「Define a material」→「ECRO_ASYM_LINE」を設定し、「STAT_NON_LINE」→「COMPORTMENT」→「RERATION」のところで「VMIS_ASYM_LIME」を選択するとできる、らしい。まだ調査中
--Salome-Meca usageによると、Salome-Meca2019でECRO_○○系統はバグが発生しやすいとの書き込みがあった。

**Flameshot [#pa530ee2]
-デスクトップの画像をキャプチャーするとき、デフォルトのPrintscreen機能でキャプチャー→gimpで編集、という作業をしていましたが、Flameshotっていうソフトがわりと便利で使いやすいです。

 sudo apt install flameshot

-でインストール後、「アプリケーション」→「take grapical screenshot」で右クリック→「このランチャをパネルに追加」で上に追加しておく。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210604.png

-あとはここをキャプチャーしたいときにクリックして範囲を選択するだけ。
-これのいいところは、簡単な編集が一緒にできます。(上の画像の赤い□とか、つまりわざわざgimpとかで作業する必要がない)



*木製治山ダム 解析 [#z5451548]
**11/2 [#ad0e9c88]
一面せん断試験で、鋼棒周りの仮想材料の設定を
-木材と同じヤング率
-木材の降伏点を入れる

ような感じで、木材にも弾塑性の要素を持たせたモデルを考えています。

木材の降伏点は、木材の基準材料強度(N/mm2)より、めりこみ応力を利用して降伏点を求める。


木材の基準材料強度(めりこみ応力)だと、スギの基準材料強度(N/mm2)は
,材中間部めり込み,材端めり込み,全面圧縮
,6.0,4.8,2.2
となる。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokukouhukuten.png

横軸が変位(mm)、縦軸が荷重(N)

2.2N/mm2だとこのような材料設定になる?
**11/1 [#wa8434db]
EYT理論で理論値を求めるためのエクセルを作成。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20211101_1.png

今後、木材のサイズを変えて解析を行うために、理論値を算出。
□120×120、□105×105、□90×90の3種類(野田さん打ち合わせ資料より)
打ち込み深さは、今の所100mmの場合を考える。
その他、計算に必要な値はとりあえず実験時と同じものを使用する。
,,降伏せん断耐力,降伏モード
,□120×120,4.22,Ⅲb
,□105×105,4.22,Ⅲb
,□90×90,4.22,Ⅲb

ここで、降伏せん断耐力は木材のサイズによらず一定なのではないかと仮定し、木材サイズと打ち込み深さごとのグラフを作成。

打ち込み深さは、実験で行った100,150,200mmと、木材サイズが小さくなることを考慮し比較用で50mmを想定。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/mokuhikaku.png

 縦軸:降伏せん断耐力(N)
 横軸:木材サイズ(□mm×mm)
 moku○○.datの○○は打ち込み深さ

これを見ると、木材のサイズにかかわらず同じ接合具(今回は異型棒鋼D16を想定)、その他の条件が同じの場合、降伏せん断耐力はほぼ変わらないと仮定できる。
今後、木材サイズを変えた一面せん断試験の解析を行う予定。

**10/28 [#r2035c2b]
-せん断試験 100mmモデルを改良。
,,降伏せん断耐力,降伏モード
,理論値,4.23,Ⅲb
,実験値,5.65,Ⅲb
,解析値,5.82,Ⅲb
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20211028_100.png

だいぶ実験値に近づいた。

**8/30 モデルの妥当性 [#cedce76b]
-打ち込み深さ100mm、150mm、200mmのモデルについて、理論値および実験値とSalome-Meca2019で行った解析結果の比較を行った。
***理論値 [#sf2ef3b4]
-ヨーロッパ型降伏理論(EYT理論)による
 P=CFdl
 P:降伏せん断耐力
 C:接合形式定数
 d:接合具の直径(異型棒鋼の場合公称径)
 l:主材厚
 木質構造設計基準・同解説 日本建築学会 P.222~参照

***実験値 [#b116de50]
-先端変位が30mmを超えたときに異型棒鋼にかかる応力を計算したもの
***解析値 [#l4ffa751]
-Salome-Meca2019で行った解析結果
-応力の一番高い箇所の結果を抽出したもの
***降伏モード [#t2fac85a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210830_kouhuku.png
***打ち込み深さ100mm [#w5c9f7be]
,,降伏せん断耐力,降伏モード
,理論値,4.23,Ⅲb
,実験値,5.65,Ⅲb
,解析値,4.91,Ⅲb
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210830_100.png
***打ち込み深さ150mm [#w5c9f7be]
,,降伏せん断耐力,降伏モード
,理論値,4.66,Ⅲb
,実験値,7.55,Ⅳ
,解析値,2.93,Ⅳ
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210830_150.png

-他よりも応力の値は低いが、応力の高い箇所(赤い部分)が広範囲に広がっている→応力が分散している可能性

-解析は間違っていない?


***打ち込み深さ200mm [#w5c9f7be]
,,降伏せん断耐力,降伏モード
,理論値,5.34,Ⅲb
,実験値,8.70,Ⅳ
,解析値,8.50,Ⅳ
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210830_200.png

***追記 打ち込み深さ150mmのモデルについて [#kf9d9742]
-応力が集中する箇所(せん断面近く)を、サブメッシュを利用して細かくメッシュを切ってみた。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20211005_1.png

,,降伏せん断耐力,降伏モード
,理論値,4.66,Ⅲb
,実験値,7.55,Ⅳ
,解析値,4.78,Ⅲa
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20211005_150.png

-こっちのほうが理論値に近く、応力が集中していることが確認できた。

***10/19 追記 [#r2fa8353]
-仮想材料の設定を見直して、少しだけヤング率を上げてみた。
,理論値,4.66,Ⅲb
,実験値,7.55,Ⅳ
,解析値,7.52,Ⅳ
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20211019_150.png

-これが今のところ実験値に一番近い。



**6/21 現在の途中経過 [#c171bd87]
-打ち込み深さ100mm、150mm、200mmのモデルを作成→解析
-100mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210621_100.png
-150mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210621_150.png

-200mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210621_200.png
-野田さんの実験結果との比較を行う。
-現段階で確認できたこと
--打ち込み深さごとの降伏モードのパターンは実験値(論文P49)とほぼ一致
-これから検証すること
--降伏耐力の実験値(論文P46)との比較
--降伏せん断耐力の実験値(論文P49)との比較
--引き抜き試験用モデル作成



**6/11 疑問 [#ja4506b5]
-異型棒鋼の木材接続面の任意の点を取り、横軸:変位,縦軸:荷重(N/mm2)のグラフを作成。


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210611_1.png

-まだ降伏しないので、もっと荷重をかけたほうがよいでしょうか?


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210611_2.png

-降伏モードは野田さんの実験と近いものになっている(写真は打ち込み深さ200mm,異型棒鋼)

**5/14 モデルver.2.01 [#t7d60310]



http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210517_1.png

モデルを改良して木材同士のめり込みを抑えるようにしました。

こちらの方が実際(野田さんの実験)の挙動に近いかと思います。

この先、面圧応力を実験の数値に近づけてみて、どのような動きになるか確かめてみます。

**5/14 モデルver.2作成 [#t7d60310]

モデルver.2を作成。


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210514_1.png


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210514_2.png

木材同士の間とラグスクリュー周りはゴム(仮想材料)を挟んでます。

設定については[[仮想材料:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E4%BB%AE%E6%83%B3%E6%9D%90%E6%96%99]]参照。

前よりよくなりましたが、まだ少し木材同士のめり込みがあるので修正&引き抜き試験に耐えれるかどうかってところでしょうか。







**5/7 後藤さんからの課題3 [#t7d60310]

いろいろなメッシュの細かさで比較してみた結果こうなりました。
(理論値もちょっと間違ってたので修正)

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210507_1.png

**4/27 後藤さんからの課題2 [#t7d60310]

メッシュの粗さを2倍に荒くして比較してみた。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210428_1.png

変位が2.352とのこと。

半分くらいになった。なぜ?

**4/27 後藤さんからの課題 [#t7d60310]
-異型棒鋼のみの片持ち梁みたいな感じにして、手計算と比較してみては?という課題です。



http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210427_1.png


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210427_2.png

計算上はZ変位=4.730に対して解析は最大変位=4.98964でした。

微妙ですがこれくらいの誤差はsalomeではよくあるのでしょうか?

もちろんモデル作成時に設定をミスってる可能性も否定できませんが・・・


**4/26 解析結果と課題 [#t7d60310]
-salomeの本を見たら、BC and loadでどこかに固定しなければうまく計算してくれないらしい。

(実際にはその物体が固定されていなくても、とのこと)

というわけで、仮で下面を固定してみたらうまく計算いきました。


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210426_1.png

次の課題として、

・ラグスクリュー(異型棒鋼)の周りに応力が分散してるように見えない→LIAISON_MAILの設定がおかしい?要確認

・仮想固定とばね固定(実際には固定されていない物体を固定させる)やり方
(下面&接続面をローラー支承にするのが理想)


・モデルはそのままだとして、引き抜き試験をどうすればよいか

そんなところでしょうか・・・


**4/23 解析失敗 [#t7d60310]
-及川さん、ありがとうございます。

とりあえずLiaison-Mailの設定はなんとかなりそうです。

ただ、めちゃくちゃ長いエラーが出てます。
(翻訳サイトに入れたらここがおかしいかも?確認してみましょう!みたいな釈然としないエラーでした。)

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210423_2.png

この辺の材料特性の値をいじればなんとかならないかな、と思ったりしてます。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/20210423_1.png

**及川です。コンパウンドで結合させる方法 [#t7d60310]


***部材をブーリアン結合で結合する方法とコンパウンドで一つのモデルにする方法 [#yf523fd6]
-[[https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ketu1.png]]
-[[ketugou.hdf:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ketugou.hdf]]
-2つの部材を結合するにはジオメトリ上で結合する方法とコンパウンドで部材ごとを一つのモデルとしてAster studyで結合を指定する方法がある。
-新しいエンティティ ービルド ーコンパウンドで一つのモデルにする方法はジオメトリで合わせる部材を選択しコンパウンドを選択する。高度なオプションのグループとサブオブジェクトを継承する。にチェックをする。そのまま適用して閉じる。一つのモデルになる。この時点では結合されていないため、メッシュを切ってそのまま計算かけても部材ごとの結果になるか、失敗する。
-ジオメトリの段階で結合する部材の結合面をグループにしておく。
-Aster study BC and load Liaison-Mail で設定を行う。
-group MA は結合されるソリッド group-MA-Eは結合するソリットの結合面 あとは適用して終了

***詳細設定をしてメッシュをきる方法 [#d90e8e48]
-通常通りメッシュを作成するを選択する。ジオメトリは要素分割したいモデル、メッシュタイプはAny、

3DのアルゴリズムをNETGEN3D > 2DのアルゴリズムをNETGEN1D-2Dにする。

-https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/mesh/3D.png
-https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/mesh/2D.png

2Dの画面の詳細設定の右にある歯車マークをクリックしてNETGEN 2D Paranetersを選択する。

-https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/mesh/syousai.png


最大サイズと最小サイズを設定する。(ジオメトリを選択して詳細設定を開くと最大最小の数値が表示される、最小サイズは選択したモデルのメッシュを切ることができる最小のサイズなためその数値を超えないようにする。)

その設定が終わったらOKで戻って適用して閉じる。あとはいつも通りメッシュを切る。


**4月21日 モデル作成 [#idfc56f8]
□90×90のモデル作成をしました。(とりあえず異型棒鋼はφ16、その周りは□20×20)



*めも [#idfc56f8]
**Liaison_UNIFコマンドについて [#d037e641]
-面の変形について制約をつけることができる
 LIAISON_UNIF=_F(
    DDL=('DX', 'DY', 'DZ'), 
    GROUP_MA=('hippari1', 'hippari2')
  ), 
-上のコマンドについて、「hippari1」「hippari2」の面の変形が固定される(DDLは変化の方向、DX、DY、DZを入力すると全く変化しない)
-Enforce DOFとは違い固定はされないため、例えば大きい構造物の小さな面に力をかけたときにその面が変形するのを防ぐときなどに使う

**野田さん打ち合わせ [#w79a1c86]
-木材

□90×90、105×105、120×120

-接合具

異型棒鋼φ16、14、12
ラグスクリュー

-検討すること

接合具種類、径、先穴、長さ

-とりあえず木材を3つ組み合わせたときの真ん中から分けた状態でモデル作成→解析



**STAT_NON_LINEの設定 [#oc79a2e5]

COMPORTEMENT

RELATION(構成式)=VMIS_ISOT_TRAC(Von Mises降伏条件、等方硬化則、多直線近似)

DEFORMATION(大変形解析)=SIMO_MIEHE(座標変換モデル)、PETIT(微小変形解析)

NEWTON(ニュートン法)

MATRICE=TANGENTE(接線剛性)

REAC_ITER(ニュートン法の接線の傾きを何ステップごとに更新するか)=1(1にしといたほうがいい、デフォルトは0)


CONVERGENCE(収束判定)

RESI_GLOB_RELA(相対誤差)=1.0e-6(デフォルト)

ITER_GLOB_MAXI(最大反復数)=100(デフォルト)



**土木学会認定2級土木技術者 過去問(土木学会誌) [#c2342817]
https://www.jsce.or.jp/journal/student/associate/index.shtml

**バイリニアモデル、トリリニアモデルとはどんなモデルですか? [#f253c18e]
http://www.jssc.or.jp/share/pdf/080703answer.pdf
 構造物や構造要素(部材)に荷重が作用し変形していく過程で示す性質を、復元力特性という。具体的には、弾性域では剛性、非弾性域ではひびわれの発生や各部が順次降伏する等が原因で生じる剛性の低下や最大耐力と最大耐力後の変形による耐力低下の性状などを表す荷重と変形の関係をいう。荷重には、一方向に漸増するものや、地震動のように正負に大きく変化しながら繰返し作用するものもあり、作用する荷重によっても復元力特性は異なってくる。一方向に漸増する荷重に対して構造物や構造要素が示す荷重-変形関係(復元力特性)を、骨格曲線という。また、地震のような繰返し荷重に対して構造物や構造要素する荷重-変形関係(復元力特性)を履歴曲線という。実際に動的解析等に用いられる復元力特性は、図-1に示すように骨格曲線(赤色)と履歴曲線から構成される。図-1の(1)に示すように、荷重と変形の関係を2本の直線(弾性域の剛性と非線形域の剛性の2つで表現)で表したものをバイリニアモデル(bilinear:2本の線の)、図-1の(2)に示すように、荷重と変形の関係を3本の直線(弾性域の剛性と非線形域を2つの領域に分けた2つの剛性の3つで表現)で表したものをトリリニアモデル(trilinear:三つの線の)という。繰返し荷重に対するバイリニアモデルとトリリニアモデルの履歴曲線については、道路橋示方書V耐震設計編,pp.338-339を参照してください。

**プロキシ設定の確認 [#r529a809]
 printenv http_proxy https_proxy

**ウインドウスナップを有効にする [#h3da5b7d]
 メニュー→MATE Tweak→ウインドウ→ウインドウスナップを有効にする

**FEAを用いた鋼部材の応答抽出法に関する検討 [#u4f6486a]
https://www.hit.or.jp/aboutus/pdf/paper_2020_04_steel.pdf
**FORC_NODA [#k2120f76]
-FORC_NODAは「ガウス点の拘束により求めた節点力」らしい。後藤さんに相談したところ、いわゆる節点力という扱いでいいんじゃないか、とのこと。

**gnuplot [#rfac8ebb]
***X軸、Y軸の範囲変更 [#qcc05d73]
 set xrange [x0:x1] 
 set yrange [y0:y1]
***文字の大きさ変更 [#id871d0c]
 set tics font "Arial,20"
-Arialはフォント、その後の数字が文字の大きさ

**libffiとctypesの関係 [#oa018076]
https://notemite.com/python/no-module-named-ctypes/


トップ   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS