特許 7538671 粘性モデルの決定方法、見かけ粘度の推定方法、およびトルクの推定方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-08-14

(45)【発行日】2024-08-22

(54)【発明の名称】粘性モデルの決定方法、見かけ粘度の推定方法、およびトルクの推定方法

(51)【国際特許分類】

   G01N 11/00 20060101AFI20240815BHJP

   G01N 19/00 20060101ALI20240815BHJP

【ＦＩ】

G01N11/00 A

G01N11/00 C

G01N19/00 A

【請求項の数】 6

(21)【出願番号】P 2020159064

(22)【出願日】2020-09-23

(65)【公開番号】P2022052592

(43)【公開日】2022-04-04

【審査請求日】2023-03-24

【前置審査】

(73)【特許権者】

【識別番号】000102692

【氏名又は名称】ＮＴＮ株式会社

(74)【代理人】

【識別番号】100174090

【弁理士】

【氏名又は名称】和気 光

(74)【代理人】

【識別番号】100100251

【弁理士】

【氏名又は名称】和気 操

(74)【代理人】

【識別番号】100205383

【弁理士】

【氏名又は名称】寺本 諭史

(72)【発明者】

【氏名】小畑 智彦

(72)【発明者】

【氏名】吉野 真人

【審査官】森口 正治

(56)【参考文献】

【文献】特開２０１８－１５５５６４（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０１Ｎ １１／００

Ｇ０１Ｎ １９／００

(57)【特許請求の範囲】

【請求項1】

非ニュートン流体として、基油と増ちょう剤を含むグリースの見かけ粘度を推定する粘性モデルの決定方法であって、
せん断速度と粘度との関係式を示す推定モデルを得る工程（ａ）と、
前記非ニュートン流体を固定板と可動板との間に挟み、前記可動板を前記固定板から離間するように移動させて前記非ニュートン流体を伸ばす引張試験において、前記非ニュートン流体の変形時の引張荷重を測定する工程（ｂ）と、
前記工程（ａ）で得た前記推定モデルに対して数値流体解析を行い、前記引張試験における前記非ニュートン流体の変形時の引張荷重を予測する工程（ｃ）と、
前記工程（ｂ）で測定した前記引張荷重と、前記工程（ｃ）で予測した前記引張荷重とを比較して、これらの荷重が一致するように前記推定モデルのパラメータを調整して粘性モデルを決定する工程（ｄ）とを有し、
前記推定モデルは、Ｃｒｏｓｓ Ｐｏｗｅｒｌａｗモデル、Ｈｅｒｃｈｅｌ－Ｂｕｌｋｌｅｙモデル、Ｐａｐａｎａｓｔａｓｉｏｕモデル、または下記式（１）のモデルに基づくものであって、
前記工程（ａ）は、回転式レオメータを用いて所定のせん断速度における前記非ニュートン流体の粘度を測定し、その測定結果に基づいて、前記推定モデルを得ることを特徴とする粘性モデルの決定方法。

【数1】

ただし、式中の記号は、η：見かけ粘度［Ｐａ・ｓ］、η_ｏｉｌ：流体成分の粘度［Ｐａ・ｓ］、τ_０：降伏応力［Ｐａ］、γ：せん断速度［ｓ^－１］、ｎ：定数、ｍ：定数である。

【請求項2】

前記工程（ｂ）において、前記引張試験は、前記非ニュートン流体を、前記固定板と前記可動板の間に初期厚さ１ｍｍ、直径７．４ｍｍで挟み、その状態から前記可動板を０．００５ｍ／ｓの速度で引き下げる条件で実施され、前記可動板の移動距離が０．０１ｍｍ～１．０ｍｍにおける引張荷重が測定されることを特徴とする請求項１記載の粘性モデルの決定方法。

【請求項3】

前記工程（ｄ）は、前記工程（ｂ）で測定した前記引張荷重を中心とした所定範囲内に、前記工程（ｃ）で予測した前記引張荷重が含まれる場合に前記推定モデルを前記粘性モデルに決定し、含まれない場合に前記推定モデルの前記パラメータを調整することを特徴とする請求項１または請求項２記載の粘性モデルの決定方法。

【請求項4】

非ニュートン流体の見かけ粘度を推定する方法であって、
請求項１から請求項３までのいずれか１項記載の粘性モデルの決定方法によって決定された粘性モデルを用いて、１０ｓ^－１以下の任意のせん断速度における前記非ニュートン流体の見かけ粘度を推定することを特徴とする見かけ粘度の推定方法。

【請求項5】

機械要素に発生する非ニュートン流体のせん断によるトルクを推定する方法であって、
請求項４記載の推定方法によって推定された前記非ニュートン流体の見かけ粘度を用いて、前記任意のせん断速度で発生するトルクを推定する、または、複数の非ニュートン流体間において前記任意のせん断速度で発生するトルクの大小を推定することを特徴とするトルクの推定方法。

【請求項6】

前記機械要素が転がり軸受であることを特徴とする請求項５記載のトルクの推定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、粘性モデルの決定方法、見かけ粘度の推定方法、およびトルクの推定方法に関し、特に、転がり軸受に封入されるグリースの見かけ粘度を推定する粘性モデルの決定方法、見かけ粘度の推定方法、および転がり軸受のトルクの推定方法に関する。

【背景技術】

【0002】

転がり軸受の内部には、転がり摩擦や滑り摩擦の軽減などを目的として、潤滑用のグリースが封入されている。グリースを封入してなるグリース封入軸受は、長寿命で外部の潤滑ユニットなどが不要かつ安価であるため、自動車や産業用機器などの汎用用途によく利用される。ここで、軸受に封入されたグリースは、せん断作用を受けて流動状態になり潤滑に寄与する。

【0003】

このような軸受の潤滑状態として、チャーニング期とチャンネリング期がある。チャーニング期は、グリースが転動体や保持器によって大きく撹拌される状態であり、定常状態でのトルクが高くなる傾向がある。一方、チャンネリング期は、グリースがシールや保持器などに付着し、ほぼ撹拌を受けない状態であり、定常状態でのトルクが低くなる傾向がある。チャーニング期のトルクを算出するためには、転動体と保持器の間や、転動体と軌道面の間でせん断を受けるグリースの粘度（見かけ粘度）を知る必要がある。

【0004】

ここで、グリースは、チキソトロピー性を有する非ニュートン流体である。このような非ニュートン流体の見かけ粘度とせん断速度との関係は、レオメータを用いて測定することができ、これを関数として表す方法として、Ｃｒｏｓｓ Ｐｏｗｅｒｌａｗモデル（非特許文献１）などの計算モデルが知られている。

【先行技術文献】

【非特許文献】

【0005】

【文献】Ｍ．Ｍ．Ｃｒｏｓｓ、「Ｒｈｅｏｌｏｇｙ ｏｆ ｎｏｎ－ｎｅｗｔｏｎｉａｎ ｆｌｕｉｄｓ：ａ ｎｅｗ ｆｌｏｗ ｅｑｕａｔｉｏｎ ｆｏｒ ｐｓｅｕｄｏｐｌａｓｔｉｃ ｓｙｓｔｅｍｓ」、Ｊ．Ｃｏｌｌｏｉｄ Ｓｃｉ.、１９６５、２０、ｐ．４１７－４３７

【発明の概要】

【発明が解決しようとする課題】

【0006】

ところで、軸受内で発生しうるせん断速度域（０．００１～１００００００ｓ^－１）は非常に広域であり、軸受内のせん断速度ごとのグリースの見かけ粘度を測定することは困難である。一般に、回転式レオメータを用いた見かけ粘度の測定は、軸受内の上記せん断速度域に対して、中速域（約１０～１００００ｓ^－１）のみと限定的である。高速域（例えば１０００ｓ^－１以上）では、遠心力により実験装置からグリースがはみ出すため、見かけ粘度の測定が困難である。

【0007】

一方、低速域（例えば１０ｓ^－１以下）では、グリースの見かけ粘度が増ちょう剤の配向などの影響を受けるためか安定しない。図８には、推定モデルとレオメータの実測値との関係を示す。図８のプロットは、回転式レオメータで測定した見かけ粘度であり、実線および破線は、回転式レオメータの結果に一致するように計算モデルのパラメータを調整した推定モデルである。推定モデル１および２は、せん断速度が約１０～１００００ｓ^－１では、両者ともにレオメータの実測値とよく一致しているが、せん断速度１０ｓ^－１以下では、両者の見かけ粘度に大きな差がある。このように、回転式レオメータによる測定では、低せん断速度域の正確なグリースモデルを得ることが困難である。また、低せん断速度域に関する知見は知られていない。そのため、例えば、グリースが封入された転がり軸受において、正確に軸受トルクを推定することや軸受内のグリース挙動をシミュレートすることが困難であった。

【0008】

本発明はこのような事情に鑑みてなされたものであり、低せん断速度域において、非ニュートン流体の見かけ粘度を精度よく推定するための粘性モデルの決定方法、その粘性モデルを用いた見かけ粘度の推定方法、およびトルクの推定方法を提供することを目的とする。

【課題を解決するための手段】

【0009】

本発明の粘性モデルの決定方法は、非ニュートン流体の見かけ粘度を推定する粘性モデルの決定方法であって、せん断速度と粘度との関係式を示す推定モデルを得る工程（ａ）と、上記非ニュートン流体を固定板と可動板との間に挟み、上記可動板を上記固定板から離間するように移動させて上記非ニュートン流体を伸ばす引張試験において、上記非ニュートン流体の変形時の引張荷重を測定する工程（ｂ）と、上記工程（ａ）で得た上記推定モデルに対して数値流体解析を行い、上記引張試験における上記非ニュートン流体の変形時の引張荷重を予測する工程（ｃ）と、上記工程（ｂ）で測定した上記引張荷重と、上記工程（ｃ）で予測した上記引張荷重とを比較して、これらの荷重が一致するように上記推定モデルのパラメータを調整して粘性モデルを決定する工程（ｄ）とを有することを特徴とする。

【0010】

本発明において、「荷重が一致する」とは、両者が完全に一致するだけでなく、ほぼ一致するものも含む。例えば、工程（ｂ）で測定した引張荷重を中心とした所定範囲内に、工程（ｃ）で予測した引張荷重が含まれる場合に一致するとすることができる。

【0011】

上記推定モデルは、Ｃｒｏｓｓ Ｐｏｗｅｒｌａｗモデル、Ｈｅｒｃｈｅｌ－Ｂｕｌｋｌｅｙモデル、Ｐａｐａｎａｓｔａｓｉｏｕモデル、または下記式（１）のモデルに基づくことを特徴とする。

【数1】

ただし、式中の記号は、η：見かけ粘度［Ｐａ・ｓ］、η_ｏｉｌ：流体成分の粘度［Ｐａ・ｓ］、τ_０：降伏応力［Ｐａ］、γ：せん断速度［ｓ^－１］、ｎ：定数、ｍ：定数である。

【0012】

上記工程（ａ）は、回転式レオメータを用いて所定のせん断速度における上記非ニュートン流体の粘度を測定し、その測定結果に基づいて、上記推定モデルを得ることを特徴とする。

【0013】

上記工程（ｂ）において、上記引張試験は、上記非ニュートン流体を、上記固定板と上記可動板の間に初期厚さ１ｍｍ、直径７．４ｍｍで挟み、その状態から上記可動板を０．００５ｍ／ｓの速度で引き下げる条件で実施され、上記可動板の移動距離が０．０１ｍｍ～１．０ｍｍにおける引張荷重が測定されることを特徴とする。

【0014】

上記工程（ｄ）は、上記工程（ｂ）で測定した上記引張荷重を中心とした所定範囲内に、上記工程（ｃ）で予測した上記引張荷重が含まれる場合に上記推定モデルを上記粘性モデルに決定し、含まれない場合に上記推定モデルの上記パラメータを調整することを特徴とする。

【0015】

上記非ニュートン流体は、基油と増ちょう剤を含むグリースであることを特徴とする。

【0016】

本発明の見かけ粘度の推定方法は、非ニュートン流体の見かけ粘度を推定する方法であって、本発明の決定方法によって決定された粘性モデルを用いて、１０ｓ^－１以下の任意のせん断速度における上記非ニュートン流体の見かけ粘度を推定することを特徴とする。

【0017】

本発明のトルクの推定方法は、機械要素に発生する非ニュートン流体のせん断によるトルクを推定する方法であって、本発明の見かけ粘度の推定方法によって推定された上記非ニュートン流体の見かけ粘度を用いて、上記任意のせん断速度で発生するトルクを推定する、または、複数の非ニュートン流体間において上記任意のせん断速度で発生するトルクの大小を推定することを特徴とする。

【0018】

上記機械要素が転がり軸受であることを特徴とする。

【発明の効果】

【0019】

本発明の粘性モデルの決定方法は、上記工程（ａ）～工程（ｄ）を有しており、せん断速度と粘度との関係式を示す推定モデルを、実験と解析の合わせこみによって補正するので、精度に優れた粘性モデルになる。特に、引張試験は低速で非ニュートン流体を引張る試験であり、その試験結果は低速せん断域での非ニュートン流体の挙動を反映していると考えられる。そのため、引張試験での非ニュートン流体の変形時の引張荷重に一致するように推定モデルのパラメータを調整することで、低速せん断域での見かけ粘度の推定に適した粘性モデルになる。その結果、非ニュートン流体の低速せん断域での見かけ粘度を精度良く推定でき、また、トルク推定やグリースの挙動解析に利用できる。

【0020】

本発明の見かけ粘度の推定方法は、本発明の決定方法によって決定された粘性モデルを用いて、１０ｓ^－１以下の任意のせん断速度における非ニュートン流体の見かけ粘度を推定するので、レオメータで測定できる範囲よりも低速せん断域の見かけ粘度を精度よく推定できる。

【0021】

本発明のトルクの推定方法は、本発明の推定方法によって推定された非ニュートン流体の見かけ粘度を用いて、任意のせん断速度における非ニュートン流体のトルクを推定する、または、任意のせん断速度における複数の非ニュートン流体間のトルクの大小を推定するので、トルクの推定をより精度よく行うことができる。

【図面の簡単な説明】

【0022】

【図1】グリースが封入された転がり軸受を示す図である。

【図2】本発明に係る粘性モデルの決定方法の概略を示す工程図である。

【図3】レオメータの一例を示す概要図である。

【図4】引張試験機の概要を示す図である。

【図5】引張試験におけるグリースの伸びを示す図である。

【図6】実施例と比較例の見かけ粘度特性を示す図である。

【図7】引張試験における移動距離と引張荷重の関係を示す図である。

【図8】推定モデルとレオメータの実測値との関係を示す図である。

【発明を実施するための形態】

【0023】

本発明における粘性モデルは、相対運動する２部材間の隙間に存在する非ニュートン流体の見かけ粘度を推定する粘性モデルである。ここで、非ニュートン流体とは、粘度が「せん断速度」に依存する流体であり、一般に非ニュートン流体の粘度は、「ずり応力」を「せん断速度」で除した見かけ粘度で表される。本発明において、非ニュートン流体には、擬塑性流体、ビンガム流体、ダイラタント流体などが用いられる。これらの中でも、擬塑性流体が好ましい。擬塑性流体として具体的には、転がり軸受に封入されるグリースや潤滑油が挙げられる。なお、以下には、非ニュートン流体としてグリースを用いた形態について説明する。

【0024】

転がり軸受に封入されるグリースは、基油と増ちょう剤とを含む。基油は、特に限定されず、通常グリースの分野で使用される一般的なものを使用できる。例えば、高度精製油、鉱油、エステル油、エーテル油、合成炭化水素油（ＰＡＯ油）、シリコーン油、フッ素油、およびこれらの混合油などを使用できる。

【0025】

基油の４０℃における動粘度としては、特に限定されないが、２０ｍｍ^２／ｓ～４００ｍｍ^２／ｓが好ましい。なお、基油として混合油を用いる場合は、該混合油の動粘度がこの範囲内であることが好ましい。

【0026】

増ちょう剤は、特に限定されず、通常グリースの分野で使用される一般的なものを使用できる。例えば、金属石けん、複合金属石けんなどの石けん系増ちょう剤、ベントン、シリカゲル、ウレア化合物、ウレア・ウレタン化合物などの非石けん系増ちょう剤を使用できる。金属石けんとしては、ナトリウム石けん、カルシウム石けん、アルミニウム石けん、リチウム石けんなどが、ウレア化合物、ウレア・ウレタン化合物としては、ジウレア化合物、トリウレア化合物、テトラウレア化合物、他のポリウレア化合物、ジウレタン化合物などが挙げられる。

【0027】

また、グリースには、必要に応じて他の公知の添加物を含有させることができる。この添加物としては、アミン系やフェノール系の酸化防止剤、塩素系、イオウ系、りん系化合物、有機モリブデンなどの極圧剤、石油スルホネート、ジノニルナフタレンスルホネート、ソルビタンエステルなどのさび止剤などが挙げられる。

【0028】

グリースが封入された転がり軸受の一例を図１に基づいて説明する。図１は深溝玉軸受の断面図である。転がり軸受１は、外周面に内輪軌道面２ａを有する内輪２と内周面に外輪軌道面３ａを有する外輪３とが同心に配置され、内輪軌道面２ａと外輪軌道面３ａとの間に複数個の玉４が配置される。この玉４は、保持器５により保持される。また、内・外輪の軸方向両端開口部８ａ、８ｂがシール部材６によりシールされ、少なくとも玉４の周囲にグリース７が封入される。内輪２、外輪３および玉４は鉄系金属材料からなり、グリース７が玉４との軌道面に介在して潤滑される。

【0029】

本発明における粘性モデルは、例えば、所定のせん断速度におけるグリースの見かけ粘度の推定に用いられる。本発明に係る粘性モデルの決定方法は、図２に示す（ａ）～（ｄ）の４つの工程を少なくとも備える。すなわち、（ａ）推定モデル２を取得する取得工程と、（ｂ）引張試験においてグリースの変形時の引張荷重３を測定する測定工程と、（ｃ）推定モデル２に対して解析を行い、引張試験におけるグリースの変形時の引張荷重４を予測する解析工程と、（ｄ）引張荷重３と引張荷重４とを比較して、これらの荷重が一致するように推定モデル２のパラメータを調整して粘性モデル５を決定する決定工程とを備える。以下に、各工程について説明する。

【0030】

＜工程（ａ）＞
この工程は、グリースの見かけ粘度を推定する推定モデル２を取得する工程である。図２において、工程（ａ）は、任意の計算モデル１を取得する一次取得工程（ａ１）と、回転式レオメータを用いたレオロジー測定結果に基づき、計算モデル１から推定モデル２を取得する二次取得工程とを有する。なお、工程（ａ）は、必ずしも回転式レオメータを用いた測定（実測）を行う必要はない。例えば、工程（ａ）において、経験則などに基づき、予めパラメータに所定の値が代入された推定モデルを取得してもよい。なお、以下では、図２に示す形態の工程（ａ）を説明する。

【0031】

グリースは非ニュートン流体であり、半固体状の高粘性流体でもある。グリースの粘度は、せん断速度の増加に対して減少する。この特性を表現する計算モデル１として、工程（ａ１）では、Ｃｒｏｓｓ Ｐｏｗｅｒｌａｗモデル（下記式（２））や、Ｈｅｒｃｈｅｌ－Ｂｕｌｋｌｅｙモデル（下記式（３））、Ｐａｐａｎａｓｔａｓｉｏｕモデル（下記式（４））などを取得する。なお、下記式（２）～（４）中のηはグリースの見かけ粘度を表す。

【0032】

【数2】

【0033】

ただし、式中の記号は、η_０：ゼロせん断粘度（せん断速度γ＝０での見かけ粘度）［Ｐａ・ｓ］、η_ｏｉｌ：基油粘度［Ｐａ・ｓ］、γ：せん断速度［ｓ^－１］、ｍ：定数、ｎ：定数である。上記式（２）では、せん断速度が無限大の時、ηはη_ｏｉｌに漸近すると仮定している。ｎは、せん断速度の変化に対する見かけ粘度の変化を表しており、工程（ａ２）のレオロジー測定によって決定される。

【0034】

【数3】

【0035】

ただし、式中の記号は、τ_０：降伏応力［Ｐａ］、γ：せん断速度［ｓ^－１］、Ｋ：定数、ｎ：定数である。

【0036】

【数4】

【0037】

ただし、式中の記号は、η_ｏｉｌ：流体成分の粘度［Ｐａ・ｓ］、τ_０：降伏応力［Ｐａ］、γ：せん断速度［ｓ^－１］、ｍ：定数である。

【0038】

工程（ａ１）で取得される計算モデル１は、上記のような公知の計算モデルに限られない。例えば、計算モデル１として、下記式（１）のモデルを用いることができる。下記式（１）は、高速せん断域での実験データに基づいてＰａｐａｎａｓｔａｓｉｏｕモデルを改良した式である。

【数5】

【0039】

ただし、式中の記号は、η_ｏｉｌ：流体成分の粘度［Ｐａ・ｓ］、τ_０：降伏応力［Ｐａ］、γ：せん断速度［ｓ^－１］、ｎ：定数、ｍ：定数である。

【0040】

グリースの場合、流体成分の粘度η_ｏｉｌには、実温度（例えば軸受内部の温度や機械装置全体の温度）の基油の動粘度に対して基油の密度を掛けた値（粘度）を用いる。また、グリースに粘度向上剤が含まれる場合は、流体成分の粘度η_ｏｉｌには、基油および粘度向上剤の粘度を用いる。

【0041】

続く工程（ａ２）では、計算モデル１における降伏応力や各定数などのパラメータを求めて推定モデル２を取得する。パラメータは、回転式レオメータを用いたグリースのレオロジー特性の評価などに基づき特定できる。グリースのレオロジー特性の評価では、レオメータとして、コーンプレート型のセルを有するものを用いることが好ましい。このようなレオメータの概要を図３に示す。図３に示すように、回転式レオメータ１１は、コーンプレート型のセル１２と、水平円盤プレート１３とから構成されており、セル１２とプレート１３とは１点で接する（僅かなギャップあり）ように配置され、これらの間に試料であるグリース１４を配置する。このレオメータでは、グリース１４に加わるせん断速度が、セル中心からの距離に依存せずに、どの位置においても同一となる。レオロジー測定の条件としては、（１）一定温度・一定方向回転での回転速度依存性、（２）一定温度・一定せん断ひずみにおける振動周波数依存性、（３）一定周波数における動的粘弾性のせん断応力依存性などがあるが、本発明では主に（１）の条件で測定を行なう。

【0042】

具体的なレオロジー測定条件としては、回転式レオメータ（Ｔｈｅｒｍｏ Ｆｉｓｈｅｒ Ｓｃｉｅｎｔｉｆｉｃ社製ＨＡＡＫＥ ＲｈｅｏＷｉｎ ＭＡＲＳ１）に、直径２０ｍｍ、先端角度１７８°のコーンプレート型のセルを用い、一定温度・一定方向回転で、温度２０℃などで行なう。この場合、せん断速度を増速し、各せん断速度で定常状態になった時の粘度を測定する。例えばＣｒｏｓｓ Ｐｏｗｅｒｌａｗモデルの場合には、このレオロジー測定結果からη_０、ｍ、ｎが求められる。求めた各値を上記式（２）に代入して、推定モデル２が得られる。この時点でｎは決定され、η_０、ｍは後続の工程を実施することで同定される。

【0043】

＜工程（ｂ）＞
この工程は、引張試験においてグリースの変形時の引張荷重３を測定する工程である。引張試験では、所定の面積および所定の間隔に調整した平行な平板間（２面間）にグリースを充填し、平板間の距離を広げてグリースを破断させる。

【0044】

図４を用いて引張試験について説明する。図４には、引張試験機の概略図を示す。引張試験機２１は、天壁に固定された円柱状の固定板２２と、円柱状の可動板２３と、固定板２２に対して可動板２３を上下に移動させる駆動装置２４とを備える。固定板２２の下面と可動板２３の上面は、互いに平行な面である。駆動装置２４は、可動板２３が載置されるステージ２５と、ステージ２５の移動を規制するリニアガイド２６と、モータ２７と、モータ２７の駆動によって回転するボールねじ２８と備える。ステージ２５は、ボールねじ２８のねじ溝に螺合しており、ボールねじ２８の回転によって軸上を摺動しながら相対的に移動する。また、固定板２２には、グリースの変形時に固定端側にかかる荷重を計測するロードセル３０が設けられている。

【0045】

以下に、引張試験の手順について説明する。

【0046】

（１）まず、可動板２３の上面に測定対象のグリース２９を所定量充填する。グリース２９の充填量は、例えば１ｃｍ^３～５ｃｍ^３である。

【0047】

（２）次に、グリース２９が固定板２２の下面に接触する位置まで可動板２３を上昇させる。このとき、グリース２９が所定の初期厚さＡ（図５（ａ）参照）となるように、可動板２３の位置を調整する。図５（ａ）の初期状態において、固定板２２と可動板２３の間隔が、グリース２９の初期厚さＡになる。初期厚さＡは、例えば０．５ｍｍ～２ｍｍに設定され、具体的には１ｍｍに設定される。この初期状態のグリース２９の直径φは、例えば５ｍｍ～１０ｍｍに設定され、具体的には７．４ｍｍに設定される。図５（ａ）に示すように、グリース２９の直径φと同じ径サイズの固定板２２および可動板２３を用いて、固定板２２および可動板２３の全体にグリース２９が拡がるように可動板２３の位置を調整してもよい。なお、固定板２２および可動板２３の径サイズや形状は、図５（ａ）に限らない。例えば、固定板２２や可動板２３を円柱状以外の形状とし、一方の対向面に円形状のグリース充填領域を設けてもよい。

【0048】

（３）次に、図４の黒矢印で示すように、可動板２３を所定速度で下降させる。可動板２３の下降速度は、例えば０．０００１ｍ／ｓ～０．１ｍ／ｓに設定され、具体的には０．００５ｍ／ｓに設定される。図４に示すように、可動板２３の下降によってグリース２９が変形する。工程（ｂ）では、可動板２３の移動に伴うグリース２９の変形時の引張荷重を経時的に測定し、可動板２３の任意の移動距離における引張荷重を測定する。なお、後続の工程（ｄ）において比較の対象となる引張荷重は、可動板２３の移動初期（例えば１秒以内）の引張荷重であることが好ましい。また、可動板２３の移動距離で言えば、好ましくは０．０１ｍｍ～１．０ｍｍ、より好ましくは０．０１ｍｍ～０．５ｍｍ、さらに好ましくは０．０５ｍｍ～０．３ｍｍの範囲内の移動距離における引張荷重で比較する。

【0049】

可動板２３の移動距離は、図５に示すように、初期厚さＡと、固定板２２と移動時の可動板２３との間隔Ｂを用いて、下記の式により算出される。
可動板２３の移動距離＝Ｂ－Ａ

【0050】

（４）図４に戻り、可動板２３の下降が更に進むと、最終的にグリース２９が破断する。グリース２９の破断の有無は、目視や高速度カメラや引張荷重を計測することによって確認できる。

【0051】

なお、図４に示す引張試験機の構成は、これに限らない。例えば、可動板２３をステージ２５と一体に形成してもよい。また、可動板２３を鉛直下向きに引き下げる構成に代えて、可動板を鉛直上向きに引き上げる構成にしてもよい。

【0052】

＜工程（ｃ）＞
この工程では、数値流体解析を行う。数値流体解析には、有限体積法によって数値解析を行う汎用の熱流体解析ソフトウェアなどが用いられる。この解析では、所定の流体物性のグリースを上述の引張試験に適用した場合、引張荷重がどのように変化するかを再現する。解析の際の計算条件には、工程（ｂ）で実施した引張試験の試験条件などを用いる。

【0053】

また、解析に適用するグリースの流体物性は、係数の数を少なくできることから、回転式レオメータの実測値から求められる推定モデルを用いることが好ましい。なお、経験則などに基づいて取得された推定モデルを用いてもよい。後述の実施例で示すように、この数値流体解析によって、引張試験におけるグリースの変形時の引張荷重４がシミュレーションされる。この場合、工程（ｂ）で取得された可動板の移動距離に対応した引張荷重がシミュレーションされる。

【0054】

＜工程（ｄ）＞
この工程では、まず、工程（ｂ）で測定した可動板の所定の移動距離における引張荷重３と、工程（ｃ）でシミュレーションした可動板の当該移動距離における引張荷重４とを比較して、これらの荷重が一致しているか否かを判定する。具体的には、各引張荷重が所定範囲内（例えば、±１０％）であるか否を判定する。そして、各引張荷重が所定範囲内である場合に一致したと判定し、各引張荷重が所定範囲内でない場合に一致していないと判定する。例えば、引張荷重３を中心とした所定範囲内（例えば、±１０％）に、引張荷重４が含まれる場合に一致したと判定し、含まれない場合に一致していないと判定する。なお、比較する引張荷重は可動板の１点の移動距離に限らず、複数点の移動距離における引張荷重を用いて比較してもよい。複数の移動距離における引張荷重を用いる場合、両者の各移動距離における引張荷重がすべて所定範囲内の場合に、引張荷重３と引張荷重４が一致したと判定する。

【0055】

両者の荷重が一致すると判定した場合は、工程（ａ）で取得した推定モデル２は正確であるとして、該推定モデル２を粘性モデル５に決定する。一方、両者の荷重が一致していないと判定した場合は、推定モデル２のパラメータを調整する。例えば、Ｃｒｏｓｓ Ｐｏｗｅｒｌａｗモデルの場合には、パラメータとして係数η_０とｍを変更する。なお、Ｈｅｒｃｈｅｌ－Ｂｕｌｋｌｅｙモデルの場合には、τ_０とＫ、Ｐａｐａｎａｓｔａｓｉｏｕモデルの場合には、τ_０とｍ、上記式（１）のモデルの場合には、τ_０とｍを変更する。

【0056】

その後、数値流体解析（工程（ｃ））を再度行い、その解析で得られた引張荷重４と引張荷重３を比較して、これらの荷重が一致しているか否かを再度判定する。そして、２つの引張荷重が一致するまで、推定モデル２のパラメータの変更を繰り返す。２つの引張荷重が一致する場合には、推定モデル２で推定した粘度特性が、実際のグリースの粘度を表せているといえる。

【0057】

このように、引張試験の解析と実験によりそれぞれ得られたグリースの引張荷重を比較することで、推定モデルの妥当性を検証することができる。そして、引張試験時の引張荷重を実験と解析で一致させることで、実際のグリースの粘度をより反映した精度の高い粘性モデルになる。

【0058】

本発明に係る見かけ粘度の推定方法は、上記の決定方法によって決定された粘性モデルを用いて、任意のせん断速度における見かけ粘度を推定する方法である。せん断速度は、チャーニング期の回転中にポケット隙間に存在するグリースに掛かるせん断速度であり、設定する軸受回転数などから算出できる。例えば、玉が保持器のポケット部中心に位置すると仮定すれば、６２０４軸受で内輪回転速度１０^４ｍｉｎ^－１の場合、ポケット部のグリースのせん断速度は１０^５ｓ^－１オーダーに相当する。任意のせん断速度を得られた粘性モデルに代入することで、チャーニング期の回転中にポケット隙間に存在するグリースのそのせん断速度における見かけ粘度を、予測値として求めることができる。

【0059】

ここで、一般的にグリースの見かけ粘度が大きくなるほど、軸受トルクも大きくなる。そのため、上述の方法で推定した見かけ粘度を用いることで、軸受トルクを推定する、または、複数のグリース間で軸受トルクの大小を推定することができる。軸受トルクの推定は、例えば、軸受を含む機械装置全体の温度を軸受内部の温度と仮定しグリースの見かけ粘度を算出することで行える。これにより、要求される軸受トルクからグリースの選定（レオメータで評価）を行なうことも可能となる。また、複数のグリース間の軸受トルクの大小の推定は、各グリースの見かけ粘度の大小を比較することで行える。例えば、複数のグリースの中から最も見かけ粘度が低いグリースを、軸受トルクが最も低いグリースとして選定することができる。

【0060】

図１では、グリースが封入される転がり軸受として玉軸受について例示したが、玉軸受以外にも、円筒ころ軸受、円すいころ軸受、自動調心ころ軸受、針状ころ軸受、スラスト円筒ころ軸受、スラスト円すいころ軸受、スラスト針状ころ軸受、スラスト自動調心ころ軸受などにも適用できる。

【0061】

上記では、グリースおよび該グリースが封入された転がり軸受について述べたが、本発明の推定方法の対象となる非ニュートン流体は、グリースに限らず、ゲルやゾル、その他のチキソトロピー性を有する流体や半固体であってもよい。また、非ニュートン流体がせん断を受ける相対運動する２部材は、転がり軸受の軸受部材に限らず、他の機械要素に適用してもよい。他の機械要素としては、例えば、リニアガイドや、クラッチ、カム、ジョイント、チェーン、歯車、滑り軸受などが挙げられる。

【0062】

本発明の推定方法は、低速せん断域でも見かけ粘度を精度よく推定できることから、１０ｓ^－１以下のせん断速度域、より好ましくは１ｓ^－１以下のせん断速度領域で使用される非ニュートン流体の推定に適している。具体的には、低速回転用の転がり軸受に封入されるグリースの見かけ粘度の推定に適している。なお、本発明における粘性モデルによれば、高速せん断域での見かけ粘度も十分推定することができる。

【実施例】

【0063】

グリースについて、レオロジー測定を行なった。レオメータ（Ｔｈｅｒｍｏ Ｆｉｓｈｅｒ Ｓｃｉｅｎｔｉｆｉｃ社製ＨＡＡＫＥ ＲｈｅｏＷｉｎ ＭＡＲＳ１）において、直径２０ｍｍ、先端角度１７８°のコーンプレート型のセルを用いて測定した。せん断速度を１０から５０００（単位：１/ｓ）まで増速し、各せん断速度で定常状態になった時の見かけ粘度を計測した。図６には、その実測値をプロットした。

【0064】

計算モデルとしては、上記式（２）に示すＣｒｏｓｓ Ｐｏｗｅｒｌａｗモデルを用いた。上述のレオロジー測定結果より、上記式（２）における係数（η_０、ｍ、ｎ）を求めて推定モデルを取得した。なお、基油の粘度η_ｏｉｌは、Ｗａｌｔｈｅｒの式によって算出した。

【0065】

得られた各係数を表１に示す。表１には、実施例と比較例の２つの推定モデルの係数の値を示している。なお、基油の粘度η_ｏｉｌおよびｎは、実施例と比較例で同じ値であり、これらの値はこの時点で決定される。

【0066】

【表1】

【0067】

図６には、上記式（２）に表１の各値を代入して得られた近似式を、実線（実施例）と破線（比較例）で示した。図６に示すように、せん断速度が約１０～１００００ｓ^－１では、実施例と比較例のいずれも回転式レオメータで測定した結果と一致している。一方、せん断速度１０ｓ^－１以下では両者に差がみられ、レオロジー測定結果だけでは、どちらの粘度特性が正確であるか判別が困難である。

【0068】

次に、図４および図５で示した引張試験を行ない、グリースの変形時の引張荷重を測定した。試験開始前のグリースの厚さ（初期厚さＡ）は１ｍｍで、直径φは７．４ｍｍであった。また、固定板および可動板には、直径７．４ｍｍの円柱状の平板を用いた。可動板を０．００５ｍ／ｓの速度で引き下げ、可動板の移動距離が０．０５ｍｍ、０．２５ｍｍ、０．５ｍｍの引張荷重を測定した（図７参照）。

【0069】

また、熱流体解析ソフトウェアを用いて数値流体解析を行い、引張試験において可動板の移動距離が０．０５ｍｍ、０．２５ｍｍ、０．５ｍｍの時の引張荷重をシミュレーションした。解析の計算条件には、実際に行った引張試験の試験条件、実施例の推定モデルおよび比較例の推定モデルに基づく流体物性などを用いた。この解析によって得られた引張荷重を図７に示す。

【0070】

図７に示すように、実施例の推定モデルからシミュレーションされた引張荷重は、実験で得られた引張荷重とよく一致した。具体的には、可動板の移動距離が０．０５ｍｍ、０．２５ｍｍ、０．５ｍｍの引張荷重の実験値との誤差がいずれも±１０％以内であった。一方、比較例の推定モデルからシミュレーションされた引張荷重は、実験で得られた引張荷重と一致しなかった。具体的には、可動板の移動距離が０．０５ｍｍ、０．２５ｍｍ、０．５ｍｍの引張荷重の実験値との誤差が２０％以上であった。実際のグリースの粘度と、推定モデルの見かけ粘度が異なっていれば、グリース内に作用する応力も異なる。この場合、比較例の見かけ粘度特性は実際のグリース粘度を表現できておらず、実施例の見かけ粘度特性は実際のグリース粘度を表現できるといえる。

【0071】

また、図７に示すように、引張荷重は、可動板の移動距離が長くなるほど小さくなっている。グリースの変形時の初期に、より大きな応力が働くためと考えられる。また、比較例と実験データの乖離具合は、可動板の移動距離が長くなるほど小さくなっている。つまり、実験データとの比較において、その判定結果に移動距離が影響する可能性がある。そのため、引張荷重の一致・不一致をより正確に判定するため、移動距離が短い領域の引張荷重を用いて判定することが好ましいと考えられる。

【0072】

上記の引張荷重の一致の判定によって、粘性モデルが決定される。例えば、実施例の推定モデルの場合には、その推定モデルを粘性モデルに決定する。一方、比較例の推定モデルの場合には、その推定モデルのパラメータの調製を行なう。具体的には、η_０、ｍの値を大きくしたり小さくしたりして変更する。そして、変更された推定モデルを用いて、数値流体解析によりグリースの変形時の引張荷重を再度シミュレーションし、引張試験の引張荷重と一致するまでパラメータの調製を行い、最終的に粘性モデルを決定する。このように、グリースの見かけ粘度は、回転式レオメータのみの同定では、実際の粘度特性と一致しない場合もあるが、引張試験の結果とを組み合わせることで一致させることができる。

【0073】

このようにして得られた粘度特性を利用すれば、実際の軸受トルクやグリース挙動を正確に推定することができる。例えば、軸受が低速で回転する場合、低速せん断時のグリースの見かけ粘度がトルクの主要因となるため、該見かけ粘度を正確に推定できなければ、推定したトルクが実際のトルクと乖離してしまう。また、軸受内のグリース挙動は、粘度が異なっていれば大きく異なる。これらのことから本発明の決定方法で決定された粘性モデルをトルク推定やグリースの挙動解析などに用いることが望ましい。

【産業上の利用可能性】

【0074】

本発明の決定方法により決定される粘性モデルによれば、低せん断速度域において、非ニュートン流体の見かけ粘度を精度よく推定でき、ひいては非ニュートン流体が存在する状態で２部材を相対運動させる際のトルクを精度よく推定できるので、例えば、低速回転用の転がり軸受のグリース選定などに好適に利用できる。

【符号の説明】

【0075】

１ 転がり軸受
２ 内輪
３ 外輪
４ 玉
５ 保持器
６ シール部材
７ グリース
８ 開口部
１１ レオメータ
１２ コーンプレート型セル
１３ 水平円盤プレート
１４ グリース
２１ 引張試験機
２２ 固定板
２３ 可動板
２４ 駆動装置
２５ ステージ
２６ リニアガイド
２７ モータ
２８ ボールねじ
２９ グリース
３０ ロードセル

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】