(57)【特許請求の範囲】
【請求項１】
連続的に粉末層を与え、前記粉末層の選択された区域を融合することによる三次元体の製造方法であって、この区域は、前記三次元体の連続した横断面に対応し、
前記方法が前記層のうちの少なくとも１層のための以下のステップ、
作業領域上に前記少なくとも１層の粉末層を塗布するステップと、
前記少なくとも１層の粉末層の選択された区域を、放射銃から前記選択された区域にエネルギーを供給することにより融合するステップと、
を備える三次元体の製造方法であって、
前記少なくとも１層の粉末層の前記選択された区域を融合する際に使用される意図されるビーム経路を確立するステップと、
前記意図されるビーム経路に沿って移動すると仮定される仮想のビームの特定のエネルギー付与の関数として、前記意図されるビーム経路に沿った前記少なくとも１層の粉末層内の温度を計算するステップと、
前記計算された温度と、前記選択された区域を前記融合するステップのために設定された条件とに因って、前記意図されるビーム経路に沿った前記仮想のビームの前記特定のエネルギー付与を調整するステップと、
前記計算および前記調整に基づき、前記少なくとも１層の前記選択された区域を融合する際に前記意図されるビーム経路において使用される実際のビームの前記特定のエネルギー付与のための操作スキームを提供するステップと、
前記少なくとも１層の粉末層の前記選択された区域を融合する際に、前記特定のエネルギー付与のための操作スキームを用いるステップと、
を備えることを特徴とする、方法。
【請求項２】
前記特定のエネルギー付与は、単位時間および単位面積あたりに前記ビームにより付与された前記エネルギーを前記ビームの速さで除したものであり、また、前記特定のエネルギー付与は、ビームの速さ、ビーム出力またはビームスポットサイズのいずれか１つ、またはこれらの組み合わせを変化させることによって変化できることを特徴とする、請求項１に記載の方法。
【請求項３】
前記方法は、融合される材料に関連する既定のデータセットの使用を備え、前記データセットは、計算された温度および設定された条件の関数として選択される前記特定のエネルギー付与の値を備えることを特徴とする、請求項１または２に記載の方法。
【請求項４】
前記融合するステップのための前記設定された条件は、前記少なくとも１層の粉末層に対して次の条件、最大温度、作業温度、溶融深度および溶融幅のうちの１つまたは複数の条件を含むことを特徴とする、請求項１ないし３のいずれか一項に記載の方法。
【請求項５】
前記温度を計算するステップは、時間依存熱伝導方程式を解くステップを含むことを特徴とする、請求項１ないし４のいずれか一項に記載の方法。
【請求項６】
前記温度を計算するステップは、前記意図されるビーム経路に沿った局所的な温度を計算するステップを含むことを特徴とする、請求項１ないし５のいずれか一項に記載の方法。
【請求項７】
前記温度を計算するステップは、前記意図されるビーム経路に沿って配置された多数の点の中または近傍において実行されるいくつかの計算を含むことを特徴とする、請求項１ないし６のいずれか一項に記載の方法。
【請求項８】
隣接する計算点の間の最大距離は、前記隣接する点の間の前記特定のエネルギー付与の許容される変化のための制限値を設定することにより設定されることを特徴とする、請求項７に記載の方法。
【請求項９】
前記意図されるビーム経路を確立するステップは、
複数の取り得るビーム経路に沿って前記温度の計算を行うステップと、
前記複数の取り得るビーム経路の中から前記意図されるビーム経路を選択するステップと、
を含むことを特徴とする、請求項１ないし８のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、連続的に粉末層を与え、この粉末層の選択された区域を融合することによる三次元体の製造方法に関する。この区域は、三次元体の連続した横断面に対応する。
【背景技術】
【０００２】
高エネルギーの電磁放射線または電子線を照射することによって融合され、かつ、凝固される粉末状の材料を用いて、層の上に層を重ねながら三次元の物体を製造する機器は、例えば、米国特許第４，８６３，５３８号、米国特許第５，６４７，９３１号およびスウェーデン特許第５２４，４６７号から知られている。このような装置は、例えば、粉末の供給源と、上下に調整可能な台すなわち作業領域の上に粉末層を連続的に塗布する手段と、作業領域全体にビームを向ける手段とを含む。粉末は、ビームが一層ずつ作業領域上を移動するとき焼結または溶融して凝固する。
【０００３】
高エネルギービームを用いて粉末を溶融または焼結する場合、被照射材料の温度を徹底的に制御して、対象物に適切な材料特性を与えること、また、幾何学的な変形を防止することが重要である。例えば、局所的な温度が高すぎれば、生成中の対象物を破壊することもあるし、温度分布があまりに不均質であれば、クラックを引き起こすこともある。さらに、完全に融合させるためには、溶融する工程の間における粉末床の上層の温度は、通常、最小値よりも高く保たなければならない。温度を制御し続ける他、製造時間を減少するように努めること、すなわち、選択された区域にわたってビームを可能な限り効率的に掃過するように努めることが、通常は重要である。
【０００４】
各粉末層の選択された部分または区域のみが、融合される。ビームは、選択された区域にわたって、この区域を完全に融合させる走査パターンまたはハッチングパターンで、特定の経路内を掃過する。多くの場合、この走査パターンは、選択された区域にわたって等間隔で配置された平行線という形態を有する。これらの選択された区域のそれぞれは、いくつかの部分区域を含んでいてもよいが、粉末床内において構築中の対象物の横断面に対応する。
【０００５】
ビームを平行線状の走査パターンで掃過するには、線を順番に走査すればよい。先に走査された線に沿った加熱された材料からの熱伝導のために、これから走査される特定の線に沿った材料内の温度は、開始温度よりも高くなる（すなわち、最初の線が走査されたときの材料内の温度よりも高くなる）。少なくとも高エネルギービームを用いる際には、材料内の局所的な温度を適切に維持するために、この温度の蓄積を考慮に入れなければならない。
【０００６】
これを考慮に入れる一つの方法は、温度の蓄積に応じてビームのエネルギー入力を調整することである。このことは、例えば、ビーム出力を変化させることにより、または、粉末層にわたってビームが移動する速さを変化させることによりなされてもよい。１つの例は、第１の走査線の終端が第２の走査線の始端に近接するビームの方向転換位置において、ビームの速さを増加することである。しかしながら、これを適切に行うには、材料内の温度の情報を入手する必要がある。この温度、さらに正確には、粉末床の表面温度は、サーマルカメラを用いて計測できる。しかしながら、ビームのリアルタイム補正または制御をそのようなカメラからの入力に基づいて行うことは、システムの応答時間が長いために、適切に行うことが難しい（温度上昇が検出されたときに温度を下げるように直ちに作用し始めたとしても、温度はしばらくの間増加し続ける可能性が高い）。サーマルカメラは、製造後、製造プロセス中になにか不具合が発生していないかの検査に役立つこともあろう。
【０００７】
米国特許第５，９０４，８９０号は、ビームの走査速さを平行線状の走査パターン内の走査線の長さの関数として変化させる方法を開示している。ビームの速さは、より長い走査線に対してはより低く、また、より短い線に対してはより高くして、特定の区域からビームが離れているときに、冷却のされ方が多様となるのを防止する。製造された製品内の密度分布を均質にすることが目的である。この方法は、ビームの速さが走査線の長さの割に高い場合には、上述の温度の蓄積に関しては役立つかもしれない。しかしながら，走査線が長い場合、ビームの速さは、走査線の終端部分のみで調整されるべきであり、また、線が同じ粉末層のいくつかの選択された区域にわたって配置される場合、または、異なるパターンの場合、温度の蓄積は、区域のあらゆる部分で同様にはならない。さらに、ビームのエネルギーが高い場合、より複雑な走査パターンが求められることもある。そのような場合には、温度の蓄積は、ビームの速さを走査線の長さに対して単に変化させることでは適切に考慮されない。
【０００８】
国際公開特許公報第２００８／０１３４８３号は、平行な走査線が特定の順番で走査され、連続して走査される線どうしの間に最小安全距離を確立できるようにする方法を開示している。このように、連続して走査される線どうしの間の熱伝導の干渉の発生を妨げることにより、走査線の間での温度（および荷電粒子）の蓄積が考慮されている。この方法は、主として高いビームの速さおよび高いビーム出力での粉末層の予熱を意図しているが、粉末を溶融するステップの間、熱伝導の干渉を防止するために用いることもできる。しかしながら、これはいくぶん時間がかかり過ぎる製造プロセスとなる。
【０００９】
したがって、徹底的な温度制御の他に時間効率の良い製造もまた可能にする、より精巧な走査方法が求められている。
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本発明の目的は、三次元体の製造のために上に説明したタイプの方法を提供することにあり、この方法は、温度の制御および製造の高速化の発展性が高められたことを呈する。この目的は、独立請求項に含まれた技術的特徴により規定された方法により達成される。従属請求項は、本発明の有利な実施形態、さらなる発展例および変形例を含む。
【課題を解決するための手段】
【００１１】
本発明は、連続的に粉末層を与え、この粉末層の選択された区域を融合することによる三次元体の製造方法に関し、この区域は、三次元体の連続した横断面に対応する。この方法は、前記層のうちの少なくとも１層のために以下のステップ、作業領域上に少なくとも１層の粉末層を塗布するステップと、放射銃から選択された区域にエネルギーを供給することにより少なくとも１層の粉末層の選択された区域を融合するステップとを備える。
【００１２】
本発明は、本方法が、少なくとも１層の粉末層の選択された区域を融合する際に使用される意図されるビーム経路を確立するステップと、意図されるビーム経路に沿って移動すると仮定される仮想のビームの特定のエネルギー付与の関数として、意図されるビーム経路に沿った少なくとも１層の粉末層内の温度を計算するステップと、計算された温度と、選択された区域を融合するステップのために設定された条件とに因って、意図されるビーム経路に沿った仮想のビームの特定のエネルギー付与を調整するステップと、計算および調整に基づき、少なくとも１層の選択された区域を融合する際に意図されるビーム経路において使用される実際のビームの特定のエネルギー付与のための操作スキームを提供するステップとを備えることを特徴とする。
【００１３】
「意図されるビーム経路」の用語は、選択された区域にわたって配置される走査パターンまたは線パターンに関し、また、選択された区域にわたって粉末を溶融／融合する目的でその区域をビームが掃過される際に、ビームスポットが辿るように意図される経路の少なくとも一部分のことを指す。原理上は、選択された区域内において粉末を完全に融合させる限りにおいて、意図されるビーム経路はあらゆる形態を有することができ、換言すると、例えば、分割されていても連続であってもよく、直線部分および曲線部分の両方を含むこともできる。さらに、ビーム経路は、例えば、線が異なる順番で走査される場合、または、単一の線が反対方向に走査される場合など、線パターンが同じであっても変化し得るものである。
【００１４】
「意図されるビーム経路に沿って移動すると仮定される仮想のビームの特定のエネルギー付与の関数として、意図されるビーム経路に沿った少なくとも１層の粉末層内の温度を計算する」ステップとは、意図されるビーム経路の中または近傍において、その延長部に沿って、局所的な温度または局所的な温度分布が計算されることを意味している。これは例えば、意図されるビーム経路に沿って移動する間に特定のエネルギー付与を発生すると仮定される仮想のビームにより材料に付与されたエネルギーを考慮に入れて、意図されるビーム経路に沿って配置された多数の点において局所的な温度（分布）を計算することなどによりなされる。
【００１５】
意図されるビーム経路沿った特定の点における（すなわち、特定の時点における）局所的な粉末層の温度は、例えば、開始時の材料層の温度分布、材料の熱的特性（熱伝導率など）、仮想のビームの特定のエネルギー付与の履歴（ビームの現在位置、および現在位置に至るまでの間のその経路において材料層に付与されたエネルギーまたは出力の量を含む）、およびビーム経路の幾何パターンに因る。
【００１６】
「ビームの特定のエネルギー付与」の用語は、単位時間および層の単位面積あたりに（仮想または実際の）ビームにより付与されたエネルギー（ビーム出力およびスポットサイズ）、すなわち、単位面積あたりに付与された出力をビームの速さで除したものを指す。したがって、特定のエネルギー付与を変化させることは、層の表面上を移動するビームの速さを変化させる、またはビームの出力を変化させる、またはビームのスポットサイズ（すなわち、特定の時点においてビームに直接曝される層の表面区域）を変化させる、またはこれらの組み合わせによってなされる。計算上、仮想のビームの特定のエネルギー付与の履歴もまたこのように、速さ、出力またはスポットサイズの変化を含む。また、ビームの形状およびビーム内におけるエネルギー／出力分布も変化してもよく、計算に含まれていてもよい。
【００１７】
計算は複雑で時間がかかり過ぎることもあるため、今までどおり特定のエネルギー付与の履歴（これは、ビームは未到達であるが、先の意図されるビーム経路の既に「融合された」部分からの熱伝導がある意図されるビーム経路の点における温度に強い影響を与えることもある）を考慮にいれながら、十分に正確な温度を計算できるような簡単化を様々に行うことができる。
【００１８】
「計算された温度と、選択された区域を融合するステップのために設定された条件とに因って、意図されるビーム経路に沿った仮想のビームの特定のエネルギー付与を調整する」ステップとは、例えば、最大温度のために設定された条件よりも、特定の点において温度が高くなったことを計算が示す場合などに、意図されるビーム経路の特定の部分にわたって、少なくとも１つのビームのパラメータ、すなわち、ビームの速さ、出力およびスポットサイズ、またはいずれか１つが調整されることを意味している（これは、例えば、特定の点の近傍のビームの速さの増加もしくはビーム出力の低減を要求する、または、先のビーム経路の部分からの間接的な、熱伝導によるその点の加熱を低減するように、特定のエネルギー付与の履歴の変更を要求することになろう）。
【００１９】
意図されるビーム経路に沿った仮想のビームの特定のエネルギー付与の調整は、経路（の部分）に沿った温度の再計算がその他のビームのパラメータを用いて行われるように扱われてもよい。別の方法として、あるいは補助として、融合される材料に関連する既定のデータセットを利用でき、前記データセットは、計算された温度および設定された条件の関数として特定のエネルギー付与の適切な値を備える。このような規定のデータは、時間がかかり過ぎる再計算を防ぐのに役立ち、例えば、意図されるビーム経路に沿って配置された多数の点において温度が計算される際に使用できる。仮想のビームの現在位置に対応する点の前方に比較的近接して位置する「次の」点において計算された温度に応じて、現在位置から「次の」点に到達するまでビームを移動する際に使用する特定のエネルギー付与の適切な値は、既定のデータから直接取得できる。意図されるビーム経路沿って配置された残りの点に対して、この手順は繰り返される。このようにして、特定のエネルギー付与は、意図されるビーム経路沿って段階的に調整される。
【００２０】
（特定のエネルギー付与のための）「操作スキーム」の用語は、実際のビームの特定のエネルギー付与、すなわち、速さ、出力およびスポットサイズのそれぞれを、粉末を融合するステップの間、時間（または、時間に関係するビーム経路に沿った位置）に対してどのように変化させるべきか、その方法を指す。したがって、操作スキームは、選択された区域を融合する際にビームの速さ、出力およびスポットサイズがどのように変化すべきかについての情報を含む。この操作スキームを提供または決定／確立するステップは、先のステップからの結果の減算および加算の形である。上記の特定のエネルギー付与の段階的な調整の例において、操作スキームは、ビームのパラメータの段階的な変動を含む。操作スキームは、意図されるビーム経路の初期部分のためなど、温度計算および特定のエネルギー付与の調整が必要ないと考えられる意図されるビーム経路の部分のためのビームのパラメータの設定の情報を含むこともできる。
【００２１】
材料内の温度は、そのエネルギー含量に関係する。したがって、真の温度を計算するかわりに、別のエネルギー関連パラメータおよび温度関連パラメータを計算し利用することができる。計算された温度の用語は、そのような関連のパラメータも含む。
【００２２】
意図されるビーム経路を確立するステップ、意図されるビーム経路に沿って温度を計算するステップ、仮想の特定のエネルギー付与を調整するステップ、および、操作スキームを決定するステップは、一度に実行されなくてもよいし、厳密に与えられた順番で実行されなくてもよい。例えば、計算および調整は、反復するやり方で実行でき、また、操作スキームは、ビーム経路全体の断片に対して段階的に決定できる。さらに、意図されるビーム経路を確立するステップは、所与の走査方向の等間隔に配置されたまっすぐな平行線からなる既定の線パターンが選択可能であり、いくぶん単純であるが、このステップは、好ましい線パターンおよび好ましい、最終的に選択された、意図されるビーム経路を獲得するための計算および調整を備えていてもよい。
【００２３】
結果、つまり、本発明は、粉末を融合する際に使用されるビームの特定のエネルギー付与が予め調整されて、異なる特定のエネルギー付与および条件におけるビーム経路に沿って結果として生じる温度を計算することにより、使用される特定の走査パターンの温度の蓄積に応じて変化するようにする方法について言及する。換言すると、本発明の方法は、計算および調整により、経路パターンに沿って通過して粉末を溶融する際に、ビームの特定のエネルギー付与が時間（または選択された区域上の位置）に対してどのように変化すべきかを予め決定できるようにする。
【００２４】
特定のエネルギー付与の操作スキームを最適化するために、例えば、製造時間を最小化する、特定の最大温度を超えないようにする、特定の時間間隔の間、特定の温度を超えないようにする、取得される温度の最大値を最小化する、ビーム経路に沿った溶融した材料の幅を均等にする、および、例えば、製造時間の最小化と取得される温度の最大値との折衷案などのこれらの様々な組み合わせなど、様々な条件を計算において使用してもよい。意図されるビーム経路を選択する前に、取り得る様々なビーム経路を評価できる。
【００２５】
計算を簡単化し速度を上げるために、条件は、ビームのパラメータ（速さ、出力およびスポットサイズ）のうちの１つまたは２つの既定の（予め計算された）値、および、例えば、互いに等間隔に配置された平行線のセットなどの既定のビーム経路の両方、またはいずれか一方を含むことができる。
【００２６】
本発明の方法は、包括的であり、選択された区域のあらゆる幾何形状に適用可能である。粉末層は、同様または異なる幾何形状を有する複数の選択された区域を備えていてもよいことに注意されたい。
【００２７】
適切な特定のエネルギー付与の操作スキームが決定されると、このスキームは、問題の層の選択された区域（の部分）の実際の溶融／融合に使用される。本発明の方法は、好ましくは、対象物が形成される層の全て、または、少なくともほとんどにおいて使用される。
【００２８】
本発明の効果は、徹底的な温度制御および選択された区域の温度分布制御を提供して、洗練された方法で融合ステップを計画可能とすることにある。同様に、これは、温度が高くなり過ぎる（これは、構築中の製品を破壊することもある）のを防止するため、温度分布を均質にする（これは、歪みおよびクラック形成を減らして製品の特性を改善する）ため、および、製造を速める（これは、製造の費用効率を高める）ために使用できる。
【００２９】
本発明の有利な実施形態において、本方法は、少なくとも１層の粉末層の選択された区域を融合する際に、特定のエネルギー付与のための操作スキームを用いるステップを備える。
【００３０】
本発明の別の有利な実施形態において、特定のエネルギー付与は、単位時間および単位面積あたりにビームにより付与されたエネルギーをビームの速さで除したものであり、また、特定のエネルギー付与は、ビームの速さ、ビーム出力またはビームスポットサイズのいずれか１つ、またはこれらの組み合わせを変化させることによって変化できる。
【００３１】
本発明の別の有利な実施形態において、本方法は、融合される材料に関連する既定のデータセットの使用を備え、前記データセットは、計算された温度および設定された条件の関数として選択される特定のエネルギー付与の値を備える。
【００３２】
本発明の別の有利な実施形態において、融合するステップのために設定された条件は、少なくとも１層の粉末層に対して次の条件、最大温度、作業温度、溶融深度および溶融幅のうちの１つまたはいくつかの条件を含む。
【００３３】
本発明の別の有利な実施形態において、温度を計算するステップは、時間依存熱伝導方程式を解くステップを含む。
【００３４】
本発明の別の有利な実施形態において、温度を計算するステップは、意図されるビーム経路に沿った局所的な温度を計算するステップを含む。
【００３５】
本発明の別の有利な実施形態において、温度を計算するステップは、意図されるビーム経路に沿って配置された多数の点の中または近傍において実行されるいくつかの計算を含む。
【００３６】
この実施形態の変形において、隣接する計算点の間の最大距離は、隣接する点の間の特定のエネルギー付与の許容される変化のための制限値を設定することにより設定される。例えば、ビームの速さだけが変化される場合、ビームの速さの許容される変化の最大値が設定される。
【００３７】
本発明の別の有利な実施形態において、意図されるビーム経路を確立するステップは、複数の取り得るビーム経路に沿って温度の計算を行うステップと、前記複数のビーム経路の中から意図されるビーム経路を選択するステップとを含む。
【００３８】
後述の本発明の詳細な説明において、次の図面を参照する。
【図面の簡単な説明】
【００３９】
【図１】本発明の方法を適用可能である、公知の三次元製品製造装置の例を示す概略図である。
【図２】ビームがｘ軸の正の方向に進行中の場合の表面温度プロファイル、ならびに、矩形の囲み線内にて対応する溶融深度および溶融幅を示す概略図である。
【図３】ＦＥＭにより計算された、いくつかの温度分布プロファイルを、式３におけるガウス関数の級数により近似された分布とともに示すグラフである。
【図４】ＦＥＭにより計算された、いくつかの温度分布プロファイルを、式３におけるガウス関数の級数により近似された分布とともに示すグラフである。
【図５】ＦＥＭにより計算された、いくつかの温度分布プロファイルを、式３におけるガウス関数の級数により近似された分布とともに示すグラフである。
【図６】点と線との距離および点と点との距離、それぞれ
【数１】
0005712306.tif 000002
および
【数２】
0005712306.tif 000003
を示す図である。ここで、
【数３】
0005712306.tif 000004
は、指数項のためのグローバル座標系における位置であり、
【数４】
0005712306.tif 000005
および
【数５】
0005712306.tif 000006
は、線ｊの線分ｋのためのグローバル座標系における座標である。
【図７】等脚台形の形状を有する選択された区域のための意図されるビーム経路の例を示す図である。ここで、意図されるビーム経路とは、ビームが、左から右へ向かうのに続いて右から左へと方向を変えながら、底部から頂部へ向かって線を走査し始めるというようなものである。
【図８】図７に示す意図されるビーム経路において使用されるビームの特定のエネルギー付与に対して決定された操作スキームを示す図である。ここで、この例における特定のエネルギー付与は、ビームの速さを変化させることにより変化される。
【発明を実施するための形態】
【００４０】
図１は、三次元製品を製造するための公知装置１の例を示す。装置１は、三次元製品３をその上で組み立てる上下に調整可能な作業テーブル２と、１つまたは２つ以上の粉末分配器４と、作業テーブル２上に粉末の薄い層を連続的に分布させて粉末床５を形成するように配置された手段２８と、粉末床５にエネルギーを送給して、粉末床５の部分を融合させる電子銃の形態の放射銃６と、前記作業テーブル２にわたって、放射銃６により照射される電子ビームを動かし、形状を決定するための偏向コイルおよびビーム形状を決定するコイル７と、装置１の様々な部品を制御するように配置された制御ユニット８とを備える。
【００４１】
一般的な作業サイクルでは、作業テーブル２が下げられ、新しい粉末層が粉末床５の上の作業領域上に塗布されて、粉末床５の上部層５´の選択された部分にわたって電子ビームが走査される。原則として、このサイクルは製品が完成されるまで繰り返される。当業者であれば、図１に概説されたタイプ、および、電子銃のかわりにレーザ銃を備える装置の両方に関して、三次元製品の製造装置の一般的な機能および構成に精通する。
【００４２】
従来、電子銃を有する装置は、通常、少なくとも１０−２ｍｂａｒ未満の真空で作動して、電子ビームが、電子銃と作業領域との間に存在する電子または分子と相互作用するのを防ぐ。
【００４３】
等脚台形の形状を有する、粉末層の選択された区域の例を図７に示す。意図されるビーム経路も示す。
【００４４】
ここで本発明の方法の実施形態を説明する。この実施形態の例において、意図されるビーム経路は、互いに等間隔に配置された複数の平行かつまっすぐな線（走査またはハッチングの線）を辿る。この例において調整されるビームのパラメータは、ビームの速さである。計算において、ビームの速さは、特定の深さにおける溶融された材料の幅（図２の溶融幅および溶融深度を参照）が全ビーム経路沿いで等しくなるように調整される。このことにより、ビーム経路の平行部分の間の距離を一定にして用いることができる。残るパラメータは、予め決定されている（または、その他の既定のパラメータから計算される）。
【００４５】
概要として、本方法の実施形態は以下のように説明できる。
１．異なる組み合わせの材料特性、材料の温度およびビーム出力のための温度プロファイルおよび関連するビームのパラメータ（スポットサイズおよびビームの速さ）から構成されるデータは、データベースに生成されて格納される。これらのデータは、図２の試験用の直方体と同様の単純な幾何形状に基づくＦＥＭ計算により得ることができる。
２．三次元体を製造するために使用される機械装置は、時間依存熱伝導方程式を解くことにより、ビームの軌道（経路）に沿って配置された多数の点の局所的な温度分布をリアルタイムで計算する。方程式の解は、先に融合された（すなわち仮想的に融合された）ハッチング線の温度プロファイルをガウス包絡線を用いて展開することにより得られる。使用されたビームおよび材料のパラメータに対応する温度プロファイルは、データベースから取得される。
３．特定の点におけるビームのパラメータは、計算された局所的な温度分布に応じて選択され、これらは、データベース内の事前に計算されたデータから取得される（これは、計算された温度分布を、事前に計算された使用される材料の温度プロファイルと比較し、計算された分布に最適にあてはまるプロファイルに対応するビームのパラメータを選択することによりなされる）。
４．１本のハッチング線が完了した時点で、線の終端における温度プロファイルもまたガウス関数により近似され、ステップ２および３（すなわち先の２つのステップ）が次のハッチング線に対して繰り返される。
【００４６】
リアルタイムで計算を行うという表現は、粉末の融合が計算と同時に行われることを意味する。一般的に、後続の層のビームのパラメータの操作スキームの計算は、先の層が融合される間に行う。原則として、第１の層の融合プロセスを開始する前に、全ての層に対する操作スキームの計算および決定を全て行うことができるが、そうすると通常、製造を開始するまでに待ち時間を生ずることになる。逆に極端な例において、実際のビームが配置される場所に非常に近接するビーム経路に沿った点に対して、操作スキームの計算および決定を行うと、計算または融合において何か不具合が起こった場合に補正や再計算をなす余裕がごくわずかになる。
【００４７】
本方法の実施形態への導入
説明した方法にしたがって、溶融プロセスを制御するために要する適切なデータを得るために、熱源が無く、その領域が−∞＜ｘ＜∞，−∞＜ｙ＜∞および−∞＜ｚ＜０である一様な材料における時間依存熱伝導方程式を考える。
式１ａ．
【数６】
0005712306.tif 000007
ここで、Ｔ（ｘ，ｙ，ｚ，ｔ）は時間依存性の温度分布、λは熱伝導率、ｃｐは熱容量およびρは材料の密度である。
【００４８】
境界条件は、次のように要約される。
式１ｂ
【数７】
0005712306.tif 000008
仮想のエネルギービームを記述するために、ｚ＝０上をｘ方向に移動するガウス分布型のソース項を使用する。同じ頂面を介する放射は、シュテファン・ボルツマンの法則に従うと仮定される。
式１ｃ
【数８】
0005712306.tif 000009
ここで、Ｐｉｎは吸収されたビーム出力、ｖｘはビームの速さ、σは分散（ビームスポットサイズ）、ｒａｄｃｏｅｆｆは表面からの放射係数、および、Ｔｓｕｒは表面上方の周辺温度である。
【００４９】
Ｔｏは、作業温度、すなわち、溶融／融合する前における材料の所望の温度である。
【００５０】
データ生成に要する時間を減らすため、移動するスポットの周囲の温度分布が定常状態（ｘ＝ｘ−ｔｖｘ，ｄｔ＝−ｄｘ／ｖｘ）に達したと仮定することにより、時間依存性を除去するのが妥当であることがある。
式２ａ
【数９】
0005712306.tif 000010
式２ｂ
【数１０】
0005712306.tif 000011
式２ｃ
【数１１】
0005712306.tif 000012
上述の熱伝導方程式は、例えば、材料特性、Ｔｏおよびビーム設定のいくつかの異なる組み合わせに対してＦＥＭ技術により解くこともできる。この手順のしくみは、図２に例示されている。
【００５１】
図２は、「試験用の直方体」を表しており、ビームがｘ軸の正方向に移動している。表面における温度プロファイルは、溶融体積が材料の溶融温度に対応する等温曲線により表される断面とともに示されている。ここで、ビームのパラメータｖｘおよびσは、特定の溶融体積プロファイルが得られるように、溶融深度および溶融幅に関して最適化されている。さらに、材料内の最大温度は、Ｔｍａｘに制限されている。もちろん、他にもビームのパラメータを最適化するために用いる条件はあり得る。例えば、溶融体積内の温度勾配を最小化することは、そのような条件の１つであり得る。
【００５２】
ハッチング線の終端におけるエネルギー入力の記述に要する温度プロファイルは、式２ａのＴ（ｘ，ｙ，ｚ）をガウス関数の級数で近似することにより得ることができる。そうすることで、任意の数のハッチング線についても半無限領域における温度分布の解析解が後に取得可能となる。級数Ｔ´（ｘ，ｙ，ｚ）は次のとおりである。
式３
【数１２】
0005712306.tif 000013
パラメータＡｉ，ｘｐｏｓｉ，σｘｉ，σｙｉ，σｚｉおよびａは、点別の非線形最小二乗法によるＴ（ｘ，ｙ，ｚ）とＴ´（ｘ，ｙ，ｚ）との間の曲線あてはめから取得可能である。ここで、ｘｐｏｓｉは、指数項ｉのビーム経路に沿ったｘの位置である。ビームの座標系において、ビームはｘ軸の正方向へ移動すると仮定され、ｘ＝０に位置されるため、これは負の値になる。
【００５３】
図３−図５において、ＦＥＭにより計算されたいくつかの温度分布が、式３により近似された分布とともに示されている。
【００５４】
適合度は、使用されるガウス関数の数により主に決定される。以下の例において、Ｎの値は１０から１２であり、各温度プロファイルに対して３０から３６のガウス関数が使用されたことを意味する。
【００５５】
材料内の時間依存性の温度分布
ビームが１本の線を走査した後の材料内の時間依存性の温度分布Ｔ´（ｘ，ｙ，ｚ，ｔ）は、グリーン関数と、式３から得られる初期条件Ｔ´（ｘ´，ｙ´，ｚ´）との畳み込みにより得られる。
式４
【数１３】
0005712306.tif 000014
上式で、
【数１４】
0005712306.tif 000015
ここで、材料の温度は、Ｔｓｕｒｆに等しく、Ｔｏとは異なると仮定する。ここでは表面を介する熱損失はゼロとする。
【数１５】
0005712306.tif 000016
ビームがＭ本の線を走査した際、式４の右辺は和で置き換えられる。
式５
【数１６】
0005712306.tif 000017
上式で、ｔｊは線ｊが完了した際の時間に等しく、Ｔｏｊは線ｊが完了した際のスポット周辺の温度、Ｔｊ´（ｘ´，ｙ´，ｚ´）は線ｊの式３による温度分布、および、Ｈ（ｔ−ｔｊ）はヘヴィサイドの階段関数であり、次のように定義される。
【数１７】
0005712306.tif 000018
【００５６】
Ｔ´（ｘ，ｙ，ｚ）の式（式４）を式５へ代入する際には、式５のｘ´およびｙ´座標がその部分の表面により決定されるグローバル座標系に関することに対して、式３における座標ｘ，ｘｐｏｓｉおよびｙは、線ｊの端点を軸とし、ｘ軸がこの線のビームの移動方向を指すローカル座標系に関することを考慮すべきである。さらに、線ｊのビーム経路が、それぞれ異なる方向のいくつかの線分による記述を要する場合、式３のｙは、線ｊの線分ｋと点（ｘ´，ｙ´）との間の距離である
【数１８】
0005712306.tif 000019
により置き換えなければならず、また、ｘ´−ｘｐｏｓｉは、線分上の指数項ｉのｘの位置と同じ線分の（ｘ´，ｙ´）の射影点との間の距離である
【数１９】
0005712306.tif 000020
により置き換えなければならない（図６参照）。
【００５７】
この方法により、あらゆる種類のビーム経路を考慮できる。しかしながら、式３における温度分布は直線におけるシミュレーションから得られることを忘れてはならない。したがって、ビーム経路の曲率が非常に著しい場合、ｘｐｏｓｉ値により決定されたものと同じ距離でこの経路に沿って式３の項を単に配置するだけでは、うまく近似できないであろう。このような場合には、曲線形状上のＦＥＭ解が必要となることもある。
【００５８】
図６は、点と線との距離および点と点との距離、それぞれ
【数２０】
0005712306.tif 000021
および
【数２１】
0005712306.tif 000022
を示す。
【数２２】
0005712306.tif 000023
は、指数項のためのグローバル座標系における位置である。
【数２３】
0005712306.tif 000024
および
【数２４】
0005712306.tif 000025
は、線ｊの線分ｋのためのグローバル座標系における座標である。
【００５９】
【数２５】
0005712306.tif 000026
に配置される指数項を少なくとも１つ含む線分ｋｊのそれぞれについて、距離の二乗である
【数２６】
0005712306.tif 000027
および
【数２７】
0005712306.tif 000028
は、（ａｘ´＋ｂｙ´＋ｃ）２項の線形結合として記述されなければならない。そうでなければ、式５の積分を解析的に解くことができない。これは、以下の式でなされる。
【数２８】
0005712306.tif 000029
ここで、ビームは点１から点２へ移動し、線分１は線ｊの最後の線分であると仮定する。このように、線分は逆向きに加算される。
【００６０】
ここで、
【数２９】
0005712306.tif 000030
は、指数項ｉのｘの位置の絶対値、すなわち、線ｊのビーム経路座標系におけるｘの位置である。
【００６１】
全てまとめると、ビームがＭ本の線を走査した際の時間依存性の温度分布のための以下の式が得られる。
式６
【数３０】
0005712306.tif 000031
上式で、
Ｋｊは、ハッチング経路ｊのまっすぐな線分の数である。
【数３１】
0005712306.tif 000032
は、各線分ｋｊ上の指数項の数である。
【００６２】
後続のセクションにおいて、和の中の項のための解析的表現を導く。しかしながら、上述のＴ´（ｘ，ｙ，ｚ，ｔ）の式を用いることにより、あらゆる種類のビーム経路に対して事実上温度を計算でき、また、計算は複数のＣＰＵ構成において効果的に行うことができて、これは、リアルタイムで計算を行うことができることを意味していることを前置きしておきたい。
【００６３】
ガウス関数の性質と積分
式６の式を解くために、ガウス関数のいくつかの性質を知っておく必要がある。
１．２つのガウス関数の積は、別のガウス関数となる。
【数３２】
0005712306.tif 000033
２．１つのガウス関数の積分
【数３３】
0005712306.tif 000034
【００６４】
Ｔ´（ｘ，ｙ，ｚ，ｔ）を計算するための積分
まず、ｚ方向の積分を考える。
【数３４】
0005712306.tif 000035
上式で、
【数３５】
0005712306.tif 000036
次に、ｘおよびｙの積分を考える。
【数３６】
0005712306.tif 000037
【００６５】
全ての線分が平行である場合において、座標系を変換してハッチング線に位置を合せることは容易であるため、ｘとｙとを区別する必要はない。したがって、以下の例において、全ての線はｘ軸に対して平行であると仮定される。
【数３７】
0005712306.tif 000038
上式で、
【数３８】
0005712306.tif 000039
【００６６】
線分が平行でなく任意の方向を有する場合、代数表現がもう少し複雑になる。このような場合、まずｘの積分を考える。
【数３９】
0005712306.tif 000040
上式で、
【数４０】
0005712306.tif 000041
0005712306.tif 000042
したがって、
【数４１】
0005712306.tif 000043
ここでｙの積分を考える。
【数４２】
0005712306.tif 000044
上式で、
【数４３】
0005712306.tif 000045
【００６７】
全ての式のまとめ
【数４４】
0005712306.tif 000046
指数項の位置
【数４５】
0005712306.tif 000047
上式で、
【数４６】
0005712306.tif 000048
平行線
【数４７】
0005712306.tif 000049
上式で、
【数４８】
0005712306.tif 000050
ハッチング線の任意の方向
【数４９】
0005712306.tif 000051
上式で、
【数５０】
0005712306.tif 000052
【数５１】
0005712306.tif 000053
【００６８】
ハッチング線沿ったビームのパラメータの計算
（仮想の）ビームがハッチング経路に沿って走査すると、スポット周辺の温度は、ここで式６の式から、および、先のハッチング線の温度プロファイルを予め計算したガウス関数を代入することにより計算できる。
【００６９】
温度を知っておくことにより、また、異なる条件のビームのパラメータのために最適化されたデータにアクセス手段を有することにより、ビームのエネルギー入力（すなわち特定のエネルギー付与）を適切に調整できる。
【００７０】
実施例
以下のハッチング例において（図７参照）、一定のビーム出力により台形部分が溶融され、ビームの速さが変化されて、一定の溶融深度および溶融幅を有するようにする。意図されるビーム経路とは、ビームが、図７において、左から右へ向かうのに続いて右から左へと方向を変えながら、底部から頂部へ向かって線を走査し始めるというようなものである。
【００７１】
スポットサイズは、融合前のその部分の温度であるＴｓｕｒｆに最適化されて、溶融たまりの最大温度がＴｍａｘに制限されるようにする。このことは、第１のハッチング線が、一定の速さおよびスポットサイズ固定で走査されることを意味する。その他の線は全て、同じスポットサイズおよび出力で走査されるが、速さは異なり、変化する。意図されるビーム経路に沿って配置された各計算点における速さは、まず点周辺の温度分布を計算することにより、次にデータベース内の速さ対温度データから、得ることができる。
データベース内の速さは、特定のビーム設定（出力およびスポットサイズ）および温度のために最適化されて、溶融深度および溶融幅が全ての線において同じとなるようにする。各ハッチング線の終端において、仮想のビームによって生成された温度プロファイルは、データベースから取得されたガウス関数によりモデル化される。データベースの温度の範囲は、ＴｓｕｒｆからＴｍｅｌｔまでとされ、予め計算されたデータの温度ステップは、２０Ｋに設定された。参照テーブル手順が使用されて、最も近い速さ、および、計算された温度のためのガウス関数を選んだ。
【００７２】
結果として得られる各ハッチング線に沿った速さプロファイルは、図８に示されている。これらのプロファイルは、意図されるビーム経路沿った局所的な温度分布の計算に基づいており、少なくとも１層の選択された区域を融合する際に意図されるビーム経路において使用されるビームの特定のエネルギー付与に対して決定された操作スキームに対応する。ここで、この例における特定のエネルギー付与は、ビームの速さを変化させることにより変化される。
【００７３】
上記の例において、段階的な手順は、線に沿った温度および速さを取得するために使用される。このことはまず、線に沿った特定の点における温度は平行線のための式６を使用して計算されることを意味する。次に、速さは、データベースを参照テーブルとして用いることにより温度から取得される。ハッチング線に沿った次の点は、固定距離Δｒにて計算され得て、ここに、時間ステップはΔｒ／Ｓｐｅｅｄに等しい。しかしながら、温度勾配は、時間および空間座標に対していくぶん劇的に変化するため、固定距離手順は十分ではなかった。ある場所において、小さいステップを要する一方で、その他の場所において、いくぶん長いステップでも十分性格であることもある。代わりに、速さにおいて許容される変化の最大値を用いた。これにより、温度の差の許容される最大値が取得でき、時間および空間の両方に関して温度の微分係数を数値的に計算することにより、空間ステップの許容される最大値が取得できる。
【００７４】
導出されたアルゴリズムは非常に効率的であり、リアルタイム計算において最大数千本のハッチング線を含めるのに問題はない。リアルタイムの用語は、ハッチング線に沿った速さを計算する時間が、実際の溶融時間未満となる計算を指す。
【００７５】
本発明は上で述べられた実施例に限定されないが、特許請求の範囲内で様々な方法で変更できる。例えば、ビームのパラメータを最適化し、データベースを作成する際に、より詳細かつ複雑な溶融プロセスの記述を用いることもできる。粉末は、溶融エンタルピーおよび溶融−凝固プロセスの詳細なモデルとあわせて、不均質材料としてモデル化できる。
【００７６】
式６による計算は、例えば、最小の溶融／融合時間に関してハッチングの方策を最適化するために使用できる。そのような最適化のために、実用上の見地から可能である限りは、全ての計算をリアルタイムで行う必要はない。しかしながら、実行できるハッチングの方策のための計算がリアルタイムで行うことが可能ならば役立つこともある。例として、最適化ステップから取得したデータを全て保存する必要はない。かわりに、最適化ステップの間に保存すべき情報は、例えば、ハッチング角度、ハッチング線間の距離、部分に対するハッチングの位置などに制限され得る。
【００７７】
上述の方法は、例えば、国際公開特許公報第２００４／０５６５１１号に記述されたような、構築中の部分を特定の温度に維持するために必要なビーム出力の計算に使用される方法と組み合わせることができる。こうして、ここに記載された方法は融合中にビームにより提供される局所的なエネルギーまたは出力付与の制御に使用されるが、全体のエネルギー入力は、その部分の幾何形状を含むエネルギー平衡の計算から計算できる。
【００７８】
上述の方法は、ハッチング線に沿った局所的な温度（分布）を取得するために、一様な材料のモデルを使用する。しかしながら、局所的な材料特性の差は、異なる場所に異なるＤの値を用いることによりモデル化できる。例えば、非常に薄い断面は、より低い熱伝導率を有するとしてモデル化できる。そのような断面であっても、最適化されたデータによりデータベースを拡張するための方法に制限はない。同様に、粉末床の下方の層が、その熱的特性が粉末床のものとは異なることがある調整可能な作業テーブルの近くに配置されていることを考慮に入れることができる。
【００７９】
操作スキームを選択された区域の当該部分に対して計算し決定する前に、選択された区域の一部分にだけ意図されるビーム経路を確立することができる。さらに、完全に確立された意図されるビーム経路の一部分のみに対して操作スキームを計算し決定できる。少なくとも１つの層の選択された区域を融合するステップは、意図されるビーム経路を確立するステップ、温度を計算するステップなどの間、選択された区域の融合されていない部分が今も進行中であることを考慮して、開始されてもよい。さらに、少なくとも１層の粉末層は、２つ以上の選択された区域を備えていてもよい。これらの選択された（部分）区域は、異なる形状を有していてもよく、別々に取り扱うことができる。
【００８０】
上に説明したように、意図されるビーム経路に沿った温度の計算において、同じ時点の温度計算を行う際に、特定の時点まで経路に沿って（仮想の）ビームにより付与されたエネルギーが考慮されている。このようにして、温度の蓄積が適切に考慮に入れられている。
【００８１】
記述の例において、意図されるビーム経路沿って配置された多数の位置において温度の計算が行われ、これらの位置のそれぞれにおいて局所的な温度分布が計算される。さらに、局所的な温度分布は、仮想のビームの位置の１つステップ先の位置において計算される。この単一ステップだけいくらか先の次の位置へビームを移動する際に使用される特定のエネルギー付与は、データベースから取得される。データベースは、異なる局所的な温度分布（粉末状の材料を使用および特定の融合条件用）に対する多数の予め計算された特定のエネルギー付与（すなわち上述の例におけるビームの速さ）を含む。次の位置における計算された局所的な温度分布は、適切な値またはデータベースからの値を選択するために使用される。
【図１】
0005712306.tif 000054
【図２】
0005712306.tif 000055
【図３】
0005712306.tif 000056
【図４】
0005712306.tif 000057
【図５】
0005712306.tif 000058
【図６】
0005712306.tif 000059
【図７】
0005712306.tif 000060
【図８】
0005712306.tif 000061