大黒屋の修論日誌 -

atコマンド

http://discypus.jp/wiki/?Linux%2Fat%A5%B3%A5%DE%A5%F3%A5%C9

構造工学作業用

09/1/13後藤メモ

深さ25%, 荷重140の側面のせん断ひずみ

sが側面、tが対称面。横軸は節点番号。データは、
側面:eyz25p140s.gnu
対称面:eyz25p140s.gnu
$ε_{y z}$ 全データ:eyz25p140.d
側面の中立軸だけ抜き出すプログラムeyz.f

次に梁側面の中立軸部9箇所でロゼットゲージにより測定した代表的なせん断ひずみ分布のグラフを，2章の手法によるFEM解析の結果とともに図-10に示す．ロゼットゲージによる測定値は，FEM解析の結果と比較的一致している．ちなみにFEM解析のひずみ分布は，初等梁の範疇で求まるせん断力図分布とも対応している。

FEM(対称面,x=0,y=0)
FEM(側面, x=40mm,y=0)

09/1/9後藤メモ

とりあえずまとめ。s4試験体(No.6)について、

2p/3載荷点(1/4点)と1/8点の側面

gnuplot用のプロットデータは、1/4点が、

s4s14.gnu

gnuplot用のプロットデータは、1/8点が、

s4s14.gnu

sgmzz.f入力用の

s4.d

gnuplotで

setxrange[-2100:2100]
setyrange[100:-100]
pl 's4s14.gnu' w l, 's4s18.gnu' w l

みたいに。

単軸ゲージにより測定した $\frac{2 P}{3}$ 載荷位置での代表的な軸ひずみ分布を， 3章?の手法によるFEM解析の結果とともに図--9に示す．いずれの試験体でも，縁部付近のひずみは合成断面としての三角形分布に近いが，縁部より少し中立軸側に離れた部分では合成断面としての三角形分布よりも小さくなっている。この傾向はFEM解でも認められる。ただしFEM解の場合，載荷部から離れた位置(1/8点)での応力分布は中立軸付近でひずみが特に小さくはならず三角形分布に近い。ということは，FEM解析における載荷節点付近の応力集中などによる影響が，試験体における不整(例えば，鋼板と集成材とが十分に一体化していないなど)による影響とたまたま近い結果を与えたという可能性もある．次にロゼットゲージにより測定した代表的なせん断ひずみ分布のグラフを図--10に示す．梁側面の中立軸部9箇所で測定した....

みたいな感じか。これでいくとすれば、図-9に、No.6試験体のピーク荷重での、 $\frac{2 P}{3}$ 載荷点(1/4点)、1/8点の側面の応力分布を併記しますか。

09/1/8後藤メモ

10*40*202の逆対称4点曲げで、左載荷点側面の軸方向応力分布を桁高に対して見るデータを取り出すには、viのビジュアルモードとかを駆使して、 ezz25p140.d のようにfrdファイルのひずみ出力部分を整形して、 sokumen.fを使って、./sokumen<ezz25p140.d>ezz25p140.gnuみたいに実行して、 gnuplotで描画。任意の要素数で任意の桁高線の軸応力分布を見るには sgmzz.f を使う。

載荷部の断面内節点に荷重を分散させて断面載荷した場合
2p/3載荷点(1/4点)の側面では、

割と測定値に近いひずみ分布なのだが、

2p/3載荷点(1/4点)の対称面(鋼板挿入部)では、

こんな感じ。なんかおかしい気がするんで、

載荷点から離れた(1/8点)の側面では、

比較的三角形分布に近い。そして

載荷点から離れた(1/8点)の対称面(鋼板挿入部)ですら、

比較的三角形分布に近い。

ちなみに、載荷を梁の上縁に線載荷した場合は、かなり応力が乱れる。
2p/3載荷点(1/4点)の側面では、

2p/3載荷点(1/4点)の対称面(鋼板挿入部)では、

載荷点から離れた(1/8点)の側面では、

載荷点から離れた(1/8点)の対称面(鋼板挿入部)では、

09/01/05 大黒屋
12/28 大黒屋

スギ

まつ

構造工学の修正意見投稿ページが、書き込み中に頻繁に落ちるので、修正意見は、まず、学内限定wiki にメモを取ってから投稿した方が安全。

dvipdfmxのエラー

pstexがあるやつを変換しようとするとエラー。いつもってわけじゃないけど。

[1Error: /undefinedresource in --findresource--
Operand stack:
   --nostringval--   Ryumin-Light-EUC-H   Font   (HeiseiMin?-W3-EUC-H)   12   Ryumin-Light   EUC-H
  EUC-H   --dict:0/10(G)--   FDepVector   --nostringval--   --dict:11/13(ro)(G)--   HeiseiMin?-W3-EUC-H   
Font   HeiseiMin?-W3-EUC-H   (HeiseiMin?-W3-EUC-H)   12   HeiseiMin?-W3   EUC-H   EUC-H   
--dict:0/10(G)--   HeiseiMin?-W3   false   HeiseiMin?-W3   HeiseiMin?-W3   CIDFont   --nostringval--   
CMap   CIDFont   CMap   CIDFont
Execution stack:
   %interp_exit   .runexec2   --nostringval--   --nostringval--   --nostringval--   2   %stopped_push   
--nostringval--   --nostringval--   --nostringval--   false   1   %stopped_push   1   3 %oparray_pop   
1   3   %oparray_pop   --nostringval--   1   3   %oparray_pop   1   3   %oparray_pop   .runexec2   
--nostringval--   --nostringval--   --nostringval--   2   %stopped_push   --nostringval--   
--nostringval--   13   4   %oparray_pop   14   4   %oparray_pop   --nostringval--   --nostringval--   
27   5   %oparray_pop   --nostringval--   30   6   %oparray_pop   --nostringval--   --nostringval--   
--nostringval--   --nostringval--
Dictionary stack:
   --dict:1072/1123(ro)(G)--   --dict:0/20(G)--   --dict:114/200(L)--   --dict:6/16(L)--   
--dict:17/17(ro)(G)--   --dict:0/10(L)--   --dict:16/24(ro)(G)--
Current allocation mode is local
Last OS error: 2
Current file position is 5866
GNU Ghostscript 7.07: Unrecoverable error, exit code 1

** WARNING ** Filtering file via command -->gs -q -dNOPAUSE -dBATCH -sPAPERSIZE=a0 -sDEVICE=pdfwrite 
-dCompatibilityLevel=1.3 -dAutoFilterGrayImages=false -dGrayImageFilter=/FlateEncode? 
-dAutoFilterColorImages=false -dColorImageFilter=/FlateEncode? -dUseFlateCompression=true 
-sOutputFile=/tmp/dvipdfmx.YqddTW -dNOKANJI -c '.setpdfwrite << /NeverEmbed? [/Courier /Courier-Bold 
/Courier-Oblique /Courier-BoldOblique? /Helvetica /Helvetica-Bold /Helvetica-Oblique /Helvetica-
BoldOblique? /Times-Roman /Times-Bold /Times-Italic /Times-BoldItalic? /Symbol /ZapfDingbats? /Ryumin-Light 
/GothicBBB-Medium] >> setdistillerparams' -f ./us.pstex -c quit<-- failed.

** WARNING ** Image format conversion for "./us.pstex" failed...
** ERROR ** pdf_ref_obj(): passed invalid object.

Output file removed.

12/24

http://tutimura.ath.cx/ptetex/?%C6%B0%BA%EE%CA%F3%B9%F0%2F41

CalculiX

逆対称の解析

曲げ応力とせん断応力

荷重P=140kN
E木は各供試体の平均値を使った。

カラマツ

スギ

後藤ちゃちゃ08/9/26

これが本当なら、実験の通りだよね。というか、この図の方が説得力があるよね。

08/11/06

応力をみてみたが、理論値と大幅にちがう。要素数を変えてもう一度やってみる。広内にも同じ計算をやってもらう。

08/11/10

要素数を変えてみたけどどうもあわない。変位はあってるっぽいけど。

08/11/11

スギ材の応力を表にしてみた。

深さ	軸方向(2P/3載荷部)	せん断応力(中央部)MPa
0	67.6	4.32
5	48.4	3.87
10	41.4	3.78
15	37.6	3.79
20	35.3	3.88
25	33.5	4.00
30	31.9	4.17
35	30.5	4.47
40	29.0	5.20
45	27.3	7.85
50	23.6	52.8

08/11/13

スギせん断支間中央

マツせん断支間中央

軸方向応力左載荷点

軸方向応力の断面図(スギ)左載荷点
せん断応力の断面図(スギ)支間中央

08/11/17

せん断は鋼板の間で応力集中している。

軸方向は計算とあわない。

軸方向応力(右支点部)

支点部も拘束の影響がありそう?。左載荷と支間中央の間、要は左支点から0.75m右のとこで応力をみてみる。

軸方向応力(左支点から0.75m右)

08/11/19

応力を強度で無次元化したやつ。強度はこっち(pdf)が最新?

スギ材

カラマツ材

断面のグラフは明日だな。

08/11/21

スギ材の支間中央のせん断応力のx方向の分布。左側が鋼板挿入部分。

08/11/26

実験のピーク荷重からだした応力を平均したの使ってで無次元化したやつ。

スギ材

カラマツ材

08/11/27

千田さんの論文を参考に、上下に1要素いれて、線載荷のやつをやってみる。それで応力集中は軽減できるらしい。

08/12/02

ngyaku1p.fは、 y方向に1要素増えるから、要素数オーバー?で ny=40だと走らないから、ny=30(32)で。

比較のために、 gyaku1p.fで ny=30でもやってみる。

08/12/05

支間中央のせん断応力。面積載荷も線載荷もほぼ同じ。線が重なっている。

スギ材

カラマツ材

支点拘束と載荷を中立軸にもってくる。cgyaku.f完成?。
cgyaku.f

マツ25%の軸方向の応力比較

上のやつに、中立軸拘束、面積載荷を追加。重なって分かりにくいけど

08/12/08

稲荷くんからもらった

sankaku.inp

08/12/11

鋼材を引っ張って降伏について調べている。

CalculiXでひずみの出力もできることが判明。

ここ

08/12/12

マツ25%の左載荷点のひずみ分布と応力分布

下拘束、上鉄板載荷

下拘束、上載荷

中立拘束、上載荷

中立拘束、中立載荷

塑性判定の降伏応力は、実験値から読み取って、

0%のとき、鋼材に集成材の材料定数をいれてたけど、

08/12/15

塑性の計算が収束していない。

深さ25%で

スギ材

b=80mm
h=200mm
L=0.5m
ヤング率7.49GPa
せん断弾性係数442MPa
P=140kN
nx=10
ny=40
nz=202

載荷位置の変位

深さ	FEM解(m)	ティモシェンコ解(m)
0	0.014663	0.011769
5	0.010722	0.0094970
10	0.0090355	0.0085735
15	0.0080776	0.0080935
20	0.0073199	0.0078158
25	0.0066528	0.0076483
30	0.0060291	0.0075479
35	0.0053904	0.0074910
40	0.0047474	0.0074629
45	0.0040351	0.0074527
50	0.0025827	0.0074512

マツ材

b=80mm
h=200mm
L=0.5m
ヤング率16.6GPa
せん断弾性係数538MPa
P=140kN
nx=10
ny=40
nz=202

載荷位置の変位

深さ	FEM解	ティモシェンコ解
0	0.0087715	0.0072654
5	0.0073216	0.0067142
10	0.0065374	0.0064028
15	0.0060294	0.0062116
20	0.0055612	0.0060896
25	0.0051041	0.0060115
30	0.0046496	0.0059631
35	0.0041681	0.0059350
40	0.0036776	0.0059209
45	0.0031344	0.0059158
50	0.0020342	0.0059151

大黒メモ08/10/09

FEM解が、ティモシェンコ理論値の2倍程度ででてくる。荷重の制御の問題か、計算を間違っているかな。ちなみに、

スギ材
b=80mm
h=200mm
L=0.5m
ヤング率7.53GPa
せん断弾性係数427MPa
P=140kN

0%
初等梁理論=6.5E-3(m)
ティモ理論=1.2E-2(m)
FEM解=2.37E-2(m)
5%
初等梁理論=4.2E-3(m)
ティモ理論=9.7E-3(m)
FEM解=2.14E-2(m)

逆対称4点曲げのたわみの式はここ使ったプログラムはこれを改造。

逆対称4点曲げの鋼板挿入のたわみの計算プログラム

soutai2.f

大黒メモ08/10/24

どうやら、要素分割のせいみたい。ちゃんとした計算できたっぽい。非線形とも数%の違い。来週から、応力もみてみる。

非線形を入れたときと入れないときの比較

b=80mm
h=200mm
L=2m
木材のヤング率:7.49GPa
木材のせん断弾性係数:442MPa
鋼材のヤング率:206GPa
鋼材のせん断弾性係数:79.2GPa
拘束条件:片持ち梁
挿入深さ:25%

初等梁理論

v = \frac{P L^{3}}{3 E I}

ティモシェンコ梁理論

v = \frac{P L^{3}}{3 E I} + \frac{P L}{k G A}

FEM解と梁理論を比較。曲げ剛性は $E_{木} I_{木} + E_{鋼} I_{鋼}$ せん断剛性は $k G_{木} A$ を使用。

荷重kN	初等(m)	ティモ(m)	FEM非線型(m)	FEM線形(m)
1	0.00240619224	0.00274555874	0.00275490	0.00275490
5	0.0120309601	0.0137277925	0.0137743	0.0137745
10	0.0240619201	0.0274555851	0.027547	0.0275490
20	0.0481238402	0.0549111702	0.055088	0.0550979
30	0.0721857622	0.0823667571	0.082614	0.082647
40	0.0962476805	0.10982234	0.11012	0.11020
50	0.120309606	0.137277931	0.13759	0.13774
60	0.144371524	0.164733514	0.16503	0.16529
70	0.168433443	0.192189097	0.19244	0.19284
80	0.192495361	0.219644681	0.21980	0.22039
90	0.216557294	0.247100279	0.2.4711	0.24794
100	0.240619212	0.274555862	0.27439	0.27549
110	0.264681131	0.30201146	0.30162	0.30304
120	0.288743049	0.329467028	0.32883	0.33059
130	0.312804937	0.356922597	0.35602	0.35814
140	0.336866885	0.384378195	0.38324	0.38569
150	0.360928804	0.411833793	0.41064	0.41323

非線形でP=150kN以降は計算が収束しないようで、

マツ集成材実験

ヤング率

挿入深さ	供試体	ヤング率(Pa)	せん断弾性係数(Pa)	EwIw?(kN^2)	EwIw?+EsIs?(kN^2)	kGw(GPa)
25%	m1	1.75841792E+10	562772608.
25%	m2	1.76567071E+10	520052352.
15%	m3	1.54035507E+10	526195808.
35%	m4	1.53642168E+10	666155968.
0%	m5	1.69810125E+10	507468256.
0%	m6	1.6761002E+10	447247200.

応力

試験体	鋼板深さ	Ew(G)	合成EI	ピーク	σ(MPa)	τ(MPa)
m6	0	16.8	906k	98k	30.3	3.03
m5	0	17.0	894k	114k	36.1	3.61
m3	15	15.4	1228k	113k	23.6	3.48
m2	25	17.7	1362k	163k	35.3	5.64
m1	25	17.6	1518k	154k	29.8	4.77
m4	35	15.4	1766k	177k	25.7	4.75
	平均	16.6		30.1	4.21

実験写真

http://www.str.ce.akita-u.ac.jp/~daikoku/jikken/

スギ集成材実験

ヤング率

挿入深さ	供試体	ヤング率(Pa)	せん断弾性係数(Pa)	EwIw?(kN^2)	EwIw?+EsIs?(kN^2)	kGw(GPa)
15%	s1	7.0131637E+9	395843641	374	897
15%	s2	7.34500209E+9	365065651	392	914
25%	s3	6.45406961E+9	540401788	344	1042
25%	s4	7.68432781E+9	446546520	410	1104
35%	s5	8.12242113E+9	400955610	433	1203
35%	s6	8.25382151E+9	377674406	440	1210
0%	s7	7.85441064E+9	450457214	419	419
0%(×)	s8	7.2322791E+9	562026668	386	386

応力

試験体	鋼板深さ	Ew	合成EI	ピーク	σ(MPa)	τ(MPa)
s7	0	7.85	7.85k	115k	36.0	3.60
s1	15	7.01	868k	137k	18.4	3.85
s2	15	7.35	786k	141k	22.0	4.47
s3	25	6.45	976k	170k	18.7	5.13
s4	25	7.68	977k	124k	16.2	3.98
s5	35	8.12	1154k	150k	17.6	4.68
s6	35	8.25	1297k	135k	14.3	3.77
平均		7.53			25.4(0%と15%)	4.39(25%と35%)

写真

http://www.str.ce.akita-u.ac.jp/~daikoku/picture/

adobereader

adobeからpdfファイルを開くと

(acroread:2948): GLib-GObject-WARNING **:

 plugin '/usr/lib/gtk-2.0/immodules/im-scim.so'

 failed to register type 'GtkIMContextSCIM'


(acroread:2948): Gtk-WARNING **: Loading IM context type 'scim' failed

ていうエラーがでる。印刷とかも問題なさそうだけど。

大黒link

有限要素法のフリーのソフトウエア CalculiX を試す http://freecaetester.blog62.fc2.com/

Gmsh

google2001 http://www.google.com/search2001.html
木質構造研究所 http://e-sunagawa.co.jp/tei01/index.html
専門用語(英語)http://ubweb.info/cetet/index.html
単位変換http://hp.vector.co.jp/authors/VA018451/javascript/jdoryoko.htm
いろんな定数や公式 http://ebw.eng-book.com/heishin/vfs/
便利帳 http://www.fujimfg.co.jp/benri/benri.htm
工学便利帳 http://hp.vector.co.jp/authors/VA016605/HANDBOOK/index.htm
tex http://www.biwako.shiga-u.ac.jp/sensei/kumazawa/tex.html
tex2 http://grace.ceser.hyogo-u.ac.jp/tex/index.html
tex3 http://hooktail.org/computer/index.php?TeX
tex4 http://www-an.acs.i.kyoto-u.ac.jp/~fujiwara/tex/index.htm
解凍 http://uguisu.skr.jp/Windows/tar.html

Xfig(小林:見たら消してください)

どうやらundo回数変更できないみたいです。違うcadとかのソフトと勘違いしていたのかも。すみません;;

undo

このソフト

大黒08/01/09

そっか。残念。わざわざありがとう。

後藤メモ07/6/12

千田さんのプログラムに河原さんが注釈を入れたものは

論文関係

構造工学委員会

構造工学論文集情報

〆切は9月18日。

木橋シンポ

08木橋シンポ

せん断弾性の変動係数を求める。

	平均	標準偏差	変動係数
kGw	0.355	0.0452769257	0.127540636
合成kG	0.5	0.0993730346	0.198746069

発表概要

(1)プログラム：以下のホームページをご覧下さい。
   http://tbl.tec.fukuoka-u.ac.jp/jsce-mokkyo-sympo/info.html
(2)発表時間：発表10分(第1鈴8分)、質疑5分、合計15分
　　発表方法の詳細は、上記のホームページ「発表に関するご案内」をご覧下さい。
(3)表彰：優秀な論文・研究発表に対し、論文賞、発表賞を授与いたします。
(4)参加申込：発表者も参加申し込みが必要です。申込方法は上記ホームページを
　　参照して下さい。(申込締切=8月10日)
(5)懇親会：シンポジウム終了後の17:30より懇親会を開催いたします。
  是非ご参加ください。
　参加費：5000円（当日受付）、会場：主婦会館プラザエフ。詳細は上記のホームページをご覧下さい。）

論文、要旨のスタイル

http://tbl.tec.fukuoka-u.ac.jp/jsce-mokkyo-sympo/info.html

〆切が一週間延期になっている。

最終原稿〆切

7月16日まで。査読内容を反映させて提出。オフセット原稿は郵送なので、7月9日の週には遅くとも出すようにする。

後藤ちゃちゃ07/4/11

投稿前の論文要旨を一般公開してしまうのはやや問題があるかも知れないので、念のため、以下に http://g2.ce.akita-u.ac.jp/cgi-bin/g2wiki/wiki.cgi?%c2%e7%b9%f5%cd%d1 (学内のみアクセス可)

授業＆TA計画

大学院シラバス

http://syllabus.eng.akita-u.ac.jp/syllabus/syllabus_graduate.php

2008後期分

	時間	月	火	水	木	金
1.2時限	8:50-10:20
3.4時限	10:30-12:00
5.6時限	12:50-14:20			後藤班ゼミ
7.8時限	14:30-16:00			情報処理の技法	構造力学特論

2008前期分

	時間	月	火	水	木	金
1.2時限	8:50-10:20
3.4時限	10:30-12:00
5.6時限	12:50-14:20			パソコンゼミ
7.8時限	14:30-16:00				構造実験(TA)

2007後期分

	時間	月	火	水	木	金
1.2時限	8:50-10:20		構造設計学特論
3.4時限	10:30-12:00			環境材料学
5.6時限	12:50-14:20	物理学実験(TA)	都市システム計画特論		特別講義Ⅰ	木構造工学
7.8時限	14:30-16:00	〃		情報処理の技法(TA)

情報処理センターの登録について

登録方法(07/04/06)

登録用のパソコンに、学生証を入れる。
IDが読み込まれるので、生年月日を入力する。
次にパスワードを設定する画面になるので6文字以上14文字未満?で入力(詳細は画面参照)
登録完了。webメールも情報処理センターの使用申請も同時に行われ、同一ID、パスワードでok.
登録後すぐに使用可。

基礎物理学実験

後期は11月開始予定?

シラバス

情報処理の技法

http://www.str.ce.akita-u.ac.jp/~gotou/zyouhou/

VBA

構造実験(アクリルの引張試験)

http://www.str.ce.akita-u.ac.jp/~gotou/zikken/

大黒屋の修論日誌 -

目次

atコマンド

構造工学作業用

09/1/13後藤メモ

09/1/9後藤メモ

09/1/8後藤メモ

dvipdfmxのエラー

CalculiX

逆対称の解析

曲げ応力とせん断応力

後藤ちゃちゃ08/9/26

08/11/06

08/11/10

08/11/11

08/11/13

08/11/17

08/11/19

08/11/21

08/11/26

08/11/27

08/12/02

08/12/05

マツ25%の軸方向の応力比較

08/12/08

08/12/11

08/12/12

下拘束、上鉄板載荷

下拘束、上載荷

中立拘束、上載荷

中立拘束、中立載荷

08/12/15

スギ材

載荷位置の変位

マツ材

載荷位置の変位

大黒メモ08/10/09

大黒メモ08/10/24

非線形を入れたときと入れないときの比較

マツ集成材実験

ヤング率

応力

実験写真

スギ集成材実験

ヤング率

応力

写真

adobereader

大黒link

Xfig(小林:見たら消してください)

大黒08/01/09

後藤メモ07/6/12

論文関係

構造工学委員会

木橋シンポ

08木橋シンポ

発表概要

論文、要旨のスタイル

最終原稿〆切

後藤ちゃちゃ07/4/11

授業＆TA計画

大学院シラバス

2008後期分

2008前期分

2007後期分

情報処理センターの登録について

登録方法(07/04/06)

基礎物理学実験

情報処理の技法

VBA

構造実験(アクリルの引張試験)