JP2003037462A - Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator - Google Patents

Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator

Info

Publication number
JP2003037462A
JP2003037462A JP2001226772A JP2001226772A JP2003037462A JP 2003037462 A JP2003037462 A JP 2003037462A JP 2001226772 A JP2001226772 A JP 2001226772A JP 2001226772 A JP2001226772 A JP 2001226772A JP 2003037462 A JP2003037462 A JP 2003037462A
Authority
JP
Japan
Prior art keywords
piezoelectric vibrator
vibration
spurious
designing
main vibration
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.)
Withdrawn
Application number
JP2001226772A
Other languages
Japanese (ja)
Inventor
Masako Tanaka
雅子 田中
Tsutomu Imai
勉 今井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2001226772A priority Critical patent/JP2003037462A/en
Publication of JP2003037462A publication Critical patent/JP2003037462A/en
Withdrawn legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To solve the problem, with the high-frequencies of a crystal oscillator, an optimum design range which averts spurious oscillations becoming narrowed, so that design becomes very difficult. SOLUTION: An area of points, representing main oscillation and spurious oscillations of a conventional mode chart is plotted by a value, in proportion to the total charge volume included in each oscillation, and a region which averts the spurious oscillation of a large area is selected as optimum design region.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【発明の属する技術分野】本発明は、圧電振動子の設計
方法、モードチャートの作成方法、周波数温度特性曲線
の作成方法、圧電振動子の設計装置および圧電振動子の
シミュレーションプログラムに関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a piezoelectric vibrator designing method, a mode chart creating method, a frequency temperature characteristic curve creating method, a piezoelectric vibrator designing apparatus, and a piezoelectric vibrator simulation program.

【0002】[0002]

【従来の技術】ATカット水晶振動子の主振動周波数
は、矩形平板の形状を持つ振動子の結晶軸x1、x2、
x3を図1のように定義し、それぞれの寸法を長さ2
a、厚み2b、幅2cとした場合、式(1)により決定
される。
2. Description of the Related Art The main vibration frequency of an AT-cut crystal oscillator is the crystal axes x1, x2 of the oscillator having the shape of a rectangular plate.
x3 is defined as in Fig. 1 and each dimension is length 2
In the case of a, thickness 2b, and width 2c, it is determined by the equation (1).

【0003】[0003]

【数1】 [Equation 1]

【0004】これは、水晶ブランクが図2に示したよう
に、厚み2bの方向x2には滑りながら、x1方向に変
位を起こす厚み滑り振動と呼ばれるもので、周波数は、
式(1)で示されたように、厚み2bにほぼ反比例す
る。この厚み滑り振動は、温度特性において、図3
(a)で示すように、3次曲線の温度特性を持ってお
り、機器類の使用範囲である−50℃〜100℃におい
て、±10ppm程度の周波数偏差しか示さないため、
高安定な周波数源として広く使用されている。
As shown in FIG. 2, this is called thickness shear vibration in which the crystal blank is displaced in the x1 direction while sliding in the direction x2 of the thickness 2b, and the frequency is
As shown in the equation (1), it is almost inversely proportional to the thickness 2b. This thickness shear vibration is shown in FIG.
As shown in (a), since it has a temperature characteristic of a cubic curve and shows only a frequency deviation of about ± 10 ppm in the use range of the equipment of −50 ° C. to 100 ° C.,
Widely used as a highly stable frequency source.

【0005】しかし、図1のような形状を持つ水晶ブラ
ンクには、主振動以外に様々な機械振動も同時に存在
し、一般にスプリアス振動(不要振動)と呼ばれる。図
4に代表的なスプリアス振動のモードを示し、それらの
モードに対応する運動方程式から導かれた周波数の算出
式を以下に示す。ここで、式(2)は図4(a)の場合
に相当し、式(3)は図4(b)の場合に相当し、式
(4)は図4(c)の場合に相当する。
However, in a crystal blank having a shape as shown in FIG. 1, various mechanical vibrations simultaneously exist in addition to the main vibration, and it is generally called spurious vibration (unnecessary vibration). FIG. 4 shows typical spurious vibration modes, and the formulas for calculating the frequencies derived from the equations of motion corresponding to these modes are shown below. Here, the equation (2) corresponds to the case of FIG. 4A, the equation (3) corresponds to the case of FIG. 4B, and the equation (4) corresponds to the case of FIG. 4C. .

【0006】[0006]

【数2】 [Equation 2]

【0007】[0007]

【数3】 [Equation 3]

【0008】[0008]

【数4】 [Equation 4]

【0009】式(2)〜式(4)により、長さ2aや幅
2c方向に定在波が立つモードでは、周波数はブランク
の長さ2aや、幅2cに大きく影響されることがわか
る。このスプリアス振動が主振動の近くに存在する場合
には、主振動とのエネルギー的な結合が生じ、その結
果、主振動のエネルギーが取られて、インピーダンスが
上昇したり、周波数がシフトしたりする。
From the equations (2) to (4), it is understood that the frequency is greatly influenced by the blank length 2a and the blank 2c in the mode in which the standing wave stands in the length 2a and the width 2c directions. When this spurious vibration exists near the main vibration, energy coupling with the main vibration occurs, and as a result, the energy of the main vibration is taken and the impedance rises or the frequency shifts. .

【0010】特に、スプリアス振動は、主振動に比べは
るかに大きな温度係数を持っているため、図3(b)で
示したように、主振動の温度依存性(以下温度特性と称
する。) を測定する中で、急速に主振動に近づき去って
いく。特に、主振動とスプリアス振動が交差する温度で
は、主振動がジャンプしてしまい、発振器に組み込まれ
た場合には、不発現象を起こしてしまうため、振動子ブ
ランクの設計に際しては、スプリアス振動を避けるよう
に長さ2a、幅2cを決定する。このため、正確なスプ
リアス振動予測は、振動子の設計にとって重要な課題で
ある。
Particularly, since spurious vibration has a temperature coefficient much larger than that of the main vibration, the temperature dependence of the main vibration (hereinafter referred to as the temperature characteristic) is shown in FIG. 3B. While measuring, it rapidly approaches the main vibration. Especially, at the temperature where the main vibration and the spurious vibration intersect, the main vibration jumps, and if it is incorporated in the oscillator, a non-firing phenomenon occurs, so avoid spurious vibration when designing the oscillator blank. Thus, the length 2a and the width 2c are determined. Therefore, accurate spurious vibration prediction is an important issue for oscillator design.

【0011】一般に、振動子の設計には、モードチャー
トと呼ばれる、主振動とスプリアス振動のそれぞれの周
波数位置をブランク形状に対してプロットしたものが用
いられる。図5に示したのは、従来の低周波モードチャ
ートの一例であるが、縦軸は周波数、横軸は振動子の幅
2cである。本図では、その他のブランク長さ2aや電
極の寸法などは一定とした。なお、縦軸には、周波数が
電極無しの水晶ブランクの周波数で規格化された規格化
周波数を、横軸には、幅2cや長さ2aが厚み2bで規
格化された値を使うこともある。
Generally, for designing a vibrator, what is called a mode chart, in which the frequency positions of the main vibration and the spurious vibration are plotted with respect to the blank shape, is used. FIG. 5 shows an example of a conventional low frequency mode chart, in which the vertical axis represents frequency and the horizontal axis represents the width 2c of the vibrator. In this figure, the other blank lengths 2a and the dimensions of the electrodes are constant. It should be noted that the vertical axis may be a standardized frequency whose frequency is standardized with the frequency of a crystal blank without electrodes, and the horizontal axis may be a value standardized with a width 2c or a length 2a being a thickness 2b. is there.

【0012】図中、線上に位置する点は主振動で、それ
以外の点はスプリアス振動である。図中に丸で囲まれた
領域では、主振動に交わるスプリアス振動が存在しな
い。この領域に対応する2c寸法をブランク幅として設
計することにより、周波数の安定した振動子を得ること
が出来る。実際の振動子では、電極が付加されたり、振
動子が支持されたりするため、式(1)〜式(4)に算
出されたとおりに周波数が観測されることはない。
In the figure, points on the line are main vibrations, and other points are spurious vibrations. In the area surrounded by a circle in the figure, there is no spurious vibration that intersects with the main vibration. By designing the 2c dimension corresponding to this region as the blank width, it is possible to obtain a vibrator with a stable frequency. In an actual oscillator, an electrode is added or the oscillator is supported, so that the frequency is not observed as calculated by the equations (1) to (4).

【0013】従来、このようなモードチャートは、形状
毎に実験によって作成することが多かったが、近年、有
限要素法(以後 FEM(Finite Element
Method)と称する。) による高精度なATカッ
ト振動子解析が確立され (特開2000−183415
号公報) 、計算機によっても、モードチャートが作成で
きるようになった。図5もFEMにより算出された結果
である。
In the past, such a mode chart was often created by experiment for each shape, but in recent years, the finite element method (hereinafter FEM (Finite Element) is used.
Method). ) Has established a highly accurate AT-cut transducer analysis (Japanese Patent Laid-Open No. 2000-183415).
It has also become possible to create mode charts using a computer. FIG. 5 also shows the result calculated by FEM.

【0014】このように、水晶振動子のスプリアス振動
予測を効率的に行い、最適設計値を探すには、モードチ
ャートは必須のものである。
As described above, the mode chart is indispensable for efficiently performing the spurious vibration prediction of the crystal unit and searching for the optimum design value.

【0015】[0015]

【発明が解決しようとする課題】最近、水晶振動子に対
する高周波化の要求が高くなってきた。式(1)で示し
たように、高周波では、式(2)〜式(4)で示される
長さ2aや幅2cに依存する輪郭振動が高次波になって
くる。長さ2a方向や幅2c方向に高次波が存在する
と、それらが複合したスプリアス振動も数多く発生する
ため、スプリアス振動の数がさらに増す。
Recently, there has been an increasing demand for higher frequencies for crystal oscillators. As shown in Expression (1), at high frequencies, the contour vibration depending on the length 2a and the width 2c shown in Expressions (2) to (4) becomes a higher-order wave. If higher-order waves exist in the length 2a direction and the width 2c direction, many spurious vibrations that combine them are generated, and the number of spurious vibrations further increases.

【0016】図6に高周波領域における従来のモードチ
ャートを示す。図5と比較して明らかなように、高周波
になると、スプリアス振動の間隔が密になり、設計に適
した領域は非常に狭くなる。しかも、量産においては、
加工誤差もあり、スプリアス振動の温度による周波数シ
フトも考慮に入れなければならないため、最適領域を見
つけることは、不可能に近い。
FIG. 6 shows a conventional mode chart in the high frequency region. As is clear from comparison with FIG. 5, at high frequencies, the intervals of spurious vibrations become close, and the region suitable for design becomes very narrow. Moreover, in mass production,
Since there is a processing error and the frequency shift due to the temperature of spurious vibration must be taken into consideration, it is almost impossible to find the optimum region.

【0017】しかし、経験上、温度特性において算出さ
れるすべてのスプリアス振動が周波数ジャンプを引き起
こすのではなく、その中の一部が主振動と結合して温度
特性を悪化させることは知られていた。スプリアス振動
の主振動に与える程度を評価するために、表面に発生す
る電荷を用いることも提案された(信学技報 US97-
45(1997−09))。 水晶は圧電性を持つため、振
動により発生した結晶内のひずみは電荷に変換される。
その電荷量は式(5)で表わされる。
However, it has been empirically known that not all spurious vibrations calculated in the temperature characteristics cause frequency jumps, but some of them splice with the main vibrations to deteriorate the temperature characteristics. . In order to evaluate the degree to which the spurious vibration is applied to the main vibration, it has been proposed to use the electric charge generated on the surface.
45 (1997-09)). Since crystal has piezoelectricity, strain in the crystal generated by vibration is converted into electric charge.
The charge amount is represented by the equation (5).

【0018】[0018]

【数5】 [Equation 5]

【0019】FEM計算では、各スプリアス振動の変位
量が算出されるため、それらを電極付加された部分で積
分することにより、表面に発生する総電荷量も算出する
ことが出来る。しかし、スプリアス振動の表面電荷は、
他のスプリアス振動、又は主振動との相互作用により、
刻々と変化するものであるのに対して、従来の評価法は
連続的な動きを捉えるようなものではなかった。
Since the displacement amount of each spurious vibration is calculated in the FEM calculation, the total amount of electric charge generated on the surface can also be calculated by integrating the displacement amount of each spurious vibration. However, the surface charge of spurious vibration is
By interaction with other spurious vibrations or the main vibration,
While it changes moment by moment, the conventional evaluation method is not one that captures continuous movement.

【0020】そこで、本発明の目的は、スプリアス振動
と主振動との相互作用を考慮しつつ、スプリアス振動が
主振動に与える影響を評価することが可能な圧電振動子
の設計方法、モードチャートの作成方法、周波数温度特
性曲線の作成方法、圧電振動子の設計装置および圧電振
動子のシミュレーションプログラムを提供することであ
る。
Therefore, an object of the present invention is to design a piezoelectric vibrator capable of evaluating the influence of spurious vibration on the main vibration while considering the interaction between the spurious vibration and the main vibration. A method for creating a frequency-temperature characteristic curve, a device for designing a piezoelectric vibrator, and a simulation program for a piezoelectric vibrator are provided.

【0021】[0021]

【課題を解決するための手段】本発明では、スプリアス
振動の強弱を判断する指標として、式(5)を用いて算
出した各スプリアス振動の各寸法条件における総表面電
荷量を用いた。水晶は圧電結晶であるので、電荷量を多
く発生するということは、大きな変位で振動を起こして
いるということになり、当然、主振動付近に存在する時
には、厚み滑り振動を阻害する要因になると考えたから
である。具体的には、従来のモードチャートの主振動及
びスプリアス振動を表わす点の面積を、各振動が結合し
た際に持つ総電荷量に比例した値でプロットすることと
した。
In the present invention, the total surface charge amount under each dimensional condition of each spurious vibration calculated by using the equation (5) is used as an index for judging the strength of the spurious vibration. Since quartz is a piezoelectric crystal, generating a large amount of charge means vibrating due to a large displacement, and of course, when it exists near the main vibration, it becomes a factor that hinders thickness shear vibration. Because I thought. Specifically, the area of the points representing the main vibration and the spurious vibration of the conventional mode chart is plotted with a value proportional to the total amount of electric charge when each vibration is coupled.

【0022】[0022]

【発明の実施の形態】以下、本発明の実施形態に係る圧
電振動子の設計方法について、図面を参照しながら説明
する。 (実施例1)図5、6に示したモードチャートに代え
て、有限要素法により前記圧電振動子の主振動及びスプ
リアス振動の変位を求め、この変位により電極表面に発
生する電荷量を算出し、この電荷量を重み付けに採用し
た方式でプロットしたモードチャートを図7、8に示
す。なお、図7は、低周波モードチャートを示し、図8
は、高周波モードチャートを示す。また、図7、8中
で、白丸で表わしたのが主振動で、総電荷量が他のスプ
リアス振動に比べて非常に大きいため、電荷量の1/2
0の表示にしてある。その他の丸はスプリアス振動を表
わす。
BEST MODE FOR CARRYING OUT THE INVENTION A method of designing a piezoelectric vibrator according to an embodiment of the present invention will be described below with reference to the drawings. (Embodiment 1) Instead of the mode charts shown in FIGS. 5 and 6, displacements of the main vibration and spurious vibrations of the piezoelectric vibrator are obtained by the finite element method, and the amount of charge generated on the electrode surface by this displacement is calculated. 7 and 8 show mode charts plotted by a method in which this charge amount is used for weighting. 7 shows a low frequency mode chart, and FIG.
Shows a high frequency mode chart. The white circles in FIGS. 7 and 8 are the main vibrations, and since the total charge amount is much larger than other spurious vibrations, it is 1/2 of the charge amount.
It is displayed as 0. Other circles represent spurious vibrations.

【0023】図8から明らかになったことは、x3軸幅
滑り輪郭振動のように、主振動と結合して温度特性にお
いて周波数ジャンプを引き起こしやすいものは、図中矢
印で示されているように、スプリアス振動自体の電荷量
が多く、主振動に大きな影響を与えている。又、それ以
外のスプリアス振動についても、丸印で示したもののよ
うに、主振動から大きく離れた領域では、電荷量は大き
くないが、主振動に近づくと、主振動と結合して電荷量
が増えるものがある。
What has been clarified from FIG. 8 is that, as in the case of the x3 axis width sliding contour vibration, that which is likely to cause a frequency jump in the temperature characteristic in combination with the main vibration is indicated by an arrow in the figure. , The spurious vibration itself has a large amount of electric charge, which greatly affects the main vibration. Also, for other spurious vibrations, as shown by the circles, the charge amount is not large in the region greatly apart from the main vibration, but when it approaches the main vibration, the charge amount is combined with the main vibration and the charge amount becomes large. There is something to increase.

【0024】スプリアス振動の電荷量が増えるというこ
とは、主振動からスプリアス振動への電荷漏洩が起こっ
ていることを意味し、主振動の電荷量が減る。即ち、主
振動が阻害され、インピーダンス増加を起こしたり、周
波数ジャンプを起こしたりする。一方では、主振動に近
づいても、ほとんど、電荷量の変化を示さず、又、離れ
てしまうスプリアス振動も存在する。これらの現象を確
認するために、いくつかの実験を行なった。
The increase in the charge amount of the spurious vibration means that charge leakage from the main vibration to the spurious vibration occurs, and the charge amount of the main vibration decreases. That is, the main vibration is hindered, causing an impedance increase and a frequency jump. On the other hand, there are spurious vibrations that show almost no change in the amount of electric charge even when they approach the main vibration, and that they move away from each other. Several experiments were conducted to confirm these phenomena.

【0025】図8に示されたx3軸幅滑り輪郭振動に最
も影響されやすいと思われるA地点の幅(2c)寸法
(1296μm)において試作された振動子の温度特性を
図9に示す。図9(a)は主振動周波数偏差の温度依存
性、図9(b)に示すのは、図9(a)中に描かれた主
振動周波数偏差の3次曲線からのずれ量(カーブフィッ
トエラー)をプロットしたものである。図9(c)は、
インピーダンスの温度依存性を示す。図8のモードチャ
ートで推測できたように、高温側において大きな周波数
ジャンプが複数存在していることが分かる。この事か
ら、モードチャート上で電荷量が大きいスプリアス振動
の周辺では、実デバイスにおいても、スプリアス振動が
存在していることが確認できた。
Width (2c) dimension of point A that is most likely to be affected by the x3 axis width sliding contour vibration shown in FIG.
FIG. 9 shows the temperature characteristics of the vibrator prototyped at (1296 μm). FIG. 9A shows the temperature dependence of the main vibration frequency deviation, and FIG. 9B shows the deviation amount (curve fit) from the cubic curve of the main vibration frequency deviation depicted in FIG. 9A. Error) is plotted. FIG. 9C shows
The temperature dependence of impedance is shown. As can be inferred from the mode chart of FIG. 8, it can be seen that there are a plurality of large frequency jumps on the high temperature side. From this, it was confirmed on the mode chart that spurious vibrations exist even in the actual device around the spurious vibrations where the charge amount is large.

【0026】このように、主振動とスプリアス振動との
結合によって発生する総電荷量に注目し、主振動の電荷
量の減少に大きな影響を与えるスプリアス振動を相対的
に大きく表示し、主振動の電荷量の減少に与える影響の
小さなスプリアス振動を相対的に小さく表示する。これ
により、どのスプリアス振動が、主振動の温度特性を悪
化させるかを容易に判別することができ、周波数温度偏
差を考慮した圧電振動子の設計を容易に行うことが可能
となる。 (実施例2)図10は、同じA点における条件で、FE
M解析により算出された総電荷量を用いた温度特性を示
したものである。すなわち、スプリアス振動の周波数位
置のみでなく、有限要素法により圧電振動子の主振動及
びスプリアス振動の結合変位を求め、この結合変位によ
り電極表面に発生する電荷量を算出する方法を採用して
スプリアス振動に重み付けを行なうことにより、温度特
性曲線を作成した。
As described above, paying attention to the total charge amount generated by the coupling of the main vibration and the spurious vibration, the spurious vibration that has a large effect on the decrease in the charge amount of the main vibration is displayed relatively large, and the main vibration The spurious vibration, which has a small effect on the decrease in the amount of electric charge, is displayed relatively small. Accordingly, it is possible to easily determine which spurious vibration deteriorates the temperature characteristic of the main vibration, and it is possible to easily design the piezoelectric vibrator in consideration of the frequency temperature deviation. (Embodiment 2) FIG. 10 shows FE under the same condition at point A.
It shows the temperature characteristics using the total amount of charges calculated by M analysis. That is, not only the frequency position of the spurious vibration but also the coupling displacement of the main vibration and the spurious vibration of the piezoelectric vibrator is obtained by the finite element method, and the method of calculating the charge amount generated on the electrode surface due to this coupling displacement is adopted. The temperature characteristic curve was created by weighting the vibration.

【0027】実施例1におけるモードチャートと同様、
温度変化によるスプリアス振動の周波数偏差と主振動と
の結合の状態が、明確に確認できる。図9(a)〜
(c)の実験結果で確認した通り、高温領域において、
3つのスプリアス振動SP1、SP2、SP3が存在
し、特に、高温側の2つのスプリアス振動SP2、SP
3は、電荷量が大きく、主振動周波数の大きなジャンプ
を誘発することを示唆している。このように、温度特性
においても、電荷量を考慮することにより、より正確な
温度特性予測を行なうことが可能となる。 (実施例3)同様に、主振動近くに存在する電荷量が小
さいスプリアス振動についての現象を、実験により調べ
た。同じ条件で試作した振動子の温度特性測定結果を、
図12(a)〜(c)に示す。なお、図12(a)は、
B点の周波数温度特性、図12(b)は、B点のカーブフ
ィットエラー、図12(c)は、B点のインピーダンス
温度特性を示す。
Similar to the mode chart in the first embodiment,
The frequency deviation of spurious vibration due to temperature change and the state of coupling with the main vibration can be clearly confirmed. 9 (a)-
As confirmed by the experimental result of (c), in the high temperature region,
There are three spurious vibrations SP1, SP2, SP3, and in particular, two spurious vibrations SP2, SP on the high temperature side.
No. 3 suggests that the charge amount is large and induces a large jump of the main vibration frequency. As described above, even in the temperature characteristic, more accurate temperature characteristic prediction can be performed by considering the charge amount. (Example 3) Similarly, the phenomenon of spurious vibration having a small amount of charge existing near the main vibration was examined by an experiment. The temperature characteristic measurement results of the oscillator prototyped under the same conditions are
It is shown in FIGS. In addition, FIG.
12B shows the frequency-temperature characteristic at point B, FIG. 12B shows the curve-fit error at point B, and FIG. 12C shows the impedance-temperature characteristic at point B.

【0028】測定結果では、特に目立つスプリアス振動
を観測することは出来なかった。このように、電荷量の
評価で、主振動と大きな結合が無いと判断されたもの
は、実デバイスにおいても大きなスプリアス振動が観測
されず、電荷量を指標とすることは、有効な手段である
ことが確認できた。一方、図11は、計算より算出した
電荷量を考慮した温度特性を示す。図中に線で示したよ
うに、主振動に対して2つのスプリアス振動SP4、S
P5が交差しているが、総電荷量は大きく増加していな
い。 (実施例4)このように、電荷量を考慮に入れたモード
チャートを採用することは、主振動と他のスプリアス振
動との周波数位置関係および電荷の授受の程度を見るの
に有効であるが、電荷量の多いスプリアス振動の周辺は
隠れてしまい、細かい点まで見ることができないという
欠点もある。
From the measurement results, it was not possible to observe particularly noticeable spurious vibrations. In this way, in the evaluation of the charge amount, when it is determined that there is no large coupling with the main vibration, no large spurious vibration is observed even in the actual device, and using the charge amount as an index is an effective means. I was able to confirm that. On the other hand, FIG. 11 shows temperature characteristics in consideration of the calculated charge amount. As shown by the line in the figure, two spurious vibrations SP4 and S
Although P5 intersects, the total charge amount has not increased significantly. (Embodiment 4) As described above, adopting the mode chart in consideration of the amount of electric charge is effective in observing the frequency positional relationship between the main vibration and other spurious vibrations and the degree of charge transfer. However, there is a drawback that the area around the spurious vibration with a large amount of electric charge is hidden, and it is not possible to see even small points.

【0029】そこで、特に主振動に注目して、その電荷
量と周波数の変化を、図13に示すように、プロットし
てもよい。図13は、主振動と電荷量の幅(2c)寸法
依存性を示す。これは、図8と同様の条件で計算した結
果である。図13においては、電荷量と周波数の変化が
同時に表現されるため、スプリアス振動に主振動の電荷
が漏洩し、電荷量が減少した時に、その影響で周波数も
シフトしてしまう様子を良く観測することができる。こ
のグラフをもとに、周波数と電荷量の双方が安定した領
域を選ぶことができ、振動子の最適設計化が実現でき
る。
Therefore, paying particular attention to the main vibration, the change in the charge amount and the frequency may be plotted as shown in FIG. FIG. 13 shows the width (2c) size dependence of the main vibration and the amount of electric charge. This is the result of calculation under the same conditions as in FIG. In FIG. 13, since changes in the charge amount and the frequency are expressed at the same time, it is well observed that the charge of the main vibration leaks to the spurious vibration and the frequency shifts due to the influence of the decrease in the charge amount. be able to. Based on this graph, it is possible to select a region in which both frequency and charge amount are stable, and it is possible to realize optimum design of the vibrator.

【0030】図14は、本発明の一実施形態に係るAT
解析プログラムの動作を示すフローチャートである。図
14において、解析対象となるモデルを決定する(ステ
ップS1)。具体的には、ブランク形状、電極の有無、
電極や圧電体の材質、温度、境界条件などを決定する。
FIG. 14 shows an AT according to an embodiment of the present invention.
It is a flow chart which shows operation of an analysis program. In FIG. 14, a model to be analyzed is determined (step S1). Specifically, the blank shape, the presence or absence of electrodes,
The material, temperature, boundary conditions, etc. of the electrodes and piezoelectric body are determined.

【0031】次に、有限要素法における解析用メッシュ
を作成し(ステップS2)、固有値解析を行う(ステッ
プS3)。具体的には、メッシュ点を要素とする行列を
作成した後、固有値および固有ベクトルを算出し、共振
周波数を求める。次に、各共振周波数の解析を行い、変
位量、歪みおよび表面電荷を算出する(ステップS
4)。
Next, an analysis mesh in the finite element method is created (step S2), and eigenvalue analysis is performed (step S3). Specifically, after creating a matrix having mesh points as elements, eigenvalues and eigenvectors are calculated to obtain the resonance frequency. Next, each resonance frequency is analyzed to calculate the amount of displacement, strain, and surface charge (step S
4).

【0032】次に、表面電荷に対応したモードチャート
を作成するとともに(ステップS5)、温度特性曲線を
表示する(ステップS6)。このように、上述した実施
形態によれば、主振動により発生する電荷量およびスプ
リアス振動により発生する電荷量をそれぞれ別個に求め
るのではなく、主振動とスプリアス振動との結合により
発生する総電荷量を求める。
Next, a mode chart corresponding to the surface charge is created (step S5) and a temperature characteristic curve is displayed (step S6). As described above, according to the above-described embodiment, the total charge amount generated by the combination of the main vibration and the spurious vibrations is not obtained separately from the charge amount generated by the main vibrations and the charge amount generated by the spurious vibrations. Ask for.

【0033】これにより、様々のスプリアス振動のう
ち、主振動と結合して温度特性を悪化させるスプリアス
振動を特定したり、主振動と結合せず、温度特性を悪化
させることのないスプリアス振動を特定したりすること
が可能となる。このため、図6の高周波モードチャート
のように、全ての領域において、スプリアス振動で埋も
ってしまい、温度特性の良好な領域を見つけることが困
難な場合においても、主振動と結合しないスプリアス振
動を特定することにより、温度特性の良好な領域を見つ
けることが容易となり、余裕を持って設計を行うことが
可能となる。
As a result, of the various spurious vibrations, the spurious vibrations that are combined with the main vibrations to deteriorate the temperature characteristics are specified, and the spurious vibrations that are not combined with the main vibrations and do not deteriorate the temperature characteristics are specified. It becomes possible to do. Therefore, as shown in the high-frequency mode chart of FIG. 6, spurious vibrations that are buried in spurious vibrations in all regions and where it is difficult to find a region with good temperature characteristics identify spurious vibrations that are not coupled to the main vibration. By doing so, it becomes easy to find a region having a good temperature characteristic, and it becomes possible to design with a margin.

【0034】[0034]

【発明の効果】本発明は、上述したような方法を取るこ
とにより、以下に記載するような格別の効果を得ること
が出来る。本発明によれば、FEM解析により水晶振動
子の周波数特性中に存在する様々なスプリアス振動の変
位量と、それより算出されるスプリアス振動の電荷量を
算出し、その値をモードチャート上に同時にプロットす
ることにより、より正確にスプリアス振動の主振動に対
する影響を予測し、最適設計領域を正確に把握すること
が出来る。これにより、設計余裕度が増し、量産性が上
がるだけでなく、効率的な開発を行なうことが出来るた
め、試作回数の低減、開発期間の短縮を行なうことが出
来る。
According to the present invention, the following effects can be obtained by taking the method as described above. According to the present invention, displacement amounts of various spurious vibrations existing in the frequency characteristic of the crystal unit and the amount of charges of the spurious vibrations calculated from the displacement amounts are calculated by FEM analysis, and the values are simultaneously displayed on the mode chart. By plotting, the influence of the spurious vibration on the main vibration can be predicted more accurately, and the optimum design area can be accurately grasped. As a result, not only the design margin is increased and the mass productivity is improved, but also efficient development can be performed, so that the number of prototypes can be reduced and the development period can be shortened.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明における振動子の座標系と各寸法を表わ
す記号の定義を示す。
FIG. 1 shows a coordinate system of a vibrator according to the present invention and definitions of symbols representing respective dimensions.

【図2】主振動の振動モードを示す。FIG. 2 shows vibration modes of a main vibration.

【図3】(a)は一般的なATカット水晶振動子の周波
数温度特性を示し、(b)はスプリアス振動が存在する
ことにより周波数温度特性に発生する周波数ジャンプを
示す。
FIG. 3A shows a frequency-temperature characteristic of a general AT-cut crystal resonator, and FIG. 3B shows a frequency jump that occurs in the frequency-temperature characteristic due to the presence of spurious vibrations.

【図4】代表的なスプリアス振動の振動モードを示し、
(a)は厚み滑り副振動を示し、(b)はx3軸幅滑り
輪郭振動を示し、(c)はx1軸輪郭振動を示す。
FIG. 4 shows a vibration mode of a typical spurious vibration,
(A) shows thickness slip secondary vibration, (b) shows x3 axis width slip contour vibration, (c) shows x1 axis contour vibration.

【図5】従来のモードチャートで横軸は幅(2c)、縦
軸は周波数を表わす。
FIG. 5 is a conventional mode chart in which the horizontal axis represents width (2c) and the vertical axis represents frequency.

【図6】高周波帯における従来のモードチャートを示
す。
FIG. 6 shows a conventional mode chart in a high frequency band.

【図7】本発明の一実施形態による、各スプリアス振動
に対して電荷の考慮がされている低周波モードチャート
を示す。
FIG. 7 shows a low frequency mode chart with charge consideration for each spurious oscillation according to one embodiment of the invention.

【図8】本発明の一実施形態による、各スプリアス振動
に対して電荷の考慮がされている高周波モードチャート
を示す。
FIG. 8 shows a high frequency mode chart in which charge is taken into account for each spurious vibration according to one embodiment of the present invention.

【図9】(a)はA点における周波数温度特性を示し、
(b)はA点におけるカーブフィットエラーの温度特性
を示し、(c)はA点におけるインピーダンス温度特性
を示す。
FIG. 9A shows frequency-temperature characteristics at point A,
(B) shows the temperature characteristic of the curve fit error at the point A, and (c) shows the impedance temperature characteristic at the point A.

【図10】電荷を考慮した温度特性曲線を示す。FIG. 10 shows a temperature characteristic curve considering electric charges.

【図11】B点における温度特性 (算出値) を示す。FIG. 11 shows temperature characteristics (calculated values) at point B.

【図12】(a)はB点における周波数温度特性を示
し、(b)はB点におけるカーブフィットエラーの温度
を示し、(c)はB点におけるインピーダンス温度特性
を示す。
FIG. 12A shows frequency temperature characteristics at point B, FIG. 12B shows temperature of curve fit error at point B, and FIG. 12C shows impedance temperature characteristics at point B.

【図13】主振動と電荷量の幅(2c)寸法依存性を示
す。
FIG. 13 shows the width (2c) size dependence of the main vibration and the charge amount.

【図14】本発明の一実施形態に係るAT解析プログラ
ムの動作を示すフローチャートである。
FIG. 14 is a flowchart showing an operation of an AT analysis program according to an embodiment of the present invention.

───────────────────────────────────────────────────── フロントページの続き Fターム(参考) 5B046 AA07 BA01 JA02 JA04 5J108 AA01 AA04 BB02 CC04 DD02 MM11 (54)【発明の名称】 圧電振動子の設計方法、モードチャートの作成方法、周波数温度特性曲線の作成方法、圧電振動 子のシミュレーション方法、圧電振動子の設計装置及び圧電振動子のシミュレーションプログラ ム   ─────────────────────────────────────────────────── ─── Continued front page    F-term (reference) 5B046 AA07 BA01 JA02 JA04                 5J108 AA01 AA04 BB02 CC04 DD02                       MM11    (54) [Title of Invention] Piezoelectric vibrator design method, mode chart creation method, frequency temperature characteristic curve creation method, piezoelectric vibration                     Child simulation method, piezoelectric vibrator design device, and piezoelectric vibrator simulation program                     Mu

Claims (19)

【特許請求の範囲】[Claims] 【請求項1】 有限要素法により圧電振動子の主振動及
びスプリアス振動の変位を求め、この変位により電極表
面に発生する電荷量を算出することを特徴とする圧電振
動子の設計方法。
1. A method of designing a piezoelectric vibrator, which comprises obtaining a displacement of a main vibration and a spurious vibration of a piezoelectric vibrator by a finite element method, and calculating an electric charge amount generated on an electrode surface by this displacement.
【請求項2】 圧電振動子の設計をする際に、主振動及
びスプリアス振動の寸法依存性を表わすモードチャート
を作成し、そのモードチャートから最適設計値を決定す
る手法において、有限要素法により前記圧電振動子の主
振動及びスプリアス振動の変位を求め、この変位により
電極表面に発生する電荷量を算出し、この電荷量を重み
付けに採用したことを特徴とするモードチャートの作成
方法。
2. When designing a piezoelectric vibrator, in a method of creating a mode chart showing the dimensional dependence of main vibration and spurious vibration and determining an optimum design value from the mode chart, a finite element method is used. A method of creating a mode chart, characterized in that displacements of main vibration and spurious vibrations of a piezoelectric vibrator are obtained, an amount of electric charges generated on an electrode surface by the displacements is calculated, and the amount of electric charges is used for weighting.
【請求項3】 圧電振動子の設計をする際に、主振動及
びスプリアス振動の寸法依存性を表わすモードチャート
を作成し、そのモードチャートから最適設計値を決定す
る手法において、有限要素法により前記圧電振動子の主
振動及びスプリアス振動の変位を求め、この変位により
電極表面に発生する電荷量を算出し、この電荷量を重み
付けに採用したモードチャートから圧電振動子の最適設
計領域を決定することを特徴とする圧電振動子の設計方
法。
3. When designing a piezoelectric vibrator, in a method of creating a mode chart showing dimensional dependence of main vibration and spurious vibration, and determining an optimum design value from the mode chart, the finite element method is used. Obtain the displacement of the main vibration and spurious vibration of the piezoelectric vibrator, calculate the amount of charge generated on the electrode surface due to this displacement, and determine the optimum design area of the piezoelectric vibrator from the mode chart that uses this amount of charge for weighting. A method for designing a piezoelectric vibrator.
【請求項4】 有限要素法により算出した圧電振動子の
周波数温度特性曲線の作成方法において、スプリアス振
動の周波数位置のみでなく、有限要素法により前記圧電
振動子の主振動及びスプリアス振動の変位を求め、この
変位により電極表面に発生する電荷量を算出する方法を
採用してスプリアス振動に重み付けを行なうことを特徴
とする周波数温度特性曲線の作成方法。
4. A method of creating a frequency-temperature characteristic curve of a piezoelectric vibrator calculated by the finite element method, wherein not only the frequency position of the spurious vibration but also the displacement of the main vibration and the spurious vibration of the piezoelectric vibrator is calculated by the finite element method. A method of creating a frequency-temperature characteristic curve, characterized in that a method of calculating the amount of electric charge generated on the electrode surface by this displacement is employed to weight the spurious vibration.
【請求項5】 有限要素法により算出した圧電振動子の
周波数温度特性曲線を用いた圧電振動子の設計方法にお
いて、スプリアス振動の周波数位置のみでなく、有限要
素法により前記圧電振動子の主振動及びスプリアス振動
の変位を求め、この変位により電極表面に発生する電荷
量を算出する方法を採用してスプリアス振動に重み付け
を行なうことにより作成した温度特性曲線を用いて、設
計最適領域を決定することを特徴とする圧電振動子の設
計方法。
5. A method for designing a piezoelectric vibrator using a frequency-temperature characteristic curve of the piezoelectric vibrator calculated by the finite element method, wherein not only the frequency position of spurious vibration but also the main vibration of the piezoelectric vibrator by the finite element method is used. And determine the optimum design area by using the temperature characteristic curve created by weighting the spurious vibrations by calculating the displacement of spurious vibrations and calculating the amount of charge generated on the electrode surface by this displacement. A method for designing a piezoelectric vibrator.
【請求項6】 有限要素法により圧電振動子の主振動及
びスプリアス振動の変位を求め、この変位により電極表
面に発生する電荷量を算出する方法を用いて主振動の電
荷量を算出し、その形状依存性から最適設計形状を決定
することを特徴とする圧電振動子の設計方法。
6. A charge amount of the main vibration is calculated by a method of calculating displacements of the main vibration and spurious vibrations of the piezoelectric vibrator by the finite element method, and calculating a charge amount generated on the electrode surface by this displacement. A method for designing a piezoelectric vibrator, characterized in that an optimum design shape is determined from shape dependence.
【請求項7】 請求項6の圧電振動子の設計方法におい
て、前記電荷量に加えて、周波数の形状依存性も考慮に
入れ、前記電荷量及び前記周波数ともに安定した領域を
最適設計領域とすることを特徴とする圧電振動子の設計
方法。
7. The method of designing a piezoelectric vibrator according to claim 6, wherein, in addition to the charge amount, shape dependency of frequency is also taken into consideration, and a region in which both the charge amount and the frequency are stable is set as an optimum design region. A method for designing a piezoelectric vibrator, comprising:
【請求項8】 請求項1、3、5、6又は7において、
前記圧電振動子が水晶振動子であることを特徴とする圧
電振動子の設計方法。
8. The method according to claim 1, 3, 5, 6, or 7,
A method for designing a piezoelectric vibrator, wherein the piezoelectric vibrator is a crystal vibrator.
【請求項9】 請求項2において、前記圧電振動子が水
晶振動子であることを特徴とするモードチャートの作成
方法。
9. The method for producing a mode chart according to claim 2, wherein the piezoelectric vibrator is a crystal vibrator.
【請求項10】 請求項4において、前記圧電振動子が
水晶振動子であることを特徴とする周波数温度特性曲線
の作成方法。
10. The method according to claim 4, wherein the piezoelectric vibrator is a crystal vibrator.
【請求項11】 主振動とスプリアス振動との結合によ
って発生する総電荷量に基づいて、圧電振動子のブラン
ク寸法を決定することを特徴とする圧電振動子の設計方
法。
11. A method for designing a piezoelectric vibrator, characterized in that a blank size of the piezoelectric vibrator is determined based on a total amount of charges generated by coupling of main vibration and spurious vibration.
【請求項12】 前記ブランク寸法は、前記スプリアス
振動による前記主振動の電荷量の減少が所定値以下とな
るように決定することを特徴とする請求項11記載の圧
電振動子の設計方法。
12. The method of designing a piezoelectric vibrator according to claim 11, wherein the blank size is determined so that a decrease in the charge amount of the main vibration due to the spurious vibration is equal to or less than a predetermined value.
【請求項13】 主振動とスプリアス振動との結合によ
って発生する総電荷量を算出することを特徴とする圧電
振動子のシミュレーション方法。
13. A method of simulating a piezoelectric vibrator, which comprises calculating a total amount of charges generated by coupling of a main vibration and a spurious vibration.
【請求項14】 圧電振動子のブランク寸法をパラメー
タとして入力するパラメータ入力手段と、 前記ブランク寸法について、主振動とスプリアス振動と
の結合によって発生する総電荷量を算出する総電荷量算
出手段とを備えることを特徴とする圧電振動子の設計装
置。
14. A parameter input means for inputting a blank size of a piezoelectric vibrator as a parameter, and a total charge amount calculating means for calculating a total charge amount generated by coupling of a main vibration and a spurious vibration with respect to the blank size. A device for designing a piezoelectric vibrator, comprising:
【請求項15】 前記スプリアス振動を示す点を、前記
総電荷量に対応してモードチャート上に拡大表示するモ
ードチャート表示手段をさらに備えることを特徴とする
請求項14記載の圧電振動子の設計装置。
15. The design of the piezoelectric vibrator according to claim 14, further comprising a mode chart display means for enlarging and displaying a point indicating the spurious vibration on a mode chart corresponding to the total charge amount. apparatus.
【請求項16】 前記モードチャート表示手段は、前記
スプリアス振動を示す点を、前記総電荷量に対応する面
積を有する円で表示することを特徴とする請求項15記
載の圧電振動子の設計装置。
16. The device for designing a piezoelectric vibrator according to claim 15, wherein the mode chart display means displays a point indicating the spurious vibration with a circle having an area corresponding to the total charge amount. .
【請求項17】 前記ブランク寸法に対する前記主振動
の電荷量および周波数の変化を表示する主振動表示手段
をさらに備えることを特徴とする請求項14または15
記載の圧電振動子の設計装置。
17. The main vibration display means for displaying changes in the charge amount and frequency of the main vibration with respect to the blank size is further provided.
The piezoelectric vibrator design apparatus described.
【請求項18】 圧電振動子のブランク寸法をパラメー
タとして、主振動とスプリアス振動との結合によって発
生する電荷量を算出するステップをコンピュータに実行
させることを特徴とする圧電振動子のシミュレーション
プログラム。
18. A computer program for simulating a piezoelectric vibrator, which causes a computer to execute a step of calculating a charge amount generated by coupling of a main vibration and a spurious vibration using a blank size of the piezoelectric vibrator as a parameter.
【請求項19】 指定されたブランク寸法に対し、温度
をパラメータとして、主振動とスプリアス振動との結合
によって発生する電荷量を算出するステップをコンピュ
ータに実行させることを特徴とする圧電振動子のシミュ
レーションプログラム。
19. A simulation of a piezoelectric vibrator, comprising causing a computer to execute a step of calculating an amount of electric charge generated by coupling of a main vibration and a spurious vibration with a temperature as a parameter for a specified blank dimension. program.
JP2001226772A 2001-05-14 2001-07-26 Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator Withdrawn JP2003037462A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001226772A JP2003037462A (en) 2001-05-14 2001-07-26 Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2001-143839 2001-05-14
JP2001143839 2001-05-14
JP2001226772A JP2003037462A (en) 2001-05-14 2001-07-26 Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator

Publications (1)

Publication Number Publication Date
JP2003037462A true JP2003037462A (en) 2003-02-07

Family

ID=26615069

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001226772A Withdrawn JP2003037462A (en) 2001-05-14 2001-07-26 Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator

Country Status (1)

Country Link
JP (1) JP2003037462A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7157983B2 (en) 2003-03-19 2007-01-02 Seiko Epson Corporation Cut of piezoelectric oscillator, piezoelectric oscillator, and piezoelectric device

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7157983B2 (en) 2003-03-19 2007-01-02 Seiko Epson Corporation Cut of piezoelectric oscillator, piezoelectric oscillator, and piezoelectric device

Similar Documents

Publication Publication Date Title
JP4306668B2 (en) Lamb wave type high frequency resonator
US7878044B2 (en) Sensor, system, and method, for measuring fluid properties using multi-mode quasi-shear-horizontal resonator
CN101682309B (en) MEMS resonators
Heikkola et al. Multiobjective optimization of an ultrasonic transducer using NIMBUS
Kang et al. A high-sensitivity pressure sensor based on surface transverse wave
Hammad Natural frequencies and mode shapes of mechanically coupled microbeam resonators with an application to micromechanical filters
JP3885893B2 (en) Method for generating a tuning fork member (DETF) having two ends based on a suspect mode map
Lin et al. Vibration modeling of disk resonator gyroscope by wave propagation
JP2003037462A (en) Method for designing piezoelectric oscillator, method for preparing mode chart, method for preparing frequency/temperature characteristics curve, simulation method of the piezoelectric oscillator, designing device of the piezoelectric oscillator and simulation program of the piezoelectric oscillator
Zhu et al. A new model for acoustic attenuation of GHz waveguide induced by metal electrodes in piezoelectric composites
JPWO2008114715A1 (en) Lamb wave type elastic wave device
Hao et al. Vibration displacement on substrate due to time-harmonic stress sources from a micromechanical resonator
JP4465464B2 (en) Lamb wave type elastic wave device
US6654711B1 (en) Correction factors for the analysis of piezoelectric devices
Wakatsuki et al. Temperature‐frequency characteristics simulation of piezoelectric resonators and their equivalent circuits based on three‐dimensional finite element modelling
Lin Coupled vibration and natural frequency analysis of isotropic cylinders or disks of finite dimensions
Lu et al. Energy-trapping analysis for the bi-stepped mesa quartz crystal microbalance using the finite element method
US6658376B1 (en) Determination of vibration frequency characteristics of electroded crystal plate
Kidambi Frequency and damping of non-axisymmetric surface oscillations of a viscous cylindrical liquid bridge
Sekimoto et al. Frequency-temperature behavior of spurious vibrations of rectangular AT-cut quartz plates
Lindenbauer et al. Fully three-dimensional analysis of tsm quartz sensors immersed in viscous liquids
CN116429281B (en) Resonator based on array structure and temperature measurement method
JP7218839B1 (en) Elastic matrix determination method and vibration analysis method for laminated core of transformer
Gufflet et al. Isochronism defect for various doubly rotated cut quartz resonators
Li et al. A Comsol analysis of vibrations of quartz crystal plates with Mindlin plate equations

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040323

A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20040430