(57)【特許請求の範囲】
【請求項１】
プラズマジェット、プラズマＣＶＭ、ＰＡＣＥ、イオンビーム、またはガスクラスターイオンビームによる局所加工ツールを用い、前記局所加工ツールと、ワークの複数の加工点とを相対的に移動させて前記ワークを加工する加工方法であって、
前記加工点のそれぞれにおける第一の加工量を求める工程と、
前記加工点のそれぞれにおける第一の加工時間を求める工程と、
前記第一の加工時間に従い加工した場合の、前記局所加工ツールが前記加工点にある時の前記加工点のワークの温度を求める工程と、
前記求めたワークの温度に対する加工レートを、予め求めておいた前記ワークの温度に対する加工レートから求め、前記求めた加工レートと前記第一の加工時間とから、前記加工点のそれぞれにおける第二の加工量を求める工程と、
前記第一の加工量と前記第二の加工量との差である加工誤差量に応じて、前記加工点のそれぞれにおける第二の加工時間を設定する工程と、
前記加工点のそれぞれにおける前記第二の加工時間に基づき、前記ワークを加工する工程と、
を有することを特徴とする加工方法。
【請求項２】
前記第二の加工量を求める工程は、ワークの温度変化による熱膨張も考慮されることを特徴とする請求項１記載の加工方法。
【請求項３】
前記ワークの温度は、前記ワークの温度の時間変化Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）であることを特徴とする請求項１または２記載の加工方法。
【請求項４】
前記第一の加工時間に従い加工した場合の、前記局所加工ツールが前記加工点にある時の前記加工点のワークの温度を求める工程は、シミュレーションにより求めることを特徴とする請求項１乃至３いずれか１項記載の加工方法。
【請求項５】
請求項１乃至４いずれか一項記載の加工方法を用いて、半導体ウェハー、レンズまたはミラーである素子、または金型による素子を製造することを特徴とする素子の製造方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は半導体ウェハー、レンズやミラー、金型等の高精度素子を製造するための、局所加工方法に関する。
【背景技術】
【０００２】
被加工物を所望の形状に加工する方法として、局所加工ツールの滞留時間を、所望加工量に応じて変動させながら走査する方法が知られている。つまり、所望加工量が多い加工点では局所加工ツールの滞留時間を相対的に長く、所望加工量が少ない加工点では局所加工ツールの滞留時間を相対的に短くして加工する。このような局所加工方法を滞留時間制御加工と呼ぶ。
【０００３】
この方法は次の工程からなる。予め局所加工ツールの加工レートを、ワークと同質のダミーワークを加工することで求めておく。次に、ワークの加工前形状を測定し、加工前形状と設計形状との差から、所望加工量の分布を算出する。そして、所望加工量の分布と加工レートから、局所加工ツールの滞留時間の分布データを作成する。この滞留時間分布を求める演算はデコンボリューションと呼ばれており、シミュレーション法や汎関数フィルターＦＦＴ法（非特許文献１）、Ｖａｎ−Ｃｉｔｔｅｒｔ法、Ｌｕｃｙ−Ｒｉｃｈａｒｄｓｏｎ法など様々な手法がある。そして最後に、滞留時間の分布データに従い、局所加工ツールとワークとを相対運動させて加工を行う。実際の加工時間は滞留時間と、加工点間を移動するのに要する時間である移動時間との和となるので、実際の加工形状はあらかじめ定めた設計形状と平行な形状となる。実際の加工形状と設計形状との差分は、加工レートと移動時間が一定であれば容易に求めることができる。従って、移動中の加工量を予め考慮し、加工後のワークが所望形状となるように加工量を算出し、その上でワークを加工することで実際の加工形状と設計形状との差分が少ないワークを得ることができる。
【０００４】
近年、高速かつ低ダメージな局所加工が可能な技術として、局所プラズマによる化学的な除去加工技術が注目されている。方式は種々あるが、プラズマジェット、プラズマＣＶＭ（Ｃｈｅｍｉｃａｌ  Ｖａｐｏｒｉｚａｔｉｏｎ  Ｍａｃｈｉｎｉｎｇ）、ＰＡＣＥ（Ｐｌａｓｍａ  Ａｓｓｉｓｔｅｄ  Ｃｈｅｍｉｃａｌ  Ｅｔｃｈｉｎｇ）などが代表的である。
【０００５】
これらのうち、プラズマジェットを利用して所望の加工形状を得る方法が非特許文献２に開示されている。プラズマジェット加工等の局所加工における加工レートはワーク表面の温度に依存する。さらにプラズマジェット加工ツールの放出する熱量が大きいためワークの温度は、ワークに対するプラズマジェット加工ツールの滞留時間に依存する。
【０００６】
ワークの加工レートは、図６（Ｓｉの加工レートの温度依存性を示した図）のように、ある温度まではワーク温度の上昇に伴って上昇するが、特定の温度以降は温度に対して加工レートは大きく変動せずに飽和する傾向がある。
【０００７】
そこで、加工レートが滞留時間に対して一定とみなせる温度まで、ワークをヒーター等で加熱することで、プラズマジェットによる局所加熱の効果を小さくし、所望の加工形状を得る方法が従来とられていた（非特許文献２）。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】根岸ほか  精密工学会誌  第６２巻  第３号（１９９６）４０８−４１２
【非特許文献２】Ｇ．Ｂｏｅｈｍ  ｅｔ．ａｌ．，Ｐｒｏｃ．ｏｆ  ｔｈｅ  ９ｔｈ  ＩＣＰＥ（１９９９）２３１−２３６
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、（非特許文献２）の方法では、プラズマジェット加工ツールなど放出する熱量の大きい局所加工ツールを使用すると、加工を進めるうちにワークの温度が大きく変動してしまう。従って、ワークの温度をヒーターなどで加熱しておいても加工レートが加工点毎に無視できないほど異なるため、形状精度が悪いという問題があった。
【００１０】
そこで本発明では、加工に伴うワークの温度変化を考慮した加工精度の高い局所加工方法を提供することを目的とする。
【課題を解決するための手段】
【００１１】
前述の課題を鑑みて本発明が提供する加工方法は、プラズマジェット、プラズマＣＶＭ、ＰＡＣＥ、イオンビーム、またはガスクラスターイオンビームによる局所加工ツールを用い、前記局所加工ツールと、ワークの複数の加工点とを相対的に移動させて前記ワークを加工する加工方法であって、前記加工点のそれぞれにおける第一の加工量を求める工程と、前記加工点のそれぞれにおける第一の加工時間を求める工程と、前記第一の加工時間に従い加工した場合の、前記局所加工ツールが前記加工点にある時の前記加工点のワークの温度を求める工程と、前記求めたワークの温度に対する加工レートを、予め求めておいた前記ワークの温度に対する加工レートから求め、前記求めた加工レートと前記第一の加工時間とから、前記加工点のそれぞれにおける第二の加工量を求める工程と、前記第一の加工量と前記第二の加工量との差である加工誤差量に応じて、前記加工点のそれぞれにおける第二の加工時間を設定する工程と、前記加工点のそれぞれにおける前記第二の加工時間に基づき、前記ワークを加工する工程と、を有することを特徴とする。
また、本発明が提供する素子の製造方法は、上記加工方法を用いて、半導体ウェハー、レンズまたはミラーである素子、または金型による素子を製造することを特徴とする。
【発明の効果】
【００１２】
本発明の局所加工方法によれば、加工に伴うワークの温度変化を考慮して滞留時間分布を算出するので、ワークの温度変化により加工レートが変化する場合でも、ワークを所望の形状へ高精度に加工ができる。
【００１３】
また、本発明の局所加工方法によれば、加工に伴うワークの温度変化による、ワークの熱変形も考慮して滞留時間分布を算出するので、更に高精度な加工ができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態１において、滞留時間分布を求めるフローを示す図である。
【図２】本発明に基づく、加工前形状と目標加工形状との関係を示す図である。
【図３】本発明の実施形態２において、滞留時間分布を求めるフローを示す図である。
【図４】本発明に基づく加工を行う加工装置の概要を示す図である。
【図５】本発明に基づく加工の実施例と比較例とを比較する表である。
【図６】Ｓｉの加工レートの温度依存性を表す図である。
【発明を実施するための形態】
【００１５】
（実施形態１）
本発明に基づく局所加工方法の第一の実施形態を、図１を用いて説明する。
【００１６】
（工程１）まず、加工前のワーク表面形状を加工前形状として、触針式や光干渉式などの形状測定器を用いて測定する。測定した加工前形状と設計形状との差を、所望の加工量Ｖ０（ｘ，ｙ）として求める。ｘ、ｙは定められた原点からのワーク上の座標をしめしており、したがって、例えばＶ０（ｘ１，ｙ１）はある特定の座標（ｘ１，ｙ１）での加工量を示している。したがってＶ０（ｘ，ｙ）は、ワーク上の特定の位置での加工量のいわば“分布”を示している。
【００１７】
尚、ここで求めた加工量Ｖ０（ｘ，ｙ）は全ての加工点で０以上の値をとる暫定値である。この加工量Ｖ０（ｘ，ｙ）を第一の加工量とする。
【００１８】
暫定値とする理由を図２に基づき説明する。上述の加工量Ｖ０（ｘ，ｙ）は加工前形状２０１と設計形状２０２との差分として定められる。局所加工ツールのワークに対する滞留時間が０であっても加工点間を移動するのに要する時間、加工は行われるので、目標加工形状は設計形状２０２の曲面に平行な形状２０３を設定する。ここに、目標加工形状２０３と設計形状２０２との差分をベース加工量と呼ぶこととする。加工レートが一定であれば、移動時間と加工レートからベース加工量を求め、目標加工形状２０３を一意的に設定することができる。ところが、加工に伴うワークの温度変化によって加工レートが変動する場合、加工点間を移動する間の加工量は、当該移動を行うまでの加工履歴によってワーク温度が異なるので、加工点毎に異なる。従って、滞留時間分布が決定していない段階で、ベース加工量を一意的に設定することができない。つまり、滞留時間分布と目標加工形状２０３は、同時に決定しなければならないので、加工量Ｖ０（ｘ，ｙ）は暫定値とする（以下の工程では計算の初期値）。
【００１９】
加工量Ｖ０（ｘ，ｙ）は暫定値であるので、加工量Ｖ０（ｘ，ｙ）に含まれるベース加工量の初期値は０としておいて良い。あるいは、予め加工領域全域の滞留時間を０として、移動時間だけの加工時間で加工を実際に行い、加工形状を測定して得た最大の加工深さを、ベース加工量の初期値としても良い。この場合には、滞留時間分布を求める計算の所要時間を短縮できる。
【００２０】
ワーク上の座標、温度、及び加工レートとの相関データである温度・加工レートＵ（ｘ，ｙ，Ｔ）は、予め実験により取得しておく。加工レートがワーク上の座標によらず一様とみなせる場合は、ワーク温度と加工レートとの相関データとして温度・加工レートＵ（Ｔ）を取得して、以下の工程にて用いても良い。
【００２１】
具体的な温度・加工レートの取得方法は、ワークの温度を放射温度計などの温度測定手段により測定しながら、局所加工ツールを一方向に一定速度で走査したり、同じ場所で所定の時間加工したりして加工を行い、加工形状を測定して求める。
【００２２】
（工程２）続いて、上述のＶ０（ｘ，ｙ）を初期値とする加工量Ｖ（ｘ，ｙ）と温度Ｔｋ１における温度・加工レートＵｋ１（ｘ，ｙ，Ｔｋ１）から、ワークの加工時間分布τ（ｘ，ｙ）を公知のデコンボリューション手法により求める。尚、加工時間は移動時間ｔ０未満にはできないので、デコンボリューションの過程で、全ての加工点において移動時間ｔ０以上になるように修正する。ここに温度Ｔｋ１は任意であるが、加工レートの極端に大きい、または小さい温度を選択すると、滞留時間を求める計算の所要時間が長くなる。そこで、温度を考慮した加工レートを取得した際に得られた温度の範囲から選択すると良い。このようにして第一の加工量とワークの温度・加工レートから算出されたワークの加工時間分布τ（ｘ，ｙ）を、第一の加工時間分布とする。
【００２３】
（工程３）そして、第一の加工時間分布であるワークの加工時間分布τ（ｘ，ｙ）に従って加工した際の、ワーク温度またはワーク温度の時間変化を次のように求める。一つは、熱源が移動する非定常の熱シミュレーションを行う方法である。もう一つは、加工中の温度を測定しながら、加工時間τ（ｘ，ｙ）に従ってダミー加工を行い、各加工点の温度変化を関数化して求める方法である。このように求めた温度は、加工ツールが加工点（ｕ，ｖ）にあるときの、時刻ｔにおける座標（ｘ，ｙ）の温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）である。尚、時刻ｔは０以上、τ（ｕ，ｖ）以下の値をとる。このようにして第一の加工時間分布に従い加工した際のワーク温度を求める。
【００２４】
（工程４）次に、局所加工ツールによる実際のワーク加工量にあたる、第二の加工量Ｖ’（ｘ，ｙ）は、ワーク温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）と、ワーク温度と加工レートの関係データである温度・加工レートＵ（ｘ，ｙ，Ｔ）と、上述のように暫定的に定めた第一の加工時間分布τ（ｘ，ｙ）から、次の（式１）で求める。尚、Ａで示した積分領域は加工領域である。
【００２５】
この場合ワーク温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）はワークの温度の時間変化である。
【００２６】
【数１】
0005814570.tif 000002
【００２７】
より簡便に実際の加工量である第二の加工量Ｖ’（ｘ，ｙ）を求めるには、局所加工ツールが加工点（ｕ，ｖ）にあるときの、時刻ｔにおける加工点（ｕ，ｖ）の温度をＴ（ｕ，ｖ，ｔ）として、次の（式２）で求める。これは、近似計算により計算を高速化する場合だけでなく、ワーク自体の大きさが局所加工ツールに対して比較的小さいなど、ワークの温度が一様で局所加工ツールの相対位置と加工時間のみでワーク温度が決まる場合にも適用可能である。
【００２８】
【数２】
0005814570.tif 000003
【００２９】
更に、簡便に第二の加工量Ｖ’（ｘ，ｙ）を求めるには、加工ツールが加工点（ｕ，ｖ）の加工を開始する時の温度をＴ（ｕ，ｖ）として、次の（式３）で求める。
【００３０】
【数３】
0005814570.tif 000004
【００３１】
（工程５）次に、暫定的に定めたＶ０（ｘ，ｙ）を初期値とする第一の加工量Ｖ（ｘ，ｙ）と、実際の加工量にあたる第二の加工量Ｖ’（ｘ，ｙ）との差を加工誤差量ε（ｘ，ｙ）として求める（式４）。
【００３２】
【数４】
0005814570.tif 000005
【００３３】
加工誤差量ε（ｘ，ｙ）が予め設定した目標誤差量δ以下であれば、決定された第一の加工時間分布τ（ｘ，ｙ）から移動時間ｔ０を引いて、実際の滞留時間分布に相当する第二の加工時間分布Ｄ（ｘ，ｙ）を求める。目標誤差量δを超えている場合には、加工誤差量ε（ｘ，ｙ）が収斂回数に対して収束したかを判定する。収束した場合には、加工誤差量ε（ｘ，ｙ）の負値のうち、絶対値が最大となる値の絶対値を、所望加工量Ｖ（ｘ，ｙ）に加えて所望加工量を修正し、工程２に戻る。収束していない場合には、加工誤差量ε（ｘ，ｙ）と温度Ｔｋ２における温度・加工レートＵｋ２（ｘ，ｙ，Ｔｋ２）から、修正時間ｄ（ｘ，ｙ）を公知のデコンボリューション手法により求める。そして、修正時間ｄ（ｘ，ｙ）を加工時間τ（ｘ，ｙ）に加え、更に、全ての加工点で移動時間ｔ０以上となるように、加工時間を修正して、実加工時間を設定する。その修正時間は、加工量Ｖが実際の加工量である第二の加工量よりも多い、加工過多の加工点では負値を取り、逆に加工不足の加工点では正値を取る。従って修正時間は加工誤差量を０に近づける作用をする。そのため温度Ｔｋ２は、加工誤差量を最小化する収斂計算が、適正に収束するよう適宜選択する。加工レートの小さい温度を選択すると、収束は早くなるが、加工誤差量が振動しやすくなる。加工レートが大きい温度を選択すると、収束は遅くなる。加工時間の修正が完了したら、工程３に戻る。
【００３４】
以上に述べた方法で第二の加工時間分布Ｄを求め、該加工時間分布に基づき局所加工ツールがワークを加工する。局所加工ツールの形態の一例は、ワークの構成元素と反応して、揮発性物質を生成する元素を含む、粒子をワークに照射することで除去加工を行うものである。具体的には、プラズマジェット、プラズマＣＶＭ、ＰＡＣＥなど、ラジカルを利用するものである。または、イオンビーム、ＧＣＩＢ（Ｇａｓ  Ｃｌｕｓｔｅｒ  Ｉｏｎ  Ｂｅａｍ）など荷電粒子を利用するものである。無論、本発明の局所加工ツールはこれらに限定されるものではなく、加工レートがワーク温度によって変化する如何なる形態の局所加工ツールであっても適用可能である。
【００３５】
また、加工点の移動と滞留を繰り返す、断続的な相対運動による加工方法を述べたが、本発明の加工方法はこれに限定されるものではない。すなわち、加工点の通過時間が当該加工点の所望加工時間と一致するように、走査速度を連続的に変調して加工する連続的な相対運動による加工方法であっても良い。
【００３６】
（実施形態２）
本発明に基づく局所加工方法の第二の実施形態を、図３を用いて説明する。第一の実施形態との差異は、加工レートと実際の加工形状をワークの熱変形を考慮して求める点にある。
【００３７】
第一の実施形態と同様の方法で、ワークの温度Ｔと加工レートの関係データである温度・加工レートＵ（ｘ，ｙ，Ｔ）を予め求めておく。ここで測定した温度・加工レートＵ（ｘ，ｙ，Ｔ）は、ある温度でのワークの形状をもとにしている。しかし、実際の加工中のワークの形状は、熱変形のため、これとは異なっている場合がある。そこで、温度・加工レートＵ（ｘ，ｙ，Ｔ）を、加工中に測定した温度とワークの熱膨張係数を基に、実際の加工中の温度・加工レートＵ’（ｘ，ｙ，Ｔ）を予め求めておき、上述のＵに換えてＵ’に変換する（工程１’）。熱膨張による変形がワークに、生じるとワークが膨張した部分を追加的に加工する必要があり、たとえば、あるワークが等方的な熱膨張特性を有してその熱膨張係数をαとすると、Ｕ（αｘ，αｙ，Ｔ）＝Ｕ’（ｘ，ｙ，Ｔ）として加工を行う。
【００３８】
以後、上記熱変形の分を加算した温度・加工レートＵ’（ｘ，ｙ，Ｔ）を利用して、（工程２）、（工程４’）および（工程５）を行う。
【００３９】
実際の加工量を求める（工程４’）では、熱変形を考慮した温度・加工レートＵ’（ｘ，ｙ，ｔ）を重ね合わせてなる加工量Ｖ”（ｘ，ｙ）を求める。
【００４０】
その上で、（工程３）で求めたワーク温度と、ワークの熱膨張係数を基に、ある温度での（例えば室温）実際の加工量にあたる第二の加工量Ｖ’（ｘ，ｙ）を求める。
【実施例】
【００４１】
（実施例１）
本実施例では、（実施形態１）に基づき、プラズマジェットで平板の石英基板に、有効径４０ｍｍ・曲率半径１０ｍの凹球面を形成した例を示す。
【００４２】
使用した加工装置の動作について図４に基づき説明する。プラズマジェット加工ツール４０２はフランジ４１６を介してチャンバー４１５内に設置されている。流量制御装置４１０により流量を制御したアルゴン、酸素、六フッ化硫黄をガス供給設備４０９からプラズマジェット加工ツール４０２に供給する。これらのガスを高周波電源４０８より印加した高周波電力により電離し、Ｆラジカルなどの活性種を含むプラズマジェット４０１を生成する。そして、該プラズマジェットをワーク４０３に照射して加工を行う。ワークホルダー４０４により保持したワークは、Ｘ軸駆動機構４０５、Ｙ軸駆動機構４０６、Ｚ軸駆動機構４０７を連接してなるワーク駆動機構４１８に搭載されている。該ワーク駆動機構によりワークとプラズマジェットとの相対走査が可能となっている。ワーク駆動機構はチャンバー外に設置の制御用コンピュータ４１４とフランジ４１７を介して接続されており、制御用コンピュータより発せられる制御信号に従って動作する。チャンバーは真空ポンプ４１３により排気が可能で、コンダクタンス可変バルブ４１２を併用することでチャンバー内の圧力を制御することができるようになっている。また、チャンバーにはビューポート４１１が設けてあり、ワークやプラズマジェットを観察することができる構成となっている。
【００４３】
凹面形状の加工を行うに先立ち、ワーク温度と加工レートの関係データである温度・加工レートＵ（ｘ，ｙ，Ｔ）を求めるため、ダミーの石英基板に対して加工を行った。加工はプラズマジェット直下のワーク表面温度を不図示の放射温度計でビューポート４１１から測定しながら、０．５ｍｍ／ｓｅｃ以上２ｍｍ／ｓｅｃ以下の範囲から選択した種々の速度で１方向に走査して行った。加工時の温度は走査速度に依存し、およそ３００℃から５００℃であった。加工前後の基板表面形状を光干渉式の形状測定器で測定し、加工断面形状から加工レートを求めた。このようにして測定した加工時の温度Ｔと温度・加工レートＵの関係をフィッティングし、Ｕ＝ｆ（ｘ，ｙ，Ｔ）なる関数で表した。
【００４４】
次に、光干渉式の形状測定器で測定したワークの加工前形状と、設計形状との差をとり、第一の加工量Ｖ０（ｘ，ｙ）を求めた（工程１）。
【００４５】
（工程２）及び（工程５）で用いる温度・加工レートＵｋ１およびＵｋ２は４８０℃における加工レートを採用した。また、デコンボリューションは汎関数フィルターＦＦＴ法により行った（工程２）。加工時間分布τ（ｘ，ｙ）からワーク温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）を求めるため、熱源が移動する非定常の熱シミュレーションを行った（工程３）。実際の加工形状Ｖ’（ｘ，ｙ）は（式２）を基に求めた（工程４）。（式２）において使用するデータは、ワーク温度と加工レートの関係データである前述のＵ＝ｆ（ｘ，ｙ，Ｔ）を、温度Ｔの時間変化として前記熱シミュレーションの結果を、それぞれ入力した。実際の積分計算は区分求積法で行った。
【００４６】
以上に述べた条件で収斂計算を実施し、第二の加工時間分布である滞留時間分布Ｄ（ｘ，ｙ）を求めた。
【００４７】
ここに求めた滞留時間分布のデータを制御用コンピュータ４１４に入力し、各加工点の座標にプラズマジェットが滞留する時間がちょうど第二の加工時間分布Ｄ（ｘ，ｙ）のデータと一致するよう、プラズマジェットを発生させつつワーク駆動機構を制御してワークに凹面形状を形成する加工を行った。
【００４８】
加工後の石英基板の形状を光干渉式の形状測定器で測定し、予想加工形状と比較したところ、ＰＶ（Ｐｅａｋ  ｔｏ  Ｖａｌｌｅｙ）にて８％の加工誤差であった。
【００４９】
（実施例２）
本実施例では、（実施形態２）に基づき、プラズマジェットで平板の石英ウェハーをワークの熱変形を考慮して平坦化した例を示す。
【００５０】
まず（工程１’）において、加工中のワークの熱変形を考慮して実際の加工中の温度・加工レートＵ’（ｘ，ｙ，Ｔ）を求めるため、ダミーの石英ウェハーに対して加工を行った。加工はプラズマジェット直下のワーク表面温度分布を不図示のサーモグラフィーでビューポート４１１から測定しながら、０．５ｍｍ／ｓｅｃ以上２ｍｍ／ｓｅｃ以下の範囲から選択した種々の速度でワークを１方向に走査して行った。加工前後の基板表面形状を光干渉式の形状測定器で測定し、加工断面形状から加工レートを求めた。ここに求めた加工断面形状および加工時の温度と、石英の熱膨張係数から加工時の温度における加工断面形状を計算した。そして、加工時の温度における加工断面形状を基に実際の加工中の温度・加工レートＵ’（ｘ，ｙ，Ｔ）を求めた。そして、加工時の温度Ｔと実際の加工中の温度・加工レートＵ’との関係をＵ’＝ｇ（ｘ，ｙ，Ｔ）なる関数で表した。
【００５１】
（工程２）と（工程５）で用いる温度・加工レートＵ’ｋ１およびＵ’ｋ２は４８０℃における実際の加工中の加工レートを採用した。また、デコンボリューションは汎関数フィルターＦＦＴ法により行った。第一の加工時間分布からワーク温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）を求めるため、ワークに対して熱源が移動するように設定された非定常の熱シミュレーションを行った。
【００５２】
（工程４’）において、まず、先に求めた実際の加工中の温度・加工レートＵ’＝ｇ（ｘ，ｙ，Ｔ）、ワーク温度Ｔ（ｕ，ｖ，ｘ，ｙ，ｔ）を（式２）に入力し、加工形状Ｖ”（ｘ，ｙ）を得た。実際の積分計算は区分求積法で行った。次に、加工形状Ｖ”（ｘ，ｙ）とワーク温度とワークの熱膨張係数を基に室温まで冷却後の加工形状を計算し、実際の加工量にあたる第二の加工量Ｖ’（ｘ，ｙ）を求めた。このときのワーク温度は各加工点での滞留が終了した時刻ｔｆにおける加工点の温度Ｔ（ｕ，ｖ，ｕ，ｖ，ｔｆ）を採用した。
【００５３】
以上に述べた条件で収斂計算を実施し、第二の加工時間分布である滞留時間分布Ｄ（ｘ，ｙ）を求めた。更に、更新された所望加工量Ｖ（ｘ，ｙ）と加工前形状から、予想加工形状を求めた。
【００５４】
ここに求めた滞留時間分布のデータを制御用コンピュータ４１４に入力し、各加工点の座標にプラズマジェットが滞留する時間がちょうど第二の加工時間分布である滞留時間分布Ｄ（ｘ，ｙ）のデータと一致するよう、ワーク駆動機構４１８を制御してワークを平坦化する加工を行った。
【００５５】
加工後の石英基板の形状を光干渉式の形状測定器で測定し、予想加工形状と比較したところ、ＰＶ（Ｐｅａｋ  ｔｏ  Ｖａｌｌｅｙ）にて６％の加工誤差であった。
【００５６】
（比較例）
４００℃に加熱した平板の石英基板に、プラズマジェットにより有効径４０ｍｍ・曲率半径１０ｍの凹球面を形成した例を示す。
【００５７】
まず、加工レートを求めるため、４００℃に加熱したダミーの石英基板の中央を２秒間照射した。加工完了後に加工形状を測定し、これを加工時間で除して加工レートを求めた。ここに求めた加工レートはワーク温度に依存せず、加工中に不変であるとみなした加工レートである。そして、加工前形状と設計形状との差から求めた第一の加工量を前記加工レートによりデコンボリュートする演算を行い、滞留時間分布を求めた。デコンボリューション手法は汎関数フィルターＦＦＴ法を適用した。また、予想加工形状は該滞留時間分布と前記加工レートとのコンボリューション演算により求めた。
【００５８】
このようにして求めた滞留時間分布のデータを制御用コンピュータ４１４に入力し、各加工点の座標にプラズマジェットが滞留する時間がちょうど前記滞留時間分布のデータと一致するよう、ワーク駆動機構を制御してワークに凹面形状を形成する加工を行った。
【００５９】
加工後の石英基板の形状を光干渉式の形状測定器で測定し、予想加工形状と比較したところ、ＰＶ（Ｐｅａｋ  ｔｏ  Ｖａｌｌｅｙ）にて１５％の加工誤差であった。
【００６０】
以上に述べた実施例と比較例とを比較する表を図５に示す。加工に伴うワーク温度の変化と、ワーク温度の変化による加工レートの変化とを考慮して求めた滞留時間分布に従い加工を行うことで、大幅に加工精度が向上した。
【産業上の利用可能性】
【００６１】
本発明の加工方法は、プラズマジェット加工等の放射する熱量の大きな局所加工ツールを大口径レンズなど比較的大きいワークを高精度に加工する際に好適に用いられる。
【符号の説明】
【００６２】
２０１  加工前形状
２０２  設計形状
２０３  目標加工形状
４０１  プラズマジェット
４０２  プラズマジェット加工ツール
４０３  ワーク
４１４  制御用コンピュータ
４１８  ワーク駆動機構
【図１】
0005814570.tif 000006
【図２】
0005814570.tif 000007
【図３】
0005814570.tif 000008
【図４】
0005814570.tif 000009
【図５】
0005814570.tif 000010
【図６】
0005814570.tif 000011