(57)【特許請求の範囲】
【請求項１】
断熱容器の種類、寸法、初期温度及び配置と、該断熱容器の周囲の外気温度と、該断熱容器に収納する内容物の種類、重量、初期温度及び配置と、該内容物の温度保持のための熱媒体の種類、重量及び配置とに関するデータを入力する入力手段と、前記入力手段により入力されたデータに基づいて対応する前記内容物の温度保持性能に関するデータを演算し出力する制御手段と、を備えることを特徴とする断熱容器の設計システム。
【請求項２】
前記内容物の温度保持性能に関するデータは、前記内容物の鮮度劣化に特有な係数をA、前記内容物の活性化エネルギーをE、気体定数をR、絶対温度をT、生産時の鮮度をF(0)として、
F(t)=F(0)EXP(-∫A・EXP(-E/R/T)dt)
によって得られる、所定時間経過後の内容物の鮮度F(t)であることを特徴とする請求項１記載の断熱容器の設計システム。
【請求項３】
断熱容器の種類、寸法、初期温度及び配置と、該断熱容器の周囲の外気温度と、該断熱容器に収納する内容物の種類、重量、初期温度及び配置と、該内容物の温度保持のための熱媒体の種類、重量及び配置とに関するデータを受信する受信手段と、前記受信手段により入力されたデータに基づいて対応する前記内容物の温度保持性能に関するデータを演算し送信する送信手段と、を備えることを特徴とする断熱容器の設計システム。
【請求項４】
コンピュータを、断熱容器の種類、寸法、初期温度及び配置と、該断熱容器の周囲の外気温度と、該断熱容器に収納する内容物の種類、重量、初期温度及び配置と、該内容物の温度保持のための熱媒体の種類、重量及び配置とに関するデータを入力する入力手段と、前記入力手段により入力されたデータに基づいて対応する前記内容物の温度保持性能に関するデータを演算し出力する制御手段と、を備えることを特徴とする断熱容器の設計システムとして機能させるためのプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
【請求項５】
コンピュータを、断熱容器の種類、寸法、初期温度及び配置と、該断熱容器の周囲の外気温度と、該断熱容器に収納する内容物の種類、重量、初期温度及び配置と、該内容物の温度保持のための熱媒体の種類、重量及び配置とに関するデータを受信する受信手段と、前記受信手段により入力されたデータに基づいて対応する前記内容物の温度保持性能に関するデータを演算し送信する送信手段と、を備えることを特徴とする断熱容器の設計システムとして機能させるためのプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は断熱容器の設計システムに関し、特に、ユーザニーズに応じた断熱容器の温度保持性能に関するデータを提供し、容器製造に係る開発期間を大幅に短縮することができる断熱容器の設計システム、断熱容器の設計プログラムを記録した記録媒体に関する。
【０００２】
【従来の技術】
食品、医薬品等の物流においては、内容物を輸送する過程で十分な鮮度が確保されることが求められる。
断熱容器の設計は、容器の断熱性能の他、内容物、輸送条件により容器内部の状態が大きく変動するために試作品による実用化レベルでの評価が欠かせず、多くの作業工数と費用が発生する。このため、断熱容器は、短時間での製品化が非常に困難な商品である。
【０００３】
【発明が解決しようとする課題】
従来の断熱容器は、過去の製作経験に基づいて設計されている。この結果、安全性を過剰に含めた設計となり、保冷容器の資源、エネルギー面での無駄が多い。
本発明の目的は、ユーザニーズに応じた断熱容器の温度保持性能を提供することができ、容器製造に係る開発期間を大幅に短縮することができる断熱容器の設計システム、及び断熱容器の設計プログラムを記録した記録媒体を提供することにある。
【０００４】
【課題を解決するための手段】
本発明の断熱容器の設計システムは、断熱容器の性能、該断熱容器の周囲の外気温度、該断熱容器に収納する内容物の熱容量及び初期温度、並びに該内容物の温度保持のための熱媒体の性能に関するデータを入力する入力手段と、該入力手段により入力されたデータに対応する前記内容物の温度保持性能に関するデータを出力する制御手段と、を備えるものである。
【０００５】
また、前記制御手段は、前記入力手段により入力されたデータに基づいて、温度保持性能を演算することで、入力条件に応じたシミュレーション結果をユーザに提供することができる。
また、前記断熱容器の性能に関するデータには、断熱容器の種類、寸法及び温度が含まれることで、断熱容器に関するより詳細な入力条件に応じたシミュレーション結果をユーザに提供することができる。
【０００６】
また、前記内容物の熱容量に関するデータには、内容物の種類及び重量が含まれることで、内容物に関するより詳細な入力条件に応じたシミュレーション結果をユーザに提供することができる。
また、前記熱媒体の性能に関するデータには、熱媒体の種類及び重量が含まれることで、保冷材又は保温材に関するより詳細な入力条件に応じたシミュレーション結果をユーザに提供することができる。
【０００７】
また、前記内容物の温度保持性能に関するデータは、所定時間経過後の内容物の温度であることで、ユーザは表示された温度分布図が示す色又は諧調に基づいて内容物の温度変化を知ることができる。
また、前記内容物の温度保持性能に関するデータは、所定時間経過後の内容物の鮮度であることで、ユーザは計算結果に基づいて内容物の鮮度を知ることができる。
【０００８】
他の観点において、本発明の断熱容器の設計システムは、断熱容器の性能、該断熱容器の周囲の外気温度、該断熱容器に収納する内容物の熱容量及び初期温度、並びに該内容物の温度保持のための熱媒体の性能に関するデータを受信する受信手段と、該受信手段により入力されたデータに対応する前記内容物の温度保持性能に関するデータを送信する送信手段と、を備えるものである。
【０００９】
他の観点において本発明は、コンピュータを、断熱容器の性能、該断熱容器の周囲の外気温度、該断熱容器に収納する内容物の熱容量及び初期温度、並びに該内容物の温度保持のための熱媒体の性能に関するデータを入力する入力手段と、該入力手段により入力されたデータに対応する前記内容物の温度保持性能に関するデータを出力する制御手段と、を備える断熱容器の設計システムとして機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体である。
【００１０】
他の観点において本発明は、コンピュータを、断熱容器の性能、該断熱容器の周囲の外気温度、該断熱容器に収納する内容物の熱容量及び初期温度、並びに該内容物の温度保持のための熱媒体の性能に関するデータを受信する受信手段と、該受信手段により入力されたデータに対応する前記内容物の温度保持性能に関するデータを送信する送信手段と、を備える断熱容器の設計システムとして機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体である。
【００１１】
本発明は、ユーザから予め要求される内容物の必要保持温度及び時間等が入力された場合、コンピュータシミュレーション又は過去の温度保持性能データを利用して、たとえば発泡プラスチックやダンボール素材からなる箱形状容器を設計提案するものである。これにより、必要に応じて試作に必要な素材や試作品を提供することにも事業発展することができ、断熱容器を必要とするユーザに代わり開発支援することができる。
【００１２】
【発明の実施の形態】
以下、本発明の実施の形態を添付図面と対応して詳細に説明する。
図１は、本発明の実施の形態による断熱容器の設計システムの構成を説明するブロック図である。本システムは、クライアント端末１、インターネット３及びサーバ５により実現される。
【００１３】
図１において、クライアント端末１、サーバ５の各々とインターネット３との間には、インターネットサービスプロバイダ（図示せず）が配置されている。このインターネットサービスプロバイダとクライアント端末１及びサーバ５の各々とは、ＩＳＤＮ回線等により接続されている。クライアント端末１は、ＴＣＰ／ＩＰ等のインターネット用の通信プロトコルを用いて、サーバ５と相互に通信可能なように接続される。
【００１４】
クライアント端末１は、断熱容器の設計支援を求めるユーザによって操作される。本実施の形態では、食品会社の物流、包装分野の技術開発者等により操作され、断熱容器の性能（種類、寸法及び温度）、該断熱容器の周囲の外気温度、該断熱容器に収納する内容物の熱容量に関するデータ（種類及び重量）及び初期温度、並びに該内容物の温度保持のための熱媒体（保冷材又は保温材）の性能（種類、重量及び温度）が入力される。
【００１５】
サーバ５は、インターネット３を介してホームページ等を窓口としてクライアント端末１に断熱容器の温度保持性能に関するデータを提供する。ここで、温度保持性能データを定義する。
クライアント端末１に提供される温度保持性能データとは、主に温度シミュレーション、鮮度シミュレーションにより得られるデータである。
【００１６】
具体的には、解析条件入力画面（図１５参照）のボタン“解析実行”をクリックして表示される温度分布（図１９参照）、及び食材鮮度パラメータ計算条件画面（図１２参照）のボタン“計算”をクリックして表示される鮮度Ｆである。
一方、サーバ５には、上記温度保持性能データと、シミュレーション画面（図６〜図１１、図１３及び図１５参照）でユーザが条件入力したデータ及び各種ボタンをクリックして表示されるデータとが、解析モデル（図１５参照）毎に格納される。ここで、サーバ５の詳細構成を説明する。
【００１７】
サーバ５は、制御部７、送信／受信部９、ＲＡＭ１１、ＲＯＭ１３、温度保持性能データベース（性能ＤＢ）１５、パラメータデータベース（パラメータＤＢ）１６、入力操作部１７、表示部１９、ＣＤ−ＲＯＭドライブ２１から構成される。
制御部７はＣＰＵであり、サーバ５全体の動作を制御する。特に、送信／受信部９の通信制御、ＲＡＭ１１、ＲＯＭ１３、性能ＤＢ１５及びパラメータＤＢ１６等からのデータ読出し／書き込み処理、後述する温度及び鮮度シミュレーションにおける温度保持性能データの演算処理等である。
【００１８】
送信／受信部９は、制御部７の指示に基づいて、クライアント端末１との間でデータの送信及び受信処理を行う。特に、送信／受信部９は、外気温度及び時間、断熱容器の種類、寸法及び温度、内容物の種類、重量及び温度、並びに熱媒体（保冷材又は保温材）の種類、重量及び温度をクライアント端末１から受信し、また、制御部７で演算された温度保持性能データをクライアント端末１に送信する。
【００１９】
ＲＡＭ１１は、サーバ５の処理に必要なデータを一時的に格納する。特に、ＲＡＭ１１は、クライアント端末１にて表示されるユーザインタフェースとなる画面データ等を格納する。ＲＯＭ１３は、本発明の処理プログラム等を固定的に格納する。
性能ＤＢ１５は、光磁気ディスクライブラリ等の大容量の記憶装置であり、断熱容器の過去のシミュレーション事例の温度保持性能データを解析モデル毎に格納する。
【００２０】
パラメータＤＢ１６は、光磁気ディスクライブラリ等の大容量の記憶装置であり、後述する温度伝導率等、シミュレーションに必要なパラメータを記憶する。入力操作部１７は、後述する性能ＤＢ１５、パラメータＤＢ１６の記憶内容の更新に応じてその都度操作される。
【００２１】
表示部１９はＬＣＤ等であり、制御部７で演算された温度保持性能データ（計算結果、グラフ）等のコードデータをその都度表示データに変換して表示処理を行う。
ＣＤ−ＲＯＭドライブ２１は、制御部７の指示に基づいてＣＤ−ＲＯＭ２３に格納されているプログラム等をＲＡＭ１１等に書き込む。
【００２２】
図２は、本実施の形態による断熱容器の設計システムの動作を説明するフローチャートである。ここでは、（１）クライアント端末１における表示画面及び該画面への入力操作（図３〜図１８）、（２）サーバ５の制御部７での温度保持性能の演算（温度・鮮度シミュレーション）を中心に説明する。
【００２３】
はじめに、ステップＳ１０１では、クライアント端末１は、インターネット３を介してホームページ等の情報提供の窓口と通信回線を確立し、次いでステップＳ１０２で、初期画面（図３参照）が表示される。
ステップＳ１０３では、この初期画面のボタン“開始”をクリックすることで、容器数入力画面（図４参照）が表示される。本実施の形態における輸送例として、被保冷食品としてケーキ（約３００ｇ）をダンボール（小箱）に包装し、パン（約３００ｇ）をダンボール（小箱）に包装し、これら２つの容器を１つのダンボール（物流容器）に収容して輸送する場合を説明する。そこで、容器数“３”を入力する。
【００２４】
ステップＳ１０３での容器数“３”の入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０４で、容器構成指定画面（図５参照）が表示される。この画面は、先に入力した容器数“３”の構成を指定する画面であり、図５では、容器を３段重ね（３重容器）に構成するか、２段重ね（２重容器）に構成するかを指定する。ここで、“容器１⊃容器２”は、容器２が容器１に収容されることを意味している。上述したことから、容器２、容器３が、ケーキとパンとを各々包装して容器１に収容される構成とするため、ここでは容器２と容器３とが同格であることを意味する記号“＝”を入力する。
【００２５】
ステップＳ１０４での記号“＝”の入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０５で、シミュレーション画面（外気条件入力画面、図６参照）に移る。
シミュレーション画面は、外気条件、容器条件、食材条件、冷媒条件、解析条件の５種類の条件入力画面からなる。図６は、外気条件を設定する画面であり、食品の輸送過程で想定される外気温度と該温度の元での輸送時間とが設定される。本実施の形態では、全体として９時間の輸送時間での温度設定に関して、外気温度２５℃での輸送が３時間、外気温度３０℃での輸送が３時間、外気温度３５℃での輸送が３時間に設定される。
【００２６】
なお、サーバ５のパラメータＤＢ１６には、パラメータ「空気の温度伝導率」、「動粘性係数v」、「重力加速度ｇ」、「体積膨張率β」、「空気の熱伝導率λ」、「プラントル数Pr」が予め格納されている。制御部７は、上記外気に関する条件入力に応答して該当するパラメータ値をパラメータＤＢ１６から読み出し、この画面に表示する。
【００２７】
また、ここで表示されている「空気の温度伝導率」は、後述する温度シミュレーションの“非定常二次元熱伝達式”で用いられ、「空気の温度伝導率」を除くパラメータ値は、自然対流を考慮する場合に、後述する温度シミュレーションの“自然対流熱伝達式”で用いられる（図１５参照）。
【００２８】
ステップＳ１０５での外気条件入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０６で、容器１の条件入力画面（図７参照）が表示される。この容器条件入力画面は、図５での容器構成指定にリンクしており、容器No.“容器１”が表示されている。
【００２９】
ここで、上述した輸送例に合わせて入力を行う。具体的には、分類への入力にはボタン“検索”をクリックし、“発泡容器”でボタン“登録”をクリックして入力する。同様に、容器名への入力にはボタン“検索”をクリックし、“物流容器”でボタン“登録”をクリックして入力する。
【００３０】
また、容器特性として横（Ｘ軸）には“３５（ｃｍ）”、縦（Ｚ軸）には“３０（ｃｍ）”及び厚さには“２．５（ｃｍ）”を入力する。さらに、食材数“２”、冷媒数“１”、容器１内の初期温度“１５（℃）”、内気温度“１５（℃）”を入力する。容器の種類（発泡容器、ダンボール容器等）、寸法（外寸、内寸、厚み）は任意に設定できる。
【００３１】
図７に示す容器１の条件入力後、ボタン“ＯＫ”をクリックすると容器２の条件入力画面（図８参照）が表示される。この容器条件入力画面は、図５での指定内容にリンクしており、容器No.“容器２”が表示されている。
ここでの入力内容は、具体的には、分類への入力にはボタン“検索”をクリックし、“ダンボール”でボタン“登録”をクリックして入力する。同様に、容器名への入力にはボタン“検索”をクリックし、“小箱１”でボタン“登録”をクリックして入力する。また、図７での入力内容に加えて、容器１における容器２の位置として、所定の位置からのＸ座標“５（ｃｍ）”、Ｚ座標“２．５（ｃｍ）”を入力する。
【００３２】
図８に示す容器２の条件入力後、ボタン“ＯＫ”をクリックすると容器３の条件入力画面（図９参照）が表示される。この容器条件入力画面は、図５での指定内容にリンクしており、容器No.“容器３”が表示されている。
ここでの入力内容は、具体的には、分類への入力にはボタン“検索”をクリックし、“ダンボール”でボタン“登録”をクリックして入力する。同様に、容器名への入力にはボタン“検索”をクリックし、“小箱２”でボタン“登録”をクリックして入力する。また、図８での入力内容と同様に、容器１における容器３の位置として、容器２と共通の所定の位置からのＸ座標“２０（ｃｍ）”、Ｚ座標“２．５（ｃｍ）”を入力する。
【００３３】
なお、サーバ５のパラメータＤＢ１６には、容器の種類に対応してパラメータ「容器温度伝導率」、「熱伝導率」、「内壁表面の境膜係数」及び「外壁表面の境膜係数」が予め格納されている。新しいものはその都度追加されていく。制御部７は、上記容器１、２、３に関する条件入力に応答して該当するパラメータ値をパラメータＤＢ１６から読み出し、図７、図８、図９の画面に表示する。
【００３４】
また、図７、図８、図９で各々表示されている「容器温度伝導率」は、後述する温度シミュレーションの“非定常二次元熱伝達式”で用いられ、「熱伝導率」、「内壁表面の境膜係数」、「外壁表面の境膜係数」及び「（容器１、２、３の）厚さ」は、後述する温度シミュレーションの“非定常二次元熱伝達式”の前提となる“冷媒の融解時間及び融解速度式”で用いられる。
【００３５】
以上により、容器１〜容器３の条件入力が完了し、図９に示す容器３の条件入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０７で、食材１の条件入力画面（図１０参照）が表示される。この容器条件入力画面は、図５での容器構成指定及び図８での容器２の条件入力内容にリンクしており、容器名“小箱１”が表示されている。
【００３６】
上述した輸送例及び入力操作に従い、ボタン“検索”及び“登録”をクリックして食材No.“食材１”、分類“加工品”、食材名“ケーキ”を入力し、該ケーキの内容（ここでは重量）“３００ｇ”、及び初期温度“５（℃）”を入力する。
つぎに、ケーキの形状として、ここでは“直方型”をクリックし、横ａ“５（ｃｍ）”、高さｂ“５（ｃｍ）”を入力する。つぎに、容器２における設置位置として、所定位置からの距離として横（Ｘ軸）“２．５（ｃｍ）”、高さ（Ｚ軸）“０．５（ｃｍ）”を入力する。
【００３７】
図１０に示す食材１の条件入力後、ボタン“ＯＫ”をクリックすると食材２の条件入力画面（図１１参照）が表示される。この容器条件入力画面は、図５での容器構成指定及び図９での容器３の条件入力内容に連動し、容器名“小箱２”が表示されている。
上述した輸送例及び入力操作に従い、ボタン“検索”及び“登録”をクリックして食材No.“食材２”、分類“加工品”、食材名“パン”を入力し、該パンの内容（ここでは重量）“３００ｇ”、及び初期温度“５（℃）”を入力する。
【００３８】
図１０の入力内容と同様に、パンの形状として、ここでは“楕円型”をクリックし、短半径ａ“２．５（ｃｍ）”、長半径ｂ“３（ｃｍ）”を入力する。つぎに、容器３における設置位置として、所定位置からの距離として横（Ｘ軸）“２（ｃｍ）”、高さ（Ｚ軸）“０．５（ｃｍ）”を入力する。
【００３９】
なお、サーバ５のパラメータＤＢ１６には、図１０、図１１に示すパラメータ「（食材）温度伝導率」が食品の種類に対応して予め格納されている。新しいものはその都度追加されていく。制御部７は、上記食材１、２に関する条件入力に応答して該当するパラメータ値をパラメータＤＢ１６から読み出し、図１０、図１１の画面に表示する。
また、図１０、図１１で各々表示されている「（食材）温度伝導率」は、後述する温度シミュレーションの“非定常二次元熱伝達式”で用いられる。
【００４０】
ここで、鮮度シミュレーションのためのパラメータの決定方法について説明する。図１０又は図１１で、「鮮度予測値」のボタン“計算”をクリックすると食材鮮度パラメータ計算条件画面（図１２参照）が各々表示される。この画面は、食材鮮度パラメータＡ、Ｅを計算する画面である。「鮮度劣化に特有な定数Ａ」及び「活性化エネルギーＥ」は食品に固有な値であり、温度と保存日数から求められる。
【００４１】
図１２では、特定の食材の温度２点での実測データが入力されている。温度“３０（℃）”で保存する場合、初期の鮮度を“１００”とすると保存日数“１（日）”経過後の鮮度は“４０”に低下し、温度“５（℃）”で保存する場合、初期の鮮度を“１００”とすると保存日数“７（日）”経過後の鮮度は“４０”に低下したことを示すデータが入力されている。そこで、ボタン“計算”をクリックすると、上記食材の入力条件に基づいて「鮮度に特有な定数Ａ」及び「活性化エネルギーＥ」が各々計算され、図１０（又は図１１）及び図１２に表示される。
【００４２】
また、ここで表示される「鮮度に特有な定数Ａ」及び「活性化エネルギーＥ」は、後述する鮮度シミュレーションの“鮮度予測式”で用いられる。
以上により、食材１〜食材２の条件入力が完了し、図１１に示す食材２の条件入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０８で、冷媒条件入力画面（図１３参照）が表示される。本実施の形態では、熱媒体として保冷材（氷、ドライアイス等の蓄冷材）を用いる場合を説明する。この冷媒条件入力画面は、図５での容器構成指定及び図７での容器１の条件入力内容（冷媒数“１”）にリンクしており、容器名“物流容器”が表示されている。
【００４３】
上述した輸送例及び入力操作に従い、ボタン“検索”及び“登録”をクリックして冷媒No.“冷媒１”、分類“蓄冷材”、冷媒名“ケーキ用”を入力し、該蓄冷材の内容（ここでは重量）“１．４ｋｇ”、及び初期温度“−１５（℃）”を入力する。冷媒の融解温度は、任意に指定できる（蓄冷タイプ０℃以下〜蓄熱タイプ０℃以上）。
【００４４】
また、蓄冷材の形状として、ここでは“直方型”をクリックし、高さａ“２（ｃｍ）”、横ｂ“２５（ｃｍ）”を入力する。つぎに、容器３における設置位置として、所定位置からの横（Ｘ軸）“２（ｃｍ）”、高さ（Ｚ軸）“０．５（ｃｍ）”を入力する。
【００４５】
なお、サーバ５のパラメータＤＢ１６には、図１３に示すパラメータ「（冷媒の）温度伝導率」及び「融解潜熱」が予め格納されている。新しいのものはその都度追加されていく。制御部７は、上記冷媒に関する条件入力に応答して該当するパラメータ値をパラメータＤＢ１６から読み出し、この画面に表示する。
【００４６】
また、ここで表示されている「（食材）温度伝導率」は、後述する温度シミュレーションの“非定常二次元熱伝達式”で用いられる。
【００４７】
以上により、容器１〜容器３、食材１〜食材２及び冷媒１の条件入力が完了する。ここで、図７〜図１１、図１３及び後述する図１５におけるボタン“配置図表示”をクリックすると条件入力レベルでの配置図が表示される。たとえば、図１３での冷媒１の条件入力を経てボタン“配置図表示”をクリックすると図１４のような配置図が表示される。容器２，３、食材１，２、及び冷媒１の容器内の配置は、任意に設定できる。また、２重容器、３重容器において、食材、冷媒及び容器の配置を任意に設定でき、図１４の配置図を表示することで全体の配置を容易に把握することができる。
【００４８】
図１３での蓄冷材の条件入力後、ボタン“ＯＫ”をクリックすると、次いでステップＳ１０９で、解析条件入力画面（図１５参照）が表示される。この画面には、外気条件、容器条件、食材条件、及び冷媒条件の内容に基づいてシミュレーションを行うための条件を入力する。
【００４９】
ここで、解析名には“解析モデル１”、解析の時間ステップ幅ｎには“１８０（sec）”、及び温度経過を観測する時間である計算実時間には“９（H）”を入力する。また、自然対流項として、ここでは“考慮しない”をクリックする。ちなみに、“考慮する”をクリックしてボタン“補正定数”をクリックすると、自然対流補正入力画面（図１６参照）が表示される。
【００５０】
この自然対流補正定数は、各容器（発泡容器、段ボール容器等）、各食材、各冷媒とも０．４２にしているが、実測値に合わせて補正する場合、０．３８〜０．４６の範囲で補正定数を各々変更しその数値で解析を行う（プラントル数Pr＝０．７〜１０）。上述のように、この各補正定数と図６で表示されている各パラメータを用いて、自然対流を考慮した温度シミュレーションが実行される。
なお、サーバ５のパラメータＤＢ１６には、各容器、各食材、各冷媒に対応して各補正定数が予め格納されている。新しいものは補正の都度追加されていく。
【００５１】
また、図１５において、ボタン“メッシュ分割指定”をクリックすると、メッシュ分割指定画面（図１７参照）が表示される。ここでは、解析モデル配置図に対して、Ｘ座標、Ｚ座標に関してスクロールバーをクリックすることによりメッシュ分割する。
【００５２】
このメッシュ分割による各領域は、後述する温度シミュレーションの実行単位となるものである。図１７に示す例では、食材１、２の高さ及び幅、容器２、３の高さ及び幅等に合わせて、メッシュ分割されている。メッシュ分割を不等間隔にすることにより、計算時間を早くすることができる。高い精度が要求される食材の領域は、細かいメッシュ分割を行うことが望ましい。
【００５３】
また、図１５において、ボタン“融解表示”をクリックすると、融解情報画面（冷媒−融解計算一覧、図１８参照）が表示される。ここでは、温度シミュレーション（“非定常二次元熱伝達式”、更には“自然対流熱伝達式”）を実行するに際して、図１３での入力条件（重量、融解潜熱）に基づいて、冷媒が融解する時間を計算により予測する。
【００５４】
ステップＳ１０９での解析条件入力後、図１７又は図１８におけるボタン“解析実行”をクリックすると、次いでステップＳ１１０で、ステップＳ１０１〜１０９での条件に基づいて、温度・鮮度シミュレーション（温度保持性能の解析）が実行される。このシミュレーションによりクライアント端末１には温度分布図（図１９参照）が表示される。ここでは、解析モデル配置図での温度変化がカラーマップ表示される。
ステップＳ１１０でのシミュレーションは、（１）温度シミュレーションと、（２）鮮度シミュレーションとに分けられる。これらについて以下に説明する。
【００５５】
はじめに、温度シミュレーションについて説明する。温度シミュレーション式は、任意の外気温度に対して、任意の経過時間における、容器内の温度（内容物、空気、等）を算出する計算式であり、（温度予測式）＝（非定常二次元熱伝導式）+（自然対流熱伝達式）で表される。以下、各式について説明する。
【００５６】
▲１▼非定常二次元熱伝導式
（基礎式）δT/δｔ＝κ・（δ2T/δｘ2＋δ2T/δｚ2）
ここで、Tは絶対温度、ｔは時間、κは温度伝導率、ｘ,ｚは座標を示す。
【００５７】
上記基礎式を数値計算するために差分化する。基礎式の時間の項に前進差分、空間の項に中心差分すると、基礎式は以下の様に差分化される。
【００５８】
Tijn+1＝Δt・κ{（Ti+1jn−2Tijn＋Ti-1jn）/(Δx)2＋（Tij+1n−2Tijn＋Tij-1n）/(Δz)2}＋Tijn
ここで、ｎは時間ステップ幅、i,jは、座標z,x方向の計算点を示す。
境界条件として、異物質が接している境界上では調和平均の温度伝導率を用いる。
【００５９】
κ-1＝（1/Ｎ）・Σκｍ-1
ここで、Ｎは境界で接している物質の数、κｍは各物質の温度伝導率を示す。
【００６０】
上述のように、この熱伝導式に基づくシミュレーションは、図１７でメッシュ分割された領域毎に実行される。このため、ここでの温度伝導率κは、図６〜図１１及び図１３で表示されている温度伝導率のうち、シミュレーション実行対象となる領域に含まれるものの温度伝導率が用いられる。
【００６１】
たとえば、図１７に示す食材“ケーキ（直方型にて図示）”が含まれるメッシュ領域に対するシミュレーションでは、ケーキの「温度伝導率」（図１０参照）が用いられ、ケーキと空気の境界上では、上記境界条件に基いて調和平均の温度伝導率が用いられる。
【００６２】
▲２▼自然対流熱伝達式
この熱伝達式に関しても、自然対流を考慮する場合には、メッシュ領域毎に実行される。物体表面から流体への対流伝熱量は流れ方向をX軸、表面に垂直な方向をY軸にとると、フーリエの法則により以下の式で与えられる。
【００６３】
【数１】
0003835980.tif 000002
ここで、
【００６４】
【数２】
0003835980.tif 000003
は熱流束、λは流体の熱伝導率である。自然対流がある場合は、熱伝導と同時に次式で与えられる熱伝達がある。
【００６５】
【数３】
0003835980.tif 000004
ここで、Tw，T∞は壁と遠方流体の温度、hは熱伝達率(W/(m2K))を示す。温度がゆっくり変化する自然対流の温度場では、以下のエネルギー方程式が成り立つ。
【００６６】
【数４】
0003835980.tif 000005
ここで、κは温度伝導率(m2/s)を示す。
【００６７】
いま、鉛直平板自然対流境界層熱伝達において、実験によりy=0でu=0，y=δでdu/dy=0の速度境界条件のもとで速度分布形
【数５】
0003835980.tif 000006
及び、y=0でT=Tw，y=δtでT=T∞，dT/dy=0の温度境界条件の下で温度分布形
【００６８】
【数６】
0003835980.tif 000007
を仮定できるものとすると、鉛直平板の層流自然対流熱伝達についてのEdeの式が導かれる。
【００６９】
【数７】
0003835980.tif 000008
ここで、
【００７０】
【数８】
0003835980.tif 000009
は局所ヌッセルト数、
【００７１】
【数９】
0003835980.tif 000010
はプラントル数、
【００７２】
【数１０】
0003835980.tif 000011
はグラスホフ数を示す。また、Pr=0.7〜10に対して、
【００７３】
【数１１】
0003835980.tif 000012
と近似できる。したがって、
【００７４】
【数１２】
0003835980.tif 000013
となる。ここで、平均熱伝達率は
【００７５】
【数１３】
0003835980.tif 000014
である。ここで、λは空気の熱伝導率、Prはプラントル数、ｇは重力加速度、βは体積膨張率、vは動粘性係数（図６参照）、ｘは密閉された空間の下からの位置までの距離を示す。この平均熱伝達率ｈｍは、対応するメッシュ領域における非定常二次元熱伝達式の結果に加算される。
【００７６】
▲３▼冷媒の融解時間及び融解速度式
図１５のボタン“融解表示”の操作で説明したように、ここでの計算は、温度シミュレーションの前段として、冷媒が融解する時間を計算により予測する。予測される融解時間は、温度シミュレーションに反映される。具体的には、予測される融解時間後に冷媒温度が変化し始めるように、非定常二次元熱伝達式での冷媒温度に反映される。
【００７７】
保冷容器への侵入熱量Qは
0003835980.tif 000015
ここで、A，ΔT，Hはそれぞれ伝熱面積、内外温度差、経過時間を示し、ｋは次式で表される熱貫流率を示す。
【００７８】
【数１４】
0003835980.tif 000016
ここで、λbox，ｄは保冷容器の熱伝導率及び厚さであり、α1，α2はそれぞれ壁面の外側及び内側の境膜係数である（図７〜９参照）。
【００７９】
氷、ドライアイス等の畜冷材が保冷容器に存在する間はこの侵入熱量Qが全てこれら畜冷材に吸収されるという条件の下で、氷、ドライアイス等の畜冷材が融解し終わるまでの融解時間を以下のように算出する。氷、ドライアイス等の畜冷材の全体がその周りから単位時間当たりに吸収する熱量が単位時間当たりの侵入熱量より多ければ周りの空気は冷やされる。すなわち、蓄冷材の熱伝達係数をhiで表すと、
0003835980.tif 000017
つまり、
【００８０】
【数１５】
0003835980.tif 000018
のとき容器内部の温度Tinは上式の両辺が等しくなるまで、すなわち平衡になるまで下がろうとする。ここで、蓄冷材の周りの熱伝達係数hiは次式で表される。
【００８１】
【数１６】
0003835980.tif 000019
結局、熱平衡になったときの容器内部の温度Tinは
【００８２】
【数１７】
0003835980.tif 000020
を満足する。いま、平衡状態を保ちつつ氷、ドライアイス等の畜冷材が融解する場合を考える。氷、ドライアイス等の畜冷材が完全に融解するまでの間に吸収する総熱量は、それらが氷点下から融解温度まで上昇するのに要する熱量と、潜熱としてもつ熱量の合計熱量に等しくなるまでの時間である。
【００８３】
２次元解析として、初期温度T0、周囲長さλ1、断面積Si、体積Viの氷を考え、その潜熱をL、密度をρとし、伝熱面積Aiを通過する熱流量をQiとすると、
0003835980.tif 000021
(δは奥行き単位長さ)となり、氷等物体の融解時間tは
【００８４】
【数１８】
0003835980.tif 000022
となる。ここで、融解温度以下のときでもqi=hi(T-273)と近似した。
同じ雰囲気に熱流束qが同じで大きさの異なる
【００８５】
【数１９】
0003835980.tif 000023
【００８６】
【数２０】
0003835980.tif 000024
で定義される氷iの融解速度vxiは
【００８７】
【数２１】
0003835980.tif 000025
により氷jの融解速度vxjに等しく、融解時間は
【００８８】
【数２２】
0003835980.tif 000026
によりtj>tiである。
【００８９】
さらに、初期温度がそれぞれTi0，Tj0の氷i,jの元体積をVi，Vjとし、氷iが全て融ける時間をti，その間に氷jが融ける体積をVjtiとすると
0003835980.tif 000027
となる。氷i,jについてそれぞれの熱収支を考えると
0003835980.tif 000028
【００９０】
【数２３】
0003835980.tif 000029
(5),(6)より
【００９１】
【数２４】
0003835980.tif 000030
を得る。一方、氷i,jへの熱流束qi，qjが異なる場合を考える。
【００９２】
0003835980.tif 000031
(8)より(7)と同様に
【００９３】
【数２５】
0003835980.tif 000032
を得る。(9)(10)より
【数２６】
0003835980.tif 000033
以上の考察では、氷i,jの初期温度Ti0，Tj0が融解温度以下の時でも熱流束は融解温度のときと同じと仮定している。外部からの侵入熱量が全て氷等の融解に使われるとする仮定の下では氷の周辺温度は変化しない。従って、氷が融解温度にあるときは、熱流束qi=hi(T-Ti)は一定である。
【００９４】
因みに、周辺温度10℃、氷の初期温度−13℃、氷の高さ10cm、熱伝導率λ=0.02305W/(mK)、動粘性係数ν=1.22×10-5m2/sのとき、熱伝達係数hは、
【００９５】
【数２７】
0003835980.tif 000034
である。氷の初期温度が0℃の場合、熱伝達係数hはh=2.64×10^0.25=4.6946となる。すなわち、熱伝達係数hは温度差に関して
【００９６】
【数２８】
0003835980.tif 000035
に比例して変化する。
【００９７】
氷の初期温度−13℃の場合は0℃の場合より熱伝達係数は1.23倍大きくなる。氷の初期温度−5℃の場合は1.5^0.25=1.1倍になるだけである。したがって、熱流束は氷の初期温度−13℃の場合は2.829倍、−5℃の場合は1.65倍大きい。
また、融解速度vxを
【００９８】
【数２９】
0003835980.tif 000036
で定義すると
【００９９】
【数３０】
0003835980.tif 000037
ゆえ、
【０１００】
【数３１】
0003835980.tif 000038
を得る。したがって、融解速度vxは上式に融解時間を代入すれば求めることができる。以下に、実際に例をとって融解時間を計算する。
【０１０１】
外気温度25℃のとき、容器内に高さ2cm幅47cm、初期温度-10℃の氷1.3 kgを置いたときの融解時間は次のように計算される。ただし、氷の潜熱80kcal/kg、密度917kg/m3、比熱2.0kJ/kg・Kとする。
【０１０２】
二分法により、熱平衡になったときの内部の温度は279.2K(6.2℃)となる。
このとき、熱伝達率は
【０１０３】
【数３２】
0003835980.tif 000039
２次元計算のため重量から奥行き長さを換算する。すなわち、
【０１０４】
【数３３】
0003835980.tif 000040
したがって、融解時間は
【０１０５】
【数３４】
0003835980.tif 000041
となる。熱伝達率は自然対流熱伝達とし、左右から熱を吸収するとして２倍した。
【０１０６】
図２０は、本実施の形態による断熱容器の設計システムのシミュレーションの精度を例示する図である。ここでは、外気温度３０℃の環境で１０時間輸送する場合を例に、段ボール容器内の温度の解析値と実測値、発泡プラスチック（ＥＰＳ）容器内の温度の解析値と実測値とを示している。
【０１０７】
段ボール容器に関して、輸送３時間後から解析値と実測値の温度差が少なくなり輸送終了時には解析値は実測値に実質的に一致している。また、ＥＰＳ容器に関して、輸送６時間迄は、解析値は実測値と約±２℃以内の精度を保っているが、輸送６時間経過後は実測値が解析値を次第に上回り、輸送終了時には約５℃程度の温度差が生じることがわかる。この温度差を調節するために、自然対流補正定数を必要に応じて変更する（図１６参照）。
【０１０８】
以上、図１９に示す容器、食材、冷媒の寸法及び配置に対応する温度分布データ（メッシュ図、初期温度分布図、温度伝導率分布図及び温度分布図）が温度保持性能データとして性能ＤＢ１５に格納されると共に、クライアント端末１に表示される（図２、ステップＳ１１１）。
【０１０９】
つぎに、鮮度シミュレーションについて説明する。鮮度シミュレーション式は、温度シミュレーション式により予測した、食材の温度変化に基いて食品の鮮度劣化を予測する計算式であり、（鮮度予測式）＝（温度予測式）+（化学反応式）で表される。
【０１１０】
鮮度を予測するにあたり、食品の保存期間と温度いわゆるＴＴＴ（Time Temperature Tolerance）に基づいて、温度シミュレーションから予測した食品の温度における鮮度予測式F(t)は、化学反応の速度を示すアレニウス則より
(1)鮮度劣化速度係数K=A・EXP(−E/R/T)
(2)時間ｔにおける定量的鮮度をF(t)で表すとF(t)=F(0)EXP(-∫Kｄｔ)で表現できる。
【０１１１】
ここで、微生物汚染は無視できることを仮定しており、Aは鮮度劣化に特有な係数、Eは活性化エネルギー（図１２参照）、Rは気体定数、Tは絶対温度を示す。F(0)=100を生産時の値として、F(t)>80を優鮮度、60<F(t)<80を良鮮度、F(t)=40を限界鮮度とする。ここで、算出される鮮度は、温度シミュレーションでの経過時間と対応して、グラフ（図示せず）にてクライアント端末１に表示される。
【０１１２】
以上、この発明の実施の形態を図面を参照して詳述してきたが、具体的な構成はこれらの実施の形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計の変更等があってもよい。
たとえば、上記実施の形態では、図２のステップＳ１０６〜Ｓ１０８で容器１における容器２及び３、食材１及び２、冷媒１の配置条件（図８〜図１１及び図１３参照）をユーザが入力し、解析を実行しているが、この配置条件のみを変更して解析を繰り返し行うこともできる。この配置条件の変更の度に、温度保持性能データは別な解析モデル（図１５参照）のデータとして格納される。容器１に対する容器２及び３、食材１及び２、並びに冷媒１の位置を変更してシミュレーションを行うことにより、保冷温度を維持し、容器１内の温度の均一化（最適配置）を図ることができる。
【０１１３】
図２１は、容器の最適配置を説明する図である。図２１（ａ）は、容器内温度の均一化を図るための最適配置を示す配置図表示画面であり、図２１（ｂ）は、該最適配置に基づく温度分布図である。
図２１（ａ）に示す最適配置は、食品と、該食品の周囲を覆うように配置される保冷材と、該食品及び保冷材を収容する容器とがＸ軸上で等間隔で配置されている。すなわち、食品の左端から保冷材（容器）の左端までの距離と、食品の右端から保冷材（容器）の右端までの距離とが等しく配置されている。図２１（ｂ）に示す温度分布は、濃い色が温度の低さを、薄い色が温度の高さを示している。
【０１１４】
図２２は、図２１に示した最適な容器配置と通常の容器配置との温度比較を示す図である。図２２は、外気温度の変化に応じた、容器内の下部に配置される食材の温度変化を示している。外気温度が２５℃から３０℃に変化する３時間経過後から次第に最適な容器配置による容器１内の温度と、通常の容器配置による容器１内の温度との間に差が生じており、輸送終了時点で約３℃程度の温度差が生じることがわかる。
【０１１５】
図２３は、図２１に示した最適な容器配置に対応する通常の容器配置の温度分布を示す図である。経過時間毎の温度分布図は図２２での外気温度に対応しており、図２３（ａ）は外気２５℃での環境で３時間経過した場合の温度分布、図２３（ｂ）は外気３０℃での環境で３時間経過（全体で６時間経過）した場合の温度分布、図２３（ｃ）は外気３５℃での環境で３時間経過（全体で９時間経過）した場合の温度分布を示している。
【０１１６】
図２３（ａ）で上部（冷媒が多く存在する位置）及び中央部（食材が存在する位置）が、図２３（ｂ）から図２３（ｃ）への時間変化に応じて濃い色（温度が低い）から薄い色（温度が高い）に変化していることがわかる。
図２４は、図２１に示した最適な容器配置の温度分布を示す図である。図２４（ａ）〜図２４（ｃ）は、図２３（ａ）〜図２３（ｃ）に対応し、外気２５℃での環境で３時間経過した場合の温度分布、外気３０℃での環境で３時間経過（全体で６時間経過）した場合の温度分布、外気３５℃での環境で３時間経過（全体で９時間経過）した場合の温度分布を各々示している。上部及び中央部の温度上昇のレベルが、図２３のそれと比べて緩和されており、容器１内の温度の均一化が図られている。
【０１１７】
また、上記実施の形態では、クライアント端末１のユーザにより図６〜図１１の条件入力画面に必要なデータを入力しているが、本発明はこれに限定されない。いずれにしても、所定の輸送時間経過後における容器内の温度分布を知ることができればよく、断熱容器の性能（図７〜図９参照）、該断熱容器の周囲の外気温度（図６参照）、該断熱容器に収納する内容物の熱容量及び初期温度（図１０及び図１１参照）、並びに該内容物の温度保持のための熱媒体の性能（図１３参照）が条件として入力されていればよい。
【０１１８】
また、上記実施の形態では、サーバ５は、クライアント端末１からの条件入力に応答して温度及び鮮度シミュレーションを行い、温度保持性能データを該端末１に表示させている。他の形態として、サーバ５は、クライアント端末１からの条件入力に応答して、過去に実行された解析モデル毎の性能データを利用してユーザに提供してもよい。
【０１１９】
この場合、サーバ５の制御部７は、クライアント端末１からの条件入力に応答して性能ＤＢ１５を参照し、該ＤＢ１５に格納されている温度保持性能データのうち、クライアント端末１からの条件に近い条件を有する温度保持性能データを読み出して、クライアント端末１に送信する。いずれにしても本発明は、クライアント端末１からの条件を受信し、該受信手段により受信されたデータに基づいて、温度保持性能データをクライアント端末１に送信するものであればよい。
【０１２０】
また、上記他の実施の形態では、過去に実行された解析モデル毎の性能データを利用してユーザに提供するが、その後条件に合致した温度保持性能データの要求をクライアント端末１から受付け、申し込みがあった場合に限り上述した温度及び鮮度シミュレーションを実行してもよい。これにより、ユーザが該シミュレーションによる温度保持性能データの評価後に、実際の試作品の作成依頼を受付け、作成依頼があった場合に限り試作品を作成して提供又は販売することも可能となる。
【０１２１】
また、上記実施の形態では、サーバ５は、インターネット３を介してクライアント端末１からの入力条件を受信し、該入力条件に応じた温度保持性能データを該端末１に送信している。更に他の形態として、サーバ５の入力操作部１７からの直接の条件入力に応じて、制御部７は温度／鮮度シミュレーションを行い、温度保持性能データを表示部１９に表示する構成としてもよい。又は、制御部７は、性能ＤＢ１５に格納されている温度保持性能データのうち、入力操作部１７からの条件に近い条件を有する温度保持性能データを読み出し、表示部１９に表示する構成としてもよい。
【０１２２】
同様にして、この場合のサーバ５の構成及び機能をクライアント端末１に設けることで、クライアント端末１は、ユーザからの条件入力に応じて、上述した温度／鮮度シミュレーションを行い、温度保持性能データを表示することができる。これにより、断熱容器の開発期間を短縮することができる。
また、本発明の断熱容器の設計システムは、本断熱容器の設計システムを機能させるためのプログラムでも実現される。このプログラムは、コンピュータで読み取り可能な記録媒体に格納されている。
【０１２３】
本システムの動作を実行するプログラムを記録した記録媒体は、図１に示されるサーバ５のＲＯＭ１３又はクライアント端末１のＲＯＭ（図示せず）そのものであっても良いし、また、外部記憶装置としてＣＤ−ＲＯＭドライブ２１等のプログラム読み取り装置が設けられ、そこに記録媒体を挿入することで読み取り可能なＣＤ−ＲＯＭ２３等であっても良い。
また、上記記録媒体は、磁気テープ、カセットテープ、フロッピーディスク、ハードディスク、ＭＯ／ＭＤ／ＤＶＤ等、又は半導体メモリであっても良い。
【０１２４】
【発明の効果】
本発明により、容器内の食品の温度及び鮮度シミュレーションに基づいた蓄冷材の配置（上部、下部又は側面）及び個数、重量、容器の材質及び肉厚（上部、下部又は側面）並びに形状（凸凹）解析により、対策コストを最小とする最適容器設計が可能となる。
また、各種パラメータの違いによる性能比較が可能となる。さらには、容器の製造に係る開発期間を大幅に短縮することが可能となる。
【図面の簡単な説明】
【図１】本発明の実施の形態による断熱容器の設計システムの構成を説明するブロック図である。
【図２】本発明の実施の形態による断熱容器の設計システムの動作を説明するフローチャートである。
【図３】初期画面を説明する図である。
【図４】容器数入力画面を説明する図である。
【図５】容器構成指定画面を説明する図である。
【図６】外気条件入力画面を説明する図である。
【図７】容器（容器１）条件入力画面を説明する図である。
【図８】容器（容器２）条件入力画面を説明する図である。
【図９】容器（容器３）条件入力画面を説明する図である。
【図１０】食材（食材１）条件入力画面を説明する図である。
【図１１】食材（食材２）条件入力画面を説明する図である。
【図１２】食材鮮度パラメータ計算条件画面を説明する図である。
【図１３】冷媒（冷媒１）条件入力画面を説明する図である。
【図１４】配置図表示画面を説明する図である。
【図１５】解析条件入力画面を説明する図である。
【図１６】自然対流補正画面を説明する図である。
【図１７】メッシュ分割指定画面を説明する図である。
【図１８】融解情報画面を説明する図である。
【図１９】解析実行時の温度分布を説明する図である。
【図２０】本発明の実施の形態による断熱容器の設計システムのシミュレーションの精度を例示する図である。
【図２１】図２１（ａ）は、容器内温度の均一化を図るための最適配置を示す配置図表示画面であり、図２１（ｂ）は、該最適配置に基づく温度分布図である。
【図２２】図２１に示した最適な容器配置と通常の容器配置との温度比較を示す図である。
【図２３】図２１に示した最適な容器配置に対応する通常の容器配置との温度分布を示す図である。
【図２４】図２１に示した最適な容器配置の温度分布を示す図である。
【符号の説明】
１  クライアント端末
３  インターネット
５  サーバ
７  制御部
９  送信／受信部
１１  ＲＡＭ
１３  ＲＯＭ
１５  性能データベース
１６  パラメータデータベース
１７  入力操作部
１９  表示部
２１  ＣＤ−ＲＯＭドライブ
２３  ＣＤ−ＲＯＭ
【図１】
0003835980.tif 000042
【図２】
0003835980.tif 000043
【図３】
0003835980.tif 000044
【図４】
0003835980.tif 000045
【図５】
0003835980.tif 000046
【図６】
0003835980.tif 000047
【図７】
0003835980.tif 000048
【図８】
0003835980.tif 000049
【図９】
0003835980.tif 000050
【図１０】
0003835980.tif 000051
【図１１】
0003835980.tif 000052
【図１２】
0003835980.tif 000053
【図１３】
0003835980.tif 000054
【図１４】
0003835980.tif 000055
【図１５】
0003835980.tif 000056
【図１６】
0003835980.tif 000057
【図１７】
0003835980.tif 000058
【図１８】
0003835980.tif 000059
【図１９】
0003835980.tif 000060
【図２０】
0003835980.tif 000061
【図２１】
0003835980.tif 000062
【図２２】
0003835980.tif 000063
【図２３】
0003835980.tif 000064
【図２４】
0003835980.tif 000065