(57)【特許請求の範囲】
【請求項１】
第１の歯車の歯面上の各計測点で計測された３次元座標データと、上記第１の歯車に噛合する第２の歯車の歯面上の各計測点で計測された３次元座標データとを歯車対の組立諸元を用いて所定の噛合回転位置で互いに関連付け、上記各３次元座標データを上記第１の歯車を基準とする円筒座標系上の３次元座標データに変換する座標変換手段と、
上記第２の歯車の歯面上の各３次元座標データに基づいて当該第２の歯車の歯面上の点を表す関数を作成する関数作成手段と、
上記第１の歯車の歯面上の点に対応して上記円筒座標系上で同一周上に存在する上記第２の歯車の歯面上の点を上記関数を用いて検索し、これら２点間の隙間を示す相対角度情報を算出する歯面間隙間情報演算手段とを備えたことを特徴とする歯車対の評価装置。
【請求項２】
上記座標変換手段は、歯車対の組立諸元を用いて互いに関連付けた上記各３次元座標データを、歯車対のデフレクション値で補正することを特徴とする請求項１記載の歯車対の評価装置。
【請求項３】
上記第１の歯車と上記第２の歯車の各噛合回転位置で演算される上記各相対角度情報の中から最小の上記相対角度情報を抽出し、該抽出した最小の相対角度情報を基準とし、上記各相対角度情報に基づいて、上記第１，第２の歯車の歯面の噛合開始から終了までの相対的な隙間距離を示す３次元の包絡面を演算する包絡面演算手段を備えたことを特徴とする請求項１または請求項２に記載の歯車対の評価装置。
【請求項４】
上記包絡面演算手段は、上記抽出した最小の相対角度情報に対応する上記包絡面上の点を歯車対の最凸点として抽出することを特徴とする請求項３記載の歯車対の評価装置。
【請求項５】
上記包絡面演算手段は、演算した３次元の包絡面の情報を２次元の等高線の情報に変換することを特徴とする請求項３または請求項４記載の歯車対の評価装置。
【請求項６】
上記関数作成手段は、上記第２の歯車の歯面上に設定した２次元の媒介変数を用い、上記第２の歯車の歯面上の各３次元座標データに基づいて当該第２の歯車の歯面上の点を表す半径座標の関数、軸座標の関数、及び角度座標の関数を作成し、
上記歯面間隙間情報演算手段は、上記第１の歯車の歯面上の点と上記円筒座標系上で同一周上に存在する上記第２の歯車の歯面上の点を示す上記媒介変数を、上記半径座標の関数及び上記軸座標の関数から演算し、当該演算した媒介変数に基づき、上記第１の歯車の歯面上の点と当該点に対応する上記第２の歯車の歯面上の点との間の隙間を示す相対角度情報を上記角度座標の関数を用いて算出することを特徴とする請求項１乃至請求項５の何れか１項に記載の歯車対の評価装置。
【請求項７】
上記関数作成手段は、第２歯面上の各格子点の座標データに基づいて、互いに近接する３点の格子点の組み合わせに基づいて歯面を近似する平面式にて設定されることを特徴とする請求項１に記載の歯車対の評価装置。
【請求項８】
第１の歯車の歯面上の各計測点で計測された３次元座標データと、上記第１の歯車に噛合する第２の歯車の歯面上の各計測点で計測された３次元座標データとを歯車対の組立諸元を用いて所定の噛合回転位置で互いに関連付け、上記各３次元座標データを上記第１の歯車を基準とする円筒座標系上の３次元座標データに変換する座標変換手順と、
上記第２の歯車の歯面上の各３次元座標データに基づいて当該第２の歯車の歯面上の点を表す関数を作成する関数作成手順と、
上記第１の歯車の歯面上の点に対応して上記円筒座標系上で同一周上に存在する上記第２の歯車の歯面上の点を上記関数を用いて検索し、これら２点間の隙間を示す相対角度情報を算出する歯面間隙間情報演算手順とを備えたことを特徴とする歯車対の評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、特に、傘歯車やハイポイドギヤ等の歯当り解析に好適な歯車対の評価装置、歯車対の評価プログラム、及び、これらを用いて歯面評価した歯車対に関する。
【背景技術】
【０００２】  一般に、傘歯車やハイポイドギヤ等の歯車対において、実用に供される各歯車の歯面形状（実歯面形状）は、理論的に共役な数学的に表現される歯面ではなく、加工可能な近似歯面形状となる。従って、特に、この種の歯車対では、実歯面の歯当り解析に基づいて、その強度評価や、振動騒音評価等を行うことが重要となってくる。【０００３】  このような歯車対の実歯面の歯当りを解析する技術として、例えば、非特許文献１（久保愛三・他５名「ハイポイドギヤ・ベベルギヤの運転性能解析法（第１報，歯面形状定義基準面のとり方と解析法）」、日本機械学会論文集（Ｃ編）６２巻５９９号（１９９６−７）、論文Ｎｏ．９５−１５４７、Ｐ２８３３〜Ｐ２８４１）には、実歯面のデータと歯面加工時に歯切り盤にセッティングされる基準歯面（理論歯面形状）のデータとの偏差を計算し、この偏差を基準歯面に対する実歯面の形状偏差として捉えることにより、基準歯面の噛合い状態を基礎として実歯面の噛合い状態を解析する技術が開示されている。【０００４】  しかしながら、上述の非特許文献１に開示された技術は、基準歯面の情報を必要とするため、例えば、各種諸元の設計から実際の歯面加工までの一連の作業をユーザ自らが行った歯車対のように、基準歯面の情報が明確な歯車対に対しては有効であるものの、基準歯面の情報が不明な歯車対に対しては適用が困難である。【０００５】  すなわち、上述の非特許文献１に開示された技術は、互いに噛合する歯車の各基準歯面を把握する必要があり、これらの基準歯面は加工機によって異なるため、加工機のセッティングを把握しない限りは、実際の歯車対の歯当りを評価することが困難であった。【０００６】
本発明は上記事情に鑑みてなされたもので、基準歯面の情報を用いることなく、実歯面の計測情報に基づいて精度のよい歯面解析を実現することができる歯車対の評価装置、及び、歯車対の評価プログラムを提供することを目的とする。
【発明の開示】
【課題を解決するための手段】
【０００７】
本発明の一例による歯車対の評価装置は、第１の歯車の歯面上の各計測点で計測された３次元座標データと、上記第１の歯車に噛合する第２の歯車の歯面上の各計測点で計測された３次元座標データとを歯車対の組立諸元を用いて所定の噛合回転位置で互いに関連付け、上記各３次元座標データを上記第１の歯車を基準とする円筒座標系上の３次元座標データに変換する座標変換手段と、上記第２の歯車の歯面上の各３次元座標データに基づいて当該第２の歯車の歯面上の点を表す関数を作成する関数作成手段と、上記第１の歯車の歯面上の点に対応して上記円筒座標系上で同一周上に存在する上記第２の歯車の歯面上の点を上記関数を用いて検索し、これら２点間の隙間を示す相対角度情報を算出する歯面間隙間情報演算手段とを備えたものである。
【０００８】
本発明の一例による歯車対の評価プログラムは、第１の歯車の歯面上の各計測点で計測された３次元座標データと、上記第１の歯車に噛合する第２の歯車の歯面上の各計測点で計測された３次元座標データとを歯車対の組立諸元を用いて所定の噛合回転位置で互いに関連付け、上記各３次元座標データを上記第１の歯車を基準とする円筒座標系上の３次元座標データに変換する座標変換手順と、上記第２の歯車の歯面上の各３次元座標データに基づいて当該第２の歯車の歯面上の点を表す関数を作成する関数作成手順と、上記第１の歯車の歯面上の点に対応して上記円筒座標系上で同一周上に存在する上記第２の歯車の歯面上の点を上記関数を用いて検索し、これら２点間の隙間を示す相対角度情報を算出する歯面間隙間情報演算手順とを備えたものである。
【図面の簡単な説明】
【００１１】
【図１】本発明の第１の実施形態に係わり、歯車対の評価装置の概略構成図
【図２】同上、歯車対の評価装置を実現するためのコンピュータシステムの一例を示す概略構成図
【図３】同上、歯当り解析ルーチンを示すフローチャート
【図４】同上、歯面間隙間演算サブルーチンを示すフローチャート
【図５】同上、包絡面演算サブルーチンを示すフローチャート
【図６】同上、ハイポイドギヤの斜視図
【図７】同上、ギヤ及びピニオンの各歯面上の各格子点を規定する円柱座標系を示す説明図
【図８】同上、ギヤ歯面上の格子点とピニオン歯面上の収束点との関係を示す説明図
【図９】同上、曲面座標の計算方法を示す説明図
【図１０】同上、ギヤ歯面を基準とした相対歯面を示す説明図
【図１１】同上、各ピニオン回転ステップでの相対歯面を示す説明図
【図１２】同上、図１１の各相対歯面を合成した包絡面を示す説明図
【図１３】同上、歯車対の歯面距離分布を示す説明図
【図１４】本発明の第２の実施形態に係わり、歯面間隙間演算サブルーチンを示すフローチャート
【図１５】同上、ギヤ歯面上の格子点とピニオン歯面に設定した平面上の点との関係を示す説明図
【図１６】同上、ギヤ歯面上の格子点に対応するピニオン歯面上の点が平面上に存在するときの条件を示す説明図
【図１７】同上、包絡面演算サブルーチンを示すフローチャート
【図１８】同上、ギヤとピニオンの歯当り領域及び接触点の軌跡を示す説明図
【図１９】同上、ギヤとピニオンとの間の伝達誤差を示す説明図
【発明を実施するための最良の形態】
【００１２】  以下、図面を参照して本発明の形態を説明する。図１乃至図１３は本発明の第１の実施形態に係わり、図１は歯車対の評価装置の概略構成図、図２は歯車対の評価装置を実現するためのコンピュータシステムの一例を示す概略構成図、図３は歯当り解析ルーチンを示すフローチャート、図４は歯面間隙間演算サブルーチンを示すフローチャート、図５は包絡面演算サブルーチンを示すフローチャート、図６はハイポイドギヤの斜視図、図７はギヤ及びピニオンの各歯面上の各格子点を規定する円柱座標系を示す説明図、図８はギヤ歯面上の格子点とピニオン歯面上の収束点との関係を示す説明図、図９は曲面座標の計算方法を示す説明図、図１０はギヤ歯面を基準とした相対歯面を示す説明図、図１１は各ピニオン回転ステップでの相対歯面を示す説明図、図１２は図１１の各相対歯面を合成した包絡面を示す説明図、図１３は歯車対の歯面距離分布を示す説明図である。【００１３】  図６において、符号１００は歯車対を示し、本実施形態において、歯車対１００は、例えば、大径をなす一方の歯車である第１の歯車（以下、ギヤともいう）１０１Ｇと、小径をなす他方の歯車である第２の歯車（以下、ピニオンともいう）１０１Ｐとが互いに噛合するハイポイドギヤである。【００１４】  この歯車対１００は、例えば図１に示す評価装置１を用いて評価される。評価装置１は、実際の歯車対１００の情報としてギヤ１０１Ｇ及びピニオン１０１Ｐの各歯面の３次元座標データ（実測値）やディメンジョンデータ等を入力するための入力部５と、入力された歯車対情報に基づいて各種演算等を行う演算部６と、演算部６で実行される各種プログラムを格納するとともに、入力された歯車対情報や演算部６での演算結果等を適宜記憶する記憶部７と、演算部６での演算結果等を出力する出力部８とを有して構成されている。【００１５】  具体的に説明すると、評価装置１には、ギヤ１０１Ｇの実歯面情報として、例えば、ギヤ１０１Ｇの注目する歯面（ギヤ歯面）１０２Ｇ上に設定したｊ×ｉ個（例えば、歯丈方向に１５個×歯筋方向に１５個）の各格子点でそれぞれ計測された３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）が入力される。また、評価装置１には、ピニオン１０１Ｐの実歯面情報として、例えば、ピニオン１０１Ｐの注目する歯面（ピニオン歯面）１０２Ｐ上に設定したｊ×ｉ個（例えば、歯丈方向に１５個×歯筋方向に１５個）の各格子点でそれぞれ計測された３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）が入力される。ここで、ギヤ１０１Ｇの実歯面情報として評価装置１に入力される各３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）は、ギヤ１０１Ｇの回転軸上に原点ＯＧを持つＸ−Ｙ−Ｚの直交座標系で規定される座標データであり、本実施形態において、この直交座標系のＹ軸は、ギヤ歯面１０２Ｇの中心を貫通するよう設定され、Ｚ軸はギヤ１０１Ｇの回転軸と同軸上に設定されている。同様に、ピニオン１０１Ｐの実歯面情報として評価装置１に入力される各３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）は、ピニオン１０１Ｐの回転軸上に原点ＯＰを持つＸ−Ｙ−Ｚの直交座標系で規定される座標データであり、本実施形態において、この直交座標系の��軸はピニオン歯面１０２Ｐの中心を貫通するよう設定され、Ｚ軸はピニオン１０１Ｐの回転軸と同軸上に設定されている。尚、本実施例では歯面中央を貫通する座標系を設定しているが、これには限定されず適宜設定することができる。【００１６】  また、評価装置１には、歯車対１００のディメンジョンデータとして、ギヤ比ｒａｔｉｏ、組立諸元（オフセットＥ、交差角Σ）、及び、各デフレクション値δＥ，δΣ，δＧ，δＰ等が入力される。ここで、各デフレクション値は歯車対１００に所定トルクが加わった際の撓み等に起因するずれ量であり、δＥはオフセットＥのずれ量、δΣは交差角Σのずれ量、δＧはギヤ１０１Ｇの回転軸方向のずれ量、δＰはピニオン１０１Ｐの回転軸方向のずれ量である（図６参照）。【００１７】  そして、評価装置１は、ギヤ歯面１０２Ｇ上の各３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）と、ピニオン歯面１０２Ｐ上の各３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）とを歯車対１００の組立諸元を用いて所定の噛合回転位置で互いに関連付け、これらを、ギヤ１０１Ｇを基準とするＲ−Ｚ−Θの円筒座標系上の３次元座標データ（ｒＧｊｉ，ｚＧｊｉ，θＧｊｉ）及び（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に変換する（図７参照）。ここで、評価装置１は、ギヤ歯面１０２Ｇ上の各３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）と、ピニオン歯面１０２Ｐ上の各３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）とを互いに関連付ける際に、デフレクション値を用いた補正を行う。【００１８】  さらに、評価装置１は、ピニオン歯面１０２Ｐ上に、該ピニオン歯面１０２Ｐ上の各格子点の番号と関連付けて２次元の媒介変数（ｊ，ｉ）を設定し、これら媒介変数（ｊ，ｉ）を用い、各３次元座標データ（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に基づいて、ピニオン歯面１０２Ｐ上の点を表す半径座標の関数ｆＲ（ｊ，ｉ）、軸座標の関数ｆＺ（ｊ，ｉ）、及び角度座標の関数ｆΘ（ｊ，ｉ）を作成する。【００１９】  さらに、評価装置１は、Ｒ−Ｚ−Θの円筒座標系上において、ギヤ歯面１０２Ｇ上の各格子点と同一の周上に存在するピニオン歯面１０２Ｐ上の各点を示す各媒介変数（ｊ，ｉ）を、関数ｆＲ（ｊ，ｉ）及びｆＺ（ｊ，ｉ）からニュートン法を用いて演算し、演算した各媒介変数（ｊ，ｉ）を用いて関数ｆΘ（ｊ，ｉ）から得られる角度情報θＰｊｉを基に、所定の噛合回転位置においてギヤ歯面１０２Ｇ上の点（格子点）と当該格子点に対応するピニオン歯面１０２Ｐ上の点との隙間を示す相対角度情報（歯面間角度）を算出する。【００２０】  ここで、評価装置１は、ピニオン１０１Ｐの１ピッチを所定の分割数で分割した回転角（ピニオン回転角）θｓ毎に（すなわち、ピニオン回転角θｓで規定されるギヤ１０１Ｇとピニオン１０１Ｐの噛合回転位置毎に）、ギヤ歯面１０２Ｇ上の各格子点を基準とする相対歯面情報を演算する。そして、評価装置１は、各噛合回転位置で演算した相対歯面情報を合成することにより、ギヤ歯面１０２Ｇとピニオン歯面１０２Ｐとの噛合開始から終了までの相対的な隙間距離を示す包絡面を演算する。【００２１】  すなわち、評価装置１の記憶部７には、上述の各種演算等を行うためのプログラムが格納されており、演算部６は、これらのプログラムを実行することにより、座標変換手段、関数作成手段、歯面間隙間情報演算手段、及び、包絡面演算手段としての各機能を実現する。【００２２】  なお、本実施形態の評価装置１は、例えば図２に示すコンピュータシステム１０で実現される。コンピュータシステム１０は、例えば、コンピュータ本体１１に、キーボード１２と、表示手段の一例としてのディスプレイ装置１３と、プリンタ１４とがケーブル１５を介して接続されて要部が構成されている。そして、このコンピュータシステム１０において、例えば、コンピュータ本体１１に配設された各種ドライブ装置やキーボード１２等が入力部１５として機能するとともに、コンピュータ本体１１に内蔵されたＣＰＵ，ＲＯＭ，ＲＡＭ等が演算部６として機能する。また、コンピュータ本体１１に内蔵されたハードディスク等が記憶部７として機能するとともに、ディスプレイ装置１３やプリンタ１４等が出力部８として機能する。【００２３】  次に、演算部６で実行される歯当り解析について、図３に示す歯当り解析ルーチンのフローチャートに従って説明する。このルーチンがスタートすると、演算部６は、先ず、ステップＳ１０１において、ギヤ１０１Ｇの歯面１０２Ｇ上の各格子点で計測された３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）、ピニオン１０１Ｐの歯面１０２Ｐ上の各格子点で計測された３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）、ギヤ比ｒａｔｉｏ、オフセットＥ、交差角Σ、各デフレクション値δＥ，δΣ，δＧ，δＰ等の歯車対情報を読み込む。【００２４】  続くステップＳ１０２において、演算部６は、ピニオン１０１Ｐの１ピッチ当たりの分割数、及び、１ステップ当りのピニオン回転角θｓを、以下の（１）〜（４）式を用いて演算する。
0004823897.tif 000002
ここで、Ｓｗａｇｌは噛合い始点−終点間のピニオン回転角度、ｎＰはピニオン歯数、Ｓｗｎは同時噛合い歯数（噛合い率）、ｏｐｎは１噛合いピッチ当りの角度ステップ数（整数）、Ｃｈｎはグリッド間隔当りの角度ステップ数、Ｍαは補正値である。【００２５】  また、Ｍｍａｘは計算すべき角度の範囲を示すものであり、歯車回転最大ステップ数Ｍｍａｘとする。なお、（２）、（４）式において、「Ｒｏｕｎｄ」とは（  ）内の値を四捨五入するためのソフト関数である。【００２６】  そして、ステップＳ１０２からステップＳ１０３に進むと、演算部６は、歯車回転ステップ数Ｍ＝１をセットし、続くステップＳ１０４において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したか否かを調べる。【００２７】  そして、演算部６は、ステップＳ１０４において、歯車回転ステップ数Ｍが未だ歯車回転最大ステップ数Ｍｍａｘに達していないと判定した場合にはステップＳ１０５に進み、以下のＳ１１１までの処理をＭｍａｘ回繰り返すことにより、ピニオン１０１Ｐがθｓ（ｒａｄ）回転する噛合回転位置毎の相対歯面データを演算する。一方、ステップＳ１０４において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したと判定すると、演算部６は、ステップＳ１１２に進む。【００２８】  ステップＳ１０４からステップＳ１０５に進むと、演算部６は、現在の歯車回転ステップ数Ｍに対応するピニオン１０１Ｐ及びギヤ１０１Ｇの回転角度ｒｏｔＰ，ｒｏｔＧを、以下の（５）、（６）式を用いて演算する。
0004823897.tif 000003
【００２９】  続くステップＳ１０６において、演算部６は、ピニオン１０１Ｐの直交座標系上の各格子点の座標をＺ軸回りにｒｏｔＰ（ｒａｄ）移動させるとともに、ギヤ１０１Ｇの直交座標系上の各格子点の座標をＺ軸回りにｒｏｔＧ（ｒａｄ）移動させる。さらに、演算部６は、ピニオンの直交座標系を、ギヤ１０１Ｇの直交座標系と一致させた状態から（クロスポイントの基準位置から）、オフセットＥだけ移動させるとともに、交差角Σだけ回転させることにより、現ステップ回転位置における両歯面１０２Ｇ，１０２Ｐの関係をギヤ１０１Ｇの直交座標系上に設定する。すなわち、演算部６は、現ステップ回転位置におけるピニオン歯面１０２Ｐ上の各格子点の座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）を、ギヤ１０１Ｇの直交座標系上の座標データ（ｘＰ’ｊｉ，ｙＰ’ｊｉ，ｚＰ’ｊｉ）に変換する。その際、演算部６は、ピニオン１０１Ｐ及びギヤ１０１Ｇの各座標系の原点ＯＰ，ＯＧをデフレクション値δＰ，δＧで補正するとともに、オフセットＥをデフレクション値δＥで補正し、交差角Σをデフレクション値δΣで補正する。【００３０】  続くステップＳ１０７において、演算部６は、ステップＳ１０６でギヤ１０１Ｇ基準の直交座標系上に表したギヤ１０１Ｇ及びピニオン１０１Ｐの全計測点（全ギヤ格子点及びピニオン格子点）の座標データを、以下の（７）〜（１２）式を用いて、ギヤ１０１Ｇを基準とするＲ−Ｚ−Θの円筒座標系の座標データにそれぞれ変換する。【００３１】
0004823897.tif 000004
【００３２】  続くステップＳ１０８において、演算部６は、ピニオン歯面１０２Ｐ上に各格子点の番号と関連付けて２次元の媒介変数ｊ，ｉを設定し、これら媒介変数ｊ，ｉを用い、ピニオン歯面１０２Ｐの各格子点の座標データ（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に基づいて、ピニオン歯面１０２Ｐ上の格子点間を補間する関数ｆＲ（ｉ，ｊ）、ｆＺ（ｉ，ｊ）、ｆθ（ｉ，ｊ）を演算する。ここで、各関数ｆＲ（ｉ，ｊ）、ｆＺ（ｉ，ｊ）、ｆθ（ｉ，ｊ）は、例えば、それぞれスプライン関数で構成される。【００３３】  そして、ステップＳ１０８からステップＳ１０９に進むと、演算部６は、図４に示す歯面間隙間演算サブルーチンのフローチャートに従って、ギヤ歯面１０２Ｇ上の各格子点とこれら各格子点にピニオン歯面１０２Ｐ上で対応する点との間の歯面間角度（相対歯面情報）を演算する。ここで、以下の説明において、ギヤ歯面１０２Ｇ上の各格子点を規定するグリッド数（ｊ，ｉ）を、ピニオン歯面１０２Ｐ上のものと区別するため、これらを（ｊＧ，ｉＧ）と表記する。【００３４】  このサブルーチンがスタートすると、演算部６は、ステップＳ２０１において、媒介変数の初期値（ｊ＝ｊｉｎｉ，ｉ＝ｉｉｎｉ）を設定し、ステップＳ２０２において、初期値ｊｉｎｉ，ｉｉｎｉで（例えば、ｊｉｎｉ＝ｉｉｎｉ＝８）規定されるピニオン歯面１０２Ｐ上の基準点の座標ｒ１（＝Ｒ（ｊｉｎｉ，ｉｉｎｉ）），ｚ１（＝Ｚ（ｊｉｎｉ，ｉｉｎｉ））を求める。【００３５】  続くステップＳ２０３において、演算部６は、後述する歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算を、ギヤ歯面１０２Ｇ上の全格子点ｊＧ，ｉＧに対して行ったか否かを調べる。そして、ギヤ歯面１０２Ｇ上の全格子点に対する歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の演算が終了していないと判定すると、演算部６は、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算対象とする格子点を新たな格子点に更新した後、ステップＳ２０４に進み、以下のステップＳ２０９までの処理により、ニュートン法を用いて歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を演算する。一方、ステップＳ２０３において、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算をギヤ歯面１０２Ｇ上の全格子点（ｊＧ，ｉＧ）に対して行ったと判定すると、演算部６は、サブルーチンを抜け、メインルーチンに戻る。【００３６】  ステップＳ２０３からステップＳ２０４に進むと、演算部６は、今回歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算対象として選択したギヤ歯面１０２Ｇ上の（ｊＧ，ｉＧ）番目の格子点での３次元座標データを（ｒ０，ｚ０，θ０）とし、このｊＧ，ｉＧ番目の格子点とＲ−Ｚ−Θの円筒座標系上で同一周上に存在するピニオン歯面１０２Ｐ上の点（収束点）での３次元座標データを（ｒｃ，ｚｃ，θｃ）とし、さらに、（ｒｃ，ｚｃ，θｃ）＝（ｒ０，ｚ０，θｃ）であると仮定して、ピニオン歯面１０２Ｐ上に設定した基準点の３次元座標データ（ｒ１，ｚ１，θ１）と、収束点の３次元座標データ（ｒ０，ｚ０，θｃ）とから、以下の（１３），（１４）式を作成する。
0004823897.tif 000005
ここで、ピニオン歯面１０２Ｐ上の任意の点（ｊ，ｉ）の円筒座標系上での半径座標（Ｒ座標）成分ｒは、以下の（１５）式を用いて演算することが可能であり、座標ｒでのｉ方向への傾き（∂ｒ／∂ｉ）、及びｊ方向への傾き（∂ｒ／∂ｊ）は、以下の（１６），（１７）式を用いて演算することが可能である。【００３７】
0004823897.tif 000006
なお、（１５）〜（１７）式中において、ｆＲｊ０（ｉ）、ｆＲｊ１（ｉ）、ｆＲｊ０（ｊ）、ｆＲｉ１（ｊ）は、図９に示すように、点（ｊ，ｉ）を囲繞する関数上の任意点での半径座標成分であり、これらは上述のステップＳ１０８で作成した関数に基づいて演算される。Ｒ（ｊ０，ｉ０）、Ｒ（ｊ０，ｉ１）、Ｒ（ｊ１，ｉ０）、Ｒ（ｊ１，ｉ１）は点（ｊ，ｉ）を囲繞する各格子点での半径座標成分である。【００３８】  同様に、ピニオン歯面１０２Ｐ上の任意の点（ｊ，ｉ）の円筒座標系上での軸座標（Ｚ座標）成分ｚは、以下の（１８）式を用いて演算することが可能であり、座標ｚでのｉ方向への傾き（∂ｚ／∂ｉ）、及びｊ方向への傾き（∂ｚ／∂ｊ）は、以下の（１９），（２０）式を用いて演算することが可能である。【００３９】
0004823897.tif 000007
なお、（１８）〜（２０）式中において、ｆＺｊ０（ｉ）、ｆＺｊ１（ｉ）、ｆＺｉ０（ｊ）、ｆＺｉ１（ｊ）は、点（ｊ，ｉ）を囲繞する関数上の任意点での軸座標成分であり、これらは上述のステップＳ１０８で作成した関数に基づいて演算される。Ｚ（ｊ０，ｉ０）、Ｚ（ｊ０，ｉ１）、Ｚ（ｊ１，ｉ０）、Ｚ（ｉ１，ｊ１）は点（ｊ，ｉ）を囲繞する各格子点での軸座標成分である。【００４０】  そして、演算部６は、（１３）、（１４）式による連立方程式を解くことにより、基準点から収束点までの媒介変数の偏差Δｉ，Δｊを求める。【００４１】  続くステップＳ２０５において、演算部６は、ステップＳ２０４で求めた偏差Δｉ，Δｊを用い、以下の（２１）、（２２）式を用いて基準点の媒介変数（ｊ，ｉ）を更新する。
0004823897.tif 000008
【００４２】  そして、ステップＳ２０６に進むと、演算部６は、ステップＳ２０５で更新した媒介変数（ｊ，ｉ）を基に、上述の（１５）、（１８）式、及び以下の（２３）式を用いて基準点での３次元座標データ（ｒ１，ｚ１，θ１）を更新する。
0004823897.tif 000009
なお、（２３）式中において、ｆθｊ０（ｉ）、ｆθｊ１（ｉ）、ｆθｉ０（ｊ）、ｆθｉ１（ｊ）は、点（ｊ，ｉ）を囲繞する各点での角度座標成分であり、これらは上述のステップＳ１０８で作成した関数に基づいて演算される。θ（ｊ０，ｉ０）、θ（ｊ０，ｉ１）、θ（ｊ１，ｉ０）、θ（ｊ１，ｉ１）は点（ｊ，ｉ）を囲繞する各格子点での角度座標成分である。【００４３】  そして、ステップＳ２０７に進むと、演算部６は、ステップＳ２０６で演算した基準点のＲ軸座標成分ｒ１及びＺ軸座標成分ｚ１が収束点のＲ座標成分ｒｃ及びＺ軸座標成分ｚｃに収束しているか否か（すなわち、ｒ１及びｚ１とｒ０及びｚ０とが予め設定した範囲内でそれぞれ一致するか否か）を調べ、収束していないと判定した場合にはステップＳ２０８に進み、収束していると判定した場合にはステップＳ２０８に進む。【００４４】  ステップＳ２０７からステップＳ２０８に進むと、演算部６は、現在選択されているギヤ歯面１０２Ｇ上の格子点に対して行われたステップＳ２０４〜Ｓ２０６の演算回数が例えば１０回以上であるか否かを調べ、演算回数が１０回よりも少ない場合にはステップＳ２０４に戻り、演算回数が１０回以上である場合にはステップＳ２０９に進む。【００４５】  そして、ステップＳ２０７或いはステップＳ２０８からステップＳ２０９に進むと、演算部６は、相対角度情報として、現在のギヤグリッドｊＧ，ｉＧ、歯車回転ステップ数Ｍにおける歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を演算した後、ステップＳ２０３に戻る。ここで、ステップＳ２０７で基準点のＲ軸座標成分ｒ１及びＺ軸座標成分ｚ１が収束点のＲ座標成分ｒｃ及びＺ軸座標成分ｚｃに収束したと判定してステップＳ２０９に進んだ場合（すなわち、現在選択されているギヤ歯面１０２Ｇ上の格子点に対応する点がピニオン歯面１０２Ｐ内に存在する場合）、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）は、以下の（２４）式によって算出される。【００４６】
0004823897.tif 000010
一方、ステップＳ２０８からステップＳ２０９に進んだ場合、現在選択されているギヤ歯面１０２Ｇ上の格子点に対応する点がピニオン歯面１０２Ｐ外に存在することを示す判定用の角度値（例えば、θ（ｊＧ，ｉＧ，Ｍ）＝２０００）が設定される。【００４７】  なお、演算部６は、本サブルーチンにおいて、ステップＳ２０９からステップＳ２０３に戻る際には、例えば、上述のステップＳ２０２で設定した値を、再度、ピニオン歯面１０２Ｐ上の基準点の座標として設定してもよい。【００４８】  メインルーチンにおいて、ステップＳ１０９からステップＳ１１０に進むと、演算部６は、ステップＳ１０９で演算した各歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を基に、ギヤ歯面１０２Ｇを基準とする相対歯面データ（例えば、図１０参照）を作成し、続くステップＳ１１１において、歯車回転ステップＭをインクリメント（Ｍ＝Ｍ＋１）した後、ステップＳ１０４に戻る。【００４９】  また、ステップＳ１０４からステップＳ１１２に進むと、演算部６は、図５に示す包絡面演算サブルーチンのフローチャートに従って、ギヤ歯面１０２Ｇとピニオン歯面１０２Ｐの噛合開始から終了までの歯面間の相対的な隙間距離を示す包絡面の演算を行う。【００５０】  このサブルーチンがスタートすると、演算部６は、ステップＳ３０１において、歯車回転ステップＭ＝１をセットし、続くステップＳ３０２において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したか否かを調べる。【００５１】  そして、演算部６は、ステップＳ３０２において、歯車回転ステップ数Ｍが未だ歯車回転最大ステップ数Ｍｍａｘに達していないと判定した場合には、ステップＳ３０３に進む。一方、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達していると判定した場合には、演算部６は、ステップＳ３０７に進む。【００５２】  ステップＳ３０２からステップＳ３０３に進むと、演算部６は、現歯車回転ステップ数Ｍにおいて、ギヤ歯面１０２Ｇ上の全格子点での歯面間角度θ（ｊＧ，ｉＧ，Ｍ）に対して後述する歯面間角度最小値θＳｍｉｎ（Ｍ）の抽出計算が終了したか否かを調べる。【００５３】  そして、ステップＳ３０３において、全格子点での歯面間角度θ（ｊＧ，ｉＧ，Ｍ）に対して歯面間角度最小値θＳｍｉｎの抽出計算が終了していないと判定すると、演算部６は、ステップＳ３０４に進み、以下の（２５）式により、現歯車回転ステップ数Ｍでの歯面間角度最小値θＳｍｉｎ（Ｍ）を、現在選択されている格子点の歯面間角度θ（ｊＧ，ｉＧ，Ｍ）で適宜更新する。【００５４】
0004823897.tif 000011
さらに、演算部６は、現在選択されている格子点の歯面間角度θ（ｊＧ，ｉＧ，Ｍ）で歯面間角度最小値θＳｍｉｎ（Ｍ）が更新された場合には、現歯車回転ステップＭにおける相対歯面の最凸点の座標（ＰＯＣｊ（Ｍ），ＰＯＣｉ（Ｍ））を現格子点の（ｊＧ，ｉＧ）に更新した後、ステップＳ３０３に戻る。【００５５】  一方、ステップＳ３０３において、全格子点での歯面間角度θ（ｊＧ，ｉＧ，Ｍ）に対して歯面間角度最小値θＳｍｉｎの抽出計算が終了したと判定すると、演算部６は、ステップＳ３０５に進む。【００５６】  そして、ステップＳ３０５において、演算部６は、以下の（２６）式により、全歯車回転ステップでの歯面間角度最小値θｍｉｎｍｉｎを、現歯車回転ステップＭでの歯面間角度最小値θＳｍｉｎで適宜更新する。【００５７】
0004823897.tif 000012
さらに、演算部６は、現歯車回転ステップＭの歯面間角度最小値θＳｍｉｎで全歯車回転ステップの歯面間角度最小値θｍｉｎｍｉｎが更新された場合には、全歯車回転ステップにおける相対歯面の最凸点（すなわち、歯車対１００の最凸点）の座標（ＡＰＥＸｊ，ＡＰＥＸｉ）を現歯車回転ステップＭにおける最凸点の座標（ＰＯＣｊ（Ｍ），ＰＯＣｉ（Ｍ））に更新する。【００５８】  そして、ステップＳ３０５からステップＳ３０６に進むと、演算部６は、歯車回転数ステップＭを更新（Ｍ＝Ｍ＋１）した後、ステップＳ３０２に戻る。【００５９】  また、ステップＳ３０２からステップＳ３０７に進むと、演算部６は、歯車回転ステップ数Ｍ毎にそれぞれ演算された全歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を用い、ステップＳ３１２までの処理により、ギヤ歯面１０２Ｇとピニオン歯面１０２Ｐの噛合開始から終了までの相対的な隙間距離をギヤ歯面１０２Ｇ上の格子点（ｊＧ，ｉＧ）毎に角度値で示す相対歯面データθＥＯ（ｊＧ，ｉＧ）を演算する。【００６０】  具体的に説明すると、演算部６は、ステップＳ３０７において、歯車回転ステップ数Ｍ＝１をセットし、続くステップＳ３０８において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したか否かを調べる。【００６１】  そして、ステップＳ３０８において、歯車回転ステップ数Ｍが未だ歯車回転最大ステップ数Ｍｍａｘに達していないと判定した場合には、演算部６は、ステップＳ３０９に進む。一方、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達していると判定した場合には、演算部６は、ステップＳ３１３に進む。【００６２】  ステップＳ３０８からステップＳ３０９に進むと、演算部６は、現歯車回転ステップ数Ｍにおいて、ギヤ歯面１０２Ｇ上の全格子点（ｊＧ，ｉＧ）に対して、後述する相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了したか否かを調べる。【００６３】  そして、ステップＳ３０９において、全格子点（ｊＧ，ｉＧ）に対して相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了していないと判定すると、演算部６は、ステップＳ３１０に進み、以下の（２７）式により、歯面間角度最小値θｍｉｎｍｉｎを基準とする相対歯面データθＥＯＭを演算する。
0004823897.tif 000013
【００６４】  続くステップＳ３１１において、演算部６は、以下の（２８）式により、現在選択されている格子点（ｊＧ，ｉＧ）の相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）を用いて、対応する相対歯面データθＥＯ（ｊＧ，ｉＧ）を適宜更新した後、ステップＳ３０９に戻る。
0004823897.tif 000014
【００６５】  一方、ステップＳ３０９において、現歯車回転ステップ数Ｍでの全格子点（ｊＧ，ｉＧ）に対して相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了したと判定すると、演算部６は、ステップＳ３１２に進み、歯車回転ステップ数Ｍをインクリメント（Ｍ＝Ｍ＋１）した後、ステップＳ３０８に戻る。【００６６】  また、ステップＳ３０８からステップＳ３１３に進むと、演算部６は、以下の（２９）式により、ギヤ歯面１０２Ｇ上の格子点（ｊＧ，ｉＧ）毎の相対歯面データθＥＯ（ｊＧ，ｉＧ）を、距離情報（相対歯面データＥＯ（ｊＧ，ｉＧ）に変換した後、サブルーチンを抜け、メインルーチンに戻る。
0004823897.tif 000015
これにより、歯車回転ステップ数Ｍ毎の各相対歯面（例えば、図１１参照）を合成した包絡面（例えば、図１２参照）が生成される。【００６７】  そして、ステップＳ１１２からステップＳ１１３に進むと、演算部６は、ステップＳ１１２で生成した３次元の包絡面データを２次元の等高線データに変換し（例えば、図１３参照）、当該等高線データをディスプレイ装置１３等の出力部８を通じて出力した後、ルーチンを抜ける。【００６８】  このような実施形態によれば、ギヤ歯面１０２Ｇ上の各３次元座標データ（ｘＧｊｉ，ｙＧｊｉ，ｚＧｊｉ）とピニオン歯面１０２Ｐ上の各３次元座標データ（ｘＰｊｉ，ｙＰｊｉ，ｚＰｊｉ）とを歯車対１００の組立諸元を用いて所定の噛合回転位置で互いに関連づけてギヤ１０１Ｇを基準とする円筒座標系の３次元座標データ（ｒＧｊｉ，ｚＧｊｉ，θＧｊｉ）及び（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に変換するとともに、ピニオン歯面１０２Ｐ上に設定した２次元の媒介変数（ｊ，ｉ）を用い各３次元座標データ（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に基づいてピニオン歯面１０２Ｐ上の点を表す関数ｆＲ（ｊ，ｉ）、ｆＺ（ｊ，ｉ）、ｆΘ（ｊ，ｉ）を作成することにより、ギヤ歯面１０２Ｇ上の各点（格子点）に対応するピニオン歯面１０２Ｐ上の各点を示す媒介変数（ｊ，ｉ）を、関数ｆＲ（ｊ，ｉ）、ｆＺ（ｊ，ｉ）からニュートン法を用いて容易且つ精度よく演算することができる。そして、演算した各媒介変数（ｊ，ｉ）を用いて関数ｆΘ（ｊ，ｉ）から得られる角度情報θＰｊｉを基に、所定の噛合回転位置においてギヤ歯面１０２Ｇ上の点（格子点）と当該格子点に対応するピニオン歯面１０２Ｐ上の点との隙間を示す歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を求めることにより、歯車対の基準歯面に関する情報を指標とすることなく、実歯面の計測情報に基づいて精度のよい歯面解析を実現することができる。【００６９】  その際、歯車対１００の組立諸元を用いて互いに関連付けた各３次元座標データを、歯車対１００のデフレクション値で補正することにより、歯面解析の精度をより向上させることができる。【００７０】  また、ギヤ１０１Ｇとピニオン１０１Ｐの各噛合回転位置（各歯車回転ステップ数Ｍ）で演算される各格子点（ｊＧ，ｉＧ）での歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の中から最小の歯面間角度（歯面間角度最小値θｍｉｎｍｉｎ）を抽出し、抽出した歯面間角度最小値θｍｉｎｍｉｎを基準とし、各歯面間角度θ（ｊＧ，ｉＧ，Ｍ）に基づいて、ギヤ１０１Ｇとピニオン１０１Ｐの噛合開始から終了までの相対的な隙間距離を示す３次元の包絡面を演算することにより、ユーザに対し、歯面解析に有用な歯車対情報を提供することができる。【００７１】  また、歯面間角度最小値θｍｉｎｍｉｎに対応するギヤ歯面１０２Ｇ上の座標（ＡＰＥＸｊ，ＡＰＥＸｉ）を歯車対１００の噛合の最凸点として抽出することにより、ユーザに対し、歯面解析に有用な歯車対情報を提供することができる。【００７２】  さらに、演算した３次元の包絡面のデータを２次元の等高線のデータに変換することにより、ユーザに対し、歯面解析により有用な歯車対情報を提供することができる。【００７３】  次に、図１４乃至図１９は本発明の第２の実施形態に係わり���図１４は歯面間隙間演算サブルーチンを示すフローチャート、図１５はギヤ歯面上の格子点とピニオン歯面に設定した平面上の点との関係を示す説明図、図１６はギヤ歯面上の格子点に対応するピニオン歯面上の点が平面上に存在するときの条件を示す説明図、図１７は包絡面演算サブルーチンを示すフローチャート、図１８はギヤとピニオンの歯当り領域及び接触点の軌跡を示す説明図、図１９はギヤとピニオンとの間の伝達誤差を示す説明図である。なお、本実施形態においては、例えば、ニュートン法とスプライン関数を用いて行う歯面間隙間の演算方法に代えて、ピニオン歯面を近似する複数の平面を用いて歯面間隙間を演算する点が上述の第１の実施形態に対して主として異なる。また、本実施形態は、上述の第１の実施形態で行った包絡面演算サブルーチンにおいて、１歯噛合いに対する演算に加え、複数歯噛合いに対する演算を行うものである。その他、上述の第１の実施形態と同様の構成については、説明を省略する。【００７４】  図１４に示すサブルーチンがスタートすると、演算部６は、ステップＳ１２０１において、ピニオン歯面１０２Ｐの各格子点の座標データ（ｒＰｊｉ，ｚＰｊｉ，θＰｊｉ）に基づいて、ピニオン歯面１０２Ｐを近似する平面を作成する。すなわち、演算部６は、ステップＳ１２０１において、互いに近接する３点の格子点の組合せを順次抽出し、各組の格子点の座標データに基づいて例えば（３０）式に示す平面式を設定することにより、ピニオン歯面１０２Ｐを複数の平面で近似する（図１５参照）。
0004823897.tif 000016
ここで、（３０）式において、ａ，ｂ，ｃは係数であり、ｄは定数である。【００７５】  続くステップＳ１２０２において、演算部６は、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算をギヤ歯面１０２Ｇ上の全格子点ｊＧ，ｉＧに対して行ったか否かを調べる。そして、ギヤ歯面１０２Ｇ上の全格子点に対する歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の演算が終了していないと判定すると、演算部６は、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算対象とする格子点を新たな格子点に更新した後、ステップＳ１２０３に進み、以下のステップＳ１２０６までの処理により、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を演算する。一方、ステップＳ１２０２において、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算をギヤ歯面１０２Ｇ上の全格子点（ｊＧ，ｉＧ）に対して行ったと判定すると、演算部６は、サブルーチンを抜け、メインルーチンに戻る。【００７６】  ステップＳ１２０２からステップＳ１２０３に進むと、演算部６は、ステップＳ１２０１で作成した平面の中から、今回、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算対象として選択したギヤ歯面１０２Ｇ上の（ｊＧ，ｉＧ）番目の格子点のＲ座標成分及びＺ座標成分を含む平面を検索する。【００７７】  ここで、例えば、図１６に示すように、ピニオン歯面１０２Ｐ上のある平面を規定する３つの格子点Ａ，Ｂ，Ｃと、ギヤ歯面１０２Ｇ上の格子点Ｄと、をＲ−Ｚの２次元座標上で表した場合、格子点Ｄから格子点Ａまでのベクトル、格子点Ｄから格子点Ｂまでのベクトル、格子点Ｄから格子点Ｃまでのベクトルは、各格子点の基準点（０，０）から各点Ａ〜Ｄまでの各ベクトルを用いて以下の（３１）〜（３３）式で表される。なお、以下の説明において、〔  〕ｖは、ベクトルを示す。
0004823897.tif 000017
また、これら各ベクトルを用いると、〔ＤＡ〕ｖと〔ＤＢ〕ｖとのなす角θ１、〔ＤＢ〕ｖと〔ＤＣ〕ｖとのなす角θ２、〔ＤＣ〕ｖと〔ＤＡ〕ｖとのなす角θ３は、以下の（３４）〜（３６）式で表される。
0004823897.tif 000018
そして、Ｒ−Ｚの２次元座標系において、θ１＋θ２＋θ３＝２πであるとき、格子点Ｄは、格子点Ａ，Ｂ，Ｃで囲まれた領域内に存在する。【００７８】  そこで、演算部６は、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）の計算対象として選択したギヤ歯面１０２Ｇ上の（ｊＧ，ｉＧ）番目の格子点での３次元座標データを（ｒ０，ｚ０，θ０）とし、この格子点のＲ座標成分ｒ０及びＺ座標成分ｚ０と、各平面を規定するピニオン歯面１０２Ｐ上の各３個の格子点のＲ座標成分及びＺ座標成分とを用いた上述のベクトル演算を順次行うことにより、該当する平面を検索する。【００７９】  ステップＳ１２０３からステップＳ１２０４に進むと、演算部６は、Ｒ−Ｚの２次元座標系において、今回選択したギヤ歯面１０２Ｇ上の格子点を含む平面が存在したか否かを調べ、平面が存在したと判定した場合には、ステップＳ１２０５に進み、平面が存在しないと判定した場合にはステップＳ１２０６に進む。【００８０】  ステップＳ１２０４からステップＳ１２０５に進むと、演算部６は、Ｒ−Ｚ−Θの円筒座標系において、ギヤ歯面１０２Ｇ上の格子点と同一周上に存在するピニオン歯面１０２Ｐ上の点のΘ座標成分θ１を求めた後、ステップＳ１２０６に進む。ここで、ギヤ歯面１０２Ｇ上の格子点と同一周上に存在するピニオン歯面１０２Ｐ上の点のＲ座標成分とＺ座標成分は、ギヤ歯面１０２Ｇ上の格子点と一致する。そこで、演算部６は、ステップＳ１２０３で検索した平面の式を変形してｒ０及びｚ０を代入した以下の（３７）式によって、θ１を求める。
0004823897.tif 000019
そして、ステップＳ１２０４或いはステップＳ１２０５からステップＳ１２０６に進むと、演算部６は、相対角度情報として、現在のギヤグリッドｊＧ，ＩＧ、歯車回転ステップ数Ｍにおける歯面間角度θ（ｊＧ，ｉＧ，Ｍ）を演算した後、ステップＳ１２０２に戻る。ここで、ステップＳ１２０４で、今回選択したギヤ歯面１０２Ｇ上の格子点を含む平面が存在したと判定してステップＳ１２０５を経てステップＳ１２０６に進んだ場合、歯面間角度θ（ｊＧ，ｉＧ，Ｍ）は上述の（２４）式によって算出される。【００８１】  一方、ステップＳ１２０４から１２０６にそのまま進んだ場合、現在選択されているギヤ歯面１０２Ｇ上の格子点に対応する点がピニオン歯面１０２Ｐ外に存在することを示す判定用の角度値（例えば、θ（ｊＧ，ｉＧ，Ｍ）＝２０００）が設定される。【００８２】  次に、図１７に示すサブルーチンがスタートすると、演算部６は、ステップＳ１３０１〜ステップＳ１３０６において、上述の第１の実施形態で説明したステップＳ３０１〜ステップＳ３０６の処理と略同様の処理を行う。【００８３】  そして、ステップＳ１３０２からステップＳ１３０７に進むと、演算部６は、歯車回転ステップ数Ｍ＝１をセットし、続くステップＳ１３０８において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したか否かを調べる。【００８４】  そして、演算部６は、ステップＳ１３０８において、歯車回転ステップ数Ｍが未だ歯車回転最大ステップ数Ｍｍａｘに達していないと判定した場合には、ステップＳ１３０９に進む。一方、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達していると判定した場合には、演算部６はステップＳ１３１２に進む。【００８５】  ステップＳ１３０８からステップＳ１３０９に進むと、演算部６は、歯車回転ステップ数Ｍでの複数歯噛合いにおける歯面間角度最小値θＭｍｉｎ（Ｍ）を、以下の（３８）式を用いて求める。
0004823897.tif 000020
ここで、（３８）式中において、ＯＰＮは１噛合いピッチ当りの角度ステップ数である。
また、ｋは、現在注目する歯の前後に連なる歯の番号であり、以下の（３９）式を満たす範囲内の整数とする。
0004823897.tif 000021
なお、（３９）式中において、「ｉｎｔ」とは（  ）内の小数点以下を切り捨てるためのソフト関数である。【００８６】  そして、ステップＳ１３０９からステップＳ１３１０に進むと、演算部６は、現在の歯車回転ステップ数Ｍにおける１歯噛合いの伝達誤差ＭＣＳ（Ｍ）、及び、複数歯噛合いの伝達誤差ＭＣＭ（Ｍ）を、以下の（４０）式、及び、（４１）式を用いてそれぞれ求める。【００８７】
0004823897.tif 000022
そして、ステップＳ１３１１において、歯車回転ステップ数Ｍをインクリメント（Ｍ＝Ｍ＋１）した後、ステップＳ１３０８に戻る。【００８８】  また、ステップＳ１３０８からステップＳ１３１２に進むと、演算部６は、上述の処理で求めた各歯車回転ステップ数Ｍの伝達誤差ＭＣＭ（Ｍ）を用い、以下の（４２）式から、複数歯噛合いにおけるギヤ歯面１０２Ｇとピニオン歯面１０２Ｐとの噛合い開始から終了までの伝達誤差の変動幅ＭＥを演算する。【００８９】
0004823897.tif 000023
これにより、例えば、図１９に示すように、複数歯噛合いによるギヤ歯面１０２Ｇとピニオン歯面１０２Ｐとの伝達誤差の特性図を得ることができる。なお、１歯噛合いによる伝達誤差の特性を図１９中に一点鎖線で示す。【００９０】  また、ステップＳ１３１２からステップＳ１３１３に進むと、演算部６は、ステップＳ１３１８までの処理により、上述の第１の実施形態で説明したステップＳ３０７〜ステップＳ３１２の処理と略同様の処理により、歯数比で回転する１対のギヤ歯面１０２Ｇとピニオン歯面１０２Ｐの噛合い開始から終了までの相対的な隙間距離をギヤ歯面１０２Ｇ上の格子点（ｊＧ，ｉＧ）毎に角度値で示す相対歯面データθＥＯ（ｊＧ，ｉＧ）を演算する。これに併せて、演算部６は、複数歯噛合いによる相対歯面データθＭＣＰ（ｊＧ，ｉＧ）を演算する。【００９１】  具体的に説明すると、演算部６は、ステップＳ１３１３において、歯車回転ステップ数Ｍ＝１をセットし、続くステップＳ１３１４において、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達したか否かを調べる。【００９２】  そして、ステップＳ１３１４において、歯車回転ステップ数Ｍが未だ歯車回転最大ステップ数Ｍｍａｘに達していないと判定した場合には、演算部６は、ステップＳ１３１５に進む。一方、歯車回転ステップ数Ｍが歯車回転最大ステップ数Ｍｍａｘに達していると判定した場合には、演算部６は、ステップＳ１３１９に進む。【００９３】  ステップＳ１３１４からステップＳ１３１５に進むと、演算部６は、現歯車回転ステップ数Ｍにおいて、ギヤ歯面１０２Ｇ上の全格子点（ｊＧ，ｉＧ）に対して、１歯噛合いによる相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）、及び、複数歯噛合いによる相対歯面データθＭＣＰＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了したか否かを調べる。【００９４】  そして、ステップＳ１３１５において、全格子点（ｊＧ，ｉＧ）に対して相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）、θＭＣＰＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了していないと判定すると、演算部６は、ステップＳ１３１６に進み、上述の（２７）式により、歯面間角度の最小値θｍｉｎｍｉｎを基準とする相対歯面データθＥＯＭ、および以下の（４３）式により、複数歯噛合いによるステップ毎歯面間角度の最小値θＭｍｉｎ（Ｍ）を基準とする相対歯面データθＭＣＰＭを演算する。
0004823897.tif 000024
【００９５】  続くステップＳ１３１７において、演算部６は、上述の（２８）式、及び、以下の（４４）式により、現在選択されている格子点（ｊＧ，ｉＧ）の相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）、θＭＣＰＭ（ｊＧ，ｉＧ，Ｍ）を用いて、対応する相対歯面データθＥＯ（ｊＧ，ｉＧ）、θＭＣＰ（ｊＧ，ｉＧ）を適宜更新した後、ステップＳ１３１５に戻る。
0004823897.tif 000025
【００９６】  一方、ステップＳ１３１５において、現歯車回転ステップ数Ｍでの全格子点（ｊＧ，ｉＧ）に対して相対歯面データθＥＯＭ（ｊＧ，ｉＧ，Ｍ）、θＭＣＰＭ（ｊＧ，ｉＧ，Ｍ）の計算が終了したと判定すると、演算部６は、ステップＳ１３１８に進み、歯車回転ステップ数Ｍをインクリメント（Ｍ＝Ｍ＋１）した後、ステップＳ１３１４に戻る。【００９７】  また、ステップＳ１３１４からステップＳ１３１９に進むと、演算部６は、上述の（２９）式、及び、以下の（４５）式により、１歯噛合い及び複数歯噛合いによるギヤ歯面１０２Ｇ上の格子点（ｊＧ，ｉＧ）毎の相対歯面データθＥＯ（ｊＧ，ｉＧ）、θＭＣＰ（ｊＧ，ｉＧ）を、距離情報（噛合い全域での相対歯面データＥＯ（ｊＧ，ｉＧ）、歯当り分布データＭＣＰ（ｊＧ，ｉＧ））に変換した後、サブルーチンを抜け、メインルーチンに戻る。
0004823897.tif 000026
これにより、歯数比で回転する１対の歯面および、実働回転角で運動する複数歯歯面の噛合いによる、歯車回転ステップ数Ｍ毎の各相対歯面を合成した包絡面が生成される。そして、複数歯噛合いによる包絡面から、例えば、図１８に示すように、ギヤとピニオンの歯当り領域及び接触点の軌跡を得ることができる。【００９８】  このような実施形態によれば、上述の第１の実施形態でえら得る効果に加え、歯面間隙間演算の演算量を飛躍的に軽減することができるという効果を奏する。【００９９】  また、ギヤとピニオンの歯当り領域及び接触点の軌跡、伝達誤差等を演算することにより、より実用的な歯車対の評価を行うことができる。【０１００】  なお、歯面間隙間演算等の演算方法は、上述の各実施形態で示したものに限定されるものではなく、円筒座標系を用いた他の演算方法を用いてもよいことは勿論である。また、上述の各実施形態の構成を適宜組み合わせてもよいことは勿論である。
【図１】
0004823897.tif 000027
【図２】
0004823897.tif 000028
【図３】
0004823897.tif 000029
【図４】
0004823897.tif 000030
【図５】
0004823897.tif 000031
【図６】
0004823897.tif 000032
【図７】
0004823897.tif 000033
【図８】
0004823897.tif 000034
【図９】
0004823897.tif 000035
【図１０】
0004823897.tif 000036
【図１１】
0004823897.tif 000037
【図１２】
0004823897.tif 000038
【図１３】
0004823897.tif 000039
【図１４】
0004823897.tif 000040
【図１５】
0004823897.tif 000041
【図１６】
0004823897.tif 000042
【図１７】
0004823897.tif 000043
【図１８】
0004823897.tif 000044
【図１９】
0004823897.tif 000045