#author("2020-01-20T11:56:39+09:00","default:kouzouken","kouzouken")
#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/graph/syseich35.png
--ch36
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/syseich36.png
-斜材
--ch33
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/syseich33.png
--ch34
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/syseich34.png

**予備試験 [#g8532b48]

tugite130710yのデータをゼロ点修正


縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch35yobi2.png
--ch36

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch36yobi.png
-斜材
--ch33

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch33yobi.png
--ch34

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch34yobi.png
**両方プロット [#rba2340d]

tugite130710.xlsx、tugite130710yのデータをゼロ点修正


縦軸 荷重(kN)、横軸 ひずみ(μ)、緑が予備試験
-支間中央部に近い垂直材
--ch35
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch35sy2.png
--ch36

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch36sy.png
-斜材
--ch33

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch33sy.png
--ch34

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/ch34sy.png
*目標 [#o5a46862]

修士発表時間
http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%c2%b4%cf%c0%bd%a4%cf%c0%a5%e1%a5%e2#i87

-等圧縮を受ける板の座屈を確認

*板の座屈 [#h2a0e021]

**等圧縮 [#m6184594]
***正方形板 4辺単純支持 [#j5c4a131]

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/kouitass8cfrd.png

&link(kouitass8c.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/kouitass8c.f90)

$\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/scpng/kouitass8cfrd2.png

***長方形板 4辺単純支持 [#s9e4b45f]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/kouitassk9frd.png

&link(kouitassk9.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/kouitassk9.f90)


長方形板だけど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/scpng/kouitassk9frd1zi.png

**純曲げ [#pa5aaff7]
三角形分布でも総荷重を分配で良さそう。
一つの三角形に1/4MNを分配してやる。

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/kouitafrd.png

**参考になりそうな文献 [#o47661f2]

-長方形板の分岐階層構造の境界条件による退化
https://www.jstage.jst.go.jp/article/jscej1984/1995/507/1995_507_65/_pdf

-弾性矩形板の曲げ解析手法の一般化について
http://ci.nii.ac.jp/els/110007887197.pdf?id=ART0009558779&type=pdf&lang=jp&host=cinii&order_no=&ppv_type=0&lang_sw=&no=1442300490&cp=

-等曲げを受ける弾性矩形板に生じる二次座屈現象の解明
https://www.jstage.jst.go.jp/article/jscej1984/1995/519/1995_519_67/_pdf

*Salomeで作ったモデルで座屈 [#e2f56b5a]
**応力集中部を細かく切ったモデル [#od3e26fe]

非一体化にして座屈させるところの鋼板を多く切る。
応力集中部を細かくしていくと1次モードは飛び出しがないモードが得られるが、2次モードで飛び出しが起こる。

鋼板要素数1は非一体化にしてるところと一体化しているところの境目付近で、ここをさらに細かく切ってみる。
鋼板要素数2は残りの鋼板全部。

-孔なし補剛材なし
,総節点数,木材要素数,非一体化鋼板要素数,鋼板要素数1,鋼板要素数2,結果,1次$P_{cr}$,2次$P_{cr}$
,631768,773611(l12),1389462(l1.6),347918(l2),291211(l8),2次で飛び出す,73.8kN,79.7kN

**補剛材なし [#e1ae5037]
-孔なし補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要素数,結果,1次$P_{cr}$
,389725,704146(l12),650860(l2),298061(l8),1次モードは正常、2次モードは要素が飛び出す,83.0kN
,535235,601922(l12),1389462(l1.6),331039(l8),1次モードは正常、2次モードは要素が飛び出す,77.7kN

--プログラムで作ったモデルと近い値が出る。

--1次モードでは座屈荷重77.7kN、2次モードでは座屈荷重78.5kNとなり、
2次モードでは要素が飛び出す。
2次モードでは鋼板の幅方向の変位が1次モードより小さい。


-メモ
--l12-l1-l8だとエラーで切れない

--3孔鋼材非一体化部だけをl1で切ると要素数470万くらいになり多すぎる。


--l10-l1-l8でもエラーで切れない、要素数1000万くらいになり多すぎる。

--l8-l1-l8もエラーで切れない

-孔あり補剛材なし
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要素数,結果,1次$P_{cr}$
,342007,599145(l12),577961(l2),231480(l8),1次モードは正常、2次モードは要素が飛び出す,71.5kN
,522828,1017878(l12),1109490(l1.6),234967(l8),1次モードは正常、2次モードは要素が飛び出す,67.5kN

--l12-l1.6-l8は切れる

--孔をつけたモデルでも解析がうまく行く。孔ありの方が1次の座屈荷重が少し低くなる。


**補剛材あり [#yc7315fa]

-孔なし補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要素数,結果,1次$P_{cr}$
,882260,1798392(l12),1653936(l2),596385(l8),要素が飛び出す,149kN


-孔あり補剛材あり
,総節点数,木材要素数,非一体化部鋼板要素数,一体化部鋼板要素数,結果,1次$P_{cr}$
,523050,598470(l12),1251573(l2),452306(l8),要素が飛び出す,175kN


--孔なしのほうが木材も鋼板も細かく切れているので1次座屈荷重が低くなってる。どちらも要素の飛び出しが起こるモードになり変な解が出てるようなので補剛材と鋼材をさらに切る必要がある。

**今までのをやり直す。 [#c05c0b12]
salomeで木材部を削って一部非一体化にしたモデルを分割してからccxで解析する。

-孔あり補剛材なしsalomeでモデリング
要素はC3D4(四面体)のAutoでlengthを変えて切る

,length,総節点数,木材要素数,鋼板要素数,結果
,12,17万,58万,15万,座屈させると要素が飛び出す。1次$P_{cr}$=92kNくらい
,10,,,,なぜかエラーがでて切れない
,8,50万,203万,36万,要素の飛び出しがなくなる。1次$P_{cr}$=147kN
,7,581146,2396666,386169,要素の飛び出しが見られない。1次$P_{cr}$=144kN
,6,92万,340万,104万,要素数が多くなりすぎてccxで解けない(強制終了)

鋼板をもっと多く切ってみる必要がありそう
2次要素とかで切ってみるとどうか?



-孔なし補剛材なしsalomeでモデリング

,length,総節点数,木材要素数,鋼板要素数,結果
,12,181632,765824,178789,座屈させると要素の飛び出しが起こる。1次$P_{cr}$=105kNくらい
,10,,,,なぜかエラーがでて切れない
,9,,,,エラーが出て切れない
,8,521420,2062159,376681,座屈させると要素の飛び出しが起こる。1次$P_{cr}$=92kNくらい
,7,597914,2412348,450689,座屈させると要素の飛び出しが起こる。1次$P_{cr}$=151kNくらいになる。


*プログラムで作った孔なし非一体化モデル [#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/scpng/811onsi4z2frd1zitou.png

-1次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/811onsi4z2frd1zi2bai.png

-2次モード等倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/811onsi4z2frd2zitou.png
-2次モード変位量2倍
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/811onsi4z2frd2zi2bai.png

**孔なし補剛材なし上の段だけ非一体化モデル [#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/scpng/onsi4uez1zi1bai.png
**孔なし補剛材あり非一体化モデル [#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_Aster_JPAubry_20131206.pdf
の15.2
ver10.8のcommファイルを使う。
これ見ながらcommファイルはオールグリーンにできたけど、postproまでいかない。
たぶんやり方が違うのだろう。
&link(commファイル,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/nari/nari2.comm)


-Asterのマニュアル(公式HP)
http://www.code-aster.org/V2/doc/default/en/index.php?man=U4&lang=en
左上で言語をenglishにできる

-Bucklingについて
--http://www.code-aster.org/V2/doc/v11/en/man_u/u2/u2.08.04.pdf

-MODE_ITER_SIMULTのマニュアル
--http://www.code-aster.org/V2/doc/default/en/man_u/u4/u4.52.03.pdf

MODE_ITER_SIMULTのオプションを使いそう

--lateral buckling of a beam (spill)
http://www.code-aster.org/V2/doc/default/en/man_v/v2/v2.03.504.pdf
どうやって解くかは書いてない

Aster→2材料で解くのと異方性はできる。座屈は調査中

-ccx2.7のマニュアル-*Friction
http://web.mit.edu/calculix_v2.7/CalculiX/ccx_2.7/doc/ccx/node207.html

-ccx2.5のマニュアル*Friction
http://its-cs1.its.uni-kassel.de/CalculiX-2.5/html/ccx/node198.html

*弾性床上の座屈 [#uf9a40dc]

http://freecaetester.blog62.fc2.com/category19-0.html

上の記事を参考にして一部非一体化にしたプログラムに接触解析の定義をいれてみる。

-接触入れるとWARNINGが出る。
--(
 *WARNING in statics: a nonlinear geometric analysis is requested
          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からアンカープレート100mmの半分を引いたもの=315mm

幅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/scpng/kouhanzfrd.png

-スパンを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/scpng/kouhan630frd.png


**木材と鋼板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/scpng/mokukou44frd.png
接触解析にしてないのでオイラー座屈、ccxで鋼板1枚を座屈させた時の座屈荷重 とほぼ変わらない値が出る。

**オイラー座屈応力 [#a3068bda]
-&link(zakutu,http://www.str.ce.akita-u.ac.jp/~gotou/pas/iinkai/deeta.pdf)(akusesu seigen ari)→100kNでは座屈していない
**1/4解析モデル [#y4702bc7]
まず普通の1/4解析モデルを作る。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakox.png
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakoy.png

どちらも下の 半解析-箱断面三角孔モデル のx変位、y変位と同じような分布になり、
たわみ平均なども大体同じ値になる。
*鋼板の座屈 - 一体化だけど一部だけ鋼板と木材の間を削ったモデル [#g10554e8]
とりえあず、中央部の木材と鋼板がくっついている部分の、木材部分を0.001mmだけ削り、その部分だけ木材と鋼板が非一体化しているとみなして計算してみる。

実物は上弦材も下側木材もPC鋼棒間で非一体化なので、FEMモデルでは上弦材と下側木材両方の中央部PC鋼棒間木材部分を0.001mmだけ削る。


**三角孔あり補剛材なし [#df7f9762]

-面載荷しているモデル。length12。
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakozakutu2x.png
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakozakutu2y.png

-木材部分に線載荷したモデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanzakutuanime.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanzakutuL2all.png

***座屈荷重(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の計算値とオイラー座屈の固定支持の値がかなり近い結果になった。

-モード1
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakozakutu2m1.png

何かが飛び出す現象がみられる。
アニメーションで詳しく見てみる。

--モード1アニメーション、Amplitude+100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanza100anime.png

--モード1アニメーション、Amplitude-100%
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakoza-100anime.png

この2つのモード1アニメーションを見る限り木材部ではなく鋼板部が飛び出していると考えられる。この挙動は線載荷と面載荷のどちらでも同様だった。

-モード2
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/4hanhakozakutu2m2.png

**三角孔あり補剛材あり [#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/scpng/4hanhakohogouza1g.png

・補剛材なしとありでは飛び出す方向が違う
→単に方向が違うだけなので補剛材のあるなしでモードが変わるかも
・またはメッシュ分割の違い
(鋼板部メッシュ数は補剛あり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]
オイラーの座屈荷重による計算結果と解析による補剛材なしのモード2の値が近いことから、
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/~gotouhan/saito/2014/ortho/2mnz100.hdf)

--&link(commファイル,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/ortho/o13.comm)

寸法をm単位で入力してるので注意。
このやり方で木材を短くしたり長くしたりして解析したが、
ほぼすべてsalomeとccxの解が同じになった。

-&link(マニュアル,http://www.code-aster.org/V2/doc/v11/fr/man_u/u4/u4.43.01.pdf)
のP13,14にLTN系の座標と例が書いてある。

-Asterでの異方性の設定について書いてあるもの
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04.131.pdf
--http://www.code-aster.org/V2/doc/v10/en/man_v/v3/v3.04.120.pdf
--http://www.code-aster.org/forum2/viewtopic.php?id=4026
--http://www.code-aster.org/V2/doc/v11/fr/man_u/u4/u4.43.01.pdf
の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/scpng/tugitehan.png

-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/tugitehanx.png

-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/tugitehany.png


継手があると支間中央部の鋼板の三角孔が継手鋼板で塞がれるのでせん断変形しづらくなる

*木材利用研究発表会 [#qc0a738b]

2014年 &link(発表スライド,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/pdf/saitomoku14.pdf)
*継手なしモデルまとめ [#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/graph/grafi-3anasaikamen.png


*3点曲げ木材上部分載荷 [#s252a518]
-全てC3D4要素
-支間中央木材部分に線載荷
-載荷上たわみ平均は、線載荷部からz軸方向に5mm離れたところの点のたわみを取る。
**補剛材ありモデル [#od9f0ef7]

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/3tenhogou5mml20mgsaikafrd.png

,載荷上たわみ平均(mm)
,9.98675
**箱断面モデル [#sc88916e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/3tenhako5mml18mgfrd.png
,載荷上たわみ平均(mm)
,11.919

**鋼板上部欠けモデル [#f148c60a]
箱断面モデルに対して上部60mm部分の鋼板がないモデル

-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/3tenkakel18mgfrdd2.png

-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/3tenkakel18mgsaikafrdd1.png

*半解析 [#lae9b1ad]
半解析で実験と同じく面に載荷してみる。
**補剛材あり三角孔モデル [#be3d440a]
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hogouxmen.png

-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hogoumenhenkeizu.png
**箱断面三角孔モデル [#z2b67ee4]
-載荷面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hanhakoload14_8mgsaikameninp.png

-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hanhakoload14_8mgsaikafrd.png

-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hakoxmen.png

-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/hanhakoloadhenkeizu.png

**箱断面(補剛材なし)孔なしモデル [#ge8b74b8]

-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/onsi14hanyhenni.png

-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/onsi14hanxheni.png

-Animation
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/onsi14hananimete.png

-プログラムから作ったモデルと、Salomeから作ったモデルでは、荷重をかける方向(主にy方向)が逆になっているので注意。
-プログラムでは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/scpng/on6m3tentaikfrd.png
-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/on6m3tentaikx.png
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/on6m3tentaiky.png

**対傾構あり三角孔あり半解析モデル [#a0fd5c34]
R20
,木材下部たわみ,6.42216mm

-x変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/sankaku2011anataikxfrd.png
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/sankaku2011anataik.png

**対傾構なし三角孔あり半解析モデル(比較用) [#f4e8ea2a]
R20
,木材下部たわみ,5.64861mm
-x変位
対傾構があるのと比べ、支間中央部の鋼板のx変位が大きい
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/sankaku2011anaxfrd.png
-y変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/sankaku2011anayfrd.png


*C3D6とC3D4の混在 [#fa4a71a2]
メッシュの接合性の問題で、C3D8とC3D4はうまく行かないけど、
C3D6とC3D4ならAutoで切っても混在させれそう。


-木材部分C3D6,鋼板部分C3D4
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/R20katac3d6frd.png

*三角孔片持ち修正版 [#v28ed3fb]
スパンを短くした。したにある古いものはスパンが端部分だけ長いので、それを修正。
木材部分上部に線載荷し、たわみは木材部分下部の平均を取ったものをsaikasita、木材部分上部の平均をとったものをsaikaueとする。3点曲げの値とも比較してみる。

孔ありは全てC3D4要素、孔なしは全てC3D8要素とする。

-面載荷木材部分上部のたわみと実験値
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/graph/FEMmensaikaue-zikken.png

*三角孔補剛材あり [#e7059fff]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m),面載荷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

-3点曲げ
,length,要素数,saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03

-孔なし
,断面,saikasitaたわみ平均(m),saikaueたわみ平均(m)
,補剛材あり,6.39474E-03,8.87783E-03

*三角孔補剛材なし(ただの箱断面) [#d8c7fc64]
-片持ち
,length,要素数,saikasitaたわみ平均(m),saikaueたわみ平均(m),面載荷saikasitaたわみ平均(m),面載荷saikaueたわみ平均(m)
,20,80万,8.41664E-03,1.18110E-02,8.10568E-03,1.26531E-02

-3点曲げ
,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) [#reeba5f1]
*ティモ孔の鋼板の換算幅算出 [#sfea2d8e]
$(元々の鋼板体積-孔部分体積)/(高さ×長さ)$
-三角孔(大きいタイプ)
,元々の鋼板の幅b, 9.0mm
,ならした鋼板の幅b2, 7.68123mm
-補剛材あり用ティモ孔理論値計算f90&link(keitimohana.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/data/keitimohana.f90)
bankの修正せん断補正係数では、補剛材部分を考慮していない。

-補剛材なし箱断面用ティモ孔理論値計算f90&link(keianahako.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/data/keianahako.f90)
今までは孔を簡略化して(長方形孔)モデル化、ティモ(孔)の理論値計算をしていたため、
ならした鋼板はもともとよりかなり薄くなるので、
ティモ(孔)と元々のティモシェンコ理論値は開きが出ていたが、
図面の三角孔の面積を使ってならした鋼板の厚さを使うと、
長方形孔をならしたときほどは薄くならない。

そのためティモ(孔)にしても、
元々の値とはあまり変わらなくなる。

-昔の長方形孔
縦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/3tenhogou125kousoku.png
-補剛材あり3点曲げファイルは3thlxxx.inp (xxxはlength)

-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3tenhogou125kousokudmen.png


,length,要素数(C3D4),saikasitaたわみ平均値(m)
,25,1262482,6.80896E-03

**補剛材あり片持ち [#na7dec11]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/R20katahogou.png
-補剛材あり片持ちファイルは3lxxxh.*

FEM解析値では、片持ち梁で
,length,要素数,saikasitaたわみ平均値(m)
,14,946976,9.66198E-03

*補剛材なし [#r9d897b0]
**補剛材なし3点曲げ解析 [#rf63141a]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3ten3kakuR20kousoku.png
-補剛材なし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/3kakur20katamesh.png
-補剛材なし片持ちファイルは3lxxx.* (xxxはlengthの数)

,length,要素数,片持ちsaikasitaたわみ平均(m),荷重(P=1.d-1 MN) /2
,20,717448,1.15251E-02
,14,920757,1.15435E-02

*鋼板欠け [#ibb126e0]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3tenkouhankake.png

-断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3tenkouhankakedanmen.png
*三角孔部分 [#ad299fb8]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kakuR20-kansei2.png
-完成

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kakuR20-sakusei2-70.png

-作成中(14/5/9)。より正確にモデリングする。
--各種寸法から直角三角形の角度θを求めて頂点からの距離xだけ
離れた点に半径20mmの円が触れるようにする。

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kakuR20-2sakuseityu2.png

-作成中(14/5/8)。
--なんとなく概算で作る。円の半径は20mm
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/3kakuR20.png

*滝田さんのプログラム [#jd30da7e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/takita3kakuinp.png

*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/content/calculix/doc/GettingStarted.pdf)
を見ながら手順通りにやった例。
-inpファイル&link(beam.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/beam.inp)
-cgxで作るモデル &link(all.msh,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/all.msh)
-拘束条件&link(fix1.nam,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/fix1.nam),&link(fix2.nam,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/fix2.nam),&link(fix3.nam,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/fix3.nam)
-載荷条件&link(load.dlo,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/load.dlo)
-frdファイル&link(beam.frd,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/beam.frd)
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/~gotouhan/saito/dload/dlokou4.f90)

--分割数nz=100の時読み込む&link(kou100.dlo,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/kou100.dlo)
--分割数nz=500の時読み込む&link(kou500.dlo,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/kou500.dlo)

5cm$\times$5cmのスパン50cmの鋼材を片持ちでDLOAD
めんどいので端一面にkousoku,1,3を与えて.dloだけcgxで吐き出させる。この時設定したpresは1.0(presは応力)
-ある荷重をかけたいなら、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,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/dlokou4.f90),&link(zd.dlo,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/zd.dlo),&link(dlokou4.png,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/dlokou4.png)
--CLOAD  min:-2.425040e-04 at node:801
DLOADにあわせて、固定してない方の端の節点4つに、1つ1/4MNずつかけている。&link(clokou4.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/clokou4.f90),&link(clokou4.png,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/clokou4.png)

-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.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/ckounx14.f90)

--理論式のたわみ$d=\frac{PL}{EA}$ 2.427184e-04  変位の単位は全てm:メートル
-CLOADでは、分配された節点荷重が重なるところは足し算する。ちゃんと割り振ってやると、
全く同じ結果になった。このやり方では、節点数が変わってしまったら面にかかる荷重の総和が変わってしまうため、注意が必要。(分割数を変えたりすると、載荷節点の荷重を変えなければいけないが。下にあるcmoku14.f90なら分割数を変えてもp=の荷重をちゃんと振り分けてくれるようにした。)

***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.png,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/dmoku500.png)
--CLOAD  min:-2.066020e-02 at node:112711 &link(cmoku14.png,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/cmoku14.png), &link(cmoku14.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/cmoku14.f90)
-ティモ片持ち等分布単位荷重(2kGA) 2.066545E-002
-ティモ片持ち等分布境界値(kGA) 2.183512E-002


幅8cm$\times$高さ4cm 長さ1mの木材 片持ち等分布荷重(pres0.03)
理論値では$k=\frac{5}{6}$として計算を行う。
--DLOAD  min:-9.210080e-02 at node:32476 , &link(png,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/dmokus8cm-8cm-1m.png)
--ティモ片持ち等分布単位荷重(2kGA) 9.226345E-002 , &link(keidmokus.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/dload/keidmokus.f90)

***DLOAD現状 [#qde8aec6]
等分布荷重の時に便利かも。cgx経由じゃないと使えないため、モデルの寸法・分割数決め打ちになってしまう。
cgxで載荷拘束を設定する場合は、
拘束節点を.nam、載荷面情報を.dloに出してそれぞれ.inpでinputしてやるように書き、
DLOADにするとできる。
-プログラムで拘束だけ与えて載荷だけcgxの.dloをインプットするやり方が楽。
-普通のCLOADで等分布したときとDLOADとで比べてDLOADが
きれいに応力分布出るのか確認。←今のところおんなじになると思われる。そういった部分でのDLOAD特有の利点は、ない(と私は思う)。
-等分布したい時なら、プログラムに慣れてない人でもcgxで視覚的に載荷条件設定できるのでよいかも?分割を多くしてくとcgxでの載荷節点設定が重くなってこわいのが欠点。
-presの値の設定(presの意味)は間違えないでね。

-後藤ちゃちゃ
DLOADは、&link(図114、115みたいな感じで,http://www.bconverged.com/calculix/doc/ccx/html/node138.html)等価節点力に変換されるということだと思います。DLOADの等価節点力で載荷した方が応力分布が一様になるかも。

**片持ち梁とティモシェンコ理論値比較 [#de1d34be]
簡単な片持ち梁の長方形断面を解いてfrdやdatで見てみたらちゃんと解けているようです。
-ifの練習モデルに載荷した&link(ifyousos.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ifyousos.f90),計算は&link(keiifyouso.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/keiifyouso.f90)
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/if-tawami-timoc.png

--y軸:たわみ(m)、x軸:nz分割数、赤:ccx、緑:ティモシェンコ理論値

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/ccxtyousa/gosa-iftawahz.png

--y軸:相対誤差、x軸:nz分割数

,ティモシェンコ理論値,1.52391E-003(m)

**要素の調査 C3D6とC3D8 [#kc919dcb]
卒論日誌にもっと詳しく書いてたんだけど、間違えて消してしまったので簡略に。
直方体でメッシュを生成するプログラムから、その直方体を2分割して三角柱メッシュを生成するプログラムを作り、収束性を比較する。
-c3d6はこんな感じで作られる
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/cgxc3d6tinp.png

-木材でc3d6とc3d8を比較,理論値は初等梁とティモ
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/c3d6c3d8.png
--y軸:たわみ(m),x軸:nz分割数

-鋼材でc3d6とc3d8を比較,理論値は初等梁
http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/sotu/koutyouc3d6c3d8.png
--y軸:たわみ(m),x軸:nz分割数

***雑感 [#nc537c35]
nz分割数が少ない時は(nz<50くらい)、単純にc3d6を作るプログラムのほうが軸方向の要素数が
2倍になるので(直方体を2分割してc3d6をつくってるから)、c3d6のほうが理論値に近いくれています。
しかし、nzが大きくなり、解が収束してくると(nz>500)、ごくわずかだがc3d8の方が理論値に近くなった。c3d8のほうが1要素あたりの節点数が多いからなのかもしれない。

*ccx_2.5 ,cgx_2.5覚書 [#ba5c6862]
**RESTART [#d09567e6]
柴田さんの修論日誌に書いてあるとおりにやると、ccx2.5でも.routファイルができた。
大体こんな感じで*STEPと*STATICの間あたりに*RESTARTを挿入した。
--(
*STEP [#e9412fb8]
*RESTART [#ec504b7c]
*RESTART,WRITE,FREQUENCY=1 [#p14d1e17]
*STATIC,SOLVER=SPOOLES [#g0f640d1]
--)
**座屈解析 [#g753cdc2]
-テクニック
総荷重を1MNかけるときなら、datに出てきた座屈荷重をそのままMN単位で読める


**エラー [#v8b95f9f]
**jacobian問題 : ERROR in e_c3d: nonpositive jacobian [#q4b8b3dc]

&link(江村さんの修論日誌,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%b9%be%c2%bc%a4%ce%bd%a4%cf%c0%c6%fc%bb%ef#i50)にも書いてくれています。有難いです。

節点座標が出力フォーマットより細かった時に、その節点を使う要素を含むモデルを解くことで発生することがある。
そんな時は出力を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を、要素の節点番号のくるみ方を$4!$=24通り試してみたが、(私のミスでなければ)すべて全く同じjacobianエラーとなった。

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 spooles solver
 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の所の解説にある、メッシュの切り方(分割数の指定 など)が非常に使える、ParaViewの解説もある
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_Aster_JPAubry_20131206.pdf

**マウスとキーボードを使った操作方法 [#mf12ea58]

,Ctrl+左クリックしながらドラッグ ,拡大&縮小
,Ctrl+真ん中クリックしながらドラッグ ,移されている領域を移動させる
,Ctrl+右クリック ,モデル回転

**salomeで.unvをつくりプログラムで.inpにする方法(解決) [#reb367d4]
***&link(Salomeメモ,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?Salome%a5%e1%a5%e2#i0)にあります。 [#w0f8c75c]
材料が2つの時もSalomeメモに書いてます。
**体積を出す [#x603c285]
Geometryで
**unv→inp用の載荷拘束テク [#deb31ce0]
-Salomeで、C3D4などでメッシュを切るとき有効
-きれいに載荷節点、拘束節点をいれたいとき有効
unv→inp変換する際は、Salome上で先に載荷拘束用の区切りをPartitionで入れておくと、
きれいに載荷拘束用の節点が入る。
unvにするf90では、
if文で
--(
区切りを入れた部分の座標-0.000001d0.and.区切りを入れた部分の座標+0.000001d0
--)
みたいに設定していくだけできれいに載荷拘束が入る。
手動で(CLOADで)面載荷したいの時も使える。

Partition入れまくったモデルを、Auto Hex(C3D8)とかでいっぺんに全部切ろうとすると、
当然エラーが出るので注意。
そういう意味ではC3D4は便利。(いっぺんに全部切っても大丈夫)

**unv→inp変換時の注意 prism [#q81eb19f]
salomeでメッシュを切った時にprismができるとcalculixでは対応していない要素となるので、計算しようとするとエラーが起きるか、要素が抜けたまま計算されてしまう。
prismができないようにメッシュを切らなければいけない。

prismは例えばSalomeでAutoで切ったときに、C3D4メッシュとC3D8メッシュの接続部などで生成されてしまう。

**メッシュをC3D8要素とかで切っていて、軸方向に多く切りたい [#ree29a8e]
-submeshを使う。
-Geometryで軸方向(多く切りたい方向)の辺すべて(例えば1個のただの箱なら軸方向に伸びている辺4つ)をグループでGroup_1とかに指定しておかなければいけない。

-まずBox_1(作りたいモデル)を作る。
-New Entity-Group-Createをクリックし、ラジオボタンで/(辺)を選択。
-Main Shapeには、Box_1を選択。するとBox_1がワイヤーフレーム表示になる。
-軸方向(z方向)に伸びている辺をクリックし、Add をクリック。
-これを軸方向(z方向)に伸びている辺すべて行う。一つの直方体なら、すべてAddすれば、
4つの辺番号が表示される。
-すべてAddしたら、Apply and Close 。この操作でできるGroupを、Group_1とする。
-Mesh画面に行き、&link(メッシュをC3D8で切って,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?Salome%a5%e1%a5%e2#i4)、例えばMesh_1ができる。
-その後、create submeshで、MeshにMesh_1を指定、GeometryにはGroup_1を指定する。
1D-AlgorithmでWire Distcretisation,HypothesisでNoneを指定し、その右にある歯車をクリックし、Nb.Segmentsをクリック。
-すると、新しいNb.Segmentsを作れるので、Number of Segmentsを軸方向に切りたい数だけ設定。

(デフォルトでは、名前が被らないようになっていると思うが、一応、この時はMesh_1を切るときに設定したNumber of Segmentsの設定のNameと、このとき設定するNumber of Segmentsの名前が被らないように注意する。)

-設定し終わったら、Applyする。するとツリーのMesh_1とあるところの絵に!がつく。
-Mesh_1を再びcomputeすると、submeshの設定が反映されて再びメッシュが切られる。

これを行うと、直方体とかなら普通にf90プログラムでnx=,ny=,nz=,と
設定するのと同じような感覚でメッシュが切れるし、軸方向分割数だけを多く稼げる。
unvからinpにしてccxで解いても、寸法・各種数値と分割数が全く同じなら全く同じ解になる(当然だが)。

**材料AとBからなるモデルを作りたいが、材料AとBの間をまたいでメッシュができてしまう [#r3e4a05b]
材料AとBの接合面を、Partitionで区切る。
Partitionでは接合面となるFaceがGeometryにある必要があるが、Geometryで、
ツリーの材料Aか材料Bを右クリックして、CreateGroup→ラジオボタンで□(面)を選び、
接合面をクリック→Addで接合面を作るのが効率が良い。

**非一体化モデルを作りたい [#ke26e2f0]
Geometryで、FuseのかわりにCompoundを使ってモデリングする。
Compoundは、Fuseとは違いMeshを切った時に一体化させた境界面で節点共有されなくなる。

ただし、Compoundした後に境界面をPartitionすると情報が上書きされて
その境目では節点が1ずつ入る(→つまり一体化する)ようだ。

-compound
選択した物体(複数可)をひとつのグループにまとめてくれる。
一体化してくれるわけではないので注意。
物体Aが物体Bにめり込んだ状態で2つをcompoundしても、
物体Aと物体Bが結合しないし、別々に物体Aと物体Bがあるという状態になる。

**メッシュを細かく切った時エラーみたいなのが出る [#z6874844]
elementsが500000(桁があってるか不明)とかより多いよ
みたいなエラーが出たら、Mesh画面でFile→Preferences→Mesh→Generalの
Mesh informationにあるAutomatic nodes compute limit を増やしたほうがいいかも。

ただしエラーが出ていても分割は終了しているっぽい(?)のでそのままでも問題はないかもしれない。
私はlimitを9999999にしている。

**グルーピングが遅い、もっと早くできないの? [#y84ecb45]
普通は、Create Meshでメッシュを切り、グルーピングを行うが、
Meshを切った後、一旦Saveし、Salomeを再起動させてからグルーピングを行うと、
普通にグルーピングするより早く終わる。
なんでだろうね。

**グルーピング後、unv→inp変換はうまく行くのに、inpの要素数が少なくなる。 [#y2930a8b]
丸みを持つモデルとかで発生する。
とりあえず、SalomeのMeshを右クリック→Mesh Informationでvolumeの数を確認し、
inpファイルで総要素数を見る。
inpファイルの方が要素数が少なかったら、グルーピングで要素数が抜け落ちてしまっている可能性がある。
丸みがある部分を持つモデルなどをC3D4(角がある要素)で作ると、少なからずその丸みは
メッシュを切る際に簡略化されてしまう。(角がある要素では丸みは"完全"には再現できないため。)
そのため、三角孔を持つ木橋などは、メッシュを切った後に、
鋼板部分をグルーピングする際はフィルターに三角孔がある鋼板を選択してしまうと、
すべての要素をグルーピングできなくなる。
こういった場合は孔を開ける前の鋼板をフィルターに設定してあげるとよい。

**C3D4とC3D6を混在させて、inpにして解くとjacobianエラーが出る [#wf7e1374]
-ccxエラー - jacobian問題 参照

**Submeshで部位ごとにメッシュを切る [#e1693575]

切り分けたい部位ごとに面などでpartitionを入れておき、
最終的なモデル(たとえばPartiton_1とする)
をExplode-Solid し、Partiton_1から小さいツリーで部位ごとに表示できるようにしておく。
Mesh画面でPartition_1を選択しCreate Meshするが、Algorithmなどは何も選択せずにApply and Closeし空のメッシュ(Mesh_1)を作る。

このMesh_1を選択し、Create sub-meshする。
GeometryでPartition_1の下の細かいツリーの中から切りたい部位を選択できるので、
それぞれ部位ごとにSubmeshを作ったら、Mesh_1をComputeすると、すべてのSubmeshの設定が反映されてメッシュが切れる。
Submeshごとに優先順位を変えることもできる。(Mesh_1を右クリックし、Change submesh priority)


-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~ [#te1e4500]
メモリが足りませんよエラー。
--(
Warning 
It seems that there is not enoght memory to show the mesh
so that the application may crash. Do you wish to continue visualization ?
--)
meshをcomputeしたとき、
Vision 2GのマシンでC3D4で(たぶん、およそ)200万要素くらいに分割しようとしたら発生。


**材料が2つのとき(を考えてる時に考えたこと) [#v86ece3e]
-例えば部材がbox_1の上にbox_2が乗っかっているモデルだとする。
-box_1とbox_2の境界面をpartitionで区切らなければいけないため、create groupを行う。
-ツリーのbox_1を右クリックし、Create Groupを選ぶ。
-ラジオボタンで、□(面)を選び、Main Shapeに面を指定する物体を選ぶ。(box_1とする)
-Geometryにbox_1だけを表示させて、box_2と接している面をクリックして、Addし、Apply and Close 。ツリーのbox_1の下に、Group_1ができる。
-これでpartitionのための設定が終わる。
-box_1とbox_2を一体化させてメッシュを切るので、ここからFuseする。
-Operations-Boolean-fuse で、object 1 にbox_1 、object 2にbox_2を選び、apply and Closeする。Fuse1ができる。
-Fuse1をpartitionで区切る。Operations-Partition を選び、ObjectsにFuse1を指定し、
Tool Objectsには、ツリーのbox_1の下できているGroup_1を選び、Apply and Close。
-Partition_1ができる。

-Mesh画面にいき、Partition_1を選択してメッシュを切る。(ここでMesh_1ができるとする)
-メッシュを切ったら、グルーピングを行う。Fuseしたモデルを切ったMesh_1を右クリックしてCreate Groupをクリック。
-MeshにMesh_1、Element Type にVolume、Nameは(かぶらなければ)好きなようにつける。まあ、Group_1mとする。(Groupの名前はGeometryにあるものとかぶっても問題はなさそうだが。)
-Group type ではGroup on filterにチェックを入れる。すると、Set Filterボタンが現れるので、クリックする。
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-今回は、Geometry のモデルでフィルターをかけるので、Criterion にBelong to Geom を
選び、Compare にはEqual to を選ぶ。
-そして、Threshold の中の空白をクリックしたら、ツリーの中からbox_1を選択。良ければApply&closeしてGroup設定画面もApply&close。
-これで、box_1の部分のメッシュがすべてGroup_1mに入れられたことになる。
-もう一度Create Groupする。今度はbox_2の部分のメッシュをグルーピングする。
-MeshにMesh_1、Element Type にVolume、NameはGroup_2mとかにする。
-Group type ではGroup on filterにチェックを入れる。すると、Set Filterボタンが現れるので、クリックする。
-Filterの設定画面が出るので、とりあえずAdd をクリック。
-Criterion にBelong to Geom を選び、Compare にはEqual to を選ぶ。
-そして、Threshold の中の空白をクリックしたら、ツリーの中からbox_2を選択。Apply&closeしてGroup設定画面もApply&close。
-これで、box_2の部分のメッシュがすべてGroup_2mに入れられたことになる。
-Mesh_1を右クリックして、unvファイルを保存する。

メッシュをC3D4で切り、2つにグループ分けしたものなら、ccxc3d4unv2.f90で解く。

partitionを使わない方法だと、C3D4とかだとbox_1とbox_2の間の部分で作られた要素が選択できない。
となると、メッシュの切り方を変えてBox_1とBox_2で別々に切って、
Compound Meshを作ってその後グループ分け、という作業が必要になるのではないか。

更に、Compoundするならその接合面で節点がくっつくほうがいいから、
Box_1の上にBox_2が乗っかってるモデルとするなら、FaceのBox_1と上端面とBox_2の下端面を2Dで全く同じように切って節点を作り、その節点を使うようにして(その面の節点は増やさないで)なんとかC3D4要素でBox_1と2を別々に切る(ここらへんはできるか不明)
という手順が必要かもしれない。調査中。
partition使ってないからこんなことになるのか?partition使ってもう一度やってみます。
-Compoundとか使うとunvの要素の入り方が変わってくる?。普通にc3d4unvのようにやると
最初が-864とかから始まる。

-partitionを使う方法
予め面を作っておき、partitionを使って区切ってからメッシュを作成、グルーピングしたら、普通に要素数の抜けなしでinpになった。
下の3角孔木箱桁橋でやったら要素の抜けが起こった。鋼板と木材が接している面を
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というコマンドを使った履歴を検索できる。
番号とその時使ったコマンドがずらっとでて、いちいち打ち直したくない場合は!番号 でそのコマンドが入る。後藤さんもaliasに入れてるらしい。
**inpファイルを見やすくする [#s9042b53]
&link(vimrc_ubuntu,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/.vimrc_ubuntu)を.vimrcの元々の中身を全部消して貼り付ける。

**メモリとかCPUの情報を見る [#zdf21144]
--(
more /proc/cpuinfo
more /proc/meminfo
--)

**外付けHDにバックアップをとる [#td205260]
--(
rsync -auv コピー元 コピー先
--)

**jobの最後にrm jobを入れとく [#racb5007]
jobの使い方は後藤資料の中でぐぐればあります。で、jobの最後に

--(
rm job
--)
を入れとくと、!./を使ったときjobが再実行されて、慌ててCtrl+cを押して強制終了するも、
解いた結果の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?%ba%ee%b6%c8%c9%f4%b2%b0


*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すると、グラフが表示されずに、.pngの中に保存される。
この2つの操作でグラフを.pngにして保存できる。

**1つのテキストファイルからいくつもグラフを書きたい!って人は、 [#xb60b973]
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 ,'hoge.txt' using 1:4 w l
--)
みたいにしてusingを使う。
using1:2 で、txtデータの1列目がxに入り、2列目がyに入るって意味なんですよ。


**線の太さをもっと太くしたい!って人は、 [#y04eee9e]
gnuplotでそのままグラフを作成すると、
スライド作成とかの時、グラフの線が細くて見にくいと言われたりする。
グラフの線を太くするには、linewidth を使う。 省略形でlwも使える。
lw 2みたいに、その後に半角スペースを空けて、整数を入れる。(数字が大きいと太い,デフォルトはlw 1の大きさ)
--(
pl 'hoge.txt' w l lw 2
--)

みたいにすると、太くなる。(だいたい2が見やすいと思う)

**グラフの線のタイトルを消してクールにしたい!って人は、 [#z308392a]
notitleを使う。
発表の時は、スライド作成ソフトやxfigとかで、グラフに矢印をつけて、FEM→ 初等梁→ ティモ→ みたいに書くほうがわかりやすい。(過去の卒論修論スライド参照)
そういう時には右上に自動生成されるグラフの線のタイトルがいらなくなるのでこれを使おう。

--(
pl 'hoge.txt' w l notitle
--)

*構造実験検証 [#x67cf198]
一人で(龍さんに手伝ってもらったりして)こっそりやったデータ。
-4/24午後7時 二回とも成功(変な測定値は出なかった)

--2回目データ:&link(140424a_2s.csv,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/zikken/140424a_2s.csv)
--1回目データ:&link(140424a_1s.csv,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/zikken/140424a_1s.csv)

*3D [#d2eb7ba5]
**引張試験 試験体 [#h211cd60]
引張試験のための供試体を作成する。
ものづくりセンターの3Dプリンタでは造形情報をstlファイルとして読み込むので、モデリングツールでモデリングした後にstlファイルに変換しなければいけない。
SalomeのGeometryでモデルを選択→export→STLAⅡfilesで吐き出す。


**引張試験 [#f5580cac]
***メモ [#a44cbbb9]
-ロガーにつながってるPcのPCD-30Aの集録データ数 の欄で何秒測定するか決められる。
-doループの終わりをデータ終わりの秒数にするだけでヘッダーがついたままで初期不整が入れられる&link(syokignu2.f,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/prog/syokignu2.f)

*回転らせん [#eeb2f178]

遊びで適当に節点だけ打ってみたけど、
中間節点を作るのが難しそう

http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/scpng/raseninp.png



*プログラム類↓ [#wa71b442]
&link(卒論日誌,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%c0%c6%c6%a3%b5%b1%a4%ce%c2%b4%cf%c0%c6%fc%bb%ef)とかにも色々転がってます。

*プレストレス木箱桁橋の角材の間にゴムを挟めるモデル [#y896aff0]
接触解析の代わりになるかも という事で滝田さんのプログラムを改変して作成しました。

ゴムの寸法や、分割数、ヤング率(ELgom)とかは各自で変えてください。

節点座標の出力桁を多くして、
ゴムを小さくした時に節点座標が出力桁以下になり、
四捨五入されてしまいメッシュの大きさが均一にならないことがある不具合を修正したgom2.f90

&link(gom2.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/model/gom2.f90)

*B32用のunv→inp変換 [#v1d92ed0]
2次元のトラスにしか対応してません。
3次元には非対応。

主な手順

-まずSalomeでyz平面上にトラスを作る
-点を打ってってそれをLineにして
-トラスができたら、全部のLineを選択して1つのWireをつくる
-その後,作ったWireを,
-1Dでメッシュを切るが、number of segmentsで分割数を3以上の偶数に指定して切る.
-Meshを.unvに吐き出して&link(b32unv.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/saito/2014/prog/b32unv.f90)の中の必要事項を書き換えてinp化

*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
,月,,構造力学特論,,パソコンゼミ
,火,,都市システム計画特論,土質工学特論,
,水,,材料設計特論,,
,木,,,交通システム計画特論(英語ゼミ),構造実験
,金,,,,


トップ   編集 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS