JP5046671B2 - 共形接触下の接触面圧および表面下応力計算方法 - Google Patents
共形接触下の接触面圧および表面下応力計算方法 Download PDFInfo
- Publication number
- JP5046671B2 JP5046671B2 JP2007026577A JP2007026577A JP5046671B2 JP 5046671 B2 JP5046671 B2 JP 5046671B2 JP 2007026577 A JP2007026577 A JP 2007026577A JP 2007026577 A JP2007026577 A JP 2007026577A JP 5046671 B2 JP5046671 B2 JP 5046671B2
- Authority
- JP
- Japan
- Prior art keywords
- contact
- surface pressure
- calculation
- cell
- calculating
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims description 122
- 238000000034 method Methods 0.000 claims description 125
- 230000008569 process Effects 0.000 claims description 62
- 238000006073 displacement reaction Methods 0.000 claims description 43
- 238000004458 analytical method Methods 0.000 claims description 37
- 230000009471 action Effects 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 9
- 238000005096 rolling process Methods 0.000 claims description 8
- 238000013459 approach Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims 2
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 description 38
- 230000006870 function Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012887 quadratic function Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Description
この発明は、接触領域が平坦と見なせない接触物体間の接触面圧および表面下応力の計算を、有限要素法や境界要素法を使用することなく行う計算方法および装置に関する。
接触面圧の計算問題は下記の二つの分野に分けられる。
一つは非共形接触問題（Counter-formal（またはNon-Conforming）Problem ）であり、接触領域が接触部の物体半径に比べて小さい場合の問題である。
他の一つは、共形接触問題（Conformal （またはConforming）Problem ）であり、接触領域が接触部の物体半径に比べて大きい場合の問題である。
一つは非共形接触問題（Counter-formal（またはNon-Conforming）Problem ）であり、接触領域が接触部の物体半径に比べて小さい場合の問題である。
他の一つは、共形接触問題（Conformal （またはConforming）Problem ）であり、接触領域が接触部の物体半径に比べて大きい場合の問題である。
非共形接触問題（Counter-formal Problem）では、接触領域はほぼ平面上にあると見なせるため、比較的解き易い。
特に接触面が二次関数で表される曲面の場合は、ヘルツ（Hertz ）接触理論などが広く用いられており、数表を用いて電卓レベルで計算が可能である。二次関数で表すことができない曲面の場合は、半無限体に集中荷重が作用する場合の解析解（ブジネスク（Boussinesq）の解や、セルッティ（Cerruti ）の解）の重ね合わせにより、数値的に面圧分布を求める手法が知られている。また、これらと類似の解法で各種文献（例えば非特許文献１や非特許文献２）が出ており、この手法をベースとした面圧計算ソフトウェアも市販されている。
特開昭５９−６９５１９号公報
ハーネット（Hartnett, M. J. ）著 ,「転がり要素軸受の接触応力解析」（"The Analysis of Contact Stresses in Rolling Element Bearings"）, ASME J. of Lubr., 101, Jan.p105,(1979),
栫井著, 「三次元弾性接触問題の境界要素解析」, トライボロジスト, 36,6,(1991), 429-433
特に接触面が二次関数で表される曲面の場合は、ヘルツ（Hertz ）接触理論などが広く用いられており、数表を用いて電卓レベルで計算が可能である。二次関数で表すことができない曲面の場合は、半無限体に集中荷重が作用する場合の解析解（ブジネスク（Boussinesq）の解や、セルッティ（Cerruti ）の解）の重ね合わせにより、数値的に面圧分布を求める手法が知られている。また、これらと類似の解法で各種文献（例えば非特許文献１や非特許文献２）が出ており、この手法をベースとした面圧計算ソフトウェアも市販されている。
しかし、共形接触問題（Conformal Problem ）については、上記のような解析手法は適用できないため、有限要素法や境界要素法を用いて解く以外に計算方法がない。
有限要素法や境界要素法を使用して計算する場合、計算時間が多大になるという課題がある。
有限要素法や境界要素法を使用して計算する場合、計算時間が多大になるという課題がある。
この発明の目的は、有限要素法や境界要素法を使用することなく、非共形接触問題で行われている数値計算方法を拡張して共形接触問題にも適用できるようにすることにより、比較的簡単に短時間で、かつ精度良く、共形接触における接触面圧の計算、さらに表面下の応力の計算をできるようにした方法および装置を提供することを目的とする。
この発明の共形接触下の接触面圧および表面下応力計算方法は、半無限体に集中荷重が作用したときの変位や表面下応力の解析解を共形接触問題に適用して、面圧および表面下応力を計算する手法であって、解析解として、法線方向荷重作用時はブジネスク（Boussinesq）の解、接線方法荷重作用時はセルッティ（Cerruti)の解を用い、法線方向および接線方法の両方向の荷重作用時は、上記両方の解を重ね合わせる手法である。
具体的には、この発明の共形接触下の接触面圧計算方法は、外力が作用したときに互いに共形接触（接触領域がほぼ平面上にあると見なせない接触）を行う２つの接触物体Ｍ１，Ｍ２の接触面における面圧を、コンピュータを用いて計算する方法であって、
前記接触物体の接触面付近の形状、弾性係数、ポアソン比、および作用する外力を入力して記憶手段に記憶させる入力過程（Ｓ１）と、入力された値を用いて前記面圧を計算する解析過程（Ｓ２）と、計算した結果を出力する出力過程（Ｓ３）とでなり、
前記解析過程（Ｓ２）として、次の各過程（Ｔ１，Ｔ２，Ｔ４，Ｔ５）を含む。
前記接触物体の接触面付近の形状、弾性係数、ポアソン比、および作用する外力を入力して記憶手段に記憶させる入力過程（Ｓ１）と、入力された値を用いて前記面圧を計算する解析過程（Ｓ２）と、計算した結果を出力する出力過程（Ｓ３）とでなり、
前記解析過程（Ｓ２）として、次の各過程（Ｔ１，Ｔ２，Ｔ４，Ｔ５）を含む。
すなわち、前記２つの接触物体Ｍ１，Ｍ２の互いに接する面の近傍に接触面形状に近い形状の基準曲面Ｓを設定してこの基準曲面ＳをセルＣが縦横に並ぶメッシュに分割し、接触物体Ｍ１，Ｍ２間の距離を計算する点である計算点Ａをセル上に定めるメッシュ分割過程（Ｔ１）と、
各セルＣに均一な単位分布荷重が作用したときの前記接触物体Ｍ１，Ｍ２の表面の前記各計算点Ａに対応する点の変位への影響量を算出する変位影響算出過程（Ｔ２）と、
各セルＣの面圧Ｐの合計が外力Ｆと釣り合い、かつ各セルＣにつき面圧Ｐと前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧Ｐを算出する接触面圧計算過程（Ｔ３）（実施形態では「接触面圧・接触領域計算過程）とを含む。
なお、上記基準曲面Ｓは、２つの接触物体Ｍ１，Ｍ２の互いに接する面の接触面形状に近い形状であることが好ましい。
各セルＣに均一な単位分布荷重が作用したときの前記接触物体Ｍ１，Ｍ２の表面の前記各計算点Ａに対応する点の変位への影響量を算出する変位影響算出過程（Ｔ２）と、
各セルＣの面圧Ｐの合計が外力Ｆと釣り合い、かつ各セルＣにつき面圧Ｐと前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧Ｐを算出する接触面圧計算過程（Ｔ３）（実施形態では「接触面圧・接触領域計算過程）とを含む。
なお、上記基準曲面Ｓは、２つの接触物体Ｍ１，Ｍ２の互いに接する面の接触面形状に近い形状であることが好ましい。
前記単位分布荷重による変位影響量算出過程（Ｔ２）では、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとして、
前記外力Ｆが法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解であるブジネスク（Boussinesq）の解を積分することで前記影響量を計算する。
前記外力Ｆが法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解であるブジネスク（Boussinesq）の解を積分することで前記影響量を計算する。
前記外力Ｆが接線方向の力である場合は、次式（ａ）〜（ｃ）で表される単位集中荷重作用時の変位解析解であるセルッティ（Cerruti)の解を積分することで前記影響量を計算する。
セルに対して法線方向の力および接線方法の力を含む場合は、前記変位影響量算出過程として前記法線方向の力である場合の影響量の計算および接線方向の力である場合の影響量の計算の両方を行い、かつ前記接触面圧計算過程（Ｔ３）として、各々の方向の力の場合の各セルＣの面圧Ｐを算出し、両方の場合の計算結果を重ね合わせる。
この発明方法によると、半無限体に集中荷重が作用したときの変位の解析解を共形接触問題に適用して、解析解として、法線方向荷重作用時はブジネスクの解、接線方法荷重作用時はセルッティの解を用い、法線方向および接線方法の両方向の荷重作用時は、上記両方の解を重ね合わせるようにしたため、比較的簡単に、したがって短時間で、共形接触における接触面圧の計算が行える。また、有限要素法を用いて接触面圧を計算した場合と同程度に精度の良い計算結果が得られることが、この発明方法と有限要素法との計算結果を比較することで、確認できた。
この発明は、このように、従来では共形接触には半無限体が適用できないと考えられていた共形接触問題に、ブジネスクやセルッティの半無限体近似の解を適用し、良好な面圧の解析結果を得ることができる。
なお、この発明方法は、非共形接触の場合にも適用可能であるが、非共形接触の場合はより簡単な計算で面圧を求めることができるため、共形接触の場合に効果的となる。
この発明は、このように、従来では共形接触には半無限体が適用できないと考えられていた共形接触問題に、ブジネスクやセルッティの半無限体近似の解を適用し、良好な面圧の解析結果を得ることができる。
なお、この発明方法は、非共形接触の場合にも適用可能であるが、非共形接触の場合はより簡単な計算で面圧を求めることができるため、共形接触の場合に効果的となる。
この発明方法において、表面下応力を計算する場合は、上記の方法で各セルＣの面圧Ｐを計算した後に、表面下応力を計算する。
この表面下応力計算過程において、
セルに対する荷重が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルＣの全域にわたり積分して、かつ全セルＣの影響の和をとって表面下応力を求める。
この表面下応力計算過程において、
セルに対する荷重が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルＣの全域にわたり積分して、かつ全セルＣの影響の和をとって表面下応力を求める。
セルに対する荷重が接線方向の力である場合は、次式（ｄ）〜（ｉ）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルＣの全域にわたり積分して、かつ全セルＣの影響の和をとって表面下応力を求める。
セルに対する荷重が法線方向の力および接線方向の力を含む場合は、上記法線方向の力である場合の表面下応力の計算、および上記接線方向の力である場合の表面下応力の計算の両方を行い、両方の計算結果を重ね合わせる。
このように表面下応力を計算するについても、計算が比較的簡単で短時間で行え、また精度の良い表面下応力の計算が行える。
この発明方法において、前記解析過程の途中で、前記メッシュに分割する大きさを修正してその修正結果により、メッシュ分割後の各過程を再度行うメッシュ修正過程（Ｔ６）を含むようにしても良い。
接触物体の接触面の形状によっては、接触領域の縁がギザギザになるなどして、離散化の誤差が大きくなる場合があるが、解析途中でメッシュを修正することで、離散化の誤差を低減させることができる。
接触物体の接触面の形状によっては、接触領域の縁がギザギザになるなどして、離散化の誤差が大きくなる場合があるが、解析途中でメッシュを修正することで、離散化の誤差を低減させることができる。
この発明方法は、前記２つの接触物体が、例えば、球面座と、この球面座上を滑る球面部を有する部材とであっても良い。
このような球面座による支持を行う場合、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要であるため、この発明方法を効果的に適用することができる。
また、この発明方法は、前記２つの接触物体が、例えば、転がり軸受における内輪または外輪である軌道輪と、ボール等の転動体とであっても良い。
転がり軸受の軌道輪とボール等の転動体との接触においても、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要であるため、この発明方法を効果的に適用することができる。
このような球面座による支持を行う場合、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要であるため、この発明方法を効果的に適用することができる。
また、この発明方法は、前記２つの接触物体が、例えば、転がり軸受における内輪または外輪である軌道輪と、ボール等の転動体とであっても良い。
転がり軸受の軌道輪とボール等の転動体との接触においても、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要であるため、この発明方法を効果的に適用することができる。
この発明の共形接触下の接触面圧計算装置は、外力が作用したときに互いに共形接触を行う２つの接触物体の接触面における面圧を、計算する装置であって、
前記接触物体Ｍ１，Ｍ２の接触面付近の形状、弾性係数、ポアソン比、および作用する外力Ｆを入力して記憶手段６に記憶させる入力処理手段４と、入力された値を用いて前記面圧を計算する解析手段７と、計算した結果を出力する出力処理手段５とでなり、
前記解析手段７として、次の各手段８，９，１１，１２を備える。
前記接触物体Ｍ１，Ｍ２の接触面付近の形状、弾性係数、ポアソン比、および作用する外力Ｆを入力して記憶手段６に記憶させる入力処理手段４と、入力された値を用いて前記面圧を計算する解析手段７と、計算した結果を出力する出力処理手段５とでなり、
前記解析手段７として、次の各手段８，９，１１，１２を備える。
すなわち、前記２つの接触物体Ｍ１，Ｍ２の互いに接する面の近傍に、基準曲面Ｓを設定してこの基準曲面ＳをセルＣが縦横に並ぶメッシュに分割し、接触物体Ｍ１，Ｍ２間の距離を計算する点である計算点Ａをセル上に定めるメッシュ分割手段８と、
各セルＣに均一な単位分布荷重が作用したときの前記接触物体Ｍ１，Ｍ２の表面の前記各計算点Ａに対応する点の変位への影響量を算出する変位影響算出手段９と、
各セルＣの面圧Ｐの合計が外力と釣り合い、かつ各セルＣにつき面圧Ｐと前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧Ｐを算出する接触面圧計算手段１０（実施形態では「接触面圧・接触領域計算手段」）とを含む。
この他に、各計算点Ａにおける接触の有無を判定して接触領域を定める接触領域計算手段１０を設けることが好ましい。
各セルＣに均一な単位分布荷重が作用したときの前記接触物体Ｍ１，Ｍ２の表面の前記各計算点Ａに対応する点の変位への影響量を算出する変位影響算出手段９と、
各セルＣの面圧Ｐの合計が外力と釣り合い、かつ各セルＣにつき面圧Ｐと前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧Ｐを算出する接触面圧計算手段１０（実施形態では「接触面圧・接触領域計算手段」）とを含む。
この他に、各計算点Ａにおける接触の有無を判定して接触領域を定める接触領域計算手段１０を設けることが好ましい。
前記変位影響量算出手段９は、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとして、
前記外力Ｆが法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解を積分することで前記影響量を計算する。
前記外力Ｆが法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解を積分することで前記影響量を計算する。
前記外力Ｆが接線方向の力である場合は、次式（ａ）〜（ｃ）で表される単位集中荷重作用時の変位解析解を積分することで前記影響量を計算する。
セルに対する荷重が法線方向の力および接線方法の力を含む場合は、前記法線方向の力である場合の影響量の計算および接線方向の力である場合の計算の両方を行うものとし、
かつ前記接触面圧計算手段１０は、セルに対する荷重が法線方向の力および接線方法の力を含む場合は、各々の方向の力の場合の各セルＣの面圧Ｐを算出し、両方の場合の計算結果を重ね合わせるものとする。
かつ前記接触面圧計算手段１０は、セルに対する荷重が法線方向の力および接線方法の力を含む場合は、各々の方向の力の場合の各セルＣの面圧Ｐを算出し、両方の場合の計算結果を重ね合わせるものとする。
この構成の面圧計算装置によると、この発明方法を実施して、共形接触を行う接触面における面圧を、比較的簡単に短時間で計算でき、かつ精度良く計算することができる。
この発明の面圧計算装置において、さらに表面下応力の計算機能を付加する場合は、次の表面下応力計算手段１３を設ける。
この表面下応力計算手段１３は、セルに対する荷重が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルの全域にわたり積分して、かつ全セルの影響の和をとって表面下応力を求める。
この表面下応力計算手段１３は、セルに対する荷重が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルの全域にわたり積分して、かつ全セルの影響の和をとって表面下応力を求める。
セルに対する荷重が接線方向の力である場合は、次式（ｄ）〜（ｉ）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルＣの全域にわたり積分して、かつ全セルＣの影響の和をとって表面下応力を求める。
セルに対する荷重が法線方向の力および接線方向の力を含む場合は、上記法線方向の力である場合の表面下応力の計算、および上記接線方向の力である場合の表面下応力の計算の両方を行い、両方の計算結果を重ね合わせる。
この構成の場合、この発明の表面下応力計算方法を実施して、共形接触を行う接触面における表面下応力を、比較的簡単に短時間で計算でき、かつ精度良く計算することができる。
この発明の共形接触下の接触面圧計算方法および計算装置によると、半無限体に集中荷重が作用したときの変位の解析解を共形接触問題に適用して、解析解として、法線方向荷重作用時はブジネスクの解、接線方法荷重作用時はセルッティの解を用い、法線方向および接線方法の両方向の荷重作用時は、上記両方の解を重ね合わせるようにしたため、比較的簡単に、したがって短時間で、共形接触における接触面圧の計算が行え、しかも有限要素法を用いて計算した場合と同程度に精度の良い計算結果を得ることができる。
この発明の共形接触下の接触面圧，表面下応力の計算方法および計算装置によると、接触面圧に加えて、さらに表面下応力についても、比較的簡単に短時間で、また精度良く、共形接触における表面下応力の計算を行うことができる。
この発明の共形接触下の接触面圧，表面下応力の計算方法および計算装置によると、接触面圧に加えて、さらに表面下応力についても、比較的簡単に短時間で、また精度良く、共形接触における表面下応力の計算を行うことができる。
この発明の一実施形態を図１ないし図７と共に説明する。図１は、この発明方法の各過程の流れ図を、図２はこの発明方法を実施する装置の概念構成のブロック図を、図４はこの発明方法で用いる基準曲面の例をそれぞれ示す。
この共形接触下の接触面圧および表面下応力計算方法は、外力が作用したときに互いに共形接触を行う２つの接触物体Ｍ１，Ｍ２（図４）の接触面における面圧を、コンピュータを用いて計算する方法であって、前記接触物体Ｍ１，Ｍ２の接触面付近の形状、弾性係数（ここでは、横弾性係数Ｇと、ポアソン比ν）、および作用する外力Ｆを入力して記憶手段に記憶させる入力過程（Ｓ１）と、入力された値を用いて前記面圧を計算する解析過程（Ｓ２）と、計算した結果を出力する出力過程（Ｓ３）とでなり、解析過程（Ｓ２）として、順次行われる次のメッシュ分割過程（Ｔ１）、変位影響算出過程（Ｔ２）、接触面圧・接触領域計算過程（Ｔ３）、および表面下応力計算過程（Ｔ４）を含む。
計算装置となるコンピュータ１（図２）は、中央処理装置（ＣＰＵ）およびメモリ等の記憶手段（いずれも図示せず）を有し、このコンピュータ１のハードウェア（オペレーションシステムを含む）、およびこのコンピュータ１に実行させる接触面圧・表面下応力計算プログラム（図示せず）により、図２に概念構成で示す各手段が構成される。
入力処理手段４は、入力装置２から、接触物体Ｍ１，Ｍ２の接触面付近の形状、弾性係数、ポアソン比、および作用する外力Ｆを入力して入力情報記憶手段６に記憶させる手段である。ここで言う「形状」は、大きさを含む概念である。入力処理手段４は、上記の形状等の情報の他に、後述の解析手段７でメッシュ分割する枡目の各軸方向の大きさやセル個数（ＮＡ，ＮＢ）等を、入力装置２から入力して入力情報記憶手段６に記憶させる機能を備えるものとしても良い。また、入力処理手段４は、入力装置２から入力された情報に対して、解析手段７で演算するための前処理を行う機能を有するものとしても良い。入力装置２は、キーボード等であっても、また通信手段や記憶素子等であっても良い。入力情報記憶手段６は、コンピュータ１の上記記憶手段における所定の記憶領域である。出力処理手段５は、解析手段７で計算した結果を、液晶ディスプレイ，プリンタ，あるいは通信機器等からなる出力装置３に出力する手段である。
入力処理手段４は、入力装置２から、接触物体Ｍ１，Ｍ２の接触面付近の形状、弾性係数、ポアソン比、および作用する外力Ｆを入力して入力情報記憶手段６に記憶させる手段である。ここで言う「形状」は、大きさを含む概念である。入力処理手段４は、上記の形状等の情報の他に、後述の解析手段７でメッシュ分割する枡目の各軸方向の大きさやセル個数（ＮＡ，ＮＢ）等を、入力装置２から入力して入力情報記憶手段６に記憶させる機能を備えるものとしても良い。また、入力処理手段４は、入力装置２から入力された情報に対して、解析手段７で演算するための前処理を行う機能を有するものとしても良い。入力装置２は、キーボード等であっても、また通信手段や記憶素子等であっても良い。入力情報記憶手段６は、コンピュータ１の上記記憶手段における所定の記憶領域である。出力処理手段５は、解析手段７で計算した結果を、液晶ディスプレイ，プリンタ，あるいは通信機器等からなる出力装置３に出力する手段である。
解析手段７は、入力処理手段４で入力されて入力情報記憶手段６に記憶されたデータを用いて前記面圧および表面下応力を計算する手段であり、メッシュ分割手段８、変位影響算出手段９、接触面圧・接触領域計算手段１０、および表面下応力計算手段１１を備える。解析手段７に備えられた上記各手段８〜１２は、それぞれ図１（Ｂ）の流れ図におけるメッシュ分割過程（Ｔ１）、変位影響算出過程（Ｔ２）、接触面圧・接触領域計算過程（Ｔ３）、および表面下応力計算過程（Ｔ４）を実行する手段であり、その実行に必要な計算式や設定値等を記憶している。上記各手段８〜１１の具体的な内容は、次に述べる上記各過程（Ｔ１）〜（Ｔ４）の説明の通りである。
解析手段７は、上記の他に、後に説明するメッシュ修正手段１２を有している。
解析手段７は、上記の他に、後に説明するメッシュ修正手段１２を有している。
図１（Ａ）の解析過程（Ｓ２）を、同図（Ｂ）および図４〜図７と共に説明する。
この解析仮定（Ｓ２）は、法線方向加重が作用したときのハーネット（Hartnett) の手法を拡張して共形接触問題（Conformal Problem ）を解く手法であり、接線方向荷重作用時はセルッティ（Cerruti)の解を用いて法線方向荷重の作用時と同様に解き、また、法線方向と接線方向の両方の荷重が作用したときは、両方法による計算結果を重ね合わせる。
この解析仮定（Ｓ２）は、法線方向加重が作用したときのハーネット（Hartnett) の手法を拡張して共形接触問題（Conformal Problem ）を解く手法であり、接線方向荷重作用時はセルッティ（Cerruti)の解を用いて法線方向荷重の作用時と同様に解き、また、法線方向と接線方向の両方の荷重が作用したときは、両方法による計算結果を重ね合わせる。
まず、メッシュ分割過程（Ｔ１）として、図４のように、接触物体Ｍ１，Ｍ２の間に基準曲面Ｓを設定し、この基準曲面ＳをセルＣが縦横に並ぶメッシュに分割し、接触物体Ｍ１，Ｍ２間の距離を計算する点である計算点ＡをセルＣ上に定める。上記基準曲面Ｓは、必ずしも２つの接触物体Ｍ１，Ｍ２の間でなくても良いが、前記２つの接触物体Ｍ１，Ｍ２の互いに接する面の近傍に設定する。同図では、ｘ軸とｙ軸の交点Ｏを球の極にとり、球の赤道で２つの接触物体Ｍ１，Ｍ２と接触する基準面Ｓを用いている。また、上記基準曲面Ｓは、２つの接触物体Ｍ１，Ｍ２の接触面形状に近い形状の曲面とする。
なお、基準曲面Ｓについて、接触面形状の正確な形状は本来わからないが、凸球面と凹球面の接触では接触面形状はほぼ球面とみなせる。ボールと内輪の接触では接触面形状は円筒面の一部に近いとみなせる（解析上もこの基準面を採用している）。従来の「接触面が平坦と見なせる」場合は平面である。複雑な曲面同士の接触では上記のような簡単な曲面で表現できない場合が出てくる。基準面形状をどのようにするかは解析精度に影響してくるため、できるだけ接触面形状に近い方が良いが、複雑にすると計算が難しくなるというデメリットがある。従って「基準曲面Ｓ」は解析する人がケースバイケースで適宜決定していく。
計算点Ａは、接触物体Ｍ１，Ｍ２間の距離を離散的な点で評価するものとするためにセルＣの位置を代表させる点である。計算点Ａは、セルＣ内の適宜の位置に設定すれば良いが、ここではセルＣの中心にとる。
メッシュ分割する枡目の大きさ、つまりセルＣの大きさは、適宜設定すればよいが、入力処理手段４からの入力等によって自由に設定可能としても良い。また、基準曲面Ｓの広さは、接触物体Ｍ１，Ｍ２の接触領域を全て含むように、接触領域より広めに適宜設定する。ここでは、セルＣがＸ軸方向にＮＡ個、Ｙ軸方向にＮＢ個並ぶ基準曲面Ｓとする。
なお、基準曲面Ｓについて、接触面形状の正確な形状は本来わからないが、凸球面と凹球面の接触では接触面形状はほぼ球面とみなせる。ボールと内輪の接触では接触面形状は円筒面の一部に近いとみなせる（解析上もこの基準面を採用している）。従来の「接触面が平坦と見なせる」場合は平面である。複雑な曲面同士の接触では上記のような簡単な曲面で表現できない場合が出てくる。基準面形状をどのようにするかは解析精度に影響してくるため、できるだけ接触面形状に近い方が良いが、複雑にすると計算が難しくなるというデメリットがある。従って「基準曲面Ｓ」は解析する人がケースバイケースで適宜決定していく。
計算点Ａは、接触物体Ｍ１，Ｍ２間の距離を離散的な点で評価するものとするためにセルＣの位置を代表させる点である。計算点Ａは、セルＣ内の適宜の位置に設定すれば良いが、ここではセルＣの中心にとる。
メッシュ分割する枡目の大きさ、つまりセルＣの大きさは、適宜設定すればよいが、入力処理手段４からの入力等によって自由に設定可能としても良い。また、基準曲面Ｓの広さは、接触物体Ｍ１，Ｍ２の接触領域を全て含むように、接触領域より広めに適宜設定する。ここでは、セルＣがＸ軸方向にＮＡ個、Ｙ軸方向にＮＢ個並ぶ基準曲面Ｓとする。
変位影響算出過程（Ｔ２）では、メッシュ分割過程（Ｔ１）の後、各セルＣに均一な分布荷重Ｐが作用したときの接触物体Ｍ１，Ｍ２の表面の前記各計算点Ａに対応する点の変位への影響量を算出する。この場合に、ここでは、個々のセルＣには全体に均一な面圧Ｐが作用すると仮定するが、１次関数または２次関数等の所定の関数に従いセルＣ内の位置によって変わる面圧が作用すると仮定しても良い。
上記影響量の算出は、次のように行う。法線方向荷重が作用したときは、上記面圧分布によって生じる２つの接触物体Ｍ１，Ｍ２の表面上の点（上記計算点Ａに対応する点）の変位量は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解（ブジネスク（Boussinesq）の解）を、次式（４）〜（６）のように積分することで得る。
なお、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとする。
上記影響量の算出は、次のように行う。法線方向荷重が作用したときは、上記面圧分布によって生じる２つの接触物体Ｍ１，Ｍ２の表面上の点（上記計算点Ａに対応する点）の変位量は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解（ブジネスク（Boussinesq）の解）を、次式（４）〜（６）のように積分することで得る。
なお、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとする。
図６は、ブジネスク（Boussinesq）の解の座標系を示す。同図のような半無限体平面の原点Ｏに単位荷重が作用するとき、任意の点であるＱ点のｘ，ｙ，ｚ方向の変位ｕ，ｖ，ｗをセル全域にわたって積分し、均一な面圧が作用する場合は、その積分値Ｕ，Ｖ，Ｗは次式（４）〜（６）となる。
（７）式はセルＣij全域の荷重分布ＰijがＱkl点の変位に影響を及ぼす影響ベクトルを表している。
なお、Ｃijは、ｘ軸方向にｉ番目でｙ軸方向にｊ番目にあるセルＣを示す。Ｑklは、ｘ軸方向にｋ番目でｙ軸方向にｌ番目の点Ａを示す（なお、ｉ：１〜ＮＡ、ｊ：１〜ＮＢ、ｋ：１〜ＮＡ、ｌ：１〜ＮＢ）。
また、上記（４）〜（６）式において、θおよびφは、図４のように極座標で示されるセルＣの中心である計算点Ａのｘ軸回りの角度およびｙ軸回りの角度を示す。ａ，ｂは、それぞれセルＣの中心である計算点ＡからそのセルＣの周縁までのｘ軸回りの角度およびＹ軸回りの角度を示す。
なお、Ｃijは、ｘ軸方向にｉ番目でｙ軸方向にｊ番目にあるセルＣを示す。Ｑklは、ｘ軸方向にｋ番目でｙ軸方向にｌ番目の点Ａを示す（なお、ｉ：１〜ＮＡ、ｊ：１〜ＮＢ、ｋ：１〜ＮＡ、ｌ：１〜ＮＢ）。
また、上記（４）〜（６）式において、θおよびφは、図４のように極座標で示されるセルＣの中心である計算点Ａのｘ軸回りの角度およびｙ軸回りの角度を示す。ａ，ｂは、それぞれセルＣの中心である計算点ＡからそのセルＣの周縁までのｘ軸回りの角度およびＹ軸回りの角度を示す。
上記の（７）式で示されるＱkl点の変位のうち、面の法線方向成分のみをとるために、Ｑkl点の法線方向のベクトルとの内積をとると、次式（８）のようになる。
この（８）式の値ｆijklが、セルＣijに分布荷重Ｐijが作用したときに、他のセルＣklの中心位置Ｑklに対応する点の変位への影響量である。
この（８）式の値ｆijklが、セルＣijに分布荷重Ｐijが作用したときに、他のセルＣklの中心位置Ｑklに対応する点の変位への影響量である。
次に、接触面圧・接触領域計算過程（Ｔ４）として、各セルＣijの面圧Ｐijの合計が外力Ｆと釣り合い、かつ各セルＣijにつき面圧Ｐijと前記影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧を算出する。
すなわち、以下の（９）式および（１０）式が釣り合うように各セルの面圧値を決めれば、離散化された全体の面圧分布が求まる。
なお、上記のように（９）式を解く過程で、２つの接触物体Ｍ１，Ｍ２間の接近量と初期隙間の差Ｄijが正の値となるセルＣijについては、そのセルＣijの面圧は零としており、この接近量と初期隙間の差Ｄijが零または負の値となるセルＣijの並び範囲が、上記「接触面圧・接触領域計算過程」で言う接触領域である。
すなわち、以下の（９）式および（１０）式が釣り合うように各セルの面圧値を決めれば、離散化された全体の面圧分布が求まる。
なお、上記のように（９）式を解く過程で、２つの接触物体Ｍ１，Ｍ２間の接近量と初期隙間の差Ｄijが正の値となるセルＣijについては、そのセルＣijの面圧は零としており、この接近量と初期隙間の差Ｄijが零または負の値となるセルＣijの並び範囲が、上記「接触面圧・接触領域計算過程」で言う接触領域である。
（９）式、（１０）式を、コンピュータにおける接触面圧・接触領域計算手段１０によって解くアルゴリズムの一例を図３に示す。
まず、α値（：２つの接触物体Ｍ１，Ｍ２の接近量）の初期値を設定する（ステップＨ１）。
そのα値に応じて各セルＣijの接触面圧を設定する（ステップＨ２）。この場合に、ｉ＝１〜ＮＡの範囲、およびｊ＝１〜ＮＢの範囲で全てのセルＣijの接触面圧Ｐijを設定する。 その設定した各セルＣijの接触面圧Ｐijで（９）式が成り立つか否かを判定し、成り立たない場合は、セルＣijの接触面圧の設定過程（ステップＨ２）に戻る（Ｈ３）。ステップＨ２では、各セルＣijの接触面圧Ｐijを再度設定し直し、その設定したセルＣijの接触面圧Ｐijで（９）式が成り立つか否かを判定し、成り立たない場合はセルＣijの接触面圧Ｐijを設定過程（Ｈ２）に戻る。このような処理を（９）式が成り立つまで繰替し、（９）式が成り立つと、ステップＨ４に進む。
ステップＨ４では、（１０）式が成り立つか否かを判定し、成り立たない場合は、ステップＨ１に戻ってα値を、再度設定し直す。
この後、前記と同様にしてステップＨ２〜Ｈ４，Ｈ１の処理を繰り返し、（１０）式が成り立つと、処理を終了する。
このようにして求めた、式（９）、（１０）が成り立つときの、ステップＨ２で設定した各セルＣijの接触面圧Ｐijが、接触面圧・接触領域計算手段１０の出力となる接触面圧の値である。
なお、セルＣijの接触面圧の最初の値は何でもよいが、できるだけ実際の面圧に近い値を設定した方が早く収束する。
まず、α値（：２つの接触物体Ｍ１，Ｍ２の接近量）の初期値を設定する（ステップＨ１）。
そのα値に応じて各セルＣijの接触面圧を設定する（ステップＨ２）。この場合に、ｉ＝１〜ＮＡの範囲、およびｊ＝１〜ＮＢの範囲で全てのセルＣijの接触面圧Ｐijを設定する。 その設定した各セルＣijの接触面圧Ｐijで（９）式が成り立つか否かを判定し、成り立たない場合は、セルＣijの接触面圧の設定過程（ステップＨ２）に戻る（Ｈ３）。ステップＨ２では、各セルＣijの接触面圧Ｐijを再度設定し直し、その設定したセルＣijの接触面圧Ｐijで（９）式が成り立つか否かを判定し、成り立たない場合はセルＣijの接触面圧Ｐijを設定過程（Ｈ２）に戻る。このような処理を（９）式が成り立つまで繰替し、（９）式が成り立つと、ステップＨ４に進む。
ステップＨ４では、（１０）式が成り立つか否かを判定し、成り立たない場合は、ステップＨ１に戻ってα値を、再度設定し直す。
この後、前記と同様にしてステップＨ２〜Ｈ４，Ｈ１の処理を繰り返し、（１０）式が成り立つと、処理を終了する。
このようにして求めた、式（９）、（１０）が成り立つときの、ステップＨ２で設定した各セルＣijの接触面圧Ｐijが、接触面圧・接触領域計算手段１０の出力となる接触面圧の値である。
なお、セルＣijの接触面圧の最初の値は何でもよいが、できるだけ実際の面圧に近い値を設定した方が早く収束する。
表面下応力計算過程（Ｔ４）では、上記のように求めた面圧分布に対して、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルＣの全域にわたり積分して、かつ全セルＣの影響の和をとって表面下応力が求まる。
接線方向荷重作用時の表面下応力は次式（ｄ）〜（ｉ）で求める。
法線方向と接線方向の両方の荷重が作用したときは、両者を重ね合わせる。例えば、セルに対して斜め方向の荷重が作用したときは、法線方向の荷重の成分と接線方向の荷重の成分を求め、法線方向の荷重および接線方向の荷重に対して、それぞれ面圧分布および表面下応力を計算し、これら各方向の面圧分布の計算結果および表面下応力の計算結果をそれぞれ重ね合わせる。この重ね合わせは単に加算すれば良い。
なお、セルに対して斜め方向の荷重が作用したときのその荷重の入力については、例えば入力過程（Ｓ１）で荷重の大きさと方向を入力し、この入力過程（Ｓ１）において、法線方向の荷重の成分と接線方向の荷重の成分を求める前処理を行うようにしても良く、また、予め入力過程（Ｓ１）で法線方向の荷重の成分と接線方向の荷重の成分をそれぞれ入力するようにしても良い。
なお、セルに対して斜め方向の荷重が作用したときのその荷重の入力については、例えば入力過程（Ｓ１）で荷重の大きさと方向を入力し、この入力過程（Ｓ１）において、法線方向の荷重の成分と接線方向の荷重の成分を求める前処理を行うようにしても良く、また、予め入力過程（Ｓ１）で法線方向の荷重の成分と接線方向の荷重の成分をそれぞれ入力するようにしても良い。
この実施形態の方法の計算方法の有効性を確認するため、同方法を実施した結果を他の計算方法と比較して示す。
図８のような凸球面の接触物体Ｍ１と凹球面の接触物体Ｍ２との接触において、軽荷重時（非共形問題（Counter-formal Problem））と重荷重時（共形問題（Conformal Problem ）に分けて計算結果を示す。
軽荷重時にこの実施形態の計算方法により求めた面圧分布を図９に示す。また、ヘルツの接触理論との比較を表１に示す。両者はほぼ一致していることがわかる。
図８のような凸球面の接触物体Ｍ１と凹球面の接触物体Ｍ２との接触において、軽荷重時（非共形問題（Counter-formal Problem））と重荷重時（共形問題（Conformal Problem ）に分けて計算結果を示す。
軽荷重時にこの実施形態の計算方法により求めた面圧分布を図９に示す。また、ヘルツの接触理論との比較を表１に示す。両者はほぼ一致していることがわかる。
重荷重時にこの実施形態の計算方法によって求めた面圧分布を図１０(a) に示す。この結果の妥当性を有限要素法（ＦＥＭ）との比較において示すため、図１０(b) にＦＥＭの結果を示す。両者は良好な一致が得られており、この実施形態の計算方法の有効性が証明できた。
なお、ＦＥＭの計算時間はPentium （登録商標）４（クロック周波数3.7GHz）で約６分、この実施形態の計算ではPentium 4 （クロック周波数2.2GHz）で約20秒であり、計算時間が大幅に短縮できる。
なお、ＦＥＭの計算時間はPentium （登録商標）４（クロック周波数3.7GHz）で約６分、この実施形態の計算ではPentium 4 （クロック周波数2.2GHz）で約20秒であり、計算時間が大幅に短縮できる。
図１１は現在市販されている非共形接触問題（Counter-formal Problem）用のソフトウェアにより重荷重時の解析を実施した結果であるが、全く使用できないことがわかる。
図１２ないし図１５は、この実施形態の計算方法の有効性を確認するために行った他の例を示す。この例は、図１５に示す内輪２１、外輪２２、ボール２３、および保持器２４により構成される転がり軸受であるアンギュラ玉軸受につき、第１の接触物体Ｍ１をボール２３、第２の接触物体Ｍ２を内輪２１として、各方法で面圧分布を求めた例である。軸受諸元は、ボール径φ：７．９３７５mm、ＰＣＤ：φ３４．５mm、溝曲率（ここでは、溝の曲率半径／ボールの曲率半径を言う）：１．０２、接触角１５°、である。
図１２〜図１４は、上記軸受の大荷重負荷時の内輪２１の軌道面２１ａとボール２３との接触面圧を示す。このうち、図１２はＦＥＭの図であり、ＭＮ軸回りにおけるボールと軌道面との接触箇所を示す。同図において、明るく表示されているのは非常にメッシュが細かい（12〜20μｍ）ことを示している。
図１３（ａ）は実施形態の方法で求めた面圧を、図１３（ｂ）は従来のハーネットの手法での面圧をそれぞれ示す。
図１４は、接触楕円長軸方向の面圧（短軸幅の中央の切断面での面圧）である。同図（ａ）は実施形態の方法で求めた面圧、同図（ｂ）は従来のハーネットの手法での面圧、同図（ｃ）はＦＥＭでの面圧をそれぞれ示す。
図１３（ａ）は実施形態の方法で求めた面圧を、図１３（ｂ）は従来のハーネットの手法での面圧をそれぞれ示す。
図１４は、接触楕円長軸方向の面圧（短軸幅の中央の切断面での面圧）である。同図（ａ）は実施形態の方法で求めた面圧、同図（ｂ）は従来のハーネットの手法での面圧、同図（ｃ）はＦＥＭでの面圧をそれぞれ示す。
これら図１２〜図１４を比較して、従来のハーネットの手法で求めた面圧は、ＦＥＭで求めた面圧と大きく異なっているが、この実施形態の方法で求めた面圧は、ＦＥＭで求めた面圧と良好な一致が得られていることが分かる。
この共形接触下の接触面圧および表面下応力計算方法によると、このように、半無限体に集中荷重が作用したときの変位の解析解を共形接触問題に適用して、解析解として、法線方向荷重作用時はブジネスクの解、接線方法荷重作用時はセルッティの解を用い、法線方向および接線方法の両方向の荷重作用時は、上記両方の解を重ね合わせるようにしたため、比較的簡単に、したがって短時間で、共形接触における接触面圧および表面下応力の計算が行える。しかも、有限要素法を用いて接触面圧を計算した場合と同程度に精度の良い計算結果を得ることができる。
また、この実施形態の計算方法を用いれば、有限要素法や境界要素法を用いる場合の高価なソフトウェアを導入する場合の初期費用、保守量、ライセンス本数といったコスト面の削減が図れ、またこの接触面圧・表面下応力計算装置を用いる場合、機能の専用化による使用方法の簡単化が図れ、これにより利用者の教育時間，教育コストの低減が図れるという利点も得られる。
また、この実施形態の計算方法を用いれば、有限要素法や境界要素法を用いる場合の高価なソフトウェアを導入する場合の初期費用、保守量、ライセンス本数といったコスト面の削減が図れ、またこの接触面圧・表面下応力計算装置を用いる場合、機能の専用化による使用方法の簡単化が図れ、これにより利用者の教育時間，教育コストの低減が図れるという利点も得られる。
なお、上記実施形態において、接触物体Ｍ１，Ｍ２の形状によっては、接触領域の縁がギザギザにあることがあるが、そのような場合は、解析途中で、例えば接触面圧・接触領域計算過程（Ｔ３）の後、メッシュ分割の修正、つまりセルＣの各軸ｘ，ｙ方向の大きさを修正してその修正結果により、メッシュ分割後の各過程を再度行わせるメッシュ修正過程（Ｔ５）を含むようにしても良い。
図２のメッシュ修正手段１３は、このようなメッシュ修正過程（Ｔ５）を実行する手段である。
図２のメッシュ修正手段１３は、このようなメッシュ修正過程（Ｔ５）を実行する手段である。
また、上記実施形態において、接触物体Ｍ１，Ｍ２の具体例の説明は省略したが、接触物体Ｍ１，Ｍ２は、例えば、球面座と、この球面座上を滑る球面部を有する部材とであっても良い。
このような球面座による支持を行う場合、作用荷重が大きくて、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要となるため、この発明方法を効果的に適用することができる。
このような球面座による支持を行う場合、作用荷重が大きくて、接触領域がほぼ平面上にあると見なせず、共形接触の問題として面圧等を求めることが必要となるため、この発明方法を効果的に適用することができる。
１…コンピュータ
４…入力処理手段
５…出力処理手段
６…入力情報記憶手段
７…解析手段
８…メッシュ分割手段
９…変位影響算出手段
１０…接触面圧・接触領域計算手段
１１…表面下応力計算手段
Ｍ１，Ｍ２…接触物体
Ｓ…基準曲面
Ｃ…セル
４…入力処理手段
５…出力処理手段
６…入力情報記憶手段
７…解析手段
８…メッシュ分割手段
９…変位影響算出手段
１０…接触面圧・接触領域計算手段
１１…表面下応力計算手段
Ｍ１，Ｍ２…接触物体
Ｓ…基準曲面
Ｃ…セル
Claims (7)
- 外力が作用したときに互いに共形接触を行う２つの接触物体の接触面における面圧を、コンピュータを用いて計算する方法であって、
前記接触物体の接触面付近の形状、弾性係数、ポアソン比、および作用する外力を入力して記憶手段に記憶させる入力過程と、入力された値を用いて前記面圧を計算する解析過程と、計算した結果を出力する出力過程とでなり、
前記解析過程として、
前記２つの接触物体の互いに接する面の近傍に、基準曲面を設定してこの基準曲面をセルが縦横に並ぶメッシュに分割し、接触物体間の距離を計算する点である計算点をセル上に定めるメッシュ分割過程と、
各セルに均一な単位分布荷重が作用したときの前記接触物体の表面の前記各計算点に対応する点の変位への影響量を算出する変位影響算出過程と、
各セルの面圧の合計が外力と釣り合い、かつ各セルにつき面圧と前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧を算出する接触面圧計算過程とを含み、
前記単位分布荷重による変位影響量算出過程では、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとして、
前記外力が法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解を積分することで前記影響量を計算し、
ことを特徴とする共形接触下の接触面圧の計算方法。 - 請求項１において、前記解析過程の途中で、前記メッシュに分割する大きさを修正してその修正結果により、メッシュ分割後の各過程を再度行うメッシュ修正過程を含む共形接触下の接触面圧計算方法。
- 請求項１または請求項２において、前記２つの接触物体が、球面座と、この球面座上を滑る球面部を有する部材とである共形接触下の接触面圧計算方法。
- 請求項１または請求項２において、前記２つの接触物体が、転がり軸受の軌道輪と転動体とである共形接触下の接触面圧計算方法。
- 外力が作用したときに互いに共形接触を行う２つの接触物体の接触面における面圧および表面下応力を、コンピュータを用いて計算する方法であって、
請求項１または請求項２または請求項３に記載の接触面圧の計算方法と、この計算方法で各セルの面圧を計算した後に、表面下応力を計算する過程を含み、
この表面下応力計算過程として、
前記外力が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルの全域にわたり積分して、かつ全セルの影響の和をとって表面下応力を求め、
共形接触下の接触面圧および表面下応力計算方法。 - 外力が作用したときに互いに共形接触を行う２つの接触物体の接触面における面圧を計算する装置であって、
前記接触物体の接触面付近の形状、弾性係数、ポアソン比、および作用する外力を入力して記憶手段に記憶させる入力処理手段と、入力された値を用いて前記面圧を計算する解析手段と、計算した結果を出力する出力処理手段とでなり、
前記解析手段として、
前記２つの接触物体の互いに接する面の近傍に、基準曲面を設定してこの基準曲面をセルが縦横に並ぶメッシュに分割し、接触物体間の距離を計算する点である計算点をセル上に定めるメッシュ分割手段と、
各セルに均一な単位分布荷重が作用したときの前記接触物体の表面の前記各計算点に対応する点の変位への影響量を算出する変位影響算出手段と、
各セルの面圧の合計が外力と釣り合い、かつ各セルにつき面圧と前記単位分布荷重による影響量の積の合計が接触領域内部において接触物体間の接近量と初期隙間の差に釣り合うという釣り合い条件式から、各セルの面圧を算出する接触面圧計算手段とを含み、
前記変位影響量算出手段は、互いに直交する３軸方向の座標位置をｘ，ｙ，ｚ、これら３軸方向の変位を、それぞれｕ，ｖ，ｗとして、
前記外力が法線方向の力である場合は、次式（１）〜（３）で表される単位集中荷重作用時の変位解析解を積分することで前記影響量を計算し、
ことを特徴とする共形接触下の接触面圧の計算装置。 - 外力が作用したときに互いに共形接触を行う２つの接触物体の接触面における面圧および表面下応力を計算する装置であって、
請求項６に記載の接触面圧の計算装置と、表面下応力計算手段とを含み、
この表面下応力計算手段は、
前記外力が法線方向の力である場合は、次式（１１）〜（１６）で与えられる単位集中荷重作用時の表面下の応力を用いて、個々のセルの全域にわたり積分して、かつ全セルの影響の和をとって表面下応力を求め、
共形接触下の接触面圧および表面下応力計算装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007026577A JP5046671B2 (ja) | 2007-02-06 | 2007-02-06 | 共形接触下の接触面圧および表面下応力計算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007026577A JP5046671B2 (ja) | 2007-02-06 | 2007-02-06 | 共形接触下の接触面圧および表面下応力計算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2008191025A JP2008191025A (ja) | 2008-08-21 |
| JP5046671B2 true JP5046671B2 (ja) | 2012-10-10 |
Family
ID=39751254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007026577A Active JP5046671B2 (ja) | 2007-02-06 | 2007-02-06 | 共形接触下の接触面圧および表面下応力計算方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5046671B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6085748B2 (ja) * | 2012-05-24 | 2017-03-01 | 長谷 陽夫 | ２物体間の接触面圧計算方法および２物体間の接触面圧計算コンピュータプログラム |
| CN103438860B (zh) * | 2013-08-20 | 2016-01-27 | 国家电网公司 | 电力隧道盾构施工地面最大沉降值预测方法 |
| JP6308544B2 (ja) * | 2013-09-12 | 2018-04-11 | 長谷 陽夫 | ２物体間の接触面圧計算方法および２物体間の接触面圧計算コンピュータプログラム |
| CN106951706B (zh) * | 2017-03-20 | 2019-07-02 | 浙江大学城市学院 | 一种地面堆载引起既有盾构隧道位移计算方法 |
| CN110147587B (zh) * | 2019-04-29 | 2023-05-23 | 浙江大学城市学院 | 一种隧道内集中荷载引起既有隧道位移变形的计算方法 |
| CN111829477B (zh) * | 2020-07-22 | 2022-11-25 | 湖北文理学院 | 轮毂轴承负游隙的计算方法、装置、存储介质及设备 |
-
2007
- 2007-02-06 JP JP2007026577A patent/JP5046671B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2008191025A (ja) | 2008-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5046671B2 (ja) | 共形接触下の接触面圧および表面下応力計算方法 | |
| Mirzaei et al. | Direct meshless local Petrov–Galerkin (DMLPG) method: a generalized MLS approximation | |
| Cao et al. | Laplace ℓ1 Huber based cubature Kalman filter for attitude estimation of small satellite | |
| Hees et al. | Combined Solar system and rotation curve constraints on MOND | |
| Elishakoff | Uncertain buckling: its past, present and future | |
| Samuel et al. | Evaluation of sphericity error from form data using computational geometric techniques | |
| JP2016224912A (ja) | 軸方向圧力補強リブ円筒殻における開口補強方法 | |
| US10060743B2 (en) | Method and apparatus of correcting output value of geomagnetic sensor | |
| JP5026225B2 (ja) | 共形接触下の接触面圧および表面下応力計算方法・計算装置 | |
| Bazzi et al. | Stability and predictability in human control of complex objects | |
| Koop et al. | Passive dynamic biped walking—Part II: stability analysis of the passive dynamic gait | |
| Gunia et al. | The influence of the geometrical parameters for stress distribution in wire raceway slewing bearing | |
| Inoue et al. | A new algorithm for continuous area cartogram construction with triangulation of regions and restriction on bearing changes of edges | |
| Banh et al. | Stress-driven design of incompressible multi-materials under frequency constraints | |
| Samuel et al. | Evaluation of sphericity error from coordinate measurement data using computational geometric techniques | |
| Park et al. | A four-node shell element with enhanced bending performance for springback analysis | |
| WO2023150352A1 (en) | Methods for analyses of internal structures and defects in materials using physics-informed neural networks | |
| Zhang et al. | An approach to analyze the position and orientation between two parts assembled by non-ideal planes | |
| Wong et al. | Study of the discrete shear gap technique in timoshenko beam elements | |
| Romero | On the development of the direct flux reconstruction scheme for high-order fluid flow simulations | |
| Malikov et al. | Analytical and numerical calculation of the force and power requirements for air bending of structured sheet metals | |
| Wu et al. | MAV quaternion attitude determination for accelerometer-magnetometer combination: Internal analysis | |
| CN104778325B (zh) | 基于表面单元的面载荷处理方法及装置 | |
| Zimmerling et al. | Optimisation of layup type and fibre orientation in continuous-fibre reinforced components via anisotropy analysis | |
| Cui et al. | An innovative investigation on the workpiece kinematics and its roundness generation in through-feed centreless grinding |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100127 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120717 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120717 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150727 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5046671 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |