斉藤輝の修論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*図などはsyuronフォルダの中にあります [#eed59bd8]
修論本体に使った座屈モードなどの
図は2014i5/tex/syuron/ にだいたい置いている。
*tugite130710鋼材荷重-ひずみ関係 [#p6cccb26]
tugite130710.xlsxのデータをゼロ点修正
座屈しそうな中央のトラス部分にひずみゲージをつけて破壊試...
ch33～36の斜材・垂直材ともに引張りのまま破壊している。
支間中央部は継手で補剛されているので心配なさそう。
縦軸 荷重(kN)、横軸 ひずみ(μ)
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
**予備試験 [#g8532b48]
tugite130710yのデータをゼロ点修正
縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
**両方プロット [#rba2340d]
tugite130710.xlsx、tugite130710yのデータをゼロ点修正
縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*目標 [#o5a46862]
修士発表時間
http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%c...
-等圧縮を受ける板の座屈を確認
*板の座屈 [#h2a0e021]
**等圧縮 [#m6184594]
***正方形板 4辺単純支持 [#j5c4a131]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
&link(kouitass8c.f90,http://www.str.ce.akita-u.ac.jp/~got...
$\sigma_{cr}=k\frac{D\pi^2}{hl^2}$
$N_{cr}=k\frac{D\pi^2}{l^2}$ よって
$N_{cr}=h \sigma_{cr}$
ここで、$\sigma_{cr}$は座屈応力、$N_{cr}$は座屈荷重、
Dは板の曲げ剛性、hは板厚、lは基準にとった板の辺長、
正方形板であり、4辺単純支持なので、k=4
$N_{cr}=0.009 \times 241.2955$
$ =2.17166 [MN] $
,座屈モード,ccx座屈荷重(MN) ,理論値との誤差(%)
, 1 ,0.2185645E+01,0.643977
, 2 ,0.2512923E+01
, 3 ,0.6733093E+01
, 4 ,0.7050507E+01
, 5 ,0.7651342E+01
板の厚さ方向より幅方向高さ方向分割の方が解に与える影響が...
幅方向高さ方向をもっと細かくすると理論値より若干小さい値...
-1次座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
***長方形板 4辺単純支持 [#s9e4b45f]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
&link(kouitassk9.f90,http://www.str.ce.akita-u.ac.jp/~got...
長方形板だけどk=4になる様なa,bの比なので、正方形板と理論...
$N_{cr}=0.009 \times 241.2955$
$ =2.17166 [MN] $
,座屈モード,ccx座屈荷重(MN) ,理論値との誤差(%)
, 1 , 0.2147385E+01 ,1.13735
, 2 , 0.2163955E+01 ,
, 3 , 0.3919156E+01,
, 4 , 0.5249193E+01,
, 5 , 0.5451060E+01,
-1次座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**純曲げ [#pa5aaff7]
三角形分布でも総荷重を分配で良さそう。
一つの三角形に1/4MNを分配してやる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**参考になりそうな文献 [#o47661f2]
-長方形板の分岐階層構造の境界条件による退化
https://www.jstage.jst.go.jp/article/jscej1984/1995/507/1...
-弾性矩形板の曲げ解析手法の一般化について
http://ci.nii.ac.jp/els/110007887197.pdf?id=ART0009558779...
-等曲げを受ける弾性矩形板に生じる二次座屈現象の解明
https://www.jstage.jst.go.jp/article/jscej1984/1995/519/1...
*Salomeで作ったモデルで座屈 [#e2f56b5a]
**応力集中部を細かく切ったモデル [#od3e26fe]
非一体化にして座屈させるところの鋼板を多く切る。
応力集中部を細かくしていくと1次モードは飛び出しがないモー...
鋼板要素数1は非一体化にしてるところと一体化しているところ...
鋼板要素数2は残りの鋼板全部。
-孔なし補剛材なし
,総節点数,木材要素数,非一体化鋼板要素数,鋼板要素数1,鋼板...
,631768,773611(l12),1389462(l1.6),347918(l2),291211(l8),2...
**補剛材なし [#e1ae5037]
-孔なし補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,389725,704146(l12),650860(l2),298061(l8),1次モードは正常...
,535235,601922(l12),1389462(l1.6),331039(l8),1次モードは...
--プログラムで作ったモデルと近い値が出る。
--1次モードでは座屈荷重77.7kN、2次モードでは座屈荷重78.5k...
2次モードでは要素が飛び出す。
2次モードでは鋼板の幅方向の変位が1次モードより小さい。
-メモ
--l12-l1-l8だとエラーで切れない
--3孔鋼材非一体化部だけをl1で切ると要素数470万くらいにな...
--l10-l1-l8でもエラーで切れない、要素数1000万くらいになり...
--l8-l1-l8もエラーで切れない
-孔あり補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,342007,599145(l12),577961(l2),231480(l8),1次モードは正常...
,522828,1017878(l12),1109490(l1.6),234967(l8),1次モードは...
--l12-l1.6-l8は切れる
--孔をつけたモデルでも解析がうまく行く。孔ありの方が1次の...
**補剛材あり [#yc7315fa]
-孔なし補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,882260,1798392(l12),1653936(l2),596385(l8),要素が飛び出...
-孔あり補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,523050,598470(l12),1251573(l2),452306(l8),要素が飛び出す...
--孔なしのほうが木材も鋼板も細かく切れているので1次座屈荷...
**今までのをやり直す。 [#c05c0b12]
salomeで木材部を削って一部非一体化にしたモデルを分割して...
-孔あり補剛材なしsalomeでモデリング
要素はC3D4(四面体)のAutoでlengthを変えて切る
,length,総節点数,木材要素数,鋼板要素数,結果
,12,17万,58万,15万,座屈させると要素が飛び出す。1次$P_{cr}...
,10,,,,なぜかエラーがでて切れない
,8,50万,203万,36万,要素の飛び出しがなくなる。1次$P_{cr}$=...
,7,581146,2396666,386169,要素の飛び出しが見られない。1次$...
,6,92万,340万,104万,要素数が多くなりすぎてccxで解けない(...
鋼板をもっと多く切ってみる必要がありそう
2次要素とかで切ってみるとどうか？
-孔なし補剛材なしsalomeでモデリング
,length,総節点数,木材要素数,鋼板要素数,結果
,12,181632,765824,178789,座屈させると要素の飛び出しが起こ...
,10,,,,なぜかエラーがでて切れない
,9,,,,エラーが出て切れない
,8,521420,2062159,376681,座屈させると要素の飛び出しが起こ...
,7,597914,2412348,450689,座屈させると要素の飛び出しが起こ...
*プログラムで作った孔なし非一体化モデル [#ud1a5503]
(鋼材と木材を支間中央部から最初のPC鋼棒のアンカープレート...
onsih4z.f90,onsi4z2.f90
を座屈させる。
1MNかけているので出力された座屈荷重がそのまま読める。
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,単純支持,3.701E-002
,固定支持,0.1480
**孔なし補剛材なし非一体化モデル [#x7a11643]
,ccx解析結果,
,n次モード,ccx座屈荷重(MN)
,1,0.7846011E-01
,2,0.2280486E+00
,3,0.2489195E+00
,4,0.4002716E+00
,5,0.4099577E+00
--オイラーの固定支持よりかなり低い値が出る。
--salomeで作ったモデルのような飛び出し部は確認されない。
--salomeで作った孔なしモデルは92kNくらいで座屈するが、そ...
-1次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-1次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-2次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-2次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**孔なし補剛材なし上の段だけ非一体化モデル [#x18a4a36]
下の段の木材と鋼板は一体化しているモデル
,n次モード,ccx座屈荷重(MN)
,1,0.8540278E-01
,2,0.2201420E+00
,3,0.3317080E+00
,4,0.4245944E+00
,5,0.4385326E+00
-1次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**孔なし補剛材あり非一体化モデル [#ca81a619]
--最初に木材部分の応力集中部(載荷点)が破壊し、
全てのモードで木材部の座屈挙動になるため、
鋼板部の座屈挙動が見ることができない。
--ためしに補剛材部分の寸法をものすごく薄くして補剛材なし...
断面にして解析したら上の補剛材なしプログラムで解いた結果...
,ccx解析結果,全て木材部の破壊
,n次モード,ccx座屈荷重(MN)
,1, 0.1028976E+00
,2, 0.1759758E+00
,3, 0.1906934E+00
,4, 0.2086057E+00
,5, 0.2153555E+00
**孔なし補剛材あり上の段だけ非一体化モデル [#g2e56c6b]
下の段の木材と鋼板は一体化しているモデル
--最初に木材部分の応力集中部(載荷点)が破壊し、
全てのモードで木材部の座屈挙動になる
,n次モード,ccx座屈荷重(MN)
,1,0.1049848E+00
,2,0.1809716E+00
,3,0.1925376E+00
,4,0.2122353E+00
,5,0.2144374E+00
*Asterのbuckling [#sa2492a9]
-解説pdf
http://framabook.org/docs/Code_Aster/beginning_with_Code_...
の15.2
ver10.8のcommファイルを使う。
これ見ながらcommファイルはオールグリーンにできたけど、pos...
たぶんやり方が違うのだろう。
&link(commファイル,http://www.str.ce.akita-u.ac.jp/~gotou...
-Asterのマニュアル(公式HP)
http://www.code-aster.org/V2/doc/default/en/index.php?man...
左上で言語をenglishにできる
-Bucklingについて
--http://www.code-aster.org/V2/doc/v11/en/man_u/u2/u2.08....
-MODE_ITER_SIMULTのマニュアル
--http://www.code-aster.org/V2/doc/default/en/man_u/u4/u4...
MODE_ITER_SIMULTのオプションを使いそう
--lateral buckling of a beam (spill)
http://www.code-aster.org/V2/doc/default/en/man_v/v2/v2.0...
どうやって解くかは書いてない
Aster→2材料で解くのと異方性はできる。座屈は調査中
-ccx2.7のマニュアル-*Friction
http://web.mit.edu/calculix_v2.7/CalculiX/ccx_2.7/doc/ccx...
-ccx2.5のマニュアル*Friction
http://its-cs1.its.uni-kassel.de/CalculiX-2.5/html/ccx/no...
*弾性床上の座屈 [#uf9a40dc]
http://freecaetester.blog62.fc2.com/category19-0.html
上の記事を参考にして一部非一体化にしたプログラムに接触解...
-接触入れるとWARNINGが出る。
--(
*WARNING in statics: a nonlinear geometric analysis is r...
but no time increment nor step is specified
the defaults (1,1) are used
--)
-接触解析で要素分割を細かくするとエラーになる
--(
*ERROR: increment size smaller than minimum
best solution and residuals are in the frd file
--)
マスターのヤング率を低くする(ezz/10とか)と何故か解ける
*座屈解析 [#p7cd23f7]
**鋼板のオイラー座屈 [#c1b15f5d]
プレストレス木箱桁橋 鋼板部プレストレス固定支持間315mmの...
--プレストレス鋼棒間730mmの半分365mmからアンカープレート1...
幅120mm高さ9mmスパン315mmで単純支持にしてオイラー座屈させ...
細長比121.24
オイラー座屈荷重理論値$Pcr=0.14937MN$
,ccx解析結果,,
,n次モード,ccx座屈荷重(MN),オイラー荷重理論値との相対誤差...
,1,0.1548903E+00,3.670
,2,0.6218848E+00,
,3,0.1393993E+01,
,4,0.1936985E+01,
,5,0.2383931E+01,
長方形断面だからちょっと誤差が出る。
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-スパンを630mmにしたやつ
,n次モード,ccx座屈荷重(MN),オイラー荷重理論値との相対誤差...
,1,0.1472202E+00,1.439
,2,0.3013863E+00
,3,0.5895368E+00
,4,0.8915795E+00
,5,0.1322299E+01
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**木材と鋼板1枚 上一段だけモデル [#tcbb0247]
木材と鋼板1枚で、支間中央部のPC鋼棒間だけが木材と鋼板が非...
,座屈モード,座屈荷重(MN)
,1, 0.1349740E+00
,2, 0.4743007E+00
,3, 0.5364287E+00
,4, 0.9919370E+00
,5, 0.1180486E+01
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
接触解析にしてないのでオイラー座屈、ccxで鋼板1枚を座屈さ...
**オイラー座屈応力 [#a3068bda]
-&link(zakutu,http://www.str.ce.akita-u.ac.jp/~gotou/pas/...
**1/4解析モデル [#y4702bc7]
まず普通の1/4解析モデルを作る。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
どちらも下の 半解析-箱断面三角孔モデル のx変位、y変位と同...
たわみ平均なども大体同じ値になる。
*鋼板の座屈 - 一体化だけど一部だけ鋼板と木材の間を削った...
とりえあず、中央部の木材と鋼板がくっついている部分の、木...
実物は上弦材も下側木材もPC鋼棒間で非一体化なので、FEMモデ...
**三角孔あり補剛材なし [#df7f9762]
-面載荷しているモデル。length12。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-木材部分に線載荷したモデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
***座屈荷重(1MN載荷時)とアニメーション [#odb2bae9]
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,単純支持,3.701E-002
,固定支持,0.1480
-ccxのBuckle
,モード,座屈荷重(MN),,
,,線載荷,面載荷(載荷試験装置寸法で載荷)
,1, 0.9249756E-01,0.9726263E-01
,2, 0.1472920E+00,0.1548828E+00
,3, 0.1699003E+00,0.1734040E+00
,4, 0.3865533E+00,0.3048128E+00
,5, 0.4932473E+00,0.3726462E+00
今回は固定支持であると考えられるため、ccxのモード2の計算...
-モード１
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
何かが飛び出す現象がみられる。
アニメーションで詳しく見てみる。
--モード1アニメーション、Amplitude+100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
--モード1アニメーション、Amplitude-100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
この2つのモード1アニメーションを見る限り木材部ではなく鋼...
-モード２
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**三角孔あり補剛材あり [#r1017a02]
***座屈荷重(1MN載荷時)とアニメーション [#o50a3b1a]
-ccxのBuckle(線載荷)
,モード,座屈荷重(MN)
,1, 0.9889383E-01
,2, 0.1574743E+00
,3, 0.3737713E+00
,4, 0.4575766E+00
,5, 0.5303813E+00
--モード1アニメーション、Amplitude+100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
・補剛材なしとありでは飛び出す方向が違う
→単に方向が違うだけなので補剛材のあるなしでモードが変わる...
・またはメッシュ分割の違い
（鋼板部メッシュ数は補剛あり18万、補剛なし15万でそんな変...
によって違うモードになったかも
**孔なし補剛材なし [#y0e7fbf9]
,モード,座屈荷重(MN)
, 1, 0.9200184E-01
, 2, 0.1465025E+00
, 3, 0.2212824E+00
, 4, 0.3723138E+00
, 5, 0.4482462E+00
4hanhakoananasi2
1次、2次モードの座屈荷重は鋼板に孔が無くても変わらない。
**孔なし補剛材あり [#rac8af2b]
,モード,座屈荷重(MN)
, 1, 0.1209061E+00
, 2, 0.1913690E+00
, 3, 0.3107388E+00
, 4, 0.4833085E+00
, 5, 0.5042332E+00
4hanhakoananasihogou
孔なしの場合、補剛材があると座屈に対してかなり補剛効果が...
**まとめ [#o07e4fc6]
オイラーの座屈荷重による計算結果と解析による補剛材なしの...
ccxによる解析結果のモード1は、断面に対し斜め方向からなど...
普通の弱軸よりもさらに弱い軸があるのではないかと考えられ...
どのモードを見ても解析結果では補剛材なしより補剛材ありの...
ちなみに鋼板と木材部を一体化させると先に上部木材部が破壊...
一体化させないモデルよりも高い座屈荷重が出る。
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,固定支持,0.1480
-ccxによる座屈解析結果
,モード,座屈荷重(MN),
,,三角孔あり補剛材なし,三角孔なし補剛材なし,三角孔あり補...
,1,0.9249756E-01,0.9200184E-01,0.9889383E-01,0.1209061E+00
,2,0.1472920E+00,0.1465025E+00,0.1574743E+00,0.1913690E+00
,3,0.1699003E+00,0.2212824E+00,0.3737713E+00,0.3107388E+00
,4,0.3865533E+00,0.3723138E+00,0.4575766E+00,0.4833085E+00
,5,0.4932473E+00,0.4482462E+00,0.5303813E+00,0.5042332E+00
*Aster [#me5f5ff7]
最終的な解析モデルをExplodeして
-荷重をかける面
-固定する面
を作っておく
**検証 [#t6e667dc]
簡単な片持ち梁を軸方向に押してみる。
-断面2cm×2cm、長さ40cm
-鋼材206GPA,ポアソン比0.3
asterはデフォルトの載荷方法だと面を指定してその面に対して...
その他の載荷方法は後でviでソース(.comm)をいじるか、
Aster画面でEficasを起動させてEficas経由で.commをいじるか...
,理論式,理論値(m)
,$d=\frac{PL}{EA}$,-4.85437E-07
,,変位最大値dmax(m),理論値との誤差(%)
,ccx節点載荷,-4.979490e-07,2.577%
,ccx面載荷,-4.846880e-07 ,0.1543%
,aster面載荷,-4.84655e-07,0.1611%
大体同じになる
**Asterで異方性 [#y23a9d1b]
ccx・理論値と比較し、異方性がちゃんと解けるか確認
***等方性で解く [#z4f5c2c1]
鋼材のたわみ(m)
節点載荷
,ccx,9.6939E-06
,aster,9.72292e-6(大体)
,初等梁,9.7087E-006
,ティモ,1.0146E-005
***直交異方性で解く [#z4f199fd]
-正方形断面の一辺20cm 長さ2mの木材
z軸方向ヤング率ezz=7.694GPa
exx=eyy=ezz/25d0
せん断弾性係数 一律でezz/15d0
片持ち梁とし、自由端に荷重100Nとする。
salomeではccxと同じように節点荷重として載荷する。
nx=ny=10,nz=100
,,たわみ最大値(m)
,初等梁理論値,2.5994E-004
,ティモ理論値,2.7164E-004
,ccx解析値,2.71617E-04
,aster,2.71617e-4
zがスパン(長さ)方向でモデリングしたなら、N(normal)に軸方...
Lに長手方向、Tに残りの1つでよいと思う。
当然せん断弾性係数やポアソン比も適切にいれる必要があると...
今回はccxのプログラムf90と同じになるように打ち込んだ。
-使ったファイル
--&link(hdfファイル,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(commファイル,http://www.str.ce.akita-u.ac.jp/~got...
寸法をm単位で入力してるので注意。
このやり方で木材を短くしたり長くしたりして解析したが、
ほぼすべてsalomeとccxの解が同じになった。
-&link(マニュアル,http://www.code-aster.org/V2/doc/v11/fr...
のP13,14にLTN系の座標と例が書いてある。
-Asterでの異方性の設定について書いてあるもの
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04....
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04....
--http://www.code-aster.org/forum2/viewtopic.php?id=4026
--http://www.code-aster.org/V2/doc/v11/fr/man_u/u4/u4.43....
のP13,14
**Commandファイル [#z0d028a9]
***荷重 [#oab18278]
-PRES_REP
デフォルトの設定。面に対して垂直な荷重をかける。
-FORCE_FACE
グループと方向を指定して荷重をかける。
*メモ [#x1bcae7a]
最近のモデル(スパン7m継手あり孔大きめ、補剛材ありなし)
-tugite1302～,1307～
*継手ありモデル [#tade4196]
-継手あり補剛材ありモデルの様子 継手部(一体化,半解析)
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
継手があると支間中央部の鋼板の三角孔が継手鋼板で塞がれる...
*木材利用研究発表会 [#qc0a738b]
2014年 &link(発表スライド,http://www.str.ce.akita-u.ac.jp...
*継手なしモデルまとめ [#j22fb978]
**3点曲げと半解析による比較 [#rd37bb00]
**木材に線載荷 [#v22e085e]
-ヤング率7.694GPa 荷重100kN
,モデル,解析手法,載荷上木材部たわみ平均値(mm),三角孔あり...
,補剛材あり三角孔ありlength25,3点曲げ,9.6324,,
,,半解析,10.2294
,補剛材あり三角孔なし,3点曲げ,8.99085,6.660%
,箱断面三角孔あり,3点曲げ,11.919
,,1/4解析,11.549
,箱断面長方形孔,3点曲げ,15.980
,箱断面三角孔なし,3点曲げ,10.9863,7.825%
,上部60mm鋼板欠け三角孔あり,3点曲げ,13.5324
**半解析と1/4解析で実験と同じように面に載荷 [#a9c7ca0a]
,モデル,解析手法,載荷上木材部たわみ平均(mm)
,補剛材あり三角孔あり,半解析,10.5003
,箱断面三角孔あり,半解析,12.4034
,,1/4解析,12.529
,上部60mm鋼板欠け三角孔あり,半解析,13.9751
-グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*3点曲げ木材上部分載荷 [#s252a518]
-全てC3D4要素
-支間中央木材部分に線載荷
-載荷上たわみ平均は、線載荷部からz軸方向に5mm離れたところ...
**補剛材ありモデル [#od9f0ef7]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
,載荷上たわみ平均(mm)
,9.98675
**箱断面モデル [#sc88916e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
,載荷上たわみ平均(mm)
,11.919
**鋼板上部欠けモデル [#f148c60a]
箱断面モデルに対して上部60mm部分の鋼板がないモデル
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*半解析 [#lae9b1ad]
半解析で実験と同じく面に載荷してみる。
**補剛材あり三角孔モデル [#be3d440a]
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**箱断面三角孔モデル [#z2b67ee4]
-載荷面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**箱断面(補剛材なし)孔なしモデル [#ge8b74b8]
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-プログラムから作ったモデルと、Salomeから作ったモデルでは...
-プログラムではy方向の負にたわみ、Salomeモデルではy方向の...
（プログラム、Salomeで個人的に作りやすいやり方でモデル化...
*昔のモデル(スパン6m継手なし孔小さめ対傾構あり) [#e434b209]
-sankaku2011
-ヤング率6.929GPa 荷重100kN
,理論値,たわみ(mm)
,初等梁理論,4.960
,ティモシェンコ,6.500
**まとめ(スパン6m継手なし孔小さめ対傾構あり) [#hfd6b642]
-半解析で木材下部のたわみ比較
荷重100kN載荷時のたわみ（半解析での載荷は50kN)
,モデル,たわみ平均(mm)
,対傾構あり三角孔あり,6.42216mm
,対傾構あり三角孔なし,5.9873mm
,対傾構なし三角孔あり,5.64861mm
,対傾構なし三角孔なし,5.501mm
**3点曲げ対傾構なし三角孔なし [#ddc17334]
,たわみを取る部分,たわみ平均(mm)
,木材上部,8.635
,木材下部,5.550
,断面平均,6.223
**半解析対傾構なし三角孔なし [#s0065e6f]
,木材下部,5.501mm
**支間中央に対傾構をつける [#af73c4c4]
-対傾構あり孔なし半解析モデル
半解析なので半分の厚さにする。対傾構をつけると、
上部分の木材から下部分の木材に応力が流れてきて、下部分が...
,木材下部たわみ,5.9873mm
支間中央部分に対傾構が入っているため、断面がふさがっている
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**対傾構あり三角孔あり半解析モデル [#a0fd5c34]
R20
,木材下部たわみ,6.42216mm
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**対傾構なし三角孔あり半解析モデル(比較用) [#f4e8ea2a]
R20
,木材下部たわみ,5.64861mm
-x変位
対傾構があるのと比べ、支間中央部の鋼板のx変位が大きい
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*C3D6とC3D4の混在 [#fa4a71a2]
メッシュの接合性の問題で、C3D8とC3D4はうまく行かないけど、
C3D6とC3D4ならAutoで切っても混在させれそう。
-木材部分C3D6,鋼板部分C3D4
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*三角孔片持ち修正版 [#v28ed3fb]
スパンを短くした。したにある古いものはスパンが端部分だけ...
木材部分上部に線載荷し、たわみは木材部分下部の平均を取っ...
孔ありは全てC3D4要素、孔なしは全てC3D8要素とする。
-面載荷木材部分上部のたわみと実験値
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*三角孔補剛材あり [#e7059fff]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m...
,50,12万,7.06216E-03
,25,52万,6.87534E-03
,23,57万,6.79788E-03
,18,90万,6.75301E-03
,18,101万,,9.94436E-03,6.56369E-03,1.07702E-02
-３点曲げ
,length,要素数,saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,補剛材あり,6.39474E-03,8.87783E-03
*三角孔補剛材なし(ただの箱断面) [#d8c7fc64]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m...
,20,80万,8.41664E-03,1.18110E-02,8.10568E-03,1.26531E-02
-３点曲げ
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,100,56989,8.17722E-03
,40,276870,8.50203E-03
,20,1489814,8.38163E-03
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,補剛材なし,8.24448E-03,1.08202E-02
**三角孔上部鋼板なし（上120mmなし) [#y267aef7]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,14,82万,1.16335E-02,1.80801E-02(要素数84万)
-3点曲げ
,length,要素数,saikasitaたわみ平均(m)
,14,153万6567,1.18096E-02
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,上部鋼板欠け,1.07214E-02,1.64507E-02
**三角孔上部鋼板なし(上60mmなし) [#mafe283a]
*プレストレス木箱桁橋三角孔(縦147.92mm×横143.66mm) [#reeb...
*ティモ孔の鋼板の換算幅算出 [#sfea2d8e]
$(元々の鋼板体積-孔部分体積)/(高さ×長さ)$
-三角孔(大きいタイプ)
,元々の鋼板の幅b, 9.0mm
,ならした鋼板の幅b2, 7.68123mm
-補剛材あり用ティモ孔理論値計算f90&link(keitimohana.f90,h...
bankの修正せん断補正係数では、補剛材部分を考慮していない。
-補剛材なし箱断面用ティモ孔理論値計算f90&link(keianahako....
今までは孔を簡略化して(長方形孔)モデル化、ティモ(孔)の理...
ならした鋼板はもともとよりかなり薄くなるので、
ティモ(孔)と元々のティモシェンコ理論値は開きが出ていたが、
図面の三角孔の面積を使ってならした鋼板の厚さを使うと、
長方形孔をならしたときほどは薄くならない。
そのためティモ(孔)にしても、
元々の値とはあまり変わらなくなる。
-昔の長方形孔
縦180mm $\times$ 横265mmに孔を簡略化
,元々の鋼板の幅b, 9.0mm
,ならした鋼板の幅b2, 2.443mm
,補剛材なしティモ孔たわみ,10.762mm
*補剛材あり [#hda634e4]
**補剛材あり3点曲げ解析 [#r7fb396e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-補剛材あり3点曲げファイルは3thlxxx.inp (xxxはlength)
-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03
**補剛材あり片持ち [#na7dec11]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/R20k...
-補剛材あり片持ちファイルは3lxxxh.*
FEM解析値では、片持ち梁で
,length,要素数,saikasitaたわみ平均値(m)
,14,946976,9.66198E-03
*補剛材なし [#r9d897b0]
**補剛材なし3点曲げ解析 [#rf63141a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-補剛材なし3点曲げファイルは3tlxxx.inp (xxxはlength)
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,100,56989,8.17722E-03
,40,276870,8.50203E-03
,20,1489814,8.38163E-03
**補剛材なし片持ち [#pa8c0825]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-補剛材なし片持ちファイルは3lxxx.* (xxxはlengthの数)
,length,要素数,片持ちsaikasitaたわみ平均(m),荷重(P=1.d-1 ...
,20,717448,1.15251E-02
,14,920757,1.15435E-02
*鋼板欠け [#ibb126e0]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
*三角孔部分 [#ad299fb8]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-完成
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-作成中(14/5/9)。より正確にモデリングする。
--各種寸法から直角三角形の角度θを求めて頂点からの距離xだけ
離れた点に半径20mmの円が触れるようにする。
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-作成中(14/5/8)。
--なんとなく概算で作る。円の半径は20mm
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
*滝田さんのプログラム [#jd30da7e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/taki...
*Mentatメモ [#vde19d27]
***弾塑性解析 [#se5a32fb]
-降伏応力を入れた後、
JOBS→PROPERTIES→JOB RESULTSで
Sterss,
Total Strain,
Elastic Strain,
Plastic Strainにチェック。
***面荷重(Faceload) [#s8b041f7]
-格子状に分割されたある大きな一つの面に対し、面荷重をかけ...
Face Loadを選択し、
ボックスの中に応力を入れることになる。
-この時の応力は、かけたい総荷重/断面積だが、
正確には、
$\frac{かけたい総荷重}{faceの数} \div \frac{断面積}{face...
となる。
-答えは変わらないけど、面はたくさんあるのに応力入れるとこ...
焦って考えて、変なことになることがありそうなので注意。
***inpからABAQUSで読む [#j0b3bcb3]
-inpをインポートする際はinpのヤング率の乗数をEで出力する...
mentatでヤング率を再設定する手間が省ける。
*ccx_2.5調査 [#u510ef96]
**DLOAD検証 [#j3ea62b7]
***DLOAD例 [#p9716c77]
&link(cgxで載荷拘束する方法,http://www.bconverged.com/con...
を見ながら手順通りにやった例。
-inpファイル&link(beam.inp,http://www.str.ce.akita-u.ac.j...
-cgxで作るモデル &link(all.msh,http://www.str.ce.akita-u....
-拘束条件&link(fix1.nam,http://www.str.ce.akita-u.ac.jp/~...
-載荷条件&link(load.dlo,http://www.str.ce.akita-u.ac.jp/~...
-frdファイル&link(beam.frd,http://www.str.ce.akita-u.ac.j...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload.png
***鋼材一本で分割数変えてDLOAD [#g47b86de]
-&link(dlokou4.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--分割数nz=100の時読み込む&link(kou100.dlo,http://www.str...
--分割数nz=500の時読み込む&link(kou500.dlo,http://www.str...
5cm$\times$5cmのスパン50cmの鋼材を片持ちでDLOAD
めんどいので端一面にkousoku,1,3を与えて.dloだけcgxで吐き...
-ある荷重をかけたいなら、pres hoge に、「設定したい荷重/...
,nz分割数,100,500
,FEMたわみ最大値(m),-3.625050e-03,-3.641140e-03
-presを変えていないので、z方向の分割数を変えてもたわみは...
***CLOADとDLOAD比較(鋼材、引張) [#ha1069cc]
1m$\times$1m、長さ50mの鋼材で比較。
-nxとnyは分割なし（1分割)。nzのみ分割(nz=200)
--DLOAD min:-2.425040e-04 at node:801 &link(dlokou4.f90,...
--CLOAD min:-2.425040e-04 at node:801
DLOADにあわせて、固定してない方の端の節点4つに、1つ1/4MN...
-nx=2,ny=1,nz=200 (pres 1.0)
--DLOAD min:-2.425090e-04 at node:1201
--CLOAD min:-2.425090e-04 at node:1201
-nx=1,ny=2,nz=200 (pres 1.0)
--DLOAD min:-2.425090e-04 at node:1201
--CLOAD min:-2.425090e-04 at node:1201
-nx=14,ny=14,nz=500 (pres 1.0)
--DLOAD min:-2.425610e-04 at node:112501
--CLOAD min:-2.425610e-04 at node:112501 &link(ckounx14...
--理論式のたわみ$d=\frac{PL}{EA}$ 2.427184e-04 変位の...
-CLOADでは、分配された節点荷重が重なるところは足し算する...
全く同じ結果になった。このやり方では、節点数が変わってし...
***CLOADとDLOAD比較(木材、片持ち梁等分布) [#aec8bd06]
10cm$\times$10cm 長さ1mの木材を片持ち等分布荷重で解析。
DLOADではpres 0.1とする。
理論値では$k=\frac{5}{6}$として計算を行う。
-nx=1,ny=1,nz=200
--DLOAD min:-2.026590e-02 at node:803
--CLOAD min:-2.026590e-02 at node:803
-nx=14,ny=14,nz=500
--DLOAD min:-2.066020e-02 at node:112711 &link(dmoku500....
--CLOAD min:-2.066020e-02 at node:112711 &link(cmoku14.p...
-ティモ片持ち等分布単位荷重(2kGA) 2.066545E-002
-ティモ片持ち等分布境界値(kGA) 2.183512E-002
幅8cm$\times$高さ4cm 長さ1mの木材 片持ち等分布荷重(pres0....
理論値では$k=\frac{5}{6}$として計算を行う。
--DLOAD min:-9.210080e-02 at node:32476 , &link(png,http...
--ティモ片持ち等分布単位荷重(2kGA) 9.226345E-002 , &link(...
***DLOAD現状 [#qde8aec6]
等分布荷重の時に便利かも。cgx経由じゃないと使えないため、...
cgxで載荷拘束を設定する場合は、
拘束節点を.nam、載荷面情報を.dloに出してそれぞれ.inpでinp...
DLOADにするとできる。
-プログラムで拘束だけ与えて載荷だけcgxの.dloをインプット...
-普通のCLOADで等分布したときとDLOADとで比べてDLOADが
きれいに応力分布出るのか確認。←今のところおんなじになると...
-等分布したい時なら、プログラムに慣れてない人でもcgxで視...
-presの値の設定（presの意味)は間違えないでね。
-後藤ちゃちゃ
DLOADは、&link(図114、115みたいな感じで,http://www.bconve...
**片持ち梁とティモシェンコ理論値比較 [#de1d34be]
簡単な片持ち梁の長方形断面を解いてfrdやdatで見てみたらち...
-ifの練習モデルに載荷した&link(ifyousos.f90,http://www.st...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/if-tawami...
--y軸:たわみ(m)、x軸:nz分割数、赤:ccx、緑:ティモシェンコ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ccxtyousa...
--y軸:相対誤差、x軸:nz分割数
,ティモシェンコ理論値,1.52391E-003(m)
**要素の調査 C3D6とC3D8 [#kc919dcb]
卒論日誌にもっと詳しく書いてたんだけど、間違えて消してし...
直方体でメッシュを生成するプログラムから、その直方体を2分...
-c3d6はこんな感じで作られる
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/cgxc...
-木材でc3d6とc3d8を比較,理論値は初等梁とティモ
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/c3d6...
--y軸:たわみ(m),x軸:nz分割数
-鋼材でc3d6とc3d8を比較,理論値は初等梁
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/kout...
--y軸:たわみ(m),x軸:nz分割数
***雑感 [#nc537c35]
nz分割数が少ない時は(nz<50くらい)、単純にc3d6を作るプログ...
２倍になるので(直方体を2分割してc3d6をつくってるから）、c...
しかし、nzが大きくなり、解が収束してくると(nz>500)、ごく...
*ccx_2.5 ,cgx_2.5覚書 [#ba5c6862]
**RESTART [#d09567e6]
柴田さんの修論日誌に書いてあるとおりにやると、ccx2.5でも....
大体こんな感じで*STEPと*STATICの間あたりに*RESTARTを挿入...
--(
*STEP [#e9412fb8]
*RESTART [#ec504b7c]
*RESTART,WRITE,FREQUENCY=1 [#p14d1e17]
*STATIC,SOLVER=SPOOLES [#g0f640d1]
--)
**座屈解析 [#g753cdc2]
-テクニック
総荷重を1MNかけるときなら、datに出てきた座屈荷重をそのま...
**エラー [#v8b95f9f]
**jacobian問題 : ERROR in e_c3d: nonpositive jacobian [#q...
&link(江村さんの修論日誌,http://www.str.ce.akita-u.ac.jp/...
節点座標が出力フォーマットより細かった時に、その節点を使...
そんな時は出力を1pe17.9とかにして細かい桁まで出してあげる...
***SalomeでAutoでメッシュを切った時に発生 [#nf5aaac5]
Salomeで、プレストレス木箱桁橋三角孔補剛材ありモデルを、
メッシュを全てC3D4で切り、非常に細かく分割した後unv→inp変...
要素数は400万くらい。
要素数70万のものではエラーは発生しなかった。
計算できないくらい細かくなりすぎた？ →細かくなり出力桁が...
--(
*ERROR in e_c3d: nonpositive jacobian
determinant in element 4518956
--)
SalomeでAutoでメッシュを切り、
C3D6とC3D4を混在させたモデルで、C3D4部分を細かくしようと...
全く同じjacobianエラーが発生する。
試しに1つのjacobianエラーが発生したinpを、要素の節点番号...
unv→inp変換プログラムの節点出力を1pe17.9とかにして細かい...
***プログラムで発生 [#s1951dbe]
--(
*ERROR in e_c3d: nonpositive jacobian
determinant in element 20358
*ERROR in e_c3d: nonpositive jacobian
determinant in element 20358
--)
C3D6のところでひとつの要素につき2つのjacobianエラーが出る。
こちらの方は、要素のくるみ方を正しく？(ccx_2.5で対応でき...
ccx_2.5のfortranでの正しい？要素のくるみ方はマニュアルに...
ccx2.6以降では、この問題が解決している可能性があり、
それらのverではccx_2.5でエラーが起きないくるみ方をしてい...
逆にエラーが出る可能性がある。
**f90のコンパイルエラー [#r49f170e]
gfortranでコンパイルを行う時、
--(
relocation truncated to fit: R_X86_64_32S against `.bss'
collect2: ld はステータス 1 で終了しました
--)
みたいなエラーが出てコンパイルできないことがある。
この時は、コンパイルに制限値以上のメモリを使用するために...
なので、制限を取っ払うオプションである-mcmodel=largeをつ...
ダイヤカット円筒をC3D4要素で一気にinp作るプログラムで、要...
--(
gfortran -mcmodel=large nantoka.f90 -o nantoka
--)
**強制終了 [#n2519c86]
ccxで解析中に、ccxが強制終了してしまう。
要素数が多いと発生。
--(
Factoring the system of equations using the symmetric sp...
Using up to 1 cpu(s) for spooles.
強制終了
--)
要するにスペック不足ですってこと。いいPCでやっても出るな...
**ERROR in u_calloc: error allocating memory [#pfe89c4b]
メモリが足りませんよって言われる。
--(
*ERROR in u_calloc: error allocating memory [#k243c6da]
num=566129320,size=4
--)
visionのメモリ2Gのマシンで、140万要素のものを解こうとした...
Corei5 のメモリ4GBのマシンでは、エラーが出ずに普通に解け...
要するにスペック不足ですってこと。
**メモ [#k4a5ceb9]
--.inpをcgxで開いてanimateとかf2とかでたわみ方を確認する...
frdとdatを信用しよう。
--f90で寸法をm、ヤング率と荷重はメガ
-.dloのpresもメガ
**要素 [#p6dbdcd4]
-C3D8 直方体
--C3D20 直方体(C3D8)を一辺に3つ節点ができるように増やした...
-C3D4 四面体
--C3D10 四面体(C3D4)を一辺に3つ以下略
-C3D6 三角柱
--C3D15 三角柱(C3D6)を以下略
-S6,S8 シェル
-他にもいっぱいあるよ。自分で調べて使ってあそぼう
-Rがついてるのは低減積分要素。制御が難しいので、簡単に使...
*Salome覚書 [#f915456c]
**Salome用参考資料 [#ifffd4df]
-チュートリアル的な感じでわかりやすい、超初心者向け
http://tyloroqtism.blogspot.jp/2014/01/salome.html
-salomeの所の解説にある、メッシュの切り方(分割数の指定 な...
http://www.geocities.jp/penguinitis2002/
-いろんなチュートリアル、よく読めば発見があるかも
http://caelinux.org/wiki/index.php/Contrib:Claws/Code_Aster
-マニュアル的なpdf
http://framabook.org/docs/Code_Aster/beginning_with_Code_...
**マウスとキーボードを使った操作方法 [#mf12ea58]
,Ctrl+左クリックしながらドラッグ ,拡大＆縮小
,Ctrl+真ん中クリックしながらドラッグ ,移されている領域を...
,Ctrl+右クリック ,モデル回転
**salomeで.unvをつくりプログラムで.inpにする方法(解決) [#...
***&link(Salomeメモ,http://www.str.ce.akita-u.ac.jp/cgi-b...
材料が2つの時もSalomeメモに書いてます。
**体積を出す [#x603c285]
Geometryで
**unv→inp用の載荷拘束テク [#deb31ce0]
-Salomeで、C3D4などでメッシュを切るとき有効
-きれいに載荷節点、拘束節点をいれたいとき有効
unv→inp変換する際は、Salome上で先に載荷拘束用の区切りをPa...
きれいに載荷拘束用の節点が入る。
unvにするf90では、
if文で
--(
区切りを入れた部分の座標-0.000001d0.and.区切りを入れた部...
--)
みたいに設定していくだけできれいに載荷拘束が入る。
手動で(CLOADで)面載荷したいの時も使える。
Partition入れまくったモデルを、Auto Hex(C3D8)とかでいっぺ...
当然エラーが出るので注意。
そういう意味ではC3D4は便利。(いっぺんに全部切っても大丈夫)
**unv→inp変換時の注意 prism [#q81eb19f]
salomeでメッシュを切った時にprismができるとcalculixでは対...
prismができないようにメッシュを切らなければいけない。
prismは例えばSalomeでAutoで切ったときに、C3D4メッシュとC3...
**メッシュをC3D8要素とかで切っていて、軸方向に多く切りた...
-submeshを使う。
-Geometryで軸方向(多く切りたい方向)の辺すべて(例えば1個の...
-まずBox_1(作りたいモデル)を作る。
-New Entity-Group-Createをクリックし、ラジオボタンで/(辺)...
-Main Shapeには、Box_1を選択。するとBox_1がワイヤーフレー...
-軸方向(z方向)に伸びている辺をクリックし、Add をクリック。
-これを軸方向(z方向)に伸びている辺すべて行う。一つの直方...
4つの辺番号が表示される。
-すべてAddしたら、Apply and Close 。この操作でできるGroup...
-Mesh画面に行き、&link(メッシュをC3D8で切って,http://www....
-その後、create submeshで、MeshにMesh_1を指定、Geometryに...
1D-AlgorithmでWire Distcretisation,HypothesisでNoneを指定...
-すると、新しいNb.Segmentsを作れるので、Number of Segment...
(デフォルトでは、名前が被らないようになっていると思うが、...
-設定し終わったら、Applyする。するとツリーのMesh_1とある...
-Mesh_1を再びcomputeすると、submeshの設定が反映されて再び...
これを行うと、直方体とかなら普通にf90プログラムでnx=,ny=,...
設定するのと同じような感覚でメッシュが切れるし、軸方向分...
unvからinpにしてccxで解いても、寸法・各種数値と分割数が全...
**材料AとBからなるモデルを作りたいが、材料AとBの間をまた...
材料AとBの接合面を、Partitionで区切る。
Partitionでは接合面となるFaceがGeometryにある必要があるが...
ツリーの材料Aか材料Bを右クリックして、CreateGroup→ラジオ...
接合面をクリック→Addで接合面を作るのが効率が良い。
**非一体化モデルを作りたい [#ke26e2f0]
Geometryで、FuseのかわりにCompoundを使ってモデリングする。
Compoundは、Fuseとは違いMeshを切った時に一体化させた境界...
ただし、Compoundした後に境界面をPartitionすると情報が上書...
その境目では節点が1ずつ入る（→つまり一体化する)ようだ。
-compound
選択した物体(複数可)をひとつのグループにまとめてくれる。
一体化してくれるわけではないので注意。
物体Aが物体Bにめり込んだ状態で2つをcompoundしても、
物体Aと物体Bが結合しないし、別々に物体Aと物体Bがあるとい...
**メッシュを細かく切った時エラーみたいなのが出る [#z68748...
elementsが500000(桁があってるか不明)とかより多いよ
みたいなエラーが出たら、Mesh画面でFile→Preferences→Mesh→G...
Mesh informationにあるAutomatic nodes compute limit を増...
ただしエラーが出ていても分割は終了しているっぽい(?)のでそ...
私はlimitを9999999にしている。
**グルーピングが遅い、もっと早くできないの？ [#y84ecb45]
普通は、Create Meshでメッシュを切り、グルーピングを行うが、
Meshを切った後、一旦Saveし、Salomeを再起動させてからグル...
普通にグルーピングするより早く終わる。
なんでだろうね。
**グルーピング後、unv→inp変換はうまく行くのに、inpの要素...
丸みを持つモデルとかで発生する。
とりあえず、SalomeのMeshを右クリック→Mesh Informationでvo...
inpファイルで総要素数を見る。
inpファイルの方が要素数が少なかったら、グルーピングで要素...
丸みがある部分を持つモデルなどをC3D4（角がある要素)で作る...
メッシュを切る際に簡略化されてしまう。（角がある要素では...
そのため、三角孔を持つ木橋などは、メッシュを切った後に、
鋼板部分をグルーピングする際はフィルターに三角孔がある鋼...
すべての要素をグルーピングできなくなる。
こういった場合は孔を開ける前の鋼板をフィルターに設定して...
**C3D4とC3D6を混在させて、inpにして解くとjacobianエラーが...
-ccxエラー - jacobian問題 参照
**Submeshで部位ごとにメッシュを切る [#e1693575]
切り分けたい部位ごとに面などでpartitionを入れておき、
最終的なモデル(たとえばPartiton_1とする)
をExplode-Solid し、Partiton_1から小さいツリーで部位ごと...
Mesh画面でPartition_1を選択しCreate Meshするが、Algorithm...
このMesh_1を選択し、Create sub-meshする。
GeometryでPartition_1の下の細かいツリーの中から切りたい部...
それぞれ部位ごとにSubmeshを作ったら、Mesh_1をComputeする...
Submeshごとに優先順位を変えることもできる。(Mesh_1を右ク...
-Submeshをたくさん作ったり、優先順位の設定が複雑だったり...
たまに
unv出力した際にunv内で要素の順番が変になり、
unv-inp変換プログラムを通した際にinpデータの要素番号が負...
priorityの順番に要素グループを作っていけばこのエラーは回...
**二次要素(中間節点がある要素)をつくる [#f34ce899]
mesh画面でconvert to quadratic
二次要素の方が精度がよいので使いましょう。
二次要素をそのままcalculixに流してとこうとしたところ、
節点の並び方が対応していないのでjacobianエラーが大量に出...
cgxでモデルを確認しても要素がかなり抜け落ちていた。
salomeのautoで切った四面体2次要素をccxで解くのは厳しそう...
**エラー [#b1736fd4]
***Warning:It seems that there is not enough memory~ [#te...
メモリが足りませんよエラー。
--(
Warning
It seems that there is not enoght memory to show the mesh
so that the application may crash. Do you wish to continu...
--)
meshをcomputeしたとき、
Vision 2GのマシンでC3D4で(たぶん、およそ)200万要素くらい...
**材料が2つのとき（を考えてる時に考えたこと) [#v86ece3e]
-例えば部材がbox_1の上にbox_2が乗っかっているモデルだとす...
-box_1とbox_2の境界面をpartitionで区切らなければいけない...
-ツリーのbox_1を右クリックし、Create Groupを選ぶ。
-ラジオボタンで、□(面)を選び、Main Shapeに面を指定する物...
-Geometryにbox_1だけを表示させて、box_2と接している面をク...
-これでpartitionのための設定が終わる。
-box_1とbox_2を一体化させてメッシュを切るので、ここからFu...
-Operations-Boolean-fuse で、object 1 にbox_1 、object 2...
-Fuse1をpartitionで区切る。Operations-Partition を選び、O...
Tool Objectsには、ツリーのbox_1の下できているGroup_1を選...
-Partition_1ができる。
-Mesh画面にいき、Partition_1を選択してメッシュを切る。(こ...
-メッシュを切ったら、グルーピングを行う。Fuseしたモデルを...
-MeshにMesh_1、Element Type にVolume、Nameは(かぶらなけれ...
-Group type ではGroup on filterにチェックを入れる。すると...
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-今回は、Geometry のモデルでフィルターをかけるので、Crite...
選び、Compare にはEqual to を選ぶ。
-そして、Threshold の中の空白をクリックしたら、ツリーの中...
-これで、box_1の部分のメッシュがすべてGroup_1mに入れられ...
-もう一度Create Groupする。今度はbox_2の部分のメッシュを...
-MeshにMesh_1、Element Type にVolume、NameはGroup_2mとか...
-Group type ではGroup on filterにチェックを入れる。すると...
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-Criterion にBelong to Geom を選び、Compare にはEqual to ...
-そして、Threshold の中の空白をクリックしたら、ツリーの中...
-これで、box_2の部分のメッシュがすべてGroup_2mに入れられ...
-Mesh_1を右クリックして、unvファイルを保存する。
メッシュをC3D4で切り、2つにグループ分けしたものなら、ccxc...
partitionを使わない方法だと、C3D4とかだとbox_1とbox_2の間...
となると、メッシュの切り方を変えてBox_1とBox_2で別々に切...
Compound Meshを作ってその後グループ分け、という作業が必要...
更に、Compoundするならその接合面で節点がくっつくほうがい...
Box_1の上にBox_2が乗っかってるモデルとするなら、FaceのBox...
という手順が必要かもしれない。調査中。
partition使ってないからこんなことになるのか？partition使...
-Compoundとか使うとunvの要素の入り方が変わってくる？。普...
最初が-864とかから始まる。
-partitionを使う方法
予め面を作っておき、partitionを使って区切ってからメッシュ...
下の３角孔木箱桁橋でやったら要素の抜けが起こった。鋼板と...
partitionで区切ってるのになあ…。(鋼板の方の面をグループ分...
→解決しました。フィルターの範囲を大きくしなければいけませ...
*パソコン関係覚書 [#e11ccfff]
**viコマンド [#x4db6525]
viのコマンド、特に移動系のコマンドなどを使いこなせてない...
-移動系(0と$は特に便利)
,0 ,その行のはじめに移動
,$ ,その行の最後に移動
,Ctrl e ,一行下へスクロール
,Ctrl y ,一行上へスクロール
,Ctrl f ,1ページ分進む
,Ctrl b ,1ページ分戻る
,H ,表示されている画面の一番上の行に移動
,M ,表示されている画面の中央の行に移動
,L ,表示されている画面の最後の行に移動
-入力系
,s ,カーソルの位置の文字を削除して入力モードに
,Ctrl s ,キー入力無効状態になる
,Ctrl q ,キー入力無効を解除
**aliasでコマンド短縮 [#l466b40b]
プログラムをコンパイルする際、gfortranといちいち打つのが...
~/.bashrcの中で# some more ls aliasesの後辺りに
--(
alias gf='gfortran -o'
--)
と追加しておくとgf がgfortran -o の短縮コマンドになる。た...
**前に使ったコマンドを見る [#u4a4ae2c]
--(
history|grep nantoka
--)
でnantokaというコマンドを使った履歴を検索できる。
番号とその時使ったコマンドがずらっとでて、いちいち打ち直...
**inpファイルを見やすくする [#s9042b53]
&link(vimrc_ubuntu,http://www.str.ce.akita-u.ac.jp/~gotou...
**メモリとかCPUの情報を見る [#zdf21144]
--(
more /proc/cpuinfo
more /proc/meminfo
--)
**外付けHDにバックアップをとる [#td205260]
--(
rsync -auv コピー元 コピー先
--)
**jobの最後にrm jobを入れとく [#racb5007]
jobの使い方は後藤資料の中でぐぐればあります。で、jobの最...
--(
rm job
--)
を入れとくと、!./を使ったときjobが再実行されて、慌ててCtr...
解いた結果のdatファイルが全部真っ白になっているということ...
**libre office [#mbf46ca9]
2014/4/21現在、Ubuntu12.04環境では、
Libre Officeが入っているわけだけど、どういうわけかすごく...
スライド作成などしたい場合には、事前に
verを3.なんたらから4.2.3にバージョンアップしたほうがいい...
コマンドはlibre office ubuntu12.04 とかでぐぐると出てく...
**役立つページ [#o35865fd]
-http://www.str.ce.akita-u.ac.jp/~chida/memo.html#i0
-http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%...
*gnuplotメモ [#wcc133bf]
**gnuplotで値をグラフにして見たいって人は、 [#k088c10c]
-まず、txtデータを作る。
hoge.txtとか(何でもいい) にx軸データ y軸データ を入れる。
hoge.txtの中身が、
--(
1 50 ←1列目がx軸データ、空白を開けたあと、2列目をy軸...
100 5000
--)
なら、
gnuplot で、
--(
pl 'hoge.txt' w l
--)
とすると、グラフ画像が出てきて、見れる。
**グラフをpngにして保存したいって人は、 [#u7fb64b4]
-まずpngを作る
--(
set term png
set output 'nantoka.png' ←nantokaはグラフにつけたいファイ...
--)
-そんで、グラフ出力する
--(
pl 'hoge.txt' w l
同時にいくつもグラフを載せたい場合は、
pl 'hoge.txt' w l , 'moge' w l , 'mage' w l
--)
set outputしたあとに、plotすると、グラフが表示されずに、....
この2つの操作でグラフを.pngにして保存できる。
**1つのテキストファイルからいくつもグラフを書きたい！って...
hoge.txtの中身を、
--(
1 50 100 0.2
100 5000 10000 20
--)
というふうにデータを並べてって、gnuplotで、
--(
pl 'hoge.txt' using 1:2 w l ,'hoge.txt' using 1:3 w l ,'h...
--)
みたいにしてusingを使う。
using1:2 で、txtデータの1列目がxに入り、2列目がyに入るっ...
**線の太さをもっと太くしたい！って人は、 [#y04eee9e]
gnuplotでそのままグラフを作成すると、
スライド作成とかの時、グラフの線が細くて見にくいと言われ...
グラフの線を太くするには、linewidth を使う。 省略形でlwも...
lw 2みたいに、その後に半角スペースを空けて、整数を入れる...
--(
pl 'hoge.txt' w l lw 2
--)
みたいにすると、太くなる。(だいたい2が見やすいと思う）
**グラフの線のタイトルを消してクールにしたい！って人は、 ...
notitleを使う。
発表の時は、スライド作成ソフトやxfigとかで、グラフに矢印...
そういう時には右上に自動生成されるグラフの線のタイトルが...
--(
pl 'hoge.txt' w l notitle
--)
*構造実験検証 [#x67cf198]
一人で（龍さんに手伝ってもらったりして）こっそりやったデ...
-4/24午後7時 二回とも成功(変な測定値は出なかった)
--2回目データ:&link(140424a_2s.csv,http://www.str.ce.akit...
--1回目データ:&link(140424a_1s.csv,http://www.str.ce.akit...
*3D [#d2eb7ba5]
**引張試験 試験体 [#h211cd60]
引張試験のための供試体を作成する。
ものづくりセンターの3Dプリンタでは造形情報をstlファイルと...
SalomeのGeometryでモデルを選択→export→STLAⅡfilesで吐き出...
**引張試験 [#f5580cac]
***メモ [#a44cbbb9]
-ロガーにつながってるPcのPCD-30Aの集録データ数 の欄で何秒...
-doループの終わりをデータ終わりの秒数にするだけでヘッダー...
*回転らせん [#eeb2f178]
遊びで適当に節点だけ打ってみたけど、
中間節点を作るのが難しそう
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*プログラム類↓ [#wa71b442]
&link(卒論日誌,http://www.str.ce.akita-u.ac.jp/cgi-bin/gw...
*プレストレス木箱桁橋の角材の間にゴムを挟めるモデル [#y89...
接触解析の代わりになるかも という事で滝田さんのプログラム...
ゴムの寸法や、分割数、ヤング率(ELgom)とかは各自で変えてく...
節点座標の出力桁を多くして、
ゴムを小さくした時に節点座標が出力桁以下になり、
四捨五入されてしまいメッシュの大きさが均一にならないこと...
&link(gom2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
*B32用のunv→inp変換 [#v1d92ed0]
2次元のトラスにしか対応してません。
3次元には非対応。
主な手順
-まずSalomeでyz平面上にトラスを作る
-点を打ってってそれをLineにして
-トラスができたら、全部のLineを選択して1つのWireをつくる
-その後,作ったWireを,
-1Dでメッシュを切るが、number of segmentsで分割数を３以上...
-Meshを.unvに吐き出して&link(b32unv.f90,http://www.str.ce...
*2014年度後期 [#r0745b2b]
,曜日,12,34,56,78
,時間帯,8:50-10:20,10:30-12:00,12:50-14:20,14:30-16:00
,月,,,基礎物理学実験D TA ,TA
,火,,,,
,水,,構造設計学特論,,
,木,,地盤工学特論,工学資源学特論,
,金,,,水防災工学特論,
集中講義 工学資源特論
*2014年度前期 [#z45bf2e6]
,曜日,12,34,56,78
,時間帯,8:50-10:20,10:30-12:00,12:50-14:20,14:30-16:00
,月,,構造力学特論,,パソコンゼミ
,火,,都市システム計画特論,土質工学特論,
,水,,材料設計特論,,
,木,,,交通システム計画特論(英語ゼミ),構造実験
,金,,,,
終了行:
#contents
*図などはsyuronフォルダの中にあります [#eed59bd8]
修論本体に使った座屈モードなどの
図は2014i5/tex/syuron/ にだいたい置いている。
*tugite130710鋼材荷重-ひずみ関係 [#p6cccb26]
tugite130710.xlsxのデータをゼロ点修正
座屈しそうな中央のトラス部分にひずみゲージをつけて破壊試...
ch33～36の斜材・垂直材ともに引張りのまま破壊している。
支間中央部は継手で補剛されているので心配なさそう。
縦軸 荷重(kN)、横軸 ひずみ(μ)
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
**予備試験 [#g8532b48]
tugite130710yのデータをゼロ点修正
縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
**両方プロット [#rba2340d]
tugite130710.xlsx、tugite130710yのデータをゼロ点修正
縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*目標 [#o5a46862]
修士発表時間
http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%c...
-等圧縮を受ける板の座屈を確認
*板の座屈 [#h2a0e021]
**等圧縮 [#m6184594]
***正方形板 4辺単純支持 [#j5c4a131]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
&link(kouitass8c.f90,http://www.str.ce.akita-u.ac.jp/~got...
$\sigma_{cr}=k\frac{D\pi^2}{hl^2}$
$N_{cr}=k\frac{D\pi^2}{l^2}$ よって
$N_{cr}=h \sigma_{cr}$
ここで、$\sigma_{cr}$は座屈応力、$N_{cr}$は座屈荷重、
Dは板の曲げ剛性、hは板厚、lは基準にとった板の辺長、
正方形板であり、4辺単純支持なので、k=4
$N_{cr}=0.009 \times 241.2955$
$ =2.17166 [MN] $
,座屈モード,ccx座屈荷重(MN) ,理論値との誤差(%)
, 1 ,0.2185645E+01,0.643977
, 2 ,0.2512923E+01
, 3 ,0.6733093E+01
, 4 ,0.7050507E+01
, 5 ,0.7651342E+01
板の厚さ方向より幅方向高さ方向分割の方が解に与える影響が...
幅方向高さ方向をもっと細かくすると理論値より若干小さい値...
-1次座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
***長方形板 4辺単純支持 [#s9e4b45f]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
&link(kouitassk9.f90,http://www.str.ce.akita-u.ac.jp/~got...
長方形板だけどk=4になる様なa,bの比なので、正方形板と理論...
$N_{cr}=0.009 \times 241.2955$
$ =2.17166 [MN] $
,座屈モード,ccx座屈荷重(MN) ,理論値との誤差(%)
, 1 , 0.2147385E+01 ,1.13735
, 2 , 0.2163955E+01 ,
, 3 , 0.3919156E+01,
, 4 , 0.5249193E+01,
, 5 , 0.5451060E+01,
-1次座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**純曲げ [#pa5aaff7]
三角形分布でも総荷重を分配で良さそう。
一つの三角形に1/4MNを分配してやる。
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**参考になりそうな文献 [#o47661f2]
-長方形板の分岐階層構造の境界条件による退化
https://www.jstage.jst.go.jp/article/jscej1984/1995/507/1...
-弾性矩形板の曲げ解析手法の一般化について
http://ci.nii.ac.jp/els/110007887197.pdf?id=ART0009558779...
-等曲げを受ける弾性矩形板に生じる二次座屈現象の解明
https://www.jstage.jst.go.jp/article/jscej1984/1995/519/1...
*Salomeで作ったモデルで座屈 [#e2f56b5a]
**応力集中部を細かく切ったモデル [#od3e26fe]
非一体化にして座屈させるところの鋼板を多く切る。
応力集中部を細かくしていくと1次モードは飛び出しがないモー...
鋼板要素数1は非一体化にしてるところと一体化しているところ...
鋼板要素数2は残りの鋼板全部。
-孔なし補剛材なし
,総節点数,木材要素数,非一体化鋼板要素数,鋼板要素数1,鋼板...
,631768,773611(l12),1389462(l1.6),347918(l2),291211(l8),2...
**補剛材なし [#e1ae5037]
-孔なし補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,389725,704146(l12),650860(l2),298061(l8),1次モードは正常...
,535235,601922(l12),1389462(l1.6),331039(l8),1次モードは...
--プログラムで作ったモデルと近い値が出る。
--1次モードでは座屈荷重77.7kN、2次モードでは座屈荷重78.5k...
2次モードでは要素が飛び出す。
2次モードでは鋼板の幅方向の変位が1次モードより小さい。
-メモ
--l12-l1-l8だとエラーで切れない
--3孔鋼材非一体化部だけをl1で切ると要素数470万くらいにな...
--l10-l1-l8でもエラーで切れない、要素数1000万くらいになり...
--l8-l1-l8もエラーで切れない
-孔あり補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,342007,599145(l12),577961(l2),231480(l8),1次モードは正常...
,522828,1017878(l12),1109490(l1.6),234967(l8),1次モードは...
--l12-l1.6-l8は切れる
--孔をつけたモデルでも解析がうまく行く。孔ありの方が1次の...
**補剛材あり [#yc7315fa]
-孔なし補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,882260,1798392(l12),1653936(l2),596385(l8),要素が飛び出...
-孔あり補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要...
,523050,598470(l12),1251573(l2),452306(l8),要素が飛び出す...
--孔なしのほうが木材も鋼板も細かく切れているので1次座屈荷...
**今までのをやり直す。 [#c05c0b12]
salomeで木材部を削って一部非一体化にしたモデルを分割して...
-孔あり補剛材なしsalomeでモデリング
要素はC3D4(四面体)のAutoでlengthを変えて切る
,length,総節点数,木材要素数,鋼板要素数,結果
,12,17万,58万,15万,座屈させると要素が飛び出す。1次$P_{cr}...
,10,,,,なぜかエラーがでて切れない
,8,50万,203万,36万,要素の飛び出しがなくなる。1次$P_{cr}$=...
,7,581146,2396666,386169,要素の飛び出しが見られない。1次$...
,6,92万,340万,104万,要素数が多くなりすぎてccxで解けない(...
鋼板をもっと多く切ってみる必要がありそう
2次要素とかで切ってみるとどうか？
-孔なし補剛材なしsalomeでモデリング
,length,総節点数,木材要素数,鋼板要素数,結果
,12,181632,765824,178789,座屈させると要素の飛び出しが起こ...
,10,,,,なぜかエラーがでて切れない
,9,,,,エラーが出て切れない
,8,521420,2062159,376681,座屈させると要素の飛び出しが起こ...
,7,597914,2412348,450689,座屈させると要素の飛び出しが起こ...
*プログラムで作った孔なし非一体化モデル [#ud1a5503]
(鋼材と木材を支間中央部から最初のPC鋼棒のアンカープレート...
onsih4z.f90,onsi4z2.f90
を座屈させる。
1MNかけているので出力された座屈荷重がそのまま読める。
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,単純支持,3.701E-002
,固定支持,0.1480
**孔なし補剛材なし非一体化モデル [#x7a11643]
,ccx解析結果,
,n次モード,ccx座屈荷重(MN)
,1,0.7846011E-01
,2,0.2280486E+00
,3,0.2489195E+00
,4,0.4002716E+00
,5,0.4099577E+00
--オイラーの固定支持よりかなり低い値が出る。
--salomeで作ったモデルのような飛び出し部は確認されない。
--salomeで作った孔なしモデルは92kNくらいで座屈するが、そ...
-1次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-1次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-2次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-2次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**孔なし補剛材なし上の段だけ非一体化モデル [#x18a4a36]
下の段の木材と鋼板は一体化しているモデル
,n次モード,ccx座屈荷重(MN)
,1,0.8540278E-01
,2,0.2201420E+00
,3,0.3317080E+00
,4,0.4245944E+00
,5,0.4385326E+00
-1次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**孔なし補剛材あり非一体化モデル [#ca81a619]
--最初に木材部分の応力集中部(載荷点)が破壊し、
全てのモードで木材部の座屈挙動になるため、
鋼板部の座屈挙動が見ることができない。
--ためしに補剛材部分の寸法をものすごく薄くして補剛材なし...
断面にして解析したら上の補剛材なしプログラムで解いた結果...
,ccx解析結果,全て木材部の破壊
,n次モード,ccx座屈荷重(MN)
,1, 0.1028976E+00
,2, 0.1759758E+00
,3, 0.1906934E+00
,4, 0.2086057E+00
,5, 0.2153555E+00
**孔なし補剛材あり上の段だけ非一体化モデル [#g2e56c6b]
下の段の木材と鋼板は一体化しているモデル
--最初に木材部分の応力集中部(載荷点)が破壊し、
全てのモードで木材部の座屈挙動になる
,n次モード,ccx座屈荷重(MN)
,1,0.1049848E+00
,2,0.1809716E+00
,3,0.1925376E+00
,4,0.2122353E+00
,5,0.2144374E+00
*Asterのbuckling [#sa2492a9]
-解説pdf
http://framabook.org/docs/Code_Aster/beginning_with_Code_...
の15.2
ver10.8のcommファイルを使う。
これ見ながらcommファイルはオールグリーンにできたけど、pos...
たぶんやり方が違うのだろう。
&link(commファイル,http://www.str.ce.akita-u.ac.jp/~gotou...
-Asterのマニュアル(公式HP)
http://www.code-aster.org/V2/doc/default/en/index.php?man...
左上で言語をenglishにできる
-Bucklingについて
--http://www.code-aster.org/V2/doc/v11/en/man_u/u2/u2.08....
-MODE_ITER_SIMULTのマニュアル
--http://www.code-aster.org/V2/doc/default/en/man_u/u4/u4...
MODE_ITER_SIMULTのオプションを使いそう
--lateral buckling of a beam (spill)
http://www.code-aster.org/V2/doc/default/en/man_v/v2/v2.0...
どうやって解くかは書いてない
Aster→2材料で解くのと異方性はできる。座屈は調査中
-ccx2.7のマニュアル-*Friction
http://web.mit.edu/calculix_v2.7/CalculiX/ccx_2.7/doc/ccx...
-ccx2.5のマニュアル*Friction
http://its-cs1.its.uni-kassel.de/CalculiX-2.5/html/ccx/no...
*弾性床上の座屈 [#uf9a40dc]
http://freecaetester.blog62.fc2.com/category19-0.html
上の記事を参考にして一部非一体化にしたプログラムに接触解...
-接触入れるとWARNINGが出る。
--(
*WARNING in statics: a nonlinear geometric analysis is r...
but no time increment nor step is specified
the defaults (1,1) are used
--)
-接触解析で要素分割を細かくするとエラーになる
--(
*ERROR: increment size smaller than minimum
best solution and residuals are in the frd file
--)
マスターのヤング率を低くする(ezz/10とか)と何故か解ける
*座屈解析 [#p7cd23f7]
**鋼板のオイラー座屈 [#c1b15f5d]
プレストレス木箱桁橋 鋼板部プレストレス固定支持間315mmの...
--プレストレス鋼棒間730mmの半分365mmからアンカープレート1...
幅120mm高さ9mmスパン315mmで単純支持にしてオイラー座屈させ...
細長比121.24
オイラー座屈荷重理論値$Pcr=0.14937MN$
,ccx解析結果,,
,n次モード,ccx座屈荷重(MN),オイラー荷重理論値との相対誤差...
,1,0.1548903E+00,3.670
,2,0.6218848E+00,
,3,0.1393993E+01,
,4,0.1936985E+01,
,5,0.2383931E+01,
長方形断面だからちょっと誤差が出る。
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-スパンを630mmにしたやつ
,n次モード,ccx座屈荷重(MN),オイラー荷重理論値との相対誤差...
,1,0.1472202E+00,1.439
,2,0.3013863E+00
,3,0.5895368E+00
,4,0.8915795E+00
,5,0.1322299E+01
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**木材と鋼板1枚 上一段だけモデル [#tcbb0247]
木材と鋼板1枚で、支間中央部のPC鋼棒間だけが木材と鋼板が非...
,座屈モード,座屈荷重(MN)
,1, 0.1349740E+00
,2, 0.4743007E+00
,3, 0.5364287E+00
,4, 0.9919370E+00
,5, 0.1180486E+01
1次モード座屈挙動
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
接触解析にしてないのでオイラー座屈、ccxで鋼板1枚を座屈さ...
**オイラー座屈応力 [#a3068bda]
-&link(zakutu,http://www.str.ce.akita-u.ac.jp/~gotou/pas/...
**1/4解析モデル [#y4702bc7]
まず普通の1/4解析モデルを作る。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
どちらも下の 半解析-箱断面三角孔モデル のx変位、y変位と同...
たわみ平均なども大体同じ値になる。
*鋼板の座屈 - 一体化だけど一部だけ鋼板と木材の間を削った...
とりえあず、中央部の木材と鋼板がくっついている部分の、木...
実物は上弦材も下側木材もPC鋼棒間で非一体化なので、FEMモデ...
**三角孔あり補剛材なし [#df7f9762]
-面載荷しているモデル。length12。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-木材部分に線載荷したモデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
***座屈荷重(1MN載荷時)とアニメーション [#odb2bae9]
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,単純支持,3.701E-002
,固定支持,0.1480
-ccxのBuckle
,モード,座屈荷重(MN),,
,,線載荷,面載荷(載荷試験装置寸法で載荷)
,1, 0.9249756E-01,0.9726263E-01
,2, 0.1472920E+00,0.1548828E+00
,3, 0.1699003E+00,0.1734040E+00
,4, 0.3865533E+00,0.3048128E+00
,5, 0.4932473E+00,0.3726462E+00
今回は固定支持であると考えられるため、ccxのモード2の計算...
-モード１
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
何かが飛び出す現象がみられる。
アニメーションで詳しく見てみる。
--モード1アニメーション、Amplitude+100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
--モード1アニメーション、Amplitude-100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
この2つのモード1アニメーションを見る限り木材部ではなく鋼...
-モード２
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**三角孔あり補剛材あり [#r1017a02]
***座屈荷重(1MN載荷時)とアニメーション [#o50a3b1a]
-ccxのBuckle(線載荷)
,モード,座屈荷重(MN)
,1, 0.9889383E-01
,2, 0.1574743E+00
,3, 0.3737713E+00
,4, 0.4575766E+00
,5, 0.5303813E+00
--モード1アニメーション、Amplitude+100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
・補剛材なしとありでは飛び出す方向が違う
→単に方向が違うだけなので補剛材のあるなしでモードが変わる...
・またはメッシュ分割の違い
（鋼板部メッシュ数は補剛あり18万、補剛なし15万でそんな変...
によって違うモードになったかも
**孔なし補剛材なし [#y0e7fbf9]
,モード,座屈荷重(MN)
, 1, 0.9200184E-01
, 2, 0.1465025E+00
, 3, 0.2212824E+00
, 4, 0.3723138E+00
, 5, 0.4482462E+00
4hanhakoananasi2
1次、2次モードの座屈荷重は鋼板に孔が無くても変わらない。
**孔なし補剛材あり [#rac8af2b]
,モード,座屈荷重(MN)
, 1, 0.1209061E+00
, 2, 0.1913690E+00
, 3, 0.3107388E+00
, 4, 0.4833085E+00
, 5, 0.5042332E+00
4hanhakoananasihogou
孔なしの場合、補剛材があると座屈に対してかなり補剛効果が...
**まとめ [#o07e4fc6]
オイラーの座屈荷重による計算結果と解析による補剛材なしの...
ccxによる解析結果のモード1は、断面に対し斜め方向からなど...
普通の弱軸よりもさらに弱い軸があるのではないかと考えられ...
どのモードを見ても解析結果では補剛材なしより補剛材ありの...
ちなみに鋼板と木材部を一体化させると先に上部木材部が破壊...
一体化させないモデルよりも高い座屈荷重が出る。
-オイラーの座屈荷重
,条件,座屈荷重(MN)
,固定支持,0.1480
-ccxによる座屈解析結果
,モード,座屈荷重(MN),
,,三角孔あり補剛材なし,三角孔なし補剛材なし,三角孔あり補...
,1,0.9249756E-01,0.9200184E-01,0.9889383E-01,0.1209061E+00
,2,0.1472920E+00,0.1465025E+00,0.1574743E+00,0.1913690E+00
,3,0.1699003E+00,0.2212824E+00,0.3737713E+00,0.3107388E+00
,4,0.3865533E+00,0.3723138E+00,0.4575766E+00,0.4833085E+00
,5,0.4932473E+00,0.4482462E+00,0.5303813E+00,0.5042332E+00
*Aster [#me5f5ff7]
最終的な解析モデルをExplodeして
-荷重をかける面
-固定する面
を作っておく
**検証 [#t6e667dc]
簡単な片持ち梁を軸方向に押してみる。
-断面2cm×2cm、長さ40cm
-鋼材206GPA,ポアソン比0.3
asterはデフォルトの載荷方法だと面を指定してその面に対して...
その他の載荷方法は後でviでソース(.comm)をいじるか、
Aster画面でEficasを起動させてEficas経由で.commをいじるか...
,理論式,理論値(m)
,$d=\frac{PL}{EA}$,-4.85437E-07
,,変位最大値dmax(m),理論値との誤差(%)
,ccx節点載荷,-4.979490e-07,2.577%
,ccx面載荷,-4.846880e-07 ,0.1543%
,aster面載荷,-4.84655e-07,0.1611%
大体同じになる
**Asterで異方性 [#y23a9d1b]
ccx・理論値と比較し、異方性がちゃんと解けるか確認
***等方性で解く [#z4f5c2c1]
鋼材のたわみ(m)
節点載荷
,ccx,9.6939E-06
,aster,9.72292e-6(大体)
,初等梁,9.7087E-006
,ティモ,1.0146E-005
***直交異方性で解く [#z4f199fd]
-正方形断面の一辺20cm 長さ2mの木材
z軸方向ヤング率ezz=7.694GPa
exx=eyy=ezz/25d0
せん断弾性係数 一律でezz/15d0
片持ち梁とし、自由端に荷重100Nとする。
salomeではccxと同じように節点荷重として載荷する。
nx=ny=10,nz=100
,,たわみ最大値(m)
,初等梁理論値,2.5994E-004
,ティモ理論値,2.7164E-004
,ccx解析値,2.71617E-04
,aster,2.71617e-4
zがスパン(長さ)方向でモデリングしたなら、N(normal)に軸方...
Lに長手方向、Tに残りの1つでよいと思う。
当然せん断弾性係数やポアソン比も適切にいれる必要があると...
今回はccxのプログラムf90と同じになるように打ち込んだ。
-使ったファイル
--&link(hdfファイル,http://www.str.ce.akita-u.ac.jp/~goto...
--&link(commファイル,http://www.str.ce.akita-u.ac.jp/~got...
寸法をm単位で入力してるので注意。
このやり方で木材を短くしたり長くしたりして解析したが、
ほぼすべてsalomeとccxの解が同じになった。
-&link(マニュアル,http://www.code-aster.org/V2/doc/v11/fr...
のP13,14にLTN系の座標と例が書いてある。
-Asterでの異方性の設定について書いてあるもの
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04....
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04....
--http://www.code-aster.org/forum2/viewtopic.php?id=4026
--http://www.code-aster.org/V2/doc/v11/fr/man_u/u4/u4.43....
のP13,14
**Commandファイル [#z0d028a9]
***荷重 [#oab18278]
-PRES_REP
デフォルトの設定。面に対して垂直な荷重をかける。
-FORCE_FACE
グループと方向を指定して荷重をかける。
*メモ [#x1bcae7a]
最近のモデル(スパン7m継手あり孔大きめ、補剛材ありなし)
-tugite1302～,1307～
*継手ありモデル [#tade4196]
-継手あり補剛材ありモデルの様子 継手部(一体化,半解析)
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
継手があると支間中央部の鋼板の三角孔が継手鋼板で塞がれる...
*木材利用研究発表会 [#qc0a738b]
2014年 &link(発表スライド,http://www.str.ce.akita-u.ac.jp...
*継手なしモデルまとめ [#j22fb978]
**3点曲げと半解析による比較 [#rd37bb00]
**木材に線載荷 [#v22e085e]
-ヤング率7.694GPa 荷重100kN
,モデル,解析手法,載荷上木材部たわみ平均値(mm),三角孔あり...
,補剛材あり三角孔ありlength25,3点曲げ,9.6324,,
,,半解析,10.2294
,補剛材あり三角孔なし,3点曲げ,8.99085,6.660%
,箱断面三角孔あり,3点曲げ,11.919
,,1/4解析,11.549
,箱断面長方形孔,3点曲げ,15.980
,箱断面三角孔なし,3点曲げ,10.9863,7.825%
,上部60mm鋼板欠け三角孔あり,3点曲げ,13.5324
**半解析と1/4解析で実験と同じように面に載荷 [#a9c7ca0a]
,モデル,解析手法,載荷上木材部たわみ平均(mm)
,補剛材あり三角孔あり,半解析,10.5003
,箱断面三角孔あり,半解析,12.4034
,,1/4解析,12.529
,上部60mm鋼板欠け三角孔あり,半解析,13.9751
-グラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*3点曲げ木材上部分載荷 [#s252a518]
-全てC3D4要素
-支間中央木材部分に線載荷
-載荷上たわみ平均は、線載荷部からz軸方向に5mm離れたところ...
**補剛材ありモデル [#od9f0ef7]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
,載荷上たわみ平均(mm)
,9.98675
**箱断面モデル [#sc88916e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
,載荷上たわみ平均(mm)
,11.919
**鋼板上部欠けモデル [#f148c60a]
箱断面モデルに対して上部60mm部分の鋼板がないモデル
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*半解析 [#lae9b1ad]
半解析で実験と同じく面に載荷してみる。
**補剛材あり三角孔モデル [#be3d440a]
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**箱断面三角孔モデル [#z2b67ee4]
-載荷面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**箱断面(補剛材なし)孔なしモデル [#ge8b74b8]
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-プログラムから作ったモデルと、Salomeから作ったモデルでは...
-プログラムではy方向の負にたわみ、Salomeモデルではy方向の...
（プログラム、Salomeで個人的に作りやすいやり方でモデル化...
*昔のモデル(スパン6m継手なし孔小さめ対傾構あり) [#e434b209]
-sankaku2011
-ヤング率6.929GPa 荷重100kN
,理論値,たわみ(mm)
,初等梁理論,4.960
,ティモシェンコ,6.500
**まとめ(スパン6m継手なし孔小さめ対傾構あり) [#hfd6b642]
-半解析で木材下部のたわみ比較
荷重100kN載荷時のたわみ（半解析での載荷は50kN)
,モデル,たわみ平均(mm)
,対傾構あり三角孔あり,6.42216mm
,対傾構あり三角孔なし,5.9873mm
,対傾構なし三角孔あり,5.64861mm
,対傾構なし三角孔なし,5.501mm
**3点曲げ対傾構なし三角孔なし [#ddc17334]
,たわみを取る部分,たわみ平均(mm)
,木材上部,8.635
,木材下部,5.550
,断面平均,6.223
**半解析対傾構なし三角孔なし [#s0065e6f]
,木材下部,5.501mm
**支間中央に対傾構をつける [#af73c4c4]
-対傾構あり孔なし半解析モデル
半解析なので半分の厚さにする。対傾構をつけると、
上部分の木材から下部分の木材に応力が流れてきて、下部分が...
,木材下部たわみ,5.9873mm
支間中央部分に対傾構が入っているため、断面がふさがっている
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**対傾構あり三角孔あり半解析モデル [#a0fd5c34]
R20
,木材下部たわみ,6.42216mm
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
**対傾構なし三角孔あり半解析モデル(比較用) [#f4e8ea2a]
R20
,木材下部たわみ,5.64861mm
-x変位
対傾構があるのと比べ、支間中央部の鋼板のx変位が大きい
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*C3D6とC3D4の混在 [#fa4a71a2]
メッシュの接合性の問題で、C3D8とC3D4はうまく行かないけど、
C3D6とC3D4ならAutoで切っても混在させれそう。
-木材部分C3D6,鋼板部分C3D4
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*三角孔片持ち修正版 [#v28ed3fb]
スパンを短くした。したにある古いものはスパンが端部分だけ...
木材部分上部に線載荷し、たわみは木材部分下部の平均を取っ...
孔ありは全てC3D4要素、孔なしは全てC3D8要素とする。
-面載荷木材部分上部のたわみと実験値
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/grap...
*三角孔補剛材あり [#e7059fff]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m...
,50,12万,7.06216E-03
,25,52万,6.87534E-03
,23,57万,6.79788E-03
,18,90万,6.75301E-03
,18,101万,,9.94436E-03,6.56369E-03,1.07702E-02
-３点曲げ
,length,要素数,saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,補剛材あり,6.39474E-03,8.87783E-03
*三角孔補剛材なし(ただの箱断面) [#d8c7fc64]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m...
,20,80万,8.41664E-03,1.18110E-02,8.10568E-03,1.26531E-02
-３点曲げ
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,100,56989,8.17722E-03
,40,276870,8.50203E-03
,20,1489814,8.38163E-03
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,補剛材なし,8.24448E-03,1.08202E-02
**三角孔上部鋼板なし（上120mmなし) [#y267aef7]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,14,82万,1.16335E-02,1.80801E-02(要素数84万)
-3点曲げ
,length,要素数,saikasitaたわみ平均(m)
,14,153万6567,1.18096E-02
-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,上部鋼板欠け,1.07214E-02,1.64507E-02
**三角孔上部鋼板なし(上60mmなし) [#mafe283a]
*プレストレス木箱桁橋三角孔(縦147.92mm×横143.66mm) [#reeb...
*ティモ孔の鋼板の換算幅算出 [#sfea2d8e]
$(元々の鋼板体積-孔部分体積)/(高さ×長さ)$
-三角孔(大きいタイプ)
,元々の鋼板の幅b, 9.0mm
,ならした鋼板の幅b2, 7.68123mm
-補剛材あり用ティモ孔理論値計算f90&link(keitimohana.f90,h...
bankの修正せん断補正係数では、補剛材部分を考慮していない。
-補剛材なし箱断面用ティモ孔理論値計算f90&link(keianahako....
今までは孔を簡略化して(長方形孔)モデル化、ティモ(孔)の理...
ならした鋼板はもともとよりかなり薄くなるので、
ティモ(孔)と元々のティモシェンコ理論値は開きが出ていたが、
図面の三角孔の面積を使ってならした鋼板の厚さを使うと、
長方形孔をならしたときほどは薄くならない。
そのためティモ(孔)にしても、
元々の値とはあまり変わらなくなる。
-昔の長方形孔
縦180mm $\times$ 横265mmに孔を簡略化
,元々の鋼板の幅b, 9.0mm
,ならした鋼板の幅b2, 2.443mm
,補剛材なしティモ孔たわみ,10.762mm
*補剛材あり [#hda634e4]
**補剛材あり3点曲げ解析 [#r7fb396e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-補剛材あり3点曲げファイルは3thlxxx.inp (xxxはlength)
-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03
**補剛材あり片持ち [#na7dec11]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/R20k...
-補剛材あり片持ちファイルは3lxxxh.*
FEM解析値では、片持ち梁で
,length,要素数,saikasitaたわみ平均値(m)
,14,946976,9.66198E-03
*補剛材なし [#r9d897b0]
**補剛材なし3点曲げ解析 [#rf63141a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-補剛材なし3点曲げファイルは3tlxxx.inp (xxxはlength)
,length,要素数(C3D4),saikasitaたわみ平均値(m)
,100,56989,8.17722E-03
,40,276870,8.50203E-03
,20,1489814,8.38163E-03
**補剛材なし片持ち [#pa8c0825]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-補剛材なし片持ちファイルは3lxxx.* (xxxはlengthの数)
,length,要素数,片持ちsaikasitaたわみ平均(m),荷重(P=1.d-1 ...
,20,717448,1.15251E-02
,14,920757,1.15435E-02
*鋼板欠け [#ibb126e0]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten...
*三角孔部分 [#ad299fb8]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-完成
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-作成中(14/5/9)。より正確にモデリングする。
--各種寸法から直角三角形の角度θを求めて頂点からの距離xだけ
離れた点に半径20mmの円が触れるようにする。
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
-作成中(14/5/8)。
--なんとなく概算で作る。円の半径は20mm
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kak...
*滝田さんのプログラム [#jd30da7e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/taki...
*Mentatメモ [#vde19d27]
***弾塑性解析 [#se5a32fb]
-降伏応力を入れた後、
JOBS→PROPERTIES→JOB RESULTSで
Sterss,
Total Strain,
Elastic Strain,
Plastic Strainにチェック。
***面荷重(Faceload) [#s8b041f7]
-格子状に分割されたある大きな一つの面に対し、面荷重をかけ...
Face Loadを選択し、
ボックスの中に応力を入れることになる。
-この時の応力は、かけたい総荷重/断面積だが、
正確には、
$\frac{かけたい総荷重}{faceの数} \div \frac{断面積}{face...
となる。
-答えは変わらないけど、面はたくさんあるのに応力入れるとこ...
焦って考えて、変なことになることがありそうなので注意。
***inpからABAQUSで読む [#j0b3bcb3]
-inpをインポートする際はinpのヤング率の乗数をEで出力する...
mentatでヤング率を再設定する手間が省ける。
*ccx_2.5調査 [#u510ef96]
**DLOAD検証 [#j3ea62b7]
***DLOAD例 [#p9716c77]
&link(cgxで載荷拘束する方法,http://www.bconverged.com/con...
を見ながら手順通りにやった例。
-inpファイル&link(beam.inp,http://www.str.ce.akita-u.ac.j...
-cgxで作るモデル &link(all.msh,http://www.str.ce.akita-u....
-拘束条件&link(fix1.nam,http://www.str.ce.akita-u.ac.jp/~...
-載荷条件&link(load.dlo,http://www.str.ce.akita-u.ac.jp/~...
-frdファイル&link(beam.frd,http://www.str.ce.akita-u.ac.j...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload.png
***鋼材一本で分割数変えてDLOAD [#g47b86de]
-&link(dlokou4.f90,http://www.str.ce.akita-u.ac.jp/~gotou...
--分割数nz=100の時読み込む&link(kou100.dlo,http://www.str...
--分割数nz=500の時読み込む&link(kou500.dlo,http://www.str...
5cm$\times$5cmのスパン50cmの鋼材を片持ちでDLOAD
めんどいので端一面にkousoku,1,3を与えて.dloだけcgxで吐き...
-ある荷重をかけたいなら、pres hoge に、「設定したい荷重/...
,nz分割数,100,500
,FEMたわみ最大値(m),-3.625050e-03,-3.641140e-03
-presを変えていないので、z方向の分割数を変えてもたわみは...
***CLOADとDLOAD比較(鋼材、引張) [#ha1069cc]
1m$\times$1m、長さ50mの鋼材で比較。
-nxとnyは分割なし（1分割)。nzのみ分割(nz=200)
--DLOAD min:-2.425040e-04 at node:801 &link(dlokou4.f90,...
--CLOAD min:-2.425040e-04 at node:801
DLOADにあわせて、固定してない方の端の節点4つに、1つ1/4MN...
-nx=2,ny=1,nz=200 (pres 1.0)
--DLOAD min:-2.425090e-04 at node:1201
--CLOAD min:-2.425090e-04 at node:1201
-nx=1,ny=2,nz=200 (pres 1.0)
--DLOAD min:-2.425090e-04 at node:1201
--CLOAD min:-2.425090e-04 at node:1201
-nx=14,ny=14,nz=500 (pres 1.0)
--DLOAD min:-2.425610e-04 at node:112501
--CLOAD min:-2.425610e-04 at node:112501 &link(ckounx14...
--理論式のたわみ$d=\frac{PL}{EA}$ 2.427184e-04 変位の...
-CLOADでは、分配された節点荷重が重なるところは足し算する...
全く同じ結果になった。このやり方では、節点数が変わってし...
***CLOADとDLOAD比較(木材、片持ち梁等分布) [#aec8bd06]
10cm$\times$10cm 長さ1mの木材を片持ち等分布荷重で解析。
DLOADではpres 0.1とする。
理論値では$k=\frac{5}{6}$として計算を行う。
-nx=1,ny=1,nz=200
--DLOAD min:-2.026590e-02 at node:803
--CLOAD min:-2.026590e-02 at node:803
-nx=14,ny=14,nz=500
--DLOAD min:-2.066020e-02 at node:112711 &link(dmoku500....
--CLOAD min:-2.066020e-02 at node:112711 &link(cmoku14.p...
-ティモ片持ち等分布単位荷重(2kGA) 2.066545E-002
-ティモ片持ち等分布境界値(kGA) 2.183512E-002
幅8cm$\times$高さ4cm 長さ1mの木材 片持ち等分布荷重(pres0....
理論値では$k=\frac{5}{6}$として計算を行う。
--DLOAD min:-9.210080e-02 at node:32476 , &link(png,http...
--ティモ片持ち等分布単位荷重(2kGA) 9.226345E-002 , &link(...
***DLOAD現状 [#qde8aec6]
等分布荷重の時に便利かも。cgx経由じゃないと使えないため、...
cgxで載荷拘束を設定する場合は、
拘束節点を.nam、載荷面情報を.dloに出してそれぞれ.inpでinp...
DLOADにするとできる。
-プログラムで拘束だけ与えて載荷だけcgxの.dloをインプット...
-普通のCLOADで等分布したときとDLOADとで比べてDLOADが
きれいに応力分布出るのか確認。←今のところおんなじになると...
-等分布したい時なら、プログラムに慣れてない人でもcgxで視...
-presの値の設定（presの意味)は間違えないでね。
-後藤ちゃちゃ
DLOADは、&link(図114、115みたいな感じで,http://www.bconve...
**片持ち梁とティモシェンコ理論値比較 [#de1d34be]
簡単な片持ち梁の長方形断面を解いてfrdやdatで見てみたらち...
-ifの練習モデルに載荷した&link(ifyousos.f90,http://www.st...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/if-tawami...
--y軸:たわみ(m)、x軸:nz分割数、赤:ccx、緑:ティモシェンコ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ccxtyousa...
--y軸:相対誤差、x軸:nz分割数
,ティモシェンコ理論値,1.52391E-003(m)
**要素の調査 C3D6とC3D8 [#kc919dcb]
卒論日誌にもっと詳しく書いてたんだけど、間違えて消してし...
直方体でメッシュを生成するプログラムから、その直方体を2分...
-c3d6はこんな感じで作られる
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/cgxc...
-木材でc3d6とc3d8を比較,理論値は初等梁とティモ
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/c3d6...
--y軸:たわみ(m),x軸:nz分割数
-鋼材でc3d6とc3d8を比較,理論値は初等梁
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/kout...
--y軸:たわみ(m),x軸:nz分割数
***雑感 [#nc537c35]
nz分割数が少ない時は(nz<50くらい)、単純にc3d6を作るプログ...
２倍になるので(直方体を2分割してc3d6をつくってるから）、c...
しかし、nzが大きくなり、解が収束してくると(nz>500)、ごく...
*ccx_2.5 ,cgx_2.5覚書 [#ba5c6862]
**RESTART [#d09567e6]
柴田さんの修論日誌に書いてあるとおりにやると、ccx2.5でも....
大体こんな感じで*STEPと*STATICの間あたりに*RESTARTを挿入...
--(
*STEP [#e9412fb8]
*RESTART [#ec504b7c]
*RESTART,WRITE,FREQUENCY=1 [#p14d1e17]
*STATIC,SOLVER=SPOOLES [#g0f640d1]
--)
**座屈解析 [#g753cdc2]
-テクニック
総荷重を1MNかけるときなら、datに出てきた座屈荷重をそのま...
**エラー [#v8b95f9f]
**jacobian問題 : ERROR in e_c3d: nonpositive jacobian [#q...
&link(江村さんの修論日誌,http://www.str.ce.akita-u.ac.jp/...
節点座標が出力フォーマットより細かった時に、その節点を使...
そんな時は出力を1pe17.9とかにして細かい桁まで出してあげる...
***SalomeでAutoでメッシュを切った時に発生 [#nf5aaac5]
Salomeで、プレストレス木箱桁橋三角孔補剛材ありモデルを、
メッシュを全てC3D4で切り、非常に細かく分割した後unv→inp変...
要素数は400万くらい。
要素数70万のものではエラーは発生しなかった。
計算できないくらい細かくなりすぎた？ →細かくなり出力桁が...
--(
*ERROR in e_c3d: nonpositive jacobian
determinant in element 4518956
--)
SalomeでAutoでメッシュを切り、
C3D6とC3D4を混在させたモデルで、C3D4部分を細かくしようと...
全く同じjacobianエラーが発生する。
試しに1つのjacobianエラーが発生したinpを、要素の節点番号...
unv→inp変換プログラムの節点出力を1pe17.9とかにして細かい...
***プログラムで発生 [#s1951dbe]
--(
*ERROR in e_c3d: nonpositive jacobian
determinant in element 20358
*ERROR in e_c3d: nonpositive jacobian
determinant in element 20358
--)
C3D6のところでひとつの要素につき2つのjacobianエラーが出る。
こちらの方は、要素のくるみ方を正しく？(ccx_2.5で対応でき...
ccx_2.5のfortranでの正しい？要素のくるみ方はマニュアルに...
ccx2.6以降では、この問題が解決している可能性があり、
それらのverではccx_2.5でエラーが起きないくるみ方をしてい...
逆にエラーが出る可能性がある。
**f90のコンパイルエラー [#r49f170e]
gfortranでコンパイルを行う時、
--(
relocation truncated to fit: R_X86_64_32S against `.bss'
collect2: ld はステータス 1 で終了しました
--)
みたいなエラーが出てコンパイルできないことがある。
この時は、コンパイルに制限値以上のメモリを使用するために...
なので、制限を取っ払うオプションである-mcmodel=largeをつ...
ダイヤカット円筒をC3D4要素で一気にinp作るプログラムで、要...
--(
gfortran -mcmodel=large nantoka.f90 -o nantoka
--)
**強制終了 [#n2519c86]
ccxで解析中に、ccxが強制終了してしまう。
要素数が多いと発生。
--(
Factoring the system of equations using the symmetric sp...
Using up to 1 cpu(s) for spooles.
強制終了
--)
要するにスペック不足ですってこと。いいPCでやっても出るな...
**ERROR in u_calloc: error allocating memory [#pfe89c4b]
メモリが足りませんよって言われる。
--(
*ERROR in u_calloc: error allocating memory [#k243c6da]
num=566129320,size=4
--)
visionのメモリ2Gのマシンで、140万要素のものを解こうとした...
Corei5 のメモリ4GBのマシンでは、エラーが出ずに普通に解け...
要するにスペック不足ですってこと。
**メモ [#k4a5ceb9]
--.inpをcgxで開いてanimateとかf2とかでたわみ方を確認する...
frdとdatを信用しよう。
--f90で寸法をm、ヤング率と荷重はメガ
-.dloのpresもメガ
**要素 [#p6dbdcd4]
-C3D8 直方体
--C3D20 直方体(C3D8)を一辺に3つ節点ができるように増やした...
-C3D4 四面体
--C3D10 四面体(C3D4)を一辺に3つ以下略
-C3D6 三角柱
--C3D15 三角柱(C3D6)を以下略
-S6,S8 シェル
-他にもいっぱいあるよ。自分で調べて使ってあそぼう
-Rがついてるのは低減積分要素。制御が難しいので、簡単に使...
*Salome覚書 [#f915456c]
**Salome用参考資料 [#ifffd4df]
-チュートリアル的な感じでわかりやすい、超初心者向け
http://tyloroqtism.blogspot.jp/2014/01/salome.html
-salomeの所の解説にある、メッシュの切り方(分割数の指定 な...
http://www.geocities.jp/penguinitis2002/
-いろんなチュートリアル、よく読めば発見があるかも
http://caelinux.org/wiki/index.php/Contrib:Claws/Code_Aster
-マニュアル的なpdf
http://framabook.org/docs/Code_Aster/beginning_with_Code_...
**マウスとキーボードを使った操作方法 [#mf12ea58]
,Ctrl+左クリックしながらドラッグ ,拡大＆縮小
,Ctrl+真ん中クリックしながらドラッグ ,移されている領域を...
,Ctrl+右クリック ,モデル回転
**salomeで.unvをつくりプログラムで.inpにする方法(解決) [#...
***&link(Salomeメモ,http://www.str.ce.akita-u.ac.jp/cgi-b...
材料が2つの時もSalomeメモに書いてます。
**体積を出す [#x603c285]
Geometryで
**unv→inp用の載荷拘束テク [#deb31ce0]
-Salomeで、C3D4などでメッシュを切るとき有効
-きれいに載荷節点、拘束節点をいれたいとき有効
unv→inp変換する際は、Salome上で先に載荷拘束用の区切りをPa...
きれいに載荷拘束用の節点が入る。
unvにするf90では、
if文で
--(
区切りを入れた部分の座標-0.000001d0.and.区切りを入れた部...
--)
みたいに設定していくだけできれいに載荷拘束が入る。
手動で(CLOADで)面載荷したいの時も使える。
Partition入れまくったモデルを、Auto Hex(C3D8)とかでいっぺ...
当然エラーが出るので注意。
そういう意味ではC3D4は便利。(いっぺんに全部切っても大丈夫)
**unv→inp変換時の注意 prism [#q81eb19f]
salomeでメッシュを切った時にprismができるとcalculixでは対...
prismができないようにメッシュを切らなければいけない。
prismは例えばSalomeでAutoで切ったときに、C3D4メッシュとC3...
**メッシュをC3D8要素とかで切っていて、軸方向に多く切りた...
-submeshを使う。
-Geometryで軸方向(多く切りたい方向)の辺すべて(例えば1個の...
-まずBox_1(作りたいモデル)を作る。
-New Entity-Group-Createをクリックし、ラジオボタンで/(辺)...
-Main Shapeには、Box_1を選択。するとBox_1がワイヤーフレー...
-軸方向(z方向)に伸びている辺をクリックし、Add をクリック。
-これを軸方向(z方向)に伸びている辺すべて行う。一つの直方...
4つの辺番号が表示される。
-すべてAddしたら、Apply and Close 。この操作でできるGroup...
-Mesh画面に行き、&link(メッシュをC3D8で切って,http://www....
-その後、create submeshで、MeshにMesh_1を指定、Geometryに...
1D-AlgorithmでWire Distcretisation,HypothesisでNoneを指定...
-すると、新しいNb.Segmentsを作れるので、Number of Segment...
(デフォルトでは、名前が被らないようになっていると思うが、...
-設定し終わったら、Applyする。するとツリーのMesh_1とある...
-Mesh_1を再びcomputeすると、submeshの設定が反映されて再び...
これを行うと、直方体とかなら普通にf90プログラムでnx=,ny=,...
設定するのと同じような感覚でメッシュが切れるし、軸方向分...
unvからinpにしてccxで解いても、寸法・各種数値と分割数が全...
**材料AとBからなるモデルを作りたいが、材料AとBの間をまた...
材料AとBの接合面を、Partitionで区切る。
Partitionでは接合面となるFaceがGeometryにある必要があるが...
ツリーの材料Aか材料Bを右クリックして、CreateGroup→ラジオ...
接合面をクリック→Addで接合面を作るのが効率が良い。
**非一体化モデルを作りたい [#ke26e2f0]
Geometryで、FuseのかわりにCompoundを使ってモデリングする。
Compoundは、Fuseとは違いMeshを切った時に一体化させた境界...
ただし、Compoundした後に境界面をPartitionすると情報が上書...
その境目では節点が1ずつ入る（→つまり一体化する)ようだ。
-compound
選択した物体(複数可)をひとつのグループにまとめてくれる。
一体化してくれるわけではないので注意。
物体Aが物体Bにめり込んだ状態で2つをcompoundしても、
物体Aと物体Bが結合しないし、別々に物体Aと物体Bがあるとい...
**メッシュを細かく切った時エラーみたいなのが出る [#z68748...
elementsが500000(桁があってるか不明)とかより多いよ
みたいなエラーが出たら、Mesh画面でFile→Preferences→Mesh→G...
Mesh informationにあるAutomatic nodes compute limit を増...
ただしエラーが出ていても分割は終了しているっぽい(?)のでそ...
私はlimitを9999999にしている。
**グルーピングが遅い、もっと早くできないの？ [#y84ecb45]
普通は、Create Meshでメッシュを切り、グルーピングを行うが、
Meshを切った後、一旦Saveし、Salomeを再起動させてからグル...
普通にグルーピングするより早く終わる。
なんでだろうね。
**グルーピング後、unv→inp変換はうまく行くのに、inpの要素...
丸みを持つモデルとかで発生する。
とりあえず、SalomeのMeshを右クリック→Mesh Informationでvo...
inpファイルで総要素数を見る。
inpファイルの方が要素数が少なかったら、グルーピングで要素...
丸みがある部分を持つモデルなどをC3D4（角がある要素)で作る...
メッシュを切る際に簡略化されてしまう。（角がある要素では...
そのため、三角孔を持つ木橋などは、メッシュを切った後に、
鋼板部分をグルーピングする際はフィルターに三角孔がある鋼...
すべての要素をグルーピングできなくなる。
こういった場合は孔を開ける前の鋼板をフィルターに設定して...
**C3D4とC3D6を混在させて、inpにして解くとjacobianエラーが...
-ccxエラー - jacobian問題 参照
**Submeshで部位ごとにメッシュを切る [#e1693575]
切り分けたい部位ごとに面などでpartitionを入れておき、
最終的なモデル(たとえばPartiton_1とする)
をExplode-Solid し、Partiton_1から小さいツリーで部位ごと...
Mesh画面でPartition_1を選択しCreate Meshするが、Algorithm...
このMesh_1を選択し、Create sub-meshする。
GeometryでPartition_1の下の細かいツリーの中から切りたい部...
それぞれ部位ごとにSubmeshを作ったら、Mesh_1をComputeする...
Submeshごとに優先順位を変えることもできる。(Mesh_1を右ク...
-Submeshをたくさん作ったり、優先順位の設定が複雑だったり...
たまに
unv出力した際にunv内で要素の順番が変になり、
unv-inp変換プログラムを通した際にinpデータの要素番号が負...
priorityの順番に要素グループを作っていけばこのエラーは回...
**二次要素(中間節点がある要素)をつくる [#f34ce899]
mesh画面でconvert to quadratic
二次要素の方が精度がよいので使いましょう。
二次要素をそのままcalculixに流してとこうとしたところ、
節点の並び方が対応していないのでjacobianエラーが大量に出...
cgxでモデルを確認しても要素がかなり抜け落ちていた。
salomeのautoで切った四面体2次要素をccxで解くのは厳しそう...
**エラー [#b1736fd4]
***Warning:It seems that there is not enough memory~ [#te...
メモリが足りませんよエラー。
--(
Warning
It seems that there is not enoght memory to show the mesh
so that the application may crash. Do you wish to continu...
--)
meshをcomputeしたとき、
Vision 2GのマシンでC3D4で(たぶん、およそ)200万要素くらい...
**材料が2つのとき（を考えてる時に考えたこと) [#v86ece3e]
-例えば部材がbox_1の上にbox_2が乗っかっているモデルだとす...
-box_1とbox_2の境界面をpartitionで区切らなければいけない...
-ツリーのbox_1を右クリックし、Create Groupを選ぶ。
-ラジオボタンで、□(面)を選び、Main Shapeに面を指定する物...
-Geometryにbox_1だけを表示させて、box_2と接している面をク...
-これでpartitionのための設定が終わる。
-box_1とbox_2を一体化させてメッシュを切るので、ここからFu...
-Operations-Boolean-fuse で、object 1 にbox_1 、object 2...
-Fuse1をpartitionで区切る。Operations-Partition を選び、O...
Tool Objectsには、ツリーのbox_1の下できているGroup_1を選...
-Partition_1ができる。
-Mesh画面にいき、Partition_1を選択してメッシュを切る。(こ...
-メッシュを切ったら、グルーピングを行う。Fuseしたモデルを...
-MeshにMesh_1、Element Type にVolume、Nameは(かぶらなけれ...
-Group type ではGroup on filterにチェックを入れる。すると...
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-今回は、Geometry のモデルでフィルターをかけるので、Crite...
選び、Compare にはEqual to を選ぶ。
-そして、Threshold の中の空白をクリックしたら、ツリーの中...
-これで、box_1の部分のメッシュがすべてGroup_1mに入れられ...
-もう一度Create Groupする。今度はbox_2の部分のメッシュを...
-MeshにMesh_1、Element Type にVolume、NameはGroup_2mとか...
-Group type ではGroup on filterにチェックを入れる。すると...
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-Criterion にBelong to Geom を選び、Compare にはEqual to ...
-そして、Threshold の中の空白をクリックしたら、ツリーの中...
-これで、box_2の部分のメッシュがすべてGroup_2mに入れられ...
-Mesh_1を右クリックして、unvファイルを保存する。
メッシュをC3D4で切り、2つにグループ分けしたものなら、ccxc...
partitionを使わない方法だと、C3D4とかだとbox_1とbox_2の間...
となると、メッシュの切り方を変えてBox_1とBox_2で別々に切...
Compound Meshを作ってその後グループ分け、という作業が必要...
更に、Compoundするならその接合面で節点がくっつくほうがい...
Box_1の上にBox_2が乗っかってるモデルとするなら、FaceのBox...
という手順が必要かもしれない。調査中。
partition使ってないからこんなことになるのか？partition使...
-Compoundとか使うとunvの要素の入り方が変わってくる？。普...
最初が-864とかから始まる。
-partitionを使う方法
予め面を作っておき、partitionを使って区切ってからメッシュ...
下の３角孔木箱桁橋でやったら要素の抜けが起こった。鋼板と...
partitionで区切ってるのになあ…。(鋼板の方の面をグループ分...
→解決しました。フィルターの範囲を大きくしなければいけませ...
*パソコン関係覚書 [#e11ccfff]
**viコマンド [#x4db6525]
viのコマンド、特に移動系のコマンドなどを使いこなせてない...
-移動系(0と$は特に便利)
,0 ,その行のはじめに移動
,$ ,その行の最後に移動
,Ctrl e ,一行下へスクロール
,Ctrl y ,一行上へスクロール
,Ctrl f ,1ページ分進む
,Ctrl b ,1ページ分戻る
,H ,表示されている画面の一番上の行に移動
,M ,表示されている画面の中央の行に移動
,L ,表示されている画面の最後の行に移動
-入力系
,s ,カーソルの位置の文字を削除して入力モードに
,Ctrl s ,キー入力無効状態になる
,Ctrl q ,キー入力無効を解除
**aliasでコマンド短縮 [#l466b40b]
プログラムをコンパイルする際、gfortranといちいち打つのが...
~/.bashrcの中で# some more ls aliasesの後辺りに
--(
alias gf='gfortran -o'
--)
と追加しておくとgf がgfortran -o の短縮コマンドになる。た...
**前に使ったコマンドを見る [#u4a4ae2c]
--(
history|grep nantoka
--)
でnantokaというコマンドを使った履歴を検索できる。
番号とその時使ったコマンドがずらっとでて、いちいち打ち直...
**inpファイルを見やすくする [#s9042b53]
&link(vimrc_ubuntu,http://www.str.ce.akita-u.ac.jp/~gotou...
**メモリとかCPUの情報を見る [#zdf21144]
--(
more /proc/cpuinfo
more /proc/meminfo
--)
**外付けHDにバックアップをとる [#td205260]
--(
rsync -auv コピー元 コピー先
--)
**jobの最後にrm jobを入れとく [#racb5007]
jobの使い方は後藤資料の中でぐぐればあります。で、jobの最...
--(
rm job
--)
を入れとくと、!./を使ったときjobが再実行されて、慌ててCtr...
解いた結果のdatファイルが全部真っ白になっているということ...
**libre office [#mbf46ca9]
2014/4/21現在、Ubuntu12.04環境では、
Libre Officeが入っているわけだけど、どういうわけかすごく...
スライド作成などしたい場合には、事前に
verを3.なんたらから4.2.3にバージョンアップしたほうがいい...
コマンドはlibre office ubuntu12.04 とかでぐぐると出てく...
**役立つページ [#o35865fd]
-http://www.str.ce.akita-u.ac.jp/~chida/memo.html#i0
-http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%...
*gnuplotメモ [#wcc133bf]
**gnuplotで値をグラフにして見たいって人は、 [#k088c10c]
-まず、txtデータを作る。
hoge.txtとか(何でもいい) にx軸データ y軸データ を入れる。
hoge.txtの中身が、
--(
1 50 ←1列目がx軸データ、空白を開けたあと、2列目をy軸...
100 5000
--)
なら、
gnuplot で、
--(
pl 'hoge.txt' w l
--)
とすると、グラフ画像が出てきて、見れる。
**グラフをpngにして保存したいって人は、 [#u7fb64b4]
-まずpngを作る
--(
set term png
set output 'nantoka.png' ←nantokaはグラフにつけたいファイ...
--)
-そんで、グラフ出力する
--(
pl 'hoge.txt' w l
同時にいくつもグラフを載せたい場合は、
pl 'hoge.txt' w l , 'moge' w l , 'mage' w l
--)
set outputしたあとに、plotすると、グラフが表示されずに、....
この2つの操作でグラフを.pngにして保存できる。
**1つのテキストファイルからいくつもグラフを書きたい！って...
hoge.txtの中身を、
--(
1 50 100 0.2
100 5000 10000 20
--)
というふうにデータを並べてって、gnuplotで、
--(
pl 'hoge.txt' using 1:2 w l ,'hoge.txt' using 1:3 w l ,'h...
--)
みたいにしてusingを使う。
using1:2 で、txtデータの1列目がxに入り、2列目がyに入るっ...
**線の太さをもっと太くしたい！って人は、 [#y04eee9e]
gnuplotでそのままグラフを作成すると、
スライド作成とかの時、グラフの線が細くて見にくいと言われ...
グラフの線を太くするには、linewidth を使う。 省略形でlwも...
lw 2みたいに、その後に半角スペースを空けて、整数を入れる...
--(
pl 'hoge.txt' w l lw 2
--)
みたいにすると、太くなる。(だいたい2が見やすいと思う）
**グラフの線のタイトルを消してクールにしたい！って人は、 ...
notitleを使う。
発表の時は、スライド作成ソフトやxfigとかで、グラフに矢印...
そういう時には右上に自動生成されるグラフの線のタイトルが...
--(
pl 'hoge.txt' w l notitle
--)
*構造実験検証 [#x67cf198]
一人で（龍さんに手伝ってもらったりして）こっそりやったデ...
-4/24午後7時 二回とも成功(変な測定値は出なかった)
--2回目データ:&link(140424a_2s.csv,http://www.str.ce.akit...
--1回目データ:&link(140424a_1s.csv,http://www.str.ce.akit...
*3D [#d2eb7ba5]
**引張試験 試験体 [#h211cd60]
引張試験のための供試体を作成する。
ものづくりセンターの3Dプリンタでは造形情報をstlファイルと...
SalomeのGeometryでモデルを選択→export→STLAⅡfilesで吐き出...
**引張試験 [#f5580cac]
***メモ [#a44cbbb9]
-ロガーにつながってるPcのPCD-30Aの集録データ数 の欄で何秒...
-doループの終わりをデータ終わりの秒数にするだけでヘッダー...
*回転らせん [#eeb2f178]
遊びで適当に節点だけ打ってみたけど、
中間節点を作るのが難しそう
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpn...
*プログラム類↓ [#wa71b442]
&link(卒論日誌,http://www.str.ce.akita-u.ac.jp/cgi-bin/gw...
*プレストレス木箱桁橋の角材の間にゴムを挟めるモデル [#y89...
接触解析の代わりになるかも という事で滝田さんのプログラム...
ゴムの寸法や、分割数、ヤング率(ELgom)とかは各自で変えてく...
節点座標の出力桁を多くして、
ゴムを小さくした時に節点座標が出力桁以下になり、
四捨五入されてしまいメッシュの大きさが均一にならないこと...
&link(gom2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/...
*B32用のunv→inp変換 [#v1d92ed0]
2次元のトラスにしか対応してません。
3次元には非対応。
主な手順
-まずSalomeでyz平面上にトラスを作る
-点を打ってってそれをLineにして
-トラスができたら、全部のLineを選択して1つのWireをつくる
-その後,作ったWireを,
-1Dでメッシュを切るが、number of segmentsで分割数を３以上...
-Meshを.unvに吐き出して&link(b32unv.f90,http://www.str.ce...
*2014年度後期 [#r0745b2b]
,曜日,12,34,56,78
,時間帯,8:50-10:20,10:30-12:00,12:50-14:20,14:30-16:00
,月,,,基礎物理学実験D TA ,TA
,火,,,,
,水,,構造設計学特論,,
,木,,地盤工学特論,工学資源学特論,
,金,,,水防災工学特論,
集中講義 工学資源特論
*2014年度前期 [#z45bf2e6]
,曜日,12,34,56,78
,時間帯,8:50-10:20,10:30-12:00,12:50-14:20,14:30-16:00
,月,,構造力学特論,,パソコンゼミ
,火,,都市システム計画特論,土質工学特論,
,水,,材料設計特論,,
,木,,,交通システム計画特論(英語ゼミ),構造実験
,金,,,,
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
