FrontPage

岡田の卒論日誌

研究

岡田の卒論日誌に引き続き行う

鋼板挿入集成材梁

参考

・集成材便覧 https://www.syuseizai.com/databox/doc/%E9%9B%86%E6%88%90%E6%9D%90%E4%BE%BF%E8%A6%A7_%E5%85%A8%E4%BD%93_20240711.pdf

・日本集成材組合サイズ表 http://www.syuseizai.com/ko-item/list.php?k1=0&k3=E75-F240

・集成材について https://www.eu-plan.co.jp/factory/2019/07/09/mokukikaku1/

・新日鉄の木鋼ハイブリット構造 https://www.eng.nipponsteel.com/steelstructures/wp/wp-content/uploads/2021/06/SShybrid_202106.pdf

・新日鉄エンジニアリング技報2013年 https://www.eng.nipponsteel.com/files_publish/page/76/vol04_10.pdf

・鋼材の厚さ https://www.kumagai-steel.co.jp/product/stock/2536/

鋼板挿入集成材梁に関する研究

・支間中央に継手を持つハイブリッド木桁の性能試験及びFEM解析(2011) https://www.jstage.jst.go.jp/article/jscejseee/67/1/67_1_108/_pdf/-char/ja

・鋼板挿入集成材梁のせん断応力に関する実験とFEM解析(2009) https://www.jstage.jst.go.jp/article/structcivil/55A/0/55A_0_932/_pdf

・鋼板挿入集成材梁のせん断性能(2009) https://www.jstage.jst.go.jp/article/structcivil/55A/0/55A_0_925/_article/-char/ja/

・鋼板挿入集成材梁におけるせん断変形評価(2007) http://library.jsce.or.jp/jsce/open/00322/2007/44-01-0018.pdf

・鋼板挿入集成材梁の有限要素弾塑性解析(2007) https://www.jstage.jst.go.jp/article/structcivil/53A/0/53A_0_880/_pdf

・鋼板挿入集成材梁の曲げ解析(2006) http://library.jsce.or.jp/jsce/open/00322/2006/2006-0062.pdf

・鋼材-木材接着の衝撃接着強さに関する基礎的研究(2005) http://library.jsce.or.jp/jsce/open/00035/2005/60-1/60-1-0287.pdf

曲げ挙動を接触解析で再現できないか

仮想剤でやっていた挙動を接触解析を使ってやってみる

摩擦なし接触

harigaikei.png

等分布載荷

・高さ方向,幅方向にずれている

test2yoko.png test2ue.png

強制変位

・幅方向のズレがわかりやすい

heniyoko.png heniue.png

使用メッシュ

・ズレはメッシュが均等に切れていないためか

testMesh.png

四角形メッシュ

sikakuMesh.png

・幅方向のズレはなくなった

SikakuMenkajuMacpyoko.png sikakuMenkajuMacpue.png

・接触の設定でMAIT(接触受ける側)とESCL(接触しにいく側)を入れ替えたらズレ方が変わった (今までは左側が下がっていたが,変更後は右側が下がった)

sikakuEscpyoko.png ・ESCL側が下がるみたい

解析中にでたエラー

・線1本で線固定しようとしたら出たエラーいつもの

<EXCEPTION> <FACTOR_11>                                                                        ║
║ Problème : la matrice est singulière ou presque singulière :                                   ║ 問題:行列が特異またはほぼ特異である:                                   
║   Lors de la factorisation de la matrice, on a rencontré un problème                           ║ 行列の因数分解中に、問題が発生しました。                           ║
║   (pivot nul ou presque nul) à la ligne 62 qui correspond au degré de liberté donné ci-dessus. ║(ピボットがゼロまたはほぼゼロ)行 62 で、これは上記の自由度に対応します。║
║                                                                                                ║
║ Risques et conseils :                                                                          ║リスクとアドバイス:                                                                   
║* Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué.  ║ * その行が物理的な自由度に対応している場合、おそらくは  剛体の動き  が正しく固定されていないことが原因です。 
║  Vérifiez les conditions aux limites.                                                          ║境界条件を確認してください。        
║  Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact.    ║ 接触を行う場合、構造物が接触によってのみ「保持」されてはいけません。  

○固定線を増やしたら解決

・四角形で解析できた設定で三角形メッシュを使って解析しようとしたらでたエラー

<ContactError> <MECANONLINE9_9>                                                                
║ Arrêt par échec dans la boucle de point fixe sur la géométrie.                          ║幾何学上の固定点ループの失敗による停止。                      
║    La base globale est sauvegardée. Elle contient les pas archivés avant l'arrêt.       ║グローバルベースは保存されています。停止前にアーカイブされたステップが含まれています。

○固定方法を面固定から2本の線固定にしたら解析できた

原因不明:線固定でも完全固定にすると上のエラーが出る,収束エラーみたいな計算結果を出すがどこでエラーを出しているかわからない

Instant de calcul:  1.666666666667e-01 - Niveau de découpe: 1
-----------------------------------------------------------------------------------------------------------------------------------------
|     CONTACT    |     NEWTON     |     RESIDU     |     RESIDU     |     OPTION     |     CONTACT    |     CONTACT    |     CONTACT    |
|    BCL. GEOM.  |    ITERATION   |     RELATIF    |     ABSOLU     |   ASSEMBLAGE   |   NEWTON GENE  |    PRESSURE    |     CRITERE    |
|    ITERATION   |                | RESI_GLOB_RELA | RESI_GLOB_MAXI |                |   VARI. CONT.  |    ERROR       |    VALEUR      |
-----------------------------------------------------------------------------------------------------------------------------------------
|     1        X |     0        X | 3.49757E-01  X | 4.93544E+03    |TANGENTE        |    12        X | 4.36303E-02    |                |
|     1        X |     1          | 2.75670E-12    | 3.89000E-08    |TANGENTE        |     0          | 4.08285E-13    | 1.00000E-01    |
-----------------------------------------------------------------------------------------------------------------------------------------
|     2        X |     0        X | 3.49596E-01  X | 4.93527E+03    |TANGENTE        |     0          | 2.09236E-07    |                |
|     2          |     1          | 1.84952E-12    | 2.61098E-08    |TANGENTE        |     0          | 4.78506E-14    | 7.78345E-03    |
-----------------------------------------------------------------------------------------------------------------------------------------
 Critère(s) de convergence atteint(s)
Instant de calcul:  2.000000000000e-01 - Niveau de découpe: 1
-----------------------------------------------------------------------------------------------------------------------------------------
|     CONTACT    |     NEWTON     |     RESIDU     |     RESIDU     |     OPTION     |     CONTACT    |     CONTACT    |     CONTACT    |
|    BCL. GEOM.  |    ITERATION   |     RELATIF    |     ABSOLU     |   ASSEMBLAGE   |   NEWTON GENE  |    PRESSURE    |     CRITERE    |
|    ITERATION   |                | RESI_GLOB_RELA | RESI_GLOB_MAXI |                |   VARI. CONT.  |    ERROR       |    VALEUR      |
-----------------------------------------------------------------------------------------------------------------------------------------
|     1        X |     0        X | 3.49678E-01  X | 5.92254E+03    |TANGENTE        |    12        X | 4.36327E-02    |                |
|     1        X |     1          | 1.22889E-12    | 2.08138E-08    |TANGENTE        |     0          | 3.88745E-13    | 1.00000E-01    |
-----------------------------------------------------------------------------------------------------------------------------------------
|     2        X |     0        X | 3.49498E-01  X | 5.92224E+03    |TANGENTE        |     0          | 9.00728E-07    |                |
|     2        X |     1          | 2.27623E-12    | 3.85708E-08    |TANGENTE        |     0          | 8.01581E-14    | 1.13776E-02    |
-----------------------------------------------------------------------------------------------------------------------------------------
|    10        X |     0        X | 3.49496E-01  X | 5.92219E+03    |TANGENTE        |     0          | 3.89100E-06    |                |
|    10        E |     1          | 1.28049E-12    | 2.16978E-08    |TANGENTE        |     0          | 4.10394E-13    | 1.22578E-02    |
-----------------------------------------------------------------------------------------------------------------------------------------
<Erreur> On dépasse le nombre de boucles de point fixe de géométrie.

・線1本の単純支持にしたら出たエラー

<A> <CONTACT3_98> 
║Contact et détection de collision.                                                            ║ 接触と衝突検出。   
║     Le contact semble "rasant",                                                              ║ 接触は「かすり接触」のように見えます。
║c'est-à-dire que le jeu est nul mais que la pression de contact est quasiment nulle également.║ つまり、遊びはゼロですが、接触圧力もほぼゼロです。
║     -> Risque & Conseil :                                                                    ║  -> リスクとアドバイス: 
║       La découpe automatique du pas de temps ne sera probablement pas efficace.              ║ 自動的な時間ステップの分割は、おそらく効果的ではないでしょう。
║Si vous constatez de fortes oscillations de la vitesse ou du déplacement au point de contact pour cet instant de calcul malgré l'activation du mode de traitement automatique de la collision, ║衝突の自動処理モードを有効にしているにもかかわらず、この計算時点で接触点における速度や変位に大きな変動が見られる場合は、
║vous pouvez essayer d'éviter instant pour lequel le contact "rasant" apparaît de deux manières différentes: ║ 接触が「かすっている」状態が発生する瞬間を、2つの方法で回避することができます。                                              ║
║        - changer la discrétisation en temps initiale                                         ║ - 初期時間離散化を変更する
║        - changer la valeur de SUBD_INST                                                      ║ - SUBD_INST の値を変更する 

○固定線を増やしたら解析できた

互い違い継ぎ手

梁構造

梁断面

集成材の間に鋼板を挟み、接着剤で固定する

集成材:\( 420\times90 \)、鋼板:\( 410\times9 \)、接着剤:\( 410\times1 \)

KousounyHariDanmen.png HariDanmen.png

継ぎ手部分

TagaiTugiteGaikei.png

TugiteYousu.png

鋼板を互い違いに挿入することで継ぎ手を行う(500mm)

 木材のみの継ぎ手モデル

MokuTugiUe.png

 継ぎ手なしモデル

TugiNasiUe.png

接着の不良を接着剤のヤング率の低減でモデル化

接着剤のヤング率を接着率100%(7500MPa)から低減させる

梁中央鋼板の引張側

載荷荷重を220kNにした

TugiteOurHizHipKou220kN.png

TugiNasiOurHizHipKou220kN.png

継ぎ手なし TugiNasiSe750EPSI0.78.png TugiNasiSe750EPSI0.86.png

継ぎ手あり TugiteSe750EPSI0.78.png TugiteSe750EPSI0.86.png

TawaKajuKou220kN.png

載荷荷重150kN

OuHizuHippaP150.png KajuHenHippaP150.png

OuKaHippaP150.png

降伏点周り OuKaHippaP75.png

梁中央鋼板の圧縮側

降伏点周り OuKaAsshuP75.png

応力分布

接着率100%、45kN載荷時の鋼材下縁の継ぎ手周りの応力分布

HipOuBunpuKou45kN.png

梁中央を5000mmとした

継ぎ手1,2

Tugite1Sieq_VMIS_t0.6.png

Tugite2Sieq_VMIS_t0.6.png

木材のみ

MokuTugiSieq_VMIS_t0.6.png

継ぎ手なし

TugiNasiSieq_VMIS_t0.6.png

接触解析

岐阜高専の接触解析を見本に練習

岐阜高専SALOME資料http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2

摩擦あり接触解析

円柱を押し付けながら角材をスライドさせる EnModeru.png

・円柱の上面を押して押し付ける

・強制変位と荷重載荷で押し付けたときに違いがみられた

・強制変位 EnYawaHeni.png ・荷重載荷 EnYawaKaju.png

摩擦ありとなしの比較

・摩擦あり

YawaPenaHeniT1.png YawaPenaHeniT20.png

・摩擦なし

YawaManaHeniT1.png YawaManaHeniT20.png

摩擦なし接触解析

弾塑性を入れても解析できた

降伏していそうに見える DansoChartView.png

千代岡くんからParavisでGenerateVectorにチェックを入れるとWarp By Vectorが使えることを教えてもらって変形を見れた

GenerateVector.png SesshokDanso.gif

SesshokuTest.png TestSita.png

ParavisでWarp By Vectorが使えなかった

メッシュはコンパウンドしないと解析できなかった

・一つのメッシュで解析した時のエラーメッセージ

Contact méthodes maillées.                                                               !
  ! La zone de contact numéro 1 contient 57 noeuds communs aux surfaces maîtres et esclaves. !
  ! Vérifiez la définition de vos surfaces de contact ou bien renseignez un des mots-clés    !
  ! SANS_NOEUD/SANS_GROUP_NO/SANS_MAILLE/SANS_GROUP_MA.
接触メッシュ法 !
  ! 接触ゾーン番号1には、マスターとスレーブの両方のサーフェスに共通する57個のノードが含まれています!
  ! 接触面の定義を確認するか、キーワードを入力してください!
  ! no_node/no_group/no_mesh/no_group_ma.

フランジを用いた継ぎ手

鋼材の板高を480mmから500mmにしたモデル

9K5000LDan.png 9K500H5000LDan.png

フランジを挿入鋼板の中に埋め込んだ [#j39c9232]

5000LFuraBetu.png 480H5000LFuraHitotu.png

 フランジが外に出ているモデルと中に埋め込んだモデルの応力分布 [#h8d27c9f]

ミーゼス応力 9K10F0_70MizeBunp.png 9K10F24_25MizeBunp.png

9K10FMizeOuKa.png

軸方向垂直応力 9K10F0_70SuiBunp.png 9K10F24_25SuiBunp.png

9K10FSuiOuKa.png

 梁中央を0とした応力分布 [#radee60a]

・フランジの切れ目で応力集中していたが、4m積雪時の荷重を載荷しても降伏応力に達していなかった

9K500H5000L89kNBunpu.png 9K500H5000LT0.5AshBunpu.png

フランジを梁と同じ幅にし、長さを5000mmにしたモデル

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/9Kou5000L.png

・これでもフランジの切れ目で降伏応力に達していた

引張側フランジを下から見た応力図 [#v20c251f]

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/5000LXZOuryoku.png https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/9Kou5000Lkire.png

横から見たフランジ切れ目の応力図

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/5000LYZOuryoku.png

梁中央を0とした応力分布

9K5000mmBunpu.png 圧縮側 9K5000mmAshBunpu.png

フランジ形状の比較(矩形、円形、先細り)

9mmH480FuraTanKeiBunpu.png 9mmH480FuraTanKeiOuKa.png

フランジを梁と同じ幅にし、長さを1100mmにしたモデル

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/9K1100L.png

・これでもフランジの切れ目で降伏応力に達していた

引張側フランジを下から見た応力図

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/9K1100LXZOuryoku.png

横から見たフランジ切れ目の応力図

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/9K1100LYZOuryoku.png

梁中央を0とした応力分布 9K1100mmBunpu.png 圧縮側 9K1100mmAshBunpu.png

フランジを先細らせたフランジ1100mmモデル

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/1100LSuboFura.png

フランジの切れ目で降伏応力に達していた

引張側フランジを下から見た応力図 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/1100LSuboXZOuryoku.png

横から見たフランジ切れ目の応力図 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/1100LSuboYZOuryoku.png

梁中央を0とした応力分布 9KSaki1100mmBunpu.png 圧縮側 9KSaki1100mmAshBunpu.png

フランジ端部を円形にしたフランジ1100mmモデル

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/EnFura1100.png

フランジの切れ目で降伏応力に達していた

引張側フランジを下から見た応力図 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/1100LEnXZOuryoku.png

横から見たフランジ切れ目の応力図 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/okada/tugite/B200/Futo/1100LEnYZOuryoku.png

梁中央を0とした応力分布 9KEn1100mmBunpu.png 圧縮側 9KEn1100mmAshBunpu.png

日誌

2025/10/14 ・四角形メッシュで解析

2025/10/10 ・接触解析で曲げ

2025/9/24 ・木材工学特集号投稿

2025/8/25 ・第24回木材工学研究発表会

2025/8/1 ・値を取る点を仮想剤の位置からずらすと応力の変化が継ぎ手なしに近づく

2025/7/22 ・継ぎ手周りの応力分布

2025/7/15 ・載荷荷重220kNで解析 ・降伏後の挙動が過去の研究と似ている?

2025/7/7 ・木材だけ切れたモデルの解析

2025/7/2 ・過去文献

2025/6/27 ・wiki更新

2025/6/23 ・引張側データ抽出

2025/6/20 ・木材提出

2025/5/26 ・MeshのHypothsesを増やしたら勝手にメッシュが作られなくなった

2025/5/22 ・解析でメッシュが勝手に作られる ・パーティションで鋼材を分割

2025/5/19 ・Paravisでcsvをエクスポートしようとすると多分有効数字を聞かれる

2025/5/14 ・ピアサポーター

2025/5/9 ・shaper部材を区別してつくるのが面倒そう

2025/5/5 ・shaperひとつの複雑なオブジェクトを作るのはGeometryより簡単(ボックスとか点を再編集できる)

2025/5/2 ・真ん中の接着層を分ける ・

2025/4/30 ・継ぎ手ありの弾塑性解析 ・降伏しにくい?

2025/4/28 ・継ぎ手ありの線形解析50kN ・salomeでは-ではなく_を使う(解析でエラーになる) ・パーテションを作るときは高度なオプションは付けないほうがいい

2025/4/25 ・継ぎ手なしの線形解析50kN・継ぎ手なしの弾塑性解析100kN ・2021ではメッシュを作ると勝手にグループが作られている

2025/4/23 ・梁モデル

2025/4/21 ・集成材便覧

メモ

salome

ParaVisでDEPLとかに名前がつくアラーム

║ Le nom utilisateur n'a pas été trouvé, on utilise le nom automatique à la place 0000000f.       ║ユーザー名が見つからなかったので、代わりに自動名を使用します 0000000f.      ║
║                                                                                                 ║
║ Ceci peut engendrer des problèmes lors du post-traitement avec PARAVIS.                         ║PARAVISで後処理をする際に問題が発生する可能性があります。 
║ Il est possible de changer manuellement le nom avec NOM_CHAM_MED et NOM_RESU_MED.               ║NOM_CHAM_MED と NOM_RESU_MED を使って、手動で名前を変更することが可能です。
║                                                                                                 ║                                                                          ║
║ Ceci est une alarme.                                                       ║これはアラームです。
║Si vous ne comprenez pas le sens de cette alarme, vous pouvez obtenir des résultats inattendus ! ║このアラームの意味を理解していないと、予期せぬ結果を招くことがあります!  

NameDepl.png 解析はできるけどParaVisでこんな感じに名前の前につく

outputのところで名前を設定してやるとParaVisにも反映される

NOM_CHAM_MED_DEPL.png NOM_CHAM_MED_unnamed.png

NameDeplKai.png

ただこれをやっても上のアラームは消えなかった(メッシュが勝手に作られるのにも関係してそう)

ParaVisでcsvファイルをエクスポートしようとすると出てくる

ParavisExportOption.png csvファイルの有効数字を決められる(画像では有効数字6桁)

Paravisでのgif画像の作り方

https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=read&page=ParaVIS%E3%83%A1%E3%83%A2&word=paravis#yfe85607

wiki上の画像のサイズ変更

https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?PukiWiki/1.4/Manual/Plugin/O-R#content_1_10


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