(57)【要約】【課題】実際に線状体を配策した際の配策形状を精度よく再現することが可能な線状体の配策形状設計方法および装置を提供する。【解決手段】線状体を固定する位置であるクリップ位置を設定し、パラメトリック曲線を用いて各クリップ位置間の線状体の配策形状を作成する線状体の配策形状設計方法であって、パラメトリック曲線のコントロールポイントを変更しつつ、線状体の全体の歪エネルギーを演算することを繰り返し、歪エネルギーが最小となるコントロールポイントを決定し、各クリップ位置間の線状体の配策形状を作成する。【選択図】図１
2014092974.tif 000002
【特許請求の範囲】
【請求項１】
線状体を固定する位置であるクリップ位置を設定し、パラメトリック曲線を用いて前記各クリップ位置間の前記線状体の配策形状を作成する線状体の配策形状設計方法であって、
前記パラメトリック曲線のコントロールポイントを変更しつつ、前記線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となる前記コントロールポイントを決定し、前記各クリップ位置間の線状体の配策形状を作成する
ことを特徴とする線状体の配策形状設計方法。
【請求項２】
前記線状体の主要な形状決定要因が弾性的性質からなる
請求項１記載の線状体の配策形状設計方法。
【請求項３】
前記パラメトリック曲線が、隣接する前記クリップ位置の一方を時間ｔ＝０の始点、他方を時間ｔ＝１の終点とした時間ｔを媒介変数とした関数であり、［数１］に示す式（１）
【数１】
2014092974.tif 000011
で表される
請求項１または２記載の線状体の配策形状設計方法。
【請求項４】
線状体を固定する位置であるクリップ位置を設定し、パラメトリック曲線を用いて前記各クリップ位置間の前記線状体の配策形状を作成する線状体の配策形状設計装置であって、
前記パラメトリック曲線のコントロールポイントを変更しつつ、前記線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となる前記コントロールポイントを決定する最適コントロールポイント演算部と、
前記最適コントロールポイント演算部が決定した前記歪エネルギーが最小となるコントロールポイントに基づき、各クリップ位置間の線状体の配策形状を作成し出力する出力部と、
を備えたことを特徴とする線状体の配策形状設計装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、大型トラックのエアブレーキ配管用のナイロンチューブなどの線状体の配策形状を設計する線状体の配策形状設計方法および装置に関するものである。
【背景技術】
【０００２】
従来より、大型トラックのエアブレーキ配管用のナイロンチューブ等の線状体の配策形状を設計する際には、ＣＡＤ（Computer Aided Design、コンピュータ支援設計）システムを用い、ＣＡＤシステムに搭載されたスプライン曲線により配策形状設計を行うことが一般的である。なお、スプライン曲線とは、設定した複数の制御点を通る滑らかな曲線をいう。
【０００３】
この出願の発明に関連する先行技術文献情報としては、特許文献１，２がある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００１−２５０４３８号公報
【特許文献２】特開２００７−１０２６７６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述のスプライン曲線を用いた設計では、線状体を固定する位置であるクリップ位置間の線状体の配策形状を設計する際に、線状体の空間の経由点を設計者が感覚で設定するしかなく、設計者により異なる配策形状となってしまい、実際に線状体を配策した際の配策形状とも大きく異なってしまっていた。
【０００６】
本発明は上記事情に鑑み為されたものであり、実際に線状体を配策した際の配策形状を精度よく再現することが可能な線状体の配策形状設計方法および装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は上記目的を達成するために創案されたものであり、線状体を固定する位置であるクリップ位置を設定し、パラメトリック曲線を用いて前記各クリップ位置間の前記線状体の配策形状を作成する線状体の配策形状設計方法であって、前記パラメトリック曲線のコントロールポイントを変更しつつ、前記線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となる前記コントロールポイントを決定し、前記各クリップ位置間の線状体の配策形状を作成する線状体の配策形状設計方法である。
【０００８】
前記線状体の主要な形状決定要因が弾性的性質からなってもよい。
【０００９】
前記パラメトリック曲線が、隣接する前記クリップ位置の一方を時間ｔ＝０の始点、他方を時間ｔ＝１の終点とした時間ｔを媒介変数とした関数であり、［数１］に示す式（１）で表されてもよい。
【００１０】
【数１】
2014092974.tif 000003
【００１１】
また、本発明は、線状体を固定する位置であるクリップ位置を設定し、パラメトリック曲線を用いて前記各クリップ位置間の前記線状体の配策形状を作成する線状体の配策形状設計装置であって、前記パラメトリック曲線のコントロールポイントを変更しつつ、前記線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となる前記コントロールポイントを決定する最適コントロールポイント演算部と、前記最適コントロールポイント演算部が決定した前記歪エネルギーが最小となるコントロールポイントに基づき、各クリップ位置間の線状体の配策形状を作成し出力する出力部と、を備えた線状体の配策形状設計装置である。
【発明の効果】
【００１２】
本発明によれば、実際に線状体を配策した際の配策形状を精度よく再現することが可能な線状体の配策形状設計方法および装置を提供できる。
【図面の簡単な説明】
【００１３】
【図１】本発明の一実施の形態に係る線状体の配策形状設計方法のフロー図である。
【図２】Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線を説明する図である。
【図３】図１の線状体の配策形状設計方法で用いるパラメトリック曲線を説明する図である。
【図４】図１の線状体の配策形状設計方法における歪エネルギーの演算を説明する図である。
【図５】本発明の一実施の形態に係る線状体の配策形状設計装置の概略構成図である。
【発明を実施するための形態】
【００１４】
以下、本発明の実施の形態を添付図面にしたがって説明する。
【００１５】
図１は、本実施の形態に係る線状体の配策形状設計方法のフロー図である。
【００１６】
図１に示すように、本実施の形態に係る線状体の配策形状設計方法では、まず、ステップＳ１にて、線状体を固定する位置であるクリップ位置の設定を行う。具体的には、クリップ位置のＸＹＺ座標と線状体の向きを示す４座標を入力する。以下のステップでは、ここで設定した各クリップ位置間の線状体の配策形状をパラメトリック曲線を用いて作成することになる。
【００１７】
ステップＳ２では、線状体の全長、外径、内径（線状体が中空構造である場合）、ヤング率を入力する。詳細は後述するが、本実施の形態では、線状体が完全弾性体であると仮定して演算を行うため、線状体はその形状の決定の主要因が弾性的性質によるものを用いることが望ましい。本実施の形態では、線状体として、大型トラックのエアブレーキ配管用のナイロンチューブを用いた。なお、線状体が弾性体でない場合にも本発明は適用可能であるが、線状体が塑性変形の顕著な材料や粘性による形状変化の顕著な材料からなる場合には、実際の配策形状とは無視できない差異が発生するおそれがある。
【００１８】
その後、ステップＳ３にて、パラメトリック曲線を用いて各クリップ位置間の線状体の配策形状を作成する。ステップＳ３で用いるパラメトリック曲線の詳細については、後述する。
【００１９】
ステップＳ４では、ステップＳ３で作成した線状体の配策形状における線状体全体（クリップ位置間の線状体全体）の歪エネルギーを演算し、歪エネルギーを、各コントロールポイントのＸＹＺ座標についてそれぞれ微分する。なお、コントロールポイントとは、つまり、パラメトリック曲線の形状を決定づけるパラメータである。ステップＳ４での歪エネルギーの演算の詳細については、後述する。
【００２０】
ステップＳ５では、最急降下法により、パラメトリック曲線の各コントロールポイントを全体の歪エネルギーが減少する方向へ移動する。また、最急下降法は、関数の最小値を探索する勾配法のアルゴリズムであり公知であるため、ここでは詳細な説明を省略する。
【００２１】
その後、ステップＳ６にて、全方向（Ｘ軸方向、Ｙ軸方向、Ｚ軸方向）の歪エネルギーの微分値が０になったかを判定する。ステップＳ６でＮＯと判定された場合ステップＳ４に戻る。ステップＳ６でＹＥＳと判定された場合、その時点でのコントロールポイントが歪エネルギーが最小となるコントロールポイントとなり、このコントロールポイントで作成される線状体の配線形状が、完全弾性体の取り得る自然な配線形状であることになる。
【００２２】
よって、ステップＳ６でＹＥＳと判定された場合、ステップＳ７にて、当該コントロールポイントでの線状体の配策形状を作成し、モニタ等の表示器に表示し、処理を終了する。
【００２３】
次に、図１のステップＳ３で用いるパラメトリック曲線について説明する。
【００２４】
本実施の形態では、パラメトリック曲線には、線状体の自然な配策形状をよく表現できること、コントロールポイント数ができるだけ少なく、計算負荷が低いこと、が要求される。
【００２５】
公知のパラメトリック曲線として、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線、ベジェ曲線、Ｂスプライン曲線、ＮＵＲＢＳ（Non-Uniform Rational B-Spline）曲線などが知られている。このうち、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線は、ベジェ曲線やその発展であるＢスプライン曲線、ＮＵＲＢＳ曲線とは異なり、物体の運動の軌跡から着想を得たものである。自然界の運動の軌跡は最小作用の原理に従っていることから、本発明者らは線状体の自然な配策形状をよく表現できると考え、このＦｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線をパラメトリック曲線として用いることを検討した。
【００２６】
Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線は、時間ｔを媒介変数とした関数であり、時間ｔ＝０の始点での位置ｓ（０）と速度ｖ（０）、および時間ｔ＝１の終点での位置ｓ（１）と速度ｖ（１）を境界条件として微分形との連立方程式を解いて得られる式であり、［数２］に示す式（２）で表される。これら位置ｓ（０），ｓ（１）および速度ｖ（０），ｖ（１）の４点が、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線のコントロールポイントとなる。
【００２７】
【数２】
2014092974.tif 000004
【００２８】
しかしながら、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線では、コントロールポイントが４点と少なく、図２に示すように、始点と終点との間で速度ｖ（ｔ）がｖ（０）からｖ（１）に徐々に近づくように補完されてしまうため、２点以上の曲げが発生するような境界条件では実際の配策形状を再現することができない。具体的には、例えば、電球型（略Ω字状）となるような配策形状は、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線で表現することはできない。
【００２９】
本発明者らは、ベジェ曲線におけるＢスプライン曲線と同様に、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線を複数本接続してコントロールポイントを増やす試みを行ったが、コントロールポイントの操作対象が局所的になってしまい、数値計算の振動や発散の要因となってしまった。
【００３０】
そこで、本発明者らは、鋭意検討の結果、図３に示すように、始点と終点の加速度ａ（０），ａ（１）を定義して新たなコントロールポイントとし、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線における定義式を時間ｔ＝０〜１を媒介変数とする５次関数に拡張することで、２点以上の曲げが発生するような境界条件においても実際の配策形状を精度よく再現することが可能となり、かつ、数値計算の安定性を確保できることを見出した。
【００３１】
すなわち、本実施の形態では、隣接するクリップ位置の一方を時間ｔ＝０の始点、他方を時間ｔ＝１の終点とした時間ｔを媒介変数とした関数であり、［数３］に示す式（１）で表されるものをパラメトリック曲線として用いた。このパラメトリック曲線では、始点の位置ｓ（０）、終点の位置ｓ（１）、始点の速度ｖ（０）、終点の速度ｖ（１）、始点の加速度ａ（０）、および終点の加速度ａ（１）の６点がコントロールポイントとなる。
【００３２】
【数３】
2014092974.tif 000005
【００３３】
ここで、式（１）の導出方法を説明しておく。
【００３４】
時間ｔでの位置ｓ（ｔ）は、係数ａ〜ｆを用いて下式（３）
ｓ（ｔ）＝ａｔ5＋ｂｔ4＋ｃｔ3＋ｄｔ2＋ｅｔ＋ｆ    ・・・（３）
で表される。式（３）より、時間ｔでの速度ｖ（ｔ）は、下式（４）
ｖ（ｔ）＝ｄｓ（ｔ）／ｄｔ
＝５ａｔ4＋４ｂｔ3＋３ｃｔ2＋２ｄｔ＋ｅ    ・・・（４）
となり、時間ｔでの加速度ａ（ｔ）は、下式（５）
ａ（ｔ）＝ｄｖ（ｔ）／ｄｔ
＝２０ａｔ3＋１２ｂｔ2＋６ｃｔ＋２ｄ    ・・・（５）
となる。
【００３５】
式（３）〜（５）にｔ＝０をそれぞれ代入すると、下式（６）〜（８）
ｆ＝ｓ（０）    ・・・（６）
ｅ＝ｖ（０）    ・・・（７）
ｄ＝１／２・ａ（０）    ・・・（８）
となり、係数ｄ〜ｆが求まる。
【００３６】
式（３）〜（５）にｔ＝１、係数ｄ〜ｆを代入すると、下式（９）〜（１１）
ｓ（１）＝ａ＋ｂ＋ｃ＋１／２・ａ（０）＋ｖ（０）＋ｓ（０）
・・・（９）
ｖ（１）＝５ａ＋４ｂ＋３ｃ＋ａ（０）＋ｖ（０）    ・・・（１０）
ａ（１）＝２０ａ＋１２ｂ＋６ｃ＋ａ（０）    ・・・（１１）
となる。
【００３７】
式（９），（１０），（１１）の連立方程式を解くと、係数ａ〜ｃは、下式（１２）〜（１４）
ａ＝−６ｓ（０）＋６ｓ（１）−３ｖ（０）−３ｖ（１）
−１／２・ａ（０）＋１／２ａ（１）  ・・・（１２）
ｂ＝１５ｓ（０）−１５ｓ（１）＋８ｖ（０）＋７ｖ（１）
＋３／２・ａ（０）−ａ（１）  ・・・（１３）
ｃ＝−１０ｓ（０）＋１０ｓ（１）−６ｖ（０）−４ｖ（１）
−３／２・ａ（０）＋１／２ａ（１）  ・・・（１４）
となる。
【００３８】
求めた係数ａ〜ｆを式（３）に代入し、行列式で表現すると、上述の式（１）が得られる。
【００３９】
次に、図１のステップＳ４における歪エネルギーの演算について説明する。
【００４０】
ステップＳ３でパラメトリック曲線を用いて作成した線状体の配策形状において、図４に示すように、ｔ＝０からｔ＝１までの間を任意の区間（区間番号をｉとする）に区切って得られる座標群を（ｘi，ｙi、ｚi）とすると、引張圧縮歪εsは［数４］に示す式（１５）、曲歪εbは［数５］に示す式（１６）で表される。なお、ここでは線状体が完全弾性体であると仮定している。
【００４１】
【数４】
2014092974.tif 000006
【００４２】
【数５】
2014092974.tif 000007
【００４３】
式（１５）の引張圧縮歪εsと、線状体のヤング率Ｅ、断面積Ａ，無負荷区間長さＬnliより、下式（１７）
Ｐs＝（ＥＡεs）／Ｌnli      ・・・（１７）
により負荷力Ｐsが求まる。蓄えられる歪エネルギーは、負荷力Ｐsがなした仕事Ｗsに等しいので、引張圧縮歪（垂直歪）による歪エネルギーΔＵsは、［数６］に示す式（１８）で得られる。なお、式（１８）にて負荷力Ｐsに１／２を乗じているのは、負荷力Ｐsが引張圧縮歪εsに対して定数でなく一次比例で増加するためである。
【００４４】
【数６】
2014092974.tif 000008
【００４５】
他方、式（１６）の曲歪εbと、線状体のヤング率Ｅ、断面二次モーメントＩ、無負荷区間長さＬnliより、負荷トルクＭを求める。ミオソテスの式より、下式（１９），（２０）
εb＝（ＭＬnli2）／（２ＥＩ）    ・・・（１９）
θ＝（ＭＬnli）／（ＥＩ）    ・・・（２０）
が成り立つ。式（１９）を変形すると、下式（２１）
Ｍ＝（２ＥＩεb）／（Ｌnli2）    ・・・（２１）
が得られ、式（１９），（２０）より、下式（２２）
θ＝（２εb）／Ｌnli    ・・・（２２）
が得られる。
【００４６】
蓄えられる歪みエネルギーは負荷トルクＭがなした仕事Ｗbに等しいので、曲歪による歪エネルギーΔＵbは、［数７］に示す式（２３）となる。
【００４７】
【数７】
2014092974.tif 000009
【００４８】
全体の歪エネルギーＵは、引張圧縮歪による歪エネルギーΔＵsと曲歪による歪エネルギーΔＵbを全区間で足し合わせたものであるから、［数８］に示す式（２４）により、全体の歪エネルギーＵを求めることができる。
【００４９】
【数８】
2014092974.tif 000010
【００５０】
線状体を配策する際には、この全体の歪エネルギーＵが最小となるような形状に自然に移行する。よって、パラメトリック曲線のコントロールポイントを変更しつつ、線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となるコントロールポイントを決定することで、最も自然な線状体の配策形状を作成することが可能になる。
【００５１】
各クリップ位置間で線状体の配策形状を作成し、全体を繋ぎ合わせれば、線状体の全体の配策形状が完成する。
【００５２】
次に、本実施の形態に係る線状体の配策形状設計方法を実施する線状体の配策形状設計装置について説明する。
【００５３】
図５に示すように、線状体の配策形状設計装置５１は、クリップ位置の情報（座標、線状体の向き）や、線状体の全長、外径、内径（線状体が中空構造である場合）、ヤング率等の情報を入力する入力部５２と、入力部５２で入力された情報を基に、パラメトリック曲線を用いて各クリップ位置間の線状体の配策形状を作成する配策形状作成部５３と、配策形状作成部５３で作成した配策形状における歪エネルギーを演算する歪エネルギー演算部５４と、配策形状作成部５３および歪エネルギー演算部５４を制御し、パラメトリック曲線のコントロールポイントを変更しつつ、線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となるコントロールポイントを決定する最適コントロールポイント演算部５５と、最適コントロールポイント演算部が決定した歪エネルギーが最小となるコントロールポイントに基づき、各クリップ位置間の線状体の配策形状を作成し出力する出力部５６と、モニタ等の表示器５７と、を備えている。
【００５４】
入力部５２、配策形状作成部５３、歪エネルギー演算部５４、最適コントロールポイント演算部５５、出力部５６は、パーソナルコンピュータ等の演算装置に搭載され、ＣＰＵ、ソフトウェア、インターフェイス、メモリ等を適宜組み合わせて実現される。
【００５５】
線状体の配策形状設計装置５１は、ＣＡＤシステムに組み込まれてもよい。この場合、線状体の配策形状設計装置５１は、ＣＡＤシステムとハードウェアを共通としつつその機能がプログラム化されＣＡＤシステムに組み込まれることとなる。線状体の配策形状設計装置５１をＣＡＤシステムに組み込む場合、線状体の配策形状設計装置５１は、線状体の配策形状設計時にＣＡＤシステムを自動的に操作して演算結果として線状体の配策形状を表示するよう構成されることになる。
【００５６】
以上説明したように、本実施の形態に係る線状体の配策形状設計方法では、パラメトリック曲線のコントロールポイントを変更しつつ、線状体の全体の歪エネルギーを演算することを繰り返し、当該歪エネルギーが最小となるコントロールポイントを決定し、各クリップ位置間の線状体の配策形状を作成している。
【００５７】
すなわち、本実施の形態では、線状体が取りうる形状を歪エネルギー最小化問題として数値計算により解いているため、従来のように設計者により異なる配策形状となることがなく、実際に線状体を配策した際の配策形状を精度よく再現することが可能になり、実際に線状体を配策する際の不具合の発生を抑制できる。
【００５８】
また、本実施の形態では、従来のように線状体の空間の経由点を手作業で設定する必要がなく、自動的に線状体の配策形状を作成可能となるため、線状体の配策形状の設計時間を大幅に短縮することが可能になる。
【００５９】
さらに、本実施の形態では、上述の式（１）で表される新規なパラメトリック曲線を用いているため、始点と終点の加速度ａ（０），ａ（１）を新たなコントロールポイントとして使用することが可能となり、２点以上の曲げ点が発生するような境界条件であっても、実際の線状体の配策形状を精度よく再現することが可能になる。また、新たに定義した始点と終点の加速度ａ（０），ａ（１）は、線状体の配策形状全体に影響するため、Ｆｅｒｇｕｓｏｎ／Ｃｏｏｎｓ曲線を直列接続した計算モデルにおいて現れた数値計算の不安定性も解消することができた。
【００６０】
本発明は上記実施の形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の変更を加え得ることは勿論である。
【００６１】
例えば、上記実施の形態では言及しなかったが、クリップ位置間に線状体が通過する中間点を設定可能とし、当該中間点を通る線状体の配策形状を作成するように構成してもよい。具体的には当該中間点を接続点とする二本のパラメトリック曲線となる。この場合、例えば、タイバンド（結束バンド）を用いて線状体を固定する場合など線状体の向きが可変となる中間点を設定する場合、その中間点での速度（速度のベクトルは線状体の向きを表す）の向きも変数となり、数値計算による歪エネルギー最小化問題探索計算での微分対象となる。このように適宜コントロールポイントを定数あるいは変数として条件設定することで実際の線状体拘束条件を計算モデルに再現し、パラメトリック曲線のコントロールポイントとして用いることで、中間点での線状体通過向きを非拘束とした配策条件での最適な配策形状を作成するよう構成することも可能である。
【００６２】
また、上記実施の形態では、歪エネルギーを最小とするコントロールポイントを探索するアルゴリズムとして最急降下法を用いたが、これに限定されるものではない。
【符号の説明】
【００６３】
５１  線状体の配策形状設計装置
５２  入力部
５３  配策形状作成部
５４  歪エネルギー演算部
５５  最適コントロールポイント演算部
５６  出力部
５７  表示器
【図１】
2014092974.tif 000012
【図２】
2014092974.tif 000013
【図３】
2014092974.tif 000014
【図４】
2014092974.tif 000015
【図５】
2014092974.tif 000016