![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
岡田の卒論日誌に引き続き行う
・日本集成材組合サイズ表 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
スパン中央の応力_歪
接触解析の場合
仮想剤の場合
応力分布
圧縮側
引張側
仮想剤でやっていた挙動を接触解析を使ってやってみる
○参考
・接触と摩擦スライドhttps://www.opencae.or.jp/wp-content/uploads/2019/12/12-contact-friction.pdf
・接触解析スレッドhttps://groups.google.com/g/code_aster-salome-meca/c/8EBQ6Cop5bc
解析メッシュ（コピーして反転）
DISCRETEの場合
CONTINUEの場合
COULOMB＝摩擦係数＝0.001にした場合
等分布載荷
・高さ方向，幅方向にずれている
強制変位
・幅方向のズレがわかりやすい
使用メッシュ
・ズレはメッシュが均等に切れていないためか
・幅方向のズレはなくなった
・接触の設定でMAIT（接触受ける側）とESCL（接触しにいく側）を入れ替えたらズレ方が変わった （今までは左側が下がっていたが，変更後は右側が下がった）
・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 \)
継ぎ手部分
鋼板を互い違いに挿入することで継ぎ手を行う(500mm)
木材のみの継ぎ手モデル
継ぎ手なしモデル
接着剤のヤング率を接着率100%(7500MPa)から低減させる
梁中央鋼板の引張側
載荷荷重を220kNにした
継ぎ手なし
継ぎ手あり
載荷荷重150kN
降伏点周り
梁中央鋼板の圧縮側
降伏点周り
応力分布
接着率100％、45kN載荷時の鋼材下縁の継ぎ手周りの応力分布
梁中央を5000mmとした
継ぎ手1,2
木材のみ
継ぎ手なし
岐阜高専SALOME資料http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2
円柱を押し付けながら角材をスライドさせる
・円柱の上面を押して押し付ける
・強制変位と荷重載荷で押し付けたときに違いがみられた
・強制変位
・荷重載荷
・摩擦あり
・摩擦なし
弾塑性を入れても解析できた
降伏していそうに見える
千代岡くんからParavisでGenerateVectorにチェックを入れるとWarp By Vectorが使えることを教えてもらって変形を見れた
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.
フランジを挿入鋼板の中に埋め込んだ [#j39c9232]
フランジが外に出ているモデルと中に埋め込んだモデルの応力分布 [#h8d27c9f]
ミーゼス応力
軸方向垂直応力
梁中央を0とした応力分布 [#radee60a]
・フランジの切れ目で応力集中していたが、4m積雪時の荷重を載荷しても降伏応力に達していなかった
・これでもフランジの切れ目で降伏応力に達していた
引張側フランジを下から見た応力図 [#v20c251f]
横から見たフランジ切れ目の応力図
梁中央を0とした応力分布
圧縮側
・これでもフランジの切れ目で降伏応力に達していた
引張側フランジを下から見た応力図
横から見たフランジ切れ目の応力図
梁中央を0とした応力分布
圧縮側
フランジを先細らせたフランジ1100mmモデル
フランジの切れ目で降伏応力に達していた
引張側フランジを下から見た応力図
横から見たフランジ切れ目の応力図
梁中央を0とした応力分布
圧縮側
フランジ端部を円形にしたフランジ1100mmモデル
フランジの切れ目で降伏応力に達していた
引張側フランジを下から見た応力図
横から見たフランジ切れ目の応力図
梁中央を0とした応力分布
圧縮側
2025/10/30 ・摩擦あり解析
2025/10/22 ・接触面を細分化した三角形メッシュで解析
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 ・集成材便覧
║ 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 ! ║このアラームの意味を理解していないと、予期せぬ結果を招くことがあります！
解析はできるけどParaVisでこんな感じに名前の前につく
outputのところで名前を設定してやるとParaVisにも反映される
ただこれをやっても上のアラームは消えなかった（メッシュが勝手に作られるのにも関係してそう）
csvファイルの有効数字を決められる（画像では有効数字6桁）
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?PukiWiki/1.4/Manual/Plugin/O-R#content_1_10
