JP2005208735A - 歯車形状のシミュレーション方法 - Google Patents
歯車形状のシミュレーション方法 Download PDFInfo
- Publication number
- JP2005208735A JP2005208735A JP2004011953A JP2004011953A JP2005208735A JP 2005208735 A JP2005208735 A JP 2005208735A JP 2004011953 A JP2004011953 A JP 2004011953A JP 2004011953 A JP2004011953 A JP 2004011953A JP 2005208735 A JP2005208735 A JP 2005208735A
- Authority
- JP
- Japan
- Prior art keywords
- tooth surface
- gear
- simulation
- shape
- tooth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Abstract
【課題】 歯車形状のシミュレーション前後における形状差が、実際の歯車の焼き入れ前後を測定した形状差とよく一致する歯車形状のシミュレーション方法を提供する。
【解決手段】 歯車形状モデルを作成し(S1)、焼き入れのシミュレーションを実行して(S3)、その結果から計算後歯面モデルを作成する(S4、5)と共に、シミュレーションを実行しない標準歯面モデルを作成する(S14、15)。計算後歯面モデルと標準歯面モデルの歯面上の重心位置を位置合わせして(S6)、その断面から歯面形状の差を求める(S7、8)。
【選択図】図2
【解決手段】 歯車形状モデルを作成し(S1)、焼き入れのシミュレーションを実行して(S3)、その結果から計算後歯面モデルを作成する(S4、5)と共に、シミュレーションを実行しない標準歯面モデルを作成する(S14、15)。計算後歯面モデルと標準歯面モデルの歯面上の重心位置を位置合わせして(S6)、その断面から歯面形状の差を求める(S7、8)。
【選択図】図2
Description
本発明は、歯車形状のシミュレーション方法に関する。
歯車形状のシミュレーションとして、有限要素法を用いたシミュレーション方法が数多く開発されている。
このような歯車形状のシミュレーションの一つとして、従来、歯車の焼き入れに起因する形状変化をシミュレーションする技術がある。
一例を挙げれば、焼き入れ前の歯車の着目する歯面を基準歯面として、焼き入れ後のシミュレーションされた歯車形状を、回転対象中心を一致させて、さらに、歯面測定領域端の一点の位置合わせを空間並進にて行い、この位置合わせ後の測定領域の対辺上の点を結んだ直線の勾配量を測定することで、焼き入れによる歯面の圧力角変化量をシミュレーションにより求める技術がある(特許文献1参照)。
特開2003−194754号公報
このような従来のシミュレーション方法では、基準歯面において基準点とも言うべき位置合わせを空間並進により行って、焼き入れ前の形状と、シミュレーションによる焼き入れ後の形状との比較を行っている。
しかしながら、実際の歯車では、歯車の測定器が、実際の歯車を、回転軸中心に回転させるとともに回転軸方向への並進により位置合わせを行って、焼き入れ前後の形状測定を行っているため、上記従来のシミュレーションにより得られた焼き入れ前後の形状差と実際の歯車を測定した焼き入れ前後の形状差がうまく一致しないと言う問題があった。
そこで本発明の目的は、歯車形状のシミュレーション前後における形状差が、実際の歯車を測定した結果による形状差とよく一致するようになる歯車形状のシミュレーション方法を提供することである。
上記目的を達成するための本発明は、歯車の歯面形状からなる標準歯面形状のデータを用いて熱処理後の形状を有限要素法によるシミュレーションによって算出する段階と、前記シミュレーションによって算出された算出後歯面形状と前記標準歯面形状のそれぞれの歯面上における重心位置の位置合わせを行う段階と、重心位置を合わせた前記算出後歯面形状と前記標準歯面形状の差を求める段階と、を有することを特徴とする歯車形状のシミュレーション方法である。
本発明によれば、シミュレーション前後の歯面の重心位置を合わせて、それらの誤差を求めることとしたので、実際の歯車を測定した結果と良く一致する歯車の形状差を得ることができるようになる。
以下、図面を参照して本発明を実施するための最良の形態を説明する。
図1は、本発明による歯車形状のシミュレーションを行うためのシステム構成を示す図面である。
このシステムは、シミュレーション対象とする歯車の形状データを入力する形状データ入力部1と、歯車の形状データと入力された諸条件から有限要素法によるシミュレーションを実行するシミュレーション部2と、シミュレーション部2に接続されたディスプレイ3および記憶装置4とからなる。
形状データ入力部1は、シミュレーションを行う歯車の形状データとして、たとえば、歯車の設計形状を用いる場合にはCAD装置からデータを入力し、また、実際に歯車の形状を測定したデータ(実測値)を用いる場合には、実測された歯車の形状データが入力される。
シミュレーション部2は、歯車の形状が熱処理が加わることによりどのように形状変化するか算出する部分であり、実質的には有限要素法を用いたシミュレーションを実行するためのシミュレーション装置である。このようなシミュレーション装置は、周知の通り、形状データ入力部1から入力された形状データと、シミュレーションの諸条件を設定するプリプロセッサ21と、有限要素法による数値解析を用いて形状変化を解析するソルバ22と、計算結果をビジュアルに表示するためのポストプロセッサ23とからなる。
なお、このようなシミュレーション装置は実際には、コンピュータに導入されたプログラムによって実行されるもので、このようなプログラムとしては、たとえば、SFTC社製のDEORM−HTなどが利用できる。
ディスプレイ3は、必要に応じて形状モデルの3次元画像および後述する歯形誤差の数値などが表示される。
記憶装置4は、入力された形状データや歯形誤差の算出結果などが記憶される。
また、このシステムには、図示しないプリンタなどが接続されて、必要に応じて形状モデルの3次元画像および後述する歯形誤差の数値などが印刷されるようにしてもよい。
次に、このシステムを用いた歯車形状のシミュレーション方法について説明する。図2は、歯車形状のシミュレーション方法の手順を示すフローチャートである。なお、この手順は、上記システムにおけるシミュレーション部2によって実行される手順を示すものである。
まず、形状データ入力部1から入力された歯車の形状データをもとに、シミュレーションを実行する歯車の形状モデルを作成する(S1)。ここで作成される形状モデルは、湯減要素法の解析に用いるメッシュモデルであり、プリプロセッサ21により作成される。
図3は、一例としてシミュレーションを実行する歯車の形状を示す図面である。図示する歯車101は、ハイポイド歯車である。
この図3に示す歯車101の全体形状データから、ここでは説明をわかりやすくするために、シミュレーションを実行する一つの歯面について説明する。図4は、この歯車101の一つの歯分の形状102を示す図面である。ここでは、歯車の軸部分102aも示しているデータである。このデータから、図5に示すように、一つの歯分のメッシュモデル103を形成する。
続いて、プリプロセッサ21から境界条件の設定が行われる(S2)。
境界条件はあらかじめ入力されたもので、ここでは、焼き入れの前後における歯面形状の変化をシミュレーションするものであるため、主に、歯車101の熱処理に関係する条件を設定することになる。このような熱処理にかかわる条件としては、たとえば、歯車の素材、その素材および/または焼き入れ環境の熱伝達係数、焼き入れの温度(熱処理温度)などである。ただし、このような境界条件は、使用するシミュレーションプログラムによって異なるので、それらに合わせて条件を設定することになる。
なお、本実施形態では、この境界条件として、後述するように、実際に行われている歯車の焼き入れとそれに伴う歯形誤差など実測値を用いて決定している。
続いて、形成された形状モデルと入力された境界条件などからソルバ22によってシミュレーションが実行される(S3)。
続いて、ポストプロセッサ23において、シミュレーション後の歯面形状から節点情報を抽出し、歯面のスプライン補完を行う(S4)。スプライン補間は、抽出した接点情報から、歯面に沿ってベジェ曲線を求めてその曲線に載る節点を結んだ線としている。
図6に、このようにして抽出されたスプライン曲線104によって表示された歯面を示す。
続いて、このようなスプライン曲線から歯面形状モデルを作成する(S5)。図7は、作成した歯面形状モデルを示す図面である。この歯面形状モデルを計算後歯面モデル110と称する。
同様にして、ポストプロセッサ23を用いて、シミュレーションを行わない歯面の形状モデルも作成する(S14、S15)。シミュレーションを行わない歯面形状モデルを標準歯面モデルと称する。
ここで、歯面形状モデルを求めているのは、そもそも本シミュレーションの目的が、歯車の焼き入れに伴う形状変化したかを求めるものであるので、このような形状変化は歯車の歯面の位置変化として求めればよい。このため、焼き入れ(熱処理)による形状変化のシミュレーション自体は歯面の立体形状がどのように変形するかを算出するが、実際にその変化量を見るためには、歯面形状のみあればよい。したがって、本実施形態では、上記ステップS4、S5、同様にS14、S15において、歯面形状モデルを形成しているものである。
続いて、ポストプロセッサ23により、計算後歯面モデルと標準歯面モデルの位置合わせを行う(S6)。
この位置合わせは、標準歯面モデルに対して、計算後歯面モデルの歯面上の重心位置を合わせるようにする。なお、歯面上の重心位置とは、たとえば、図8に示すように、歯面モデル100上にメッシュを展開し(実際にはシミュレーションで使用したメッシュモデルのメッシュをそのまま利用してもよい)、そのメッシュ個数が歯面の端から縦、横ともに50:50となる位置を歯面上の重心とするとよい。
図9は、この重心の位置合わせを説明するための図面である。
重心位置合わせは、歯車101の軸中心A(図1参照)を中心に歯車101全体を回転させると共に、軸そのものを上下(図1における矢印Z方向)に移動させ、さらに歯面そのものを放射方向(歯車の軸中心から外側方向またはその逆方向、図9における矢印Y方向)に移動させて、計算後歯面モデル110と標準歯面モデル111の重心位置Oを合わせる。なお、この操作は、ポストプロセッサ23によって表示されている計算後歯面モデル110と標準歯面モデル111の3次元表示図形を手動により上記の通り移動させて重心の位置合わせを行ってもよいし、上記のような移動をあらかじめプログラムして、重心位置が合うように自動的に行ってもよい。
続いて、位置合わせを行った2つの歯面モデル(すなわち、計算後歯面モデル110と標準歯面モデル111)について、図10に示すように、重心Oを通る歯面内側から外側方向(この歯車101での圧力角方向)の断面(図示B−B断面)と、同じく重心Oを通り圧力角に直交する方向(ねじれ角方向)の断面(図示C−C断面)を作成する(S7)
続いて、これら2つの歯面モデルの断面における差を算出し、算出された値を誤差曲線を作成する(S8)。この誤差曲線は、互いの歯面の法線距離の値である。ここで得られた誤差曲線が、焼き入れによる形状変化を示すものとなる。
続いて、これら2つの歯面モデルの断面における差を算出し、算出された値を誤差曲線を作成する(S8)。この誤差曲線は、互いの歯面の法線距離の値である。ここで得られた誤差曲線が、焼き入れによる形状変化を示すものとなる。
その後は、必要に応じて誤差曲線を数値化して出力する(S9)。このときディスプレイ上に誤差曲線そのものと数値をあわせて出力してもよいし、またプリンタにより印刷してもよい。その後、処理を終了する。
図11は、圧力角誤差(図10におけるB−B線断面)における誤差曲線の一例を示す図面である。このようにして得られた誤差曲線は、実際の歯車において計測した焼き入れ前後の形状差と良く一致したものとなる。
次に、上記シミュレーションにおいて用いる境界条件の設定方法について説明する。
図12は、シミュレーションの境界条件の設定手順を示す図面である。
歯形誤差の触れ幅を管理項目として使用する場合について説明する。
この歯形誤差の触れ幅は実際の歯車製造工程における管理項目の中から、歯形誤差に効いてくる要因を探す出すことにより、それを境界条件として設定すれば、実際の工程に合わせたシミュレーションを実行することができる。
これには、まず、工程管理データを取得する(S21)。工程管理データは、たとえば、図13に示すように、実際の製造工程における歯車の形状測定日、品番(同一歯車は同一番号)、チャージ(材料に添加する成分)、歯筋誤差、歯形誤差、焼き入れ保持温度、などのデータである。
続いて、周知の決定木手法により、歯形誤差に最も効いてくる要因を検索する(S22)。
図14は、決定木手法により、歯形誤差に最も効いてくる要因(項目)を検索する例を示す図面である。ここでは、要因として焼き入れ保持温度とチャージについて決定木によりどの程度不良品が発生したかにより、焼き入れ保持温度とチャージのいずれがより歯車の歯形誤差に対して効いているかを決定する。なお、図においては、良品の場合をOK、不良品発生の場合をNGと記している。
続いて、所得した工程管理データの中から、決定木で決定した最も歯形誤差に効く要因の範囲を抽出し(S23)、処理終了となる。そしてここで抽出した最も歯形誤差に効く要因の範囲を前述のシミュレーションにおける境界条件として用いる。
図15は、決定木で決定した最も歯形誤差に効く要因を、たとえば焼き入れ保持温度となった場合における焼き入れ保持温度の変化に対する歯形誤差の範囲を示す図面である。
図示するように、焼き入れ保持温度の変化範囲が約838〜865℃で歯形誤差が振れている。したがって、この焼入れ保持温度の変化範囲を前述した境界条件設定の際の熱伝達係数の一要因の振り幅の範囲として用いることで、シミュレーションを行う際に、実際に発生する歯形誤差をシミュレーションすることが可能となる。
これにより、これまでたとえばシミュレーションを実行する技術者の経験や勘に頼っていたシミュレーション条件の設定を容易にし、しかも、実際に起こり得る範囲で設定することができるようになり、シミュレーションの無駄を減らすことができる。
以上説明したように、本実施の形態によれば、シミュレーション前後の歯面の形状モデルにおける歯面の重心位置を合わせて、それらの誤差を求めることとしたので、実際の歯車を測定した結果と良く一致する歯形誤差を得ることができるようになる。また、歯面の重心を合わせる際には、歯車101の軸を中心とした回転、軸方向への移動、および歯車101の放射方向へ移動だけであるので、比較的簡単に位置合わせを行うことが可能で、また、計算機により行う場合でも移動させる方向が少ないので計算機負荷を小さくすることができる。
また、本実施の形態では、実際の工程管理データを元にシミュレーションに用いる境界条件を設定しているので、実際に起こり得る歯形誤差をシミュレーションすることができる。したがって、より実際に近い歯形誤差の品質変動をシミュレーションにより見ることができる。
算出後歯面形状と標準歯面形状は、有限要素法に用いるメッシュモデルの歯面上の節点をベジェ曲線を利用してスプラインを描いて作成しているので、ベジェ曲線の使用により変化の大きな歯面の曲面を誤差が少なく表すことができる。また、ベジェ曲線のような関数を用いることで、複雑な歯面も容易に作成することができる。
なお、以上説明した発明を実施するための最良の形態は、あくまでも本発明を用いた一実施形態に過ぎず、本発明はこのような形態に限定されるものではない。たとえば、シミュレーションを行う歯車は、当然、本形態で示した歯車に限定されず、さまざまな形状の歯車において実施可能である。また、シミュレーションを行う際の境界条件は、工程管理データから得ているが、これについてもシミュレーションを行う際に、任意に設定してもよい。その他、当業者においてさまざまな改変が行えることはいうまでもない。
本発明は、歯車101の焼き入れの前後における形状変化のシミュレーション、歯車101を使用することにより熱履歴が加わった場合の形状変化のシミュレーションなどに使用することができる。
1…形状データ入力部、
2…シミュレーション部、
3…ディスプレイ、
4…記憶装置、
21…プリプロセッサ、
22…ソルバ、
23…ポストプロセッサ、
101…歯車、
102…一歯分の形状、
103…メッシュモデル、
110…計算後歯面モデル、
111…標準歯面モデル。
2…シミュレーション部、
3…ディスプレイ、
4…記憶装置、
21…プリプロセッサ、
22…ソルバ、
23…ポストプロセッサ、
101…歯車、
102…一歯分の形状、
103…メッシュモデル、
110…計算後歯面モデル、
111…標準歯面モデル。
Claims (4)
- 歯車の歯面形状からなる標準歯面形状のデータを用いて熱処理後の形状を有限要素法によるシミュレーションによって算出する段階と、
前記シミュレーションによって算出された算出後歯面形状と前記標準歯面形状のそれぞれの歯面上における重心位置の位置合わせを行う段階と、
重心位置を合わせた前記算出後歯面形状と前記標準歯面形状の差を求める段階と、
を有することを特徴とする歯車形状のシミュレーション方法。 - 前記位置合わせを行う段階は、前記標準歯面形状に対して前記算出後歯面形状を、歯車の軸を中心にした回転、前記軸の上下移動、および歯車の放射方向への歯面の移動により行うことを特徴とする請求項1記載の歯車形状のシミュレーション方法。
- 前記算出後歯面形状と前記標準歯面形状は、有限要素法に用いるメッシュモデルの歯面上の節点をベジェ曲線を利用して接続することにより形成していることを特徴とする請求項1または2記載の歯車形状のシミュレーション方法。
- 前記シミュレーションは、実際に製造された歯車の歯形誤差の測定結果から、当該歯形誤差に起因する要因を抽出して、当該要因の範囲を前記シミュレーションにおける条件範囲とすることを特徴とする請求項1〜4のいずれか一つに記載の歯車形状のシミュレーション方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004011953A JP2005208735A (ja) | 2004-01-20 | 2004-01-20 | 歯車形状のシミュレーション方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004011953A JP2005208735A (ja) | 2004-01-20 | 2004-01-20 | 歯車形状のシミュレーション方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005208735A true JP2005208735A (ja) | 2005-08-04 |
Family
ID=34898483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004011953A Pending JP2005208735A (ja) | 2004-01-20 | 2004-01-20 | 歯車形状のシミュレーション方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005208735A (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2012275A1 (en) | 2007-07-06 | 2009-01-07 | Toyota Jidosha Kabusiki Kaisha | Design support system, method and program |
CN106021721A (zh) * | 2016-05-19 | 2016-10-12 | 山东科技大学 | 一种渗碳圆柱齿轮参数化仿真分析方法及cae系统 |
CN109033669A (zh) * | 2018-08-08 | 2018-12-18 | 中南大学 | 基于万能运动参数驱动的螺旋锥齿轮仿真加工建模方法 |
KR20200036961A (ko) * | 2018-09-28 | 2020-04-08 | 주식회사 태웅 | 대형 조미니 시험편을 이용한 구조 안정성 예측방법 |
CN111390299A (zh) * | 2020-03-25 | 2020-07-10 | 中国北方车辆研究所 | 一种浮动支撑摩擦片加工变形预测方法 |
CN114818181A (zh) * | 2022-04-21 | 2022-07-29 | 哈尔滨工程大学 | 一种基于齿廓线的直齿圆柱齿轮有限元网格自动生成方法和计算机设备 |
-
2004
- 2004-01-20 JP JP2004011953A patent/JP2005208735A/ja active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2012275A1 (en) | 2007-07-06 | 2009-01-07 | Toyota Jidosha Kabusiki Kaisha | Design support system, method and program |
US8274528B2 (en) | 2007-07-06 | 2012-09-25 | Toyota Jidosha Kabushiki Kaisha | Support system, method and program |
CN106021721A (zh) * | 2016-05-19 | 2016-10-12 | 山东科技大学 | 一种渗碳圆柱齿轮参数化仿真分析方法及cae系统 |
CN109033669A (zh) * | 2018-08-08 | 2018-12-18 | 中南大学 | 基于万能运动参数驱动的螺旋锥齿轮仿真加工建模方法 |
CN109033669B (zh) * | 2018-08-08 | 2022-12-27 | 中南大学 | 基于万能运动参数驱动的螺旋锥齿轮仿真加工建模方法 |
KR20200036961A (ko) * | 2018-09-28 | 2020-04-08 | 주식회사 태웅 | 대형 조미니 시험편을 이용한 구조 안정성 예측방법 |
KR102117347B1 (ko) * | 2018-09-28 | 2020-06-02 | 주식회사 태웅 | 대형 조미니 시험편을 이용한 구조 안정성 예측방법 |
CN111390299A (zh) * | 2020-03-25 | 2020-07-10 | 中国北方车辆研究所 | 一种浮动支撑摩擦片加工变形预测方法 |
CN111390299B (zh) * | 2020-03-25 | 2021-03-05 | 中国北方车辆研究所 | 一种浮动支撑摩擦片加工变形预测方法 |
CN114818181A (zh) * | 2022-04-21 | 2022-07-29 | 哈尔滨工程大学 | 一种基于齿廓线的直齿圆柱齿轮有限元网格自动生成方法和计算机设备 |
CN114818181B (zh) * | 2022-04-21 | 2023-08-25 | 哈尔滨工程大学 | 基于齿廓线直齿圆柱齿轮有限元网格自动生成方法和设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dilberoglu et al. | Shrinkage compensation approach proposed for ABS material in FDM process | |
Leser et al. | A digital twin feasibility study (Part II): Non-deterministic predictions of fatigue life using in-situ diagnostics and prognostics | |
KR101355839B1 (ko) | 덴드라이트 암 스페이싱 및 기공도 기반 모델을 이용한 피로 수명 계산을 위한 재료 물성 분포 결정 | |
Kaveh et al. | Optimization of the printing parameters affecting dimensional accuracy and internal cavity for HIPS material used in fused deposition modeling processes | |
CN111274671A (zh) | 一种基于数字孪生的复杂产品装配过程精准修配方法及其运行系统 | |
Vivet et al. | On the modelling of gear alignment errors in the tooth contact analysis of spiral bevel gears | |
Kruth et al. | Uncertainty determination for CMMs by Monte Carlo simulation integrating feature form deviations | |
EP1060037A1 (en) | Formation method and device for curved plates | |
US10318675B2 (en) | Post-processing system for finite element analysis | |
Haefner et al. | Function-oriented measurements and uncertainty evaluation of micro-gears for lifetime prognosis | |
US11982998B2 (en) | Manufacturing monitoring assistance device, manufacturing monitoring assistance method, and manufacturing monitoring assistance program | |
Jaramillo et al. | Fast dimensional inspection of deformable parts from partial views | |
Lepore et al. | A FEM based methodology to simulate multiple crack propagation in friction stir welds | |
Forte et al. | A CAE approach for the stress analysis of gear models by 3D digital photoelasticity | |
EP3581318A1 (en) | System, method and control unit for friction stir welding apparatus | |
JP2005208735A (ja) | 歯車形状のシミュレーション方法 | |
Li et al. | Error measurements for flank milling | |
Hashemian et al. | Surface fairness: a quality metric for aesthetic assessment of compliant automotive bodies | |
Zheng et al. | Physics-informed deep Monte Carlo quantile regression method for interval multilevel Bayesian Network-based satellite circuit board reliability analysis | |
US10242921B2 (en) | Method of forming pattern of semiconductor device from which various types of pattern defects are removed | |
Li et al. | A novel variable selection approach for redundant information elimination purpose of process control | |
Khodaygan | Manufacturing error compensation based on cutting tool location correction in machining processes | |
Marciniec et al. | Comparative analysis of numerical methods for the determination of contact pattern of spiral bevel gears | |
Ren et al. | Stochastic modeling and diagnosis of leak areas for surface assembly | |
Lipman et al. | Testing implementations of geometric dimensioning and tolerancing in CAD software |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061127 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091104 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091224 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100209 |