江村の修論日誌 -

目次

• みんな
• 補足
• ダイヤカット円筒が内外圧に強いと言われるわけ
• 雪平とか半円筒とか泥船渡河
• ダイヤカット半円筒
• 自由端にヤング率$10^6$倍のビーム要素を付けたケース
• 軸方向:8
• 軸方向:12
• 軸方向:16
• 軸方向:20
• ダイヤカット半円筒(全方向変位・回転拘束)
• (有)ビーム要素ダイヤカット半円筒の曲げ剛性
• (有)ビーム要素ダイヤカット半円筒の曲げ剛性(回転も拘束)
• ダイヤカット半円筒座屈解析(ビーム要素つき)
• 半円筒
• 自由端にヤング率$10^6$倍のビーム要素を付けたケース
• 雪平半円筒
• 自由端にヤング率$10^6$倍のビーム要素を付けたケース
• 曲げ試験(FullCure720)
• 試験方法
• 横方向
• スパン20cm(たわみ測定点は左支点より9cm)
• スパン18cm(たわみ測定点は左支点より8cm)
• スパン16cm(たわみ測定点は左支点より7cm)
• スパン14cm(たわみ測定点は左支点より6cm)
• スパン12cm(たわみ測定点は左支点より5cm)
• スパン10cm(たわみ測定点は左支点より4cm)
• スパン8cm(たわみ測定点は左支点より3cm)
• スパン6cm(たわみ測定点は左支点より2cm)
• スパン4cm(たわみ測定点は左支点より1cm)
• グラフ
• 引張試験(FullCure720)
• 試験方法
• 横方向
• 造形直後(10月9日15:30造形)
• 養生2日(10月9日15:30造形)
• 養生7日(10月9日15:30造形)
• 0014 1/2
• ひずみ比較
• ポアソン比比較
• 孔削り油差し試験体
• 予備引張試験(FullCure720)
• viコマンド
• 論文
• 英語論文のタイトル(発声練習-英語論文のタイトルで大文字にするのはどこか-より)
• 雑感
• なぜ広まった？ 「『訊く』が正しい」という迷信より
• 「大学の広報誌にカバーガールを載せるのはいかがなものか」を読んで
• anet「【入試課】 秋田大学「入試・イベント情報」 facebookについて」
• アルミ缶の寸法
• jacobianエラーについて
• g77をインストールしてccxs6.f90をg77でコンパイル(拡張子変更後)してinpファイルをつくった。
• Gmailのhtmlファイル非添付
• エラー
• メッシュの比較
• 静的解析
• m=8,n=8
• 座屈解析
• m=8,n=8
• S6に変換
• ダイヤカットをSALOMEでモデリング
• Elmer
• tex.vim??
• viで複数行に文字を挿入
• 座屈解析の覚書
• 無次元化
• 曲げ剛性
• 座屈荷重
• 半ドーナツ型の断面二次モーメントを任意の軸周りで求められるプログラム
• niji9.f
• 剛性を求めるプログラム
• gousei.f
• 半ドーナツ型を縦分割して断面二次モーメントを求めるプログラム
• niji7.f
• 矩形断面を縦横分割して断面二次モーメントを求めるプログラム
• niji3.f
• 矩形を縦分割して断面二次モーメントを求めるプログラム
• niji1.f
• 2011年度前期
• 2011年度後期

高野豆腐

みんな

• 後藤さんが編集してくれた卒論・修論発表のisoファイルいいぞ!!DVDに焼くべし。
実家に送ればきっと喜んでくれるはず。
• 裏謝辞は書いたほうが面白いぞ!!
• 改めて修論発表を見直すとめちゃくちゃおもしろい。
限られた時間の中でいかにわかりやすく説明するかってむずかしいなあ。 早口になっちゃうしなあ。いやあむずかしい。

補足

• CalculiXによるダイヤカット圧縮時の変位(ヤング率102MPa)と実験値
daiya102mpa.txt

ダイヤカット円筒が内外圧に強いと言われるわけ

雪平とか半円筒とか泥船渡河

• サイズ

ダイヤカット半円筒

gou.txt

自由端にヤング率$10^6$倍のビーム要素を付けたケース

軸方向:8

• 周方向:8,軸方向:8

計算結果(載荷部分中央の節点番号は3)から $EI=24.18$N$\times$m$^{2}$

• ccxs6daiyahan.fをで1要素を更に4分割。

• 周方向:8,軸方向:8

計算結果(載荷部分中央の節点番号は3)から$EI=18.47$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan88h.frd

• 周方向:16,軸方向:8
12月3日
• ccxs6daiyahan.fで1要素を更に4分割。

計算結果(載荷部分中央の節点番号は5)から$EI=47.09$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan168h.frd

• 周方向:32,軸方向:8
12月3日
• ccxs6daiyahan.fで1要素を更に4分割。

計算結果(載荷部分中央の節点番号は9)から$EI=59.01$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan328h.frd

• 周方向:64,軸方向:8
12月3日
• ccxs6daiyahan.fで1要素を更に4分割。

計算結果(載荷部分中央の節点番号は17)から$EI=59.62$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan648h.frd

軸方向:12

• 周方向:8,軸方向:12
12月3日

計算結果(載荷部分中央の節点番号は3)から$EI=14.01$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan812h.frd

• 周方向:16,軸方向:12
12月4日

計算結果(載荷部分中央の節点番号は5)から$EI=38.27$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan1612h.frd

• 周方向:32,軸方向:12
12月4日

計算結果(載荷部分中央の節点番号は9)から$EI=58.47$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan3212h.frd

• 周方向:64,軸方向:12
12月4日

計算結果(載荷部分中央の節点番号は17)から$EI=59.37$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan6412h.frd

軸方向:16

• 周方向:8,軸方向:16
12月3日

計算結果(載荷部分中央の節点番号は3)から$EI=10.66$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan816h.frd

• 周方向:16,軸方向:16
12月4日

計算結果(載荷部分中央の節点番号は5)から$EI=31.85$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan1616h.frd

• 周方向:32,軸方向:16
12月4日

計算結果(載荷部分中央の節点番号は9)から$EI=57.69$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan3216h.frd

• 周方向:64, 軸方向:16
12月4日

計算結果(載荷部分中央の節点番号は17)から$EI=59.53$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan6416h.frd

軸方向:20

• 周方向:8, 軸方向:20
12月4日

計算結果(載荷部分中央の節点番号は3)から$EI=8.422$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan820h.frd

• 周方向:16, 軸方向:20
12月4日

計算結果(載荷部分中央の節点番号は5)から$EI=27.46$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan1620h.frd

• 周方向:32,軸方向:20
12月4日

計算結果(載荷部分中央の節点番号は9)から$EI=56.71$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan3220h.frd

• 周方向:64,軸方向:20
12月4日

計算結果(載荷部分中央の節点番号は17)から$EI=47.31$N$\times$m$^{2}$

• 変形図(倍率200)

ccxs6daiyahan6420h.frd

ダイヤカット半円筒(全方向変位・回転拘束)

• 周方向パターン数:8~64、軸方向パターン数:8
• ccxs6daiyahan88r.dat

• ccxs6daiyahan168r.dat

• ccxs6daiyahan328r.dat

• ccxs6daiyahan648r.dat

• 周方向パターン数:8~64、軸方向パターン数:12
• ccxs6daiyahan812r.dat

• ccxs6daiyahan1612r.dat

• ccxs6daiyahan3212r.dat

• ccxs6daiyahan3212r.dat

• ccxs6daiyahan6412r.dat

• 周方向パターン数:8~64、軸方向パターン数:16
• ccxs6daiyahan816r.dat

• ccxs6daiyahan1616r.dat

• ccxs6daiyahan3216r.dat

• ccxs6daiyahan6416r.dat

• 周方向パターン数:8~64、軸方向パターン数:20
• ccxs6daiyahan820r.dat

• ccxs6daiyahan1620r.dat

• ccxs6daiyahan3220r.dat

• ccxs6daiyahan6420r.dat

(有)ビーム要素ダイヤカット半円筒の曲げ剛性

ccxs6daiyahangraph.txt

• 周方向パターン数に対して敏感
• 周方向パターン数が多いほど、軸方向パターン数が少ないほど高い値になる傾向
• 軸方向パターン数が多いほど表面積が増える分、硬くなりそうだが

(有)ビーム要素ダイヤカット半円筒の曲げ剛性(回転も拘束)

ccxs6daiyahanRgraph.txt

• 拘束の差は周方向パターン数が増えるほど小さくなり、軸方向パターン数が増えるほど小さくなる。

ダイヤカット半円筒座屈解析(ビーム要素つき)

ccxs6daiyahanRgraph.txt

倍率0.0005

ccxs6daiyahan88z.frd

ccxs6daiyahan168z.frd

ccxs6daiyahan328z.frd

ccxs6daiyahan648z.frd

ccxs6daiyahan812z.frd

ccxs6daiyahan1612z.frd

ccxs6daiyahan3212z.frd

ccxs6daiyahan6412z.frd

ccxs6daiyahan816z.frd

ccxs6daiyahan1616z.frd

ccxs6daiyahan3216z.frd

ccxs6daiyahan6416z.frd

ccxs6daiyahan820z.frd

ccxs6daiyahan1620z.frd

ccxs6daiyahan3220z.frd

ccxs6daiyahan6420z.frd

• 軸方向:40
• 周方向:8,軸方向:40
12月4日

計算結果(載荷部分中央の節点番号は3)から$EI=3.986$N$\times$m$^{2}$

• 軸方向:100
• 周方向:8,軸方向:100
12月4日

計算結果(載荷部分中央の節点番号は3)から$EI=1.879$N$\times$m$^{2}$

• 周方向:36,軸方向:8

計算結果(載荷部分中央の節点番号は10)から $EI=59.59$N$\times$m$^{2}$

• ccxs6daiyahan.fで1要素を更に4分割。

計算結果(載荷部分中央の節点番号は10)から$EI=59.23$N$\times$m$^{2}$

• 周方向:100,軸方向:2

計算結果(載荷部分中央の節点番号は26)から $EI=60.48$N$\times$m$^{2}$

• 周方向:8,軸方向:40

計算結果(載荷部分中央の節点番号は3)から $EI=8.393$N$\times$m$^{2}$

• 周方向:12,軸方向:14

計算結果(載荷部分中央の節点番号は4)から $EI=35.24$N$\times$m$^{2}$

半円筒

• 周方向分割数:100, 高さ方向分割数:80
計算結果(載荷部分中央の節点番号は26)から $EI=-41.85$N$\times$m$^{2}$

自由端にヤング率$10^6$倍のビーム要素を付けたケース

計算結果(載荷部分中央の節点番号は26)から $EI=59.66$N$\times$m$^{2}$

雪平半円筒

• 周方向:4, 軸方向(数え方ややこしい):9

計算結果(最下部部中央の節点番号は440)から $EI=3.129$N$\times$m$^{2}$

自由端にヤング率$10^6$倍のビーム要素を付けたケース

計算結果(載荷部分中央の節点番号は440)から $EI=7.3331$N$\times$m$^{2}$

• 以下ccxs6s3yuki.f使用

• 周方向:2, 軸方向:5 ビーム要素つき

計算結果(最下部部中央の節点番号は126)から $EI=0.4697$N$\times$m$^{2}$

• 変形図(倍率50)

yuki25fh.frd

• 周方向:4, 軸方向:9 ビーム要素つき(ccxs6s3yuki.f使用)

計算結果(最下部部中央の節点番号は440)から $EI=7.722$N$\times$m$^{2}$

• 変形図(倍率200)

yuki49h.frd

• 周方向:8, 軸方向:19 ビーム要素つき

計算結果(最下部部中央の節点番号は1832)から $EI=24.10$N$\times$m$^{2}$

• 変形図(倍率200)

yuki819h.frd

• 周方向:16, 軸方向:39 ビーム要素つき

計算結果(最下部部中央の節点番号は7496)から $EI=48.96$N$\times$m$^{2}$

• 変形図(倍率200)

yuki1639h.frd

曲げ試験(FullCure720)

試験方法

• 単純梁にU字フック?を使って中央載荷。
• スパン中央のたわみの測定は難しいので少し位置をずらして測定。
• 試験体は厚さ5mm, 幅10mm, 長さ250mmで造形
• こちらの方法より軸方向ヤング率を導き、せん断弾性係数を算出する
• 荷重はフック・バケツ等を含め204g

横方向

スパン20cm(たわみ測定点は左支点より9cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.99^3}{12}=104$mm$^4$

• 曲げヤング率
こちらの式 $v=\frac{P}{48EI}(3L^2z-4z^3)$にそれぞれ代入して求める。

• 表裏平均:$1.62+1.66=1.64$GPa
• $(\frac{h}{\ell})^2=(\frac{4.99}{200})^2=6.23\times10^{-4}$

スパン18cm(たわみ測定点は左支点より8cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.97^3}{12}=102$mm$^4$

• 曲げヤング率

• 表裏平均:$1.58+1.62=1.60$GPa
• $(\frac{h}{\ell})^2=(\frac{4.97}{180})^2=7.62\times10^{-4}$

スパン16cm(たわみ測定点は左支点より7cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.96^3}{12}=102$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{1.54+2.18}{2}=1.86$GPa
• $(\frac{h}{\ell})^2=(\frac{4.96}{160})^2=9.61\times10^{-4}$

スパン14cm(たわみ測定点は左支点より6cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.01\times4.95^3}{12}=101$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{1.78+1.71}{2}=1.75$GPa
• $(\frac{h}{\ell})^2=(\frac{4.95}{140})^2=1.25\times10^{-3}$

スパン12cm(たわみ測定点は左支点より5cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.01\times4.94^3}{12}=101$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{1.38+1.16}{2}=1.27$GPa
• $(\frac{h}{\ell})^2=(\frac{4.94}{120})^2=1.69\times10^{-3}$

スパン10cm(たわみ測定点は左支点より4cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.01\times4.94^3}{12}=101$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{1.58+2.00}{2}=1.79$GPa
• $(\frac{h}{\ell})^2=(\frac{4.94}{100})^2=2.44\times10^{-3}$

スパン8cm(たわみ測定点は左支点より3cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.94^3}{12}=100$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{0.613+0.607}{2}=0.610$GPa
• $(\frac{h}{\ell})^2=(\frac{4.94}{80})^2=3.81\times10^{-3}$

スパン6cm(たわみ測定点は左支点より2cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.94^3}{12}=100$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{0.261+0.275}{2}=0.268$GPa
• $(\frac{h}{\ell})^2=(\frac{4.94}{80})^2=6.78\times10^{-3}$

スパン4cm(たわみ測定点は左支点より1cm)

• 日時

• 寸法(mm)

断面二次モーメント$I=\frac{10.00\times4.94^3}{12}=100$mm$^4$

• 曲げヤング率

• 表裏平均:$\frac{0.0791+0.0759}{2}=0.0775$GPa
• $(\frac{h}{\ell})^2=(\frac{4.94}{40})^2=1.53\times10^{-2}$

グラフ

縦軸:曲げヤング率逆数 横軸:$(\frac{h}{\ell})^2$ ※magehazureでスパン20cm〜12cmまでを外れ値とした。以下も外れ値とした値を除いて算出

引張試験(FullCure720)

試験方法

• 学生実験とほぼ同じ。
• 開始30秒後に荷重を載荷して、24時間計測。
• 24時間後に一旦計測停止、ファイル名を変えて再度計測開始。
• 開始30秒後に荷重を除荷して、24時間計測。
• 2ファイル間には若干の時間があるが、グラフをつなげて見たいので両ファイルをひとつにまとめた。

横方向

造形直後(10月9日15:30造形)

• 日時

• 寸法(mm)

• 荷重(g)

• 開始0秒のひずみで初期化(yoko1009s.csv)

• 開始30秒のひずみで初期化(yoko1009_wo29s.csv)

養生2日(10月9日15:30造形)

• 日時

• 寸法(mm)

• 荷重(g)

• 開始0秒のひずみで初期化(yoko1011s.csv)

• 開始30秒のひずみで初期化(yoko1011_wo29s.csv)

養生7日(10月9日15:30造形)

• 日時

• 寸法(mm)

• 荷重(g)

• 開始0秒のひずみで初期化(yoko1016s.csv)

• 開始30秒のひずみで初期化(yoko1016_wo29s.csv)

0014 1/2

• 開始0秒のひずみで初期化(yoko0112s.csv)

ひずみ比較

• 各試験体の時間毎のひずみ($\mu\epsilon$)比較
※時間は載荷後の時間。載荷は試験開始30秒後に行なっているので実際の計測時間は表示時間+30秒

• 各たてひずみ

• たてひずみ平均

• 各よこひずみ

• よこひずみ平均

ポアソン比比較

孔削り油差し試験体

• 横方向10/9造形試験体
• 日時

• 寸法(mm)

• 荷重(g)

• 開始121秒のひずみで初期化(yoko1022aburas.csv)

予備引張試験(FullCure720)

• およそ24時間1513g載荷後、除荷して24時間のひずみ(斜め方向)
• naname24_24s.csv

• 1kgを載荷して100000秒(27.777...時間)のひずみ(横方向)
• yoko27hs.csv

• 1kgを載荷して2時間のひずみ(横方向)
• yoko1kg2hs.csv

• 縦方向(3Dプリンターのヘッドが動く向きに対して試験体の軸が直角)試験体
載荷試験結果(csv)より各荷重の載荷1秒前のひずみを抜き取って描いた応力-ひずみ曲線ざっくりグラフ(引張正)
• gnuplotの最小二乗法による各ゲージのざっくりヤング率

viコマンド

• 文字数カウント

:%s/./&/g

論文

英語論文のタイトル(発声練習-英語論文のタイトルで大文字にするのはどこか-より)

• 単語の1文字目は大文字
• 接続詞,冠詞,前置詞は小文字

雑感

なぜ広まった？ 「『訊く』が正しい」という迷信より

最近ネット上で、小説から実用書まで幅広く書く人が「何故」「訳」「訊く」などを平仮名に直されるということについて愚痴を言っているのを読んだ。

「難しい」「堅い」「読めない」という理由を言われたという。

情けない話だ。

それが理由なら、ルビを振ればいいという話になってしまう。


それらを仮名書きするのには、原則がある。

和語の副詞や形式名詞、助詞に相当するものは平仮名にするというものだ。

「何故」を平仮名で書くのは和語の副詞だからで、「訳（わけ）」は「事」「為」と同じ形式名詞だから。

また、「訊く」は実用文にはふさわしくないスタイルで、何よりも「聞く」が伝統的な書き方だからだ。

どうして、それらの背景が言えないのか。

「大学の広報誌にカバーガールを載せるのはいかがなものか」を読んで

この広報誌を見かける(しっかり読んでない)たびに後藤さんが言うような

私が見た限りでは、男の学生が「表紙モデル」なったことはないし、 女の学生にしても、
どうもカバーガールとして見栄えのする かわいい外見の人を選択的に採用しているのではないかという 疑いを私は抱く。

カバーガールとして見栄えのする かわいい外見の人を選択的に採用

ひとつ思い出した。某ミスコンでいつからか外国人の候補者が登場し出した。ひとりだけ。 毎年チェックしてるわけじゃないからわかんないけど、もしやと思って某外国人候補者に きいてみたら、案の定、大学から頼まれたとのこと。まあ嘘ってことは考えにくい。 なんとか資源学部とかこういう事例とか、気持ち悪い。

anet「【入試課】 秋田大学「入試・イベント情報」 facebookについて」

アルミ缶の寸法

jacobianエラーについて

• 一行の長さを短くしてccx2.3を実行したところjacobianエラーは消えず。
短くした上でqflpコマンドで要素を反転させると反転させた要素番号のjacobianエラーは消えた。

• 一行の長さを短くしないでqflpコマンドで要素を反転させてもjacobianエラーは消えず。

g77をインストールしてccxs6.f90をg77でコンパイル(拡張子変更後)してinpファイルをつくった。

• 両者をccx_2.3で実行
ccxs6をgfortranでコンパイルしてつくったccxs6gfortran.inpはjacobianエラーが出るのに対し、 g77でコンパイルしてつくったccxs6g77.inp はjob finishした。

単純に両者の座標の桁数が違う。

Gmailのhtmlファイル非添付

• (新規作成/返信共に)本文欄すぐ上の「<<書式なしのテキスト」をクリックで自分が送るメールがテキスト形式のみになる。
戻したい場合は「リッチテキスト形式>>」をクリック。 参考

エラー

• ダイヤカットをgmshでメッシュをだんだん細かくしていくと、ある大きさからcalculixの
計算過程で

*ERROR in e_c3d: nonpositive jacobian
        determinant in element    要素番号

一応解決策はみつけたけど、メッシュが複雑すぎていちいち直すのは非現実的。 じゃあと思いnetgenを試してみた。 メッシュを切るまではうまくいく。

しかし、inpファイルに吐き出すと節点の途中までしか書いてない不完全なファイルが 出来上がってしまい計算に使えない。

メッシュの比較

静的解析

m=8,n=8

• gmsh法
• gmshで通常のメッシュ
ccxs688.dat

• Element size at pointsで全節点を選択->Value「10」
ccxs6g88.dat

• 旧法
daiyahanfc88.dat

メッシュを細かくしたら旧法の値に近づいた。 メッシュの細かさはかなり値に影響を与えるらしい。

座屈解析

m=8,n=8

• gmsh法
• element size10
ccxs6g88zaku.dat

• 旧法
daiya88zaku.dat

S6に変換

• 以下のリンクにgmshでメッシュを切ったあとcalculixで使うためにS6にinpファイルを変換するプログラムがある
http://dip28p.web.fc2.com/calculix/gmshinpconvert/index.html
• 上のプログラムを実行後、余計な要素を削除した後に
要素番号を付け直すプログラム ccxs6.f90

ダイヤカットをSALOMEでモデリング

• daiyahanf.f90をコンパイル・実行
• できたdaiyahanf.objの座標をSALOMEにコピペで節点を打っていく(7桁?)
• 節点を線で結ぶ(節点も線も順番は関係なし)
• New Entity->Build->Faceで三角形の面を作る
三角形を成す三つの線を選択するとObjects欄に3_objectsと表示されるので、Apply
• すべての面を作ったらOperations->Boolean->Fuseで2つの三角形をひとつにまとめる
• それを繰り返してひとつの図形にまとめる

Elmer

$ElmerGrid? 5 1 daiyahanf.inp

Starting program Elmergrid
Elmergrid reading in-line arguments
Output will be saved to file daiyahanf.
input of type 5 can't create output of type 1

Thank you for using Elmergrid!
Send bug reports and feature wishes to peter.raback@csc.fi

tex.vim??

"elsej.tex" 334L, 12997C
/usr/share/vim/vim73/syntax/tex.vim の処理中にエラーが検出されました:
行  793:
E115: 引用符 (') がありません: 'ς')
E116: Invalid arguments for function s:Greek('texGreek','\\varsigma\>'^I^I,'  �)
E15: 無効な式です: s:Greek('texGreek','\\varsigma\>'^I^I,'ς')
行  823:
E115: 引用符 (') がありません: '????)
E116: Invalid arguments for function s:SuperSub?('texSuperscript','\^','0','????)

E15: 無効な式です: s:SuperSub?('texSuperscript','\^','0','????)
行  827:
E115: 引用符 (') がありません: '????)
E116: Invalid arguments for function s:SuperSub?('texSuperscript','\^','4','????)

E15: 無効な式です: s:SuperSub?('texSuperscript','\^','4','????)
行  828:
E115: 引用符 (') がありません: '????)
E116: Invalid arguments for function s:SuperSub?('texSuperscript','\^','5','????)

E15: 無効な式です: s:SuperSub?('texSuperscript','\^','5','????)
行  829:

と表示された.

文字コードの問題だったみたい。 こちらを 参考にしたら改善された。

viで複数行に文字を挿入

選択を開始する行の一番上にカーソルを置き、 Ctrl+v 範囲を選択、 Shift+iで挿入モードへ # 入力(1文字だけ入力されます) Esc これで行頭に#が入力される。 コメント文にしたい時とか便利かも

座屈解析の覚書

• inpファイルのCLOADで載荷節点の倍率bを書き込む。daiyahanfは10$^{-8}$くらい
• ccxで計算して.datファイルのBUCKLING FACTOR を読む。
• BUCKLING FACTORの値に倍率をかけた値が座屈荷重。

無次元化

曲げ剛性

• 波形構造
この曲げ剛性 を無次元化してみた

niji9.fとniji13.fで半円筒と四分の一円筒の中立軸周りの断面二次をだす。それで daiyad3d.txtとdaiyad3d2.txtの各曲げ剛性を割った。

• 半円筒を無次元化してみた

座屈荷重

• 波形構造(l=0)
※円筒は周60高さ70

• 波形構造(l=1)
※円筒は周60高さ40

半ドーナツ型の断面二次モーメントを任意の軸周りで求められるプログラム

niji9.f

剛性を求めるプログラム

gousei.f

半ドーナツ型を縦分割して断面二次モーメントを求めるプログラム

niji7.f

矩形断面を縦横分割して断面二次モーメントを求めるプログラム

niji3.f

矩形を縦分割して断面二次モーメントを求めるプログラム

niji1.f

2011年度前期

2011年度後期