2月16-21日

今週もまた，どたばたした1週間だった。
・月曜日は，中村先生，国枝先生とともに道路局の報告会にった。委員長は，大物感たっぷりの人でした。G大学の先生がその前に発表されていて，ずいぶんとへりくだっている人でへんなプレゼンだなあと勝手におもっていたが，実は，ちょー大物ということで，これも驚いた。

・17日はJASS5NWGの報告書で半日つぶれた。

・19日は，名大の博士の公聴会，午後はファカルティーミーティング，その後は内緒の委員会，構造系で行った,非常勤講師の先生方との懇親会でした。
博士の公聴会，シナリオ不十分もあって，個人的に納得できない発表があった。審査委員会に入っているわけではないが，名大の博士ということでもう少しクオリティを考えてもよいと思う。

・20日は，広大でS氏とK氏の博士公聴会があって，出席した。S氏については，佐藤先生の配慮で副査に入れていただいた。どちらも社会人で5年以上の年月がかかっていた。公聴会には田澤先生をはじめ50名を超える方が聞きにいらっしゃっており，これぞ博士の公聴会，晴れの場という雰囲気があった。
毎度のことであるが，佐藤先生の出す博士のブランド力というのは，相当に高い感服する。

黄表紙3月号

黄表紙3月号は，超高強度の初期応力，初期欠陥に関する論文のその２が掲載される。
部材断面で，柱軸と直角方向の応力を検討した論文で，古くはJCI時間依存性の委員会のシンポジウムで発表した論文がスタートだった。
発見，の論文で示唆したものを解析・実験の両面から明らかにしたもの。

これで，自己収縮の部材への影響については一段落で，計6本の大作になった。広大で佐藤先生のご教授を得て以降の成果である。科研の若手(A)「超高強度コンクリートを用いた実大RC部材の検証」の成果ということで，審査していただいた方々にご理解いただければと思う。

・ひび割れ発見
・圧縮主筋の降伏問題
・初期応力・初期欠陥その１（夏・冬のひび割れリスクの相違）
・収縮の冬期が大きいことの詳細な検討（寺本筆頭）
・膨張材・収縮低減剤による鉄筋周囲のひび割れ抑制（鈴木氏筆頭）
・初期応力・初期欠陥その２（軸垂直方向応力，縦ひび割れ問題）

英文化もそろそろ検討しなくてはいけないな。

2月は，論文を黄表紙に2編投稿した。収縮関連。これは，材料的な奴で私が博士論文にとりかかってから，ずっと構想をねっていた論文なのでマニアックです。新しいシリーズの幕開けになります。
これも結構な大作になる予定。


超高強度関連の論文を1本提出してあったのだが，提出時にそうとうひっかかっていた。
結局，今月投稿した論文で自分の論文を否定することになってしまった。こんなに早くものごとがひもとけるとは思っていなかった。論文は取り下げをお願いしようと思う。査読委員会の方々にはお詫びしないと。


水和関連を早くとりまとめなくては，と思いつつもう１年が過ぎた。いい加減にしないと風化してしまうなあ。住宅局助成の報告書とりまとめもあるので，今月中には目途をつけたい。

採用決定

Elsevierの

Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment

に採用が決定しました。2月5日です。

含水人口軽量骨材から，低水セメント比のセメントペーストマトリクスに水分が移動する挙動を中性子らジオグラフィで可視化したもので，セメント硬化体の水和にともなう乾燥より前に，セメントの反応によるイオン濃度増大にともなう浸透圧で水分がでていそうだ，ということを実験的に確認した論文です。

いわれればみな当たり前と思うんですけどね。

でも，気づく人もすくなかったんですよ。

ちなみに自己養生効果は，骨材からというよりは骨材の周囲で水セメント比がリッチになった箇所から養生されるかんじです。

中性子の実験も，もう少しハンドリングが良いとよいんですけど，なかなか東海村は遠いし，原発のハードルはちょっと高い。

理科大，東大，原研４者の共同研究です。

最近のヒット : 文献整理

文献整理は非常に重要である。
論文を書くとき，モデルを検討するとき，実験結果を評価するとき，実験方法を考えるとき，すべての段階（なんか順番ぎゃくだけど）で必要。
で，最近は便利で論文を検索できるようになってはいるんだけど，それでも自分の中で，特に一度利用したことのある文献というのの整理は重要。

いろいろ，専用ソフトをかったり，出力してバインダーにしたりしたんだけど，結局完全版には至っていない。
現在，はじめた取り組みはPDFにして，メールソフトの中に入れ込む（あるいはG-Mailに放り込む）という手法。
メールに参考文献形式の記載と，キーワード，何の論文で引用されていたか，自分で何を引用したか，を記載すると，検索が容易。
ユーザーを研究室で共有する試みを行っているが，これはうまくいくかどうかわからない。

Beckyは検索機能も迅速なので，G-Mailからダウンロードして，検索してもよいしね。

というかんじで，新しい論文整理を初めてみた。しかし，文献名を入力するだけでも面倒くさいんだよなー。

昨日のつづき

あのあと考えてみました。
日本の個人資産総額1500兆円のうち，8割が高齢者のものだそうだ。
無利子国債の形にもっていくということのメリットは，それを短期的に吐き出すための施策である。
つまり，短期的に吐き出させるために無税という形で対応するということになるけれども，この部分のメリットが浮き出ていない点が問題なんでしょう。
利率が無い＋相続税が減額＋短期的にフローが出る

というときのフローが，どこにどう使われるから，どうなりますか，という形がもう少し明示的だといいんじゃないでしょうか。

昨日のリンク先のページにも，談話があるけど，「～いわれている」とか伝聞形で頼りない。
もう少し実行力のあるプランがあれば，その減税に対する効果っていうものがわかりやすくなるのだと思う。

http://news.livedoor.com/article/detail/4008551/


追記　文章少し変えました。　昨日のも同様。

無利子非課税国債

まあ，また，こんなことを書くとおこられちゃうかもしれないんだけども。
景気対策で，公共投資をするとしても，本来は将来性のあるものに利用しなくては，金の循環は行われないし，景気浮揚感もでてきこないと思う。

http://www.yomiuri.co.jp/politics/news/20090207-OYT1T00174.htm?from=main1

でも，人口が減るなかで経済成長率がどうのこうの，という議論で今回の振れ幅を見てみると，見てみると日本の経済成長率なんてほんと微々たるもので，それがアメリカバブルの影響でしたなんて哀しくて見ていられない。今の日本はその変動幅自体が小さくなっているので，結局短期的な解決というのは，意味が無い。

むしろ，お金を使うので有れば，長期的な経済成長を促す枠組みが不況であっても重要なんじゃないか，というのが持論です。

そういえば，まったく新しい取り組みとして，ベーシックインカムの議論はおもしろい。
（リンク1，2）
昨年，ホリエモンのブログから，山崎元氏のページへのリンクで少しだけ話題になった。時代がこういうものになっている，という視点は重要で，政策や政治を見る視点が70年代のつねに仕事があふれていて，常に税収が増えていって，常に人口が増えているんですというもので固定されていることから脱却することは重要でしょう。

マスコン講習会

JCIのマスコン講習会で広島に行きました。
今回は、照査事例ということで、計算事例についてのご説明。
おかげさまで、今回のもので講習会における私の役目はごめんとなりました。
札幌もいきたかったのですけども、残念ながら行けません。

しかし、広島駅は近いなあ。東広島は遠い・・・。時間を買うのなら、広島行ってから戻った方がよいですね。

Gmail オフライン機能の利用

すでに英語でオフライン機能が実装されていますが，これも日本語で利用できるようです。
http://google-mania.net/archives/1127

しかし，G-mailのβ版っていうのは，いつ正式版になるのでしょうかねえ・・・。

コンパイラ

今日は，家でずっとプログラムをいじっていた。
Intelコンパイラの最適化オプションでCPU依存のパラメータがあるのだが，大学と家でCPUが異なることをわすれていてそのままにしていた。
なんか，CCBMの値が感覚的にちょっと変だな，と思って全部チェックしていたら，最適化オプションがCPUに対応していない場合に変な値を吐き出す関数があることがわかった。

そういわれれば，そうなんだけどあまりCPUが変わった場合というのを意識していなかったので，今回は焦った。というよりも，時間を6時間ほどつかってしまってしょんぼり。

やれやれ

バッファロー「PF-50WG」

http://dc.watch.impress.co.jp/cda/item/2009/01/30/10086.html

デジタルフォトフレーム。無線LAN付。
これはよさそうだ。

ミニモニタ

http://www.amazon.co.jp/exec/obidos/ASIN/B001QVJ83C/gj_moji-22/ref=nosim/
USB接続らしい。これは便利そうだ。

ただ，モニタ自体安くなっているので，場所に困るとか，延長して別のところで見せたりとか，ちょっと変わった利用の仕方かな。
あるいはサーバーとか？

実験棟などでは，普段はPC単体でデータログして，確認のためにこれに繋いで見るっていうのはありかもしれないな。小型PCも中古だといくらでも手にはいるので，HDDを換装すれば信頼できるようになるし。

非対称，構造対称マトリクスの解法 MKL PARADISO

備忘録。
水熱連成で，非対称マトリクスになってしまう場合について。
MKLのParadisoで非対称マトリクスを解くためのサブルーチンの一例。

サンプルコード：

INCLUDE 'mkl_pardiso.f90'
subroutine MainMKL()
use MKL_PARDISO
real(8) :: A(5,5)
real(8) :: val(13)
integer :: col(13)
integer :: ptrB(5)
integer :: ptrE(5)
integer :: m
real(8) :: alpha

!MKL PARAMETERS
type(mkl_pardiso_handle) :: pt(64)
integer :: maxfct=1, mnum=1, mtype=11, phase=13, n, nrhs=1, error, msglvl=0
integer :: iparm(64)
integer :: nnz
integer :: ncol
integer :: idum
real(8) :: ddum
real(8), allocatable :: value(:)
real(8), allocatable :: b(:)
real(8), allocatable :: x(:)
integer, allocatable :: columns(:)
integer, allocatable :: rowindex(:)
integer, allocatable :: perm(:)

integer :: i, j, k
integer :: nzidx
integer :: MAX_OF_NODES

pt(:)=mkl_pardiso_handle(0)
iparm(1)=0

A=0.D0
A(1,1)=1.D0
A(1,2)=-1.D0
A(1,4)=-3.D0
A(2,1)=-2.D0
A(2,2)=5.D0
A(3,3)=4.D0
A(3,4)=6.D0
A(3,5)=4.D0
A(4,1)=-4.D0
A(4,3)=2.D0
A(4,4)=7.D0
A(5,2)=8.D0
A(5,5)=-5.D0


!n x nを定義
MAX_OF_NODES=5
n=MAX_OF_NODES


!rowIndexを定義
allocate(rowindex(1:(n+1)))
allocate(b(1:n))
allocate(x(1:n))
allocate(perm(1:n))
rowindex(1)=1
nnz=0
do i = 1, n
ncol=0
do j = 1, n
if(A(i,j)/=0.D0) nnz=nnz+1
if((A(i,j)/=0.D0)) ncol=ncol+1
enddo
rowindex(i+1)=rowindex(i)+ncol
enddo
allocate(value(1:nnz))
allocate(columns(1:nnz))

nzidx=0
do i =1, n
do j=1, n
if(A(i,j)/=0.D0)then
nzidx=nzidx+1
value(nzidx)=A(i,j)
columns(nzidx)=j
endif
enddo
enddo

do i=1, nnz
write(*,*) value(i), columns(i)
enddo
!-------------------------------------------------------------
!Factorization & Solution
!-------------------------------------------------------------
phase = 13 ! only factorization

!ここで，OutFvecを代入
do i = 1, n
b(i) = 1.D0
end do

call pardiso (pt, maxfct, mnum, mtype, phase, n, value, rowindex, columns,&
& perm, nrhs, iparm, msglvl, b, x, error)
write(*,*) 'error', error
do i = 1, n
write(*,*) ' x(',i,') = ', x(i), b(i)
end do


b=matmul(A, x)
do i = 1, n
write(*,*) ' b(',i,') = ', b(i)
end do

!-------------------------------------------------------------
!Termination and release of memory
!-------------------------------------------------------------
phase = -1 ! release internal memory
CALL pardiso (pt, maxfct, mnum, mtype, phase, MAX_OF_NODES, value, rowindex, columns,&
& perm, nrhs, iparm, msglvl, b, x, error)
return
end subroutine MainMKL

ぼー

なんというか，学生の卒論，修論データをチェックしていると，本当にいろいろなことに気づく。
まさしく，研究テーマの宝の山なんですが，あまりにおおくて，メモしていても忘れてしまう。
1年前の今頃，この実験はかならずやろう，と思っていたことを思い出したり，なんでこのデータをこのグラフで今まで整理してこなかったのか，と本当に愕然とするような事態にも遭遇する。

あと，うちでは，乾燥装置（11％，105度，冷凍）を整備したけれども，実験の材齢のデータで水和をとめておくことは本当に重要。11％RHは長期的にやるとC-S-Hが形かわっている感じがするけど，その他のデータは半年くらいなら融通が利く。東北のH先生には，ドモホルン方式だね，といわれたけど，追加的に確認実験をして，整合性がとれたり，仮説が立証されると本当に小躍りする。


しばらくは，表面関係の本を読みあさっていたけれども，また，水分移動，熱移動の話にもどってきた。今年は，JCIで熱勾配による水分移動を熱伝導率の含水率依存性から評価する論文を岸君が出した。ペーストの熱伝導率では，含水率が大きいときは，20％程度も水分移動で熱が運ばれることを実験的に示すとともに，熱勾配の物性値をある指標から定量的に予測できることも示した。
水熱連成大事だね，という前置きの論文です。

で，20日から，水熱の完全連成プログラムをつくっている。完全に連立するとマトリクスが対称にならないので，いろいろ思案していたけども，結局MKLの非対称アルゴリズムで対応することにした。来年度の科研があたれば，水・熱・応力完全連成を目指したいところです。

また，JCIでは4年生二人ががんばって，論文を投稿した。一人はセメント硬化体の温度履歴依存性について，相組成と反応率から説明する論文。収縮理論と絡めるとその大事さがわかると思うけれども，そこが記載できなかったので，ちょっと唐突な論文になってしまったかも。

もう一人は，高炉セメントの水和で温度履歴の影響。
彼の研究は，投稿には間に合わなかったけど石膏やアルミの挙動とひずみの関係について，ある程度の知見が得られた。石膏を多く添加することの意義，というのがあきらかになったと思っている。

寺本は，相変わらず超高強度のデータを示した。
私は，マスコンのひび割れ幅の式です。JCI指針のものを投稿させていただいた。

そろそろ，科研および国交省の助成のおわりが近づいているので，いままでのデータをまとめて論文・報告書を書かないといけないなあ。
頭にある奴を早くはき出したい。