近藤の修論日誌 - $雪荷重を受けるプレストレス木箱桁橋の座屈挙動$

目次

• ちょっと気になったこと（2018/02/15．もう，遅い…）
• フランジ付きモデルの材料定数の影響
• フランジ付きモデルの一体化方法の影響
• 対傾構の影響
• まとめ(橋長・桁高・孔を変化させた時)(対傾構を2つを入れてみる．もう期限的に余裕がないので，無理して高次のモードまで確認しない．)
• 桁高500mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• 桁高600mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• 桁高700mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• 桁高800mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• 桁高900mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• 桁高1000mm
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• グラフ
• 鉛直材80mm・斜材80mm
• 鉛直材95mm・斜材95mm
• 鉛直材80mm・斜材80mm(対傾構2つ)
• 鉛直材95mm・斜材95mm(対傾構2つ)
• Salome-Meca2017の座屈解析メモ
• 固有値が算定されるのに結果が見れない場合
• TRI_DIAG(Lanczos法)の精度を上げる方法
• 固有値解析
• SORENSEN
• JACOBI
• TRI_DIAG
• QZ
• 幾何剛性マトリックス
• 載荷方向と座屈荷重の正負
• CalculiXとかMarcとかのこと
• オンサイト木橋の鋼板をSalome上で自動生成するPythonスクリプト
• ねじ山を作りたくない
• 軸方向に色々な荷重が作用した時の柱の座屈
• 片持ち梁の先端に集中荷重が作用した場合（メモ）
• 片持ち梁に等分布荷重が作用した場合（メモ）
• Bessel関数を解いたら(p.97)
• エネルギー法1(p.99)
• エネルギー法2(p.101)
• 片持ち梁に等分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合
• 片持ち梁に三角形分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合
• 片持ち梁に三角形分布荷重（引張方向）と先端に集中荷重（圧縮方向）が作用した場合
• 片持ち梁に台形分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合
• 片持ち梁に台形分布荷重（引張方向）と先端に集中荷重（圧縮方向）が作用した場合
• 両端固定の問題に応用したい場合
• 等分布載荷
• 台形分布載荷
• 直交異方性材料の直応力-直ひずみ関係
• 木材利用2017メモ
• 5月初旬(発表概要)までに
• 後藤さんと相談した概要の流れをメモ
• 6月から8月(発表スライドor論文)までに
• おまけ
• 印刷用の折り紙円筒
• 支間中央のみを取り出して線形座屈解析
• 後藤ちゃちゃ(17/3/19)
• 理工学デザイン
• PC鋼棒の引張+トルク+摩擦
• CLT関連で使えそうな板の文献をメモ
• 板の座屈
• プレストレス木箱桁橋
• 孔無し一体化モデルの1/4解析
• 孔無し非一体化モデルの1/4解析
• 緊張力の低下
• 孔無し一部非一体化モデルの1/4解析(局部座屈)
• オイラー座屈から3点曲げ載荷時の座屈荷重を推定
• オイラー座屈から等分布載荷時の座屈する積雪深を推定
• datファイルのデータ整理
• 弾性床上の座屈
• 鋼材と木材の接触解析
• CLT
• 梁
• 梁(5層7プライ/7層7プライ)の3点曲げ
• 梁(5層7プライ/7層7プライ)の３点曲げを繰り返す(全体が弾塑性体)
• h=240mm(7層8プライ)のときの梁理論との比較
• b=24mm,h=240mm(7層8プライ)のときの1層と8層のFEMでの比較
• 床版
• 供試体2(7層8プライ/接続金物無し)の静的載荷
• 供試体2(7層8プライ/接続金物無し/ゴム有り)の静的載荷
• 供試体2(7層8プライ/接続金物無し)の弾塑性解析
• 供試体1(7層8プライ/接続金物有り/ゴム有り)の静的載荷
• 供試体3(7層8プライ/供試体2の半分)の静的載荷
• $\ell$=2010mm,h=240mm(7層8プライ)のときの1層モデルと8層モデルのFEMでの比較
• 1層モデルと8層モデルのたわみの誤差についての考察
• 1層モデルのポアソン比について考えてみる
• 長方形板の4辺支持の式を計算
• 供試体2(7層8プライ/ボルト有り)の1/4解析
• CLT橋
• 円筒
• 円筒(角筒)
• 折り目を剛結とした回転螺旋折り
• 折り目を剛結とした反転螺旋折り
• 板の要素数を増やす
• 鋼棒のねじれ座屈
• C3D8とC3D6の混在
• C3D8だけ
• B32
• 横ねじれ座屈の非線形解析
• 手法
• 伊藤さんのトラス橋
• CalculiX関連
• unv-->inp
• 線形座屈解析から初期たわみを与える
• frd-->vtk
• paraviewメモ
• コンター図関連(Color Map Editor)
• あるグループのみ表示
• csvに出力
• 動画化
• Salome-Meca関連
• ある断面の軸方向応力分布を描画
• Maximaリンク
• ImageMagick?メモ
• 複数のスクリーンショットを一括でppm(アスキー)に変換
• 定点観測で撮影した連続写真を重ね合わせ
• pdfを変換するときに画質を下げない
• $\LaTeX$メモ
• minipageのキャプション
• 単位
• ダブルクォーテーション
• 画像をわざと表示させない
• $\LaTeX$周りのツール群
• pdfcrop
• Ghostscriptによるpdfのグレースケール化
• おまけ（修論の背表紙）
• その他，便利なコマンドライン
• pdfimages
• grep
• 特定の文字を含むファイルを検出
• 3Dプリンター メモ
• gnuplotメモ
• 日本語のpngを出力する
• Fortran90(gfortran)メモ
• 固有値計算
• LAPACK
• FFTW
• Julia メモ
• Ubuntu上で最新版にする
• パッケージのインストール
• Python メモ
• コメント文にutf-8の日本語を使えるようにする
• Matplotlibで日本語を使えるようにする
• Rustメモ
• rust.vim
• Goメモ
• パッケージのインストール
• フォーマッタ
• Haskell(GHC)メモ
• パッケージのインストール
• METAFONTメモ
• gftodvi
• mktextfm
• Vimメモ
• ある単語の数を数える
• シンタックスハイライトの追加
• USBメモリにバックアップ
• HHKB Pro2
• 講義関連
• LibreOffice?
• 構造力学特論
• 土質工学特論
• TA関連
• 基礎物理学実験
• 土木環境工学実験(アクリルの引張試験)
• 時間割と単位
• 2016年度
• 2016年度前期
• 2016年度後期

$雪荷重を受けるプレストレス木箱桁橋の座屈挙動$

$プレストレス木箱桁橋鋼板部の局部座屈挙動（仮）$

• 桁高やスパンが変わったときに，座屈挙動がどう変わるか？

M1のとき:$汎用有限要素ツールを利用した非線形解析（案）$

ちょっと気になったこと（2018/02/15．もう，遅い…）

フランジ付きモデルの材料定数の影響

• 対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, 鉛直材95mm)

• 対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm, 斜材85mm, 鉛直材95mm)

• やはり，異方性の影響は1割ぐらいある．木材を等方性材料にすると座屈がほとんど1PC鋼棒間のみで生じるかんじで，隣のPC鋼棒間の変位があまり大きくない．座屈時にアンカープレートらへんが斜めに変形するかどうかの違いだろう．

フランジ付きモデルの一体化方法の影響

• 対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, 鉛直材95mm)

• 対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm, 斜材85mm, 鉛直材95mm)

• 一体化方法によって座屈モード(PC鋼棒間をまたぐかどうか)が違い，全体的非一体化モデルの方が座屈荷重が小さくなる．
• フランジなしモデルでは，一体化方法による影響が1割もあったけど，フランジ付きだとそれほど影響がないみたい．
• フランジ付きモデルだと，どう一体化しても，結局すり抜ける．
• フランジなしモデルの座屈モードは，部分的非一体化モデルにすれば，すり抜けがほとんど生じないため，木材利用2017当時は部分的非一体化モデルの方が実現象に近いと考えていた．
• 実は部分的非一体化モデルの方が木材と鋼板の接している面積が大きいので，（境界の要素のアスペクト比を考慮すると）要素数が多くなってしまい計算が大変．
• 全体的非一体化モデルにしても，飛び出しが生じないモデル化の方法（段差を鋼板側に設ける）が分かったので，別に部分的非一体化モデルにする必要もない．
• 木材利用2017当時は，全体的非一体化モデルにすると飛び出しばかりが生じてしまい意味のある解が得られないことが多く，場合によっては部分的非一体化モデルにせざるを得なかった．

対傾構の影響

• 大池のモデルで，対傾構の数を変化させる．
• 1.対傾構なし(0)
• 2.端部・支間中央(2)
• 3.さらに間に1枚入れる(2+1=3つ)
• 4.さらに間に1枚入れる(3+2=5つ)

• 対傾構を入れると，（境界の木材をSub Meshで分割したせいもあって）要素数がかなり増えるので，計算するのに1日1モデルが限界．対傾構2枚までなら1日2モデルも可能か…？
• モデル化を工夫する(onsitz.py)と，対傾構があってもそこまで要素数が増えないので，1日にもっと沢山計算できる．

• 寸法は大池モデルにし，1/4解析を行う．
• たわみ，軸方向応力について共に最大値のみ参考値としてメモをとっておく．

• [考察]
• 対傾構が直接効いてくるのは支点付近の座屈だと思うが，そもそも大池モデルだと支点付近の座屈が4次モード以上であるため，対傾構が増えたところで座屈荷重はさして変わらない．
• しかし，対傾構をまたいで座屈するモードは算出されていないため，現実的ではないが全てのPC鋼棒間に対傾構を設ければ，座屈荷重も上がるのではないかと思われる．
• 対傾構があると断面変形が抑えられるため，最大たわみは当然小さくなる．で，引張側にも力が伝達するので，引張側の軸方向応力は少し大きくなる．

まとめ(橋長・桁高・孔を変化させた時)(対傾構を2つを入れてみる．もう期限的に余裕がないので，無理して高次のモードまで確認しない．)

• とりあえず，橋長については，約10mから約21mmまで(虎毛と大池の橋長の約数になるように)730mm$\times$3=2.190mずつ伸ばす．
• とりあえず，桁高についても，500mmから800mmまで(虎毛と大池の橋長の約数になるように)100mmずつ増やす．
• 鉛直材の挙動を見るために，1000mmまで増やしてもいいのだけど，見栄えがあまりよくないだろうか？
• 鉛直材・垂直材の幅は80・85・90・95mmの組み合わせが16通りもあるが，そんなには計算できないので，ある程度絞る必要がある．
• ということで，一番きついであろう80mmx80mmを始めにやってみた．
• 80mmx80mmとの比較のために，一番変化が分かりやすい95mmx95mmをやってみた．
• CalculiXみたいに夜通し連続して計算できれば，パターン数が多くても何とかならないこともないと思う．ASTKでそれができるかもしれないのだが，複数のメッシュファイルを読み込んで，別々の名前で結果を出力する方法がいまいち分からない．
• ASTKがSalome-Meca2017でなくなってしまったかと思っていたが，AsterStudy? -> Tools -> Plugins -> Salome-Meca -> Run astk で使える．

• 表中の支間中央No.のNo.は，PC鋼棒間を支間中央から順に数えたもの．
• 表中の鉛直材No.のNo.は，鉛直材を端部から順に数えたもの．
• たわみと軸方向応力について，せっかく計算したので何かの足しになればとメモしておく．最大値をメモしているだけで，わざわざ平均値まで求めていない．
• 当たり前だが，橋長が短いほど断面変形が大きく，最大のたわみが梁理論の値から離れる感じ．たわみを断面で平均すれば，お馴染みのティモシェンコ梁の式とかなり合う(概算)．
• 固有値が近接しているため，何次モードまで計算するべきか迷う．せん断が支配的だと逆向きモードが出力されやすいため，スパンが短い場合は次数を多めに計算しておく必要あり．

桁高500mm

鉛直材80mm・斜材80mm

• 橋長10.220mで，いくつか飛び出しモードが出てしまった．

鉛直材95mm・斜材95mm

• 橋長10.220mで，飛び出しモードばかりが出てしまい，なかなか有意な結果が出ない．要求次数を増やすべし!

鉛直材80mm・斜材80mm(対傾構2つ)

• やはり，木材と鋼板の境目で木材側の要素のアスペクト比を小さくできると，（その代わりに鋼板側の要素のアスペクトが大きくなっても，）飛び出しが生じないようだ．

鉛直材95mm・斜材95mm(対傾構2つ)

桁高600mm

鉛直材80mm・斜材80mm

• 橋長10.220mで，飛び出し発生!(3次モードは得られたけど)

鉛直材95mm・斜材95mm

• 橋長10.220mで，飛び出し発生!　というか，このケースのほとんどが飛び出しモード．

鉛直材80mm・斜材80mm(対傾構2つ)

鉛直材95mm・斜材95mm(対傾構2つ)

桁高700mm

鉛直材80mm・斜材80mm

• 橋長10.220mで，いくつか飛び出しモードが出てしまった．
• Sub Mesh の優先順位を1:境界面 2:木材 に変えた方がいいかもしれない．（追記：Sub Mesh の優先順位を変えるとメッシュが切れない（進まないの意）し，要素数が増えすぎる）
• double mesh がどうのこうのと書かれているので，remove double mesh だったかをしておくのも手かな？

鉛直材95mm・斜材95mm

• やはり橋長10.220mは飛び出しが多い．経験的に座屈しにくい形状のものほど飛び出しやすい．

鉛直材80mm・斜材80mm(対傾構2つ)

鉛直材95mm・斜材95mm(対傾構2つ)

桁高800mm

鉛直材80mm・斜材80mm

• ここの表だけ木材の幅を420mmに間違えて(0.1mm削ると419.9mmのはず)計算していました．すみません．でも，確実に3桁を下回る程度の誤差だろうから，やり直すにしてもなるべく後回しにしたい．

鉛直材95mm・斜材95mm

• 橋長10mは飛び出しモードも算出された．

鉛直材80mm・斜材80mm(対傾構2つ)

鉛直材95mm・斜材95mm(対傾構2つ)

桁高900mm

鉛直材80mm・斜材80mm

鉛直材95mm・斜材95mm

鉛直材80mm・斜材80mm(対傾構2つ)

• ちなみに，特に橋長10mなどの短いものは，変位が最大となるところが支間中央ではないので注意．

鉛直材95mm・斜材95mm(対傾構2つ)

• 橋長10mの2次モードは，鋼板全体が波々に座屈する．

桁高1000mm

鉛直材80mm・斜材80mm

鉛直材95mm・斜材95mm

鉛直材80mm・斜材80mm(対傾構2つ)

• ちなみに，橋長10mなどの短いものは，変位が最大となるところは支間中央ではないので注意．

鉛直材95mm・斜材95mm(対傾構2つ)

• 橋長10mの1次モードは，鉛直材というよりも，支点付近の鋼板が全体する．

グラフ

鉛直材80mm・斜材80mm

• 発表や論文には使えないが，ブラウザ上でマウスで動かせたら見やすいのではないかと，（単なる興味で？）plotly.jsを試してみる．↓
• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/graph/plotly/3d_e80s80.html
• もしかすると，色はz軸の値にするのではなくて，座屈箇所によって色分けした方がわかりやすいかもと思い，やってみる．↓
• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/graph/plotly/3d2_e80s80.html
• なんだか，surfaceのx軸とy軸がよく分からないので，mesh3dというグラフタイプでやった．plotly.jsの3次元プロットやJavaScript?の配列を勉強中．

• 設計参考資料として使いやすいように，2次元グラフも作ってみる．

鉛直材95mm・斜材95mm

• ブラウザ上で表示
• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/graph/plotly/3d_e95s95.html
• ブラウザ上で表示(座屈箇所ごとに色分け)
• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/graph/plotly/3d2_e95s95.html

• 設計参考資料として使いやすいように，2次元グラフも作ってみる．

鉛直材80mm・斜材80mm(対傾構2つ)

• 設計参考資料として使いやすいように，2次元グラフも作ってみる．

鉛直材95mm・斜材95mm(対傾構2つ)

• 設計参考資料として使いやすいように，2次元グラフも作ってみる．

Salome-Meca2017の座屈解析メモ

固有値が算定されるのに結果が見れない場合

• 座屈解析の CALC_MODES を編集 -> VERI_MODES にチェックを入れ，Edit を押す -> STOP_ERREUR を'NON'にする
• ただし，計算されたモードは自己責任で評価して下さい，とのことだ．
• また，medファイルが重くて？モデリングや計算に使ったhdfファイルでは開けないので，別のhdfファイルで開いたほうがいい．

TRI_DIAG(Lanczos法)の精度を上げる方法

• 座屈解析の CALC_MODES を編集 -> SOLVEUR_MODAL の Edit を押す -> COEF_DIM_ESPACE にチェックを入れ，値を入力
• COEF_DIM_ESPACE は，デフォルトでは4なので，それを適宜上げる．目安としては，求めたい次数以上にするといいとか．
• 例えば求めたい次数が6のときに，COEF_DIM_ESPACEがデフォルトのままだと，3次で有効数字が3桁，4次で有効数字が2桁，5・6次で有効数字は1桁になるので，高次モードの信頼性がかなり低い！
• ただし，要求次数が6のとき，COEF_DIM_ESPACEをデフォルトの4から6にすると，計算時間も1.5倍に増える．
• と思いきや，要求次数6のとき，試しにCOEF_DIM_ESPACEを12にしてみたら，計算時間がデフォルトの1/10にも減った．
• さらに，COEF_DIM_ESPACEを24にしてみたら，計算時間がデフォルトと同じぐらいになった．
• ということで，精度を上げるためにはCOEF_DIM_ESPACEを大きくしたほうがいいが，それに加えて計算速度を求めるのなら職人的に調度良い値を見つけるといい．経験的には，要求次数の1.5~2.0倍ぐらいかな？

• PREC_LANCZOS(精度(precision)だと思う)を変えても，結果は何も変わらなかった．一体これは何なのか？

• CalculiXでは Lanczos vectors はデフォルトで4$\times$次数らしく，*BUCKLEの次の行の3つ目の値を変えればいい．
• デフォルトで4$\times$次数とは！なかなか多いですね．しかも，デフォルトのCode_Asterよりも倍ぐらい計算が早いのに．

*BUCKLE
6, 0.01, 24, 1000
**次数，精度，Lanczos vectors，最大反復回数…この行はただのコメント

固有値解析

• 昨年度の方法を真似てSORENSEN法？とやらを使っているが，他の（QZ・JACOBI・TRI_DIAG）も試してみたい．
• よくわからないが，アプローチ？を追加してみる．
• REELとIMAGで同じ結果になり，COMPLEXEは結果が出ない．どうやら正負とは関係ない．

SORENSEN

• 後藤さんいわく，何かの間違いで普通でない非対称の行列ができてしまった時に，そうとも知らずに勝手に計算できてしまうのは良くないのではないか，とのことなので基本的に使わないことにする．

JACOBI

• 後藤さんによると，時間がかかっているからといって必ずしも正確なわけでもないらしいので，基本的に使わないことにする．

TRI_DIAG

• 原則として，TRI_DIAGで計算することにしようかと思う．

• TRI_DIAGの精度(PREC_LANCZOS)はデフォルトでは1.0E-08だが，CalculiXの場合1.0E-02である．ちょっと必要精度が高すぎるように思える．
• CaluliXとオンサイト木橋の座屈を比較してみたところ，座屈時の積雪深はCalculiXの方が0.2~0.3%高くなり，座屈モードはほぼ一致した（2次モードまで確認）．ちなみに，CalculiXでは特にエラーや警告は出ていないし，精度を0.01から0.00001にしたところで結果は何一つ変わらないようだ．
• ところで，CalculiXでは荷重を節点数で割っただけで，等価節点力に置き換えていないので，より要素分割の細かい鋼板側に荷重が多めにかかっている筈である．なので，CalculiXの方が座屈荷重が小さくなりそうな気がするのだけど，そうはなっていない．不思議．
• 試しにPREC_LANCZOSを1.0E-02に変えてみても，結果は何一つ変わらなかった．

• messageファイルにCOEF_DIM_ESPACEを上げてみては？と書かれているので，デフォルトの4からちょっとずつ上げていってみようかと思う．
• COEF_DIM_ESPACEとはここでいうnのことだろうか，そうしたら求めたい固有値の数以上の値を与えればいいのかも．
• 一応2・3次モードぐらいまで確認しておきたいので，正負のモードが出ることを考慮して6次まで要求している．
• デフォルトのCOEF_DIM_ESPACE:4から6に上げてみたところ，2次モードまではWarnが出なくなった（誤差(norme d'erreur)が0.000001以上あると出る）が，計算時間も1.5倍ぐらいに増えた．ちなみに，残りの4モードはWarnが出ているが，逆向きのモードなので無視してもいいと思う．とはいえ，COEF_DIM_ESPACE:4の時でさえ，誤差は有効数字6桁目に影響するという程度で十分小さいと思うけど．

QZ

幾何剛性マトリックス

載荷方向と座屈荷重の正負

• 片持ち梁のオイラー座屈を色々な組み合わせて計算してみたところ，軸の向きとも関係なく，荷重の向きを引張で与えた時，座屈荷重が正で出てくる．
• トラヘアの座屈を載荷方向を変えて計算したところ，載荷方向を変えると座屈荷重の正負が逆に出てくる（つまり，軸の向きに依存する？）．また，トラヘアの場合，片方に載荷すると，絶対値がほぼ同じで正負が逆のモードが出力される．

CalculiXとかMarcとかのこと

• CalculiXは固有値解析(というかほとんどの行列計算)をするのにSPOOLESというライブラリを使っているらしいので，SPOOLESを調べればどういう方法で計算しているのかが分かるかも．また，仲山さんの調査結果を見る限り，CalculiX(?法)とMarc(Lanczos法とInverse Power法)とで結果が違うようだが，CalculiXのマニュアルには"Lanczos vectors"という単語が出てくることから，一応CalculiXの解法はLanczos法の仲間なのだと思う．ただ，CalculiXの板要素S6は立体要素C3D15に置き換えたりしている？ようなので(変なの！)，単にCalculiXとMarcで要素が違うだけなのかもしれないけど．

オンサイト木橋の鋼板をSalome上で自動生成するPythonスクリプト

• onsi_kou.py
• 溝形のフランジ付きモデル．桁高・鉛直材幅・斜材幅・橋長・フランジ幅などを与えると，自動で形状が決まる．1/4解析用に作ったので，半解析や全体解析をしたい場合は，TranslationやMirrorを使う．
• 近い将来(少なくとも2020年頃までに)，Salome-MecaもPython3系に移行するだろうから，一部書き換えないと使えなくなると思う．
• 寸法(弦材幅や斜材の角度)は，Brepファイルにエクスポート -> FreeCAD -> Part -> Measure で確認．いくつかの図面と見比べると，合っていそうだ．
• onsiz_org.py
• 座屈解析用に鋼板，木材，支点の線を作る．鋼板はonsi_kou.pyと全く同じだが，木材はアンカープレートの箇所のみ一体化させるモデル(木材利用2017の「全体的非一体化モデル」．ちょっと作るのが面倒臭い)に対応．
• 次のように要素分割すれば，飛び出しが生じにくいということが分かったため，このモデルで座屈解析をすることにしようと思う．(ちなみに以下に記したLengthはこれ以上小さいと，モデルによっては計算できなくなるという結構ギリギリの値である．座屈荷重の収束性を確認したわけではないが，オイラー座屈・トラヘア座屈の確認をした時，鋼板を2次要素でこのくらい細かく切ると座屈荷重も3桁程度収束する．)
• 鋼板 : NETGEN 1D-2D-3D -> Max.size : 20 , Min.size : 10
• 木材(SubMesh?) : NETGEN 1D-2D-3D -> Max.size : 60 , Min.size : 30
• 境界面(SubMesh?) : NETGEN 1D-2D -> Max.size : 10 , Min.size : 0.1
• Convert to Quadratic
• onsiz.py
• onsiz_org.pyの改良版．要素の飛び出しが鋼板と接しているアスペクト比が極端に大きい木材の要素で生じているため，一体化箇所の段差を鋼板側に設けることによって，要素の飛び出しを回避することを狙ったもの．
• (やり直す時間的余裕がなかったので)今回の計算には使っていない．使ったのは，onsiz_org.pyの方．今後，似たようなことをする人には，onsiz.pyの方法をおすすめしておく．
• onsitz_org.py
• onsiz_org.pyに端部と支間中央の対傾構をくっつけたもの．対傾構は，鋼板と上部の角材と一体化するようにしている．また，1/4解析では実物通りにモデル化することができないので，支間中央の対傾構の位置を対称面に移動させたり，実物よりも危険側になるように寸法を少し変更したりしている．
• onsitz.py
• onsitz_org.pyでは，どうでもいいところ(木材の隙間と対傾構の間)の要素分割が細かすぎて発表までに計算が終わらないので，問題の箇所をフィレットして要素分割を大幅に減らした．また，要素の飛び出しが，鋼板と接している木材でアスペクト比が極端に大きい要素で生じているため，一体化箇所の段差を鋼板側に設けることによって，要素の飛び出しを回避することを狙ったもの．

• 海老さんと川村さんによると，4面体に要素分割をするとき，Max Length と Min Length を指定してやると，メッシュ関係のエラーが出なくなるらしい．

ねじ山を作りたくない

軸方向に色々な荷重が作用した時の柱の座屈

片持ち梁の先端に集中荷重が作用した場合（メモ）

• $P_{cr}=\frac{\pi^{2} EI}{4\ell^2}\approx \frac{2.467EI}{\ell^2}$

片持ち梁に等分布荷重が作用した場合（メモ）

Bessel関数を解いたら(p.97)

• $(q\ell)_{cr}=\frac{7.837EI}{\ell^2}$

エネルギー法1(p.99)

• $(q\ell)_{cr}=\frac{7.89EI}{\ell^2}$

エネルギー法2(p.101)

• $(q\ell)_{cr}=\frac{7.84EI}{\ell^2}$

片持ち梁に等分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合

• $P_{cr}=\frac{(6\pi^5+(3\pi^5-12\pi^3)k)EI}{((-1536+432\pi+8\pi^3)k^2+(24\pi^3-96\pi)k+24\pi^3)\ell^2}$
• 試しに，k=0を代入すると，先端に集中荷重が載った場合の座屈荷重と一致する．

• ちなみに，等分布荷重が引張方向に作用していた場合，kの値を負にすればいいだけなので（一応，計算して確かめてみたが），
• $P_{cr}=\frac{(6\pi^5-(3\pi^5-12\pi^3)k)EI}{((-1536+432\pi+8\pi^3)k^2-(24\pi^3-96\pi)k+24\pi^3)\ell^2}$

• 計算が間違っている可能性も高いので，そのうちにFEMで確認してみよう．

片持ち梁に三角形分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合

• $P_{cr}=\frac{(30\pi^6+(5\pi^6-30\pi^4)k)EI}{((-4080+360\pi^2+6\pi^4)k^2+(40\pi^4-240\pi^2)k+120\pi^4)\ell^2}$

片持ち梁に三角形分布荷重（引張方向）と先端に集中荷重（圧縮方向）が作用した場合

• $P_{cr}=-\frac{\left(\left(5{{\pi}^{6}}-15{{\pi}^{4}}\right)k-15{{\pi}^{6}}\right)EI}{\left(\left(-7560+720{{\pi}^{2}}+8{{\pi}^{4}}\right){{k}^{2}}+\left(120{{\pi}^{2}}-40{{\pi}^{4}}\right)k+60{{\pi}^{4}}\right){{l}^{2}}}$

片持ち梁に台形分布荷重（圧縮方向）と先端に集中荷重（圧縮方向）が作用した場合

• $P_{cr}=\frac{\left(30{{\pi}^{6}}+\left(5{{\pi}^{6}}-30{{\pi}^{4}}\right)k_{1}+\left(10{{\pi}^{6}}-30{{\pi}^{4}}\right)k_{2}\right)EI}{(\left(-15120+1440{{\pi}^{2}}+16{{\pi}^{4}}\right){{k_{2}}^{2}}+\left(-240{{\pi}^{2}}+80{{\pi}^{4}}+\left(18{{\pi}^{4}}+360{{\pi}^{2}}-7680\pi+19200\right)k_{1}\right)k_{2}+\left(-4080+360{{\pi}^{2}}+6{{\pi}^{4}}\right){{k_{1}}^{2}}+\left(40{{\pi}^{4}}-240{{\pi}^{2}}\right)k_{1}+120{{\pi}^{4}}){{l}^{2}}}$
• maximaに貼り付ける用

P1(k1,k2,l1):=((30*%pi^6+(5*%pi^6-30*%pi^4)*k1+(10*%pi^6-30*%pi^4)*k2)*EI)/(((-15120+1440*%pi^2+16*%pi^4)*k2^2+(-240*%pi^2+80*%pi^4+(18*%pi^4+360*%pi^2-7680*%pi+19200)*k1)*k2+(-4080+360*%pi^2+6*%pi^4)*k1^2+(40*%pi^4-240*%pi^2)*k1+120*%pi^4)*l1^2)

• $k_{2}$に0あるいは$k_{1}$に0代入したら三角形分布と同じになる．
• $k_{1}=k_{2}$にしたら等分布と同じになる．

片持ち梁に台形分布荷重（引張方向）と先端に集中荷重（圧縮方向）が作用した場合

• ↑の式をk1とk2をマイナスにしたら同じになる.
• $P_{cr}=\frac{\left(-30{{\pi}^{6}}-\left(5{{\pi}^{6}}-30{{\pi}^{4}}\right)k1-\left(10{{\pi}^{6}}-30{{\pi}^{4}}\right)k2\right)EI)}{{(\left(-15120+1440{{\pi}^{2}}+16{{\pi}^{4}}\right){{k2}^{2}}+\left(240{{\pi}^{2}}-80{{\pi}^{4}}+\left(18{{\pi}^{4}}+360{{\pi}^{2}}-7680\pi+19200\right)k1\right)k2+\left(-4080+360{{\pi}^{2}}+6{{\pi}^{4}}\right){k1}^{2}}+\left(240{{\pi}^{2}}-40{{\pi}^{4}}\right)k1+120{{\pi}^{4}}){{l}^{2}})}$
• つまり，maximaでは，P1を貼り付けたあとに次のようにすればいい．

P2(k3,k4,l2):=P1(-k3,-k4,l2)

両端固定の問題に応用したい場合

• まず，上の片持ちのそれぞれの式（分布荷重が圧縮方向と引張方向）にkの値を代入する．
• そして，Pcrをそれぞれ$Pcr_{1},Pcr_{2}$，スパン$\ell$をそれぞれ$\ell_{1}, \ell_{2}$として，
• 連立方程式：$Pcr_{1}=Pcr_{2}, \ell_{1}+\ell_{2}=0.5$
• ここでは，具体的な値ではなく，便宜上$\ell_{1}+\ell_{2}=0.5$としたが，本来は両端固定のスパンの半分の値を入れるべき．
• より，$\ell_{1}, \ell_{2}$を求め，元の式$Pcr_{1}$か$Pcr_{2}$に代入すれば，座屈荷重Pcrが求められる．

等分布載荷

• かなり長いので，書きたくない．maximaに貼る用のみ．

sol(k1,k2):=solve(P1(k1,k1,0.5-l2)=P2(k2,k2,l2),l2)

• 0.5は仮で，本当は$\frac{\ell}{2}$をいれる．別に$\ell_{2}$でなく，$\ell_{1}$についても解いてもいい．
• 等分布荷重が載っていない場合と同じになるかの確認のためk1=k2=0を代入すると，$\ell_{2}=0.250$になる．また，それを$Pcr_{2}$に代入するとPcr=4$\frac{\pi^2 EI}{\ell^2}$とも一致する．

bfloat(P2(k2,k2,sol(k1,k2)))

• ここでは，kに具体的な数値を入れる．
• kの計算方法については考え中．多分，曲げモーメントによる軸力をzで一回微分して，zに何かを代入して，さらにそれに何かをかけるのだろうけど，k自体がl2とP2の入れ子になっていたら非常に厄介．

台形分布載荷

• かなり長いので，書きたくない．maximaに貼る用のみ．

sol(k1,k2,k3,k4):=solve(P1(k1,k2,0.5-l2)=P2(k3,k4,l2),l2)

• 0.5は仮で，本当は$\frac{\ell}{2}$をいれる．別に$\ell_{2}$でなく，$\ell_{1}$についても解いてもいい．

bfloat(P2(k3,k4,sol(k1,k2,k3,k4)))

• ここでは，k1, k2, k3, k4に具体的な値を入れる．
• k1,k2,k3,k4の計算方法については考え中．

直交異方性材料の直応力-直ひずみ関係

• 式(1)について、Maximaで求めた結果(ihou1.wxm)を(wiki用に書き換えるのが面倒臭いので)画像のまま貼り付けると、

• 確かに複雑……。

木材利用2017メモ

• アンカープレートの部分だけ木材と鋼板を一体化させる。
• アンカープレートの境界面でパーティションし、木材はLength:50,鋼材はLength:10ぐらいの四面体要素に切り、2次要素にする。
• これだと、鋼板の隣の非一体化部分が木材にすり抜けるので、従来のモデリングの方がいかにもそれらしい座屈をする。ということで、従来に近い非一体化方法でもやる。
• 堅固山さんとの違いは、
• 木材を直交異方性にする。
• 等方性と異方性にしたときの座屈荷重の比較をする?
• 3点曲げに加えて積雪荷重での座屈も計算する。
• 特に、支点付近の鉛直材に焦点を当てる。
• つまり、中央以外の非一体化部分もモデリングする。
• 大池や虎毛などの大きいモデルでもやってみる。
• 要素数の問題で部分的に取り出したモデルで計算するかも。

5月初旬(発表概要)までに

• 7mモデルでやる。
• 穴の拡大率は0,5,10,15,20,25,30%の7通りとし、7mモデルの鋼板については堅固山さん作のデータを使用。
• 何もなしモデルと補剛材+対傾構ありモデルの2通りでやってみる。
• 載荷条件は、3点曲げ(線載荷でいいかな?)と積雪荷重(積雪1mあたり$3.5kN/m^2$)の2通りでやる。

これだけで28通り。ここまでの計算自体はjobを実行すれば2〜3日で終わるけど、ファイルのサイズが結構な量で120GBもいってしまうので、保存場所に困る。

• さらに木材を等方性と異方性にしたものとで比較する?ことになったので、2倍して56通り?(済)
• 中央と支点部の鋼板と木材だけ非一体化にしたものでもやってみるので、さらに2倍して112通り?(済)

ひとまず、7mスパンについては全通りの計算が終わったが、結果が多すぎて、どうまとめたらいいものやら分からない……。

後藤さんと相談した概要の流れをメモ

• 1.端だけ非一体化・等方性で穴の大きさを変える(堅固山さんのようにする)と、穴の大きさと座屈荷重はほぼ比例関係になる。
• 2.しかし、木材を直交異方性にすると、それよりも座屈荷重が2割も下がってしまう。
• 不思議。それと、要素の跳び出しが等方性の時よりも増えるので、あまりいい気はしない。
• (3.)さらにアンカープレートのところだけ一体化にすると、座屈荷重がもう1割下がるが、これはめり込みが起きるからであり、鋼板中央部に着目する場合、先ほどの数値モデルの方が現実的である。
• しかも、木材と鋼板の境界に応力集中する箇所が増えるからか、要素の跳び出しがさらに増える。これもあまりいい気はしない。

• 4.実際に影響が大きいと予想される積雪荷重が載った場合、3点曲げよりも低次のモードで支点部付近の鉛直材の座屈や支点部の横倒れ?が発生する。しかも、その順番は穴の大きさに依存する。

• 対傾構と補剛材の影響をどこに挟むのがいいのかを訊くのを忘れてた!　座屈荷重が1桁ぐらい高くなるので、これを入れたら7mスパンの場合、かなり安心できる。

6月から8月(発表スライドor論文)までに

• 大池モデル(や虎毛モデル)でも同じことをやる。
• 7mモデルで支点近くだけを取り出した一部モデルと従来の1/4モデルで同じ結果がでることを確認できたら、一部モデルでやるかも。

こんなにやれるかな？

おまけ

• 5mスパン・板厚6mmの図面があったら(k2にはない?)、本当に70kN近くで実験のように座屈(はらみ?)するのかを確かめておきたい(済)。
• 座屈荷重は，部分的非一体化モデルで282kN, 全体的非一体化モデルで252kN.実験値と全然合わないが，これでいいのだろうか……．
• というか，そもそも実験値（tugite）のはらみ（7mに3点曲げで100kNかけたもの）が6mmもあるのに、FEMだとその1/20ぐらいしかない．まずは，座屈解析に入る前に，そのへんの原因をはっきりさせて，はらみが合うようなモデル化を見つけないけないと思う．
• ちなみに，CalculiXで座屈の*STEP〜*END STEPの前にプレストレスの*STEP〜*END STEPをいれると，プレストレスを入れた状態での座屈解析ができるが，プレストレスがあってもなくても座屈荷重はほぼ変わらない（相対誤差1%未満）．
• どうやら非線形座屈解析の出番はないかも。計算するにしても、ひとつにつき2〜3日は掛かりそうなので、厳選してやる必要があると思う。ただ、遊びでやってみたい。
• 幾何学非線形で計算すると，座屈荷重に至る前にまた別のあるところに応力が集中しすぎて，そこで飛び出しが発生する．困ったものだ．ヤング率と荷重のオーダーを一緒に下げてみるとか？

• 計算に使った積雪1mにつき$3.5kN/m^2$だが、昔の道路橋示方書(5版)によると、

印刷用の折り紙円筒

• システム-->システム管理-->プリンター
• HL5270DNを右クリックして，プロパティ
• ジョブオプション-->「用紙に合わせる」にチェック-->適用

• prjabara.f90
• 蛇腹折り。例えば、周方向6分割・高さ方向8分割の例。今までオープンキャンパスで配っていたものよりもきれい。
近藤夏実さんが折り畳まることを確認済み．
• prkrasen.f90
• 回転螺旋折り。例えば、周方向5分割・高さ方向10分割・$\beta$=30度の例。
近藤夏実さんに折って貰ったところ，あまりにも糊付けが難しそうなため，ちゃんと畳めるかどうか確認ができていない．とても悪い例であり，オープンキャンパスで配るべきでない！配るならパターンを変えるべし．
• prhrasen.f90
• 反転螺旋折り。例えば、周方向5分割・高さ方向6分割・$\beta$=30度の例。
近藤夏実さんが折り畳まることを確認済み．こちらは糊付けが楽そうだが，見た目はどうかな？
• 螺旋折り系は、紙に目一杯描くようにしていないため、パターンによっては紙からはみ出る(要微調整)。$\beta$が折り畳めない角度を入力していた場合、.texのヘッダーにエラーメッセージを出すようにしている。
• prdaiya.f90
• ダイヤカット円筒。例えば、周方向6分割・高さ方向8分割の例。
川村さんに折って貰ったところ，ちゃんと完成する．（これは折り畳めなくていい）

支間中央のみを取り出して線形座屈解析

• 拘束は
• (z)橋軸方向対称面
• (x)幅員方向対称面
• (y)切断面の中立軸
• 載荷は
• (y)輪荷重P/4
• (z)$\frac{曲げモーメント_{切断面}}{角材同士の中心間隔}$を軸方向に上下の角材(+鋼板)で逆向きに作用させる。
• ここで、角材と鋼板の剛性比($\frac{E_{鋼}I_{鋼}}{E_{鋼}I_{鋼}+E_{木}I_{木}}$or$\frac{E_{木}I_{木}}{E_{鋼}I_{鋼}+E_{木}I_{木}}$)をちゃんと考慮する。そうしないと、鋼材の負担が少なくなるため、座屈荷重が高くなってしまう。

• 孔空き鋼板の3点曲げ(載荷は輪荷重)とした。
• モデルはk2に上がっていた堅固山さんのbrepファイルをそのままお借りしました。ありがたや。
• 堅固山さんが計算したものは、$P_{cr}$=270kNなので、1/4モデル(等方性)とほぼ一緒。

• ちなみに、応力集中する要素の飛び出しなどのノイズと思われるものは除いている。
• モデルの体積比は$\frac{1/4モデル}{一部モデル}=5$なので、要素数も5倍にならないといけないと思うかもしれないが、それだと計算が大変なので、1/4モデルは全体的にそこそこ粗めに要素分割して、鋼板の非一体化部分だけ(sub meshで)一部モデルと同じぐらい細かく要素分割している。
• 変位図(変位の倍率は見やすいように調整)
• そんなことをしていいのだか分からないけど、1/4モデルの1次モードと2次モードは見にくいので、幅方向だけ変位を逆にしている。

• 相対誤差が大体2~3%ぐらいなので、なかなか良い精度だと思う。
• 分布荷重が載った場合についても、いずれやりたい。

後藤ちゃちゃ(17/3/19)

• ただの近藤メモ:堅固山さんに聞いてみたところ、純粋な鉛直材の座屈を確認したいときは、50モードぐらい出力するようにしているらしい。ちなみに、なぜここで「純粋な鉛直材の座屈」と書いたかというと、それよりも低い荷重で、非一体化部分の鋼板が開きながら鉛直材も折れ曲がるという座屈をするから。
• 後藤ちゃちゃ(17/3/23)虎毛や大池のモデルは、鋼板上下が折れ曲がっているチャネル材です。フランジ部は、ウェブ厚を含めて75mm, ウェブ厚を引いた張り出し部だけなら75-9=66mm（図面アップしてましたっけ）。

理工学デザイン

システムデザイン工学専攻	１２月２２日（木）１２時５０分〜１４時２０分
○会場：ＶＢＬ棟２階大セミナー室
○発表者は開始５分前までに準備を終えて会場に集合のこと。
○発表者は３名（教員２名、学生１名）への説明が終わった後も、
他の学生のポスターを見学したり、担当でない先生方へも対応するなりして、
終了時間まで会場に残ってください。
○当日、発表も質疑もない学生は出席しなくて構いません。

○ポスター発表当日の午前１１時からポスターの貼り付けなど準備できます。
○ポスター用ボードにポスター発表番号が記載されていますので、
指定の場所にポスターを貼り付けてください。
○各自、ピンかセロテープを用意して、ポスター用ボードに貼り付けてください。
○終了後は、忘れずにポスターを持ち帰りください。

○Ａ１サイズ（縦８４１ｍｍ×横５９４ｍｍ）縦置きでポスターを準備してください。
○ポスターの上段に、タイトル、所属（○○専攻○○コース）、
発表者氏名を明記してください。
○ポスターの作成は自由（手作りでも結構です）に行ってください。
情報統括センターを利用してポスターを作成する人は、
予めセンターに予約（前日17:00まで）してから行ってください。
　　　なお、ポスター印刷する際は、利用申請書を必ず持参してください。
○研究テーマや方向性、方法をポスターにまとめ、発表する。
指定された教員２名に説明する以外に、他専攻の学生に対して説明を行う。
○各学生は、自身の属する専攻以外の学生のポスター発表を訪れ、
質疑応答し、その学生の研究内容を十分に理解した上で研究内容に関してレポートを作成し、
1か月以内に理工（大学院）担当 に提出する。
特に、複数の質問事項を明確に記載し、質疑応答を記録する。
その質疑応答から理解した事をレポートにまとめる。

• 後藤ちゃちゃ(16/11/18)
各種の非線形問題（ねじり座屈とか接触とか）一通り解けているようだし （CalculiXも一通り解けそうですね）、線形でも回転螺旋折りとかも面白いし、 解析が困難な各種の非線形問題を工夫しながら解いてますが、 オープンソースのツールだけで、割と何でも解けてます的なオムニバス的な 話にしてもいいかもしれません。まだM1なんで色々とおもしろそうなところを やってます。これからやっていきたい本命はこれなんですが、というのを最後に 入れてもいいですが。

PC鋼棒の引張+トルク+摩擦

• toruku.f90
• PC鋼棒を引張った後、トルクを導入し、滑り出しをみる。
• PC鋼棒(梁要素)の先端に円盤(立体要素)をくっつけて、もうひとつの円盤(立体要素)と接触するようにした。しかし、円盤の大きさや摩擦係数はどうしたら良いものか。まだまだ、修正が必要。
• 菊地さんのと境界条件や荷重条件がやや異なるので、いずれそちらに合わせようかな。
• tth.f90
• gnuplot用にdatをまとめる。

CLT関連で使えそうな板の文献をメモ

• 土木学会『構造力学公式集』@研究室
• 相対2辺が単純支持された等方性長方形薄板：p.328
• 部分分布荷重：p.330
• 残りの相対2辺が弾性支持桁：p.331 .ここで、$\gamma=0$とすると、残りの相対2辺が自由辺の条件になる。
• 直交異方性板：p.367 .少しだけ長方形板の記述あり。ただし、これを見ただけではなかなか使えるようにはならないと思う。既に勉強した人の為のリファレンス的なもの？
• ティモシェンコ『板とシェルの理論(上・下)』 @図書館
• 直交異方性板は下巻に書かれている。これも具体例はあまり無い。
• 渡辺昇『橋梁工学』 @研究室
• 等方性板：p.300
• 直交異方性板：p.324（そんなに詳しくない…）
• アムバルツミャン『異方弾性板の理論』 @図書館
• 異方性の薄板、厚板、積層板について座屈、有限変位、振動といった具合に何でも書かれているが、結構読むのが大変そう。
• J.R. Vinson & R.L.Sierakowski『複合材料の構造力学』@図書館
• 異方性の積層板について：pp.42-101
• さらに、7章に「強度と破損の理論」が載っている。これは面白そう。
• 宮入裕夫『複合材料入門　ー基礎と応用ー』@図書館
• 異方性の積層板についてはあまり書かれていないのだが，何と！ハニカムパネルの曲げ剛性・せん断剛性・ねじり剛性について書かれている（p.180~）．
• ということで，ハニカムパネルについて読んでみると，面内引張・面外曲げ・ねじりについては表面材（上下の蓋）のみで抵抗し，面外せん断については心材（ハニカムコア）のみで抵抗する，という大胆な仮定を用いていた．工学的！
• ハニカムパネルについては，図書館にはないが，宮入裕夫『サンドイッチ構造』あたりにより詳しく書かれているのでは？

板の座屈

• メモ：4辺単純支持の座屈荷重
$N=k\frac{D\pi^{2}}{b^{2}} .$ ここで、 $k=(m\frac{b}{a}+\frac{1}{m}\frac{a}{b})^{2} , D=\frac{Et^3}{12(1-\nu^{2})} .$ t:板厚　a:長さ　b:幅　m:長軸方向の凹凸の数

と、色んな資料に書かれているが、Dは単位幅あたりの曲げ剛性なので、(次元的にも)上の式に幅bを掛けた$N=k\frac{D\pi^{2}}{b}$が座屈荷重なるのではないだろうか。

• itaz.f90
• 4辺単純支持を板要素(S8)で計算。幅厚比b/t=100で試してみたところ、理論値(修正した方)よりも1%ぐらい座屈荷重が小さくなる。
• hitaz.f90
• itaz.f90の結果で非線形座屈計算。4辺単純支持した板は、それほど変形が大きくないので、分かりづらい。

プレストレス木箱桁橋

孔無し一体化モデルの1/4解析

• onsiyon.f90
• プレストレス木箱桁橋プログラムの練習として一番単純なモデルを作ってみる。プレストレスも入れられるようにした。プレではないけど。
• 3点曲げか、4点曲げか選べるように改良(『木材の特徴を活かす木橋』のp.19と比べたかったから)。載荷板に載荷するように改良。変位計の位置で変位を出力するように改良。
• 1/4解析にしたのは、計算が大変な接触解析とか局部座屈とかへの移行を想定したため。

孔無し非一体化モデルの1/4解析

• onsimyon.f90
• プレストレスのかけ方は、鋼板の上下の木材と隣り合わせている部分に全面載荷。
• 載荷板に載荷するように改良。変位計の位置で変位を出力するように改良。なるべく均等に載荷するように改良。
• onsimayon.f90
• アンカープレートのみにプレストレスをかけてみたが、プレストレス時の計算が収束しない。
• せっかくなので、スパン7m($E_{zz}$=7.694GPa)の孔無し一体化モデルと孔無し非一体化モデルの断面のずれ具合を比較してみた。

• とりあえず摩擦係数は中央値をとって、鋼材と木材:0.40 , 木材と木材:0.40 にした。
• 軸方向は128分割で、上村さんの修論本体のグラフからするとこれで十分と判断。その代わり、断面の分割を細かくしてみる(上図)と、断面の変位差が大きくなった。ちなみに鋼板の厚さは4分割。
• 変位図(載荷方向)は、10倍に拡大している。
• 実験値は上村さんの修論本体の表の実験2(k2にあるuemurahontai2.pdfのp.102)から持ってきた。一体化モデルでも、十分実験値に近い結果になった。
• 濡れ具合の他に、↓摩擦面の仕上げ方によっても、摩擦係数は大きく変わるらしい。現地製材の木材と、橋梁に使うような大型の鋼板であれば、摩擦面がかなり粗いと思われる。
• ttps://www.jstage.jst.go.jp/article/aijs/76/666/76_666_1469/_pdf

緊張力の低下

• onsirelax.f90
• onsimyon.f90の計算後、荷重はそのまま一定に保ち、緊張力だけを低下させる。
• 載荷板に載荷するように改良。変位計で変位を読み取るように改良。なるべく均等に載荷するように改良。
• とりあえず、スパン6m($E_{zz}$=6.921GPa)のモデルに荷重70kN(二等林道橋の後輪荷重55kN*(1+木橋の衝撃荷重0.25))を載荷し、緊張力を70kNから40kN(根拠なし)まで減少させてみた。摩擦係数も根拠なし。
• スパン7mモデル($E_{zz}$=7.694GPa)にして、計算し直そう。

• 上図は緊張力が53kNの時の載荷方向の変異(1倍)図。落下防止材やら対傾構やらで、これほどは滑らないような気がする。剛性は低下すると思うけど。
• 破壊試験(100kN)では、鉛直変位の最大値が15mmぐらいなので、それを目安に考えると、1〜2割ぐらい緊張力が落ちていたら危ない?
• グラフには鉛直変位20mm以上は載せていない。データが途中で途切れているのは、そこでいきなり滑って、計算が止まったから。
• ところで，色々とひっくり返すようなことをいうと，緊張力はそもそも乾燥収縮やクリープやらで40%までに低下したときに丁度良くなるように設計されていたはずなので，緊張力が28kNの時に完全に鋼板と木材が滑っているのはいいのだろうか？まあ，木材同士はずれてないから，元のプレストレス木床版としてはそれでいいのかもしれないけど……．

孔無し一部非一体化モデルの1/4解析(局部座屈)

• onsikyo.f90
• プレストレスを鋼板の上下木材部分に全面載荷。もちろん、プレストレス力を0にもできる。
• onsikyoa.f90
• プレストレスをアンカープレートのみに載荷。こちらの方がやや実物に近い。
• これらをスパン7m($E_{zz}$=7.694GPa)、鋼板の板厚を9mmとして計算し、pd2onsi.f90でまとめると、

• 下の図は500kNの時の幅員方向変位の1倍。当たり前だが、そんなに派手に座屈してくれない。
プレストレス無し 全面プレストレス アンカープレートのみプレストレス
• 斉藤輝さんの線形座屈解析の結果よりも随分と結果が危険側にずれてしまった。と思っていが、どうやら輝さんのは、1/4解析なのに荷重1MNを与えている(k2にある"onsi4z2.f90")ようだ。これを荷重も0.25MNに直したら、座屈荷重が4倍になって、なかなかいい具合に非線形座屈の値と一致した(Pcr=339.1kNとかになる!)。
• ただし、鋼板の板厚方向の要素分割は、今は適当に4分割だが、輝さんに合わせるなら、8分割ぐらいに増やしたほうがいいかも。
• 理工学デザインの発表用に、線形座屈と非線形座屈を同じ要素分割でなるべく細かくして計算したものがあるが、全部書き換えるのは面倒なので、上の結果には反映させていない。ほんの少し変わるけど、大体こんな感じ。

オイラー座屈から3点曲げ載荷時の座屈荷重を推定

• onsioira_3ten.f90
• 堅固山さんが(孔空き鋼板の)オイラー座屈を計算しているので、オイラー荷重から3点曲げの荷重を逆算してみた。圧縮力を分布荷重ではなく集中荷重に置き換えているところや，断面内でのせん断力の分布を考慮していないところなどが特に怪しいが、こんなものでもオーダー的にそれなりの予測ができそう。

• 「異部分」は，木材を直交異方性材料とした，部分的非一体化モデルのこと．
• 「等部分」は，木材を等方性材料とした，部分的非一体化モデルのこと．
• 鋼板上部の座屈について，ここでは(1.直交異方性を考慮していないので座屈荷重は上がる↑)が，それと同時に(2.斜材と鉛直材による弾性支持を無視しているので座屈荷重が下がる↓)ため，うまいことたまたまFEMと近い値になったのだと思う．
• ということで，FEMで中央の鉛直材と斜材をなくして，木材を等方性材料にしたり異方性材料にしたり，はたまた異方性で，ヤング率（$E_{xx},E_{yy}$だけ）を高くしたり，せん断弾性係数（$\tau_{xy},\tau_{zx}$だけ）を高くしたりしてみると，

• となることから，座屈荷重は，異方性の影響が約2割↓，中央の鉛直材の影響が約2割↑なのだろう．
• 特に，異方性の中でも，木材の軸直角方向のヤング率の影響よりも，せん断弾性係数の影響が大きいことが分かる．つまり，アンカープレートらへんが斜め向きになるということ．
• ちなみに，ヤング率を変える（軸直交方向を硬くしても柔らかくしても）と「飛び出し」が生じるが，せん断弾性係数だけを変えても「飛び出し」は生じない．
• ただし，桁高800mmになると中央の鉛直材と斜材の補剛効果が下がるので，この推定方法よりも座屈荷重が低くなってしまう．ただ，桁高500mmモデルの中央の鉛直・斜材の効果が約2割であることが分かったので，計算結果よりも$2-2\times (\frac{桁高500mmの鉛直材の長さ}{桁高800mmの鉛直材の長さ})^3$割低下すると予測できると思う．

• 3点曲げ載荷における鉛直材の座屈は$\frac{\ell}{4}$あたりで生じる（しかも，鉛直材だけでなくて板全体が曲がる）ため，この方法では推定できない気がする．

オイラー座屈から等分布載荷時の座屈する積雪深を推定

• onsioira_yuki.f90
• onsioira_3ten.f90の精度が割といいので，等分布載荷でも同じことをやってみた．

• 鉛直材の長さの決め方が難しく，それによって結果が2倍ぐらい変わってしまう．とはいえ，座屈荷重があまりに大きすぎて現実的な値ではないため，鉛直材の座屈について議論する必要がないのではないかと思う．

• 対傾構をつけず，補剛材をつけたらどうなるかやってみたい．波打つ様な座屈でないなら，それなりにうまく行く気がする．
• Trahairでも$S(\frac{\ell}{2})$や$\tau=\frac{QS}{bI}$とかを整理したらできるかもしれない。

datファイルのデータ整理

• plonsi.f90
• 一体化モデルの1/4解析用。
• plonsim.f90
• 非一体化モデルの1/4解析用。
• plonsirelax.f90
• 緊張力の抜けの1/4解析用。

弾性床上の座屈

• yuka.f90
• 特に初期不整を与えなくても、断面が変形して弾性床と接触することで勝手に偏心していき、そのうち座屈する。
• ただし、"*SURFACE BEHAVIOR"の値をかなり上手く調整しないと、計算できない。
• やってみたのは単純梁と片持ち梁だが、どちらも床と離れて座屈した。この場合の座屈荷重(見た目で判断)は、普通の座屈とほとんど同じになる。
• 仮に床方向に座屈させたとしても、鋼材と木材(繊維直角方向)だとヤング率比が1000倍近くあるので、鋼材を座屈させる場合、弾性床の恩恵をほとんど受けないかも。

鋼材と木材の接触解析

• masatu.f90
• JIS K7125を真似る。
• 鋼材と木材の摩擦係数は0.2(湿)~0.6(乾)の中央値をとって0.4とした。
• ちなみに木材(上に載っている方)の繊維方向は滑らす方向と垂直(横向き)にしてある。
• "*SURFACE BEHAVIOR"が謎。
• Abaqusの"*SURFACE BEHAVIOR"の説明
• "*SURFACE BEHAVIOER"をあまりいじらない方がいいと言っている人もいるようだが、寸法とかヤング率によって変化させないと収束しないことが多いように感じる。
• 要は収束の「する」、「しない」に関わってくる。
• "*FRICTION"の右のstick range slopeの値が謎。
• マニュアルとfigure 119によると"stick slope"は、力/体積の次元を持っているとか何とかで……、ペナルティー法のパラメータだそうだ。
• ネットを漁っていると、食い込みのない理想的な値は∞だが、大きすぎると収束しなくなるので、適度に小さくする必要があるらしい。ひとまず、海老さんにならって、軟らかい方のヤング率ぐらいにしておくと丁度いいかな?

• pd.f90でまとめると、

• 縦軸が水平力、横軸が水平変位。変位は載荷面(横)の平均であり、試験材がせん断変形しているため、滑っていなくても緩やかな変位があるのだと思う。

CLT

梁

梁(5層7プライ/7層7プライ)の3点曲げ

• clt.f90
• 木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似る。層数および層ごとの向きを設定できるようにした。
• ひき板の強軸方向ヤング率を試しに6GPaにして計算してみると、5層7プライと7層7プライのどちらの場合も実験結果に近くなる。

梁(5層7プライ/7層7プライ)の３点曲げを繰り返す(全体が弾塑性体)

• cltso.f90
• 木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似る。ただし、時間がかかりそうなので繰り返し回数は少ししか試していない。
• CalculiXの"*PLASTIC"はミーゼス応力での降伏値。木材は圧縮と引張で降伏の仕方が異なるので、千田さんのように圧縮部のみを弾塑性体にする方法もあるようだ。
• "*PLASTIC"の行に続けて",HARDENING=〜"を入れると、硬化則(isotropic,kinematic,combined)を設定できる。何も入れないと、等方硬化則。
• 応力-ひずみ曲線(多直線近似)はかなり適当に入れてあるが、一応、ある程度の大きさの荷重を超えた後に除荷すると、永久ひずみが残ってくれる。
• "*PLASTIC"ではなく"*DEFORMATION PLASTICITY"を使うのもありかも。
• フィンガージョイントの降伏点が分かれば、そこだけ弾塑性材料にして破壊できるが、床版の場合、そこよりも先にボルトが降伏しそう。

h=240mm(7層8プライ)のときの梁理論との比較

• clt3ten.f90
• clt5ten.f90
• pclt35ten.f90
• 特に意味はないけど、大黒屋さんの図6-7をCLTでやってみた。
• 桁高は重ね方の都合上、変えられないので、全方向の寸法を大黒屋さんの1.2倍にした。
• つまり、b=24mm,h=240mm,$\ell$=240〜4800mm(240mm刻みで変化).ただし、$\ell$が全体の半分になっていることに注意。
• 大黒屋さんに合わせて、分割数はnx=5,ny=24,nz=100.(大黒屋さんに合わせると、ny=20になるが、8層もののため、nyは8の倍数でしか分割できない)

• 結果↑は大黒屋さんのものに似ている。

b=24mm,h=240mm(7層8プライ)のときの1層と8層のFEMでの比較

• 大黒屋さんの真似をしたついでに、FEMの相対誤差($\frac{v_{1層}-v_{8層}}{v_{8層}}\times100$)を一応まとめておく。

• 曲げが支配的になるほど、両者の値が近くなる。曲げ剛性から1層モデルのヤング率を求めていることから、予想がつくけど。
• 弱軸の挙動もちょっと気になる。3点曲げと5点曲げで違う形状の曲線を描いているあたりが。

床版

供試体2(7層8プライ/接続金物無し)の静的載荷

• clt1.f90
• 1層モデル。海老名さんの指摘により寸法や境界条件を直した。変位計の位置で変位も出力するようにもした。弾塑性への移行等も想定して、半解析にした。変位計の位置はExcelのデータと図面で違う？
• $\overline{E_z}=3.5$GPa, $\overline{E_x}=1.7$GPa, $E_y=0.2$GPa,
• $\nu_{xy}=\nu_{xz}=0.4\cdot\frac{\overline{E_{x}}}{\overline{E_{z}}}=0.4\cdot\frac{1.7}{3.5}\approx 0.194$, $\nu_{yx}=\nu_{yz}=0.4\cdot\frac{E_{y}}{\overline{E_{z}}}=0.4\cdot\frac{0.2}{3.5}\approx 0.0229$, $\nu_{zx}=\nu_{zy}=0.4$(支間方向が支配的だと想定してポアソン比0.4をここにもってくる),
• $G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\approx 0.33$GPa.
• ちなみに、ヤング率の平均化は、$\frac{\sum_{i=1}^{8} E_{i}\cdot\int_{A_{i}} y^2 dA_{i}}{\sum_{i=1}^{8}\int_{A_{i}} y^2 dA_{i}}$というようにして計算すると、ぴったりこの値になる。
• clt8.f90
• 8層モデル。海老名さんの指摘により寸法や境界条件を直した。変位計の位置で変位も出力するようにもした。弾塑性への移行等も想定して、半解析にした。変位計の位置はExcelのデータと図面で違う？
• $E_z=5$GPa, $E_x=E_y=0.2$GPa.
• $\nu_{xy}=\nu_{xz}=\nu_{yx}=\nu_{yz}=\frac{0.4}{25}=0.016$, $\nu_{zx}=\nu_{zy}=0.4$,
• $G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\approx 0.33$GPa.
• cltiso.f90
• 等方性材料。実験や上の2つとは全然合わないけど、構造力学公式集に載っている板のたわみの式に近くなる。

供試体2(7層8プライ/接続金物無し/ゴム有り)の静的載荷

• clt1g.f90
• clt8g.f90
• ゴムのヤング率は分からないが、作るだけ作っておいた。ちょっと動かしてみた感じだと、それほど大きな影響はないような……。
• ちなみに、*HYPERELASTICでゴムのような超弾性も扱えるが、超弾性について勉強していないのでちゃんと扱える気がせず、試しておりません。

供試体2(7層8プライ/接続金物無し)の弾塑性解析

• clt1so.f90
• clt8so.f90
• 中立軸の上下で別々の弾塑性材料にしただけ。応力ひずみ関係の値は、ちゃんとkatuyou/4shou.pdfを読んでから書きなおそう。半解析。
• とはいえ、本とか色んなページ(例えばこことか)に書かれているスギの機械等級E50(曲げヤング率:3.9〜5.9GPa)の基準強度と比べると、katuyou/4shou.pdfに書かれている基準強度は高すぎる気がする。藤田さん曰く、1000kN載荷でもほとんど塑性化しないというし。
• しかも、応力-ひずみは、この辺のをそっくりそのまま真似ただけなので、ちょっと雑な気がしないでもない(とくにひずみの値)。
• clt1soko.f90
• clt8soko.f90
• 中立軸の上下で別々の弾塑性材料にし、かつ支持部周辺を圧縮降伏になるようにした。ただし、「支持部周辺」をかなりざっくりと決めているため、ある程度の試行調整が必要。
• 藤田さんの計算結果が実験値よりも硬いのは、支持部周辺は圧縮が作用しているのに中立軸より下を全部引張降伏にしたことによるかもしれない思い、とりあえず作ってみた。しかし、まだ計算はしていない。
• plcltso.f90
• gnuplot用に荷重-変位関係をまとめる。

• ところで，入力データに特に書き込んでいないが，デフォルトでMisesの降伏条件を直交異方性材料に拡張したというHillの降伏条件になっているのだろうか？ソースをいじってみたいのだが，どこにあるのかいまいち分からない．CAELinuxでは/opt/CalculiXにあったのだけど．
• 後で後藤さんから訊いた話では，apt-getでCalculiXをインストールすると，ソースファイルをいじって使うことができないらしい．実行ファイルだけインストールされるとかかな？残念．

供試体1(7層8プライ/接続金物有り/ゴム有り)の静的載荷

• clt1jg.f90
• clt8jg.f90
• 鋼材は板要素(S8)として木材と節点共有させたので、若干寸法がずれている。鋼材と木材間の滑りを考慮していないので、実験よりもかなり変位が小さい。

供試体3(7層8プライ/供試体2の半分)の静的載荷

• clt1han.f90
• clt8han.f90
• 図面をちゃんと見ていないので、寸法は正確でないと思うけれど、一応作っておいた。実験の結果も見ていないけど。

$\ell$=2010mm,h=240mm(7層8プライ)のときの1層モデルと8層モデルのFEMでの比較

• clthaba.f90
• H鋼の間隔($\ell$=2010),床版厚(h=240mm(7層8プライ))を一定、橋軸方向だけb=400〜4000mm(400mm刻みで変化)。1/4解析。
• 載荷は、床版中央に500x200mm(の1/4)の大きさを面載荷。
• 境界条件は、面を(ボルト〜H鋼の端)高さ方向と、横方向に固定。
• clthababun.f90
• 床版の設計上で知っておきたいのは中央集中荷重よりも分布荷重時の挙動だと思ったので、clthaba.f90を分布荷重にしてみた。
• pclthaba.f90
• clthaba.f90の相対誤差($\frac{v_{1層}-v_{8層}}{v_{8層}}\times100$)をまとめる。読み取るのは、載荷面の平均値。
• 分布荷重についても同様に中央のたわみを見たが、分布荷重の場合、板の中央よりも支持していない辺の方が変形するので、これでいいのだろうか?
• たわみの収束性の確認をするのが面倒だったので、要素分割は24x100x100(高さx幅x長さ)というように計算できる限界まで細かくした。

• 8層ものを1層のモデルに曲げ剛性から近似した場合、梁では両者の誤差は少ないが、板に近づくにつれ1割ぐらいの誤差が生じることが分かる。
• あと、変位の相対誤差の大小に関わらず、1層モデルのミーゼス応力が8層モデルの約2/3ぐらいになる(単にヤング率の比3.5/5の関係?)ので、降伏についての考慮をするときはちゃんと8層モデルにした方がいいと思う。
• それから、境界条件をCLT床版に近くして計算したが、より一般的な問題を考えるのなら、中立軸でローラー支承にした方がいいかも。
• それなら、立体要素の8層モデルと1層モデル、板要素の1層モデルの3つの比較ができる。面白そう。

1層モデルと8層モデルのたわみの誤差についての考察

1層モデルのポアソン比について考えてみる

• clthabas.f90：中央に輪載荷
• clthabab.f90：等分布載荷
• pclthaba.f90：gnuplot用にデータをまとめる

長方形板の4辺支持の式を計算

• 4hen.f90
• 等方性長方形板の4辺単純支持。構造力学公式集より。
• 2辺単純支持とかも作っていたが、たわみが発散する。式が複雑なので、多分どこか書き間違えているか、do文のところを間違えている。確認中。

供試体2(7層8プライ/ボルト有り)の1/4解析

• cltbyon.f90
• 木材利用研究論文報告集15のpp.17-22の有山さんのを真似る。つまり、輪荷重を受けた後の「一般部」の破壊試験(H28/09/05)のイメージ。
• 実験では降伏するまで載荷したものの、これといってCLTの損傷は見当たらなかったので、ボルトがどうにかなってしまったものと思われる。そこでボルトもモデル化。
• 計算してみたところ、随分とボルトの下部に応力が集中しているようだったので、弾塑性解析にしたら収束しないかもと思い、線形計算しかしていない。
• 直方体(CLT)と円柱(ボルト)が組み合わさったメッシュを生成するのは大変そうなので、ボルトは梁要素(B32)にしてみた。なので、ボルトの径と本数を変えての計算が容易にできる。最適な組み合わせを探すのにいいのでは？
• メッシュをかなり細かくしないと、ボルトの断面が膨らむ(梁要素なのに?)ので1/4解析にした。

• 静的試験について、エクセルのC列目には最大荷重が約10kNと書かれているが、『木材利用研究論文報告集15』のp.20には静的試験が100kN載荷と書かれている。W列目の単位がtfになっていることからも、エクセルの荷重(C列目)は単位がkNでなく、tfの間違いだと思う。エクセルの荷重が9.85$\times$9.8$\approx$97kNだと、FEMとおおよそ合う。

CLT橋

• nclt.f90
• 3点曲げ載荷．単純に，CLTと集成材の境界面の節点を共有させる．CLTと集成材で別々に材料定数を与えられるようにした．
• 境界の節点を共有させると，実験値よりもかなり剛性が高くなる(材料定数はかなりでたらめ(CLTのラミナ強軸E=5GPa, 集成材強軸E=7GPa)で計算したら，荷重P=150kNでたわみは$\delta$=30mmぐらい)．中村先生も言っていたが，実物はどちらかというと非合成桁に近いのだろう．
• ncltj.f90
• 3点曲げ載荷．CLTと集成材の間に，接合要素みたいなものを入れて，滑りを再現してみる．具体的には，接合要素（仮）をとても薄い直交異方性材料にして，せん断弾性係数のみを低くする(極端な食い込みを防ぐために直方向ヤング率は木材と同じ)．
• 接合要素のせん断弾性係数の値によって，計算結果が大きく（2倍ぐらい）変わるのだが，どのようにしてこの値を決めればいいのか分からないので，実験に合わせて接合要素のせん断弾性係数を決めるような感じになる．例えば，厚さ0.1mmの場合，$G_{xy}=G_{yz}=G_{zx}=\frac{E_{zz}/15}{50000}$とすると，たわみとひずみは実験値に結構近くなる．ただし，材料定数はでたらめ(CLTのラミナ強軸E=5GPa, 集成材強軸E=7GPa)なので，それほど意味はない．

• 上がCLT，下が集成材．150kN載荷時の変形を1倍で描くと，滑りを割と上手く再現できていることが分かる．
• ncltj_bun.f90
• ncltj.f90を等分布載荷にしたもの．
• A活荷重（13.5kN/m$^2$）が作用時のたわみは，中村先生の計算では$\delta$=20mmぐらいになるが，接合要素の諸元を上に書いたもので計算すると，$\delta$=18mmぐらいになる．
• ちなみに，『木材利用研究発表会講演概要集16』のp.83の最後の式（A活荷重でのたわみ）は，そのすぐ上にある見かけの剛性の式で$\ell=5800$mmを使っていることから，$\ell=6000$mmではなく$\ell=5800$mmだと思う．それで計算すると，たわみは$\delta$=18mmぐらい．
• ncltjb.f90<--需要がなさそうなので，つくっていない
• ラグボルトを梁要素として，CLTと集成材の節点上に配置する(cltbyon.f90と同じ方法で)．つまり，接合要素のせん断弾性係数をとことん低くして，せん断をラグボルトに受け持ってもらう．こうすれば，ラグボルトの本数によって滑りが変わってくれるのでは？

円筒

円筒(角筒)

• entou.f90
• 要素はS8.境界条件を色々と変えられるようにした。座屈モードが境界条件にとても敏感。
• hentou.f90
• 初期不整を与えて、幾何学非線形解析をすると、途中で計算が止まるものの、座屈しているらしい所までは見れる。

折り目を剛結とした回転螺旋折り

• gkaiten.f90
• 折り目は剛結。板一枚につき1要素なので、かなりメッシュが粗い。山内さんの計算の様に上部が広がる。
• s6x4.f90で要素分割を細かくするために、要素の並び順を一部修正。
• 幾何学非線形解析にすると、高さ0.87hぐらいで計算が止まってしまう。下図は、線形計算したもの。要素を細かくしての計算はまだやっていない。
• 段の高さに応じて高さ方向変位を多点拘束(*EQUATION)によって拘束している。

• 山内さんのモデルの作り方では、一段ごとにずらす角度を計算していないように見える。これだと、知らずに折り畳みにくいモデル($\alpha\neq\frac{\pi}{m}$)を作っていた可能性があるのでは？

折り目を剛結とした反転螺旋折り

• ghanten.f90
• 折り目は剛結。板一枚につき1要素なので、かなりメッシュが粗い。
• s6x4.f90で要素分割を細かくするために、要素の並び順を一部修正。
• 幾何学非線形解析にすると、高さ0.85hぐらいで計算が止まってしまう。下図は、線形計算したもの。要素を細かくしての計算はまだやっていない。
• 上下を板に挟まれた中の段は均等に畳まるが、片方が自由になっている上面と下面は変形が大きくなり、場合に寄ってはそこが裏返ったようになる。そこで、多点拘束で1段ごとの節点の高さ方向の変位の割合を一定にすると、(線形解析でも)高さ0までわりと綺麗に折り畳まる。

• 線形解析をすると、回転螺旋折は上部が拡がってしまい、反転螺旋折はあまり拡がらない。どちらも折り目が剛結なのにこういった差が出るのは、ヒンジがどうのという話ではない気がする。恐らく、前者は上部がぐるぐる回るつくりで、板要素の面に対して水平な(自由度を持たない?)方向の剛体変位(or回転?)が大きくなってしまうからだと思う。なので、段数が多くてぐるぐる回る回転螺旋折とかの場合、折り目をヒンジにしたとしても、幾何学非線形解析でないと拡がってしまうのではないか?
• Salome-Mecaに詳しそうな人の記事を見ていて、何となくそんな気がした。
• 後に後藤さんから聞いた話によると、板要素の誤魔化してる(面内回転の)自由度のことを"drilling degrees of freedom"というらしい。

板の要素数を増やす

• s6x4.f90
• 螺旋折りプログラムでは板一枚につき要素数1つしかできないので、それを4倍ずつ要素数を増やすプログラム。しかし、多点拘束を入れるのが面倒臭いので、そこは対応させていない。
• 多分、S6要素のみのモデルなら、何にでも使えると思う。ただし、要素を構成する節点番号の並び順を、要素を二等辺三角形とみなした場合の(等辺-->等辺-->底辺)にしないと要素の形状が細長くなってしまう。
• 元々は、構造解析関係のプログラム置き場にあるccxs6s3.f90を使おうとしていたが、なぜだか重節点になるところができてしまうので、これを参考に作りなおしてみた。ハニカムだとそれでもうまくいくのだけれど。

• 折り目を重節点にして多点拘束を入れてみたいが、頂点部分の節点が6つも重なるので、なかなか複雑。

• 立体要素の方が(遠藤さんの結果によると)幾何学非線形解析もうまくいくらしいし、簡単に折り目のヒンジをモデル化できそう。折り目を剛結とヒンジでやってみたい。

鋼棒のねじれ座屈

C3D8とC3D6の混在

• nezid6.f90
• ねじって線形座屈解析。軸方向に2000分割ぐらいしないと、座屈荷重が収束しない。境界条件も難しい。何回転もするようなモーメント荷重が入れられない(?)ので、1N･mに相当するねじれ角($\theta$=$T\frac{\ell}{GJ}$=$1\frac{\ell}{GJ}$)を与えて計算し、buckling factorがそのまま座屈トルク(N･ｍ)になるようにした。
• 石坂さんの($\ell$=3m,d=13mm,E=206GPa,$\nu$=0.3)を計算すると、理論値$T_{座屈}$=$2.861\times\frac{\pi}{\ell}\times EI$=865.3N・mに対して、CCXが861.7N･mぐらい。恐らく、円周方向の分割数が少ないため、断面積が円形より小さくなり、その分座屈トルクが理論値よりも小さくなっている。

• hnezid6.f90
• 例の如く、初期不整を与えてねじる。初期不整をとても小さくしないとC3D6部分が負のヤコビアンになる。つまり、要素が変てこな形になって計算ができない。しかし、このエラーが出ないように初期不正を小さくすると、今度は逆に座屈が起こらず、計算も途中で止まってしまう。

C3D8だけ

• nezid8.f90
• 初期不正を色々と工夫していると、どうやら初期不正が大きくても、C3D8の部分は負のヤコビアンになりにくいことが分かった。そこで、限りなく内径の小さい円管として、全要素をC3D8だけで作ることにした。何となく、こちらの方がいびつな形。石坂さんの問題の座屈トルクは861.7N･mになり、孔の影響はほとんどない。
• hnezid8.f90
• nonpositive jacobian は出なくなったが、途中で計算が止まってしまう。

• 気が向いたら、C3D20でもやってみたい。

B32

• nezib32.f90
• 梁要素。$T_{座屈}$=863.1N・mで、立体要素よりも理論値(865.3N・m)に近い。
• 後藤さんが『オイラー角を用いた簡潔な有限変位解析手法』で書いているように、座屈後挙動を追うには、梁中央で$u$と$\theta_x$を拘束する必要があるかもしれないので、その条件になるように改良した。
• hnezib32.f90
• 例のごとく、非線形座屈解析。座屈する瞬間は確認できるが、座屈後挙動は少ししか追えないので、残念ながらグリーンヒル問題を解けているとはいえない。

横ねじれ座屈の非線形解析

• 後藤ちゃちゃ(16/5/17)
• FEMで初期不正のある座屈シミュレーションから$P-\delta^{2}$法だのSouthwell法で座屈点を算定している例については、k2のbunken/g01j95enko1.pdf参照。
• あと、ローマ数字や◯数字の漢字コードは、昔は機種依存文字だったけど、最近は、UTF-8で扱えるようになったものの、このページの文字コードは実はEUC-JPなので、表示はできているものの、どうだろうか。
• とりあえず、半角数字にしました。
• 平沢先生の図面をk2のmokusinにアップしました。

• katasen.f-->harisen.fでやっているような座屈の計算を、CalculiXとかで工夫してやってみる。
• 汎用有限要素ツールは、固有値解析で座屈を一気に求めているので、途中の変形していく過程を考慮していないし、変形中のことが分からない(らしい)。
• かといって、複雑な形状とかを立体要素にモデル化して、真面目に計算するというのは、どうも現実的でなさそう(ここのp.4によると)。
• そこで、プリポストプロセッサーと連携しやすく、立体要素が使えて、材料非線形も考慮できるCalculiXとかでそれができたら、とても便利。
• 過去の横ねじれ座屈の卒論等を参考にした。
• 寸法：b=2cm,h=25cm,$\ell$=2.0m.
• スパンは一つしか試していない
• 分割：幅方向8,高さ方向30,軸方向200.
• 過去の卒論等の"幅方向6,高さ方向20,軸方向100"だと、線形の座屈荷重が3桁まで収束しないので、多めに分割(たわみはそれで収束するけど）。
• 材料定数
• 等方性:E=206GPa,$\nu$=0.4 .

手法

• 1.CalculiXで線形座屈解析をする。
• yoko.f90
• 2.1/5000(例えば)の初期不整を与える。
• 初期不整の大きさはそれなりの工夫が必要。
• 3.CalculiXで幾何学非線形解析をする。
• hyoko.f90

• gifのファイルサイズが大きいので一回消したのだが，遠藤さんが折り畳み円筒でこれをやりたがっていたので，参考になればと思い，再UPする．
• cgxppm.fとかで、色の最大値を調整したほうが綺麗かも。
• 1~5次モードでもいかにもそれらしい挙動を示す。異方性を入れても同様。
• 4.座屈荷重をP-$\delta^{2}$法とかで求める。
• pd2.f90
• sw.f90
• meck.f90
• P-$\delta^2$法：終わりのの17~5点を使って、線形回帰。

• Southwell法：初めの76点を使って、線形回帰。

• 終わりの方のプロットも使うと、座屈荷重が少し大きくなってしまう。横方向の変位だけ集めていて、軸方向に引っ込んで行くのを無視しているから?
• Meck法：初めの76点を使って、線形回帰。

• 終わりの方のプロットも使うと、座屈荷重が少し大きくなってしまう。回転軸が傾いていくのを無視しているから?

• 理論値はちょっと小さい。式の意味はほとんど理解していないが、スパンが短いからだと思う。そこは、もうちょっと勉強したい。
• 計算にはtrahair.f90を使った。

• なお、結果は省くが、オイラー座屈なども同様なかんじ。ただ、荷重の載荷方向を軸圧縮方向に変えればいいだけ。

伊藤さんのトラス橋

• 金具詳細モデル(鋼材:E=206GPa,$\nu$=0.3)を2次要素にして、Code_AsterとCalculiXで計算し直してみた。ついでに、向きごとに木材の異方性を入れてみた。

• 要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v_{C3D10}}$=-1.2%
• 異方性によるたわみの相対誤差(CalculiX(多分間違えている)):$\frac{v_{等方性C3D10}-v_{直交異方性C3D10}}{v_{直交異方性C3D10}}$=-12.6% --> -12.0%
• 初めは斜め部材の角度を60度にして計算したが、ちゃんと61.763度にして計算すると、ほんの少しだけ等方性の結果に近くなった。
• 異方性によるたわみの相対誤差(Code_Aster(あっていそう)):$\frac{v_{等方性4面体2次}-v_{直交異方性4面体2次}}{v_{直交異方性4面体2次}}$=-3.2%（60度でやった)
• ついでに、樹脂材料(E=2.84GPa,$\nu$=0.3)でもやってみた。

• 要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v_{C3D10}}$=-12.4% .金具部分が全体の剛性に与える影響が結構大きい。
• 異方性によるたわみの相対誤差(CalculiX(多分間違えている)):$\frac{v_{等方性C3D10}-v_{直交異方性C3D10}}{v_{直交異方性C3D10}}$=-13.7% --> -13.5%
• 異方性によるたわみの相対誤差(Code_Aster(あっていそう)):$\frac{v_{等方性4面体2次}-v_{直交異方性4面体2次}}{v_{直交異方性4面体2次}}$=-4.3%（60度でやった)
• ちなみに、2次元のトラス(toraITO.f90)や骨組(honeITO.f90)、ティモシェンコ梁の骨組(honetimo.f90)として計算すると、

• Marcを使った3次元トラスの計算は伊藤さんによるものなので、3桁までしか分からないが、2次元で計算してもそれなりに精度がいいことが分かる。
• 金具詳細モデルの計算ほどせん断の影響が大きくない。不思議。

• (2017/7)Code_Asterでも直交異方性材料として計算したところ，CalculiXよりもそれらしい値になった（トラスの部材に曲げがほとんど作用しないという点で）．多分，CalculiXでどこかを間違えている（木材利用2016のも怪しいので参考にしないでください．すみません）．確認中．
• 伊藤さんのトラス橋をやりたい人のために，斜めの直交異方性材料に直したcommファイル（kondo.comm）をk2のgakusei/2015/ito/sasikomi_model2/syousai_modelに入れておいた．geometryで木材を向きごとに4つのグループに分ける必要があるが，それは各自やってください．
• ただ，読み取っている箇所(中央の金具上の角材同士の中心線の交点)の変位より最大変位が2〜3倍も大きいので，あまり良い載荷方法ではないと思う．菊地さん(2015)みたいに板を乗っけるとか？
• 確か，Marcのトラス要素と比較するために，そのような載荷方法を採ったのだと思う．

CalculiX関連

unv-->inp

• unvc3d10.f90
• unvc3d4.f90を2次要素用に書き換えてみた。やはり、CalculiXでも2次要素?の方がメッシュが粗くてもたわみが収束しやすい。
• ついでに、節点群(*NSET)を座標指定でなくて、Salomeでつくった節点のグループから読み取れるようにした。
• mokuc3d10.f90(木材利用2016で使ったのはmokuc3d10_org.f90)
• unvc3d10.f90を複数材料・異方性対応にした。ついでに、トラスの斜材にも異方性を設定できるようにした。
• 局所座標系の設定に不確かなところ(SYSTEM=CYLINDRICAL)があったので、より安全なもの(SYSTEM=RECTANGULAR)に改良した。とはいえ、改良前と比べても、伊藤さんのトラス橋の計算結果は3桁一致するけど。
• jacobianエラー：輝さんの日誌によると、節点座標の出力形式を1pe17.9とかにすると出なくなる。それでも、複雑で細かい形状のものを4面体に要素分割すると、変てこな形(jacobianが負値をとる?)の要素がいくつかできてしまう。そういう要素は、(コメントアウトしてもそこの要素で穴が開くわけでないようなので)**でコメントアウトするなりして、手をかける必要あり。とすると、Code_Asterはそういった要素を自動ではじいてから計算しているのだろうか?
• menc3d10.f90
• mokuc3d10.f90を接触解析対応にした。つまり、salomeで作った面情報をinpファイルに書き込める。
• ccx_2.9だと問題ないが、ccx_2.6.1だとfrdファイルの要素情報に変な数値が書きこまれてしまいcgxで開けない。
• 他にも、モデルによっては面情報を上手く変換できていないことがあるので、絶対にどこか間違っている。そのうちになおすべし。
• zac3d10.f90
• mokuc3d10.f90をオンサイト木橋の線形座屈解析用に余計なところを消して、すっきりさせた。また、こちらの方がモデリング等が楽かと思い、節点群をグループ指定でなく、原型であるunvc3d4.f90のように座標指定に戻した。
• s6.f90
• s3.f90
• Salomeで切った3角形の板要素(unv)をCalculiX用入力データ(inp)に変換し，おまけで面積を算出するプログラム．Code_Asterの板要素の法線は，デフォルトでは作成したモデル通りではなく，(関塚さんのエラーメッセージによると)全ての要素の法線がx軸方向を向いてしまうらしい．つまり，面の向きごとに全てグループを作って，法線ベクトルを設定しなければいけない訳だが，ParaVis?ではモデルの形状が元のモデルのまま表示されるため，与えた法線ベクトルが正しいかどうかを確認する術がない．そこで，SALOMEで作ったメッシュをそのままinpファイルに変換してCalculiXで計算してやれば，その確認ができるのではないだろうか？という算段．ちなみに，ParaVis?(というかParaView?)は節点変位が3成分のベクトルになっていないと変位図が描けない(Warp by vector が使えない(Warp by scalar を使うのか？))ため，板要素と梁要素の変位が視覚的に確認できないという致命的な欠点がある．しかも，medファイルはバイナリなので，加工しようにも無理があるし…．
• 面内曲げの変位にS3とS6とでさほど結果に差はないのだが，面外曲げだと明らかに差がある．しかも，面内曲げの変位は近い値が出ると言ったものの，応力は1.5倍ぐらいの開きがある．ということで，Code_Asterはどうか知らないが，CalculiXで計算させるのなら，断然S6がおすすめ．

線形座屈解析から初期たわみを与える

• syokitawami.f90
• 「線形座屈解析用のinpファイル」と「全節点変位を出力したdatファイル」から「初期たわみを与えたinpファイル」を生成するプログラム。
• 何次の座屈モードであっても、どんな要素であっても、何種材料であっても気にせずに使えるので、非線形座屈解析の決定版とでもいいたいところ。
• ただ、inpファイルからコピー＆ペーストした部分の一列目に空白が入るのが、ちょっと気に食わない。何でこうなるのだろう？分かり次第直そう。

• Code_Asterでも，Medでなく，IDEAS(unv)などのテキストファイルに出力することができるので，それを利用して初期たわみを与えたunvファイルを生成し，meshに読み込むことが可能なはず．しかし，CalculiXよりも（フォーマットの理解だけでなく，(Salome自体の動作が重いため)生成したメッシュの確認が）確実に面倒臭いだろうから，余程の需要がない限りやらないことにする．

frd-->vtk

• c3d10vtk.f90
• c3d10のfrdファイルをparaview用にvtkファイルに変換する。cgxではできなかったコンター図の色の調整、特定の点・面の情報の出力、動画化などが簡単にできるようになる.主に非線形座屈の動画化や，木+鋼などの2種材料の応力を見る際にコンター図の色の範囲を変更して見やすくするといった目的で作る．
• frdファイルには書かれていない主応力を計算してvtkファイルにベクトルとして書き込めるようにしたので，praviewのGlyphを使えば主応力の大きさと方向を描ける．
• c3d8vtk.f90
• よく使うC3D8でも同じもの作る．
• s6vtk.f90を作ろうとしたのだが，S6はfrdファイルにおいて1要素15節点の中間節点を有する三角柱に置き換えられてしまうので，vtkフォーマットではできない．やるなら，ただの三角柱要素や中間節点を有する三角形に置き換えて誤魔化すとか？
• というか，vtkフォーマット以外を採用するのもいいかも．

• paraviewに計算結果を入れられると，こんな↓遊び？(CLTの最小主応力の向きと大きさ)もできて楽しい．

paraviewメモ

コンター図関連(Color Map Editor)

• ハートマーク(Choose preset)：初期設定の(青<-->赤)から虹色などの様々な色に変えられる．
• 個人的には，デフォルトの Cool to Warm が大変見づらく、不満だった。Rainbow Desaturated が最も見分けやすいと思う．
• ここ的な見やすさを求めるのなら，"Grayscale","Linear YGB","Linear Green","Linear Blue"辺りがいいと思う．
• Number Of Table Values：分割数．適度に荒い方が見やすい．
• 32分割ぐらいにすると、後述する Maximam Number of Labels に合わせられるので丁度いいと思う。
• C(rescale to custom data range)：コンター図の範囲の上下限値を設定．
• 目(rescale to visible data range)：コンター図の範囲を丁度良くしてくれる．
• e(Edit Color Legend Parameters)：カラフルな棒の大きさや向き，文字の大きさや数，フォント等を調整．
• Maximam Number of Labels を Number Of Table Values の数に合わせると、csvなどに出力しなくても大体の値が読み取れる。

あるグループのみ表示

• Filter -> Alphabetical -> Extract Group -> 見たいグループだけにチェックを入れる

csvに出力

• (Filters)Slice(あるいは Clip ) -> Select Points On(あるいは Select Points Through ?) -> (Filters)Extract Selection -> ロロ -> Spreadsheet View -> （右上の方にある）Show only selected Elements. -> (Fileの)Export Secen
• Select Points On : 表面の節点のみを選択
• Select Points Through : 内部の節点も選択
• Extract Group と組み合わせて(作ったグループだけ表示させて)、Select Points Through を使うのが良さそう。

動画化

• (File)Save Animation
• 細かな設定は試していないから知らない．
• その後，ImageMagick?のconvertを使う．

Salome-Meca関連

ある断面の軸方向応力分布を描画

• danmen.f90
• paravisから出力したデータをgnuplot用に整える。ついでに、断面の節点変位の平均値をデータファイルの行頭に出力する。
• ちなみに、断面をグルーピングし忘れたときでも、以下の方法で断面の情報を出力できるが、若干csvの列データが無くなるのでプログラムをいじる必要あり。
• パイプブラウザーのrmedファイルを選択したまま、フィルター->Common->Slice
• 見たい断面の、原点座標(origin)を入力し、面の法線(nomal)を選択。
• sankaku.f90
• $\sigma=0$の線とか$\sigma=\frac{M}{I}y$の線とかをgnuplot用につくる。
• szz.f90
• 'MED'ではなく'RESULTAT'に出力したデータを、gnuplot用に整える。rmedファイルとresultatファイルを同時に出力したい時は、ASTKを使う。danmen.f90の方法の方が楽。

Maximaリンク

• 斉木先生のページ
• ttp://mechanics.civil.tohoku.ac.jp/saiki/strcanalysis/kozomaxima.pdf

ImageMagick?メモ

複数のスクリーンショットを一括でppm(アスキー)に変換

mogrify -format ppm -compress none *png

• cgxppm.f関連で使えないかと。

定点観測で撮影した連続写真を重ね合わせ

convert *.JPG -evaluate-sequence max gousei.png

• 場合によっては、コマンド中のmaxをminにした方がいい。
• 特に意味はないが、上のコマンドではデジカメの画像形式からpngに変えている。

• この画像をImageJ(Windowsの場合は公式ページからダウンロード)で読み込んで、
• 変位を測りたい場合は、
• Straght:距離が既知な線に沿って、線を引く。(ドラッグ&ドロップみたいにして)
• Analyze-->Measure:線の長さをpixelで測る。
• Analyze-->Set Scale-->Known distance:既知な線の実際の長さを入力する。(これで写真上の長さ[pixel]が実際の長さ[mm など]に置き換わる)
• Analyze-->Measure:線の長さを実際の長さで測る。
• たわみ角(傾き)を測りたい場合は、
• Angle tool:測りたい角度をなす3点を順にクリック
• Analyze-->Measure
• よく使いそうな保存の形式は、
• File-->Save as-->Results:結果をテキストで出力。
• File-->Save as-->XY Coordinates:現在、表示されている線の座標をテキストで出力する。

• 野田さんと3年生のたわみ試験で試してみたところ、精度がとても悪かった。(誤差が±10%以上あった気がする。ちなみに、上の合成写真はその実験前に一人で行ったものなので、また別の話。)
• 原因としては、
• 1.三脚のカメラを固定するところが壊れているので、シャッターを切るときにカメラが動く。
• これは、ちゃんとした三脚を使えばいいだけの話。
• 2.被写体と1m近く離したのだが、たわみが小さい(8mm程度)こととカメラの解像度(そんなに悪くない気がする)のせいで、ImageJ上で点を打とうと拡大すると結構ぼやける。
• 3.なるべく被写体の真正面にカメラを据え付けたつもりだが、下手だったかもしれない。被写体が画面の端に行くほど歪んでしまうが、被写体とカメラの距離が近いためその影響も大きそう。
• 原因3と原因2は、どちらかを優先させると、もう片方が疎かになるため、なかなか調度よい距離が掴めない。一番いい方法は、被写体からなるべく程度距離を離して、ものすごく解像度の高いカメラを使うとか？

• 同時にRegard3DというWindows+Mac用のオープンソースのソフトでSfMもやってみたのだが、15〜16枚程度の写真ではほとんど特徴点が抽出されなかった。アクリル板がカラフルじゃないと駄目なのかもしれない。
• 特徴点がたくさん抽出できる場合には、 1.Regard3D(SfM + メッシュの作成) -> 2.MeshLab?(余計な点の除去) -> 1〜２を変形後の写真でも行い -> 3.FreeCAD(長さを測る) が使えることを確認済み。
• ネット上にあったカエルのぬいぐるみや猿のぬいぐるみで試した。この時は、8〜10枚程度の（横480pixelくらい？の解像度がそれほど高くない）画像でできた。
• ただし、百目石橋でもそれをやろうとすると、相当いいPC(ゲーミングPCとか)でないとできない。

• 野田さんによると，tiffだと拡大してもぼやけにくいとか．例えば，
• convert -threshold 50% -edge 1 hoge.JPEG hoge.tif
• などとすると，輝度50%を境(50%でなくてもいい)に2値化し，1画素の線(2や3でもいい)でエッジ抽出し，tiffに変換されるので，以前よりは見やすいかも．ただし，後からtiffにしたところで意味があるかどうかは謎．カメラの設定でtiffにするのと違うのだろうか？

pdfを変換するときに画質を下げない

$\LaTeX$メモ

minipageのキャプション

単位

$\frac{P}{A} = \SI{100}{N/mm^2}$

ダブルクォーテーション

``hoge''

画像をわざと表示させない

\documentclass[a4j,12pt,dvipdfmx,draft]{jarticle}

$\LaTeX$周りのツール群

pdfcrop

• 余白無し

pdfcrop hoge.pdf hoge_c.pdf

• 左右・上下に余白

pdfcrop --margins "10 10" hoge.pdf hoge_c.pdf

• マージンは"左右　上下"の順．
• 左・上・右・下に余白

pdfcrop --margins "10 10 10 -50" hoge.pdf hoge_c.pdf

• マージンは"左　上　右　下"．

• ちなみに，出力ファイル名を指定しない場合，勝手にhoge-crop.pdfを出力してくれる．

Ghostscriptによるpdfのグレースケール化

gs -sOutputFile=hoge_g.pdf -sDEVICE=pdfwrite -sColorConversionStrategy=Gray -dProcessColorModel=/DeviceGray? -dCompatibilityLevel=1.4 -dAutoRotatePages=/None -dNOPAUSE -dBATCH hoge.pdf

• ここで，1つ目のオプションの"hoge_g.pdf"は出力ファイルなので好きな名前に，最後のオプションの"hoge.pdf"は変換したい入力ファイルに変える．
• また，"?"はwikiの機能で勝手に付いてしまっているだけなので，消すこと．
• 他の細かなオプションについては分かりません．
• 参考URL：
• http://takuya-1st.hatenablog.jp/entry/2014/09/22/010805

おまけ（修論の背表紙）

• se.tex
• 用紙の向きが縦書き．余白や空白を適宜調整する必要あり．
• sel.tex
• 用紙の向きが横向き．タイトルが短いため横向きにしてみた．dvipdfmx -l の"-l"のつけ忘れに要注意．
• 例
ちなみにfirefoxだと，縦書きの「ー」は横向きのままのようだが，evinceだとちゃんと|みたいに縦向きになる．

その他，便利なコマンドライン

pdfimages

pdfimages -png hoge.pdf out

grep

特定の文字を含むファイルを検出

grep ほげ -rl .

3Dプリンター メモ

• ボールジョイントも作れそう(PRN3DでなくてEDENでの話)。
• http://rayerd.sblo.jp/article/172174211.html
• 造形後に関節を折り曲げることで、隙間につまったサポート材を壊している。
• http://rayerd.sblo.jp/article/172249880.html
• http://rayerd.sblo.jp/article/172529855.html

gnuplotメモ

日本語のpngを出力する

set term png font "/usr/share/fonts/truetype/takao-gothic/TakaoPGothic.ttf"  size 500,350

nkf -guess hoge.plt

nkf -w-overwrite hoge.plt

Fortran90(gfortran)メモ

固有値計算

• jacobi.f90
• ヤコビ法によって実対称行列の固有値と固有ベクトルを計算するサブルーチン．どうやらヤコビ法で求めた固有ベクトルは正規化されているようだ．
• frdファイルをvtkファイルに変換しようとして，frdファイル中に主応力が書かれていないことを知る．主応力もparaviewで見れると色々と面白いかもしれないので，応力テンソルを自分で対角化してみることにした．というか，コンクリートの教科書にありそうな，場所ごとの最大主応力の大きさと向きの図を描いてみたかっただけ．
• scipyやoctave,julia等を使えば一行で簡単に固有値と固有ベクトルを求めることができるけど，最初から組み込み関数に頼るのは面白くないし，(固有値以外の計算も含めて)計算速度が遅い．LAPACKの整備?の仕方が分からないというのも作った理由の一つではあるけど．
• sjacobi.f90
• jacobi.f90を昇順に並び替えるようにした．

LAPACK

• Ubuntu16.04LTSでのインストール方法
• https://qiita.com/AnchorBlues/items/69c1744de818b5e045ab
• 全くこの通りで整備できた(LAPACKとLAPACK95)．
• コンパイルの仕方やお試し計算
• http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux/fortran.html#sec7

FFTW

• Ubuntu16.04LTSでのインストール方法
• http://moznion.hatenadiary.com/entry/20120718/1342626456
• 全くこの通りで整備できた(リンク先はUbuntu12.04だけど)．
• コンパイルの仕方やお試し計算
• http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux/fortran.html#sec8

Julia メモ

Ubuntu上で最新版にする

sudo add-apt-repository ppa:staticfloat/juliareleases
sudo add-apt-repository ppa:staticfloat/julia-deps
sudo apt-get update
sudo apt-get install julia

sudo apt-get autoremove julia

パッケージのインストール

Pkg.add("Plots")

• ちなみに，ここでの例のPlotsは何だかplotly.jsと同じみたい．plotlyは対話的に動かせる高機能なツールなので普段使いにはとても便利だけど，論文・スライドに使えるほど綺麗ではないと思う．
• JuliaはPGFPlotsにも対応しているようなので，やはり論文などにはそのまま$\LaTeX$に取り込めるPGFPlotsを使うのが綺麗でいいのでは……．と思いきや，なかなかPGFPlots.jlを使えるようにならない．

Python メモ

コメント文にutf-8の日本語を使えるようにする

# -*- coding: utf-8 -*-

Matplotlibで日本語を使えるようにする

mv ~/.cache/matplotlib/fontList.cache ~/.cache/matplotlib/fontList.cache_org (←fontList.chaceをリセットしたいだけなので、rmで消してもいいけど、一応残す)
sudo cp /usr/share/fonts/truetype/takao-gothic/TakaoPGothic.ttf /usr/share/matplotlib/mpl-data/fonts/ttf/ (←やる必要がなかったかも！)
sudo vi /etc/matplotlibrc

Rustメモ

rust.vim

• https://github.com/rust-lang/rust.vim/blob/master/syntax/rust.vim
• を"/usr/share/vim/vim74/syntax/"にダウンロードする．
• 次に"sudo vi /usr/share/vim/vim74/filetype.vim"で，どこか(アルファベット順に並んでいるのでrubyの下とか?)に

" Rust
au BufNewFile?,BufRead? *.rs                      setf rust

• を追加し，拡張子が被っているHerculesをいじってごまかす(本当はscripts.vimというのを書いてファイルの種類を判別するのがいいらしいが，Herculeを全く知らないので判別方法を考えるのが面倒臭かった)．

" Hercules
" au BufNewFile?,BufRead? *.vc,*.ev,*.rs,*.sum,*.errsum   setf hercules
 au BufNewFile?,BufRead? *.vc,*.ev,*.sum,*.errsum setf hercules

• どうやらvimには拡張子(言語)別に辞書機能のようなものや自動インデント機能などがあるそうで世の中の人は結構使っているらしいが，まあ，それは困っていないので別にいいや．

Goメモ

パッケージのインストール

go get github.com/〇〇

• "github.com/○○"はURL．

• 数値計算系ライブラリの有名どころを探してみると，それなりに揃っている．
• BLAS/LAPACK:github.com/gonum/gonum
• FFTW:github.com/runningwild/go-fftw

• ただし，これはこの方法でgolangをインストールした時の話なので，apt-getでインストールした場合にはライブラリもapt-getを使わないといけないのかも？

フォーマッタ

go fmt hoge.go

Haskell(GHC)メモ

パッケージのインストール

cabal update
cabal install hmatrix

METAFONTメモ

gftodvi

• mf hoge.mf と打つと，hoge.2602gf というファイルができる．
• METAFONTについては，クヌース先生が書いた本（日本語訳）が図書館にあるので，それを参考に．
• 他の参考ページ：ttp://hayato1224.web.fc2.com/MFmain.pdf

• gfファイルをdviファイルに変換するために，gftodviというコマンドを使うらしいのだが，ただ gftodvi hoge.2602gf と打っても，

gftodvi: fatal: tfm file `gray.tfm' not found.

• というエラーが出てしまい，生成されたdviファイルは0バイトとなってしまう．そこで，意味はわからないが，予め

mktextfm gray

• と打っておくと，gftodvi hoge.2602gf でgfファイルをdviファイルに変換できるようになる．
• あとは，xdvi hoge.dvi& でも開けるし，dvipdfmx hoge.dvi とすれば，pdfファイルも生成される．

mktextfm

mktextfm hoge

Vimメモ

ある単語の数を数える

:%s/hoge/&/gn

シンタックスハイライトの追加

• "/usr/share/vim/vim74/syntax"にGitHub?などからhoge.vimを追加．
• "/usr/share/vim/vim74/filetype.vim"を編集し，拡張子から言語を自動判別するように設定．
• もしかすると，同じ拡張子を使っている言語もあるかもしれないので注意!(scripts.vimで自動判別するようにもできるとか)

USBメモリにバックアップ

rsync -auv /home/kouzou/2016 /media/disk

• オプション"-auv"は，「a:パーミッションを維持，u:差分のみ，v:処理を表示」という意味．
• "disk"は、USBメモリに付けた名前。
• ちなみに，k2へのバックアップ時にも使える(パスを間違えそうで怖いけど)．

HHKB Pro2

• システム -> 設定 -> ハードウェア -> キーボード -> レイアウト -> キーボードの型式 -> Happy Hacking Keyboad
• 右上のキーボードマークを右クリック -> 設定 -> + -> US -> (△で)USを日本語よりも上にする．
間違いなく打ちやすいが，日本語キーボードを普段使っていると，記号の位置に迷う． とはいえ，じっと眺めているとUS配列の方が記号の位置がきれいに思える．

講義関連

LibreOffice?

• https://help.libreoffice.org/Main_Page/ja

構造力学特論

• torasu2.f90
• hone2.f90
• 格点を剛結にするとどのくらい曲げが入るのか気になったので、torasu.f90やhone.f90をいじって各部材ごとの節点力(全体座標系のまま)を出力するようにした。
• toraITO.f90
• honeITO.f90
• 伊藤さんの木製トラス橋を２次元の問題として解いてみたかったので、要素数を9から15まで増やせるようにした。とはいっても、ただ配列の宣言をいじっただけ。
• honetimo.f90
• さらに、伊藤さんの木製トラス橋の2次元骨組を剛性行列をティモシェンコ梁(岩熊先生のページ)のものに置き換えてみた。

土質工学特論

• 間隙水圧分布を計算:u.f90
• set zrange [] reverse : z軸を下向き正にする。
• set pm3d : z軸の値を色付けする。x軸でこれをやりたい。調べるべし。

TA関連

基礎物理学実験

• 指導書
• 基礎物理学実験TA研修システム
• Ubuntu12.04のFirefoxでWindows Media Playerを再生するには、
• gnome-mplayer
• gecko-mediaplayerSynaptic
• をパッケージマネージャーでインストールすると、エラーは出るけど見れる。

土木環境工学実験(アクリルの引張試験)

• ヤング率とポアソン比の確認:yangu.f90
• syokignu.fでデータ整形後

時間割と単位

• 「関連性のあるコースで開講する専門科目」とは：理工学研究科の専門科目なら、どれでもいいらしい（2016年4月に理工学研究科担当の人に確認）。
• しかし，2017年度入学の人は，ガイダンスで，システムデザイン工学専攻以外の科目は認めない（つまり機械の専門科目をとれということ）と言われたらしい．心配なので2017年4月に理工学研究科担当の人（2016年度とは違う人に変わった）にもう一度訊いてみると，2016年度入学の人は理工学研究科の1期生で，大学の方もよくわかってなかったからそれでいいよ，とのこと．
• 2017年度に取得予定の残り12単位の内訳は、
• システムデザイン工学課題研究：専門必修(2年間)=10単位
• システムデザイン工学演習：専門必修(通年)=2単位

2016年度

• 理工学デザイン：1共通必修(後期前半)
• システムデザイン工学課題研究：5専門必修(通年)

2016年度前期

2016年度後期