近藤の修論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
$雪荷重を受けるプレストレス木箱桁橋の座屈挙動$
$プレストレス木箱桁橋鋼板部の局部座屈挙動（仮）$
-桁高やスパンが変わったときに，座屈挙動がどう変わるか？
M1のとき:$汎用有限要素ツールを利用した非線形解析（案）$
*ちょっと気になったこと（2018/02/15．もう，遅い…） [#j306...
対称条件から1/4モデルで座屈解析を行ったが，「支点付近の座...
しかも，座屈が幅員方向に非対称なモードとなる場合は，対傾...
ただ，十分な精度が得られる要素分割で1/2解析なんて，対傾構...
小さいモデルで確認しておきたい．
*フランジ付きモデルの材料定数の影響 [#d3d3211d]
一体化方法の影響について計算してしまったので，木材を等方...
-対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, ...
, 木材の材料 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 直交異方性 , 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 2...
, 等方性 , 27.2, 支間中央1・2 , 27.7, 支間中央1~3 , 28.5 ...
-対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm,...
, 木材の材料 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 直交異方性 , 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26...
, 等方性 , 27.5, 支間中央1・2 , 27.8, 支間中央1~3 , 28.5 ...
-やはり，異方性の影響は1割ぐらいある．木材を等方性材料に...
*フランジ付きモデルの一体化方法の影響 [#b6442f90]
木材利用2017で，部分的非一体化モデルがいいと言っておきな...
-対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, ...
, 非一体化方法 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) ,...
, 全体的 , 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 25.9 ...
, 部分的 , 25.4, 支間中央1 , 27.5, 支点付近 , 30.7 , 支間...
-対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm,...
, 非一体化方法 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) ,...
, 全体的 , 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 ,...
, 部分的 , 26.5, 支間中央1 , 31.7, 支間中央1 , --- , --- ...
-一体化方法によって座屈モード(PC鋼棒間をまたぐかどうか)が...
--フランジなしモデルでは，一体化方法による影響が1割もあっ...
-フランジ付きモデルだと，どう一体化しても，結局すり抜ける．
--フランジなしモデルの座屈モードは，部分的非一体化モデル...
-実は部分的非一体化モデルの方が木材と鋼板の接している面積...
-全体的非一体化モデルにしても，飛び出しが生じないモデル化...
--木材利用2017当時は，全体的非一体化モデルにすると飛び出...
*対傾構の影響 [#n60c922d]
-大池のモデルで，対傾構の数を変化させる．
--1.対傾構なし(0)
--2.端部・支間中央(2)
--3.さらに間に1枚入れる(2+1=3つ)
--4.さらに間に1枚入れる(3+2=5つ)
-対傾構を入れると，（境界の木材をSub Meshで分割したせいも...
--モデル化を工夫する(&link(onsitz.py,http://www.str.ce.ak...
-寸法は大池モデルにし，1/4解析を行う．
--たわみ，軸方向応力について共に最大値のみ参考値としてメ...
, 対傾構の数 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 0, 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 25.9 , 支間...
, 2, 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 , 支間...
, 3, 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 , 支間...
, 5, 25.3, 支間中央1~2 , 25.9, 支間中央1~3 , 27.5 , 支間...
-[考察]
--対傾構が直接効いてくるのは支点付近の座屈だと思うが，そ...
--しかし，対傾構をまたいで座屈するモードは算出されていな...
--対傾構があると断面変形が抑えられるため，最大たわみは当...
*まとめ(橋長・桁高・孔を変化させた時)(対傾構を2つを入れて...
3.5kN/m$^2$の等分布荷重載荷時の座屈解析の結果．
-とりあえず，橋長については，約10mから約21mmまで(虎毛と大...
-とりあえず，桁高についても，500mmから800mmまで(虎毛と大...
--鉛直材の挙動を見るために，1000mmまで増やしてもいいのだ...
-鉛直材・垂直材の幅は80・85・90・95mmの組み合わせが16通り...
--ということで，一番きついであろう80mmx80mmを始めにやって...
--80mmx80mmとの比較のために，一番変化が分かりやすい95mmx9...
--CalculiXみたいに夜通し連続して計算できれば，パターン数...
---ASTKがSalome-Meca2017でなくなってしまったかと思ってい...
-表中の支間中央No.のNo.は，PC鋼棒間を支間中央から順に数え...
-表中の鉛直材No.のNo.は，鉛直材を端部から順に数えたもの．
-たわみと軸方向応力について，せっかく計算したので何かの足...
--当たり前だが，橋長が短いほど断面変形が大きく，最大のた...
-固有値が近接しているため，何次モードまで計算するべきか迷...
**桁高500mm [#jc7d20ca]
***鉛直材80mm・斜材80mm [#j968097b]
-橋長10.220mで，いくつか飛び出しモードが出てしまった．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 58.4, 支間中央1~3 , 59.0, 支間中央1~3 , ---- , ...
, 12.410, 39.2, 支間中央1~3 , 40.3, 支間中央1~3 , 41.0, ...
, 14.600, 28.4, 支間中央1~4 , 28.9, 支間中央1~4 , 29.5, ...
, 16.790, 21.4, 支間中央1~3 , 21.8, 支間中央1~4 , 22.4, ...
, 18.980, 16.7, 支間中央1~3 , 17.1, 支間中央1~5 , 17.2, ...
, 21.170, 13.4, 支間中央1~3 , 13.6, 支間中央1~4 , 13.9, ...
***鉛直材95mm・斜材95mm [#uf3b6c41]
-橋長10.220mで，飛び出しモードばかりが出てしまい，なかな...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 60.9, 支間中央1~2 , 61.4, 支間中央1~3 , ---- , ...
, 12.410, 41.3, 支間中央1~3 , 41.5, 支間中央1~3 , 42.5, ...
, 14.600, 29.8, 支間中央1~4 , 29.9, 支間中央1~4 , 30.8, ...
, 16.790, 22.5, 支間中央1~3 , 22.7, 支間中央1~4 , 23.1, ...
, 18.980, 17.5, 支間中央1~4 , 17.6, 支間中央1~4 , 17.9, ...
, 21.170, 14.1, 支間中央1~4 , 14.2, 支間中央1~4 , 14.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#t8bc4aab]
-やはり，木材と鋼板の境目で木材側の要素のアスペクト比を小...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 60.0, 支間中央1~4 , 62.1, 支間中央1~4 , 63.5 , ...
, 12.410, 39.7, 支間中央1~4 , 40.7, 支間中央1~3 , 42.3, ...
, 14.600, 28.6, 支間中央1~4 , 29.5, 支間中央1~5 , ---- , ...
, 16.790, 21.5, 支間中央1~4 , 22.0, 支間中央1~4 , 22.5, ...
, 18.980, 16.8, 支間中央1~5 , 17.2, 支間中央1~6 , 17.4, ...
, 21.170, 13.4, 支間中央1~4 , 13.7, 支間中央1~5 , 13.9, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#i74314c3]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 61.9, 支間中央1~3 , 64.2, 支間中央1~4 , 66.0, ...
, 12.410, 41.7, 支間中央1~4 , 42.5, 支間中央1~4 , 43.4, ...
, 14.600, 30.0, 支間中央1~4 , 30.8, 支間中央1~5 , 31.6, ...
, 16.790, 22.6, 支間中央1~4 , 23.0, 支間中央1~4 , 23.5, ...
, 18.980, 17.6, 支間中央1~5 , 17.9, 支間中央1~5 , 18.4, ...
, 21.170, 14.1, 支間中央1~4 , 14.3, 支間中央1~5 , 14.5, ...
**桁高600mm [#o6d177d7]
***鉛直材80mm・斜材80mm [#qaced9ef]
-橋長10.220mで，飛び出し発生!(3次モードは得られたけど)
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 61.4, 支点+鉛直材1 , 62.3, 支間中央1~3 , 64.9, ...
, 12.410, 42.5, 支間中央1~3 , 44.7, 支間中央1~4 , 46.6, ...
, 14.600, 30.8, 支間中央1~3 , 32.2, 支間中央1~4 , 32.9, ...
, 16.790, 23.2, 支間中央1~3 , 24.1, 支間中央1~4 , 25.0, ...
, 18.980, 18.1, 支間中央1~4 , 18.8, 支間中央1~5 , 19.3, ...
, 21.170, 14.5, 支間中央1~3 , 15.0, 支間中央1~4 , 15.4, ...
***鉛直材95mm・斜材95mm [#s71a8d20]
-橋長10.220mで，飛び出し発生! というか，このケースのほと...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 65.2, 支間中央1~3 , ---- , 飛び出し! , ---- , ...
, 12.410, 44.2, 支間中央1~3 , 46.7, 支間中央1~4 , 47.8, ...
, 14.600, 32.0, 支間中央1~3 , 33.4, 支間中央1~4 , 34.3, ...
, 16.790, 24.1, 支間中央1~3 , 25.2, 支間中央1~4 , 25.9, ...
, 18.980, 18.8, 支間中央1~4 , 19.5, 支間中央1~5 , 20.0, ...
, 21.170, 15.1, 支間中央1~3 , 15.6, 支間中央1~4 , 16.1, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#gbaf52d5]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 64.8, 支間中央1~4 , ---- , ---- , ---- , ---- ,...
, 12.410, 43.5, 支間中央1~3 , 45.2, 支間中央1~4 , 47.2, ...
, 14.600, 31.5, 支間中央1~4 , 32.8, 支間中央1~5 , 34.4, ...
, 16.790, 23.6, 支間中央1~4 , 24.4, 支間中央1~4 , 25.1, ...
, 18.980, 18.4, 支間中央1~5 , ---- , ---- , ---- , ---- ,...
, 21.170, 14.7, 支間中央1~5 , 15.1, 支間中央1~5 , 15.5, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#r9fef15b]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 66.9, 支間中央2~4 , 72.3, 支間中央1~4 , 72.6, ...
, 12.410, 45.5, 支間中央1~3 , 47.4, 支間中央1~4 , 48.8, ...
, 14.600, 32.8, 支間中央1~4 , 34.4, 支間中央1~5 , 35.1, ...
, 16.790, 24.6, 支間中央1~4 , 25.5, 支間中央1~5 , 26.2, ...
, 18.980, 19.2, 支間中央1~5 , 19.9, 支間中央1~6 , 20.6, ...
, 21.170, 15.3, 支間中央1~5 , 15.8, 支間中央1~5 , 16.2, ...
**桁高700mm [#b1b28363]
***鉛直材80mm・斜材80mm [#k96f9889]
-橋長10.220mで，いくつか飛び出しモードが出てしまった．
--Sub Mesh の優先順位を1:境界面 2:木材 に変えた方がいいか...
--double mesh がどうのこうのと書かれているので，remove do...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 49.3, 支点+鉛直材1 , 63.6, 支点+鉛直材1 , ---- ...
, 12.410, 40.7, 支点+鉛直材1 , 44.2, 支間中央1・2 , 46.1,...
, 14.600, 31.9, 支間中央1~3 , 33.1, 支間中央1~4 , 33.9, ...
, 16.790, 24.1, 支間中央1~3 , 25.0, 支間中央1~4 , 25.8, ...
, 18.890, 18.8, 支間中央1~3 , 19.5, 支間中央1~4 , 19.9, ...
, 21.170, 15.1, 支間中央1~3 , 15.5, 支間中央1~4 , 16.0, ...
***鉛直材95mm・斜材95mm [#t274bf45]
-やはり橋長10.220mは飛び出しが多い．経験的に座屈しにくい...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 55.0, 支点+鉛直材1 , 63.4, 支点+鉛直材1 , 66.0 ...
, 12.410, 45.0, 支間中央1・2 , 45.5, 支点+鉛直材1 , 47.3,...
, 14.600, 32.4, 支間中央1~3 , 33.9, 支間中央1~4 , 34.8, ...
, 16.790, 24.6, 支間中央1~3 , 25.6, 支間中央1~4 , 26.5, ...
, 18.890, 19.2, 支間中央1~3 , 19.9, 支間中央1~4 , 20.4, ...
, 21.170, 15.4, 支間中央1~3 , 15.9, 支間中央1~4 , 16.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#lf7e27f6]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 67.4, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 12.410, 45.5, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 14.600, 32.9, 支間中央1~4 , 34.0, 支間中央1~5 , ---- , ...
, 16.790, 24.7, 支間中央1~4 , 25.3, 支間中央1~4 , 26.0, ...
, 18.890, 19.3, 支間中央1~5 , 19.9, 支間中央1~5 , 19.8, ...
, 21.170, 15.4, 支間中央1~4 , 15.7, 支間中央1~5 , 16.1, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#mc6a2237]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 69.4, 支間中央2~3 , 73.2 , 支間中央2~3 , 76.0, ...
, 12.410, 46.6, 支間中央1~3 , 47.9, 支間中央1~4 , 50.3, ...
, 14.600, 33.7, 支間中央2~4 , 34.9, 支間中央2~5 , 36.9, ...
, 16.790, 25.3, 支間中央1~4 , 25.8, 支間中央1~4 , 26.6, ...
, 18.890, 19.7, 支間中央1~4 , 20.2, 支間中央2~5 , 20.9, ...
, 21.170, 15.8, 支間中央1~4 , 16.1, 支間中央1~5 , 16.5, ...
**桁高800mm [#ib28b7b4]
***鉛直材80mm・斜材80mm [#abad02dc]
-''ここの表だけ木材の幅を420mmに間違えて(0.1mm削ると419.9...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 39.3, 支点+鉛直材1 , 51.1, 支点+鉛直材1・2 , 54...
, 12.410, 32.5, 支点+鉛直材1 , 41.9, 支点+鉛直材1・2 , 43...
, 14.600, 27.7, 支点+鉛直材1 , 32.0, 支間中央1・2 , 33.0,...
, 16.790, 24.2, 支点+鉛直材1 , 24.5, 支間中央1~2 , 25.2, ...
, 18.980, 18.9, 支間中央1~3 , 19.5, 支間中央1~4 , 19.8, ...
, 21.170, 15.4, 支間中央1~3 , 15.7, 支間中央1~4 , 16.0, ...
***鉛直材95mm・斜材95mm [#w6ee712b]
-橋長10mは飛び出しモードも算出された．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 44.3, 支点+鉛直材1 , 62.4, 支点+鉛直材1・2 , 64...
, 12.410, 36.6, 支点+鉛直材1 , 44.8, 支点+鉛直材1・2 , 43...
, 14.600, 31.3, 支点+鉛直材1 , 32.1, 支間中央1・2 , 33.1,...
, 16.790, 24.5, 支間中央1~3 , 25.2, 支間中央2~3 , 26.0, ...
, 18.980, 19.0, 支間中央1~3 , 19.5, 支間中央1~4 , 19.9, ...
, 21.170, 15.4, 支間中央1~3 , 15.7, 支間中央1~4 , 16.1, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#f08f6c00]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 54.8, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 12.410, 44.1, 鉛直材2~3 , 45.8, 鉛直材1~5 , ---- , ----...
, 14.600, 33.3, 支間中央2~3 , ---- , ---- , ---- , ---- ,...
, 16.790, 25.2, 支間中央1~3 , 25.3, 支間中央1~3 , ---- , ...
, 18.980, 19.6, 支間中央2~3 , 19.9, 支間中央2~4 , 20.4, ...
, 21.170, 15.8, 支間中央1~3 , 15.8, 支間中央1~4 , 16.1, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#s8e16104]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 68.6, 支間中央2~3 , 71.9, 支間中央2~4 , ---- , ...
, 12.410, 46.4, 支間中央1~2 , 47.3, 支間中央1~3 , 49.4, ...
, 14.600, 33.4, 支間中央2~3 , 34.4, 支間中央2~4 , 35.9, ...
, 16.790, 25.3, 支間中央1~3 , 25.5, 支間中央1~4 , 26.2, ...
, 18.980, 19.6, 支間中央2~3 , 20.0, 支間中央2~4 , 20.6, ...
, 21.170, 15.8, 支間中央1~3 , 15.9, 支間中央1~4 , 16.2, ...
**桁高900mm [#i96aa86f]
***鉛直材80mm・斜材80mm [#f36b9e88]
桁高900mmから計算時間が5倍ぐらい長くなった．多分，要素分...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 32.4, 支点+鉛直材1 , 42.3, 支点+鉛直材1・2 , 45...
, 12.410, 26.8, 支点+鉛直材1 , 34.7, 支点+鉛直材1~3 , 36....
, 14.600, 22.9, 支点+鉛直材1 , 29.5, 支点+鉛直材1~3 , 30....
, 16.790, 19.9, 支点+鉛直材1 , 24.3, 支間中央1~3 , 25.0, ...
, 18.980, 17.6, 支点+鉛直材1 , 18.5, 支間中央1・2 , 19.2,...
, 21.170, 15.3, 支間中央1~3 , 15.6, 支間中央1~4 , 15.8, ...
***鉛直材95mm・斜材95mm [#k2393a16]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 36.7, 支点+鉛直材1 , 52.0, 支点+鉛直材1・2 , 62...
, 12.410, 30.4, 支点+鉛直材1 , 42.9, 支点+鉛直材1~2 , 43....
, 14.600, 26.0, 支点+鉛直材1 , 20.8, 支間中央1~3 , 32.3, ...
, 16.790, 22.9, 支点+鉛直材1 , 23.9, 支間中央1~3 , 24.7, ...
, 18.980, 18.2, 支間中央1~3 , 18.9, 支間中央1~5 , 19.5, ...
, 21.170, 15.0, 支間中央1~3 , 15.4, 支間中央1~5 , 15.8, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#m49094a7]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 45.7, 鉛直材2~3 , 47.7, 鉛直材1~4 , ---- , ----...
, 12.410, 36.6, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 14.600, 30.7, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 16.790, 24.8, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 18.980, 19.3, 支間中央2~4 , ---- , ---- , ---- , ---- ,...
, 21.170, 15.5, 支間中央1~4 , 15.8, 支間中央1~5 , ---- , ...
-ちなみに，特に橋長10mなどの短いものは，変位が最大となる...
***鉛直材95mm・斜材95mm(対傾構2つ) [#a5b90213]
-橋長10mの2次モードは，鋼板全体が波々に座屈する．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 66.7, 支間中央2~4 , 69.9, 鋼板全体 , ---- , ---...
, 12.410, 44.9, 支間中央1~3 , 46.9, 支間中央1~4 , 48.8, ...
, 14.600, 32.5, 支間中央2~4 , 33.8 , 支間中央2~5 , 35.4 ,...
, 16.790, 24.3, 支間中央1~3 , 25.2, 支間中央1~4 , 25.9, ...
, 18.980, 19.0, 支間中央2~4 , 19.6, 支間中央2~5 , 20.2, ...
, 21.170, 15.2, 支間中央1~4 , 15.6, 支間中央1~5 , 16.0, ...
**桁高1000mm [#ibd647a4]
***鉛直材80mm・斜材80mm [#f8c140d5]
橋長20m/18mなどは計算時間がとても長い．要素数と要求次数の...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 27.3, 支点+鉛直材1 , 35.8, 支点+鉛直材1~3 , 38....
, 12.410, 22.6, 支点+鉛直材1 , 29.3, 支点+鉛直材1~3 , 31....
, 14.600, 19.3, 支点+鉛直材1 , 24.9, 支点+鉛直材1~3 , 26....
, 16.790, 16.8, 支点+鉛直材1 , 21.5, 支点+鉛直材1~3 , 22....
, 18.980, 14.9, 支点+鉛直材1 , 17.9, 支間中央1~3 , 18.8, ...
, 21.170, 13.3, 支点+鉛直材1 , 14.9, 支間中央1~4 , 15.4, ...
***鉛直材95mm・斜材95mm [#i2fb2ddc]
鉛直材80mm・斜材80mmと比較しやすいように最も太いものを計...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 31.2, 支点+鉛直材1 , 44.3, 支点+鉛直材1~2 , 60....
, 12.410, 25.9, 支点+鉛直材1 , 36.5, 支点+鉛直材1~2 , 42....
, 14.600, 22.1, 支点+鉛直材1 , 29.5, 支間中央1~3 , 31.1, ...
, 16.790, 19.3, 支点+鉛直材1 , 23.1, 支間中央1~4 , 24.3, ...
, 18.980, 17.1, 支点+鉛直材1 , 17.4, 支間中央1~3 , 18.4, ...
, 21.170, 14.4, 支間中央1~4 , 15.0, 支間中央1~5 , 15.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#i6f228f2]
3次モードぐらいまで確認したかったけど，もう期限的にその余...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 38.8, 鉛直材2~3 , 40.6, 鉛直材1~4 , ---- , ----...
, 12.410, 31.1, 鉛直材2~4 , ---- , ---- , ---- , ---- ,, ...
, 14.600, 26.0, 鉛直材2~4 , 27.1, 鉛直材1~4 , 27.5, 鉛直...
, 16.790, 22.2, 鉛直材2~4 , 23.1, 鉛直材2~5 , ---- , ----...
, 18.980, 18.8, 支間中央2~5 , ---- , ---- , ---- , ---- ,...
, 21.170, 15.0, 支間中央1~4 , ---- , ---- , ---- , ---- ,...
-ちなみに，橋長10mなどの短いものは，変位が最大となるとこ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#p9977548]
-橋長10mの1次モードは，鉛直材というよりも，支点付近の鋼板...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 62.1, 鉛直材1~6 , 64.6, 支間中央2~5 , 66.0, 鉛...
, 12.410, 43.5, 支間中央1~4 , 46.1, 支間中央1~5 , ---- , ...
, 14.600, 31.4, 支間中央2~5 , 33.2, 支間中央2~6 , 35.0, ...
, 16.790, 23.5, 支間中央1~4 , 24.6, 支間中央1~5 , 25.8, ...
, 18.980, 18.3, 支間中央2~5 , 19.2 , 支間中央2~6 , 20.0, ...
, 21.170, 14.6, 支間中央1~5 , 15.2, 支間中央1~6 , 15.8, ...
**グラフ [#r2857111]
手始めに構造研でお馴染みの3次元プロットでもつくってみる．...
***鉛直材80mm・斜材80mm [#g5209b20]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-発表や論文には使えないが，ブラウザ上でマウスで動かせたら...
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-もしかすると，色はz軸の値にするのではなくて，座屈箇所に...
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-なんだか，surfaceのx軸とy軸がよく分からないので，mesh3d...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材95mm・斜材95mm [#o746beae]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-ブラウザ上で表示
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-ブラウザ上で表示(座屈箇所ごとに色分け)
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材80mm・斜材80mm(対傾構2つ) [#fd860625]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材95mm・斜材95mm(対傾構2つ) [#yda7b568]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
*Salome-Meca2017の座屈解析メモ [#t6be1072]
**固有値が算定されるのに結果が見れない場合 [#z3abdde0]
固有値が近くて，誤差が許容範囲以下に収まらない時，message...
-座屈解析の CALC_MODES を編集 -> VERI_MODES にチェックを...
-ただし，計算されたモードは自己責任で評価して下さい，との...
-また，medファイルが重くて？モデリングや計算に使ったhdfフ...
**TRI_DIAG(Lanczos法)の精度を上げる方法 [#kfe64a70]
Code_Asterの座屈解析では，誤差(norme d'erruer)が0.000001...
-座屈解析の CALC_MODES を編集 -> SOLVEUR_MODAL の Edit を...
--COEF_DIM_ESPACE は，デフォルトでは4なので，それを適宜上...
---例えば求めたい次数が6のときに，COEF_DIM_ESPACEがデフォ...
--ただし，要求次数が6のとき，COEF_DIM_ESPACEをデフォルト...
--と思いきや，要求次数6のとき，試しにCOEF_DIM_ESPACEを12...
--さらに，COEF_DIM_ESPACEを24にしてみたら，計算時間がデフ...
--ということで，精度を上げるためにはCOEF_DIM_ESPACEを大き...
-PREC_LANCZOS(精度(precision)だと思う)を変えても，結果は...
-CalculiXでは Lanczos vectors はデフォルトで4$\times$次数...
--デフォルトで4$\times$次数とは！なかなか多いですね．しか...
--(
*BUCKLE [#e09c3fd1]
6, 0.01, 24, 1000
**次数，精度，Lanczos vectors，最大反復回数…この行はただ...
--)
**固有値解析 [#p2357bac]
-昨年度の方法を真似てSORENSEN法？とやらを使っているが，他...
-よくわからないが，アプローチ？を追加してみる．
--REELとIMAGで同じ結果になり，COMPLEXEは結果が出ない．ど...
***SORENSEN [#jdf3999b]
IRAM(Implicitly Restarted Arnoldi Method)というらしく(Sor...
-後藤さんいわく，何かの間違いで普通でない非対称の行列がで...
***JACOBI [#kcb8e997]
ただのjacobi法かと思いきや，公式ドキュメントにBatheとWils...
-後藤さんによると，時間がかかっているからといって必ずしも...
***TRI_DIAG [#vdb3c836]
名前はただの3重対角化っぽいけど，Lanczos法で実対称行列の...
-原則として，TRI_DIAGで計算することにしようかと思う．
-TRI_DIAGの精度(PREC_LANCZOS)はデフォルトでは1.0E-08だが...
--CaluliXとオンサイト木橋の座屈を比較してみたところ，座屈...
---ところで，CalculiXでは荷重を節点数で割っただけで，等価...
-試しにPREC_LANCZOSを1.0E-02に変えてみても，結果は何一つ...
-messageファイルにCOEF_DIM_ESPACEを上げてみては？と書かれ...
--COEF_DIM_ESPACEとは&link(ここ,http://femingway.com/?p=2...
-一応2・3次モードぐらいまで確認しておきたいので，正負のモ...
--デフォルトのCOEF_DIM_ESPACE:4から6に上げてみたところ，2...
***QZ [#q84ea5c5]
そのままQZ法．大規模行列には向いていない(自由度1000以下)...
**幾何剛性マトリックス [#f14a6fd4]
昨年度の方法では，幾何剛性マトリックスを作るために，"SIEF...
**載荷方向と座屈荷重の正負 [#u13f97ff]
なぜか逆向きになり，使用者がこんがらかる問題．
-片持ち梁のオイラー座屈を色々な組み合わせて計算してみたと...
-トラヘアの座屈を載荷方向を変えて計算したところ，載荷方向...
**CalculiXとかMarcとかのこと [#m806e0a9]
-CalculiXは固有値解析(というかほとんどの行列計算)をするの...
*オンサイト木橋の鋼板をSalome上で自動生成するPythonスクリ...
オンサイト木橋の鋼板について，図面のなかには情報が十分で...
ところで，後藤さんからオンサイト木橋の鋼板の形状を決める...
-&link(onsi_kou.py,http://www.str.ce.akita-u.ac.jp/~gotou...
--溝形のフランジ付きモデル．桁高・鉛直材幅・斜材幅・橋長...
--近い将来(少なくとも2020年頃までに)，Salome-MecaもPython...
--寸法(弦材幅や斜材の角度)は，Brepファイルにエクスポート ...
-&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~goto...
--座屈解析用に鋼板，木材，支点の線を作る．鋼板は&link(ons...
--次のように要素分割すれば，飛び出しが生じにくいというこ...
---鋼板 : NETGEN 1D-2D-3D -> Max.size : 20 , Min.size : 10
---木材(SubMesh) : NETGEN 1D-2D-3D -> Max.size : 60 , Min...
---境界面(SubMesh) : NETGEN 1D-2D -> Max.size : 10 , Min....
---Convert to Quadratic
-&link(onsiz.py,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
--(やり直す時間的余裕がなかったので)今回の計算には使って...
-&link(onsitz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
--&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
-&link(onsitz.py,http://www.str.ce.akita-u.ac.jp/~gotouha...
--&link(onsitz_org.py,http://www.str.ce.akita-u.ac.jp/~go...
-海老さんと川村さんによると，4面体に要素分割をするとき，M...
*ねじ山を作りたくない [#n60cf996]
菊地翔さんのPC鋼棒の緊張力ネタで「今後はねじ山を考慮した...
&link(McMaster(機械部品屋)のページ,https://www.mcmaster.c...
*軸方向に色々な荷重が作用した時の柱の座屈 [#fac7ce18]
**片持ち梁の先端に集中荷重が作用した場合（メモ） [#t20b94...
そこかしこに載っているもの．
-$P_{cr}=\frac{\pi^{2} EI}{4\ell^2}\approx \frac{2.467EI}...
**片持ち梁に等分布荷重が作用した場合（メモ） [#bfb3d57d]
ティモシェンコ『弾性安定の理論』の上巻から引用．
***Bessel関数を解いたら(p.97) [#kbea664c]
Bessel関数とやらの解法がよく分からない．
-$(q\ell)_{cr}=\frac{7.837EI}{\ell^2}$
***エネルギー法1(p.99) [#m421c3a4]
第一近似として，たわみを$y(x)=\delta (1-\cos \frac{\pi x}...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$(q\ell)_{cr}=\frac{7.89EI}{\ell^2}$
***エネルギー法2(p.101) [#fcb9d5b1]
さらに助変数を用いて，たわみを$y(x)=\delta_{1}(1-\cos \fr...
-$(q\ell)_{cr}=\frac{7.84EI}{\ell^2}$
**片持ち梁に等分布荷重（圧縮方向）と先端に集中荷重（圧縮...
等分布荷重qは，集中荷重Pと係数kを用いて，仮に$q=\frac{kP}...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$P_{cr}=\frac{(6\pi^5+(3\pi^5-12\pi^3)k)EI}{((-1536+432\...
--試しに，k=0を代入すると，先端に集中荷重が載った場合の座...
-ちなみに，等分布荷重が引張方向に作用していた場合，kの値...
-$P_{cr}=\frac{(6\pi^5-(3\pi^5-12\pi^3)k)EI}{((-1536+432\...
-計算が間違っている可能性も高いので，そのうちにFEMで確認...
**片持ち梁に三角形分布荷重（圧縮方向）と先端に集中荷重（...
図がないと載荷条件が分かりにくいが，描くのが面倒臭いので...
三角形分布荷重qは，集中荷重Pと係数kを用いて，仮に$q(x)=\f...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
-$P_{cr}=\frac{(30\pi^6+(5\pi^6-30\pi^4)k)EI}{((-4080+360...
**片持ち梁に三角形分布荷重（引張方向）と先端に集中荷重（...
図がないと載荷条件が分かりにくいが，描くのが面倒臭いので...
三角形分布荷重qは，集中荷重Pと係数kを用いて，仮に$q(x)=\f...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
-$P_{cr}=-\frac{\left(\left(5{{\pi}^{6}}-15{{\pi}^{4}}\ri...
**片持ち梁に台形分布荷重（圧縮方向）と先端に集中荷重（圧...
台形分布の下底（固定端側）を$\frac{k_{1}P}{\ell}$，下底（...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$P_{cr}=\frac{\left(30{{\pi}^{6}}+\left(5{{\pi}^{6}}-30{...
-maximaに貼り付ける用
--(
P1(k1,k2,l1):=((30*%pi^6+(5*%pi^6-30*%pi^4)*k1+(10*%pi^6-...
--)
--$k_{2}$に0あるいは$k_{1}$に0代入したら三角形分布と同じ...
--$k_{1}=k_{2}$にしたら等分布と同じになる．
**片持ち梁に台形分布荷重（引張方向）と先端に集中荷重（圧...
-↑の式をk1とk2をマイナスにしたら同じになる.
-$P_{cr}=\frac{\left(-30{{\pi}^{6}}-\left(5{{\pi}^{6}}-30...
-つまり，maximaでは，P1を貼り付けたあとに次のようにすれば...
--(
P2(k3,k4,l2):=P1(-k3,-k4,l2)
--)
**両端固定の問題に応用したい場合 [#dc6d4761]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-まず，上の片持ちのそれぞれの式（分布荷重が圧縮方向と引張...
-そして，Pcrをそれぞれ$Pcr_{1},Pcr_{2}$，スパン$\ell$をそ...
--連立方程式：$Pcr_{1}=Pcr_{2}, \ell_{1}+\ell_{2}=0.5$
---ここでは，具体的な値ではなく，便宜上$\ell_{1}+\ell_{2}...
-より，$\ell_{1}, \ell_{2}$を求め，元の式$Pcr_{1}$か$Pcr_...
***等分布載荷 [#d19735da]
kと$\ell_{2}$を$Pcr_{2}$に代入すればいい．$\ell_{2}$は，...
-かなり長いので，書きたくない．maximaに貼る用のみ．
--(
sol(k1,k2):=solve(P1(k1,k1,0.5-l2)=P2(k2,k2,l2),l2)
--)
--0.5は仮で，本当は$\frac{\ell}{2}$をいれる．別に$\ell_{2...
--等分布荷重が載っていない場合と同じになるかの確認のためk...
--(
bfloat(P2(k2,k2,sol(k1,k2)))
--)
--ここでは，kに具体的な数値を入れる．
--kの計算方法については考え中．多分，曲げモーメントによる...
***台形分布載荷 [#tde5cd2a]
kと$\ell_{2}$を$Pcr_{2}$に代入すればいい．$\ell_{2}$は，...
-かなり長いので，書きたくない．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...
--(
bfloat(P2(k3,k4,sol(k1,k2,k3,k4)))
--)
--ここでは，k1, k2, k3, k4に具体的な値を入れる．
--k1,k2,k3,k4の計算方法については考え中．
*直交異方性材料の直応力-直ひずみ関係 [#jfd24410]
前々から&link(構造関係メモ,http://www.str.ce.akita-u.ac.j...
-式(1)について、Maximaで求めた結果(&link(ihou1.wxm,http:/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
-確かに複雑……。
*木材利用2017メモ [#k31eb3da]
-アンカープレートの部分だけ木材と鋼板を一体化させる。
--アンカープレートの境界面でパーティションし、木材はLengt...
---これだと、鋼板の隣の非一体化部分が木材にすり抜けるので...
-堅固山さんとの違いは、
--木材を直交異方性にする。
---等方性と異方性にしたときの座屈荷重の比較をする?
--3点曲げに加えて積雪荷重での座屈も計算する。
--特に、支点付近の鉛直材に焦点を当てる。
---つまり、中央以外の非一体化部分もモデリングする。
--大池や虎毛などの大きいモデルでもやってみる。
---要素数の問題で部分的に取り出したモデルで計算するかも。
**5月初旬(発表概要)までに [#z3f1729b]
-7mモデルでやる。
--穴の拡大率は0,5,10,15,20,25,30%の7通りとし、7mモデルの...
--何もなしモデルと補剛材+対傾構ありモデルの2通りでやって...
-載荷条件は、3点曲げ(線載荷でいいかな?)と積雪荷重(積雪1m...
これだけで28通り。ここまでの計算自体はjobを実行すれば2～3...
-さらに木材を等方性と異方性にしたものとで比較する?ことに...
-中央と支点部の鋼板と木材だけ非一体化にしたものでもやって...
ひとまず、7mスパンについては全通りの計算が終わったが、結...
**後藤さんと相談した概要の流れをメモ [#f0170a09]
-1.端だけ非一体化・等方性で穴の大きさを変える(堅固山さん...
-2.しかし、木材を直交異方性にすると、それよりも座屈荷重が...
--不思議。それと、要素の跳び出しが等方性の時よりも増える...
-(3.)さらにアンカープレートのところだけ一体化にすると、座...
--しかも、木材と鋼板の境界に応力集中する箇所が増えるから...
-4.実際に影響が大きいと予想される積雪荷重が載った場合、3...
-対傾構と補剛材の影響をどこに挟むのがいいのかを訊くのを忘...
**6月から8月(発表スライドor論文)までに [#h1d08971]
-大池モデル(や虎毛モデル)でも同じことをやる。
-7mモデルで支点近くだけを取り出した一部モデルと従来の1/4...
こんなにやれるかな？
**おまけ [#afc06b56]
-5mスパン・板厚6mmの図面があったら(k2にはない?)、本当に70...
--座屈荷重は，部分的非一体化モデルで282kN, 全体的非一体化...
---というか，そもそも実験値（tugite）のはらみ（7mに3点曲...
-ちなみに，CalculiXで座屈の*STEP～*END STEPの前にプレスト...
-どうやら非線形座屈解析の出番はないかも。計算するにしても...
--幾何学非線形で計算すると，座屈荷重に至る前にまた別のあ...
-計算に使った積雪1mにつき$3.5kN/m^2$だが、昔の道路橋示方...
,降りたての雪, $150kg/m^3$
,やや落着いた雪, $300kg/m^3$
,圧縮された雪または多量に水を含んだ雪, $500～700kg/m^3$
*印刷用の折り紙円筒 [#s8a60d5d]
色々な折り目パターンに変えても、$\LaTeX$を経由してA4用紙...
ところで，7月10日の高校の先生との個別相談会用に印刷してみ...
-システム-->システム管理-->プリンター
-HL5270DNを右クリックして，プロパティ
-ジョブオプション-->「用紙に合わせる」にチェック-->適用
-&link(prjabara.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--蛇腹折り。例えば、&link(周方向6分割・高さ方向8分割の例,...
近藤夏実さんが折り畳まることを確認済み．
-&link(prkrasen.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--回転螺旋折り。例えば、&link(周方向5分割・高さ方向10分割...
近藤夏実さんに折って貰ったところ，あまりにも糊付けが難し...
-&link(prhrasen.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--反転螺旋折り。例えば、&link(周方向5分割・高さ方向6分割...
近藤夏実さんが折り畳まることを確認済み．こちらは糊付けが...
---螺旋折り系は、紙に目一杯描くようにしていないため、パタ...
-&link(prdaiya.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--ダイヤカット円筒。例えば、&link(周方向6分割・高さ方向8...
川村さんに折って貰ったところ，ちゃんと完成する．（これは...
*支間中央のみを取り出して線形座屈解析 [#pc2d0104]
局部座屈の着目点と比較して全体が大きい構造物だと、十分な...
手法は以下の通り。
-拘束は
--(z)橋軸方向対称面
--(x)幅員方向対称面
--(y)切断面の中立軸
-載荷は
--(y)輪荷重P/4
--(z)$\frac{曲げモーメント_{切断面}}{角材同士の中心間隔}$...
---ここで、角材と鋼板の剛性比($\frac{E_{鋼}I_{鋼}}{E_{鋼}...
-孔空き鋼板の3点曲げ(載荷は輪荷重)とした。
-モデルはk2に上がっていた堅固山さんのbrepファイルをそのま...
-堅固山さんが計算したものは、$P_{cr}$=270kNなので、1/4モ...
,, 1/4モデル(等方性) , 一部モデル(等方性) , 1/4モデル(異...
, 要素(C3D10)数 , 175493, 92862, 175493, 92862
, 1次モード , 270.0kN, 263.6kN, 228.5kN, 232.8kN
, 2次モード , 871.2kN, 851.3kN, 728.6kN, 745.7kN
, 3次モード , 1769kN, 1729kN, 1494kN, 1536kN
--ちなみに、応力集中する要素の飛び出しなどのノイズと思わ...
--モデルの体積比は$\frac{1/4モデル}{一部モデル}=5$なので...
-変位図(変位の倍率は見やすいように調整)
--そんなことをしていいのだか分からないけど、1/4モデルの1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ooi...
--相対誤差が大体2~3%ぐらいなので、なかなか良い精度だと思...
-分布荷重が載った場合についても、いずれやりたい。
**後藤ちゃちゃ(17/3/19) [#f902c5d0]
トラスの影響線で考えると、圧縮力が最も大きくなる鉛直材は...
一番近いところになるので、実は、支点側の座屈を解析する方...
その辺は、後でお話しましょう。
ところで、きのう、岩熊先生の退職パーティーに行って、鬆徒...
&link(3.5.2.4 積層板の見かけ上のYoung率,http://mechanics....
-ただの近藤メモ:堅固山さんに聞いてみたところ、純粋な鉛直...
--後藤ちゃちゃ(17/3/23)虎毛や大池のモデルは、鋼板上下が折...
*理工学デザイン [#k0eaa856]
--(
システムデザイン工学専攻 １２月２２日（木）１２時５０分～...
○会場：ＶＢＬ棟２階大セミナー室
○発表者は開始５分前までに準備を終えて会場に集合のこと。
○発表者は３名（教員２名、学生１名）への説明が終わった後も、
他の学生のポスターを見学したり、担当でない先生方へも対応...
終了時間まで会場に残ってください。
○当日、発表も質疑もない学生は出席しなくて構いません。
○ポスター発表当日の午前１１時からポスターの貼り付けなど準...
○ポスター用ボードにポスター発表番号が記載されていますので、
指定の場所にポスターを貼り付けてください。
○各自、ピンかセロテープを用意して、ポスター用ボードに貼り...
○終了後は、忘れずにポスターを持ち帰りください。
○Ａ１サイズ（縦８４１ｍｍ×横５９４ｍｍ）縦置きでポスター...
○ポスターの上段に、タイトル、所属（○○専攻○○コース）、
発表者氏名を明記してください。
○ポスターの作成は自由（手作りでも結構です）に行ってくださ...
情報統括センターを利用してポスターを作成する人は、
予めセンターに予約（前日17:00まで）してから行ってください。
なお、ポスター印刷する際は、利用申請書を必ず持参し...
○研究テーマや方向性、方法をポスターにまとめ、発表する。
指定された教員２名に説明する以外に、他専攻の学生に対して...
○各学生は、自身の属する専攻以外の学生のポスター発表を訪れ、
質疑応答し、その学生の研究内容を十分に理解した上で研究内...
1か月以内に理工（大学院）担当 に提出する。
特に、複数の質問事項を明確に記載し、質疑応答を記録する。
その質疑応答から理解した事をレポートにまとめる。
--)
-後藤ちゃちゃ(16/11/18)
各種の非線形問題（ねじり座屈とか接触とか）一通り解けてい...
（CalculiXも一通り解けそうですね）、線形でも回転螺旋折り...
解析が困難な各種の非線形問題を工夫しながら解いてますが、
オープンソースのツールだけで、割と何でも解けてます的なオ...
話にしてもいいかもしれません。まだM1なんで色々とおもしろ...
やってます。これからやっていきたい本命はこれなんですが、...
入れてもいいですが。
*PC鋼棒の引張+トルク+摩擦 [#i2a752fc]
-&link(toruku.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--PC鋼棒を引張った後、トルクを導入し、滑り出しをみる。
--PC鋼棒(梁要素)の先端に円盤(立体要素)をくっつけて、もう...
---菊地さんのと境界条件や荷重条件がやや異なるので、いずれ...
-&link(tth.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--gnuplot用にdatをまとめる。
*CLT関連で使えそうな板の文献をメモ [#i726a43d]
CLTのたわみ(実験値、FEMの値)が、どの理論式に近くなるのか...
-土木学会『構造力学公式集』@研究室
--相対2辺が単純支持された等方性長方形薄板：p.328
---部分分布荷重：p.330
---残りの相対2辺が弾性支持桁：p.331 .ここで、$\gamma=0$と...
--直交異方性板：p.367 .少しだけ長方形板の記述あり。ただし...
-ティモシェンコ『板とシェルの理論(上・下)』 @図書館
--直交異方性板は下巻に書かれている。これも具体例はあまり...
-渡辺昇『橋梁工学』 @研究室
--等方性板：p.300
--直交異方性板：p.324（そんなに詳しくない…）
-アムバルツミャン『異方弾性板の理論』 @図書館
--異方性の薄板、厚板、積層板について座屈、有限変位、振動...
-J.R. Vinson & R.L.Sierakowski『複合材料の構造力学』@図書館
--異方性の積層板について：pp.42-101
--さらに、7章に「強度と破損の理論」が載っている。これは面...
-宮入裕夫『複合材料入門 ー基礎と応用ー』@図書館
--異方性の積層板についてはあまり書かれていないのだが，何...
--ということで，ハニカムパネルについて読んでみると，面内...
---ハニカムパネルについては，図書館にはないが，宮入裕夫『...
*板の座屈 [#zaff6db6]
-メモ：4辺単純支持の座屈荷重
$N=k\frac{D\pi^{2}}{b^{2}} .$
ここで、
$k=(m\frac{b}{a}+\frac{1}{m}\frac{a}{b})^{2} , D=\frac{Et...
t:板厚 a:長さ b:幅 m:長軸方向の凹凸の数
と、色んな資料に書かれているが、Dは単位幅あたりの曲げ剛性...
-&link(itaz.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--4辺単純支持を板要素(S8)で計算。幅厚比b/t=100で試してみ...
-&link(hitaz.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--itaz.f90の結果で非線形座屈計算。4辺単純支持した板は、そ...
*プレストレス木箱桁橋 [#ce0220fb]
**孔無し一体化モデルの1/4解析 [#b0cb7c7c]
-&link(onsiyon.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--プレストレス木箱桁橋プログラムの練習として一番単純なモ...
---3点曲げか、4点曲げか選べるように改良(『木材の特徴を活...
--1/4解析にしたのは、計算が大変な接触解析とか局部座屈とか...
**孔無し非一体化モデルの1/4解析 [#me62708c]
いわゆる接触解析。
-&link(onsimyon.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--プレストレスのかけ方は、鋼板の上下の木材と隣り合わせて...
---載荷板に載荷するように改良。変位計の位置で変位を出力す...
-&link(onsimayon.f90,http://www.str.ce.akita-u.ac.jp/~got...
--アンカープレートのみにプレストレスをかけてみたが、プレ...
-せっかくなので、スパン7m($E_{zz}$=7.694GPa)の孔無し一体...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--とりあえず摩擦係数は中央値をとって、鋼材と木材:0.40 , ...
--軸方向は128分割で、上村さんの修論本体のグラフからすると...
--変位図(載荷方向)は、10倍に拡大している。
--実験値は上村さんの修論本体の表の実験2(k2にあるuemurahon...
--濡れ具合の他に、↓摩擦面の仕上げ方によっても、摩擦係数は...
---ttps://www.jstage.jst.go.jp/article/aijs/76/666/76_666...
***緊張力の低下 [#q027a29a]
木材のクリープやらPC鋼棒のリラクセーションやらで緊張力が...
-&link(onsirelax.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(onsimyon.f90,http://www.str.ce.akita-u.ac.jp/~got...
---載荷板に載荷するように改良。変位計で変位を読み取るよう...
-とりあえず、スパン6m($E_{zz}$=6.921GPa)のモデルに荷重70k...
--スパン7mモデル($E_{zz}$=7.694GPa)にして、計算し直そう。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--上図は緊張力が53kNの時の載荷方向の変異(1倍)図。落下防止...
--破壊試験(100kN)では、鉛直変位の最大値が15mmぐらいなので...
--グラフには鉛直変位20mm以上は載せていない。データが途中...
--ところで，色々とひっくり返すようなことをいうと，緊張力...
**孔無し一部非一体化モデルの1/4解析(局部座屈) [#t76f6c01]
中央のアンカープレート間の鋼板のみ非一体化にして、かつ接...
-&link(onsikyo.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--プレストレスを鋼板の上下木材部分に全面載荷。もちろん、...
-&link(onsikyoa.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--プレストレスをアンカープレートのみに載荷。こちらの方が...
-これらをスパン7m($E_{zz}$=7.694GPa)、鋼板の板厚を9mmとし...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
,, 線形座屈(輝さん) , プレストレスなし , 全面プレストレス...
, 座屈荷重(kN) , 339.1(78.5) , 326.6 , 座屈しない? , 338.8
--下の図は500kNの時の幅員方向変位の1倍。当たり前だが、そ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--&link(斉藤輝さんの線形座屈解析の結果,http://www.str.ce....
--ただし、鋼板の板厚方向の要素分割は、今は適当に4分割だが...
-理工学デザインの発表用に、線形座屈と非線形座屈を同じ要素...
***オイラー座屈から3点曲げ載荷時の座屈荷重を推定 [#zb7e6f...
-&link(onsioira_3ten.f90,http://www.str.ce.akita-u.ac.jp/...
--堅固山さんが(孔空き鋼板の)オイラー座屈を計算しているの...
, 鋼板上部の座屈(h=163) , 鋼板上部の座屈(h=120) , 鋼板上...
, 216kN, 192kN, 218kN,, 6201kN, 2777kN
---「異部分」は，木材を直交''異''方性材料とした，''部分''...
---「等部分」は，木材を''等''方性材料とした，''部分''的非...
--鋼板上部の座屈について，ここでは(1.直交異方性を考慮して...
--ということで，FEMで中央の鉛直材と斜材をなくして，木材を...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/kai...
, 鋼板上部の座屈(h=163) , FEM(等部分) , FEM(異部分) , FEM...
, 216kN, 262kN, 218kN, 209kN, 170kN, 192kN, 215kN
--となることから，座屈荷重は，異方性の影響が約2割↓，中央...
--特に，異方性の中でも，木材の軸直角方向のヤング率の影響...
---ちなみに，ヤング率を変える（軸直交方向を硬くしても柔ら...
---ただし，桁高800mmになると中央の鉛直材と斜材の補剛効果...
--3点曲げ載荷における鉛直材の座屈は$\frac{\ell}{4}$あたり...
***オイラー座屈から等分布載荷時の座屈する積雪深を推定 [#w...
-&link(onsioira_yuki.f90,http://www.str.ce.akita-u.ac.jp/...
--&link(onsioira_3ten.f90,http://www.str.ce.akita-u.ac.jp...
, 鋼板上部の座屈(h=163) , 鋼板上部の座屈(h=120) , 鋼板上...
, 20.5m, 18.2m, 20.6m,, 658.7m, 295.0m, 278.8m
--鉛直材の長さの決め方が難しく，それによって結果が2倍ぐら...
-対傾構をつけず，補剛材をつけたらどうなるかやってみたい．...
-Trahairでも$S(\frac{\ell}{2})$や$\tau=\frac{QS}{bI}$とか...
**datファイルのデータ整理 [#u09e325b]
-&link(plonsi.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--一体化モデルの1/4解析用。
-&link(plonsim.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--非一体化モデルの1/4解析用。
-&link(plonsirelax.f90,http://www.str.ce.akita-u.ac.jp/~g...
--緊張力の抜けの1/4解析用。
*弾性床上の座屈 [#wf37705a]
-&link(yuka.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--特に初期不整を与えなくても、断面が変形して弾性床と接触...
---ただし、"*SURFACE BEHAVIOR"の値をかなり上手く調整しな...
---やってみたのは単純梁と片持ち梁だが、どちらも床と離れて...
--仮に床方向に座屈させたとしても、鋼材と木材(繊維直角方向...
*鋼材と木材の接触解析 [#o862c564]
-&link(masatu.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--&link(JIS K7125,https://www.djklab.com/service/koubunsh...
--鋼材と木材の摩擦係数は0.2(湿)~0.6(乾)の中央値をとって0....
--ちなみに木材(上に載っている方)の繊維方向は滑らす方向と...
--"*SURFACE BEHAVIOR"が謎。
---&link(Abaqusの"*SURFACE BEHAVIOR"の説明,https://engr-a...
---"*SURFACE BEHAVIOER"を&link(あまりいじらない方がいいと...
---要は収束の「する」、「しない」に関わってくる。
--"*FRICTION"の右のstick range slopeの値が謎。
---&link(マニュアル,http://web.mit.edu/calculix_v2.7/Calc...
---ネットを漁っていると、食い込みのない理想的な値は∞だが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/mas...
--&link(pd.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/mas...
--縦軸が水平力、横軸が水平変位。変位は載荷面(横)の平均で...
*CLT [#n13ea13e]
**梁 [#l6cd48cb]
***梁(5層7プライ/7層7プライ)の3点曲げ [#u0bb2155]
-&link(clt.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似...
--ひき板の強軸方向ヤング率を試しに6GPaにして計算してみる...
***梁(5層7プライ/7層7プライ)の３点曲げを繰り返す(全体が弾...
-&link(cltso.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似...
--CalculiXの"*PLASTIC"はミーゼス応力での降伏値。木材は圧...
---"*PLASTIC"の行に続けて",HARDENING=～"を入れると、硬化...
--応力-ひずみ曲線(多直線近似)はかなり適当に入れてあるが、...
--"*PLASTIC"ではなく&link("*DEFORMATION PLASTICITY",http:...
-フィンガージョイントの降伏点が分かれば、そこだけ弾塑性材...
***h=240mm(7層8プライ)のときの梁理論との比較 [#dc468976]
-&link(clt3ten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(clt5ten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(pclt35ten.f90,http://www.str.ce.akita-u.ac.jp/~got...
--特に意味はないけど、&link(大黒屋さんの,http://library.j...
--桁高は重ね方の都合上、変えられないので、全方向の寸法を...
---つまり、b=24mm,h=240mm,$\ell$=240～4800mm(240mm刻みで...
--大黒屋さんに合わせて、分割数はnx=5,ny=24,nz=100.(大黒屋...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
--結果↑は大黒屋さんのものに似ている。
***b=24mm,h=240mm(7層8プライ)のときの1層と8層のFEMでの比...
-大黒屋さんの真似をしたついでに、FEMの相対誤差($\frac{v_{...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
--曲げが支配的になるほど、両者の値が近くなる。曲げ剛性か...
--弱軸の挙動もちょっと気になる。3点曲げと5点曲げで違う形...
**床版 [#pf637764]
''注意！''以下のプログラム類で1層化したモデルについては，...
$\nu_{zx}=\sqrt{0.0064\frac{E_{z}}{E_{x}}}=\sqrt{0.0064\f...
$\nu_{xy}=\nu_{xz}=0.0558, \nu_{zy}=\nu_{zx}=0.115$, $\nu...
に書き換えてみると，8層のものにかなり近い値がでるので，一...
***供試体2(7層8プライ/接続金物無し)の静的載荷 [#k6931ef5]
-&link(clt1.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--1層モデル。海老名さんの指摘により寸法や境界条件を直した...
---$\overline{E_z}=3.5$GPa, $\overline{E_x}=1.7$GPa, $E_y...
--- $\nu_{xy}=\nu_{xz}=0.4\cdot\frac{\overline{E_{x}}}{\o...
---$G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\ap...
--ちなみに、ヤング率の平均化は、$\frac{\sum_{i=1}^{8} E_{...
-&link(clt8.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--8層モデル。海老名さんの指摘により寸法や境界条件を直した...
---$E_z=5$GPa, $E_x=E_y=0.2$GPa.
---$\nu_{xy}=\nu_{xz}=\nu_{yx}=\nu_{yz}=\frac{0.4}{25}=0....
---$G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\ap...
-&link(cltiso.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--等方性材料。実験や上の2つとは全然合わないけど、構造力学...
***供試体2(7層8プライ/接続金物無し/ゴム有り)の静的載荷 [#...
-&link(clt1g.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
-&link(clt8g.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--ゴムのヤング率は分からないが、作るだけ作っておいた。ち...
--ちなみに、*HYPERELASTICでゴムのような超弾性も扱えるが、...
***供試体2(7層8プライ/接続金物無し)の弾塑性解析 [#k1f144ec]
-&link(clt1so.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
-&link(clt8so.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--中立軸の上下で別々の弾塑性材料にしただけ。応力ひずみ関...
---とはいえ、本とか色んなページ(&link(例えばこことか,http...
---しかも、応力-ひずみは、&link(この辺,http://www.str.ce....
-&link(clt1soko.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(clt8soko.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--中立軸の上下で別々の弾塑性材料にし、かつ支持部周辺を圧...
---藤田さんの計算結果が実験値よりも硬いのは、支持部周辺は...
-&link(plcltso.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--gnuplot用に荷重-変位関係をまとめる。
-ところで，入力データに特に書き込んでいないが，デフォルト...
--後で後藤さんから訊いた話では，apt-getでCalculiXをインス...
***供試体1(7層8プライ/接続金物有り/ゴム有り)の静的載荷 [#...
-&link(clt1jg.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
-&link(clt8jg.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--鋼材は板要素(S8)として木材と節点共有させたので、若干寸...
***供試体3(7層8プライ/供試体2の半分)の静的載荷 [#t2782c36]
-&link(clt1han.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(clt8han.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--図面をちゃんと見ていないので、寸法は正確でないと思うけ...
***$\ell$=2010mm,h=240mm(7層8プライ)のときの1層モデルと8...
-&link(clthaba.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--H鋼の間隔($\ell$=2010),床版厚(h=240mm(7層8プライ))を一...
--載荷は、床版中央に500x200mm(の1/4)の大きさを面載荷。
--境界条件は、面を(ボルト～H鋼の端)高さ方向と、横方向に固...
-&link(clthababun.f90,http://www.str.ce.akita-u.ac.jp/~go...
--床版の設計上で知っておきたいのは中央集中荷重よりも分布...
-&link(pclthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(clthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
---分布荷重についても同様に中央のたわみを見たが、分布荷重...
--たわみの収束性の確認をするのが面倒だったので、要素分割...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
-8層ものを1層のモデルに曲げ剛性から近似した場合、梁では両...
--あと、変位の相対誤差の大小に関わらず、1層モデルのミーゼ...
-それから、境界条件をCLT床版に近くして計算したが、より一...
--それなら、立体要素の8層モデルと1層モデル、板要素の1層モ...
***1層モデルと8層モデルのたわみの誤差についての考察 [#k00...
支間方向にz軸，橋軸方向にx軸，載荷方向にy軸をとった（梁の...
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{z...
この逆関係は，
$$\left( \begin{array}{c} \sigma_{x}\\ \sigma_{z} \end{ar...
となるので，&link(岩熊先生のページ,http://mechanics.civil...
$$ D_{x} \equiv \frac{E_{x}}{1-{\nu_{xz}}{\nu_{zx}}}I^{(p...
と書ける．ここで，$I^{(p)}$は単位幅の断面二次モーメント（...
例えば，上述したように挽き板一枚の材料定数
$$E_z=5GPa, E_x=E_y=0.2GPa, \nu_{xy}=\nu_{xz}=\nu_{yx}=\n...
から，8層分を平均化した材料定数を
$$\overline{E_z}=3.5GPa, \overline{E_x}=1.7GPa, E_y=0.2GP...
というように決めた時，8層モデルの板の曲げ剛性は（$\sum_{i...
$$ D_{8層x} = \frac{\overline{E_{x}I^{(p)}}}{1-0.016\cdot...
1層モデルの板の曲げ剛性は
$$ D_{1層x} = \frac{\overline{E_{x}I^{(p)}}}{1-0.19428571...
ここで，1層モデルと8層モデルとで梁としてEIが同じになるよ...
とはいえ，CLTにKirchhoff-Loveの仮定を適用するのにもかなり...
***1層モデルのポアソン比について考えてみる [#b53c4aa6]
というわけで，8層モデルのたわみに合わせようとするなら，ど...
軸のとり方やら材料定数やら寸法やらを上述した7層8プライの...
$$ \nu_{zx}=\sqrt{0.0064\frac{E_{z}}{E_{x}}}=\sqrt{0.0064...
で，あまり変位に影響しなさそうなポアソン比を，仮に今求め...
になる．こんな決め方でいいのかどうか，幅/スパン比を変えて...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/clt...
となり，以前のものと比べてかなり改善される．ちなみに今回...
-&link(clthabas.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(clthabab.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(pclthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
***長方形板の4辺支持の式を計算 [#z08442e4]
-&link(4hen.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--等方性長方形板の4辺単純支持。構造力学公式集より。
--2辺単純支持とかも作っていたが、たわみが発散する。式が複...
***供試体2(7層8プライ/ボルト有り)の1/4解析 [#b2859c74]
-&link(cltbyon.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--木材利用研究論文報告集15のpp.17-22の有山さんのを真似る...
--実験では降伏するまで載荷したものの、これといってCLTの損...
---計算してみたところ、随分とボルトの下部に応力が集中して...
--直方体(CLT)と円柱(ボルト)が組み合わさったメッシュを生成...
--メッシュをかなり細かくしないと、ボルトの断面が膨らむ(梁...
--静的試験について、エクセルのC列目には最大荷重が約10kNと...
**CLT橋 [#lb2cea18]
中村先生のCLT橋（桁に集成材，床版にCLTを使い，それらをラ...
-&link(nclt.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--3点曲げ載荷．単純に，CLTと集成材の境界面の節点を共有さ...
--境界の節点を共有させると，実験値よりもかなり剛性が高く...
-&link(ncltj.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--3点曲げ載荷．CLTと集成材の間に，接合要素みたいなものを...
--接合要素のせん断弾性係数の値によって，計算結果が大きく...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/ncl...
--上がCLT，下が集成材．150kN載荷時の変形を1倍で描くと，滑...
-&link(ncltj_bun.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(ncltj.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--A活荷重（13.5kN/m$^2$）が作用時のたわみは，中村先生の計...
---ちなみに，『木材利用研究発表会講演概要集16』のp.83の最...
-ncltjb.f90<--需要がなさそうなので，つくっていない
--ラグボルトを梁要素として，CLTと集成材の節点上に配置する...
*円筒 [#f07d6080]
**円筒(角筒) [#q6b64dd3]
円筒プログラムと板要素の練習。
-&link(entou.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--要素はS8.境界条件を色々と変えられるようにした。座屈モー...
-&link(hentou.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--初期不整を与えて、幾何学非線形解析をすると、途中で計算...
**折り目を剛結とした回転螺旋折り [#l2348b95]
-&link(gkaiten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--折り目は剛結。板一枚につき1要素なので、かなりメッシュが...
---&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--幾何学非線形解析にすると、高さ0.87hぐらいで計算が止まっ...
--段の高さに応じて高さ方向変位を多点拘束(*EQUATION)によっ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ent...
-&link(山内さんのモデルの作り方,http://www.str.ce.akita-u...
**折り目を剛結とした反転螺旋折り [#o3dd0190]
-&link(ghanten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--折り目は剛結。板一枚につき1要素なので、かなりメッシュが...
---&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--幾何学非線形解析にすると、高さ0.85hぐらいで計算が止まっ...
--上下を板に挟まれた中の段は均等に畳まるが、片方が自由に...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ent...
-線形解析をすると、回転螺旋折は上部が拡がってしまい、反転...
--&link(Salome-Mecaに詳しそうな人,http://salome-meca.coco...
--後に後藤さんから聞いた話によると、板要素の誤魔化してる(...
**板の要素数を増やす [#b84ef3d6]
-&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--螺旋折りプログラムでは板一枚につき要素数1つしかできない...
--多分、S6要素のみのモデルなら、何にでも使えると思う。た...
--元々は、&link(構造解析関係のプログラム置き場,http://www...
-折り目を重節点にして&link(多点拘束,http://www.str.ce.aki...
-立体要素の方が(遠藤さんの結果によると)幾何学非線形解析も...
*鋼棒のねじれ座屈 [#ef3ac8e2]
**C3D8とC3D6の混在 [#ff707b9d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(nezid6.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--ねじって線形座屈解析。軸方向に2000分割ぐらいしないと、...
--石坂さんの($\ell$=3m,d=13mm,E=206GPa,$\nu$=0.3)を計算す...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(hnezid6.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--例の如く、初期不整を与えてねじる。初期不整をとても小さ...
**C3D8だけ [#cf73940e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(nezid8.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--初期不正を色々と工夫していると、どうやら初期不正が大き...
-&link(hnezid8.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--nonpositive jacobian は出なくなったが、途中で計算が止ま...
-気が向いたら、C3D20でもやってみたい。
**B32 [#c799c384]
-&link(nezib32.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--梁要素。$T_{座屈}$=863.1N・mで、立体要素よりも理論値(86...
---後藤さんが『オイラー角を用いた簡潔な有限変位解析手法』...
-&link(hnezib32.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--例のごとく、非線形座屈解析。座屈する瞬間は確認できるが...
*横ねじれ座屈の非線形解析 [#c3399edd]
-後藤ちゃちゃ(16/5/17)
--FEMで初期不正のある座屈シミュレーションから$P-\delta^{2...
--あと、ローマ数字や◯数字の漢字コードは、昔は機種依存文字...
---とりあえず、半角数字にしました。
--平沢先生の図面をk2のmokusinにアップしました。
-&link(katasen.f,http://www.str.ce.akita-u.ac.jp/~gotou/p...
--汎用有限要素ツールは、固有値解析で座屈を一気に求めてい...
--かといって、複雑な形状とかを立体要素にモデル化して、真...
--そこで、プリポストプロセッサーと連携しやすく、立体要素...
-過去の横ねじれ座屈の卒論等を参考にした。
--寸法：b=2cm,h=25cm,$\ell$=2.0m.
---スパンは一つしか試していない
--分割：幅方向8,高さ方向30,軸方向200.
---過去の卒論等の"幅方向6,高さ方向20,軸方向100"だと、線形...
--材料定数
---等方性:E=206GPa,$\nu$=0.4 .
**手法 [#qf64594b]
-1.CalculiXで線形座屈解析をする。
--&link(yoko.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
-2.1/5000(例えば)の初期不整を与える。
--初期不整の大きさはそれなりの工夫が必要。
-3.CalculiXで幾何学非線形解析をする。
--&link(hyoko.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---gifのファイルサイズが大きいので一回消したのだが，遠藤...
---&link(cgxppm.f,http://www.str.ce.akita-u.ac.jp/~gotou/...
---1~5次モードでもいかにもそれらしい挙動を示す。異方性を...
-4.座屈荷重を&link(P-$\delta^{2}$法とか,http://www.str.ce...
--&link(pd2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--&link(sw.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--&link(meck.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--P-$\delta^2$法：終わりのの17~5点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
--Southwell法：初めの76点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---終わりの方のプロットも使うと、座屈荷重が少し大きくなっ...
--Meck法：初めの76点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---終わりの方のプロットも使うと、座屈荷重が少し大きくなっ...
,, 固有値 , $P-\delta^2$法 , Southwell法 , Meck法 , Timos...
,座屈荷重(kN), 46.8, 46.9, 47.3, 47.3, 45.3, 45.5
--&link(理論値,http://www.str.ce.akita-u.ac.jp/~gotou/zik...
---計算には&link(trahair.f90,http://www.str.ce.akita-u.ac...
-なお、結果は省くが、オイラー座屈なども同様なかんじ。ただ...
*伊藤さんのトラス橋 [#xd333fa9]
-金具詳細モデル(鋼材:E=206GPa,$\nu$=0.3)を2次要素にして、...
, C3D4(Code_Aster) , C3D10(Code_Aster&CalculiX) , 梁要素?...
, 0.239mm, 0.242mm, 0.228mm, 0.277mm-->0.275mm , 0.250mm
--要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v...
--異方性によるたわみの相対誤差(CalculiX(多分間違えている)...
---初めは斜め部材の角度を60度にして計算したが、ちゃんと61...
--異方性によるたわみの相対誤差(Code_Aster(あっていそう)):...
-ついでに、樹脂材料(E=2.84GPa,$\nu$=0.3)でもやってみた。
, C3D4(Code_Aster) , C3D10(Code_Aster&CalculiX) , C3D10直...
, 0.586mm, 0.669mm, 0.775mm-->0.773mm, 0.699mm
--要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v...
--異方性によるたわみの相対誤差(CalculiX(多分間違えている)...
--異方性によるたわみの相対誤差(Code_Aster(あっていそう)):...
-ちなみに、2次元のトラス(&link(toraITO.f90,http://www.str...
, 3次元トラス(Marc) , 2次元トラス , 2次元骨組 , 2次元骨組...
, 0.228mm, 0.22749mm, 0.22677mm, 0.22679mm
--Marcを使った3次元トラスの計算は伊藤さんによるものなので...
--金具詳細モデルの計算ほどせん断の影響が大きくない。不思...
-(2017/7)Code_Asterでも直交異方性材料として計算したところ...
--伊藤さんのトラス橋をやりたい人のために，斜めの直交異方...
-ただ，読み取っている箇所(中央の金具上の角材同士の中心線...
--確か，Marcのトラス要素と比較するために，そのような載荷...
*CalculiX関連 [#u3310209]
**unv-->inp [#g38f0dd0]
-&link(unvc3d10.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(unvc3d4.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--ついでに、節点群(*NSET)を座標指定でなくて、Salomeでつく...
-&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(unvc3d10.f90,http://www.str.ce.akita-u.ac.jp/~got...
--局所座標系の設定に不確かなところ(SYSTEM=CYLINDRICAL)が...
---jacobianエラー：輝さんの日誌によると、節点座標の出力形...
-&link(menc3d10.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~go...
---ccx_2.9だと問題ないが、ccx_2.6.1だとfrdファイルの要素...
---他にも、モデルによっては面情報を上手く変換できていない...
-&link(zac3d10.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~go...
-&link(s6.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
-&link(s3.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
--Salomeで切った3角形の板要素(unv)をCalculiX用入力データ(...
---面内曲げの変位にS3とS6とでさほど結果に差はないのだが，...
**線形座屈解析から初期たわみを与える [#j6172625]
-&link(syokitawami.f90,http://www.str.ce.akita-u.ac.jp/~g...
--「線形座屈解析用のinpファイル」と「全節点変位を出力した...
--何次の座屈モードであっても、どんな要素であっても、何種...
---ただ、inpファイルからコピー＆ペーストした部分の一列目...
-Code_Asterでも，Medでなく，IDEAS(unv)などのテキストファ...
**frd-->vtk [#r7d8378c]
-&link(c3d10vtk.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--c3d10のfrdファイルをparaview用にvtkファイルに変換する。...
--frdファイルには書かれていない主応力を計算してvtkファイ...
-&link(c3d8vtk.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--よく使うC3D8でも同じもの作る．
-s6vtk.f90を作ろうとしたのだが，S6はfrdファイルにおいて1...
--というか，vtkフォーマット以外を採用するのもいいかも．
-paraviewに計算結果を入れられると，こんな↓遊び？(CLTの最...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
*paraviewメモ [#wff29fdc]
主に個人の好みや拘わりで変更できるところを記す．
**コンター図関連(Color Map Editor) [#zd684856]
-ハートマーク(Choose preset)：初期設定の(青<-->赤)から虹...
--個人的には，デフォルトの Cool to Warm が大変見づらく、...
---&link(ここ,http://www.str.ce.akita-u.ac.jp/~gotou/prog...
-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 の数...
**あるグループのみ表示 [#cbdee9c6]
-Filter -> Alphabetical -> Extract Group -> 見たいグルー...
**csvに出力 [#z4be0f93]
Extract Selection をしないと，どこかの節点を飛ばされるら...
-(Filters)Slice(あるいは Clip ) -> Select Points On(ある...
--Select Points On : 表面の節点のみを選択
--Select Points Through : 内部の節点も選択
--Extract Group と組み合わせて(作ったグループだけ表示させ...
**動画化 [#g65c47a7]
-(File)Save Animation
--細かな設定は試していないから知らない．
-その後，ImageMagickのconvertを使う．
*Salome-Meca関連 [#s486f8e6]
**ある断面の軸方向応力分布を描画 [#m68c736b]
-&link(danmen.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--paravisから出力したデータをgnuplot用に整える。ついでに...
--ちなみに、断面をグルーピングし忘れたときでも、以下の方...
---パイプブラウザーのrmedファイルを選択したまま、フィルタ...
---見たい断面の、原点座標(origin)を入力し、面の法線(nomal...
-&link(sankaku.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--$\sigma=0$の線とか$\sigma=\frac{M}{I}y$の線とかをgnuplo...
-&link(szz.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--'MED'ではなく'RESULTAT'に出力したデータを、gnuplot用に...
*Maximaリンク [#k95ca0eb]
とても便利な数式処理のソフトウェア．terminal上で動くmaxim...
-斉木先生のページ
--ttp://mechanics.civil.tohoku.ac.jp/saiki/strcanalysis/k...
*ImageMagickメモ [#o77bc022]
Windowsの場合は、多分&link(ImageMagickの公式ページ,https:...
**複数のスクリーンショットを一括でppm(アスキー)に変換 [#t...
--(
mogrify -format ppm -compress none *png
--)
-&link(cgxppm.f,http://www.str.ce.akita-u.ac.jp/~gotou/pr...
**定点観測で撮影した連続写真を重ね合わせ [#nedf90c1]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
--(
convert *.JPG -evaluate-sequence max gousei.png
--)
-場合によっては、コマンド中のmaxをminにした方がいい。
-特に意味はないが、上のコマンドではデジカメの画像形式から...
-この画像をImageJ(Windowsの場合は&link(公式ページ,https:/...
--変位を測りたい場合は、
---Straght:距離が既知な線に沿って、線を引く。(ドラッグ&ド...
---Analyze-->Measure:線の長さをpixelで測る。
---Analyze-->Set Scale-->Known distance:既知な線の実際の...
---Analyze-->Measure:線の長さを実際の長さで測る。
--たわみ角(傾き)を測りたい場合は、
---Angle tool:測りたい角度をなす3点を順にクリック
---Analyze-->Measure
-よく使いそうな保存の形式は、
--File-->Save as-->Results:結果をテキストで出力。
--File-->Save as-->XY Coordinates:現在、表示されている線...
-野田さんと3年生のたわみ試験で試してみたところ、精度がと...
-原因としては、
--1.三脚のカメラを固定するところが壊れているので、シャッ...
---これは、ちゃんとした三脚を使えばいいだけの話。
--2.被写体と1m近く離したのだが、たわみが小さい(8mm程度)こ...
--3.なるべく被写体の真正面にカメラを据え付けたつもりだが...
---原因3と原因2は、どちらかを優先させると、もう片方が疎か...
-同時に&link(Regard3D,http://www.regard3d.org/)というWind...
--追記：Ubuntu18.04LTSでは&link(COLMAP,https://colmap.git...
--特徴点がたくさん抽出できる場合には、 1.Regard3D(SfM + ...
---ネット上にあったカエルのぬいぐるみや猿のぬいぐるみ(&li...
--ただし、百目石橋でもそれをやろうとすると、相当いいPC(ゲ...
-野田さんによると，tiffだと拡大してもぼやけにくいとか．例...
-convert -threshold 50% -edge 1 hoge.JPEG hoge.tif
-などとすると，輝度50%を境(50%でなくてもいい)に2値化し，1...
**pdfを変換するときに画質を下げない [#u11b1304]
pdfを画像に変換したい時，ただ，"convert hoge.pdf hoge.png...
*$\LaTeX$メモ [#v3c82d4c]
**minipageのキャプション [#ff759cff]
画像をminipageで横並びさせると，互いのキャプションが窮屈...
**単位 [#p7961415]
無頓着なもので，単位が斜体だとか間に半角スペースが入って...
--(
$\frac{P}{A} = \SI{100}{N/mm^2}$
--)
texlive-scienceは，"sudo apt-get install texlive-science"...
他にも色々と便利なマクロが梱包されていると思う．
**ダブルクォーテーション [#h41d44c1]
ずっと知らなかったが，"hoge"だと残念なことになる，正しくは
--(
``hoge''
--)
**画像をわざと表示させない [#s5197781]
最終版まで画像は表示される必要はなく，画像の範囲だけ調整...
--(
\documentclass[a4j,12pt,dvipdfmx,draft]{jarticle}
--)
というように，draftを入れてやると，画像は表示されず，画像...
**$\LaTeX$周りのツール群 [#tfc9f21c]
身の周りの人があまり使っていない？けど，とても便利なツー...
***pdfcrop [#o9674473]
TeXliveに勝手についてくるtexlive-extra-utilsに入っているp...
-余白無し
--(
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-cro...
***Ghostscriptによるpdfのグレースケール化 [#u5e5e5a5]
白黒で印刷をする前に，カラーで作ったpdfがグレースケール化...
"convert -type grayscale -density 300 hoge.pdf hoge_g.pdf...
--(
gs -sOutputFile=hoge_g.pdf -sDEVICE=pdfwrite -sColorConve...
--)
-ここで，1つ目のオプションの"hoge_g.pdf"は出力ファイルな...
-また，"?"はwikiの機能で勝手に付いてしまっているだけなの...
-他の細かなオプションについては分かりません．
-参考URL：
--http://takuya-1st.hatenablog.jp/entry/2014/09/22/010805
***Inkscape [#nf9f01af]
ドロー系の描画ソフト。Xfigの上位互換でgnuplot→Inkscape→$\...
コマンドラインで形式変換が行えるので、シェルスクリプトに...
**おまけ（修論の背表紙） [#z0b8365b]
修論本体をpdfで提出するようになったものの，依然として修論...
-&link(se.tex,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
--用紙の向きが縦書き．余白や空白を適宜調整する必要あり．
-&link(sel.tex,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--用紙の向きが横向き．タイトルが短いため横向きにしてみた...
---&link(例,http://www.str.ce.akita-u.ac.jp/~gotouhan/j20...
ちなみにfirefoxだと，縦書きの「ー」は横向きのままのようだ...
*その他，便利なコマンドライン [#ta1d4bec]
**pdfimages [#x7d58f6a]
pdfから画像を抜き取るのに使える．例えば，pngとして画像を...
--(
pdfimages -png hoge.pdf out
--)
**grep [#p3818bef]
***特定の文字を含むファイルを検出 [#dc7039b5]
--(
grep ほげ -rl .
--)
これだと今作業しているところ(.)にあるファイルで「ほげ」と...
勿論，カレントディレクトリ以外のパスでもいいし，オプショ...
*その他、便利なソフトウェア [#a143fc30]
**[[Hugin:http://hugin.sourceforge.net/]] [#aa12ccdf]
パノラマ写真をつくるためのソフトウェア。aptでインストール...
**[[Resynthesizer:http://www.logarithmic.net/pfh/resynthe...
GIMPでスマート消去するためのプラグイン。
-インストール：https://github.com/bootchk/resynthesizer
-使い方：自由選択 → フィルター → 強調 → スマート消去 → 半...
*3Dプリンター メモ [#f2cd76b0]
-ボールジョイントも作れそう(PRN3DでなくてEDENでの話)。
--http://rayerd.sblo.jp/article/172174211.html
---造形後に関節を折り曲げることで、隙間につまったサポート...
--http://rayerd.sblo.jp/article/172249880.html
--http://rayerd.sblo.jp/article/172529855.html
*gnuplotメモ [#id35cce2]
$\LaTeX$でグラフを描く時は，gnuplotとxfigを使うのよりも，...
**日本語のpngを出力する [#x8e397c5]
ロードするファイルの頭（set out "hoge.png"の上）に
--(
set term png font "/usr/share/fonts/truetype/takao-gothic...
--)
を入れる（sizeは好きな値にする．特に指定しない場合は入れ...
--(
nkf -guess hoge.plt
--)
でShift-JISですよ，などといわれた場合には
--(
nkf -w-overwrite hoge.plt
--)
でutf-8に変換する必要がある．
*Fortran90(gfortran)メモ [#ee62b9d7]
**固有値計算 [#z320f171]
-&link(jacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--ヤコビ法によって実対称行列の固有値と固有ベクトルを計算...
---frdファイルをvtkファイルに変換しようとして，frdファイ...
---scipyやoctave,julia等を使えば一行で簡単に固有値と固有...
-&link(sjacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--&link(jacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
**LAPACK [#i3976944]
学習・自作が難しそうな大規模な行列演算が簡単にできる．こ...
-Ubuntu16.04LTSでのインストール方法
--https://qiita.com/AnchorBlues/items/69c1744de818b5e045ab
--全くこの通りで整備できた(LAPACKとLAPACK95)．
-コンパイルの仕方やお試し計算
--http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux...
**FFTW [#a97b2f39]
高速フーリエ変換などが簡単にできる．
野田班でフリーのfftアナライザを使っていたが，自分で作って...
&link(梅原さん,http://www.str.ce.akita-u.ac.jp/~gotou/pro...
-Ubuntu16.04LTSでのインストール方法
--http://moznion.hatenadiary.com/entry/20120718/1342626456
--全くこの通りで整備できた(リンク先はUbuntu12.04だけど)．
-コンパイルの仕方やお試し計算
--http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux...
*Julia メモ [#adc05a47]
マトリクス構造での導入案があるとかないとか？
インタープリタ言語の割に実行速度がとても早く，;でシェルモ...
強いて言えば，（グラフ描きツールなど）色々と選択肢があり...
**Ubuntu上で最新版にする [#j30a97d0]
--(
sudo add-apt-repository ppa:staticfloat/juliareleases
sudo add-apt-repository ppa:staticfloat/julia-deps
sudo apt-get update
sudo apt-get install julia
--)
既にJuliaをインストール済みの場合は，
--(
sudo apt-get autoremove julia
--)
としてから上記の手順を踏まないと，駄目でした．
**パッケージのインストール [#ced80124]
例えば，グラフ描画ツールなら，juliaの対話モードを立ち上げ...
--(
Pkg.add("Plots")
--)
と打つだけでインストールできる．そして，使うときは using ...
-ちなみに，ここでの例のPlotsは何だかplotly.jsと同じみたい...
-JuliaはPGFPlotsにも対応しているようなので，やはり論文な...
**JuliaFEM [#aaa26ba5]
こんなのもある。
-http://www.juliafem.org/
*Python メモ [#sa31cbf5]
**コメント文にutf-8の日本語を使えるようにする [#t225c8d5]
ヘッダーにこれをいれる。ただし、Python3系では入れなくてよ...
--(
# -*- coding: utf-8 -*-
--)
**Matplotlibで日本語を使えるようにする [#jf2ccabd]
端末にて、
--(
mv ~/.cache/matplotlib/fontList.cache ~/.cache/matplotlib...
sudo cp /usr/share/fonts/truetype/takao-gothic/TakaoPGoth...
sudo vi /etc/matplotlibrc
--)
と入力した後、matplotlibrc内の"#font.family : sans-serif"...
*Rustメモ [#sfcec57b]
初めて聞く概念ばかりで難しい…．そろそろLinuxもC以外の言語...
**rust.vim [#y876ed73]
まだ若い言語なので（Ubuntu16.04での）デフォルトのvimでシ...
-https://github.com/rust-lang/rust.vim/blob/master/syntax...
-を"/usr/share/vim/vim74/syntax/"にダウンロードする．
-次に"sudo vi /usr/share/vim/vim74/filetype.vim"で，どこ...
--(
" Rust
au BufNewFile,BufRead *.rs setf rust
--)
-を追加し，拡張子が被っているHerculesをいじってごまかす(...
--(
" Hercules
" au BufNewFile,BufRead *.vc,*.ev,*.rs,*.sum,*.errsum s...
au BufNewFile,BufRead *.vc,*.ev,*.sum,*.errsum setf herc...
--)
-どうやらvimには拡張子(言語)別に辞書機能のようなものや自...
*Goメモ [#k38b93fc]
シンプルで早い．試していないけど，クロスコンパイルできる...
**パッケージのインストール [#b7e8fd15]
標準ライブラリだけでも物凄い数があるが，
それ以外にもオープンソースのライブラリが沢山公開されてい...
それらを入手するときの注意点として，golangだけでなく''git...
--(
go get github.com/〇〇
--)
-"github.com/○○"はURL．
-数値計算系ライブラリの有名どころを探してみると，それなり...
--BLAS/LAPACK:github.com/gonum/gonum
--FFTW:github.com/runningwild/go-fftw
-ただし，これは&link(この方法,http://golang-jp.org/doc/in...
**フォーマッタ [#qb158e8f]
インデントやスペース数が揃ってなくてもコマンド一発で綺麗...
--(
go fmt hoge.go
--)
*Haskell(GHC)メモ [#y3b3b7bc]
純粋関数型言語とやらを触ってみたい．
**パッケージのインストール [#abab45ce]
例えば，行列演算ライブラリとかは，
--(
cabal update
cabal install hmatrix
--)
gnuplotも&link(呼び出せる,http://d.hatena.ne.jp/mmitou/20...
*METAFONTメモ [#c59d00e6]
同研究室の某学生がM冨フォントを使いたいと話していたので，...
モデルとなる字の入手が困難を極めるため，道のりは険しい……．
**gftodvi [#efee22ee]
文字の確認方法．
-mf hoge.mf と打つと，hoge.2602gf というファイルができる．
--METAFONTについては，クヌース先生が書いた本（日本語訳）...
--他の参考ページ：ttp://hayato1224.web.fc2.com/MFmain.pdf
-gfファイルをdviファイルに変換するために，gftodviというコ...
--(
gftodvi: fatal: tfm file `gray.tfm' not found.
--)
-というエラーが出てしまい，生成されたdviファイルは0バイト...
--(
mktextfm gray
--)
-と打っておくと，gftodvi hoge.2602gf でgfファイルをdviフ...
-あとは，xdvi hoge.dvi& でも開けるし，dvipdfmx hoge.dvi ...
**mktextfm [#ye59d54e]
tfmという形式にすると自作フォントを$\LaTeX$で使うことがで...
--(
mktextfm hoge
--)
しかし，普通の文字や数式の感覚で使えない？ので，使い勝手...
*Vimメモ [#h3f9ed74]
**ある単語の数を数える [#o7aeff5a]
--(
:%s/hoge/&/gn
--)
**シンタックスハイライトの追加 [#l9a91177]
デフォルトでかなり沢山の言語が，勝手にシンタックスハイラ...
-"/usr/share/vim/vim74/syntax"にGitHubなどからhoge.vimを...
-"/usr/share/vim/vim74/filetype.vim"を編集し，拡張子から...
--もしかすると，同じ拡張子を使っている言語もあるかもしれ...
*USBメモリにバックアップ [#f293b776]
--(
rsync -auv /home/kouzou/2016 /media/disk
--)
-オプション"-auv"は，「a:パーミッションを維持，u:差分のみ...
-"disk"は、USBメモリに付けた名前。
-ちなみに，k2へのバックアップ時にも使える(パスを間違えそ...
*HHKB Pro2 [#sb4fe682]
HHKB Pro2 をただPCに繋いだだけだと，\や_が打てなくて困る...
&link(この配列,https://www.pfu.fujitsu.com/hhkeyboard/lea...
-システム -> 設定 -> ハードウェア -> キーボード -> レイア...
-右上のキーボードマークを右クリック -> 設定 -> + -> US ->...
間違いなく打ちやすいが，日本語キーボードを普段使っている...
とはいえ，じっと眺めているとUS配列の方が記号の位置がきれ...
*講義関連 [#edbe1bee]
**LibreOffice [#e18c2630]
-https://help.libreoffice.org/Main_Page/ja
**構造力学特論 [#befb0708]
-&link(torasu2.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(hone2.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--格点を剛結にするとどのくらい曲げが入るのか気になったの...
-&link(toraITO.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(honeITO.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--伊藤さんの木製トラス橋を２次元の問題として解いてみたか...
-&link(honetimo.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--さらに、伊藤さんの木製トラス橋の2次元骨組を剛性行列をテ...
**土質工学特論 [#f7040b16]
-間隙水圧分布を計算:&link(u.f90,http://www.str.ce.akita-u...
--set zrange [] reverse : z軸を下向き正にする。
--set pm3d : z軸の値を色付けする。x軸でこれをやりたい。調...
*TA関連 [#v35f5e11]
**基礎物理学実験 [#ua539bd7]
-&link(指導書,http://sht.phys.akita-u.ac.jp/~buturi/)
-&link(基礎物理学実験TA研修システム,http://sht.phys.akita...
--Ubuntu12.04のFirefoxでWindows Media Playerを再生するに...
---gnome-mplayer
---gecko-mediaplayerSynaptic
--をパッケージマネージャーでインストールすると、エラーは...
**土木環境工学実験(アクリルの引張試験) [#aff52f8f]
-ヤング率とポアソン比の確認:&link(yangu.f90,http://www.st...
--syokignu.fでデータ整形後
*時間割と単位 [#i354bb5c]
, 授業科目区分 , 修了に必要な単位 , 取得単位数 , 2017年度...
,,,,,
, 共通科目(必修) , 2単位 , 2単位,,
, 専門科目(必修) , 13単位 , 1単位, 12単位,
, 共通科目および専門科目(選択) , 10単位以上 , 12単位,, (...
, 下の3つを合わせて15単位以上, 2単位以上 , 2単位,, (他専...
,,, 4単位,, (他共通)
-「関連性のあるコースで開講する専門科目」とは：理工学研究...
--しかし，2017年度入学の人は，ガイダンスで，システムデザ...
-2017年度に取得予定の残り12単位の内訳は、
--システムデザイン工学課題研究：専門必修(2年間)=10単位
--システムデザイン工学演習：専門必修(通年)=2単位
**2016年度 [#m81db7e0]
-理工学デザイン：1共通必修(後期前半)
-システムデザイン工学課題研究：5専門必修(通年)
***2016年度前期 [#y7eed7a5]
,, 月 , 火 , 水 , 木 , 金
,,,,,,
, 10:20 , 構造力学特論 , 都市システム計画特論 ,,, 材料設...
, ｜ , 218 , 総合研究棟7階 ,,, c317
, 12:00 , 2専門選択 , 2専門選択 ,,, 2専門選択
,,,,,,
, 12:50 , , 土質工学特論 , 基礎物理学実験B , 土木環境工学...
, ｜ ,, c317 , 教育4号館104 , 218 ,
, 14:20 ,, 2専門選択 , (TA) , (TA) ,
,,,,,,
, 14:30 , ゼミ ,, 基礎物理学実験B , 英語ゼミ , 地域防災学...
, ｜ , 218 ,, 教育4号館104 , 218 , 共320
, 16:00 , (おまけの人) ,, (TA) , (おまけの人) , 2共通選択
***2016年度後期 [#n42ca291]
,, 月 , 火 , 水 , 木 , 金
, 12:50 , 基礎物理学実験D , 構造設計学特論 , 木構造工学 ,...
, ｜ , 教育4号館104 , c317 , c317 , c319 , 218
, 14:20 , (TA) , 2専門選択 , 2専門選択 , 1専門必修 , (お...
,,,,,,
, 14:30 , 基礎物理学実験D ,,地域産業アントレプレナー論(前...
, ｜ , 教育4号館104 ,, 共127 ,, 共224
, 16:00 , (TA) ,, 1共通必修 ,, 2共通選択
,,,,,,
, 16:10 , 地域産業論 ,,,,
, ｜ , 2号館3階学習支援室 ,,,,
, 17:40 , 2他専門選択 ,,,,
終了行:
#contents
$雪荷重を受けるプレストレス木箱桁橋の座屈挙動$
$プレストレス木箱桁橋鋼板部の局部座屈挙動（仮）$
-桁高やスパンが変わったときに，座屈挙動がどう変わるか？
M1のとき:$汎用有限要素ツールを利用した非線形解析（案）$
*ちょっと気になったこと（2018/02/15．もう，遅い…） [#j306...
対称条件から1/4モデルで座屈解析を行ったが，「支点付近の座...
しかも，座屈が幅員方向に非対称なモードとなる場合は，対傾...
ただ，十分な精度が得られる要素分割で1/2解析なんて，対傾構...
小さいモデルで確認しておきたい．
*フランジ付きモデルの材料定数の影響 [#d3d3211d]
一体化方法の影響について計算してしまったので，木材を等方...
-対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, ...
, 木材の材料 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 直交異方性 , 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 2...
, 等方性 , 27.2, 支間中央1・2 , 27.7, 支間中央1~3 , 28.5 ...
-対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm,...
, 木材の材料 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 直交異方性 , 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26...
, 等方性 , 27.5, 支間中央1・2 , 27.8, 支間中央1~3 , 28.5 ...
-やはり，異方性の影響は1割ぐらいある．木材を等方性材料に...
*フランジ付きモデルの一体化方法の影響 [#b6442f90]
木材利用2017で，部分的非一体化モデルがいいと言っておきな...
-対傾構なしの大池モデル(橋長16.790, 桁高800mm, 斜材85mm, ...
, 非一体化方法 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) ,...
, 全体的 , 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 25.9 ...
, 部分的 , 25.4, 支間中央1 , 27.5, 支点付近 , 30.7 , 支間...
-対傾構2枚(端部と中央)の大池モデル(橋長16.790, 桁高800mm,...
, 非一体化方法 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) ,...
, 全体的 , 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 ,...
, 部分的 , 26.5, 支間中央1 , 31.7, 支間中央1 , --- , --- ...
-一体化方法によって座屈モード(PC鋼棒間をまたぐかどうか)が...
--フランジなしモデルでは，一体化方法による影響が1割もあっ...
-フランジ付きモデルだと，どう一体化しても，結局すり抜ける．
--フランジなしモデルの座屈モードは，部分的非一体化モデル...
-実は部分的非一体化モデルの方が木材と鋼板の接している面積...
-全体的非一体化モデルにしても，飛び出しが生じないモデル化...
--木材利用2017当時は，全体的非一体化モデルにすると飛び出...
*対傾構の影響 [#n60c922d]
-大池のモデルで，対傾構の数を変化させる．
--1.対傾構なし(0)
--2.端部・支間中央(2)
--3.さらに間に1枚入れる(2+1=3つ)
--4.さらに間に1枚入れる(3+2=5つ)
-対傾構を入れると，（境界の木材をSub Meshで分割したせいも...
--モデル化を工夫する(&link(onsitz.py,http://www.str.ce.ak...
-寸法は大池モデルにし，1/4解析を行う．
--たわみ，軸方向応力について共に最大値のみ参考値としてメ...
, 対傾構の数 , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , ...
, 0, 24.5, 支間中央1・2 , 25.2, 支間中央1~3 , 25.9 , 支間...
, 2, 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 , 支間...
, 3, 25.2, 支間中央1~3 , 25.4, 支間中央1~3 , 26.1 , 支間...
, 5, 25.3, 支間中央1~2 , 25.9, 支間中央1~3 , 27.5 , 支間...
-[考察]
--対傾構が直接効いてくるのは支点付近の座屈だと思うが，そ...
--しかし，対傾構をまたいで座屈するモードは算出されていな...
--対傾構があると断面変形が抑えられるため，最大たわみは当...
*まとめ(橋長・桁高・孔を変化させた時)(対傾構を2つを入れて...
3.5kN/m$^2$の等分布荷重載荷時の座屈解析の結果．
-とりあえず，橋長については，約10mから約21mmまで(虎毛と大...
-とりあえず，桁高についても，500mmから800mmまで(虎毛と大...
--鉛直材の挙動を見るために，1000mmまで増やしてもいいのだ...
-鉛直材・垂直材の幅は80・85・90・95mmの組み合わせが16通り...
--ということで，一番きついであろう80mmx80mmを始めにやって...
--80mmx80mmとの比較のために，一番変化が分かりやすい95mmx9...
--CalculiXみたいに夜通し連続して計算できれば，パターン数...
---ASTKがSalome-Meca2017でなくなってしまったかと思ってい...
-表中の支間中央No.のNo.は，PC鋼棒間を支間中央から順に数え...
-表中の鉛直材No.のNo.は，鉛直材を端部から順に数えたもの．
-たわみと軸方向応力について，せっかく計算したので何かの足...
--当たり前だが，橋長が短いほど断面変形が大きく，最大のた...
-固有値が近接しているため，何次モードまで計算するべきか迷...
**桁高500mm [#jc7d20ca]
***鉛直材80mm・斜材80mm [#j968097b]
-橋長10.220mで，いくつか飛び出しモードが出てしまった．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 58.4, 支間中央1~3 , 59.0, 支間中央1~3 , ---- , ...
, 12.410, 39.2, 支間中央1~3 , 40.3, 支間中央1~3 , 41.0, ...
, 14.600, 28.4, 支間中央1~4 , 28.9, 支間中央1~4 , 29.5, ...
, 16.790, 21.4, 支間中央1~3 , 21.8, 支間中央1~4 , 22.4, ...
, 18.980, 16.7, 支間中央1~3 , 17.1, 支間中央1~5 , 17.2, ...
, 21.170, 13.4, 支間中央1~3 , 13.6, 支間中央1~4 , 13.9, ...
***鉛直材95mm・斜材95mm [#uf3b6c41]
-橋長10.220mで，飛び出しモードばかりが出てしまい，なかな...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 60.9, 支間中央1~2 , 61.4, 支間中央1~3 , ---- , ...
, 12.410, 41.3, 支間中央1~3 , 41.5, 支間中央1~3 , 42.5, ...
, 14.600, 29.8, 支間中央1~4 , 29.9, 支間中央1~4 , 30.8, ...
, 16.790, 22.5, 支間中央1~3 , 22.7, 支間中央1~4 , 23.1, ...
, 18.980, 17.5, 支間中央1~4 , 17.6, 支間中央1~4 , 17.9, ...
, 21.170, 14.1, 支間中央1~4 , 14.2, 支間中央1~4 , 14.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#t8bc4aab]
-やはり，木材と鋼板の境目で木材側の要素のアスペクト比を小...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 60.0, 支間中央1~4 , 62.1, 支間中央1~4 , 63.5 , ...
, 12.410, 39.7, 支間中央1~4 , 40.7, 支間中央1~3 , 42.3, ...
, 14.600, 28.6, 支間中央1~4 , 29.5, 支間中央1~5 , ---- , ...
, 16.790, 21.5, 支間中央1~4 , 22.0, 支間中央1~4 , 22.5, ...
, 18.980, 16.8, 支間中央1~5 , 17.2, 支間中央1~6 , 17.4, ...
, 21.170, 13.4, 支間中央1~4 , 13.7, 支間中央1~5 , 13.9, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#i74314c3]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 61.9, 支間中央1~3 , 64.2, 支間中央1~4 , 66.0, ...
, 12.410, 41.7, 支間中央1~4 , 42.5, 支間中央1~4 , 43.4, ...
, 14.600, 30.0, 支間中央1~4 , 30.8, 支間中央1~5 , 31.6, ...
, 16.790, 22.6, 支間中央1~4 , 23.0, 支間中央1~4 , 23.5, ...
, 18.980, 17.6, 支間中央1~5 , 17.9, 支間中央1~5 , 18.4, ...
, 21.170, 14.1, 支間中央1~4 , 14.3, 支間中央1~5 , 14.5, ...
**桁高600mm [#o6d177d7]
***鉛直材80mm・斜材80mm [#qaced9ef]
-橋長10.220mで，飛び出し発生!(3次モードは得られたけど)
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 61.4, 支点+鉛直材1 , 62.3, 支間中央1~3 , 64.9, ...
, 12.410, 42.5, 支間中央1~3 , 44.7, 支間中央1~4 , 46.6, ...
, 14.600, 30.8, 支間中央1~3 , 32.2, 支間中央1~4 , 32.9, ...
, 16.790, 23.2, 支間中央1~3 , 24.1, 支間中央1~4 , 25.0, ...
, 18.980, 18.1, 支間中央1~4 , 18.8, 支間中央1~5 , 19.3, ...
, 21.170, 14.5, 支間中央1~3 , 15.0, 支間中央1~4 , 15.4, ...
***鉛直材95mm・斜材95mm [#s71a8d20]
-橋長10.220mで，飛び出し発生! というか，このケースのほと...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 65.2, 支間中央1~3 , ---- , 飛び出し! , ---- , ...
, 12.410, 44.2, 支間中央1~3 , 46.7, 支間中央1~4 , 47.8, ...
, 14.600, 32.0, 支間中央1~3 , 33.4, 支間中央1~4 , 34.3, ...
, 16.790, 24.1, 支間中央1~3 , 25.2, 支間中央1~4 , 25.9, ...
, 18.980, 18.8, 支間中央1~4 , 19.5, 支間中央1~5 , 20.0, ...
, 21.170, 15.1, 支間中央1~3 , 15.6, 支間中央1~4 , 16.1, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#gbaf52d5]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 64.8, 支間中央1~4 , ---- , ---- , ---- , ---- ,...
, 12.410, 43.5, 支間中央1~3 , 45.2, 支間中央1~4 , 47.2, ...
, 14.600, 31.5, 支間中央1~4 , 32.8, 支間中央1~5 , 34.4, ...
, 16.790, 23.6, 支間中央1~4 , 24.4, 支間中央1~4 , 25.1, ...
, 18.980, 18.4, 支間中央1~5 , ---- , ---- , ---- , ---- ,...
, 21.170, 14.7, 支間中央1~5 , 15.1, 支間中央1~5 , 15.5, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#r9fef15b]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 66.9, 支間中央2~4 , 72.3, 支間中央1~4 , 72.6, ...
, 12.410, 45.5, 支間中央1~3 , 47.4, 支間中央1~4 , 48.8, ...
, 14.600, 32.8, 支間中央1~4 , 34.4, 支間中央1~5 , 35.1, ...
, 16.790, 24.6, 支間中央1~4 , 25.5, 支間中央1~5 , 26.2, ...
, 18.980, 19.2, 支間中央1~5 , 19.9, 支間中央1~6 , 20.6, ...
, 21.170, 15.3, 支間中央1~5 , 15.8, 支間中央1~5 , 16.2, ...
**桁高700mm [#b1b28363]
***鉛直材80mm・斜材80mm [#k96f9889]
-橋長10.220mで，いくつか飛び出しモードが出てしまった．
--Sub Mesh の優先順位を1:境界面 2:木材 に変えた方がいいか...
--double mesh がどうのこうのと書かれているので，remove do...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 49.3, 支点+鉛直材1 , 63.6, 支点+鉛直材1 , ---- ...
, 12.410, 40.7, 支点+鉛直材1 , 44.2, 支間中央1・2 , 46.1,...
, 14.600, 31.9, 支間中央1~3 , 33.1, 支間中央1~4 , 33.9, ...
, 16.790, 24.1, 支間中央1~3 , 25.0, 支間中央1~4 , 25.8, ...
, 18.890, 18.8, 支間中央1~3 , 19.5, 支間中央1~4 , 19.9, ...
, 21.170, 15.1, 支間中央1~3 , 15.5, 支間中央1~4 , 16.0, ...
***鉛直材95mm・斜材95mm [#t274bf45]
-やはり橋長10.220mは飛び出しが多い．経験的に座屈しにくい...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 55.0, 支点+鉛直材1 , 63.4, 支点+鉛直材1 , 66.0 ...
, 12.410, 45.0, 支間中央1・2 , 45.5, 支点+鉛直材1 , 47.3,...
, 14.600, 32.4, 支間中央1~3 , 33.9, 支間中央1~4 , 34.8, ...
, 16.790, 24.6, 支間中央1~3 , 25.6, 支間中央1~4 , 26.5, ...
, 18.890, 19.2, 支間中央1~3 , 19.9, 支間中央1~4 , 20.4, ...
, 21.170, 15.4, 支間中央1~3 , 15.9, 支間中央1~4 , 16.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#lf7e27f6]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 67.4, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 12.410, 45.5, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 14.600, 32.9, 支間中央1~4 , 34.0, 支間中央1~5 , ---- , ...
, 16.790, 24.7, 支間中央1~4 , 25.3, 支間中央1~4 , 26.0, ...
, 18.890, 19.3, 支間中央1~5 , 19.9, 支間中央1~5 , 19.8, ...
, 21.170, 15.4, 支間中央1~4 , 15.7, 支間中央1~5 , 16.1, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#mc6a2237]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 69.4, 支間中央2~3 , 73.2 , 支間中央2~3 , 76.0, ...
, 12.410, 46.6, 支間中央1~3 , 47.9, 支間中央1~4 , 50.3, ...
, 14.600, 33.7, 支間中央2~4 , 34.9, 支間中央2~5 , 36.9, ...
, 16.790, 25.3, 支間中央1~4 , 25.8, 支間中央1~4 , 26.6, ...
, 18.890, 19.7, 支間中央1~4 , 20.2, 支間中央2~5 , 20.9, ...
, 21.170, 15.8, 支間中央1~4 , 16.1, 支間中央1~5 , 16.5, ...
**桁高800mm [#ib28b7b4]
***鉛直材80mm・斜材80mm [#abad02dc]
-''ここの表だけ木材の幅を420mmに間違えて(0.1mm削ると419.9...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 39.3, 支点+鉛直材1 , 51.1, 支点+鉛直材1・2 , 54...
, 12.410, 32.5, 支点+鉛直材1 , 41.9, 支点+鉛直材1・2 , 43...
, 14.600, 27.7, 支点+鉛直材1 , 32.0, 支間中央1・2 , 33.0,...
, 16.790, 24.2, 支点+鉛直材1 , 24.5, 支間中央1~2 , 25.2, ...
, 18.980, 18.9, 支間中央1~3 , 19.5, 支間中央1~4 , 19.8, ...
, 21.170, 15.4, 支間中央1~3 , 15.7, 支間中央1~4 , 16.0, ...
***鉛直材95mm・斜材95mm [#w6ee712b]
-橋長10mは飛び出しモードも算出された．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 44.3, 支点+鉛直材1 , 62.4, 支点+鉛直材1・2 , 64...
, 12.410, 36.6, 支点+鉛直材1 , 44.8, 支点+鉛直材1・2 , 43...
, 14.600, 31.3, 支点+鉛直材1 , 32.1, 支間中央1・2 , 33.1,...
, 16.790, 24.5, 支間中央1~3 , 25.2, 支間中央2~3 , 26.0, ...
, 18.980, 19.0, 支間中央1~3 , 19.5, 支間中央1~4 , 19.9, ...
, 21.170, 15.4, 支間中央1~3 , 15.7, 支間中央1~4 , 16.1, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#f08f6c00]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 54.8, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 12.410, 44.1, 鉛直材2~3 , 45.8, 鉛直材1~5 , ---- , ----...
, 14.600, 33.3, 支間中央2~3 , ---- , ---- , ---- , ---- ,...
, 16.790, 25.2, 支間中央1~3 , 25.3, 支間中央1~3 , ---- , ...
, 18.980, 19.6, 支間中央2~3 , 19.9, 支間中央2~4 , 20.4, ...
, 21.170, 15.8, 支間中央1~3 , 15.8, 支間中央1~4 , 16.1, ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#s8e16104]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 68.6, 支間中央2~3 , 71.9, 支間中央2~4 , ---- , ...
, 12.410, 46.4, 支間中央1~2 , 47.3, 支間中央1~3 , 49.4, ...
, 14.600, 33.4, 支間中央2~3 , 34.4, 支間中央2~4 , 35.9, ...
, 16.790, 25.3, 支間中央1~3 , 25.5, 支間中央1~4 , 26.2, ...
, 18.980, 19.6, 支間中央2~3 , 20.0, 支間中央2~4 , 20.6, ...
, 21.170, 15.8, 支間中央1~3 , 15.9, 支間中央1~4 , 16.2, ...
**桁高900mm [#i96aa86f]
***鉛直材80mm・斜材80mm [#f36b9e88]
桁高900mmから計算時間が5倍ぐらい長くなった．多分，要素分...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 32.4, 支点+鉛直材1 , 42.3, 支点+鉛直材1・2 , 45...
, 12.410, 26.8, 支点+鉛直材1 , 34.7, 支点+鉛直材1~3 , 36....
, 14.600, 22.9, 支点+鉛直材1 , 29.5, 支点+鉛直材1~3 , 30....
, 16.790, 19.9, 支点+鉛直材1 , 24.3, 支間中央1~3 , 25.0, ...
, 18.980, 17.6, 支点+鉛直材1 , 18.5, 支間中央1・2 , 19.2,...
, 21.170, 15.3, 支間中央1~3 , 15.6, 支間中央1~4 , 15.8, ...
***鉛直材95mm・斜材95mm [#k2393a16]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 36.7, 支点+鉛直材1 , 52.0, 支点+鉛直材1・2 , 62...
, 12.410, 30.4, 支点+鉛直材1 , 42.9, 支点+鉛直材1~2 , 43....
, 14.600, 26.0, 支点+鉛直材1 , 20.8, 支間中央1~3 , 32.3, ...
, 16.790, 22.9, 支点+鉛直材1 , 23.9, 支間中央1~3 , 24.7, ...
, 18.980, 18.2, 支間中央1~3 , 18.9, 支間中央1~5 , 19.5, ...
, 21.170, 15.0, 支間中央1~3 , 15.4, 支間中央1~5 , 15.8, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#m49094a7]
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 45.7, 鉛直材2~3 , 47.7, 鉛直材1~4 , ---- , ----...
, 12.410, 36.6, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 14.600, 30.7, 鉛直材2~3 , ---- , ---- , ---- , ---- ,, ...
, 16.790, 24.8, 支間中央1~3 , ---- , ---- , ---- , ---- ,...
, 18.980, 19.3, 支間中央2~4 , ---- , ---- , ---- , ---- ,...
, 21.170, 15.5, 支間中央1~4 , 15.8, 支間中央1~5 , ---- , ...
-ちなみに，特に橋長10mなどの短いものは，変位が最大となる...
***鉛直材95mm・斜材95mm(対傾構2つ) [#a5b90213]
-橋長10mの2次モードは，鋼板全体が波々に座屈する．
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 66.7, 支間中央2~4 , 69.9, 鋼板全体 , ---- , ---...
, 12.410, 44.9, 支間中央1~3 , 46.9, 支間中央1~4 , 48.8, ...
, 14.600, 32.5, 支間中央2~4 , 33.8 , 支間中央2~5 , 35.4 ,...
, 16.790, 24.3, 支間中央1~3 , 25.2, 支間中央1~4 , 25.9, ...
, 18.980, 19.0, 支間中央2~4 , 19.6, 支間中央2~5 , 20.2, ...
, 21.170, 15.2, 支間中央1~4 , 15.6, 支間中央1~5 , 16.0, ...
**桁高1000mm [#ibd647a4]
***鉛直材80mm・斜材80mm [#f8c140d5]
橋長20m/18mなどは計算時間がとても長い．要素数と要求次数の...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 27.3, 支点+鉛直材1 , 35.8, 支点+鉛直材1~3 , 38....
, 12.410, 22.6, 支点+鉛直材1 , 29.3, 支点+鉛直材1~3 , 31....
, 14.600, 19.3, 支点+鉛直材1 , 24.9, 支点+鉛直材1~3 , 26....
, 16.790, 16.8, 支点+鉛直材1 , 21.5, 支点+鉛直材1~3 , 22....
, 18.980, 14.9, 支点+鉛直材1 , 17.9, 支間中央1~3 , 18.8, ...
, 21.170, 13.3, 支点+鉛直材1 , 14.9, 支間中央1~4 , 15.4, ...
***鉛直材95mm・斜材95mm [#i2fb2ddc]
鉛直材80mm・斜材80mmと比較しやすいように最も太いものを計...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 31.2, 支点+鉛直材1 , 44.3, 支点+鉛直材1~2 , 60....
, 12.410, 25.9, 支点+鉛直材1 , 36.5, 支点+鉛直材1~2 , 42....
, 14.600, 22.1, 支点+鉛直材1 , 29.5, 支間中央1~3 , 31.1, ...
, 16.790, 19.3, 支点+鉛直材1 , 23.1, 支間中央1~4 , 24.3, ...
, 18.980, 17.1, 支点+鉛直材1 , 17.4, 支間中央1~3 , 18.4, ...
, 21.170, 14.4, 支間中央1~4 , 15.0, 支間中央1~5 , 15.4, ...
***鉛直材80mm・斜材80mm(対傾構2つ) [#i6f228f2]
3次モードぐらいまで確認したかったけど，もう期限的にその余...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 38.8, 鉛直材2~3 , 40.6, 鉛直材1~4 , ---- , ----...
, 12.410, 31.1, 鉛直材2~4 , ---- , ---- , ---- , ---- ,, ...
, 14.600, 26.0, 鉛直材2~4 , 27.1, 鉛直材1~4 , 27.5, 鉛直...
, 16.790, 22.2, 鉛直材2~4 , 23.1, 鉛直材2~5 , ---- , ----...
, 18.980, 18.8, 支間中央2~5 , ---- , ---- , ---- , ---- ,...
, 21.170, 15.0, 支間中央1~4 , ---- , ---- , ---- , ---- ,...
-ちなみに，橋長10mなどの短いものは，変位が最大となるとこ...
***鉛直材95mm・斜材95mm(対傾構2つ) [#p9977548]
-橋長10mの1次モードは，鉛直材というよりも，支点付近の鋼板...
, 橋長(m) , 1次(m) , 箇所 , 2次(m) , 箇所 , 3次(m) , 箇所...
, 10.220, 62.1, 鉛直材1~6 , 64.6, 支間中央2~5 , 66.0, 鉛...
, 12.410, 43.5, 支間中央1~4 , 46.1, 支間中央1~5 , ---- , ...
, 14.600, 31.4, 支間中央2~5 , 33.2, 支間中央2~6 , 35.0, ...
, 16.790, 23.5, 支間中央1~4 , 24.6, 支間中央1~5 , 25.8, ...
, 18.980, 18.3, 支間中央2~5 , 19.2 , 支間中央2~6 , 20.0, ...
, 21.170, 14.6, 支間中央1~5 , 15.2, 支間中央1~6 , 15.8, ...
**グラフ [#r2857111]
手始めに構造研でお馴染みの3次元プロットでもつくってみる．...
***鉛直材80mm・斜材80mm [#g5209b20]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-発表や論文には使えないが，ブラウザ上でマウスで動かせたら...
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-もしかすると，色はz軸の値にするのではなくて，座屈箇所に...
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-なんだか，surfaceのx軸とy軸がよく分からないので，mesh3d...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材95mm・斜材95mm [#o746beae]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-ブラウザ上で表示
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-ブラウザ上で表示(座屈箇所ごとに色分け)
--http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/g...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材80mm・斜材80mm(対傾構2つ) [#fd860625]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
***鉛直材95mm・斜材95mm(対傾構2つ) [#yda7b568]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
-設計参考資料として使いやすいように，2次元グラフも作って...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/gra...
*Salome-Meca2017の座屈解析メモ [#t6be1072]
**固有値が算定されるのに結果が見れない場合 [#z3abdde0]
固有値が近くて，誤差が許容範囲以下に収まらない時，message...
-座屈解析の CALC_MODES を編集 -> VERI_MODES にチェックを...
-ただし，計算されたモードは自己責任で評価して下さい，との...
-また，medファイルが重くて？モデリングや計算に使ったhdfフ...
**TRI_DIAG(Lanczos法)の精度を上げる方法 [#kfe64a70]
Code_Asterの座屈解析では，誤差(norme d'erruer)が0.000001...
-座屈解析の CALC_MODES を編集 -> SOLVEUR_MODAL の Edit を...
--COEF_DIM_ESPACE は，デフォルトでは4なので，それを適宜上...
---例えば求めたい次数が6のときに，COEF_DIM_ESPACEがデフォ...
--ただし，要求次数が6のとき，COEF_DIM_ESPACEをデフォルト...
--と思いきや，要求次数6のとき，試しにCOEF_DIM_ESPACEを12...
--さらに，COEF_DIM_ESPACEを24にしてみたら，計算時間がデフ...
--ということで，精度を上げるためにはCOEF_DIM_ESPACEを大き...
-PREC_LANCZOS(精度(precision)だと思う)を変えても，結果は...
-CalculiXでは Lanczos vectors はデフォルトで4$\times$次数...
--デフォルトで4$\times$次数とは！なかなか多いですね．しか...
--(
*BUCKLE [#e09c3fd1]
6, 0.01, 24, 1000
**次数，精度，Lanczos vectors，最大反復回数…この行はただ...
--)
**固有値解析 [#p2357bac]
-昨年度の方法を真似てSORENSEN法？とやらを使っているが，他...
-よくわからないが，アプローチ？を追加してみる．
--REELとIMAGで同じ結果になり，COMPLEXEは結果が出ない．ど...
***SORENSEN [#jdf3999b]
IRAM(Implicitly Restarted Arnoldi Method)というらしく(Sor...
-後藤さんいわく，何かの間違いで普通でない非対称の行列がで...
***JACOBI [#kcb8e997]
ただのjacobi法かと思いきや，公式ドキュメントにBatheとWils...
-後藤さんによると，時間がかかっているからといって必ずしも...
***TRI_DIAG [#vdb3c836]
名前はただの3重対角化っぽいけど，Lanczos法で実対称行列の...
-原則として，TRI_DIAGで計算することにしようかと思う．
-TRI_DIAGの精度(PREC_LANCZOS)はデフォルトでは1.0E-08だが...
--CaluliXとオンサイト木橋の座屈を比較してみたところ，座屈...
---ところで，CalculiXでは荷重を節点数で割っただけで，等価...
-試しにPREC_LANCZOSを1.0E-02に変えてみても，結果は何一つ...
-messageファイルにCOEF_DIM_ESPACEを上げてみては？と書かれ...
--COEF_DIM_ESPACEとは&link(ここ,http://femingway.com/?p=2...
-一応2・3次モードぐらいまで確認しておきたいので，正負のモ...
--デフォルトのCOEF_DIM_ESPACE:4から6に上げてみたところ，2...
***QZ [#q84ea5c5]
そのままQZ法．大規模行列には向いていない(自由度1000以下)...
**幾何剛性マトリックス [#f14a6fd4]
昨年度の方法では，幾何剛性マトリックスを作るために，"SIEF...
**載荷方向と座屈荷重の正負 [#u13f97ff]
なぜか逆向きになり，使用者がこんがらかる問題．
-片持ち梁のオイラー座屈を色々な組み合わせて計算してみたと...
-トラヘアの座屈を載荷方向を変えて計算したところ，載荷方向...
**CalculiXとかMarcとかのこと [#m806e0a9]
-CalculiXは固有値解析(というかほとんどの行列計算)をするの...
*オンサイト木橋の鋼板をSalome上で自動生成するPythonスクリ...
オンサイト木橋の鋼板について，図面のなかには情報が十分で...
ところで，後藤さんからオンサイト木橋の鋼板の形状を決める...
-&link(onsi_kou.py,http://www.str.ce.akita-u.ac.jp/~gotou...
--溝形のフランジ付きモデル．桁高・鉛直材幅・斜材幅・橋長...
--近い将来(少なくとも2020年頃までに)，Salome-MecaもPython...
--寸法(弦材幅や斜材の角度)は，Brepファイルにエクスポート ...
-&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~goto...
--座屈解析用に鋼板，木材，支点の線を作る．鋼板は&link(ons...
--次のように要素分割すれば，飛び出しが生じにくいというこ...
---鋼板 : NETGEN 1D-2D-3D -> Max.size : 20 , Min.size : 10
---木材(SubMesh) : NETGEN 1D-2D-3D -> Max.size : 60 , Min...
---境界面(SubMesh) : NETGEN 1D-2D -> Max.size : 10 , Min....
---Convert to Quadratic
-&link(onsiz.py,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
--(やり直す時間的余裕がなかったので)今回の計算には使って...
-&link(onsitz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
--&link(onsiz_org.py,http://www.str.ce.akita-u.ac.jp/~got...
-&link(onsitz.py,http://www.str.ce.akita-u.ac.jp/~gotouha...
--&link(onsitz_org.py,http://www.str.ce.akita-u.ac.jp/~go...
-海老さんと川村さんによると，4面体に要素分割をするとき，M...
*ねじ山を作りたくない [#n60cf996]
菊地翔さんのPC鋼棒の緊張力ネタで「今後はねじ山を考慮した...
&link(McMaster(機械部品屋)のページ,https://www.mcmaster.c...
*軸方向に色々な荷重が作用した時の柱の座屈 [#fac7ce18]
**片持ち梁の先端に集中荷重が作用した場合（メモ） [#t20b94...
そこかしこに載っているもの．
-$P_{cr}=\frac{\pi^{2} EI}{4\ell^2}\approx \frac{2.467EI}...
**片持ち梁に等分布荷重が作用した場合（メモ） [#bfb3d57d]
ティモシェンコ『弾性安定の理論』の上巻から引用．
***Bessel関数を解いたら(p.97) [#kbea664c]
Bessel関数とやらの解法がよく分からない．
-$(q\ell)_{cr}=\frac{7.837EI}{\ell^2}$
***エネルギー法1(p.99) [#m421c3a4]
第一近似として，たわみを$y(x)=\delta (1-\cos \frac{\pi x}...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$(q\ell)_{cr}=\frac{7.89EI}{\ell^2}$
***エネルギー法2(p.101) [#fcb9d5b1]
さらに助変数を用いて，たわみを$y(x)=\delta_{1}(1-\cos \fr...
-$(q\ell)_{cr}=\frac{7.84EI}{\ell^2}$
**片持ち梁に等分布荷重（圧縮方向）と先端に集中荷重（圧縮...
等分布荷重qは，集中荷重Pと係数kを用いて，仮に$q=\frac{kP}...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$P_{cr}=\frac{(6\pi^5+(3\pi^5-12\pi^3)k)EI}{((-1536+432\...
--試しに，k=0を代入すると，先端に集中荷重が載った場合の座...
-ちなみに，等分布荷重が引張方向に作用していた場合，kの値...
-$P_{cr}=\frac{(6\pi^5-(3\pi^5-12\pi^3)k)EI}{((-1536+432\...
-計算が間違っている可能性も高いので，そのうちにFEMで確認...
**片持ち梁に三角形分布荷重（圧縮方向）と先端に集中荷重（...
図がないと載荷条件が分かりにくいが，描くのが面倒臭いので...
三角形分布荷重qは，集中荷重Pと係数kを用いて，仮に$q(x)=\f...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
-$P_{cr}=\frac{(30\pi^6+(5\pi^6-30\pi^4)k)EI}{((-4080+360...
**片持ち梁に三角形分布荷重（引張方向）と先端に集中荷重（...
図がないと載荷条件が分かりにくいが，描くのが面倒臭いので...
三角形分布荷重qは，集中荷重Pと係数kを用いて，仮に$q(x)=\f...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
-$P_{cr}=-\frac{\left(\left(5{{\pi}^{6}}-15{{\pi}^{4}}\ri...
**片持ち梁に台形分布荷重（圧縮方向）と先端に集中荷重（圧...
台形分布の下底（固定端側）を$\frac{k_{1}P}{\ell}$，下底（...
そして，ティモシェンコのエネルギー法を真似，たわみを第一...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-$P_{cr}=\frac{\left(30{{\pi}^{6}}+\left(5{{\pi}^{6}}-30{...
-maximaに貼り付ける用
--(
P1(k1,k2,l1):=((30*%pi^6+(5*%pi^6-30*%pi^4)*k1+(10*%pi^6-...
--)
--$k_{2}$に0あるいは$k_{1}$に0代入したら三角形分布と同じ...
--$k_{1}=k_{2}$にしたら等分布と同じになる．
**片持ち梁に台形分布荷重（引張方向）と先端に集中荷重（圧...
-↑の式をk1とk2をマイナスにしたら同じになる.
-$P_{cr}=\frac{\left(-30{{\pi}^{6}}-\left(5{{\pi}^{6}}-30...
-つまり，maximaでは，P1を貼り付けたあとに次のようにすれば...
--(
P2(k3,k4,l2):=P1(-k3,-k4,l2)
--)
**両端固定の問題に応用したい場合 [#dc6d4761]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/bun...
-まず，上の片持ちのそれぞれの式（分布荷重が圧縮方向と引張...
-そして，Pcrをそれぞれ$Pcr_{1},Pcr_{2}$，スパン$\ell$をそ...
--連立方程式：$Pcr_{1}=Pcr_{2}, \ell_{1}+\ell_{2}=0.5$
---ここでは，具体的な値ではなく，便宜上$\ell_{1}+\ell_{2}...
-より，$\ell_{1}, \ell_{2}$を求め，元の式$Pcr_{1}$か$Pcr_...
***等分布載荷 [#d19735da]
kと$\ell_{2}$を$Pcr_{2}$に代入すればいい．$\ell_{2}$は，...
-かなり長いので，書きたくない．maximaに貼る用のみ．
--(
sol(k1,k2):=solve(P1(k1,k1,0.5-l2)=P2(k2,k2,l2),l2)
--)
--0.5は仮で，本当は$\frac{\ell}{2}$をいれる．別に$\ell_{2...
--等分布荷重が載っていない場合と同じになるかの確認のためk...
--(
bfloat(P2(k2,k2,sol(k1,k2)))
--)
--ここでは，kに具体的な数値を入れる．
--kの計算方法については考え中．多分，曲げモーメントによる...
***台形分布載荷 [#tde5cd2a]
kと$\ell_{2}$を$Pcr_{2}$に代入すればいい．$\ell_{2}$は，...
-かなり長いので，書きたくない．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...
--(
bfloat(P2(k3,k4,sol(k1,k2,k3,k4)))
--)
--ここでは，k1, k2, k3, k4に具体的な値を入れる．
--k1,k2,k3,k4の計算方法については考え中．
*直交異方性材料の直応力-直ひずみ関係 [#jfd24410]
前々から&link(構造関係メモ,http://www.str.ce.akita-u.ac.j...
-式(1)について、Maximaで求めた結果(&link(ihou1.wxm,http:/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
-確かに複雑……。
*木材利用2017メモ [#k31eb3da]
-アンカープレートの部分だけ木材と鋼板を一体化させる。
--アンカープレートの境界面でパーティションし、木材はLengt...
---これだと、鋼板の隣の非一体化部分が木材にすり抜けるので...
-堅固山さんとの違いは、
--木材を直交異方性にする。
---等方性と異方性にしたときの座屈荷重の比較をする?
--3点曲げに加えて積雪荷重での座屈も計算する。
--特に、支点付近の鉛直材に焦点を当てる。
---つまり、中央以外の非一体化部分もモデリングする。
--大池や虎毛などの大きいモデルでもやってみる。
---要素数の問題で部分的に取り出したモデルで計算するかも。
**5月初旬(発表概要)までに [#z3f1729b]
-7mモデルでやる。
--穴の拡大率は0,5,10,15,20,25,30%の7通りとし、7mモデルの...
--何もなしモデルと補剛材+対傾構ありモデルの2通りでやって...
-載荷条件は、3点曲げ(線載荷でいいかな?)と積雪荷重(積雪1m...
これだけで28通り。ここまでの計算自体はjobを実行すれば2～3...
-さらに木材を等方性と異方性にしたものとで比較する?ことに...
-中央と支点部の鋼板と木材だけ非一体化にしたものでもやって...
ひとまず、7mスパンについては全通りの計算が終わったが、結...
**後藤さんと相談した概要の流れをメモ [#f0170a09]
-1.端だけ非一体化・等方性で穴の大きさを変える(堅固山さん...
-2.しかし、木材を直交異方性にすると、それよりも座屈荷重が...
--不思議。それと、要素の跳び出しが等方性の時よりも増える...
-(3.)さらにアンカープレートのところだけ一体化にすると、座...
--しかも、木材と鋼板の境界に応力集中する箇所が増えるから...
-4.実際に影響が大きいと予想される積雪荷重が載った場合、3...
-対傾構と補剛材の影響をどこに挟むのがいいのかを訊くのを忘...
**6月から8月(発表スライドor論文)までに [#h1d08971]
-大池モデル(や虎毛モデル)でも同じことをやる。
-7mモデルで支点近くだけを取り出した一部モデルと従来の1/4...
こんなにやれるかな？
**おまけ [#afc06b56]
-5mスパン・板厚6mmの図面があったら(k2にはない?)、本当に70...
--座屈荷重は，部分的非一体化モデルで282kN, 全体的非一体化...
---というか，そもそも実験値（tugite）のはらみ（7mに3点曲...
-ちなみに，CalculiXで座屈の*STEP～*END STEPの前にプレスト...
-どうやら非線形座屈解析の出番はないかも。計算するにしても...
--幾何学非線形で計算すると，座屈荷重に至る前にまた別のあ...
-計算に使った積雪1mにつき$3.5kN/m^2$だが、昔の道路橋示方...
,降りたての雪, $150kg/m^3$
,やや落着いた雪, $300kg/m^3$
,圧縮された雪または多量に水を含んだ雪, $500～700kg/m^3$
*印刷用の折り紙円筒 [#s8a60d5d]
色々な折り目パターンに変えても、$\LaTeX$を経由してA4用紙...
ところで，7月10日の高校の先生との個別相談会用に印刷してみ...
-システム-->システム管理-->プリンター
-HL5270DNを右クリックして，プロパティ
-ジョブオプション-->「用紙に合わせる」にチェック-->適用
-&link(prjabara.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--蛇腹折り。例えば、&link(周方向6分割・高さ方向8分割の例,...
近藤夏実さんが折り畳まることを確認済み．
-&link(prkrasen.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--回転螺旋折り。例えば、&link(周方向5分割・高さ方向10分割...
近藤夏実さんに折って貰ったところ，あまりにも糊付けが難し...
-&link(prhrasen.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--反転螺旋折り。例えば、&link(周方向5分割・高さ方向6分割...
近藤夏実さんが折り畳まることを確認済み．こちらは糊付けが...
---螺旋折り系は、紙に目一杯描くようにしていないため、パタ...
-&link(prdaiya.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--ダイヤカット円筒。例えば、&link(周方向6分割・高さ方向8...
川村さんに折って貰ったところ，ちゃんと完成する．（これは...
*支間中央のみを取り出して線形座屈解析 [#pc2d0104]
局部座屈の着目点と比較して全体が大きい構造物だと、十分な...
手法は以下の通り。
-拘束は
--(z)橋軸方向対称面
--(x)幅員方向対称面
--(y)切断面の中立軸
-載荷は
--(y)輪荷重P/4
--(z)$\frac{曲げモーメント_{切断面}}{角材同士の中心間隔}$...
---ここで、角材と鋼板の剛性比($\frac{E_{鋼}I_{鋼}}{E_{鋼}...
-孔空き鋼板の3点曲げ(載荷は輪荷重)とした。
-モデルはk2に上がっていた堅固山さんのbrepファイルをそのま...
-堅固山さんが計算したものは、$P_{cr}$=270kNなので、1/4モ...
,, 1/4モデル(等方性) , 一部モデル(等方性) , 1/4モデル(異...
, 要素(C3D10)数 , 175493, 92862, 175493, 92862
, 1次モード , 270.0kN, 263.6kN, 228.5kN, 232.8kN
, 2次モード , 871.2kN, 851.3kN, 728.6kN, 745.7kN
, 3次モード , 1769kN, 1729kN, 1494kN, 1536kN
--ちなみに、応力集中する要素の飛び出しなどのノイズと思わ...
--モデルの体積比は$\frac{1/4モデル}{一部モデル}=5$なので...
-変位図(変位の倍率は見やすいように調整)
--そんなことをしていいのだか分からないけど、1/4モデルの1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ooi...
--相対誤差が大体2~3%ぐらいなので、なかなか良い精度だと思...
-分布荷重が載った場合についても、いずれやりたい。
**後藤ちゃちゃ(17/3/19) [#f902c5d0]
トラスの影響線で考えると、圧縮力が最も大きくなる鉛直材は...
一番近いところになるので、実は、支点側の座屈を解析する方...
その辺は、後でお話しましょう。
ところで、きのう、岩熊先生の退職パーティーに行って、鬆徒...
&link(3.5.2.4 積層板の見かけ上のYoung率,http://mechanics....
-ただの近藤メモ:堅固山さんに聞いてみたところ、純粋な鉛直...
--後藤ちゃちゃ(17/3/23)虎毛や大池のモデルは、鋼板上下が折...
*理工学デザイン [#k0eaa856]
--(
システムデザイン工学専攻 １２月２２日（木）１２時５０分～...
○会場：ＶＢＬ棟２階大セミナー室
○発表者は開始５分前までに準備を終えて会場に集合のこと。
○発表者は３名（教員２名、学生１名）への説明が終わった後も、
他の学生のポスターを見学したり、担当でない先生方へも対応...
終了時間まで会場に残ってください。
○当日、発表も質疑もない学生は出席しなくて構いません。
○ポスター発表当日の午前１１時からポスターの貼り付けなど準...
○ポスター用ボードにポスター発表番号が記載されていますので、
指定の場所にポスターを貼り付けてください。
○各自、ピンかセロテープを用意して、ポスター用ボードに貼り...
○終了後は、忘れずにポスターを持ち帰りください。
○Ａ１サイズ（縦８４１ｍｍ×横５９４ｍｍ）縦置きでポスター...
○ポスターの上段に、タイトル、所属（○○専攻○○コース）、
発表者氏名を明記してください。
○ポスターの作成は自由（手作りでも結構です）に行ってくださ...
情報統括センターを利用してポスターを作成する人は、
予めセンターに予約（前日17:00まで）してから行ってください。
なお、ポスター印刷する際は、利用申請書を必ず持参し...
○研究テーマや方向性、方法をポスターにまとめ、発表する。
指定された教員２名に説明する以外に、他専攻の学生に対して...
○各学生は、自身の属する専攻以外の学生のポスター発表を訪れ、
質疑応答し、その学生の研究内容を十分に理解した上で研究内...
1か月以内に理工（大学院）担当 に提出する。
特に、複数の質問事項を明確に記載し、質疑応答を記録する。
その質疑応答から理解した事をレポートにまとめる。
--)
-後藤ちゃちゃ(16/11/18)
各種の非線形問題（ねじり座屈とか接触とか）一通り解けてい...
（CalculiXも一通り解けそうですね）、線形でも回転螺旋折り...
解析が困難な各種の非線形問題を工夫しながら解いてますが、
オープンソースのツールだけで、割と何でも解けてます的なオ...
話にしてもいいかもしれません。まだM1なんで色々とおもしろ...
やってます。これからやっていきたい本命はこれなんですが、...
入れてもいいですが。
*PC鋼棒の引張+トルク+摩擦 [#i2a752fc]
-&link(toruku.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--PC鋼棒を引張った後、トルクを導入し、滑り出しをみる。
--PC鋼棒(梁要素)の先端に円盤(立体要素)をくっつけて、もう...
---菊地さんのと境界条件や荷重条件がやや異なるので、いずれ...
-&link(tth.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--gnuplot用にdatをまとめる。
*CLT関連で使えそうな板の文献をメモ [#i726a43d]
CLTのたわみ(実験値、FEMの値)が、どの理論式に近くなるのか...
-土木学会『構造力学公式集』@研究室
--相対2辺が単純支持された等方性長方形薄板：p.328
---部分分布荷重：p.330
---残りの相対2辺が弾性支持桁：p.331 .ここで、$\gamma=0$と...
--直交異方性板：p.367 .少しだけ長方形板の記述あり。ただし...
-ティモシェンコ『板とシェルの理論(上・下)』 @図書館
--直交異方性板は下巻に書かれている。これも具体例はあまり...
-渡辺昇『橋梁工学』 @研究室
--等方性板：p.300
--直交異方性板：p.324（そんなに詳しくない…）
-アムバルツミャン『異方弾性板の理論』 @図書館
--異方性の薄板、厚板、積層板について座屈、有限変位、振動...
-J.R. Vinson & R.L.Sierakowski『複合材料の構造力学』@図書館
--異方性の積層板について：pp.42-101
--さらに、7章に「強度と破損の理論」が載っている。これは面...
-宮入裕夫『複合材料入門 ー基礎と応用ー』@図書館
--異方性の積層板についてはあまり書かれていないのだが，何...
--ということで，ハニカムパネルについて読んでみると，面内...
---ハニカムパネルについては，図書館にはないが，宮入裕夫『...
*板の座屈 [#zaff6db6]
-メモ：4辺単純支持の座屈荷重
$N=k\frac{D\pi^{2}}{b^{2}} .$
ここで、
$k=(m\frac{b}{a}+\frac{1}{m}\frac{a}{b})^{2} , D=\frac{Et...
t:板厚 a:長さ b:幅 m:長軸方向の凹凸の数
と、色んな資料に書かれているが、Dは単位幅あたりの曲げ剛性...
-&link(itaz.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--4辺単純支持を板要素(S8)で計算。幅厚比b/t=100で試してみ...
-&link(hitaz.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--itaz.f90の結果で非線形座屈計算。4辺単純支持した板は、そ...
*プレストレス木箱桁橋 [#ce0220fb]
**孔無し一体化モデルの1/4解析 [#b0cb7c7c]
-&link(onsiyon.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--プレストレス木箱桁橋プログラムの練習として一番単純なモ...
---3点曲げか、4点曲げか選べるように改良(『木材の特徴を活...
--1/4解析にしたのは、計算が大変な接触解析とか局部座屈とか...
**孔無し非一体化モデルの1/4解析 [#me62708c]
いわゆる接触解析。
-&link(onsimyon.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--プレストレスのかけ方は、鋼板の上下の木材と隣り合わせて...
---載荷板に載荷するように改良。変位計の位置で変位を出力す...
-&link(onsimayon.f90,http://www.str.ce.akita-u.ac.jp/~got...
--アンカープレートのみにプレストレスをかけてみたが、プレ...
-せっかくなので、スパン7m($E_{zz}$=7.694GPa)の孔無し一体...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--とりあえず摩擦係数は中央値をとって、鋼材と木材:0.40 , ...
--軸方向は128分割で、上村さんの修論本体のグラフからすると...
--変位図(載荷方向)は、10倍に拡大している。
--実験値は上村さんの修論本体の表の実験2(k2にあるuemurahon...
--濡れ具合の他に、↓摩擦面の仕上げ方によっても、摩擦係数は...
---ttps://www.jstage.jst.go.jp/article/aijs/76/666/76_666...
***緊張力の低下 [#q027a29a]
木材のクリープやらPC鋼棒のリラクセーションやらで緊張力が...
-&link(onsirelax.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(onsimyon.f90,http://www.str.ce.akita-u.ac.jp/~got...
---載荷板に載荷するように改良。変位計で変位を読み取るよう...
-とりあえず、スパン6m($E_{zz}$=6.921GPa)のモデルに荷重70k...
--スパン7mモデル($E_{zz}$=7.694GPa)にして、計算し直そう。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--上図は緊張力が53kNの時の載荷方向の変異(1倍)図。落下防止...
--破壊試験(100kN)では、鉛直変位の最大値が15mmぐらいなので...
--グラフには鉛直変位20mm以上は載せていない。データが途中...
--ところで，色々とひっくり返すようなことをいうと，緊張力...
**孔無し一部非一体化モデルの1/4解析(局部座屈) [#t76f6c01]
中央のアンカープレート間の鋼板のみ非一体化にして、かつ接...
-&link(onsikyo.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--プレストレスを鋼板の上下木材部分に全面載荷。もちろん、...
-&link(onsikyoa.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--プレストレスをアンカープレートのみに載荷。こちらの方が...
-これらをスパン7m($E_{zz}$=7.694GPa)、鋼板の板厚を9mmとし...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
,, 線形座屈(輝さん) , プレストレスなし , 全面プレストレス...
, 座屈荷重(kN) , 339.1(78.5) , 326.6 , 座屈しない? , 338.8
--下の図は500kNの時の幅員方向変位の1倍。当たり前だが、そ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ons...
--&link(斉藤輝さんの線形座屈解析の結果,http://www.str.ce....
--ただし、鋼板の板厚方向の要素分割は、今は適当に4分割だが...
-理工学デザインの発表用に、線形座屈と非線形座屈を同じ要素...
***オイラー座屈から3点曲げ載荷時の座屈荷重を推定 [#zb7e6f...
-&link(onsioira_3ten.f90,http://www.str.ce.akita-u.ac.jp/...
--堅固山さんが(孔空き鋼板の)オイラー座屈を計算しているの...
, 鋼板上部の座屈(h=163) , 鋼板上部の座屈(h=120) , 鋼板上...
, 216kN, 192kN, 218kN,, 6201kN, 2777kN
---「異部分」は，木材を直交''異''方性材料とした，''部分''...
---「等部分」は，木材を''等''方性材料とした，''部分''的非...
--鋼板上部の座屈について，ここでは(1.直交異方性を考慮して...
--ということで，FEMで中央の鉛直材と斜材をなくして，木材を...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/kai...
, 鋼板上部の座屈(h=163) , FEM(等部分) , FEM(異部分) , FEM...
, 216kN, 262kN, 218kN, 209kN, 170kN, 192kN, 215kN
--となることから，座屈荷重は，異方性の影響が約2割↓，中央...
--特に，異方性の中でも，木材の軸直角方向のヤング率の影響...
---ちなみに，ヤング率を変える（軸直交方向を硬くしても柔ら...
---ただし，桁高800mmになると中央の鉛直材と斜材の補剛効果...
--3点曲げ載荷における鉛直材の座屈は$\frac{\ell}{4}$あたり...
***オイラー座屈から等分布載荷時の座屈する積雪深を推定 [#w...
-&link(onsioira_yuki.f90,http://www.str.ce.akita-u.ac.jp/...
--&link(onsioira_3ten.f90,http://www.str.ce.akita-u.ac.jp...
, 鋼板上部の座屈(h=163) , 鋼板上部の座屈(h=120) , 鋼板上...
, 20.5m, 18.2m, 20.6m,, 658.7m, 295.0m, 278.8m
--鉛直材の長さの決め方が難しく，それによって結果が2倍ぐら...
-対傾構をつけず，補剛材をつけたらどうなるかやってみたい．...
-Trahairでも$S(\frac{\ell}{2})$や$\tau=\frac{QS}{bI}$とか...
**datファイルのデータ整理 [#u09e325b]
-&link(plonsi.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--一体化モデルの1/4解析用。
-&link(plonsim.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--非一体化モデルの1/4解析用。
-&link(plonsirelax.f90,http://www.str.ce.akita-u.ac.jp/~g...
--緊張力の抜けの1/4解析用。
*弾性床上の座屈 [#wf37705a]
-&link(yuka.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--特に初期不整を与えなくても、断面が変形して弾性床と接触...
---ただし、"*SURFACE BEHAVIOR"の値をかなり上手く調整しな...
---やってみたのは単純梁と片持ち梁だが、どちらも床と離れて...
--仮に床方向に座屈させたとしても、鋼材と木材(繊維直角方向...
*鋼材と木材の接触解析 [#o862c564]
-&link(masatu.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--&link(JIS K7125,https://www.djklab.com/service/koubunsh...
--鋼材と木材の摩擦係数は0.2(湿)~0.6(乾)の中央値をとって0....
--ちなみに木材(上に載っている方)の繊維方向は滑らす方向と...
--"*SURFACE BEHAVIOR"が謎。
---&link(Abaqusの"*SURFACE BEHAVIOR"の説明,https://engr-a...
---"*SURFACE BEHAVIOER"を&link(あまりいじらない方がいいと...
---要は収束の「する」、「しない」に関わってくる。
--"*FRICTION"の右のstick range slopeの値が謎。
---&link(マニュアル,http://web.mit.edu/calculix_v2.7/Calc...
---ネットを漁っていると、食い込みのない理想的な値は∞だが...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/mas...
--&link(pd.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/mas...
--縦軸が水平力、横軸が水平変位。変位は載荷面(横)の平均で...
*CLT [#n13ea13e]
**梁 [#l6cd48cb]
***梁(5層7プライ/7層7プライ)の3点曲げ [#u0bb2155]
-&link(clt.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似...
--ひき板の強軸方向ヤング率を試しに6GPaにして計算してみる...
***梁(5層7プライ/7層7プライ)の３点曲げを繰り返す(全体が弾...
-&link(cltso.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--木材利用研究論文報告集15のpp.11-16の佐々木さんのを真似...
--CalculiXの"*PLASTIC"はミーゼス応力での降伏値。木材は圧...
---"*PLASTIC"の行に続けて",HARDENING=～"を入れると、硬化...
--応力-ひずみ曲線(多直線近似)はかなり適当に入れてあるが、...
--"*PLASTIC"ではなく&link("*DEFORMATION PLASTICITY",http:...
-フィンガージョイントの降伏点が分かれば、そこだけ弾塑性材...
***h=240mm(7層8プライ)のときの梁理論との比較 [#dc468976]
-&link(clt3ten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(clt5ten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(pclt35ten.f90,http://www.str.ce.akita-u.ac.jp/~got...
--特に意味はないけど、&link(大黒屋さんの,http://library.j...
--桁高は重ね方の都合上、変えられないので、全方向の寸法を...
---つまり、b=24mm,h=240mm,$\ell$=240～4800mm(240mm刻みで...
--大黒屋さんに合わせて、分割数はnx=5,ny=24,nz=100.(大黒屋...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
--結果↑は大黒屋さんのものに似ている。
***b=24mm,h=240mm(7層8プライ)のときの1層と8層のFEMでの比...
-大黒屋さんの真似をしたついでに、FEMの相対誤差($\frac{v_{...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
--曲げが支配的になるほど、両者の値が近くなる。曲げ剛性か...
--弱軸の挙動もちょっと気になる。3点曲げと5点曲げで違う形...
**床版 [#pf637764]
''注意！''以下のプログラム類で1層化したモデルについては，...
$\nu_{zx}=\sqrt{0.0064\frac{E_{z}}{E_{x}}}=\sqrt{0.0064\f...
$\nu_{xy}=\nu_{xz}=0.0558, \nu_{zy}=\nu_{zx}=0.115$, $\nu...
に書き換えてみると，8層のものにかなり近い値がでるので，一...
***供試体2(7層8プライ/接続金物無し)の静的載荷 [#k6931ef5]
-&link(clt1.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--1層モデル。海老名さんの指摘により寸法や境界条件を直した...
---$\overline{E_z}=3.5$GPa, $\overline{E_x}=1.7$GPa, $E_y...
--- $\nu_{xy}=\nu_{xz}=0.4\cdot\frac{\overline{E_{x}}}{\o...
---$G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\ap...
--ちなみに、ヤング率の平均化は、$\frac{\sum_{i=1}^{8} E_{...
-&link(clt8.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--8層モデル。海老名さんの指摘により寸法や境界条件を直した...
---$E_z=5$GPa, $E_x=E_y=0.2$GPa.
---$\nu_{xy}=\nu_{xz}=\nu_{yx}=\nu_{yz}=\frac{0.4}{25}=0....
---$G_{xy}=G_{yz}=G_{zx}=\frac{E_{z}}{15}=\frac{5}{15}\ap...
-&link(cltiso.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--等方性材料。実験や上の2つとは全然合わないけど、構造力学...
***供試体2(7層8プライ/接続金物無し/ゴム有り)の静的載荷 [#...
-&link(clt1g.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
-&link(clt8g.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--ゴムのヤング率は分からないが、作るだけ作っておいた。ち...
--ちなみに、*HYPERELASTICでゴムのような超弾性も扱えるが、...
***供試体2(7層8プライ/接続金物無し)の弾塑性解析 [#k1f144ec]
-&link(clt1so.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
-&link(clt8so.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--中立軸の上下で別々の弾塑性材料にしただけ。応力ひずみ関...
---とはいえ、本とか色んなページ(&link(例えばこことか,http...
---しかも、応力-ひずみは、&link(この辺,http://www.str.ce....
-&link(clt1soko.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(clt8soko.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--中立軸の上下で別々の弾塑性材料にし、かつ支持部周辺を圧...
---藤田さんの計算結果が実験値よりも硬いのは、支持部周辺は...
-&link(plcltso.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--gnuplot用に荷重-変位関係をまとめる。
-ところで，入力データに特に書き込んでいないが，デフォルト...
--後で後藤さんから訊いた話では，apt-getでCalculiXをインス...
***供試体1(7層8プライ/接続金物有り/ゴム有り)の静的載荷 [#...
-&link(clt1jg.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
-&link(clt8jg.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--鋼材は板要素(S8)として木材と節点共有させたので、若干寸...
***供試体3(7層8プライ/供試体2の半分)の静的載荷 [#t2782c36]
-&link(clt1han.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(clt8han.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--図面をちゃんと見ていないので、寸法は正確でないと思うけ...
***$\ell$=2010mm,h=240mm(7層8プライ)のときの1層モデルと8...
-&link(clthaba.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--H鋼の間隔($\ell$=2010),床版厚(h=240mm(7層8プライ))を一...
--載荷は、床版中央に500x200mm(の1/4)の大きさを面載荷。
--境界条件は、面を(ボルト～H鋼の端)高さ方向と、横方向に固...
-&link(clthababun.f90,http://www.str.ce.akita-u.ac.jp/~go...
--床版の設計上で知っておきたいのは中央集中荷重よりも分布...
-&link(pclthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(clthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
---分布荷重についても同様に中央のたわみを見たが、分布荷重...
--たわみの収束性の確認をするのが面倒だったので、要素分割...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/clt...
-8層ものを1層のモデルに曲げ剛性から近似した場合、梁では両...
--あと、変位の相対誤差の大小に関わらず、1層モデルのミーゼ...
-それから、境界条件をCLT床版に近くして計算したが、より一...
--それなら、立体要素の8層モデルと1層モデル、板要素の1層モ...
***1層モデルと8層モデルのたわみの誤差についての考察 [#k00...
支間方向にz軸，橋軸方向にx軸，載荷方向にy軸をとった（梁の...
$$\left(\begin{array}{c} \varepsilon_{x}\\ \varepsilon_{z...
この逆関係は，
$$\left( \begin{array}{c} \sigma_{x}\\ \sigma_{z} \end{ar...
となるので，&link(岩熊先生のページ,http://mechanics.civil...
$$ D_{x} \equiv \frac{E_{x}}{1-{\nu_{xz}}{\nu_{zx}}}I^{(p...
と書ける．ここで，$I^{(p)}$は単位幅の断面二次モーメント（...
例えば，上述したように挽き板一枚の材料定数
$$E_z=5GPa, E_x=E_y=0.2GPa, \nu_{xy}=\nu_{xz}=\nu_{yx}=\n...
から，8層分を平均化した材料定数を
$$\overline{E_z}=3.5GPa, \overline{E_x}=1.7GPa, E_y=0.2GP...
というように決めた時，8層モデルの板の曲げ剛性は（$\sum_{i...
$$ D_{8層x} = \frac{\overline{E_{x}I^{(p)}}}{1-0.016\cdot...
1層モデルの板の曲げ剛性は
$$ D_{1層x} = \frac{\overline{E_{x}I^{(p)}}}{1-0.19428571...
ここで，1層モデルと8層モデルとで梁としてEIが同じになるよ...
とはいえ，CLTにKirchhoff-Loveの仮定を適用するのにもかなり...
***1層モデルのポアソン比について考えてみる [#b53c4aa6]
というわけで，8層モデルのたわみに合わせようとするなら，ど...
軸のとり方やら材料定数やら寸法やらを上述した7層8プライの...
$$ \nu_{zx}=\sqrt{0.0064\frac{E_{z}}{E_{x}}}=\sqrt{0.0064...
で，あまり変位に影響しなさそうなポアソン比を，仮に今求め...
になる．こんな決め方でいいのかどうか，幅/スパン比を変えて...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/clt...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/clt...
となり，以前のものと比べてかなり改善される．ちなみに今回...
-&link(clthabas.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(clthabab.f90,http://www.str.ce.akita-u.ac.jp/~goto...
-&link(pclthaba.f90,http://www.str.ce.akita-u.ac.jp/~goto...
***長方形板の4辺支持の式を計算 [#z08442e4]
-&link(4hen.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--等方性長方形板の4辺単純支持。構造力学公式集より。
--2辺単純支持とかも作っていたが、たわみが発散する。式が複...
***供試体2(7層8プライ/ボルト有り)の1/4解析 [#b2859c74]
-&link(cltbyon.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--木材利用研究論文報告集15のpp.17-22の有山さんのを真似る...
--実験では降伏するまで載荷したものの、これといってCLTの損...
---計算してみたところ、随分とボルトの下部に応力が集中して...
--直方体(CLT)と円柱(ボルト)が組み合わさったメッシュを生成...
--メッシュをかなり細かくしないと、ボルトの断面が膨らむ(梁...
--静的試験について、エクセルのC列目には最大荷重が約10kNと...
**CLT橋 [#lb2cea18]
中村先生のCLT橋（桁に集成材，床版にCLTを使い，それらをラ...
-&link(nclt.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--3点曲げ載荷．単純に，CLTと集成材の境界面の節点を共有さ...
--境界の節点を共有させると，実験値よりもかなり剛性が高く...
-&link(ncltj.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--3点曲げ載荷．CLTと集成材の間に，接合要素みたいなものを...
--接合要素のせん断弾性係数の値によって，計算結果が大きく...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/ncl...
--上がCLT，下が集成材．150kN載荷時の変形を1倍で描くと，滑...
-&link(ncltj_bun.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(ncltj.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--A活荷重（13.5kN/m$^2$）が作用時のたわみは，中村先生の計...
---ちなみに，『木材利用研究発表会講演概要集16』のp.83の最...
-ncltjb.f90<--需要がなさそうなので，つくっていない
--ラグボルトを梁要素として，CLTと集成材の節点上に配置する...
*円筒 [#f07d6080]
**円筒(角筒) [#q6b64dd3]
円筒プログラムと板要素の練習。
-&link(entou.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--要素はS8.境界条件を色々と変えられるようにした。座屈モー...
-&link(hentou.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--初期不整を与えて、幾何学非線形解析をすると、途中で計算...
**折り目を剛結とした回転螺旋折り [#l2348b95]
-&link(gkaiten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--折り目は剛結。板一枚につき1要素なので、かなりメッシュが...
---&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--幾何学非線形解析にすると、高さ0.87hぐらいで計算が止まっ...
--段の高さに応じて高さ方向変位を多点拘束(*EQUATION)によっ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ent...
-&link(山内さんのモデルの作り方,http://www.str.ce.akita-u...
**折り目を剛結とした反転螺旋折り [#o3dd0190]
-&link(ghanten.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--折り目は剛結。板一枚につき1要素なので、かなりメッシュが...
---&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--幾何学非線形解析にすると、高さ0.85hぐらいで計算が止まっ...
--上下を板に挟まれた中の段は均等に畳まるが、片方が自由に...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/ent...
-線形解析をすると、回転螺旋折は上部が拡がってしまい、反転...
--&link(Salome-Mecaに詳しそうな人,http://salome-meca.coco...
--後に後藤さんから聞いた話によると、板要素の誤魔化してる(...
**板の要素数を増やす [#b84ef3d6]
-&link(s6x4.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--螺旋折りプログラムでは板一枚につき要素数1つしかできない...
--多分、S6要素のみのモデルなら、何にでも使えると思う。た...
--元々は、&link(構造解析関係のプログラム置き場,http://www...
-折り目を重節点にして&link(多点拘束,http://www.str.ce.aki...
-立体要素の方が(遠藤さんの結果によると)幾何学非線形解析も...
*鋼棒のねじれ座屈 [#ef3ac8e2]
**C3D8とC3D6の混在 [#ff707b9d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(nezid6.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--ねじって線形座屈解析。軸方向に2000分割ぐらいしないと、...
--石坂さんの($\ell$=3m,d=13mm,E=206GPa,$\nu$=0.3)を計算す...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(hnezid6.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--例の如く、初期不整を与えてねじる。初期不整をとても小さ...
**C3D8だけ [#cf73940e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/nez...
-&link(nezid8.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--初期不正を色々と工夫していると、どうやら初期不正が大き...
-&link(hnezid8.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--nonpositive jacobian は出なくなったが、途中で計算が止ま...
-気が向いたら、C3D20でもやってみたい。
**B32 [#c799c384]
-&link(nezib32.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--梁要素。$T_{座屈}$=863.1N・mで、立体要素よりも理論値(86...
---後藤さんが『オイラー角を用いた簡潔な有限変位解析手法』...
-&link(hnezib32.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--例のごとく、非線形座屈解析。座屈する瞬間は確認できるが...
*横ねじれ座屈の非線形解析 [#c3399edd]
-後藤ちゃちゃ(16/5/17)
--FEMで初期不正のある座屈シミュレーションから$P-\delta^{2...
--あと、ローマ数字や◯数字の漢字コードは、昔は機種依存文字...
---とりあえず、半角数字にしました。
--平沢先生の図面をk2のmokusinにアップしました。
-&link(katasen.f,http://www.str.ce.akita-u.ac.jp/~gotou/p...
--汎用有限要素ツールは、固有値解析で座屈を一気に求めてい...
--かといって、複雑な形状とかを立体要素にモデル化して、真...
--そこで、プリポストプロセッサーと連携しやすく、立体要素...
-過去の横ねじれ座屈の卒論等を参考にした。
--寸法：b=2cm,h=25cm,$\ell$=2.0m.
---スパンは一つしか試していない
--分割：幅方向8,高さ方向30,軸方向200.
---過去の卒論等の"幅方向6,高さ方向20,軸方向100"だと、線形...
--材料定数
---等方性:E=206GPa,$\nu$=0.4 .
**手法 [#qf64594b]
-1.CalculiXで線形座屈解析をする。
--&link(yoko.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
-2.1/5000(例えば)の初期不整を与える。
--初期不整の大きさはそれなりの工夫が必要。
-3.CalculiXで幾何学非線形解析をする。
--&link(hyoko.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---gifのファイルサイズが大きいので一回消したのだが，遠藤...
---&link(cgxppm.f,http://www.str.ce.akita-u.ac.jp/~gotou/...
---1~5次モードでもいかにもそれらしい挙動を示す。異方性を...
-4.座屈荷重を&link(P-$\delta^{2}$法とか,http://www.str.ce...
--&link(pd2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan...
--&link(sw.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--&link(meck.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--P-$\delta^2$法：終わりのの17~5点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
--Southwell法：初めの76点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---終わりの方のプロットも使うと、座屈荷重が少し大きくなっ...
--Meck法：初めの76点を使って、線形回帰。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/kondo/yok...
---終わりの方のプロットも使うと、座屈荷重が少し大きくなっ...
,, 固有値 , $P-\delta^2$法 , Southwell法 , Meck法 , Timos...
,座屈荷重(kN), 46.8, 46.9, 47.3, 47.3, 45.3, 45.5
--&link(理論値,http://www.str.ce.akita-u.ac.jp/~gotou/zik...
---計算には&link(trahair.f90,http://www.str.ce.akita-u.ac...
-なお、結果は省くが、オイラー座屈なども同様なかんじ。ただ...
*伊藤さんのトラス橋 [#xd333fa9]
-金具詳細モデル(鋼材:E=206GPa,$\nu$=0.3)を2次要素にして、...
, C3D4(Code_Aster) , C3D10(Code_Aster&CalculiX) , 梁要素?...
, 0.239mm, 0.242mm, 0.228mm, 0.277mm-->0.275mm , 0.250mm
--要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v...
--異方性によるたわみの相対誤差(CalculiX(多分間違えている)...
---初めは斜め部材の角度を60度にして計算したが、ちゃんと61...
--異方性によるたわみの相対誤差(Code_Aster(あっていそう)):...
-ついでに、樹脂材料(E=2.84GPa,$\nu$=0.3)でもやってみた。
, C3D4(Code_Aster) , C3D10(Code_Aster&CalculiX) , C3D10直...
, 0.586mm, 0.669mm, 0.775mm-->0.773mm, 0.699mm
--要素によるたわみの相対誤差:$\frac{v_{C3D4}-v_{C3D10}}{v...
--異方性によるたわみの相対誤差(CalculiX(多分間違えている)...
--異方性によるたわみの相対誤差(Code_Aster(あっていそう)):...
-ちなみに、2次元のトラス(&link(toraITO.f90,http://www.str...
, 3次元トラス(Marc) , 2次元トラス , 2次元骨組 , 2次元骨組...
, 0.228mm, 0.22749mm, 0.22677mm, 0.22679mm
--Marcを使った3次元トラスの計算は伊藤さんによるものなので...
--金具詳細モデルの計算ほどせん断の影響が大きくない。不思...
-(2017/7)Code_Asterでも直交異方性材料として計算したところ...
--伊藤さんのトラス橋をやりたい人のために，斜めの直交異方...
-ただ，読み取っている箇所(中央の金具上の角材同士の中心線...
--確か，Marcのトラス要素と比較するために，そのような載荷...
*CalculiX関連 [#u3310209]
**unv-->inp [#g38f0dd0]
-&link(unvc3d10.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(unvc3d4.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--ついでに、節点群(*NSET)を座標指定でなくて、Salomeでつく...
-&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~got...
--&link(unvc3d10.f90,http://www.str.ce.akita-u.ac.jp/~got...
--局所座標系の設定に不確かなところ(SYSTEM=CYLINDRICAL)が...
---jacobianエラー：輝さんの日誌によると、節点座標の出力形...
-&link(menc3d10.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~go...
---ccx_2.9だと問題ないが、ccx_2.6.1だとfrdファイルの要素...
---他にも、モデルによっては面情報を上手く変換できていない...
-&link(zac3d10.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--&link(mokuc3d10.f90,http://www.str.ce.akita-u.ac.jp/~go...
-&link(s6.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
-&link(s3.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
--Salomeで切った3角形の板要素(unv)をCalculiX用入力データ(...
---面内曲げの変位にS3とS6とでさほど結果に差はないのだが，...
**線形座屈解析から初期たわみを与える [#j6172625]
-&link(syokitawami.f90,http://www.str.ce.akita-u.ac.jp/~g...
--「線形座屈解析用のinpファイル」と「全節点変位を出力した...
--何次の座屈モードであっても、どんな要素であっても、何種...
---ただ、inpファイルからコピー＆ペーストした部分の一列目...
-Code_Asterでも，Medでなく，IDEAS(unv)などのテキストファ...
**frd-->vtk [#r7d8378c]
-&link(c3d10vtk.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--c3d10のfrdファイルをparaview用にvtkファイルに変換する。...
--frdファイルには書かれていない主応力を計算してvtkファイ...
-&link(c3d8vtk.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--よく使うC3D8でも同じもの作る．
-s6vtk.f90を作ろうとしたのだが，S6はfrdファイルにおいて1...
--というか，vtkフォーマット以外を採用するのもいいかも．
-paraviewに計算結果を入れられると，こんな↓遊び？(CLTの最...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
*paraviewメモ [#wff29fdc]
主に個人の好みや拘わりで変更できるところを記す．
**コンター図関連(Color Map Editor) [#zd684856]
-ハートマーク(Choose preset)：初期設定の(青<-->赤)から虹...
--個人的には，デフォルトの Cool to Warm が大変見づらく、...
---&link(ここ,http://www.str.ce.akita-u.ac.jp/~gotou/prog...
-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 の数...
**あるグループのみ表示 [#cbdee9c6]
-Filter -> Alphabetical -> Extract Group -> 見たいグルー...
**csvに出力 [#z4be0f93]
Extract Selection をしないと，どこかの節点を飛ばされるら...
-(Filters)Slice(あるいは Clip ) -> Select Points On(ある...
--Select Points On : 表面の節点のみを選択
--Select Points Through : 内部の節点も選択
--Extract Group と組み合わせて(作ったグループだけ表示させ...
**動画化 [#g65c47a7]
-(File)Save Animation
--細かな設定は試していないから知らない．
-その後，ImageMagickのconvertを使う．
*Salome-Meca関連 [#s486f8e6]
**ある断面の軸方向応力分布を描画 [#m68c736b]
-&link(danmen.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--paravisから出力したデータをgnuplot用に整える。ついでに...
--ちなみに、断面をグルーピングし忘れたときでも、以下の方...
---パイプブラウザーのrmedファイルを選択したまま、フィルタ...
---見たい断面の、原点座標(origin)を入力し、面の法線(nomal...
-&link(sankaku.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--$\sigma=0$の線とか$\sigma=\frac{M}{I}y$の線とかをgnuplo...
-&link(szz.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--'MED'ではなく'RESULTAT'に出力したデータを、gnuplot用に...
*Maximaリンク [#k95ca0eb]
とても便利な数式処理のソフトウェア．terminal上で動くmaxim...
-斉木先生のページ
--ttp://mechanics.civil.tohoku.ac.jp/saiki/strcanalysis/k...
*ImageMagickメモ [#o77bc022]
Windowsの場合は、多分&link(ImageMagickの公式ページ,https:...
**複数のスクリーンショットを一括でppm(アスキー)に変換 [#t...
--(
mogrify -format ppm -compress none *png
--)
-&link(cgxppm.f,http://www.str.ce.akita-u.ac.jp/~gotou/pr...
**定点観測で撮影した連続写真を重ね合わせ [#nedf90c1]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kondo/hok...
--(
convert *.JPG -evaluate-sequence max gousei.png
--)
-場合によっては、コマンド中のmaxをminにした方がいい。
-特に意味はないが、上のコマンドではデジカメの画像形式から...
-この画像をImageJ(Windowsの場合は&link(公式ページ,https:/...
--変位を測りたい場合は、
---Straght:距離が既知な線に沿って、線を引く。(ドラッグ&ド...
---Analyze-->Measure:線の長さをpixelで測る。
---Analyze-->Set Scale-->Known distance:既知な線の実際の...
---Analyze-->Measure:線の長さを実際の長さで測る。
--たわみ角(傾き)を測りたい場合は、
---Angle tool:測りたい角度をなす3点を順にクリック
---Analyze-->Measure
-よく使いそうな保存の形式は、
--File-->Save as-->Results:結果をテキストで出力。
--File-->Save as-->XY Coordinates:現在、表示されている線...
-野田さんと3年生のたわみ試験で試してみたところ、精度がと...
-原因としては、
--1.三脚のカメラを固定するところが壊れているので、シャッ...
---これは、ちゃんとした三脚を使えばいいだけの話。
--2.被写体と1m近く離したのだが、たわみが小さい(8mm程度)こ...
--3.なるべく被写体の真正面にカメラを据え付けたつもりだが...
---原因3と原因2は、どちらかを優先させると、もう片方が疎か...
-同時に&link(Regard3D,http://www.regard3d.org/)というWind...
--追記：Ubuntu18.04LTSでは&link(COLMAP,https://colmap.git...
--特徴点がたくさん抽出できる場合には、 1.Regard3D(SfM + ...
---ネット上にあったカエルのぬいぐるみや猿のぬいぐるみ(&li...
--ただし、百目石橋でもそれをやろうとすると、相当いいPC(ゲ...
-野田さんによると，tiffだと拡大してもぼやけにくいとか．例...
-convert -threshold 50% -edge 1 hoge.JPEG hoge.tif
-などとすると，輝度50%を境(50%でなくてもいい)に2値化し，1...
**pdfを変換するときに画質を下げない [#u11b1304]
pdfを画像に変換したい時，ただ，"convert hoge.pdf hoge.png...
*$\LaTeX$メモ [#v3c82d4c]
**minipageのキャプション [#ff759cff]
画像をminipageで横並びさせると，互いのキャプションが窮屈...
**単位 [#p7961415]
無頓着なもので，単位が斜体だとか間に半角スペースが入って...
--(
$\frac{P}{A} = \SI{100}{N/mm^2}$
--)
texlive-scienceは，"sudo apt-get install texlive-science"...
他にも色々と便利なマクロが梱包されていると思う．
**ダブルクォーテーション [#h41d44c1]
ずっと知らなかったが，"hoge"だと残念なことになる，正しくは
--(
``hoge''
--)
**画像をわざと表示させない [#s5197781]
最終版まで画像は表示される必要はなく，画像の範囲だけ調整...
--(
\documentclass[a4j,12pt,dvipdfmx,draft]{jarticle}
--)
というように，draftを入れてやると，画像は表示されず，画像...
**$\LaTeX$周りのツール群 [#tfc9f21c]
身の周りの人があまり使っていない？けど，とても便利なツー...
***pdfcrop [#o9674473]
TeXliveに勝手についてくるtexlive-extra-utilsに入っているp...
-余白無し
--(
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-cro...
***Ghostscriptによるpdfのグレースケール化 [#u5e5e5a5]
白黒で印刷をする前に，カラーで作ったpdfがグレースケール化...
"convert -type grayscale -density 300 hoge.pdf hoge_g.pdf...
--(
gs -sOutputFile=hoge_g.pdf -sDEVICE=pdfwrite -sColorConve...
--)
-ここで，1つ目のオプションの"hoge_g.pdf"は出力ファイルな...
-また，"?"はwikiの機能で勝手に付いてしまっているだけなの...
-他の細かなオプションについては分かりません．
-参考URL：
--http://takuya-1st.hatenablog.jp/entry/2014/09/22/010805
***Inkscape [#nf9f01af]
ドロー系の描画ソフト。Xfigの上位互換でgnuplot→Inkscape→$\...
コマンドラインで形式変換が行えるので、シェルスクリプトに...
**おまけ（修論の背表紙） [#z0b8365b]
修論本体をpdfで提出するようになったものの，依然として修論...
-&link(se.tex,http://www.str.ce.akita-u.ac.jp/~gotouhan/j...
--用紙の向きが縦書き．余白や空白を適宜調整する必要あり．
-&link(sel.tex,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
--用紙の向きが横向き．タイトルが短いため横向きにしてみた...
---&link(例,http://www.str.ce.akita-u.ac.jp/~gotouhan/j20...
ちなみにfirefoxだと，縦書きの「ー」は横向きのままのようだ...
*その他，便利なコマンドライン [#ta1d4bec]
**pdfimages [#x7d58f6a]
pdfから画像を抜き取るのに使える．例えば，pngとして画像を...
--(
pdfimages -png hoge.pdf out
--)
**grep [#p3818bef]
***特定の文字を含むファイルを検出 [#dc7039b5]
--(
grep ほげ -rl .
--)
これだと今作業しているところ(.)にあるファイルで「ほげ」と...
勿論，カレントディレクトリ以外のパスでもいいし，オプショ...
*その他、便利なソフトウェア [#a143fc30]
**[[Hugin:http://hugin.sourceforge.net/]] [#aa12ccdf]
パノラマ写真をつくるためのソフトウェア。aptでインストール...
**[[Resynthesizer:http://www.logarithmic.net/pfh/resynthe...
GIMPでスマート消去するためのプラグイン。
-インストール：https://github.com/bootchk/resynthesizer
-使い方：自由選択 → フィルター → 強調 → スマート消去 → 半...
*3Dプリンター メモ [#f2cd76b0]
-ボールジョイントも作れそう(PRN3DでなくてEDENでの話)。
--http://rayerd.sblo.jp/article/172174211.html
---造形後に関節を折り曲げることで、隙間につまったサポート...
--http://rayerd.sblo.jp/article/172249880.html
--http://rayerd.sblo.jp/article/172529855.html
*gnuplotメモ [#id35cce2]
$\LaTeX$でグラフを描く時は，gnuplotとxfigを使うのよりも，...
**日本語のpngを出力する [#x8e397c5]
ロードするファイルの頭（set out "hoge.png"の上）に
--(
set term png font "/usr/share/fonts/truetype/takao-gothic...
--)
を入れる（sizeは好きな値にする．特に指定しない場合は入れ...
--(
nkf -guess hoge.plt
--)
でShift-JISですよ，などといわれた場合には
--(
nkf -w-overwrite hoge.plt
--)
でutf-8に変換する必要がある．
*Fortran90(gfortran)メモ [#ee62b9d7]
**固有値計算 [#z320f171]
-&link(jacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotouh...
--ヤコビ法によって実対称行列の固有値と固有ベクトルを計算...
---frdファイルをvtkファイルに変換しようとして，frdファイ...
---scipyやoctave,julia等を使えば一行で簡単に固有値と固有...
-&link(sjacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--&link(jacobi.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
**LAPACK [#i3976944]
学習・自作が難しそうな大規模な行列演算が簡単にできる．こ...
-Ubuntu16.04LTSでのインストール方法
--https://qiita.com/AnchorBlues/items/69c1744de818b5e045ab
--全くこの通りで整備できた(LAPACKとLAPACK95)．
-コンパイルの仕方やお試し計算
--http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux...
**FFTW [#a97b2f39]
高速フーリエ変換などが簡単にできる．
野田班でフリーのfftアナライザを使っていたが，自分で作って...
&link(梅原さん,http://www.str.ce.akita-u.ac.jp/~gotou/pro...
-Ubuntu16.04LTSでのインストール方法
--http://moznion.hatenadiary.com/entry/20120718/1342626456
--全くこの通りで整備できた(リンク先はUbuntu12.04だけど)．
-コンパイルの仕方やお試し計算
--http://www.rcs.arch.t.u-tokyo.ac.jp/kusuhara/tips/linux...
*Julia メモ [#adc05a47]
マトリクス構造での導入案があるとかないとか？
インタープリタ言語の割に実行速度がとても早く，;でシェルモ...
強いて言えば，（グラフ描きツールなど）色々と選択肢があり...
**Ubuntu上で最新版にする [#j30a97d0]
--(
sudo add-apt-repository ppa:staticfloat/juliareleases
sudo add-apt-repository ppa:staticfloat/julia-deps
sudo apt-get update
sudo apt-get install julia
--)
既にJuliaをインストール済みの場合は，
--(
sudo apt-get autoremove julia
--)
としてから上記の手順を踏まないと，駄目でした．
**パッケージのインストール [#ced80124]
例えば，グラフ描画ツールなら，juliaの対話モードを立ち上げ...
--(
Pkg.add("Plots")
--)
と打つだけでインストールできる．そして，使うときは using ...
-ちなみに，ここでの例のPlotsは何だかplotly.jsと同じみたい...
-JuliaはPGFPlotsにも対応しているようなので，やはり論文な...
**JuliaFEM [#aaa26ba5]
こんなのもある。
-http://www.juliafem.org/
*Python メモ [#sa31cbf5]
**コメント文にutf-8の日本語を使えるようにする [#t225c8d5]
ヘッダーにこれをいれる。ただし、Python3系では入れなくてよ...
--(
# -*- coding: utf-8 -*-
--)
**Matplotlibで日本語を使えるようにする [#jf2ccabd]
端末にて、
--(
mv ~/.cache/matplotlib/fontList.cache ~/.cache/matplotlib...
sudo cp /usr/share/fonts/truetype/takao-gothic/TakaoPGoth...
sudo vi /etc/matplotlibrc
--)
と入力した後、matplotlibrc内の"#font.family : sans-serif"...
*Rustメモ [#sfcec57b]
初めて聞く概念ばかりで難しい…．そろそろLinuxもC以外の言語...
**rust.vim [#y876ed73]
まだ若い言語なので（Ubuntu16.04での）デフォルトのvimでシ...
-https://github.com/rust-lang/rust.vim/blob/master/syntax...
-を"/usr/share/vim/vim74/syntax/"にダウンロードする．
-次に"sudo vi /usr/share/vim/vim74/filetype.vim"で，どこ...
--(
" Rust
au BufNewFile,BufRead *.rs setf rust
--)
-を追加し，拡張子が被っているHerculesをいじってごまかす(...
--(
" Hercules
" au BufNewFile,BufRead *.vc,*.ev,*.rs,*.sum,*.errsum s...
au BufNewFile,BufRead *.vc,*.ev,*.sum,*.errsum setf herc...
--)
-どうやらvimには拡張子(言語)別に辞書機能のようなものや自...
*Goメモ [#k38b93fc]
シンプルで早い．試していないけど，クロスコンパイルできる...
**パッケージのインストール [#b7e8fd15]
標準ライブラリだけでも物凄い数があるが，
それ以外にもオープンソースのライブラリが沢山公開されてい...
それらを入手するときの注意点として，golangだけでなく''git...
--(
go get github.com/〇〇
--)
-"github.com/○○"はURL．
-数値計算系ライブラリの有名どころを探してみると，それなり...
--BLAS/LAPACK:github.com/gonum/gonum
--FFTW:github.com/runningwild/go-fftw
-ただし，これは&link(この方法,http://golang-jp.org/doc/in...
**フォーマッタ [#qb158e8f]
インデントやスペース数が揃ってなくてもコマンド一発で綺麗...
--(
go fmt hoge.go
--)
*Haskell(GHC)メモ [#y3b3b7bc]
純粋関数型言語とやらを触ってみたい．
**パッケージのインストール [#abab45ce]
例えば，行列演算ライブラリとかは，
--(
cabal update
cabal install hmatrix
--)
gnuplotも&link(呼び出せる,http://d.hatena.ne.jp/mmitou/20...
*METAFONTメモ [#c59d00e6]
同研究室の某学生がM冨フォントを使いたいと話していたので，...
モデルとなる字の入手が困難を極めるため，道のりは険しい……．
**gftodvi [#efee22ee]
文字の確認方法．
-mf hoge.mf と打つと，hoge.2602gf というファイルができる．
--METAFONTについては，クヌース先生が書いた本（日本語訳）...
--他の参考ページ：ttp://hayato1224.web.fc2.com/MFmain.pdf
-gfファイルをdviファイルに変換するために，gftodviというコ...
--(
gftodvi: fatal: tfm file `gray.tfm' not found.
--)
-というエラーが出てしまい，生成されたdviファイルは0バイト...
--(
mktextfm gray
--)
-と打っておくと，gftodvi hoge.2602gf でgfファイルをdviフ...
-あとは，xdvi hoge.dvi& でも開けるし，dvipdfmx hoge.dvi ...
**mktextfm [#ye59d54e]
tfmという形式にすると自作フォントを$\LaTeX$で使うことがで...
--(
mktextfm hoge
--)
しかし，普通の文字や数式の感覚で使えない？ので，使い勝手...
*Vimメモ [#h3f9ed74]
**ある単語の数を数える [#o7aeff5a]
--(
:%s/hoge/&/gn
--)
**シンタックスハイライトの追加 [#l9a91177]
デフォルトでかなり沢山の言語が，勝手にシンタックスハイラ...
-"/usr/share/vim/vim74/syntax"にGitHubなどからhoge.vimを...
-"/usr/share/vim/vim74/filetype.vim"を編集し，拡張子から...
--もしかすると，同じ拡張子を使っている言語もあるかもしれ...
*USBメモリにバックアップ [#f293b776]
--(
rsync -auv /home/kouzou/2016 /media/disk
--)
-オプション"-auv"は，「a:パーミッションを維持，u:差分のみ...
-"disk"は、USBメモリに付けた名前。
-ちなみに，k2へのバックアップ時にも使える(パスを間違えそ...
*HHKB Pro2 [#sb4fe682]
HHKB Pro2 をただPCに繋いだだけだと，\や_が打てなくて困る...
&link(この配列,https://www.pfu.fujitsu.com/hhkeyboard/lea...
-システム -> 設定 -> ハードウェア -> キーボード -> レイア...
-右上のキーボードマークを右クリック -> 設定 -> + -> US ->...
間違いなく打ちやすいが，日本語キーボードを普段使っている...
とはいえ，じっと眺めているとUS配列の方が記号の位置がきれ...
*講義関連 [#edbe1bee]
**LibreOffice [#e18c2630]
-https://help.libreoffice.org/Main_Page/ja
**構造力学特論 [#befb0708]
-&link(torasu2.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(hone2.f90,http://www.str.ce.akita-u.ac.jp/~gotouha...
--格点を剛結にするとどのくらい曲げが入るのか気になったの...
-&link(toraITO.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
-&link(honeITO.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--伊藤さんの木製トラス橋を２次元の問題として解いてみたか...
-&link(honetimo.f90,http://www.str.ce.akita-u.ac.jp/~goto...
--さらに、伊藤さんの木製トラス橋の2次元骨組を剛性行列をテ...
**土質工学特論 [#f7040b16]
-間隙水圧分布を計算:&link(u.f90,http://www.str.ce.akita-u...
--set zrange [] reverse : z軸を下向き正にする。
--set pm3d : z軸の値を色付けする。x軸でこれをやりたい。調...
*TA関連 [#v35f5e11]
**基礎物理学実験 [#ua539bd7]
-&link(指導書,http://sht.phys.akita-u.ac.jp/~buturi/)
-&link(基礎物理学実験TA研修システム,http://sht.phys.akita...
--Ubuntu12.04のFirefoxでWindows Media Playerを再生するに...
---gnome-mplayer
---gecko-mediaplayerSynaptic
--をパッケージマネージャーでインストールすると、エラーは...
**土木環境工学実験(アクリルの引張試験) [#aff52f8f]
-ヤング率とポアソン比の確認:&link(yangu.f90,http://www.st...
--syokignu.fでデータ整形後
*時間割と単位 [#i354bb5c]
, 授業科目区分 , 修了に必要な単位 , 取得単位数 , 2017年度...
,,,,,
, 共通科目(必修) , 2単位 , 2単位,,
, 専門科目(必修) , 13単位 , 1単位, 12単位,
, 共通科目および専門科目(選択) , 10単位以上 , 12単位,, (...
, 下の3つを合わせて15単位以上, 2単位以上 , 2単位,, (他専...
,,, 4単位,, (他共通)
-「関連性のあるコースで開講する専門科目」とは：理工学研究...
--しかし，2017年度入学の人は，ガイダンスで，システムデザ...
-2017年度に取得予定の残り12単位の内訳は、
--システムデザイン工学課題研究：専門必修(2年間)=10単位
--システムデザイン工学演習：専門必修(通年)=2単位
**2016年度 [#m81db7e0]
-理工学デザイン：1共通必修(後期前半)
-システムデザイン工学課題研究：5専門必修(通年)
***2016年度前期 [#y7eed7a5]
,, 月 , 火 , 水 , 木 , 金
,,,,,,
, 10:20 , 構造力学特論 , 都市システム計画特論 ,,, 材料設...
, ｜ , 218 , 総合研究棟7階 ,,, c317
, 12:00 , 2専門選択 , 2専門選択 ,,, 2専門選択
,,,,,,
, 12:50 , , 土質工学特論 , 基礎物理学実験B , 土木環境工学...
, ｜ ,, c317 , 教育4号館104 , 218 ,
, 14:20 ,, 2専門選択 , (TA) , (TA) ,
,,,,,,
, 14:30 , ゼミ ,, 基礎物理学実験B , 英語ゼミ , 地域防災学...
, ｜ , 218 ,, 教育4号館104 , 218 , 共320
, 16:00 , (おまけの人) ,, (TA) , (おまけの人) , 2共通選択
***2016年度後期 [#n42ca291]
,, 月 , 火 , 水 , 木 , 金
, 12:50 , 基礎物理学実験D , 構造設計学特論 , 木構造工学 ,...
, ｜ , 教育4号館104 , c317 , c317 , c319 , 218
, 14:20 , (TA) , 2専門選択 , 2専門選択 , 1専門必修 , (お...
,,,,,,
, 14:30 , 基礎物理学実験D ,,地域産業アントレプレナー論(前...
, ｜ , 教育4号館104 ,, 共127 ,, 共224
, 16:00 , (TA) ,, 1共通必修 ,, 2共通選択
,,,,,,
, 16:10 , 地域産業論 ,,,,
, ｜ , 2号館3階学習支援室 ,,,,
, 17:40 , 2他専門選択 ,,,,
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
