Microsoft Word - 17_SKYACTIV-Gノッキング検出性のMBD_1016改.doc
|
|
|
- ひろみ くまじ
- 7 years ago
- Views:
Transcription
1 No.1(01) 特集 : モデルベース開発 17 KYACTIV-G ノッキング検出性の MBD MBD for Knock Detection in KYACTIV-G 要約 阿部広幸 *1 森本博貴 * 住谷章 * Hiroyuki Abe Hiroki Morimoto Akira umitani 荒川博之 *4 井上淳 *5 Hiroyuki Arakawa Atsushi Inoue KYACTIV-G では熱効率を徹底的に追求し,14.0 という高圧縮比を採用した (1) この高圧縮比の採用にはノッキングに代表される異常燃焼の回避が重要であり,Knock Control ystem を用いた点火時期制御において, 従来と比べ高いノッキング検出性に基づく緻密な制御が不可欠であった 今回, ノッキング発生時のセンサ振動状態を燃焼室空洞共鳴とエンジン構造体の連成応答解析で模擬し, エンジン構造変更によるノッキング検出性への影響を定量的に予測する手法を開発したので紹介する 本技術の適用により机上検討の段階でノッキング検出性のポテンシャル見極めが可能となり, エンジンの点火時期制御開発を早期化した ummary In the pursuit of thermal efficiency improvement, KYACTIV-G has taken up a high compression ratio of The avoidance of irregular combustion such as knocking was essential to the adoption of a high compression ratio, and even for ignition timing control which utilizes the Knock Control ystem more precision was necessary, requiring high knock detection accuracy. ensor vibration during knocking was simulated by coupled response analysis of the combustion chamber resonance and the engine structure. The quantitative prediction of the effects of changes in engine structure on knock detection. This technology enables knock detection ability to be assessed during the preliminary verification stages, accelerating the selection of engine specifications which enhance the ability to control ignition timing. 1. はじめに KYACTIV-G では, 熱効率を徹底的に追求し, 高圧縮比 14.0 を採用している しかし, 高圧縮比化を行うとノッキング ( 以下, ノック ) やプリイグニッション等の異常燃焼が起こりやすくなる これを回避するための制御による無駄なトルク低下なしに, 高圧縮比化のメリットを最大限生かすため, ノックセンサで検出する振動でノック発生を判定する KC(Knock Control ystem) を開発してきた ノックによるトルク低下 () およびプリイグニッション回避技術 () については, 既に技報にて紹介しているため詳細は割愛する この KC により点火時期セットを大幅に効率化し * 1~4 エンジン性能開発部 Engine Performance Development Dept. てきた 一方でエンジン開発時にノック強度の代用特性とするため設置した燃焼室圧力センサ値と, ノックセンサ振動値の相関が弱いため, エンジン構造変更がノック検出性に与える影響が定量的でなく, エンジン制御キャリブレーションに多くの工数を要していた そこで今回, 燃焼室 -エンジンの連成応答解析により机上でノック発生時におけるノックセンサ部振動を模擬する手法を開発した これによりノックセンサ振動に対するノック加振力や, 燃焼室からセンサ部への振動伝達特性に対する多くの制御因子を明らかにし, 課題であった燃焼室圧力センサ値とノックセンサ振動値の関係を解明した 以上で得られた知見を基にノック検出性の MBD を行ったのでここで紹介する * 5 パワートレイン技術開発部 Powertrain Technology Development Dept. 91
2 No ノック発生 点火時期制御の流れ Powertrain ノックが発生すると Fig. 1 のように 燃焼室内の空 洞共鳴を励起し 増幅した圧力波が燃焼室壁面を加振 することでエンジン全体が振動する この振動 加速 度 をシリンダブロック吸気側壁面にあるノックセン サで検出している Combustion Chamber KC は上記振動を周波数分析することで 気筒 サ Fig. イクルごとに特定周波数帯での強度を算出し 振動強 Finite Element Model. 解析手法 度の頻度分布を推定している その分布の標準偏差な 計 算 は モ ー ダ ル 周 波 数 応 答 解 析 NATRAN どから算出した閾値 Kσ を振動強度が超えた場合 olution 111 で燃焼室 PT ユニットを連成させて行 ノック発生と判定する った ここでは振動の伝達特性を評価するため 体積 Resonance of Combustion Chamber 加速度 球面放射の加速度 で単位加振を行い 燃焼 Engine Vibration Due to Cylinder 室 PT ユニット構造体でノックに対する振動応答感度 を評価する 加振点は燃焼室内のノック発生点とした ノックは 主にエンドガスゾーン シリンダ外周部の未燃混合気 部 で発生すると考えられており 本解析ではエンド ガスゾーン全周でノックが発生すると想定した Fig. Knock ensor 応答評価は以下で行う (1) 燃焼室指圧計測用センサ部の圧力 Frequency (%) Vibration Intensity Distribution of Vibration Intensity Threshold (Kσ) Combustion Cycle () ノックセンサ位置の振動 加速度 Fig. 4 Knock Excitation Points Vibration Intensity Fig. 1 Knock Detection Flow (1). 解析手法.1 解析モデル 解析は燃焼室 燃焼ガス PT Powertrain ユ Fig. Load and Response Points of Combustion Chamber ニットを有限要素でモデル化して行う Fig. ここでいう燃焼室はノックが発生するクランクアン () Acceleration グルにおける容積をモデル化する 解析で使用する燃 焼ガス物性値は 燃料噴射量と燃焼室容積から求めた 質量密度と式(1)から求めた音速を使用する P α 音速 κ 比熱(1.5) ρ 質量密度 P 燃焼指圧 (1) Fig. 4 Response Point of Knock ensor. 結果処理方法 有限要素長は ノック発生周波数帯の燃焼室空洞共 鳴 エンジン振動状態を十分表現できる長さとする 結果処理は制御性評価周波数帯のノックセンサ加速 今回は平均メッシュサイズ 8mm 約 00 万自由度で 度 燃 焼 室 圧 力 値 そ れ ぞ れ で RM Root Mean モデル化を行った quare 値をとり 燃焼室圧力を横軸 ノックセンサ 加速度を縦軸として加振点ごとに ここでは 8 点 プ ロットする Fig. 5 9
3 No.1(01) Acceleration 0 Frequency (Hz) Knock ensor Output Plots per Load Frequency (Hz) ensor Output Fig. 5 Treatment of Calculated Data ここで Pmax( 時間軸での指圧最大値 ) 付近で定義された燃焼室形状において 1 次の空洞共鳴は Front-Rear モードと Intake-Exhaust モードがあり, これらはほぼ同じ周波数帯に存在する 例えば Front もしくは Rear 側でノックが発生した場合,Front-Rear モードを主に励振することから指圧計測点は応答の腹となり, 圧力値は大きくなる (Fig. 6(b)) 一方,Intake もしくは Exhaust 側でノックが発生した場合,Intake-Exhaust モードが主に励振され, 指圧計測点は応答の節となり, 圧力値は小さくなる (Fig. 6(a)) この結果より, ノック強度が一定でもノック発生位置により指圧センサ圧力はおよそ 8 倍も変わることを示している よって冒頭で述べた燃焼室圧力センサ値とノックセンサ振動値との相関が弱いメカニズムの解明ができた Node at ensor (a) Exhaust Intake Knock Point (c) approx. 8 times Front (b) Anti-Node at ensor Rear 前述のようにノック発生点 ( 加振点 ) ごとに励起される空洞共鳴が変わり, それによりセンサ振動が変化するため, 加振点をエンドガスゾーン全周として評価すると Fig. 6(c) のような点群となる 更にノックは Pmax から以後 10 あたりの特定 CA(Crank Angle) 間で発生するため, それぞれの燃焼室状態において Fig. 6 のような点群を求める必要がある Fig. 7 に例として, つの点火時期状態を評価するために行った下記 4 条件 CA の解析結果を示す 1 弱ノック発生点火時期から 1 点火進角した場合の Pmax 時の CA 上記進角時の Pmax 後 10 の CA 弱ノック発生点火時期から 5 点火遅角した場合の Pmax 時の CA 4 上記遅角時の Pmax 後 10 の CA ここで,4はノック発生タイミングばらつきの最遅角側を表している P-θ 1 4 Vibration- 1 4 Advanced 1 Retarded Crank Angle from TDC (deg) Fig. 7 Calculated Data for Each Crank Angle Fig. 7 点群分布において,1 進角状態 ( 紺色点 ) では CA5 (1) と 5 () 間に分布が存在し,5 遅角状態 ( マゼンタ点 ) では CA5 () と 45 (4) 間に分布が存在すると仮定し,No.1~4 気筒ごとに実機結果と CAE 点群分布比較を行うと, 単位加振の解析結果にもかかわらず実機同様の特徴を再現しており, 実機との相関があると判断した (Fig. 8) 0.E+00.E-04 4.E-04 6.E-04 8.E-04 1.E-0 (Pa) Fig. 6 ensor by Loaded Point 9
4 No.1(01) Cal. Cylinder No. 1 Exp. Cylinder No. 1 値 ) を超えた時をノック発生と判定するため, 応答感 Upward rise Cal. Cylinder No. ector High level Cal. Cylinder No. ector Cal. Cylinder No. 4 Exp. Cylinder No. Exp. Cylinder No. Exp. Cylinder No. 4 度ばらつきはない方がノック判定閾値は安定し, 更に純粋なノック強度を評価できるため, ノック検出性としては理想である (Fig. 9(b)) CA 5 CA 5 (a) Vibration- (b) Ideal Characteristic Fig. 9 Interpretation of Response ensitivity Oval Fig. 8 Calculation/Experiment Vibration- 従ってこの振動 - 圧力分布結果を基にノック検出性の評価指標を定めた 4. ノック検出性評価手法 KC の点火時期制御では以下 つの現象が課題となる (1) ノックを判定できないため, 必要な遅角を行わずノック音が聞こえる () ノックを誤判定し, 不要な遅角を行ってしまいトルクが低下する これら現象はノック判定閾値に対するノックセンサ値の関係で決まるが, これに対して Fig. 8 の振動 圧力分布から理想特性を考えると以下のようになる Fig. 9(a) のように振動 - 圧力分布の縦軸 ( ノックセンサ振動 ) に着目すると同一加振力にも関わらずノック発生点, ノック発生タイミング CA5,CA5 で縦軸値が の幅でばらつく つまりノック発生点 -ノックセンサ振動への応答感度が だけばらつくことを示している 章 ノック発生 ~ 点火時期制御の流れ で述べたが,KC は振動強度の頻度分布を推定し,Kσ( 閾 4.1 定常時のノック検出性評価上記理想特性を踏まえ 1 進角状態のノック検出性を評価する場合,Fig. 10 のように分布の縦軸のばらつき幅を各気筒で評価し, ばらつき小をノック検出性良, ばらつき大をノック検出性悪と定義する Acceleration Range Fig. 10 Vibration- CA5 CA5 Deviation 1 4 Cylinder Number Bad Good Evaluation of Knock Detection During teady tate 4. 過渡時のノック検出性評価ある気筒でノック判定した場合, 点火時期が一気に遅角するので, 遅角前と評価する CA, 燃焼ガス物性が変わる それによる燃焼室圧力分布の変化からノックセンサ部への応答感度も大きく変化する しかしながら遅角時において通常ノックは発生しない ( 粗悪燃料使用時, 吸気高温時など除く ) ため, 応答感度が下 94
5 No.1(01) がる場合は問題ではない 一方, 応答感度が大きく上昇する場合, 遅角切り替わりでセンサ振動が上昇し, ノック判定閾値も上昇することでノック検出性が悪化する可能性がある よって遅角切り替わり前後のノック検出性評価も必要であり, Fig. 11 における遅角前 ( 進角時 ) 分布の最小値と遅角時の最大値の幅 をばらつきとして各気筒を評価する 著者 Vibration- max. when retarded ' min. when advanced 阿部広幸森本博貴住谷章 Acceleration Range ' Deviation 1 4 Cylinder Number Fig. 11 Evaluation of Knock Detection During Transient 上記のように定常時 (4.1 節 ) と過渡時 (4. 節 ), つの場合のノック検出性を評価することで, 対象ユニット構造のノック検出性優劣を机上で評価できる 5. おわりに 今回開発した技術により, ノックによるセンサ部の振動発生メカニズムを明らかにし, 机上でノック検出性を評価する指標を確立した この指標からノック検出性に優れた PT ユニット仕様検討が可能となり, 高精度な点火時期制御を実現し, ノック音とエンジントルクを高次元でバランスしたエンジンを短期間で開発するプロセスを構築できた 今後は次世代のエンジンに対しても本技術適用を進め, 更にお客様に喜ばれる高品質な商品開発に貢献していく所存である 参考文献 Bad Good (1) 富澤和廣ほか : KYACTIV-G 新型ガソリンエンジンの紹介, 01 年自動車技術会シンポジウム (01) () 長谷川裕一ほか : CX-5 KYACTIV-G のエンジン技術,,No.0,pp. 8(01) () 江角圭太郎ほか : KYACTIV-G 制御技術の紹介,,No.9,pp.6 40(011) 荒川博之 井上淳 95
研究論文 ガソリン高圧噴射を用いた高圧縮比エンジンの燃焼技術 ( 第 3 報 ) * 混合気制御による熱効率改善の検討 神長隆史 1) 養祖隆 1) 長野高皓 1) 藤川竜也 1) 山川正尚 1) Combustion Technologies of High Compressio
研究論文 * 混合気制御による熱効率改善の検討 神長隆史 養祖隆 長野高皓 藤川竜也 山川正尚 Combustion Technologies of High Compression Ratio Engine Using High Pressure Gasoline Injection (Third Report) - Study on Improvement of Thermal Efficiency
Development of Induction and Exhaust Systems for Third-Era Honda Formula One Engines Induction and exhaust systems determine the amount of air intake
Development of Induction and Exhaust Systems for Third-Era Honda Formula One Engines Induction and exhaust systems determine the amount of air intake supplied to the engine, and as such are critical elements
FFT
ACTRAN for NASTRAN Product Overview Copyright Free Field Technologies ACTRAN Modules ACTRAN for NASTRAN ACTRAN DGM ACTRAN Vibro-Acoustics ACTRAN Aero-Acoustics ACTRAN TM ACTRAN Acoustics ACTRAN VI 2 Copyright
03マツダ技報29本文-h1-01-1案0620
特集 :SKYACTIV TECHNOLOGY 7 Introduction of SKYACTIV-G Control Technology *1 江角圭太郎 佃 *2 厚典 *3 今村友之 Keitaro Ezumi Atsunori Tsukuda Tomoyuki Imamura *4 遠藤孝次 *5 上野隆司 小森 Kouji Endou Takashi Ueno Satoshi Komori
5 11 3 1....1 2. 5...4 (1)...5...6...7...17...22 (2)...70...71...72...77...82 (3)...85...86...87...92...97 (4)...101...102...103...112...117 (5)...121...122...123...125...128 1. 10 Web Web WG 5 4 5 ²
Microsoft Word - 校了_03_新型SKYACTIV-G気筒休止システムNVH性能開発.docx
3 SKYACTIV-G 2.5 気筒休止システムの NVH 対応技術 Vehicle NVH Performance Technology for Cylinder Deactivation of SKYACTIV-G 2.5 濱野健太 *1 田中和宏 *2 陳宇 *3 Kenta Hamano 土路生修 *4 Osamu Torobu Kazuhiro Tanaka Yu Chen 要約 新型
第6章 実験モード解析
第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法
国土技術政策総合研究所 研究資料
3. 解析モデルの作成汎用ソフトFEMAP(Ver.9.0) を用いて, ダムおよび基礎岩盤の有限要素メッシュを8 節点要素により作成した また, 貯水池の基本寸法および分割数を規定し,UNIVERSE 2) により差分メッシュを作成した 3.1 メッシュサイズと時間刻みの設定基準解析結果の精度を確保するために, 堤体 基礎岩盤 貯水池を有限要素でモデル化する際に, 要素メッシュの最大サイズならびに解析時間刻みは,
Microsoft PowerPoint - ce07-13b.ppt
制御工学 3 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価キーワード : 設計手順, 性能評価 8. ID 補償による制御系設計キーワード : ( 比例 ),I( 積分 ),D( 微分 ) 8.3 進み 遅れ補償による制御系設計キーワード : 遅れ補償, 進み補償 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. ループ整形の考え方を用いて, 遅れ補償,
F7-10 エンジンの Design of F7-10 High Bypass Turbofan Engine for P-1 Maritime Patrol Aircraft 空 部 ス 部 空エンジン 部 空 部 ス 部 空エンジン 部 F7-10 エンジン の P-1 の ファン エンジン 部
F7-0 Design of F7-0 High Bypass Turbofan Engine for P- Maritime Patrol Aircraft 空 空 空 空 F7-0 P- ファ ( : ) 00 量 量 F7-0 6 tf F7-0 The F7-0, a turbofan engine for the P- maritime patrol aircraft operated by
圧電型加速度センサ Piezoelectric Acceleration sensor 特長 Features 圧電素子に圧電型セラミックを用いた加速度センサは 小型堅牢 高感度で広帯域を特長としております 従って 低い周波数の振動加速度から衝突の様な高い加速度の測定まで 各分野で 幅広く使用されて
圧電型加速度センサ 小型タイプ φ3.5 5.85 2.5(H)mm 加速度 MAX100,000m/s 2 高温タイプ MAX250 小型 3 軸タイプ 8 7 5.5(H)mm Super miniature type φ3.5 5.85 2.5(H)mm 100,000m/s 2 High temperature resistance type MAX250 and Triaxial type
技術解説_有田.indd
Acceleration / G 2 18 16 14 12 1 8 6 4 2 Damping : 1. Period / s XY.1.1 1. 6533 283 3333 423 155 15 (X) 26.12 Hz 15 12 (Y) 28.32 Hz (Z) 43.98 Hz GS Yuasa Technical Report 211 年 6 月 第8巻 水平方向 X_3G 1.7e+7
75 unit: mm Fig. Structure of model three-phase stacked transformer cores (a) Alternate-lap joint (b) Step-lap joint 3 4)
3 * 35 (3), 7 Analysis of Local Magnetic Properties and Acoustic Noise in Three-Phase Stacked Transformer Core Model Masayoshi Ishida Kenichi Sadahiro Seiji Okabe 3.7 T 5 Hz..4 3 Synopsis: Methods of local
多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典
多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典 重回帰分析とは? 重回帰分析とは複数の説明変数から目的変数との関係性を予測 評価説明変数 ( 数量データ ) は目的変数を説明するのに有効であるか得られた関係性より未知のデータの妥当性を判断する これを重回帰分析という つまり どんなことをするのか? 1 最小 2 乗法により重回帰モデルを想定 2 自由度調整済寄与率を求め
Microsoft PowerPoint - chap8.ppt
第 8 章 : フィードバック制御系の設計法 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価 キーワード : 設計手順, 性能評価 8. 補償による制御系設計 キーワード : ( 比例 ),( 積分 ),( 微分 ) 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. 学習目標 : 補償の有効性について理解し, その設計手順を習得する. 第 8 章
no15
Development of High Performance Catalyst Temperature Sensor for NOx Catalyst Control Atsushi KURANO Kaoru KUZUOKA Sotoo TAKAHASHI Itsuhei OGATA In order to meet each countrys low emission vehicle regulations
第62巻 第1号 平成24年4月/石こうを用いた木材ペレット
Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting
Microsoft Word - 14_SKYACTIV-D i-stop開発における MBD適用事例の紹介_1018改.doc
特集 : モデルベース開発 14 SKYACTIV-D i-stop 開発における MBD 適用事例の紹介 Introduction of MBD for SKYACTIV-D i-stop Development 要約 中本仁寿 *1 中島宏幸 *2 朝倉浩之 *3 Yoshihisa Nakamoto Hiroyuki Nakashima Hiroyuki Asakura 寺岡陽一 *4 田賀淳一
Table 1. Reluctance equalization design. Fig. 2. Voltage vector of LSynRM. Fig. 4. Analytical model. Table 2. Specifications of analytical models. Fig
Mover Design and Performance Analysis of Linear Synchronous Reluctance Motor with Multi-flux Barrier Masayuki Sanada, Member, Mitsutoshi Asano, Student Member, Shigeo Morimoto, Member, Yoji Takeda, Member
スライド 1
ANDROMEDA Andromeda Simulation system 背景 自動車開発プロセスにおける試験機の位置づけと業界の動向 先行開発 製品開発 製品開発時の コンポーネント開発メニューとして提供する車両テスト実走行テスト E/G T/Mなどパワートレイン耐久 性能 適合試験などを行う完成車の試験を行う完成車の実走行試験を行う ECUなどエレクトロニクス関連 エンジン単体試験 耐久 性能試験
研究論文 HCCI 機関の急峻な燃焼に関する実験的研究 * 須山謙太 1) 伊藤直也 2) 寺島昂 3) 東條智也 4) 飯島晃良 5) 吉田幸司 6) 庄司秀夫 7) Experimental Research on Rapid Combustion in HCCI Engine
研究論文 1369 * 須山謙太 1) 伊藤直也 ) 寺島昂 3) 東條智也 ) 飯島晃良 5) 吉田幸司 6) 庄司秀夫 7) Experimental Research on Rapid Combustion in HCCI Engines. Kenta Suyama Naoya Ito Akira Terashima Tomoya Tojo Akira Iijima Koji Yoshida
Microsoft Word - 校了_08_新型 SKYACTIV-G 2.5 燃焼技術の開発.docx
No.34(2017) マツダ技報 特集 : 新型 CX-5 8 新型 SKYACTIV-G 2.5 燃焼技術の開発 Gasoline Engine SKYACTIV-G 2.5 Combustion Technology 小谷佳苗 *1 上村匠 *2 平下茂行 *3 Kanae Kodani 西尾貴史 *4 Takafumi Nishio Takumi Uemura 野村健太郎 *5 Kentaro
Microsoft PowerPoint - ⑥H20説明資料(岡山大河原)配布 ppt
点火プラグ実装型燃料 残留ガス濃度 計測センサシステム 岡山大学大学院自然科学研究科産業創成工学専攻教授冨田栄二准教授河原伸幸 ガソリンエンジンへの要求 高い熱効率 ( 低燃費 ), クリーンな排気ガスノッキング燃焼のサイクル変動 エンジンシリンダ内での特性評価 点火プラグ近傍の燃料濃度 ( 空燃比 ) 赤外吸収法水素炎イオン検出器 (FID) 点火プラグ近傍の乱流特性点火プラグ埋込型 LDV 装置
No.33(216) 2. ガソリンターボエンジンの制御技術 2.1 ターボシステム SKYACTIV-G 2.5Tの制御システムとして, 電動式のウェイストゲートバルブとエアバイパスバルブをマツダで初採用した どちらのバルブも, 従来の機械式では吸気管内の圧力設定値によって作動させるが, 電動式は
特集 : 新型 CX-9 4 新型 2.5L ターボチャージャーエンジン制御技術の紹介 Introduction of New 2.5L Turbocharger Engine Control Technology 西尾貴史 *1 東尾理克 *2 砂流雄剛 *3 Takafumi Nishio Masayoshi Higashio Yugou Sunagare 足利謙介 *4 坊田亮郎 *5 Kensuke
EPA Tier4排出ガス規制対応2~3tフォークリフト向けディーゼルエンジンの開発,三菱重工技報 Vol.54 No.1(2017)
M-FET 特集技術論文 49 EPA Tier4 排出ガス規制対応 2~3t フォークリフト向けディーゼルエンジンの開発 Development of EPA Interim Tier4 Certified Diesel Engines for 2-3t Forklift Trucks *1 山本高之 *2 平岡賢二 Takayuki Yamamoto Kenji Hiraoka *2 入口信也
kubostat2017b p.1 agenda I 2017 (b) probability distribution and maximum likelihood estimation :
kubostat2017b p.1 agenda I 2017 (b) probabilit distribution and maimum likelihood estimation [email protected] http://goo.gl/76c4i 2017 11 14 : 2017 11 07 15:43 1 : 2 3? 4 kubostat2017b (http://goo.gl/76c4i)
<4D F736F F D F838D815B835E838A CC94528FC4975C91AA8B5A8F708A4A94AD5F E646F63>
論文 解説 27 ロータリエンジンの燃焼予測技術開発 Analysis Technology of Combustion Model in a Rotary Engine 要約 中島聖 *1 尹正虎 *2 本郷均 *3 Akira Nakashima Jyong-Ho Yun Hitoshi Hongou 植木信也 *4 香川良二 *5 横畑英明 *6 Shinya Ueki Ryoji Kagawa
untitled
1.0 1. Display Format 8*2 Character 2. Power Supply 3.3V 3. Overall Module Size 30.0mm(W) x 19.5mm(H) x max 5.5mm(D) 4. Viewing Aera(W*H) 27.0mm(W) x 10.5mm(H) 5. Dot Size (W*H) 0.45mm(W) x 0.50mm(H) 6.
Table 1. Assumed performance of a water electrol ysis plant. Fig. 1. Structure of a proposed power generation system utilizing waste heat from factori
Proposal and Characteristics Evaluation of a Power Generation System Utilizing Waste Heat from Factories for Load Leveling Pyong Sik Pak, Member, Takashi Arima, Non-member (Osaka University) In this paper,
< B837B B835E82C982A882AF82E991CF905593AE90AB8CFC8FE382C98AD682B782E988EA8D6C8E40>
1 / 4 SANYO DENKI TECHNICAL REPORT No.10 November-2000 一般論文 日置洋 Hiroshi Hioki 清水明 Akira Shimizu 石井秀幸 Hideyuki Ishii 小野寺悟 Satoru Onodera 1. まえがき サーボモータを使用する機械の小型軽量化と高応答化への要求に伴い サーボモータは振動の大きな環境で使用される用途が多くなってきた
RLC 共振回路 概要 RLC 回路は, ラジオや通信工学, 発信器などに広く使われる. この回路の目的は, 特定の周波数のときに大きな電流を得ることである. 使い方には, 周波数を設定し外へ発する, 外部からの周波数に合わせて同調する, がある. このように, 周波数を扱うことから, 交流を考える
共振回路 概要 回路は ラジオや通信工学 などに広く使われる この回路の目的は 特定の周波数のときに大きな電流を得ることである 使い方には 周波数を設定し外へ発する 外部からの周波数に合わせて同調する がある このように 周波数を扱うことから 交流を考える 特に ( キャパシタ ) と ( インダクタ ) のそれぞれが 周波数によってインピーダンス *) が変わることが回路解釈の鍵になることに注目する
PowerPoint プレゼンテーション
不飽和土の力学を用いた 締固めメカニズムの解明 締固めとは 土に力を加え 間隙中の空気を追い出すことで土の密度を高めること 不飽和土 圧縮性の減少透水性の減少せん断 変形抵抗の増大 などに効果あり 締固め土は土構造物の材料として用いられている 研究背景 現場締固め管理 締固め必須基準 D 値 施工含水比 施工層厚 水平まきだし ( ρdf ) 盛土の乾燥密度 D値 = 室内締固め試験による最大乾燥密度
Time of Arrival of Shogi Japanese Chessin Japan: Response to Koichi Masukawa s Criticism Yoshinori KIMURA This writing is a response to Koichi Masukawa s criticism of my assertion regarding the time of
Microsoft PowerPoint _分科会PPT資料6-1(JGA+MHI)最終.ppt
超高効率天然ガスエンジン コンバインドシステム技術開発 ( 事後評価 ) 分科会配付資料資料 6-1 超高効率天然ガスエンジン コンバインドシステム技術開発 6. プロジェクトの詳細説明 (1) 超高効率ガスエンジンの開発について (JGA MHI) 平成 20 年 6 月 22 日社団法人日本ガス協会三菱重工業株式会社 1 事業原簿 Ⅲ-7 頁 効率向上への課題 ノッキングを抑制する新燃焼技術の実現
デンソーテクニカルレビュー Vol. 13 No 特集ポート噴射エンジンにおける付着燃料の液膜厚さ解析 * Analysis of the Fuel Liquid Thickness on the Intake Port and Combustion Chamber of a Por
デンソーテクニカルレビュー Vol. 13 No. 1 8 特集ポート噴射エンジンにおける付着燃料の液膜厚さ解析 * Analysis of the Fuel Liquid Thickness on the Intake Port and Combustion Chamber of a Port Fuel Injection Engine 髙橋幸宏 中瀬善博 加藤雄一 Yukihiro TAKAHASHI
2017 (413812)
2017 (413812) Deep Learning ( NN) 2012 Google ASIC(Application Specific Integrated Circuit: IC) 10 ASIC Deep Learning TPU(Tensor Processing Unit) NN 12 20 30 Abstract Multi-layered neural network(nn) has
Microsoft Word - 4_構造特性係数の設定方法に関する検討.doc
第 4 章 構造特性係数の設定方法に関する検討 4. はじめに 平成 年度 年度の時刻歴応答解析を実施した結果 課題として以下の点が指摘 された * ) 脆性壁の評価法の問題 時刻歴応答解析により 初期剛性が高く脆性的な壁については現在の構造特性係数 Ds 評価が危険であることが判明した 脆性壁では.5 倍程度必要保有耐力が大きくなる * ) 併用構造の Ds の設定の問題 異なる荷重変形関係を持つ壁の
自動車技術会論文集-Vol48_No1/本文win-CS6_02.indd
プレス成形時の板厚変化を考慮した振動騒音解析用有限要素モデルの精度向上 * 久保川 明輝 ) 鍵山 恭彦 ) Accuracy Enhancement of FE Model for NVH Analysis by Application of Stamping Effect Haruki Kubokawa Yasuhiko Kagiyama The accuracy and reliability
23_02.dvi
Vol. 2 No. 2 10 21 (Mar. 2009) 1 1 1 Effect of Overconfidencial Investor to Stock Market Behaviour Ryota Inaishi, 1 Fei Zhai 1 and Eisuke Kita 1 Recently, the behavioral finance theory has been interested
周波数特性解析
周波数特性解析 株式会社スマートエナジー研究所 Version 1.0.0, 2018-08-03 目次 1. アナログ / デジタルの周波数特性解析................................... 1 2. 一巡周波数特性 ( 電圧フィードバック )................................... 4 2.1. 部分周波数特性解析..........................................
Introduction to System Identification
y(t) モデルベースデザイン 制御系設計のためのシステム同定入門 s 2 Teja Muppirala t s 2 3s 4 2012 The MathWorks, Inc. 1 モデルベースデザイン 正確なモデルがあることが大前提 実行可能な仕様書 シミュレーションによる設計 モデル 連続したテスト 検証 コード生成による実装 2 動的システムのモデリング モデリング手法 第一原理モデリング データドリブンモデリング
28 Horizontal angle correction using straight line detection in an equirectangular image
28 Horizontal angle correction using straight line detection in an equirectangular image 1170283 2017 3 1 2 i Abstract Horizontal angle correction using straight line detection in an equirectangular image
Microsoft Word doc
. 正規線形モデルのベイズ推定翠川 大竹距離減衰式 (PGA(Midorikawa, S., and Ohtake, Y. (, Attenuation relationships of peak ground acceleration and velocity considering attenuation characteristics for shallow and deeper earthquakes,
電子部品はんだ接合部の熱疲労寿命解析
43 Evaluation for Thermal Fatigue Life of Solder Joints in Electronic Components Haruhiko Yamada, Kazuyoshi Ogawa 2 63Sn- 37Pb 95Pb-5Sn Si Cu Si 63Sn-37Pb Since automotive electronic components are used
エラー動作 スピンドル動作 スピンドルエラーの計測は 通常 複数の軸にあるセンサーによって行われる これらの計測の仕組みを理解するために これらのセンサーの 1つを検討する シングル非接触式センサーは 回転する対象物がセンサー方向またはセンサー反対方向に移動する1 軸上の対象物の変位を測定する 計測
LION PRECISION TechNote LT03-0033 2012 年 8 月 スピンドルの計測 : 回転数および帯域幅 該当機器 : スピンドル回転を測定する静電容量センサーシステム 適用 : 高速回転対象物の回転を計測 概要 : 回転スピンドルは 様々な周波数でエラー動作が発生する これらの周波数は 回転スピード ベアリング構成部品の形状のエラー 外部影響およびその他の要因によって決定される
マツダ技法32_150216
論文 解説 35 サブタンクレス FFV(Flexible Fuel Vehicle) の開発 Development of Sub-Tank Less FFV (Flexible Fuel Vehicle) 古市展之 *1 島本祥 *2 西尾貴史 *3 Nobuyuki Furuichi Sho Shimamoto Takafumi Nishio 渡辺友巳 *4 大橋美貴典 *5 安田京平 *6
Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).
Fig. 1 The scheme of glottal area as a function of time Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels). Fig, 4 Parametric representation
平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損
平成 3 年度 JAXA 航空プログラム公募型研究報告会資料集 (3 年度採用分 1 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損傷発生位置の推定発表内容 (1 荷重同定 1:11 点衝撃荷重同定 ( 荷重同定 : 分布荷重同定 (3 今後の予定
ブローダウン過給システムを用いたガソリンHCCI機関の 運転領域拡大について
1/ 35 ブローダウン過給システムを用いたガソリン HCCI 機関の運転領域拡大 後藤俊介, 窪山達也, 森吉泰生 ( 千葉大学 ) 畑村耕一, 鈴木正剛,( 畑村エンジン研究事務所 ) 山田敏生,(CDAJ) 高梨淳一,( 本田技術研究所 ) 本日の発表内容 2/ 35 1. 背景 2. 目的 3. ブローダウン過給 (BDSCI) システム 4.EGRガイド 5. 実験装置 6. 実験結果と考察
Presentation Title
( ) 2017 The MathWorks, Inc. 1 新型乗用車の CO2 規制の比較 出典 :the International Council on Clean Transportation 2 自動車台数は今後も増加 多くの割合は内燃機関を使用 OEM 各社はコンセプト 想定ターゲットに応じて最適なパワートレインを選択 3 厳しい燃費規制に対応必至 従来の燃焼エンジンからフルバッテリのモータまで
問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた
問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた 測定データを図 1-2 に示す データから, オーステナイト系ステンレス鋼どうしの摩擦係数を推定せよ
Microsoft PowerPoint - 1.プロセス制御の概要.pptx
プロセス制御工学 1. プロセス制御の概要 京都大学 加納学 Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University [email protected] http://www-pse.cheme.kyoto-u.ac.jp/~kano/
MRI | 所報 | 分権経営の進展下におけるグループ・マネジメント
JOURNAL OF MITSUBISHI RESEARCH INSTITUTE No. 35 1999 (03)3277-0003 FAX (03)3277-0520 [email protected] 76 Research Paper Group Management in the Development of Decentralized Management Satoshi Komatsubara,
<4D F736F F F696E74202D A957A A8EC0895E8D7182C982A882AF82E EF89FC915082CC82BD82DF82CC A83808DC5934B89BB A2E >
Techno Forum 2012 実運航における燃費改善のためのトリム最適化 株式会社 MTI 技術戦略グループ上級研究員堀正寿 1 目次 1. はじめに 2. 最適トリムの評価手法 2-1. オペレーションプロファイル調査 2-2. 水槽試験とトリム影響解析 2-3. 実船検証 3. トリムチャートと運用 4. まとめ 2 1-1 トリムの定義 1. はじめに 船尾喫水 (da) と船首喫水 (df)
地球観測衛星データの保存・配布システム
Earth Remote Sensing Data Archive and Distribution System NASAEOS Earth Observing SystemASTER Advanced Spaceborne Thermal Emission and Reflection Radiometer ASTERASTER EOS ASTER DADSData Archive and Distribution
1
問題を解こう. 熱力学の基礎 問題. 容積 [m ] の密閉容器内に 温度 0[ ] 質量 0[kg] の酸素が含まれている この容器内の圧力を求めよ ただし 酸素の気体定数を R= 59.8[J/kg K] とする 解答 酸素の体積 V=m 質量 m=0kg なので 酸素の比容積 v=/0 m /kg である 式 (.) において ガス定数 R=59.8 温度 T=(0+7)K であるので 圧力
加振装置の性能に関する検証方法 Verification Method of Vibratory Apparatus DC-X デジタルカメラの手ぶれ補正効果に関する測定方法および表記方法 ( 光学式 ) 発行 一般社団法人カメラ映像機器工業会 Camera & Imaging Pr
加振装置の性能に関する検証方法 Verification Method of Vibratory Apparatus DC-X011-2012 デジタルカメラの手ぶれ補正効果に関する測定方法および表記方法 ( 光学式 ) 発行 一般社団法人カメラ映像機器工業会 Camera & Imaging Products Association 目 次 1. まえがき ------------------------------------------------------------------------------------------------------
木村の理論化学小ネタ 理想気体と実在気体 A. 標準状態における気体 1mol の体積 標準状態における気体 1mol の体積は気体の種類に関係なく 22.4L のはずである しかし, 実際には, その体積が 22.4L より明らかに小さい
理想気体と実在気体 A. 標準状態における気体 1mol の体積 標準状態における気体 1mol の体積は気体の種類に関係なく.4L のはずである しかし, 実際には, その体積が.4L より明らかに小さい気体も存在する このような気体には, 気体分子に, 分子量が大きい, 極性が大きいなどの特徴がある そのため, 分子間力が大きく, 体積が.4L より小さくなる.4L とみなせる実在気体 H :.449
塗装深み感の要因解析
17 Analysis of Factors for Paint Depth Feeling Takashi Wada, Mikiko Kawasumi, Taka-aki Suzuki ( ) ( ) ( ) The appearance and quality of objects are controlled by paint coatings on the surfaces of the objects.
Microsoft PowerPoint - 資料04 重回帰分析.ppt
04. 重回帰分析 京都大学 加納学 Division of Process Control & Process Sstems Engineering Department of Chemical Engineering, Koto Universit [email protected] http://www-pse.cheme.koto-u.ac.jp/~kano/ Outline
ファイナンスのための数学基礎 第1回 オリエンテーション、ベクトル
時系列分析 変量時系列モデルとその性質 担当 : 長倉大輔 ( ながくらだいすけ 時系列モデル 時系列モデルとは時系列データを生み出すメカニズムとなるものである これは実際には未知である 私たちにできるのは観測された時系列データからその背後にある時系列モデルを推測 推定するだけである 以下ではいくつかの代表的な時系列モデルを考察する 自己回帰モデル (Auoregressive Model もっとも頻繁に使われる時系列モデルは自己回帰モデル
京都大学博士 ( 工学 ) 氏名宮口克一 論文題目 塩素固定化材を用いた断面修復材と犠牲陽極材を併用した断面修復工法の鉄筋防食性能に関する研究 ( 論文内容の要旨 ) 本論文は, 塩害を受けたコンクリート構造物の対策として一般的な対策のひとつである, 断面修復工法を検討の対象とし, その耐久性をより
塩素固定化材を用いた断面修復材と犠牲陽極材を併用し Titleた断面修復工法の鉄筋防食性能に関する研究 ( Abstract_ 要旨 ) Author(s) 宮口, 克一 Citation Kyoto University ( 京都大学 ) Issue Date 2015-01-23 URL https://doi.org/10.14989/doctor.k18 Right Type Thesis
PowerPoint プレゼンテーション
音響解析プログラム WAON 最新開発動向と適用例のご紹介 サイバネットシステム株式会社 メカニカル CAE 事業部 WAON 推進室 アジェンダ 1. 会社紹介 2. WAON とは? 3. なぜ WAON なのか? 4. 各種適用例のご紹介 5. 最新開発動向 2 1. 会社紹介サイバネットシステム ( 株 ) メカニカル CAE 事業部 音響 構造 熱 電磁場 熱流体 衝突 板成形 樹脂流動などの各種解析
MKT-TAISEI Co.,Ltd.
MKT-TAISEI Co.,Ltd. CONTENTS はやぶさ と共に燃え尽きた アンテナ アンテナ 2.45GHz 4 パッチアンテナ UHF 帯アンテナ 2.4GHz フラットアンテナ 22.5GHz/24GHz/26.5GHz マイクロストリップアレーアンテナ 3 4 5 antenna_contents-04 アンテナ ANTENNA 製品例 2.45GHz 4 Patch Antenna
Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Enginee
Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Engineering, Kanazawa University, Kakuma-machi, Kanazawa-shi,
PowerPoint プレゼンテーション
1 2017 年度ベンチマーク試験について エネルギーフロー車両試験 (2 台 ) (1) Chevrolet Malibu (2) BMW 320i エンジンユニット単体試験 (2 ユニット ) (3) Mercedes C200 (M274) (4) Audi A4 (EA888) エンジンフリクション測定 (2 ユニット ) (5) Chevrolet Malibu(Ecotec 1.5L )
Synthesis and Development of Electric Active Stabilizer Suspension System Shuuichi BUMA*6, Yasuhiro OOKUMA, Akiya TANEDA, Katsumi SUZUKI, Jae-Sung CHO
Synthesis and Development of Electric Active Stabilizer Suspension System Shuuichi BUMA*6, Yasuhiro OOKUMA, Akiya TANEDA, Katsumi SUZUKI, Jae-Sung CHO and Masaru KOBAYASHI Chassis Engineering Management
Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷
熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI プロジェクト @ 宮崎県美郷町 熊本大学副島慶人川村諒 1 実験の目的 従来 信号の受信電波強度 (RSSI:RecevedSgnal StrengthIndcator) により 対象の位置を推定する手法として 無線 LAN の AP(AccessPont) から受信する信号の減衰量をもとに位置を推定する手法が多く検討されている
Microsoft PowerPoint - H22制御工学I-10回.ppt
制御工学 I 第 回 安定性 ラウス, フルビッツの安定判別 平成 年 6 月 日 /6/ 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し
. はじめに 本文書は, ジャイロもしくは周波数発振器などの性能評価にしばしば用いられるアラン分散 について記したものである.. 目的 本文書は, アラン分散の概念及び計算方法, そして評価方法について述べ, アラン分散を 用いた解析のノウハウを習得することを目的とする 3. 参考文書 参考文書を以
目次. はじめに.... 目的... 3. 参考文書... 4. アラン分散とは... 3 5. アラン分散の定義... 3 6. アラン分散の計算方法... 4 7. アラン分散計算プログラム... 6 8. グラフの読み取り方... 9 9. アラン分散関係の参考資料...0 . はじめに 本文書は, ジャイロもしくは周波数発振器などの性能評価にしばしば用いられるアラン分散 について記したものである..
03マツダ技報29本文-h1-01-1案0620
新 型 デミオのエンジン 技 術 No.29(2011) 特 集 : 新 型 マツダデミオ 2 新 型 デミオのエンジン 技 術 Engine Technology for New Demio *1 富 澤 和 廣 *2 松 尾 佳 朋 大 槻 Kazuhiro Tomizawa Yoshitomo Matsuo Ken Ohtsuki *4 室 谷 満 幸 後 藤 *5 剛 *3 健 *6 上 月
第 4 週コンボリューションその 2, 正弦波による分解 教科書 p. 16~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問 1. 以下の図にならって,1 と 2 の δ 関数を図示せよ δ (t) 2
第 4 週コンボリューションその, 正弦波による分解 教科書 p. 6~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問. 以下の図にならって, と の δ 関数を図示せよ. - - - δ () δ ( ) - - - 図 δ 関数の図示の例 δ ( ) δ ( ) δ ( ) δ ( ) δ ( ) - - - - - - - -
音響部品アクセサリ本文(AC06)PDF (Page 16)
Guide for Electret Condenser Microphones A microphone as an audio-electric converting device, whose audio pickup section has a structure of a condenser consisting of a diaphragm and a back plate opposite
Vol. 36, Special Issue, S 3 S 18 (2015) PK Phase I Introduction to Pharmacokinetic Analysis Focus on Phase I Study 1 2 Kazuro Ikawa 1 and Jun Tanaka 2
Vol. 36, Special Issue, S 3 S 18 (2015) PK Phase I Introduction to Pharmacokinetic Analysis Focus on Phase I Study 1 2 Kazuro Ikawa 1 and Jun Tanaka 2 1 2 1 Department of Clinical Pharmacotherapy, Hiroshima
Microsoft PowerPoint - 【最終提出版】 MATLAB_EXPO2014講演資料_ルネサス菅原.pptx
MATLAB/Simulink を使用したモータ制御アプリのモデルベース開発事例 ルネサスエレクトロニクス株式会社 第二ソリューション事業本部産業第一事業部家電ソリューション部 Rev. 1.00 2014 Renesas Electronics Corporation. All rights reserved. IAAS-AA-14-0202-1 目次 1. はじめに 1.1 モデルベース開発とは?
[2] 1. 2. 2 2. 1, [3] 2. 2 [4] 2. 3 BABOK BABOK(Business Analysis Body of Knowledge) BABOK IIBA(International Institute of Business Analysis) BABOK 7
32 (2015 ) [2] Projects of the short term increase at present. In order to let projects complete without rework and delays, it is important that request for proposals (RFP) are written by reflecting precisely
アクティブフィルタ テスト容易化設計
発振を利用したアナログフィルタの テスト 調整 群馬大学工学部電気電子工学科高橋洋介林海軍小林春夫小室貴紀高井伸和 発表内容. 研究背景と目的. 提案回路 3. 題材に利用したアクティブフィルタ 4. 提案する発振によるテスト方法 AG( 自動利得制御 ) バンドパス出力の帰還による発振 3ローパス出力の帰還による発振 4ハイパス出力の帰還による発振. 結果 6. まとめ 発表内容. 研究背景と目的.
Time Variation of Earthquake Volume and Energy-Density with Special Reference to Tohnankai and Mikawa Earthquake Akira IKAMi and Kumizi IIDA Departmen
Time Variation of Earthquake Volume and Energy-Density with Special Reference to Tohnankai and Mikawa Earthquake Akira IKAMi and Kumizi IIDA Department of Earth Sciences, Nagoya University (Received January
