JPWO2012157637A1 - 磁力特性算出方法、磁力特性算出装置及びコンピュータプログラム - Google Patents

Abstract

予め記憶してあるδ（深度）—ΔＨｃＪ（保磁力増加量）の対応を示すデータベース及びＤｙ拡散条件（Ｄｙ拡散における拡散係数と拡散束と処理時間）の情報を示すデータベースを用い、磁石の形状情報及びＤｙ導入面情報から磁石内のＤｙ導入量の分布を算出し、Ｄｙ導入量の分布から磁石内でのΔＨｃＪの分布を算出する。保磁力が均質でなく分布を有する磁石について、算出したΔＨｃＪの分布を用いてＪ−Ｈカーブを算出し、温度係数を用いて所定の温度での減磁率を算出する。

Description

本発明は、磁石の情報を記憶しておき演算手段により磁石の磁力特性を算出する方法に関し、特に、磁石表面よりジスプロシウムなどの重希土類元素を内部に拡散させ導入した磁石内部の磁力特性、特に減磁特性を精度よく算出することができる磁力特性算出方法、磁力特性算出装置及びコンピュータプログラムに関する。

高性能な永久磁石として、Ｎｄ−Ｆｅ−Ｂ系焼結磁石は、多様な機器、特にハードディスクドライブ、又は各種モータに使用されている。

Ｎｄ−Ｆｅ−Ｂ系焼結磁石は、高温に晒されたり減磁界が印加された場合、残留磁束密度が減少することがある（減磁）。この減磁には、常温に戻したときに回復する「可逆減磁」と、回復しない「不可逆減磁」とがある。上述したような多様な機器では、使用環境の温度も様々であり、高温下で減磁界が印加された場合でも不可逆減磁が発生しないことが求められる。従来より、Ｎｄ−Ｆｅ−Ｂ系焼結磁石にジスプロシウム（Ｄｙ）などの重希土類元素を添加することによって磁石の主相中のＮｄを重希土類元素で置換し、高い保磁力を得ることで高温下で減磁界が印加された場合でも不可逆減磁が発生しない磁石を製作する技術が利用されている。以降、当明細書では「不可逆減磁」を単に減磁と表現する。

しかしながら、Ｎｄ−Ｆｅ−Ｂ系焼結磁石におけるＮｄを重希土類元素で置換する方法では、保磁力が増加すると残留磁束密度が依然として低下する。また、重希土類元素は希少資源であるため、耐熱性は維持しつつも使用する重希土類元素の量を削減することが望ましい。そこで、ジスプロシウムなどの重希土類元素を拡散源から磁石表面へ供給し、粒界拡散を主体とした熱拡散によって磁石の主相外殻部に重希土類元素の濃度が高い層をつくる技術が提案されている（特許文献１）。特許文献２には、磁石表面からジスプロシウムなどの重希土類元素を磁石内へ拡散させる技術が開示されている。これにより、永久磁石全体の残留磁束密度の低下を抑制しつつ、保磁力を向上させた高性能の永久磁石を製造することができる。

磁石がモータなどに使用される場合、磁石を適切に設計するために、また、モータの特性などの仕様を決定するために、磁石の磁力特性、特に、熱などによる減磁率を知ることが重要である。特許文献３から５には、永久磁石の減磁評価を行なう方法が開示されている。

国際公開２００７／１０２３９１号パンフレット 国際公開２００６／４３３４８号パンフレット 特開２００４−１２７０５６号公報 特開２００４−２５７８７９号公報 国際公開２００８／１２３２５１号パンフレット

ジスプロシウムなどの重希土類元素を表面から導入して主相外殻部又は内部に拡散させたＮｄ−Ｆｅ−Ｂ系焼結磁石の保磁力は磁石の表面近傍にて特に高く、深さ方向に対して不均質である。そして保磁力は温度変化に対して非線形の変化をし、減磁界も各部位で異なるので、各部位で減磁率も異なる。したがって、重希土類元素を磁石表面より主相外殻部に偏在させたＮｄ−Ｆｅ−Ｂ系焼結磁石の減磁率を精度よく求めるためには、磁石内部の保磁力が異なる部位毎に減磁率を得る必要がある。

特許文献３から５に開示されている技術はいずれも、磁石内の部位毎に減磁率を求める方法であるが、初期的な値はバルクで測定された磁石全体の磁束密度、減磁率から算出する構成であった。これに対し、発明者は、拡散処理後の磁石内の保磁力の分布を、ジスプロシウムを拡散させた磁石の既知の保磁力増加量の情報を用いて推定した結果、実測値を精度よく再現できるという知見を得た。重希土類元素を拡散したＮｄ−Ｆｅ−Ｂ系焼結磁石内の減磁前の保磁力分布を推定できれば、その後の磁石内での異なる部分における減磁率を精度良く算出することが可能である。

本発明は斯かる知見に基づきなされたものであり、重希土類元素を表面から内部に拡散させ導入した磁石内部の保磁力増加量の分布、磁石内部の磁力特性、特に減磁特性を精度よく算出することができる磁力特性算出方法、磁力特性算出装置及びコンピュータプログラムを提供することを目的とする。

本発明に係る磁力特性算出方法は、重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を求める方法において、重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を予め記憶しておき、前記磁石の大きさと形状とを示す形状情報を受け付ける第１ステップ、受け付けた形状情報に対応させて、導入面情報を受け付ける第２ステップ、記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する第３ステップ、並びに、算出した導入量分布と、記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する第４ステップを含むことを特徴とする。

本発明に係る磁力特性算出方法は、前記拡散係数は、導入される重希土類元素の濃度依存性の関数で表されてあることを特徴とする。

本発明に係る磁力特性算出方法は、重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを予め記憶しておき、記憶してある磁化曲線と、第４ステップで算出した保磁力増加量の分布とに基づき、前記磁石の各部位の所定の第１温度における磁化曲線を算出する第５ステップ、算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する第６ステップ、及び、第６ステップで算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する第７ステップを更に含むことを特徴とする。

本発明に係る磁力特性算出方法は、第４ステップで算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する第８ステップ、及び、前記磁石の減磁率が所定の率以下となる減磁温度を特定する第９ステップを更に含むことを特徴とする。

本発明に係る磁力特性算出装置は、重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を求める磁力特性算出装置において、重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を記憶しておく記憶手段、前記磁石の大きさと形状とを示す形状情報を受け付ける手段、受け付けた形状情報に対応させて、導入面情報を受け付ける手段、前記記憶手段に記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する手段、並びに、算出した導入量分布と、前記記憶手段に記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する保磁力増加量分布算出手段を備えることを特徴とする。

本発明に係る磁力特性算出装置は、重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを予め記憶しておく手段、記憶してある磁化曲線と、前記保磁力増加量分布算出手段にて算出した保磁力増加量の分布とに基づき、前記磁石の所定の第１温度における磁化曲線を算出する手段、算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する手段、及び、算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する手段を更に備えることを特徴とする。

本発明に係る磁力特性算出装置は、前記保磁力増加量分布算出手段にて算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する手段、及び、前記磁石の減磁率が所定の率以下となる減磁温度を特定する手段を更に備えることを特徴とする。

本発明に係るコンピュータプログラムは、記憶手段を備えるコンピュータに、重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を、前記記憶手段に記憶させてある重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を用いて算出させるコンピュータプログラムであって、コンピュータに、前記磁石の大きさと形状とを示す形状情報を取得する第１ステップ、形状情報に対応させて、導入面情報を取得する第２ステップ、記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する第３ステップ、並びに、算出した導入量分布と、記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する第４ステップを実行させることを特徴とする。

本発明に係るコンピュータプログラムは、重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを記憶してある記憶手段を更に用い、前記コンピュータに、記憶してある磁化曲線と、第４ステップで算出した保磁力増加量の分布とに基づき、前記磁石の各部位の所定の第１温度における磁化曲線を算出する第５ステップ、算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する第６ステップ、及び、第６ステップで算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する第７ステップを更に実行させることを特徴とする。

本発明に係るコンピュータプログラムは、前記コンピュータに、第４ステップで算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する第８ステップ、及び、前記磁石の減磁率が所定の率以下となる減磁温度を特定する第９ステップを更に実行させることを特徴とする。

本発明では、磁力特性を算出する対象である磁石の形状及びジスプロシウムなどの重希土類元素の導入面の情報に対応した拡散条件（拡散係数と拡散束と処理時間）の情報に基づき拡散方程式を用いて、磁石内における重希土類元素の導入量の分布が算出される。算出された導入量の分布から磁石内における重希土類元素の拡散による保磁力増加量の分布が求められる。求められた保磁力増加量の分布を用いることにより、拡散前の均質な保磁力から、拡散後の部位毎に異なる保磁力を精度よく推定することが可能となる。

本発明では、拡散条件の一つである拡散係数は、導入された重希土類元素の濃度をパラメータとする関数で示されており、導入量の分布の算出の際に拡散方程式に使用される。これにより、算出精度を維持しつつも、記憶しておく情報量を削減することができると共に補完処理が不要となる。

本発明では、拡散前の均質な磁石（母材）の磁化曲線と、上述の発明にて求めた保磁力増加量の分布とに基づき、拡散後の磁石の各部位での所定の第１温度（例えば常温）における磁化曲線が求められ、更に、求めた磁化曲線から、異なる保磁力毎に記憶してある温度変化に対する保磁力変化率を示す温度係数の情報に基づき、拡散後の磁石の各部位での所定の第２温度（例えば加熱温度）における磁化曲線が求められ、更に、所定の第２温度で異なる減磁界が印加された後の第１温度（常温）に戻った後の磁石内の各部位の磁化曲線（磁力特性）が求められる。これにより、磁石全体としての減磁率を精度よく算出することが可能となる。

本発明では、拡散前の均質な磁石（母材）の磁化曲線と、上述の発明にて求めた保磁力増加量の分布とに基づき、拡散後の磁石の所定の第１温度（例えば常温）の磁化曲線が精度よく求められ、更に、異なる保磁力毎に記憶してある保磁力の温度変化量を示す温度係数の情報に基づき、異なる複数の第２温度での磁石内における各部位での保磁力の異なる複数の第２温度で異なる減磁界が印加されて減磁した場合の第１温度における磁石全体の減磁率が求められ、所定の率以下となる減磁温度が特定される。これにより、重希土類元素を拡散したＮｄ−Ｆｅ−Ｂ系焼結磁石の磁力特性として、高温下での減磁特性を精度よく算出することが可能である。

本発明による場合、重希土類元素を表面から内部に拡散させ導入した磁石内での部位毎に異なる保磁力の分布を計算により精度よく推定することが可能となるので、磁石全体としての磁力特性を精度よく算出することができる。

特に、モータのような高温下での使用が想定される磁石では、不可逆的な残留磁束密度の低下率所謂不可逆減磁率を精度よく算出することが必要であるが、重希土類元素を拡散させた磁石において拡散及び温度に対する保磁力の変化が分布を有していても減磁率も精度よく算出することができる。

本実施の形態における磁力特性算出装置の構成を示すブロック図である。 本実施の形態における磁力特性算出装置の演算部が実行する磁力特性算出処理の一例を示すフローチャートである。 拡散条件データベースの内容例である濃度と拡散係数との関係を示す説明図である。 拡散係数の同定処理の一例を示すフローチャートである。 Ｄｙ導入量−ΔＨｃＪデータベースの内容例であるＤｙ導入量と保磁力増加量との関係を示す説明図である。 Ｊ−Ｈカーブの修正例を示すグラフである。 記憶部に予め記憶されている保磁力の温度係数の内容例を示すグラフである。 本実施例のＩＰＭモータの略示上方斜視図である。 本実施例の磁石の母材の磁力特性を示すグラフである。 本実施例の母材に対するＤｙ導入量と保磁力増加量ΔＨｃＪとの対応を示すグラフである。 保磁力増加量ΔＨｃＪ分布の算出精度の検証に用いた磁石の形状情報を示す説明図である。 本実施例の磁石（厚さ９．５ｍｍ）における深度δでの保磁力増加量ΔＨｃＪ分布の算出結果と実際の実測結果とを示すグラフである。 本実施例の磁石（厚さ２．５ｍｍ）における深度δでの保磁力増加量ΔＨｃＪ分布の算出結果と実際の実測結果とを示すグラフである。 本実施例の磁石について算出されたΔＨｃＪ分布の例を模式的に示す模式図である。 本実施例の磁石を用いたＩＰＭモータの減磁特性の算出結果及び実測結果を比較した例を示すグラフである。 本実施例の磁石（母材）が２％減磁したときの温度におけるＢｒ減少率の分布の算出結果を模式的に示す模式図である。 本実施例の磁石（Ｄｙ拡散材）が２％減磁したときの温度におけるＢｒ減少率の分布の算出結果を模式的に示す模式図である。

以下、本発明をその実施の形態を示す図面に基づき具体的に説明する。
なお、以下の実施の形態では、本発明に係るコンピュータプログラムに基づき、磁力特性算出方法をコンピュータに実行させて磁力特性算出装置として動作させ、重希土類元素としてジスプロシウム（以下Ｄｙと表記する）を拡散させたＮｄ−Ｆｅ−Ｂ系焼結磁石の磁力特性を算出する例を説明する。

図１は、本実施の形態における磁力特性算出装置１の構成を示すブロック図である。本実施の形態における磁力特性算出装置１には、パーソナルコンピュータが用いられる。磁力特性算出装置１は、各構成部による動作を制御し、演算を実行する演算部１０と、各種情報を記憶する記憶部１１と、演算部１０の処理に利用される一時記憶部１２と、可搬型記録媒体２から情報を読み出す読取部１３と、ディスプレイ１４、キーボード１５、マウス１６等の入出力装置及び演算部１０の間を中継するインタフェース（Ｉ／Ｆ）１７とを備える。

演算部１０は、ＣＰＵ（Central Processing Unit）、ＭＰＵ（Micro Processing Unit）などを用いる。演算部１０は、記憶部１１に記憶されている磁力特性算出プログラム１Ｐを読み出して実行する。これにより、演算部１０は設計に利用するＮｄ−Ｆｅ−Ｂ系焼結磁石の磁力特性を算出する各処理を実行する。

記憶部１１は、ハードディスク（Hard Disk）又はＳＳＤ（Solid State Drive）等の外部記憶装置を用いる。記憶部１１は、上述した磁力特性算出プログラム１Ｐが記憶されるほか、演算部１０が磁力特性を算出する処理の際に参照することができるように、後述するＤｙ導入量−ΔＨｃＪデータベース（導入量−保磁力増加量特性情報）１１１及び拡散条件データベース（Ｄｙなどの重希土類元素の拡散処理における拡散係数と拡散束と処理時間とを含む拡散条件の情報）１１２が記憶されている。

一時記憶部１２は、ＤＲＡＭ（Dynamic Random Access Memory）、ＳＲＡＭ（Static RAM）等の揮発性のランダムアクセスメモリを用いる。一時記憶部１２は、記憶部１１から読み出される磁力特性算出プログラム１Ｐなど、演算部１０の処理によって発生する各種情報が一時的に記憶される。

読取部１３は、ＤＶＤ、ＣＤ−ＲＯＭ、フレキシブルディスク等の可搬型記録媒体２からデータを読み出すことが可能である。可搬型記録媒体２には、コンピュータを磁力特性算出装置１として動作させるための磁力特性算出プログラム２Ｐが記録されている。記憶部１１に記憶されている磁力特性算出プログラム１Ｐは、演算部１０が読取部１３によって可搬型記録媒体２から読み出した磁力特性算出プログラム２Ｐを複製したものであってもよい。

Ｉ／Ｆ１７は、後述するように演算部１０によって出力される画像情報などをディスプレイ１４へ出力する処理、キーボード１５により入力される情報を検知して演算部１０へ通知する処理、マウス１６により入力される情報を検知して演算部１０へ通知する処理等を行なう。磁力特性算出装置１を操作し、磁石及び磁石を用いた製品を設計しようとするオペレータ（技術者）は、キーボード１５及びマウス１６を利用し、オペレータが設計に利用する磁石についての情報を入力し、当該磁石についての特性を演算部１０に算出させることが可能である。

上述のように構成される磁力特性算出装置１において、Ｄｙなどの重希土類元素を拡散させるＮｄ−Ｆｅ−Ｂ系焼結磁石の磁力特性として、熱及び減磁界による減磁率が算出され、所定の減磁率以下となる最高温度（減磁温度）が特定される過程について説明する。ここで、拡散される重希土類元素としてＤｙを例に挙げて説明する。まず、演算部１０は、特性の算出対象の磁石内におけるＤｙ拡散後の保磁力（ＨｃＪ）の増加量（ΔＨｃＪ）の分布を求める。そして演算部１０は、拡散前の母材の保磁力に基づいて磁石内の保磁力（ＨｃＪ）分布を求め、各部位におけるＪ−Ｈカーブを特定し、磁石全体の磁力特性として所定の温度（所定の第２温度、例えば１００℃）で使用された場合に常温（所定の第１温度、例えば２０℃）に戻った後の減磁率を算出する。更に演算部１０は、所定の減磁率以下となる減磁温度を特定する。

図２は、本実施の形態における磁力特性算出装置１の演算部１０が、Ｄｙ拡散後の磁石の磁力特性を算出する処理手順の一例を示すフローチャートである。

演算部１０は、磁石の大きさと形状とを示す形状情報を入力又は選択させる画面を作成してＩ／Ｆ１７を介してディスプレイ１４に出力し、キーボード１５及びマウス１６により磁石の形状情報をＩ／Ｆ１７を介して受け付ける（ステップＳ１）。演算部１０は、磁石のいずれの面からＤｙを拡散させ導入するかを示す導入面情報を入力又は選択させる画面を作成してＩ／Ｆ１７を介してディスプレイ１４に出力し、キーボード１５及びマウス１６により導入面情報をＩ／Ｆ１７を介して受け付ける（ステップＳ２）。ステップＳ１で受け付ける磁石の形状情報とは例えば、有限要素法のメッシュ情報（節点・要素情報）である。同様にステップＳ２で受け付ける導入面情報とは、形状情報に対応する導入面の数と導入面夫々を特定する情報である。

演算部１０は、受け付けた形状情報及び導入面情報に対応させて、記憶してある拡散条件（拡散係数、拡散束、処理時間）の情報を、拡散条件データベース１１２から読み出す。演算部１０は、受け付けた形状情報及び導入面情報について、読み出した拡散条件（拡散係数、拡散束、処理時間）の情報に基づき、フィックの拡散方程式（第２法則）である拡散方程式を用いて、磁石内におけるＤｙ導入量分布を算出する（ステップＳ３）。なお、拡散方程式、拡散条件（拡散係数、拡散束、処理時間）については後述する。

演算部１０は、ステップＳ３で算出した磁石内におけるＤｙ導入量分布を、Ｄｙ導入量−ΔＨｃＪデータベース１１１に基づき、磁石内の保磁力増加量ΔＨｃＪ分布を算出する（ステップＳ４）。

次に演算部１０は、拡散前の磁石即ち磁石の母材の磁石特性と、ステップＳ４で算出した保磁力増加量ΔＨｃＪとに基づき、第１温度（例えば常温（２０℃））でのＤｙ拡散後の保磁力ＨｃＪ分布を算出する（ステップＳ５）。演算部１０は、算出した保磁力ＨｃＪ分布に基づき、異なる温度係数を考慮し、減磁評価温度の第２温度（例えば１００℃）まで磁石の温度が上昇した場合のＪ−Ｈカーブを磁力特性として算出し（ステップＳ６）、磁石の温度が減磁評価温度の第２温度まで上昇した状態で減磁界が印加されるような負荷を与えて減磁が発生した後、第１温度に戻した場合のＪ−Ｈカーブを、記憶してある温度係数に基づいて算出する（ステップＳ７）。演算部１０は、ステップＳ６及びステップＳ７の算出結果に基づき、減磁評価温度の第２温度における減磁率を算出する（ステップＳ８）。なおここで減磁率は、減磁評価温度になる前後の常温における評価対象の磁石を用いたモータ特性のトルクの減少率を用いる。

次に演算部１０は、ステップＳ８で減磁率を算出した第２温度が所定の減磁率以下となる最大の第２温度であるか否かを判断する（ステップＳ９）。最大の第２温度でない場合には（Ｓ９：ＮＯ）、ステップＳ６へ処理を戻し、第２温度に別の温度を設定して、ステップＳ６−Ｓ８の処理を繰り返す。一方、演算部１０は、最大の第２温度であると判断した場合（Ｓ９：ＹＥＳ）、その第２温度を、磁石の減磁率が所定の率以下となる減磁温度として特定して、処理を終了する。

以下、図２のフローチャートに示した各処理について詳細を説明する。まず、ステップＳ３におけるＤｙ導入量分布の算出処理について説明する。

拡散方程式として、フィックの拡散方程式（第２法則）を用いる。このフィックの拡散方程式は、時間を無限と見なさない非定常状態における拡散処理、すなわち、拡散における濃度が時間の経過とともに変化する際に使用される。なお、境界条件は、ノイマン境界条件またはディリクレ境界条件を磁石表面に設定する。

拡散係数は、Ｄｙ拡散のし易さを表した係数であり、拡散束は、単位時間当たりに拡散面の単位面積を通過するＤｙ量であり、処理時間は、Ｄｙ拡散を施した時間である。本実施の形態では拡散係数を、Ｄｙ濃度依存性の係数とする。濃度依存性を考慮した拡散係数（関数）を決定することにより、このような拡散係数は同定される、なお、拡散係数を示す関数は、保磁力増加量ΔＨｃＪ実測値と一致するように決定される。図３は、拡散条件データベース１１２の内容例である濃度と拡散係数との関係を示す説明図であり、図３では、濃度Ｃに対する拡散係数Ｄをグラフによって示している。

以下に、拡散係数の同定処理（関数の決定処理）について説明する。図４は、演算部１０が拡散係数を同定する処理手順の一例を示すフローチャートである。

図３から理解されるように、濃度Ｃが高くなるにつれて拡散係数Ｄは指数関数的に減少する。まず、演算部１０は、濃度依存性を考慮した拡散係数Ｄの近似式として下記（１）式を定義する（ステップＳ３１）。
Ｄ＝ｋ１・ＥＸＰ（−ｋ２・Ｃ）＋ｋ３ …（１）
Ｃ：濃度
ｋ１，ｋ２，ｋ３：係数

演算部１０は、係数ｋ１，ｋ２，ｋ３の値を設定して、上記（１）式に従って拡散係数Ｄを計算する（ステップＳ３２）。演算部１０は、計算した拡散係数Ｄに基づき拡散方程式を用いて、Ｄｙ導入量分布を計算する（ステップＳ３３）。次に、演算部１０は、計算した磁石内におけるＤｙ導入量の分布を、Ｄｙ導入量−ΔＨｃＪデータベース１１１に基づき、磁石内の保磁力増加量ΔＨｃＪ分布に変換する（ステップＳ３４）。

演算部１０は、変換された保磁力増加量ΔＨｃＪ分布と実測された保磁力増加量ΔＨｃＪ分布とを比較し、その差が所定範囲内に入っているか否かを判断する（ステップＳ３５）。演算部１０は、所定範囲内に入っていないと判断した場合（Ｓ３５：ＮＯ）、ステップＳ３２へ処理を戻し、係数ｋ１，ｋ２，ｋ３に別の値を設定してあらためて拡散係数Ｄを計算した後、ステップＳ３３−Ｓ３５の処理を繰り返す。

一方、演算部１０は、前記差が所定範囲内に入っていると判断した場合（Ｓ３５：ＹＥＳ）、そのときの係数ｋ１，ｋ２，ｋ３の値を用いて拡散係数Ｄを同定して（ステップＳ３６）、処理を終了する。

次に、ステップＳ４におけるＤｙ導入量−ΔＨｃＪデータベース１１１について説明する。図５は、Ｄｙ導入量−ΔＨｃＪデータベース１１１の内容例であるＤｙ導入量と保磁力増加量との関係を示す説明図であり、図５では、Ｄｙ導入量に対する保磁力増加量ΔＨｃＪをグラフによって示している。Ｄｙ導入量−ΔＨｃＪデータベース１１１は具体的には、複数の異なるＤｙ導入量毎の保磁力増加量ΔＨｃＪの情報でもよいし、図５の説明図に示されている曲線を近似した数式であってもよい。

次に、ステップＳ５、Ｓ６及びＳ７におけるＪ−Ｈカーブの算出処理について詳細を説明する。Ｊ−Ｈカーブは、磁石の磁化曲線の内、磁化Ｊ（Ｔ）と磁場Ｈ（Ａ／ｍ）との関係を表した磁力特性を特定するための重要な情報である。

磁力特性算出装置１の記憶部１１には、拡散前の母材となる磁石の磁力特性の情報が記憶されてある。磁力特性の情報には、磁化曲線（Ｊ−Ｈカーブ、Ｂ−Ｈカーブ）が含まれる。演算部１０は、母材の磁化曲線を用い、拡散による保磁力増加量ΔＨｃＪにより、部位毎の保磁力ＨｃＪを求め、求めた保磁力ＨｃＪから、部位毎にＪ−Ｈカーブを算出する。算出されたＪ−Ｈカーブを、記憶してある保磁力毎の温度係数を用いて修正し、常温時（第１温度、例えば２０℃）の磁力特性、高温時（第２温度、例えば１００℃）の磁力特性を得る。

磁石が１００℃に加熱されて負荷が印加された後に、常温（第１温度、例えば２０℃）に戻った場合のＪ−Ｈカーブの修正方法について説明する。図６は、Ｊ−Ｈカーブの修正例を示すグラフである。横軸に磁場Ｈ、縦軸に磁化Ｊを示す。図６中の細い破線にて、２０℃及び１００℃での拡散後の磁石内の任意の部位におけるＪ−Ｈカーブを示している。

演算部１０の処理により、無負荷時のＪ−Ｈカーブ上のパーミアンス係数Ｐｃ´（図６中ラインｉ）が算出される。なお、ここでいうパーミアンス係数はＪ−Ｈカーブ上のパーミアンス係数であり、以下の説明でも同様である。無負荷時のパーミアンス係数Ｐｃ´は、磁力特性算出対象の磁石の形状及び磁気回路構造により決まり、無負荷時の２０℃での動作点は、ラインｉと２０℃のＪ−Ｈカーブとの交点（Ａ）になる。演算部１０は、２０℃で負荷が印加されたときの動作点Ｂを計算し、ラインｉを動作点Ｂと重なるように平行移動させてラインｉｉとし、これにより印加された減磁界Ｈｄを算出する。演算部１０は、１００℃のＪ−Ｈカーブとラインｉｉとにより、１００℃で減磁界Ｈｄが印加された場合の動作点Ｃを算出する。動作点Ｃは、１００℃のＪ−Ｈカーブの偏曲部（クニック）より下にあるため、不可逆減磁が発生した状態である。１００℃で減磁界Ｈｄが印加され減磁が発生したときの等価的なＪ−Ｈカーブは、図６中の太い破線にて示されている。

Ｊ−Ｈカーブにて磁化Ｊがゼロのときの磁場Ｈが保磁力ＨｃＪであり、磁場Ｈがゼロのときの磁化Ｊが残留磁束密度Ｂｒである。演算部１０は、１００℃で不可逆減磁した後の２０℃でのＪ−Ｈカーブを、１００℃で減磁界Ｈｄが印加され減磁したときの等価的なＪ−Ｈカーブ（図６中太い破線）のＨｃＪとＢｒとを、記憶してある上述の温度係数の情報に基づき算出する。１００℃で不可逆減磁が発生した後の２０℃でのＪ−Ｈカーブは、図６中の太い実線にて示されている。

なお上述したように、Ｄｙ拡散した磁石の保磁力は、磁石内で均一ではなく部位毎に保磁力が異なる。したがって残留磁束密度Ｂｒの低下の程度も部位毎に異なるから、磁石の減磁特性を精度よく求めるには、部位毎に保磁力の温度変化の違いを考慮しなければならない。したがって、異なる保磁力（絶対値）に応じた温度係数が必要である。

図７は、記憶部１１に予め記憶されている保磁力の温度係数の内容例を示すグラフである。図７では、横軸に保磁力ＨｃＪ（ｋＡ／ｍ）を示し、縦軸に温度変化に対する保磁力の変化の割合β（％／℃）を示している。図７中の白丸に示す温度変化係数の実測値から、実線に示すように保磁力ＨｃＪについて例えば２次の近似式を算出しておき、Ｄｙ拡散後の磁石内の部位毎に異なる任意の保磁力ＨｃＪに対して温度係数を使用できるようにする。このようにすることにより、異なる保磁力の値を部位毎に有する、即ち保磁力分布を有するＤｙ拡散後のＮｄ−Ｆｅ−Ｂ系焼結磁石の磁力特性を精度よく算出することが可能となる。

図２のフローチャートのステップＳ６及びＳ７にて、温度上昇による減磁特性を、Ｊ−Ｈカーブの温度係数による修正にて求めた。これによって、拡散後の磁石が１００℃に加熱された状態で負荷が印加され、減磁が発生した場合の減磁率をステップＳ８にて算出することができる。なおステップＳ８における算出処理では、演算部１０は、ステップＳ７にて部位毎に算出したＪ−ＨカーブをＢ−Ｈカーブに変換し、磁石全体の減磁率を算出するための既存のプログラムに基づく処理によって算出する。

なお、ステップＳ６では１００℃におけるＪ−Ｈカーブを算出し、１００℃にて減磁界が印加され減磁した後に２０℃に低下したときのＪ−ＨカーブをステップＳ７にて算出することにより、１００℃における減磁率を算出したが、ステップＳ６における温度を１００℃以外の温度にすることにより、異なる温度での減磁率をも算出することが可能であることは勿論である。ステップＳ８にて保磁力毎に異なる減磁率を複数算出し、ステップＳ９にて減磁率が所定の率以下となる温度（減磁温度）を特定してＤｙ拡散後の磁石の減磁特性として求めることができる。モータに用いるＮｄ−Ｆｅ−Ｂ系焼結磁石は、モータの回転、モータ周辺の環境等によって高温下での使用が多く、温度による残留磁束密度の低下度合いが重要である。つまり、どれほどの温度までであれば、継続的に使用しても減磁しないのかなどの情報が必要となる。したがって、本実施の形態における磁力特性算出装置１によって精度よく求まる減磁率または減磁温度は、非常に有用である。

上述したように算出される磁力特性の例（本発明の実施例）として、ＩＰＭモータに用いられる磁石の磁力特性、特に温度に対する減磁特性（減磁率）を算出し、実測した減磁特性と比較した。以下、この算出結果、比較結果について説明する。

図８は、本実施例のＩＰＭモータの略示上方斜視図である。図８中の３は、本実施例のＩＰＭモータであり、Ｍは、ＩＰＭモータ３に用いられるＤｙ拡散がされたＮｄ−Ｆｅ−Ｂ系焼結磁石である。ＩＰＭモータ３は、Ｖ字に配置されるようにロータに嵌めこんだ構成である。磁石Ｍ夫々は平板形状をなす。Ｄｙ導入面は、磁石Ｍ，Ｍ，…夫々のＩＰＭモータ３の外周面及び外周面に垂直な面である。

本実施例では、磁石Ｍの母材として以下を用いる。
ＮＭＸ−Ｓ５２（日立金属株式会社製、Ｎｄ−Ｆｅ−Ｂ系焼結磁石）

図９は、本実施例の磁石Ｍの母材の磁力特性を示すグラフである。図９は、２０℃、６０℃、１００℃及び１４０℃における磁化曲線であり、横軸に磁場Ｈ（ｋＡ／ｍ）、縦軸に磁化Ｂ又はＪ（Ｔ）を示している。上側の曲線がＪ−Ｈカーブであり、下側の曲線がＢ−Ｈカーブである。磁力特性算出装置１の記憶部１１に、図９に示す母材の磁化曲線を記憶しておくことにより、演算部１０は上述したようにＤｙ拡散後のＪ−Ｈカーブ及び温度変化後のＪ−Ｈカーブを算出することが可能である。

図１０は、本実施例の母材に対するＤｙ導入量と保磁力増加量ΔＨｃＪとの対応を示すグラフである。横軸にＤｙ導入量（mass％）、縦軸に保磁力増加量ΔＨｃＪ（ｋＡ／ｍ）をとって、白丸にて実測値を示し実線にて近似式を示している。なお、本実施例は、Ｄｙ拡散させた母材を２．８ｍｍ×２．８ｍｍ×１．０ｍｍの個片試料に切り出し、ＩＣＰ（Inductively Coupled Plasma）分析法を用いてＤｙ導入量を測定することにより得た。また、ΔＨｃＪは、ＶＳＭ（Vibrating Sample Magnetometer）で測定した試料の保磁力ＨｃＪからの母材の保磁力ＨｃＪの値の差により求めた。図１０に示すように対応するＤｙ導入量−ΔＨｃＪデータベース１１１を記憶部１１に記憶しておくことにより、演算部１０は上述したように磁石内におけるΔＨｃＪ分布を算出することが可能である。

本実施例では、磁石Ｍとして、特許文献１に記載の方法で、処理温度を９００℃とし、４．０時間Ｄｙを供給後、拡散させて製作した。

ここで、拡散係数を同定するために、実測した深度δにおける保磁力増加量ΔＨｃＪ分布への拡散方程式の算出の精度について検証した結果を示す。深度δは、拡散した表面からの距離を指す。図１１は、保磁力増加量ΔＨｃＪ分布の算出精度の検証に用いた磁石の形状情報を示す説明図である。検証には、縦が４２．５ｍｍ、横が３２．５ｍｍ、そして厚さが２．５ｍｍまたは９．５ｍｍである扁平直方体形状の厚みが異なる磁石を用い、夫々に同条件で特許文献１に記載の方法を用いてＤｙ拡散を行なった後に、磁石の上面中央部Ｃ及び周縁部Ｒにおける保磁力増加量ΔＨｃＪを測定した。なお、上面周縁部は、長辺（４２．５ｍｍ）からの距離２．０ｍｍであって短辺（３２．５ｍｍ）からの距離が約２４ｍｍの部分である。そして図１１の白矢印に示すように、上面及び４つの側面夫々から同様に磁石へＤｙを拡散させた。ここで、Ｄｙ濃度の測定にはＩＣＰ分析法を用いた。保磁力ＨｃＪの測定は、磁石を２．８ｍｍ×２．８ｍｍ×１．０ｍｍの測定部位毎に切り出し、ＶＳＭを用いて行なった。なお、ＩＣＰ分析法及びＶＳＭの測定ピッチをいずれも０．５ｍｍにて測定できるよう、同一条件で製作した複数の磁石で測定した。

磁石中央部の磁石の深度δにおける保磁力増加量ΔＨｃＪの実測値の分布に基づき、拡散方程式に使用する拡散係数を同定した。この際、拡散束を１．０×１０^-7（mass％）、処理時間を１４４００（sec）の固定条件とした。また、拡散係数を前述した式（１）のＤ＝ｋ１・ＥＸＰ（−ｋ２・Ｃ）＋ｋ３（Ｃ：濃度、ｋ１，ｋ２，ｋ３：係数）と近似して、図４のフローチャートに従って、これらの係数ｋ１，ｋ２，ｋ３の値を求めた。求めた係数を用いて、具体的に拡散係数Ｄを下記（２）のように同定した。
Ｄ＝５．０×１０^-11 ・ＥＸＰ（−７．０・Ｃ）＋１．１×１０^-11 …（２）

図１２は、厚さ９．５ｍｍの磁石の中央部Ｃ（図１２の（ａ））及び周縁部Ｒ（図１２の（ｂ））における深度（距離）δにおける保磁力増加量ΔＨｃＪ分布の算出結果と実際の実測結果とを示すグラフである。図１２では、横軸に深度δを「ｍｍ（ミリメートル）」単位で、縦軸に保磁力増加量ΔＨｃＪを「ｋＡ／ｍ（キロアンペア パー メートル）」単位で示している。図１２に示すように、深度（距離）δが０．５ｍｍ〜５．０ｍｍの範囲内で０．５ｍｍおきの算出結果（●印）と実測結果（○印）とは、中央部及び周縁部ともに精度良く合致している。

また、図１３は、厚さ２．５ｍｍの磁石の中央部Ｃ（図１３の（ａ））及び周縁部Ｒ（図１３の（ｂ））における深度（距離）δにおける保磁力増加量ΔＨｃＪ分布の算出結果と実際の実測結果とを示すグラフである。図１３の横軸、縦軸は、図１２と同じである。図１３に示すように、深度（距離）δが０．５ｍｍ〜２．０ｍｍの範囲内で０．５ｍｍおきの算出結果（●印）と実測結果（○印）とは、中央部及び周縁部ともに精度良く合致している。

図１４は、本実施例の磁石Ｍについて算出されたΔＨｃＪ分布の例を模式的に示す模式図である。図１４に示す例では、磁石Ｍの軸方向中央断面と幅方向中央断面とにおけるΔＨｃＪ分布を示している。上述したように磁石Ｍは、ＩＰＭモータ３の外周面及び外周面に垂直な面がＤｙ導入面である。したがって、磁石Ｍにおける保磁力増加量ΔＨｃＪは、磁石Ｍ夫々のＩＰＭモータ３の外周面側及び外周面に垂直な面側ほど大きく、ＩＰＭモータ３の中心側にて減少する分布を有すると算出されている。

図１５は、図１４に示したΔＨｃＪ分布を有する本実施例の磁石Ｍを用いたＩＰＭモータ３の減磁特性の算出結果及び実測結果を比較した例を示すグラフである。図１５の横軸は、減磁率を評価するための減磁評価温度（℃）を示し、縦軸は減磁率（％）を示す。○印はＤｙ拡散処理前の母材の異なる温度に対する減磁率の実測値、実線は該母材の減磁率の算出値であり、また、△印はＤｙが拡散された磁石Ｍの減磁率の実測値、破線は該磁石Ｍの減磁率の算出値である。

なお、減磁率は、磁石Ｍを用いたＩＰＭモータ３を温度が設定された恒温槽中で動作させた後に常温に戻し、常温状態でのトルクを測定し、負荷が印加される前の常温におけるトルクに対する減少率を算出して求めた。なお恒温槽の温度は、図１５に示したように、５０〜１４０℃まで１０℃おきに設定され、夫々の設定温度にてＩＰＭモータ３を回転させた。

図１５に示した減磁率の実測結果及び算出結果の内、減磁率が２％となるときの設定温度に着目すると、以下のようになる。母材の減磁率が２％となる温度は、実測値では１０１℃、算出値では１００℃であり、誤差はプラス１℃であった。Ｄｙが拡散された磁石Ｍの減磁率が２％となる温度は、実測値では１２２℃、算出値では１２４℃であり、誤差はプラス２℃であった。これにより、解析誤差が１０℃以下であり、十分な精度で減磁特性が解析できていると言える。Ｄｙ拡散によって磁石Ｍの減磁耐熱性が向上していることが分かる。

図１６及び図１７は、図１４に示したΔＨｃＪ分布を有する本実施例の磁石Ｍが２％減磁したときの温度におけるＢｒ減少率の分布の算出結果を模式的に示す模式図である。図１６に母材の例、図１７に拡散材（Ｄｙが拡散された磁石Ｍ）の例を示す。図１４と同様に、母材及び磁石Ｍの軸方向中央断面と幅方向中央断面とにおける分布を示している。上述したように、母材は１００℃で減磁率が２％となると算出されている。母材は保磁力の分布が均質であるが、Ｂｒ減少率は分布を有し、ＩＰＭモータ３の外側の角部分が最もＢｒ減少率が高い。Ｄｙが拡散された磁石Ｍは１２４℃で減磁率が２％となることが計算により得られており、このときのＢｒ減少率は、図１４に示した保磁力の増加量が大きい部位のＢｒ減少率が低く、保磁力の増加量が小さい部位のＢｒ減少率が高い傾向にある。なお、Ｂｒ減少率の算出は、磁石Ｍを減磁評価を行なう温度まで加熱し、減磁界を印加した後、室温（２０℃）に戻したときのＢｒの減少量から求めた。

本実施例に示したように、磁力特性算出装置１の磁力特性算出プログラム１に基づく処理によって、Ｄｙなどの重希土類元素を拡散して得られるＮｄ−Ｆｅ−Ｂ系焼結磁石の磁石内での部位毎に異なる保磁力増加量ΔＨｃＪの分布についての算出結果から、減磁率を精度よく求めることができる。

なお、重希土類元素としてＤｙを用いる例を説明した。しかしながら本発明はこれに限らないし、Ｔｂ等の重希土類元素を拡散して得られる磁石の磁力特性の算出にも広く適用され得る。また、磁石の形状は平板形状を用いる例を説明した。しかしながら、本発明は、これに限らず、弓形形状、リング形状、棒形状の磁石の磁力特性の算出にも広く適用され得る。また、モータは、ＩＰＭモータに限らず、ＳＰＭモータにも適用され得る。

なお、開示された実施の形態は、全ての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上述の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内での全ての変更が含まれることが意図される。

１ 磁力特性算出装置
１０ 演算部
１１ 記憶部
１１１ Ｄｙ導入量−ΔＨｃＪデータベース（導入量−保磁力増加量特性情報）
１１２ 拡散条件データベース（重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報）
１Ｐ 磁力特性算出プログラム
２Ｐ 磁力特性算出プログラム
Ｍ 磁石

Claims (10)

1. 重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を求める方法において、
   重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を予め記憶しておき、
   前記磁石の大きさと形状とを示す形状情報を受け付ける第１ステップ、
   受け付けた形状情報に対応させて、導入面情報を受け付ける第２ステップ、
   記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する第３ステップ、並びに、
   算出した導入量分布と、記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する第４ステップ
   を含むことを特徴とする磁力特性算出方法。
2. 前記拡散係数は、導入される重希土類元素の濃度依存性の関数で表されてあることを特徴とする請求項１に記載の磁力特性算出方法。
3. 重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを予め記憶しておき、
   記憶してある磁化曲線と、第４ステップで算出した保磁力増加量の分布とに基づき、前記磁石の各部位の所定の第１温度における磁化曲線を算出する第５ステップ、
   算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する第６ステップ、及び、
   第６ステップで算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する第７ステップ
   を更に含むことを特徴とする請求項１または２に記載の磁力特性算出方法。
4. 第４ステップで算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する第８ステップ、及び、
   前記磁石の減磁率が所定の率以下となる減磁温度を特定する第９ステップ
   を更に含むことを特徴とする請求項３に記載の磁力特性算出方法。
5. 重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を求める磁力特性算出装置において、
   重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を記憶しておく記憶手段、
   前記磁石の大きさと形状とを示す形状情報を受け付ける手段、
   受け付けた形状情報に対応させて、導入面情報を受け付ける手段、
   前記記憶手段に記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する手段、並びに、
   算出した導入量分布と、前記記憶手段に記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する保磁力増加量分布算出手段
   を備えることを特徴とする磁力特性算出装置。
6. 重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを予め記憶しておく手段、
   記憶してある磁化曲線と、前記保磁力増加量分布算出手段にて算出した保磁力増加量の分布とに基づき、前記磁石の所定の第１温度における磁化曲線を算出する手段、
   算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する手段、及び、
   算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する手段
   を更に備えることを特徴とする請求項５に記載の磁力特性算出装置。
7. 前記保磁力増加量分布算出手段にて算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する手段、及び、
   前記磁石の減磁率が所定の率以下となる減磁温度を特定する手段
   を更に備えることを特徴とする請求項６に記載の磁力特性算出装置。
8. 記憶手段を備えるコンピュータに、重希土類元素を表面から内部に拡散させ導入した磁石における磁力特性を、前記記憶手段に記憶させてある重希土類元素導入量に対する拡散させ導入したことによる保磁力の増加量の特性を示す導入量−保磁力増加量特性情報、及び、重希土類元素の拡散における拡散係数と拡散束と処理時間とを含む拡散条件の情報を用いて算出させるコンピュータプログラムであって、
   コンピュータに、
   前記磁石の大きさと形状とを示す形状情報を取得する第１ステップ、
   形状情報に対応させて、導入面情報を取得する第２ステップ、
   記憶してある前記拡散条件の情報に基づき拡散方程式を用いて、導入した重希土類元素の前記磁石内における導入量分布を算出する第３ステップ、並びに、
   算出した導入量分布と、記憶してある前記導入量−保磁力増加量特性情報とに基づき、前記磁石内における重希土類元素を拡散させ導入したことによる保磁力増加量の分布を算出する第４ステップ
   を実行させることを特徴とするコンピュータプログラム。
9. 重希土類元素拡散前の磁化曲線と、保磁力毎に異なる磁石の温度変化に対する保磁力変化率を示す温度係数の情報とを記憶してある記憶手段を更に用い、
   前記コンピュータに、
   記憶してある磁化曲線と、第４ステップで算出した保磁力増加量の分布とに基づき、前記磁石の各部位の所定の第１温度における磁化曲線を算出する第５ステップ、
   算出した磁化曲線と、記憶してある前記温度係数の情報とに基づき、所定の第２温度における磁化曲線を算出する第６ステップ、及び、
   第６ステップで算出した磁化曲線に基づき、第２温度で各部位に異なる減磁界が印加され減磁した後の前記所定の第１温度で減磁率を算出する第７ステップ
   を更に実行させることを特徴とする請求項８に記載のコンピュータプログラム。
10. 前記コンピュータに、
    第４ステップで算出した保磁力増加量の分布に基づき、前記磁石の異なる温度における減磁特性を算出する第８ステップ、及び、
    前記磁石の減磁率が所定の率以下となる減磁温度を特定する第９ステップ
    を更に実行させることを特徴とする請求項９に記載のコンピュータプログラム。

Images