JP2002267960A - Method and device for vibration analysis of optical equipment - Google Patents

Method and device for vibration analysis of optical equipment

Info

Publication number
JP2002267960A
JP2002267960A JP2001071091A JP2001071091A JP2002267960A JP 2002267960 A JP2002267960 A JP 2002267960A JP 2001071091 A JP2001071091 A JP 2001071091A JP 2001071091 A JP2001071091 A JP 2001071091A JP 2002267960 A JP2002267960 A JP 2002267960A
Authority
JP
Japan
Prior art keywords
vibration
optical
calculation
polynomial
finite element
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.)
Granted
Application number
JP2001071091A
Other languages
Japanese (ja)
Other versions
JP3595775B2 (en
Inventor
Tatsuro Kawakami
竜郎 川上
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP2001071091A priority Critical patent/JP3595775B2/en
Publication of JP2002267960A publication Critical patent/JP2002267960A/en
Application granted granted Critical
Publication of JP3595775B2 publication Critical patent/JP3595775B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Mechanical Optical Scanning Systems (AREA)
  • Lenses (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a vibration analyzing method of optical equipment by which the influence of vibration in the actual operation of an optical equipment given to the performance of the optical equipment can accurately be known. SOLUTION: The vibration analyzing method of optical equipment for analyzing the influence of the vibration of an optical component constituting the optical equipment given to the performance of the optical equipment, includes a vibration measuring process (S1) of actually measuring the vibration frequency and amplitude of the optical component, a numerical model generating process (S2) of generating a simulation model of a finite element method according to information on the measured vibration frequency and amplitude of the optical component, a vibration computing process (S3) of computing the vibration by using the simulation model generated by the numerical model generating process, a curved surface fitting process of approximating the deformation state due to the vibration of the optical component computed by the vibration computing process with a shape function capable of unique shape representation in a three-dimensional space, and an optical computing process (S4) of computing characteristics related to the performance of the optical equipment according to the shape function.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、複写機やレーザー
ビームプリンタなどの電子写真方式の画像形成装置に装
備される光学機器、またはカメラのオートフォーカス機
能を司る光学機器の設計技術に係り、光学機器の実使用
状況下で生じ得る光学部品の変形が機器の性能劣化に与
える影響を詳細に解析することで、設計改善の方向性を
明瞭に見出すようにした光学機器の振動解析方法及び装
置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a technology for designing an optical device mounted on an electrophotographic image forming apparatus such as a copying machine or a laser beam printer, or an optical device for controlling an autofocus function of a camera. The present invention relates to a method and an apparatus for analyzing vibration of an optical device, in which the influence of the deformation of an optical component that may occur under the actual use condition of the device on the performance degradation of the device is analyzed in detail, so that the direction of the design improvement is clearly found. Things.

【0002】[0002]

【従来の技術】複写機やLBP(レーザービームプリン
ター)等に装着される光学ユニットの設計は、各ユニッ
トを静止した状態とし、3次元空間における配置を光学
経路の計算により決定していた。光学経路の計算は、各
ユニットの光の透過部、反射部の形状を空間変数によっ
て関数表現し、光の直進、反射、屈折といった諸原理を
用いて行っていた。
2. Description of the Related Art The design of an optical unit to be mounted on a copying machine, an LBP (laser beam printer), or the like is such that each unit is kept stationary and the arrangement in a three-dimensional space is determined by calculating an optical path. The calculation of the optical path has been performed by expressing the shapes of the light transmitting portion and the reflecting portion of each unit as a function by a spatial variable and using various principles such as straight traveling, reflection, and refraction of light.

【0003】複写機やLBP等の製品の実稼動時には、
各光学ユニットが振動し、静止時を前提にして設計され
た光学経路が変化してしまう。光学ユニットの振動によ
る光学経路の変化については、特開平11−11913
6号公報で光学ユニットを有限要素モデルで表現し、振
動状態のシミュレーションによって変化量を計算する手
法が提案されている。
[0003] During the actual operation of products such as copiers and LBP,
Each optical unit vibrates, and an optical path designed on the assumption that the optical unit is stationary changes. The change of the optical path due to the vibration of the optical unit is described in JP-A-11-11913.
Japanese Patent Laid-Open No. 6-2006 proposes a method of expressing an optical unit by a finite element model and calculating a change amount by simulation of a vibration state.

【0004】また圧電型加速度センサおよびFFTアナ
ライザにより、光学部品の振動挙動を実測により求める
技術も広く普及している。
[0004] Techniques for actually measuring the vibration behavior of an optical component using a piezoelectric acceleration sensor and an FFT analyzer have also become widespread.

【0005】[0005]

【発明が解決しようとする課題】近年プリンターの高速
化、高画質化に伴い、実稼動時に発生する光学ユニット
の振動や機内の昇温によって発生する光学ユニットの熱
歪など、光学ユニットの変形が画像に与える影響を無視
できなくなってきている。
In recent years, with the increase in speed and image quality of printers, deformation of optical units, such as vibration of optical units during actual operation and thermal distortion of optical units caused by temperature rise inside the machine, has been increasing. The effect on the image is no longer negligible.

【0006】より高画質のプリンターを提供するために
は、各光学ユニットが静止した理想的な状況の光学経路
の計算に加え、光学ユニットに発生する変形が光学経路
に与える影響を正確に考慮して、各光学ユニットの配置
を決定する必要がある。
In order to provide a printer with higher image quality, in addition to calculating the optical path in an ideal situation where each optical unit is stationary, the influence of the deformation generated in the optical unit on the optical path must be accurately considered. Therefore, it is necessary to determine the arrangement of each optical unit.

【0007】しかしながら上記従来例(特開平11−1
19136号公報)では、有限要素法による振動シミュ
レーションのみでは減衰効果で決定される各光学ユニッ
トの振幅値が容易に予測できず、正確に光学経路を計算
できる系の範囲が限定されてしまう欠点があった。
However, the above conventional example (Japanese Patent Laid-Open No. 11-1)
19136) has a drawback that the amplitude value of each optical unit determined by the damping effect cannot be easily predicted only by the vibration simulation by the finite element method, and the range of the system in which the optical path can be calculated accurately is limited. there were.

【0008】一方、圧電型加速度センサおよびFFTア
ナライザを用い、機器の振動状態を測定することは可能
であるが、振動時の光学経路を求める上では測定点数が
機器内へのセンサの配線等の物理的な問題で不十分であ
った。また光学ユニットの回転変動に伴う光学経路を測
定するには、本来並進方向の振幅しか測れない圧電型加
速度センサ2点を用い、2点での応答および2点間の距
離から回転角を測定することができる。しかし、測定す
る光学部品のサイズがセンサのサイズに対し相対的に小
さくなるほど、正確な回転角を測定することが困難にな
り、光学ユニットの微量な回転変動を測定できないとい
う欠点があった。
On the other hand, it is possible to measure the vibration state of a device by using a piezoelectric acceleration sensor and an FFT analyzer. However, when obtaining an optical path during vibration, the number of measurement points is limited by the wiring of the sensor inside the device. Physical problems were inadequate. In addition, to measure an optical path due to rotation fluctuation of the optical unit, two rotation points of the piezoelectric acceleration sensor which can measure only the amplitude in the translation direction are used, and the rotation angle is measured from the response at the two points and the distance between the two points. be able to. However, as the size of the optical component to be measured becomes smaller relative to the size of the sensor, it becomes more difficult to accurately measure the rotation angle, and there is a disadvantage that a minute rotation fluctuation of the optical unit cannot be measured.

【0009】従って、本発明は上述した課題に鑑みてな
されたものであり、その目的は、光学機器の実稼動時の
振動が光学機器の性能に与える影響を正確に知ることが
できる光学機器の振動解析方法及び装置を提供すること
である。
SUMMARY OF THE INVENTION Accordingly, the present invention has been made in view of the above-mentioned problems, and an object of the present invention is to provide an optical device capable of accurately knowing the influence of the vibration during the actual operation of the optical device on the performance of the optical device. An object of the present invention is to provide a vibration analysis method and apparatus.

【0010】[0010]

【課題を解決するための手段】上述した課題を解決し、
目的を達成するために、本発明に係わる光学機器の振動
解析方法は、光学機器の実稼動状況下における、該光学
機器を構成する光学部品の振動が前記光学機器の性能に
与える影響を解析するための光学機器の振動解析方法で
あって、前記光学部品の振動数及び振幅を実際に測定す
る振動測定工程と、 測定された前記光学部品の振動数
及び振幅の情報に基づいて有限要素法のシミュレーショ
ンモデルを作成する数値モデル作成工程と、該数値モデ
ル作成工程において生成されたシミュレーションモデル
を用いて振動計算を行う振動計算工程と、該振動計算工
程において計算された前記光学部品の振動による変形状
態を、3次元空間に対する一意の形状表現が可能な形状
関数で近似する曲面フィッティング工程と、前記形状関
数に基づいて前記光学機器の性能に関係する特性の計算
を行う光学計算工程と、を具備することを特徴としてい
る。
Means for Solving the Problems The above-mentioned problems are solved,
In order to achieve the object, a vibration analysis method for an optical device according to the present invention analyzes an influence of vibration of an optical component included in the optical device on the performance of the optical device under actual operating conditions of the optical device. A vibration analysis method for actually measuring the frequency and amplitude of the optical component; and a finite element method based on the measured information on the frequency and amplitude of the optical component. A numerical model generating step of generating a simulation model, a vibration calculating step of performing vibration calculation using the simulation model generated in the numerical model forming step, and a deformation state of the optical component calculated by the vibration calculated in the vibration calculating step Is approximated by a shape function capable of expressing a unique shape in a three-dimensional space; and And an optical calculation step of calculating characteristics related to the performance of the scientific instrument.

【0011】また、この発明に係わる光学機器の振動解
析方法において、前記数値モデル作成工程は、光学経路
の許容誤差を指定する許容誤差指定工程と、前記シミュ
レーションモデル中の有限要素を選択する要素選択工程
と、選択された有限要素の材料特性及び/又は形状特性
を変化させて振動計算を行い、前記振動測定工程で測定
された振動数及び振幅に対する前記有限要素の材料特性
及び/又は形状特性の感度を調べる感度計算工程と、計
算された感度に基づいて、前記許容誤差内に収まるま
で、有限要素の材料特性及び/又は形状特性を変更して
収束計算を行うモデル最適化工程とを備えることを特徴
としている。
Further, in the vibration analysis method for an optical device according to the present invention, the numerical model creating step includes an allowable error specifying step for specifying an allowable error of an optical path, and an element selection for selecting a finite element in the simulation model. A vibration calculation by changing the material property and / or shape property of the selected finite element, and calculating the material property and / or shape property of the finite element with respect to the frequency and amplitude measured in the vibration measurement step. A sensitivity calculation step of examining the sensitivity, and a model optimization step of performing a convergence calculation by changing a material property and / or a shape property of a finite element based on the calculated sensitivity until the calculated value falls within the allowable error. It is characterized by.

【0012】また、この発明に係わる光学機器の振動解
析方法において、前記曲面フィッティング工程では、前
記形状関数を空間変数x,y,zの多項式で表現し、前
記シミュレーションモデルによって離散的に計算された
変形状態から最小2乗法によって前記多項式の各係数を
求めて曲面フィッティングを行うことを特徴としてい
る。
In the method of analyzing vibration of an optical instrument according to the present invention, in the curved surface fitting step, the shape function is represented by a polynomial of spatial variables x, y, and z, and is discretely calculated by the simulation model. It is characterized in that the surface fitting is performed by obtaining each coefficient of the polynomial from the deformed state by the least square method.

【0013】また、この発明に係わる光学機器の振動解
析方法において、前記曲面フィッティング工程は、前記
多項式による曲面の近似精度の許容値を指定する許容値
指定工程と、前記多項式の次数を指定する次数指定工程
と、前記多項式によって計算された光学部品の変形量と
測定された変形量とを比較し、その差が前記許容値内で
あれば前記次数指定工程において定義された多項式を採
用し、許容値外であれば許容値内に収まるまで前記多項
式の次数を増加させて計算する収束計算工程とを備える
ことを特徴としている。
In the vibration analysis method for an optical instrument according to the present invention, the curved surface fitting step includes an allowable value specifying step of specifying an allowable value of approximation accuracy of the curved surface by the polynomial, and an order of specifying an order of the polynomial. Specifying step, comparing the deformation amount of the optical component calculated by the polynomial and the measured deformation amount, if the difference is within the allowable value, adopts the polynomial defined in the order specifying step, the allowable A convergence calculation step of increasing the degree of the polynomial and calculating the degree until the value falls outside the allowable value.

【0014】また、本発明に係わる光学機器の振動解析
装置は、光学機器の実稼動状況下における、該光学機器
を構成する光学部品の振動が前記光学機器の性能に与え
る影響を解析するための光学機器の振動解析装置であっ
て、前記光学部品の振動数及び振幅を実際に測定する振
動測定手段と、測定された前記光学部品の振動数及び振
幅の情報に基づいて有限要素法のシミュレーションモデ
ルを作成する数値モデル作成手段と、該数値モデル作成
手段により生成されたシミュレーションモデルを用いて
振動計算を行う振動計算手段と、該振動計算手段により
計算された前記光学部品の振動による変形状態を、3次
元空間に対する一意の形状表現が可能な形状関数で近似
する曲面フィッティング手段と、前記形状関数に基づい
て前記光学機器の性能に関係する特性の計算を行う光学
計算手段と、を具備することを特徴としている。
Further, the vibration analyzing apparatus for an optical device according to the present invention is provided for analyzing the influence of the vibration of an optical component constituting the optical device on the performance of the optical device under actual operating conditions of the optical device. An apparatus for analyzing vibration of an optical device, comprising: a vibration measuring unit for actually measuring the frequency and amplitude of the optical component; and a simulation model of a finite element method based on the measured information on the frequency and amplitude of the optical component. Numerical model creating means for creating, a vibration calculating means for performing a vibration calculation using a simulation model generated by the numerical model creating means, and the deformation state of the optical component due to vibration calculated by the vibration calculating means, A curved surface fitting unit that approximates with a shape function capable of expressing a unique shape in a three-dimensional space; It is characterized by including the optical calculation means for performing calculation of characteristics relating to performance, the.

【0015】また、この発明に係わる光学機器の振動解
析装置において、前記数値モデル作成手段は、光学経路
の許容誤差を指定する許容誤差指定手段と、前記シミュ
レーションモデル中の有限要素を選択する要素選択手段
と、選択された有限要素の材料特性及び/又は形状特性
を変化させて振動計算を行い、前記振動測定手段により
測定された振動数及び振幅に対する前記有限要素の材料
特性及び/又は形状特性の感度を調べる感度計算手段
と、計算された感度に基づいて、前記許容誤差内に収ま
るまで、有限要素の材料特性及び/又は形状特性を変更
して収束計算を行うモデル最適化手段とを備えることを
特徴としている。
Further, in the vibration analyzing apparatus for optical equipment according to the present invention, the numerical model creating means includes an allowable error specifying means for specifying an allowable error of an optical path, and an element selecting means for selecting a finite element in the simulation model. Means for calculating a vibration by changing a material property and / or a shape property of the selected finite element, and calculating a material property and / or a shape property of the finite element with respect to a frequency and an amplitude measured by the vibration measuring means. A sensitivity calculating means for examining the sensitivity, and a model optimizing means for performing a convergence calculation by changing a material property and / or a shape property of a finite element based on the calculated sensitivity until the calculated value falls within the allowable error. It is characterized by.

【0016】また、この発明に係わる光学機器の振動解
析装置において、前記曲面フィッティング手段は、前記
形状関数を空間変数x,y,zの多項式で表現し、前記
シミュレーションモデルによって離散的に計算された変
形状態から最小2乗法によって前記多項式の各係数を求
めて曲面フィッティングを行うことを特徴としている。
Further, in the vibration analyzing apparatus for optical equipment according to the present invention, the curved surface fitting means expresses the shape function by a polynomial of spatial variables x, y, z, and is discretely calculated by the simulation model. It is characterized in that the surface fitting is performed by obtaining each coefficient of the polynomial from the deformed state by the least square method.

【0017】また、この発明に係わる光学機器の振動解
析装置において、前記曲面フィッティング手段は、前記
多項式による曲面の近似精度の許容値を指定する許容値
指定手段と、前記多項式の次数を指定する次数指定手段
と、前記多項式によって計算された光学部品の変形量と
測定された変形量とを比較し、その差が前記許容値内で
あれば前記次数指定手段により定義された多項式を採用
し、許容値外であれば許容値内に収まるまで前記多項式
の次数を増加させて計算する収束計算手段とを備えるこ
とを特徴としている。
Further, in the vibration analyzing apparatus for optical equipment according to the present invention, the curved surface fitting means includes an allowable value specifying means for specifying an allowable value of approximation accuracy of the curved surface by the polynomial, and an order for specifying an order of the polynomial. Designating means, comparing the deformation of the optical component calculated by the polynomial and the measured deformation, and if the difference is within the allowable value, adopts the polynomial defined by the degree specifying means, A convergence calculation means for calculating by increasing the degree of the polynomial until the value falls outside the allowable value.

【0018】[0018]

【発明の実施の形態】以下、本発明の好適な一実施形態
について説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, a preferred embodiment of the present invention will be described.

【0019】まず、一実施形態の概要について説明す
る。
First, an outline of an embodiment will be described.

【0020】一実施形態の光学機器の設計方法は、光学
機器の実使用状況下において、機器を構成する光学部品
の振動が光学機器性能に与える影響を解析する方法であ
って、光学部品の振動数、振幅を実際に測定する振動測
定工程と、測定された前記光学部品の振動数、振幅情報
から有限要素法のシミュレーションモデルを作成する数
値モデル作成工程と、前記数値モデル作成工程において
生成されたシミュレーションモデルを用い、振動計算を
行う振動計算工程と、前記振動計算工程において計算さ
れた光学部品の振動(変形)状態を、3次元空間に対す
る一意の形状表現が可能な形状関数に変換する曲面フィ
ッティング工程と、前記形状関数に基づいて光路軌跡等
の計算を行う光学計算工程とを有している。
The method for designing an optical device according to one embodiment is a method for analyzing the influence of the vibration of an optical component constituting the device on the performance of the optical device under actual use conditions of the optical device. A vibration measurement step of actually measuring the number and amplitude, a numerical model creation step of creating a simulation model of the finite element method from the measured vibration frequency and amplitude information of the optical component, and a numerical model creation step. A vibration calculation step of performing a vibration calculation using a simulation model; and a curved surface fitting for converting a vibration (deformation) state of the optical component calculated in the vibration calculation step into a shape function capable of expressing a unique shape in a three-dimensional space. And an optical calculation step of calculating an optical path trajectory and the like based on the shape function.

【0021】この構成において、振動測定工程では、有
限要素法等の数値シミュレーションのみでは予測しきれ
ない各光学ユニットの振動数毎の減衰効果を数値モデル
作成工程に渡し、数値モデル作成工程では、各光学ユニ
ットの詳細形状を有限要素法モデルによって定義するこ
とで、測定では機器内へのセンサの配線等の物理的な問
題で不十分となる測定点を補い、振動計算工程では、実
際の振動測定では測定困難な光学ユニットの振動成分
(特に回転変動)を正確に計算し、曲面フィッティング
工程では、前記振動計算工程において計算された光学面
の振動状態が有限要素法モデルの特性上、節点上での空
間的な離散点データとして表現されるので、任意の位置
での光学計算を行うために、光学面の離散点での振動状
態をもとに3次元空間に対する一意の形状表現が可能な
形状関数に変換して光学計算工程に渡し、光学計算工程
では、各光学面の振動モードを考慮した光学経路の計算
を行うことが可能となるため、実稼動時に発生する光学
ユニットの振動変形が光学経路に与える影響を正確に考
慮して、各光学ユニットの配置を決定することを可能と
する。
In this configuration, in the vibration measurement step, the damping effect for each frequency of each optical unit, which cannot be predicted only by numerical simulation such as the finite element method, is passed to the numerical model creation step. By defining the detailed shape of the optical unit using a finite element method model, it compensates for measurement points that would be insufficient due to physical problems such as sensor wiring in the equipment during measurement, and the actual vibration measurement in the vibration calculation process Then, the vibration component (especially rotation fluctuation) of the optical unit which is difficult to measure is accurately calculated, and in the curved surface fitting step, the vibration state of the optical surface calculated in the vibration calculation step is determined by the characteristics of the finite element method model, Is represented as spatial discrete point data, so that optical calculations at arbitrary positions can be performed based on the vibrational state at discrete points on the optical surface. Is converted to a shape function capable of expressing a unique shape for the optical calculation process and passed to the optical calculation process.In the optical calculation process, it is possible to calculate the optical path in consideration of the vibration mode of each optical surface. The arrangement of each optical unit can be determined by accurately considering the influence of the generated vibration deformation of the optical unit on the optical path.

【0022】また、上記の数値モデル作成工程におい
て、光学経路の許容誤差を指定する許容誤差指定工程
と、モデル中の有限要素を選択する要素選択工程と、選
択された有限要素の材料特性、形状特性を変化させ振動
計算を行い、測定された振動数、振幅に対する有限要素
の材料特性、形状特性の感度を調べる感度計算工程と、
前記感度に基づいて、光学経路の変動量が前記許容誤差
内に収まるまで、有限要素の材料特性、形状特性を変更
して収束計算を行う、モデル最適化工程とを有してい
る。
Further, in the above numerical model creation step, an allowable error specifying step for specifying an allowable error of the optical path, an element selecting step for selecting a finite element in the model, a material property and a shape of the selected finite element A sensitivity calculation step of performing a vibration calculation by changing the characteristics, and examining the sensitivity of the material characteristics of the finite element to the measured frequency and amplitude, and the sensitivity of the shape characteristics,
A model optimization step of performing a convergence calculation by changing the material characteristics and shape characteristics of the finite element until the variation amount of the optical path falls within the allowable error based on the sensitivity.

【0023】この構成により、製品の実使用時に発生す
る振動等により、光学経路の変動量が問題になる場合、
許容誤差指定工程において、光学経路の許容量を指定
し、要素選択工程において、光学経路を決定しうる設計
変数を有限要素モデルの中で定義し、感度計算工程にお
いて、前記要素選択工程で定義された設計変数の光学経
路の変動に対する感度係数を求め、モデル最適化工程に
おいて、前記許容誤差指定工程で指定された光学経路の
変動量が許容値内に収まるまで、前記感度計算工程で抽
出された設計変数の内、感度係数が高いものを優先的に
操作し、繰り返し振動計算を行い、計算機によって最適
な設計変数の決定を可能とする。
According to this configuration, when the amount of fluctuation of the optical path becomes a problem due to vibration or the like generated during actual use of the product,
In the allowable error specifying step, the allowable amount of the optical path is specified, in the element selecting step, design variables that can determine the optical path are defined in the finite element model, and in the sensitivity calculating step, the design variables are defined in the element selecting step. The sensitivity coefficient for the variation of the optical path of the design variable obtained is obtained, and in the model optimization step, the sensitivity coefficient is extracted in the sensitivity calculation step until the variation amount of the optical path specified in the allowable error designation step falls within an allowable value. Among the design variables, the one with the higher sensitivity coefficient is preferentially operated, the vibration calculation is repeatedly performed, and the optimum design variable can be determined by the computer.

【0024】また、上記の曲面フィッティング工程にお
いて、前記形状関数を空間変数x,y,zの多項式で表
現し、数値モデルによって離散的に計算された変形状態
から最小2乗法によって前記多項式の各係数を求めて曲
面フィッティングを行う。
In the above-mentioned curved surface fitting step, the shape function is represented by a polynomial of spatial variables x, y, and z, and each coefficient of the polynomial is calculated by a least square method from a deformed state discretely calculated by a numerical model. And perform surface fitting.

【0025】この構成により、前記振動計算工程で計算
された有限要素法モデルの節点上のみでの離散的な光学
面の振動結果を、前記光学計算工程で光学ユニットの構
成上考えられる任意の光学経路を計算するために曲面の
方程式に変換する際に、様々な振動状態の表現が可能で
あり、かつ光学面の局部的な変形に対し、光学経路計算
の誤差を小さくすることが可能となる。
With this configuration, the result of the vibration of the discrete optical surface only on the node of the finite element method model calculated in the vibration calculation step can be converted into any optical structure considered in the configuration of the optical unit in the optical calculation step. When converting to a curved surface equation to calculate the path, it is possible to represent various vibration states, and it is possible to reduce the error of the optical path calculation for local deformation of the optical surface .

【0026】また、上記の曲面フィッティング工程にお
いて、曲面の近似精度の許容値を指定する許容値指定工
程と、形状関数の多項式次数を指定する次数指定工程
と、多項式によって計算された光学部品の変形量と測定
された変形量とを比較し、その差が前記許容値内であれ
ば前記次数指定工程において定義された多項式を採用
し、許容値外であれば許容値内に収まるまで多項式次数
を上げて計算する収束計算工程とを有している。
In the above-mentioned curved surface fitting step, an allowable value specifying step of specifying an allowable value of the approximation accuracy of the curved surface, an order specifying step of specifying a polynomial order of the shape function, and a deformation of the optical component calculated by the polynomial Comparing the amount and the measured deformation amount, if the difference is within the allowable value, the polynomial defined in the order designation step is adopted, and if the difference is outside the allowable value, the polynomial order is reduced to within the allowable value. And a convergence calculation step of performing calculation.

【0027】この構成により、前記曲面フィッティング
工程で、離散的な光学面の振動結果を多項式近似により
補間する際、振動モードに対し不必要に高次の多項式次
数を設定すると、かえってフィッティング精度の低下を
招くことを防止することが可能となる。以下、本発明の
一実施形態について複写機の画像形成部への書き込み装
置であるレーザースキャナーの光学設計を対象に説明す
る。
According to this configuration, when the vibration result of the discrete optical surface is interpolated by the polynomial approximation in the curved surface fitting step, if the higher-order polynomial order is set unnecessarily for the vibration mode, the fitting accuracy is reduced. Can be prevented. Hereinafter, an embodiment of the present invention will be described with respect to an optical design of a laser scanner which is a writing device for an image forming unit of a copying machine.

【0028】図1は、本発明の一実施形態の光学機器の
振動解析方法の手順を示すフローチャートである。
FIG. 1 is a flowchart showing a procedure of a method for analyzing vibration of an optical apparatus according to an embodiment of the present invention.

【0029】図1において、まず振動測定工程(ステッ
プS1)において実際に各光学ユニットの振動を測定す
る。
In FIG. 1, first, the vibration of each optical unit is actually measured in a vibration measuring step (step S1).

【0030】図2に複写機の画像形成部への書き込み装
置であるレーザースキャナーユニットを示す。光学機器
は、レーザー照射ユニット2Aと、レンズ群2B,2D
と、ポリゴンミラー2Cとそれらを保持するスキャナー
ユニット2Gと、折り返しミラー2Eとから構成され
る。
FIG. 2 shows a laser scanner unit which is a writing device for an image forming section of a copying machine. The optical device includes a laser irradiation unit 2A and lens groups 2B and 2D.
, A polygon mirror 2C, a scanner unit 2G for holding the polygon mirror 2C, and a folding mirror 2E.

【0031】画像形成の過程は、まずレーザー照射ユニ
ット2Aにより画像のデジタル信号に基づいてレーザー
が照射され、レンズ2Bを介して回転しているポリゴン
ミラー2Cで反射されることでスキャニングが行われ
る。ポリゴンミラー2Cで反射されたレーザーはレンズ
群2Dを介し、折り返しミラー2Eの反射面上でスキャ
ニングスピードが等速度になるよう補正され、押し返し
ミラー2Eで反射されたレーザーは最終的には感光ドラ
ム2F上に照射される。また感光ドラム2Fはその中心
軸まわりに回転しており、スキャンしている主走査方向
に対して垂直方向に照射位置を送り出すことによって画
像を形成する。
In the image forming process, a laser is first irradiated by a laser irradiation unit 2A based on a digital signal of an image, and is reflected by a rotating polygon mirror 2C via a lens 2B to perform scanning. The laser reflected by the polygon mirror 2C is corrected via the lens group 2D so that the scanning speed becomes uniform on the reflection surface of the return mirror 2E, and the laser reflected by the pushback mirror 2E finally becomes the photosensitive drum 2F. Irradiated on top. The photosensitive drum 2F rotates around its central axis, and forms an image by sending out an irradiation position in a direction perpendicular to the main scanning direction in which the image is scanned.

【0032】図3は、振動測定工程(ステップS1)で
用いられる振動測定装置の具体的な構成を示す図であ
る。測定方法としては、製品稼動時において圧電型加速
度センサ3Aを直接各光学部品に取付け、FFTアナラ
イザー3Bにより周波数毎の振幅データとして収集す
る。またインパクトハンマーによるモーダル試験をも実
施し、モーダルパラメータ(固有振動数、固有モード、
減衰比)を収集する。収集された実験データはコンピュ
ータ3C内のハードディスク等の記憶装置3Dに格納さ
れ、数値モデル作成工程(ステップS2)に渡される。
また測定された実稼動時の振動状態およびモーダルパラ
メータはコンピュータに接続された表示装置3H上に表
示され、採取した振動波形データ(3E)、振動モード
(3F)、モーダルパラメータ(3G)等の結果の確認
が可能になっている。
FIG. 3 is a diagram showing a specific configuration of the vibration measuring device used in the vibration measuring step (Step S1). As a measuring method, the piezoelectric acceleration sensor 3A is directly attached to each optical component when the product is operating, and the amplitude data for each frequency is collected by the FFT analyzer 3B. In addition, a modal test using an impact hammer was performed, and modal parameters (eigenfrequency, eigenmode,
Damping ratio). The collected experimental data is stored in a storage device 3D such as a hard disk in the computer 3C, and is passed to a numerical model creation step (Step S2).
The measured vibration state and modal parameters at the time of actual operation are displayed on a display device 3H connected to a computer, and the results of the collected vibration waveform data (3E), vibration mode (3F), modal parameters (3G), etc. Can be checked.

【0033】数値モデル作成工程(ステップS2)で
は、振動計算工程(ステップS3)で光学機器の振動計
算を行うための有限要素法モデル作成と、光学計算工程
(ステップS4)で機器振動時のレーザー軌跡の光学計
算を行うための光学面の定義とを行う。
In the numerical model creation step (step S2), a finite element method model for calculating the vibration of the optical device is calculated in the vibration calculation step (step S3), and the laser generated when the equipment is vibrated in the optical calculation step (step S4). An optical surface for performing the optical calculation of the trajectory is defined.

【0034】図4は、数値モデル作成工程(ステップS
2)で使用される数値モデル作成装置の具体的な動作を
示す図である。記憶装置3Dには有限要素モデルを生成
するためのプログラムが格納されており、振動計算の対
象を定義する。数値モデルは、折り返しミラー2Eやス
キャナーユニット2Gのみならず、それらユニットを保
持する枠体4Bをも含み、光学機器の振動特性に寄与す
る周辺領域が対象となる。作成された有限要素法モデル
は表示装置3H上に表示され、光学経路を決定する光学
面4Dをマウス等のポインティングデバイス4Aによっ
て定義する。振動特性に関係する材料定数(ヤング率、
ポアソン比、質量密度)はキーボード4Eにより入力し
定義すると同時に、振動測定工程(ステップS1)にお
いて測定されたモーダルパラメータの内の減衰比が記憶
装置3Dから読み込まれ、振動計算のための材料定義が
完成する。また光学機器に加わる加振力の振動波形も記
憶装置3Dから読み込まれ、ポインティングデバイス4
Aによって数値モデル内の所定の位置を指定することで
定義される。
FIG. 4 shows a numerical model creation process (step S
It is a figure which shows the specific operation | movement of the numerical model preparation apparatus used by 2). The storage device 3D stores a program for generating a finite element model, and defines a target of vibration calculation. The numerical model includes not only the folding mirror 2E and the scanner unit 2G, but also a frame 4B that holds these units, and covers a peripheral region that contributes to the vibration characteristics of the optical device. The created finite element method model is displayed on the display device 3H, and an optical surface 4D for determining an optical path is defined by a pointing device 4A such as a mouse. Material constants related to vibration characteristics (Young's modulus,
Poisson's ratio, mass density) are input and defined by the keyboard 4E, and at the same time, the damping ratio of the modal parameters measured in the vibration measurement step (step S1) is read from the storage device 3D, and the material definition for vibration calculation is obtained. Complete. Also, the vibration waveform of the excitation force applied to the optical device is read from the storage device 3D, and the pointing device 4
It is defined by designating a predetermined position in the numerical model by A.

【0035】振動計算工程(ステップS3)では、数値
モデル作成工程(ステップS2)で作成された数値モデ
ルを用いて有限要素法による振動解析を実施する。
In the vibration calculation step (step S3), a vibration analysis is performed by the finite element method using the numerical model created in the numerical model creation step (step S2).

【0036】図5は、振動計算工程(ステップS3)で
用いられる振動計算装置の構成および動作を示す図であ
る。振動計算工程(S3)は固有モード計算工程(ステ
ップS31)および周波数応答計算工程(ステップS3
2)とから成り、固有モード計算工程(ステップS3
1)が終了した時点で、振動測定工程(ステップS1)
で測定した固有モードおよび固有振動数を記憶装置3D
から呼び出し、表示装置3H上で計算によって求めた固
有モードおよび固有振動数と比較を行うことで数値モデ
ルの妥当性を確認する。所定の精度が得られない場合に
は、数値モデル作成工程(ステップS2)に戻り数値モ
デルの修正を行う。所定の精度が得られた場合には、周
波数応答計算工程(ステップS32)に進み振動計算を
行い、計算結果が表示装置3Hに表示される。ここで振
動測定工程(ステップS1)において測定された実稼動
時の振動モードが点と線から成るワイヤーフレームモデ
ル5Aとして呼び出され、周波数応答計算工程(ステッ
プS32)での計算結果5Bと重ね合わせることで比較
を行い、再度計算精度の確認が行われる(5C)。もし
所定の精度が得られない場合、数値モデル作成工程(ス
テップS2)に戻り、加振力あるいは減衰比の調整を行
う。所定の精度が得られた場合には計算結果を記憶装置
3Dに格納し光学計算工程(ステップS4)に移る。
FIG. 5 is a diagram showing the configuration and operation of the vibration calculation device used in the vibration calculation step (step S3). The vibration calculation step (S3) includes an eigenmode calculation step (step S31) and a frequency response calculation step (step S3).
2) and an eigenmode calculation step (step S3).
When 1) is completed, a vibration measuring step (step S1)
Eigenmode and eigenfrequency measured by the storage device 3D
And compares it with the eigenmode and eigenfrequency obtained by calculation on the display device 3H to confirm the validity of the numerical model. If the predetermined accuracy cannot be obtained, the process returns to the numerical model creation step (step S2) to correct the numerical model. When the predetermined accuracy is obtained, the process proceeds to the frequency response calculation step (step S32), the vibration is calculated, and the calculation result is displayed on the display device 3H. Here, the vibration mode at the time of actual operation measured in the vibration measurement step (step S1) is called as a wire frame model 5A including points and lines, and is superimposed on the calculation result 5B in the frequency response calculation step (step S32). Are performed, and the calculation accuracy is confirmed again (5C). If the predetermined accuracy cannot be obtained, the process returns to the numerical model creation step (step S2), and the excitation force or the damping ratio is adjusted. When the predetermined accuracy is obtained, the calculation result is stored in the storage device 3D, and the process proceeds to the optical calculation step (Step S4).

【0037】本実施形態では最終的な光学機器の振動状
態を得るために、周波数軸上での振動状態を計算する周
波数応答解析工程(ステップS32)を採用したが、周
波数応答解析の代用として時間軸上での振動状態を計算
する過渡応答解析を採用してもかまわない。
In this embodiment, in order to obtain the final vibration state of the optical apparatus, the frequency response analysis step (step S32) of calculating the vibration state on the frequency axis is employed. A transient response analysis for calculating the vibration state on the shaft may be adopted.

【0038】図6は、光学計算工程(ステップS4)の
構成を示す図である。
FIG. 6 is a diagram showing the configuration of the optical calculation step (step S4).

【0039】光学計算工程(ステップS4)は曲面フィ
ッティング工程(ステップS41)と光学経路計算工程
(ステップS42)とから構成される。
The optical calculation step (step S4) includes a curved surface fitting step (step S41) and an optical path calculation step (step S42).

【0040】曲面フィッティング工程(ステップS4
1)での動作を折り返しミラーの光学面への適応事例を
用いて説明する。
The curved surface fitting step (step S4)
The operation in 1) will be described using an example of application to the optical surface of the folding mirror.

【0041】図7は、振動計算工程において求めた折り
返しミラーの振動状態図である。2Eは、静止時の折り
返しミラーの形状を示している。有限要素法で計算され
た光学面7Aの振動状態は節点7B上の解として得られ
る離散点データの集合体である。従って、レーザービー
ム7Cが折り返しミラー2E上をスキャニングした際の
光学経路を計算するためには、離散データである光学面
の変形状態を連続的な形状関数に変換し、任意の位置で
の反射角等を計算できるようにしなければならない。離
散的な変形データを連続的な形状関数に変換する曲面フ
ィッティングの手法としては種々あるが、後の光学経路
計算工程(ステップS42)における計算方法の簡便さ
から、本実施形態では形状関数として式(1)で定義す
るようなyz多項式を採用している。
FIG. 7 is a vibration state diagram of the folding mirror obtained in the vibration calculation step. 2E shows the shape of the folding mirror at rest. The vibration state of the optical surface 7A calculated by the finite element method is a set of discrete point data obtained as a solution on the node 7B. Therefore, in order to calculate the optical path when the laser beam 7C scans the folding mirror 2E, the deformation state of the optical surface, which is discrete data, is converted into a continuous shape function, and the reflection angle at an arbitrary position is calculated. And so on. There are various methods of curved surface fitting for converting discrete deformation data into a continuous shape function. However, in the present embodiment, a shape function is used as a shape function in the present embodiment because of a simple calculation method in an optical path calculation step (step S42). The yz polynomial as defined in (1) is adopted.

【0042】[0042]

【数1】 この多項式は図7のように光学面の座標系をyz平面に
選び、振動による移動量をスカラーベクトルXとしてい
る。式(1)中の各yz項に対する係数Cを図8に示
す。曲面フィッティングの具体的な方法は、振動計算工
程(ステップS3)で求めた光学面の節点変位ベクトル
Xおよび節点座標yzを記憶装置より呼び出し、最小2
乗法により式(1)に対し最も誤差が小さくなるような
係数群CNを求める。フィッティングを行う際、扱う多
項式次数を何次にするかがフィッティング精度に影響を
与える。対象次数が少なすぎると十分な振動形状を表現
できなくなり、一方多すぎても高次の影響による誤差が
発生する。したがって、光学面における曲面フィッティ
ングの結果と節点上の振動計算の結果は図9のように表
示装置によって可視化され、フィッティングの様子が確
認できるようになっている。また図10のようにフィッ
ティングエラーおよび多項式係数を表示することが可能
である。
(Equation 1) In this polynomial, the coordinate system of the optical surface is selected on the yz plane as shown in FIG. FIG. 8 shows the coefficient C for each yz term in equation (1). A specific method of the curved surface fitting is that the node displacement vector X and the node coordinates yz of the optical surface obtained in the vibration calculation step (step S3) are called from the storage device, and the minimum 2
A coefficient group CN that minimizes the error with respect to equation (1) is obtained by the multiplication method. When performing fitting, the order of the polynomial order to be handled affects fitting accuracy. If the target order is too small, it is not possible to express a sufficient vibration shape, while if it is too large, an error due to higher-order effects will occur. Therefore, the result of the curved surface fitting on the optical surface and the result of the vibration calculation on the node are visualized by the display device as shown in FIG. 9 so that the state of the fitting can be confirmed. Further, it is possible to display the fitting error and the polynomial coefficient as shown in FIG.

【0043】上記の処理により、十分なフィッティング
精度が得られていない場合には、曲面フィッティング工
程(ステップS41)に戻り、式(1)中の多項式次数
を再度設定してフィッティングを行い、また十分な精度
が得られた場合には、光学経路計算工程(ステップS4
2)に移る。
If sufficient fitting accuracy is not obtained by the above processing, the process returns to the curved surface fitting step (step S41), and the polynomial degree in equation (1) is set again to perform fitting. If a high accuracy is obtained, the optical path calculation step (step S4)
Go to 2).

【0044】光学経路計算工程(ステップS42)では
光学面の変形によってレーザービームの光学経路がどの
ように変化するのかを計算し、最終的な照射位置を求め
る。光学面の形状関数fは式(2)により得られ、折り
返しミラーに照射されるレーザービームのスカラーベク
トルをSとすると、式(3)によって反射位置での方向
微分係数が容易に計算でき、レーザービームが振動して
いる光学面で反射してどの方向に進むのかが求まる。
In the optical path calculation step (step S42), how the optical path of the laser beam changes due to the deformation of the optical surface is calculated, and the final irradiation position is obtained. The shape function f of the optical surface is obtained by equation (2). If the scalar vector of the laser beam applied to the folding mirror is S, the directional differential coefficient at the reflection position can be easily calculated by equation (3). The direction in which the beam is reflected by the oscillating optical surface and travels is determined.

【0045】[0045]

【数2】 図2に示したレーザー照射ユニット2A、レンズ群2
B,2D、ポリゴンミラー2Cの光学面についても前述
の折り返しミラーと同様の光学経路計算を実施し、最終
的に感光ドラム2F上に照射されるレーザースポット位
置を計算する。計算されたレーザースポット位置は画像
表示装置によって確認できる。
(Equation 2) Laser irradiation unit 2A and lens group 2 shown in FIG.
For the optical surfaces B, 2D, and the polygon mirror 2C, the same optical path calculation as that of the above-described folding mirror is performed, and the laser spot position finally irradiated on the photosensitive drum 2F is calculated. The calculated laser spot position can be confirmed by the image display device.

【0046】図11に機器振動時において、8角形であ
るポリゴンミラー2Cの2回転分に相当する計16本の
直線を感光ドラム2F上に描画した際の画像出力例を示
す。ポリゴンミラー2Cおよび感光ドラム2Fの回転数
は一定であるため、実際には16本の直線が等間隔で引
けるはずであるが、各光学面の振動により直線間隔がま
ばらとなる様子が確認できる。もし振動による影響が与
える出力画像品質が許容できない場合、モデル最適化工
程(ステップS5)に移行する。
FIG. 11 shows an example of image output when a total of 16 straight lines corresponding to two rotations of the octagonal polygon mirror 2C are drawn on the photosensitive drum 2F when the apparatus is vibrated. Since the rotation speeds of the polygon mirror 2C and the photosensitive drum 2F are constant, 16 straight lines should actually be able to be drawn at equal intervals. However, it can be seen that the linear intervals are sparse due to the vibration of each optical surface. If the output image quality affected by the vibration is unacceptable, the process proceeds to the model optimization step (step S5).

【0047】図12に光学機器の設計フロー内でのモデ
ル最適化工程(ステップS5)の詳細な構成を示す。モ
デル最適化工程(ステップS5)は許容誤差指定工程
(ステップS51)、要素選択工程(ステップS52)
および感度計算工程(ステップS53)によって構成さ
れる。
FIG. 12 shows the detailed configuration of the model optimizing step (step S5) in the design flow of the optical equipment. The model optimization step (step S5) includes an allowable error designation step (step S51) and an element selection step (step S52).
And a sensitivity calculation step (step S53).

【0048】許容誤差指定工程(ステップS51)で
は、光学機器の設計における所定の位置でのレーザービ
ームスポットの許容変動量を指定する。
In the allowable error designating step (step S51), the allowable variation of the laser beam spot at a predetermined position in the design of the optical equipment is designated.

【0049】以下に図13を用いて本工程での動作を説
明する。許容誤差指定工程(ステップS51)では、数
値モデル作成工程(ステップS3)において作成された
有限要素モデルが表示装置3H上に表示され、振動時レ
ーザービームスポットの許容変動量をどの位置で定義す
るかをポインティングデバイス等で指示する。また同時
に許容誤差量を入力する。本実施形態においては、感光
ドラム2F上で許容誤差量0.1μmとなるように設定
する。
The operation in this step will be described below with reference to FIG. In the allowable error designating step (step S51), the finite element model created in the numerical model creating step (step S3) is displayed on the display device 3H, and the position at which the allowable variation amount of the laser beam spot during vibration is defined Is instructed by a pointing device or the like. At the same time, an allowable error amount is input. In the present embodiment, the allowable error amount is set to be 0.1 μm on the photosensitive drum 2F.

【0050】要素選択工程(ステップS52)ではレー
ザービームのスポットの変動量に対し、光学機器構成部
品の剛性や質量密度の感度を感度計算工程(ステップS
53)において計算するにあたり、感度計算の対象とす
る部品の選択を、再度有限要素モデル内の要素群を選択
することで定義する。
In the element selecting step (step S52), the sensitivity calculation step (step S52) determines the stiffness and mass density sensitivity of the components of the optical device with respect to the variation of the spot of the laser beam.
In the calculation in 53), selection of a component to be subjected to sensitivity calculation is defined by selecting an element group in the finite element model again.

【0051】以下に図14を用いて本工程の動作を説明
する。本実施形態では折り返しミラー2Eの構成要素で
あるミラー14A、ミラーホルダー14B、および取り
付けバネ14Cの材料特性である剛性(ヤング率)を感
度計算の対象として選択する。表示装置3H上に表示さ
れた有限要素モデル内の要素をポインティングデバイス
等で選択すると、同一部品(同一材質)全てが感度計算
の対象要素として認識される。同作業を部品毎に繰り返
し、感度計算の対象要素を定義する。また選択された部
品は表示装置3H上でリスト表示され確認が容易になっ
ている。
The operation of this step will be described below with reference to FIG. In the present embodiment, rigidity (Young's modulus), which is a material characteristic of the mirror 14A, the mirror holder 14B, and the mounting spring 14C, which are the components of the folding mirror 2E, is selected as an object of the sensitivity calculation. When an element in the finite element model displayed on the display device 3H is selected by a pointing device or the like, all the same parts (the same material) are recognized as the sensitivity calculation target elements. This operation is repeated for each part, and the target element of the sensitivity calculation is defined. Also, the selected parts are displayed in a list on the display device 3H to facilitate confirmation.

【0052】感度計算工程(ステップS53)では、要
素選択工程(ステップS52)において選択された部品
の材料特性(本実施形態では要素のヤング率)のスポッ
ト変動量に対する感度を計算する。本実施形態における
感度とは、レーザービームのスポット変動量を要素のヤ
ング率の増減比率で除したものである。この感度が大き
い部品程、その剛性の変更がスポット変動量の抑制に対
し有効であることを意味している。
In the sensitivity calculation step (step S53), the sensitivity to the spot variation of the material characteristics (the Young's modulus of the element in the present embodiment) of the component selected in the element selection step (step S52) is calculated. The sensitivity in the present embodiment is obtained by dividing the spot fluctuation amount of the laser beam by the increase / decrease ratio of the Young's modulus of the element. This means that the higher the sensitivity of the part, the more the change in the rigidity is more effective in suppressing the spot fluctuation amount.

【0053】以下図15を用いて本工程の動作を説明す
る。要素選択工程(ステップS52)で選択された部品
の剛性を一律10%増加(あるいは減少)させ、振動計
算工程(ステップS3)、光学計算工程(S4)を実施
し部品剛性のスポット変動量に対する感度を算出する。
本実施形態では10%としたが、増減比率は各部品とも
一定であれば、10%である必要はない。算出された感
度の結果は、感度数値の大小によって有限要素モデル中
の要素に対し色分けされ、表示装置3Hにて可視化され
るので、どの部品剛性がレーザービームのスポット変動
量に寄与しているのかが容易に判断できる。本実施形態
では感度数値が大きいほど暖かい色で示す仕様にしてお
り、取り付けバネ14Cが最も感度が高い。感度計算の
結果、光学機器を構成する各部品のうちどの部品の改善
に注力するか目安を立て、数値モデル作成工程(ステッ
プS2)に戻り、具体的な対策モデルを再構築する。本
作業を許容誤差指定工程(ステップS51)で指定した
許容誤差内に収まるまで繰り返す。
The operation of this step will be described below with reference to FIG. The rigidity of the part selected in the element selection step (step S52) is increased (or decreased) by 10% uniformly, and the vibration calculation step (step S3) and the optical calculation step (S4) are performed, and the sensitivity of the part rigidity to the spot variation is performed. Is calculated.
In the present embodiment, it is 10%. However, if the change ratio is constant for each component, it is not necessary to be 10%. The calculated sensitivity results are color-coded for the elements in the finite element model according to the magnitude of the sensitivity value and are visualized on the display device 3H, so which component rigidity contributes to the spot variation amount of the laser beam. Can be easily determined. In this embodiment, the specification is such that the larger the sensitivity value, the warmer the color, and the mounting spring 14C has the highest sensitivity. As a result of the sensitivity calculation, a guide is made as to which one of the components constituting the optical apparatus should be focused on, and the process returns to the numerical model creation step (step S2) to reconstruct a specific countermeasure model. This operation is repeated until it falls within the allowable error specified in the allowable error specifying step (step S51).

【0054】以上説明したように、本実施形態によれ
ば、以下の効果を得ることができる。
As described above, according to the present embodiment, the following effects can be obtained.

【0055】振動測定工程では、有限要素法等の数値シ
ミュレーションのみでは予測しきれない各光学ユニット
の振動数毎の減衰効果を数値モデル作成工程に渡し、数
値モデル作成工程では、各光学ユニットの詳細形状を有
限要素法モデルにて定義することで、測定では機器内へ
のセンサの配線等の物理的な問題で不十分となる測定点
を補い、振動計算工程では、実際の振動測定では測定困
難な光学ユニットの振動成分(特に回転変動)を正確に
計算し、曲面フィッティング工程では、振動計算工程に
よって計算された光学面の振動状態が有限要素法モデル
の特性上、節点上での空間的な離散点データとして表現
されるので、任意の位置での光学計算を行うために、光
学面の離散点での振動状態をもとに3次元空間に対する
一意の形状表現が可能な形状関数に変換して光学計算工
程に渡し、光学計算工程では、各光学面の振動モードを
考慮した光学経路の計算を行うことが可能となるため、
実稼動時に発生する光学ユニットの振動変形が光学経路
に与える影響を正確に考慮して、各光学ユニットの配置
を決定することが可能となる。
In the vibration measuring step, the damping effect for each frequency of each optical unit, which cannot be predicted only by the numerical simulation such as the finite element method, is passed to the numerical model creating step. By defining the shape with a finite element method model, it compensates for measurement points that are insufficient due to physical problems such as wiring of sensors into the device in the measurement, and it is difficult to measure in the vibration calculation process with actual vibration measurement The vibration component (especially rotation fluctuation) of a simple optical unit is accurately calculated, and in the curved surface fitting process, the vibration state of the optical surface calculated by the vibration calculation process is spatially limited at the nodes due to the characteristics of the finite element method model. Since it is expressed as discrete point data, in order to perform optical calculation at an arbitrary position, a unique shape expression for a three-dimensional space based on the vibration state at discrete points on the optical surface For passing the optical calculations step by converting the ability shape function, the optical calculation process, it becomes possible to perform the calculation of the optical path of which the oscillation mode of each optical surface,
The arrangement of each optical unit can be determined by accurately considering the influence of the vibration deformation of the optical unit generated during actual operation on the optical path.

【0056】また、製品の実使用時に発生する振動等に
より、光学経路の変動量が問題になる場合、許容誤差指
定工程にて光学経路の許容量を指定し、要素選択工程に
て、光学経路を決定しうる設計変数を有限要素モデルの
中で定義し、感度計算工程にて、要素選択工程で定義さ
れた設計変数の光学経路の変動に対する感度係数を求
め、モデル最適化工程にて、許容誤差指定工程で指定さ
れた光学経路の変動量が許容値内に収まるまで、感度計
算工程で抽出された設計変数の内、感度係数が高いもの
を優先的に操作し、繰り返し振動計算を行い計算機によ
って最適な設計変数の決定が可能となる。
When the amount of fluctuation of the optical path becomes a problem due to vibration or the like generated during actual use of the product, the allowable amount of the optical path is specified in the allowable error specifying step, and the optical path is specified in the element selecting step. In the finite element model, the design variables that can determine the parameters are determined, and in the sensitivity calculation step, the sensitivity coefficient to the variation of the optical path of the design variables defined in the element selection step is obtained. Until the variation amount of the optical path specified in the error specification process falls within the allowable value, the design variable extracted in the sensitivity calculation process with the higher sensitivity coefficient is preferentially operated, and the vibration is calculated repeatedly. This makes it possible to determine optimal design variables.

【0057】また、振動計算工程で計算された有限要素
法モデルの節点上のみでの離散的な光学面の振動結果
を、光学計算工程で光学ユニットの構成上考えられる任
意の光学経路を計算するために曲面の方程式に変換する
際に、様々な振動状態の表現が可能であり、かつ光学面
の局部的な変形に対し、光学経路計算の誤差を小さくす
ることが可能となる。
The result of the vibration of the discrete optical surface only on the node of the finite element method model calculated in the vibration calculation step is used to calculate any optical path conceivable in the configuration of the optical unit in the optical calculation step. Therefore, when converting into an equation of a curved surface, it is possible to represent various vibration states, and it is possible to reduce the error of the optical path calculation for local deformation of the optical surface.

【0058】また、曲面フィッティング工程で、離散的
な光学面の振動結果を多項式近似により補間する際、振
動モードに対し不必要に高次の多項式次数を設定する
と、かえってフィッティング精度の低下を招くことを防
止する効果がある。
Further, in the curved surface fitting step, when the vibration results of the discrete optical surfaces are interpolated by the polynomial approximation, if the polynomial order of the vibration mode is set unnecessarily high, the fitting accuracy is reduced. Has the effect of preventing.

【0059】[0059]

【発明の効果】以上説明した様に、本発明によれば、光
学機器の実稼動時の振動が光学機器の性能に与える影響
を正確に知ることができる光学機器の振動解析方法及び
装置を提供することが可能となる。
As described above, according to the present invention, there is provided a method and apparatus for analyzing the vibration of an optical device which can accurately determine the influence of the vibration during the actual operation of the optical device on the performance of the optical device. It is possible to do.

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

【図1】本発明の一実施形態の光学機器の振動解析方法
の手順を示すフローチャートである。
FIG. 1 is a flowchart illustrating a procedure of a vibration analysis method for an optical device according to an embodiment of the present invention.

【図2】レーザースキャナーユニットの構成を示す概念
図である。
FIG. 2 is a conceptual diagram illustrating a configuration of a laser scanner unit.

【図3】振動測定装置の具体的な構成および動作を示す
説明図である。
FIG. 3 is an explanatory diagram showing a specific configuration and operation of the vibration measuring device.

【図4】数値モデル作成装置の動作を示す説明図であ
る。
FIG. 4 is an explanatory diagram showing the operation of the numerical model creation device.

【図5】振動計算装置の構成例および動作を示す説明図
である。
FIG. 5 is an explanatory diagram showing a configuration example and operation of a vibration calculation device.

【図6】光学計算工程の構成を示す説明図である。FIG. 6 is an explanatory diagram showing a configuration of an optical calculation step.

【図7】振動解析工程において求めた折り返しミラーの
振動状態を示す説明図である。
FIG. 7 is an explanatory diagram showing a vibration state of a return mirror obtained in a vibration analysis step.

【図8】曲面フィッティングにおける多項式係数の説明
図である。
FIG. 8 is an explanatory diagram of polynomial coefficients in curved surface fitting.

【図9】曲面フィッティング工程の動作を示す説明図で
ある。
FIG. 9 is an explanatory diagram showing an operation in a curved surface fitting step.

【図10】曲面フィッティング工程の動作を示す説明図
である。
FIG. 10 is an explanatory diagram showing an operation in a curved surface fitting step.

【図11】光学計算工程での計算結果を示す説明図であ
る。
FIG. 11 is an explanatory diagram showing a calculation result in an optical calculation step.

【図12】光学機器の設計フロー内でのモデル最適化工
程(ステップS5)の構成を示す説明図である。
FIG. 12 is an explanatory diagram showing a configuration of a model optimization step (step S5) in a design flow of an optical device.

【図13】許容誤差指定工程の動作を示す説明図であ
る。
FIG. 13 is an explanatory diagram showing an operation of an allowable error designation step.

【図14】要素選択工程の動作を示す説明図である。FIG. 14 is an explanatory diagram showing an operation of an element selection step.

【図15】感度計算工程の結果を示す説明図である。FIG. 15 is an explanatory diagram showing a result of a sensitivity calculation step.

【符号の説明】[Explanation of symbols]

2A レーザー照射ユニット 2B,2D レンズ郡 2C ポリゴンミラー 2E 折り返しミラー 2F 感光ドラム 2G スキャナーユニット 3A 圧電型加速度センサ 3B FFTアナライザ 3C コンピュータ 3D 記憶装置 3E 表示装置に表示された振動波形データ 3F 表示装置に表示された振動モード 3G 表示装置に表示されたモーダルパラメータ 3H 表示装置 4A ポインティングデバイス 4B 枠体 4D 光学面 4E キーボード 5A 測定された実稼動モードのワイヤーフレーム 5B 計算された実稼動モードのFEMモデル 5C 表示装置に表示された実稼動モード 7A 光学面 7B 振動時光学面 7C レーザービーム 14A ミラー 14B ミラーホルダー 14C 取り付けバネ 2A laser irradiation unit 2B, 2D lens group 2C polygon mirror 2E folding mirror 2F photosensitive drum 2G scanner unit 3A piezoelectric acceleration sensor 3B FFT analyzer 3C computer 3D storage device 3E Vibration waveform data displayed on display device 3F Displayed on display device Vibration mode 3G Modal parameter displayed on 3G display device 3H Display device 4A Pointing device 4B Frame 4D Optical surface 4E Keyboard 5A Measured actual operation mode wireframe 5B Calculated actual operation mode FEM model 5C Display device Actual operation mode displayed 7A Optical surface 7B Vibration optical surface 7C Laser beam 14A Mirror 14B Mirror holder 14C Mounting spring

Claims (8)

【特許請求の範囲】[Claims] 【請求項1】 光学機器の実稼動状況下における、該光
学機器を構成する光学部品の振動が前記光学機器の性能
に与える影響を解析するための光学機器の振動解析方法
であって、 前記光学部品の振動数及び振幅を実際に測定する振動測
定工程と、 測定された前記光学部品の振動数及び振幅の情報に基づ
いて有限要素法のシミュレーションモデルを作成する数
値モデル作成工程と、 該数値モデル作成工程において生成されたシミュレーシ
ョンモデルを用いて振動計算を行う振動計算工程と、 該振動計算工程において計算された前記光学部品の振動
による変形状態を、3次元空間に対する一意の形状表現
が可能な形状関数で近似する曲面フィッティング工程
と、 前記形状関数に基づいて前記光学機器の性能に関係する
特性の計算を行う光学計算工程と、 を具備することを特徴とする光学機器の振動解析方法。
1. A vibration analysis method for an optical device for analyzing the influence of the vibration of an optical component constituting the optical device on the performance of the optical device in an actual operation state of the optical device, the method comprising: A vibration measuring step of actually measuring the frequency and amplitude of the component; a numerical model creating step of creating a simulation model of the finite element method based on the measured information of the frequency and amplitude of the optical component; A vibration calculation step of performing a vibration calculation using the simulation model generated in the creation step; and a shape capable of expressing a deformation state due to the vibration of the optical component calculated in the vibration calculation step in a three-dimensional space. A curved surface fitting step approximating with a function, and an optical calculation for calculating characteristics related to the performance of the optical device based on the shape function Vibration analysis method of the optical apparatus characterized by comprising the extent, the.
【請求項2】 前記数値モデル作成工程は、光学経路の
許容誤差を指定する許容誤差指定工程と、前記シミュレ
ーションモデル中の有限要素を選択する要素選択工程
と、選択された有限要素の材料特性及び/又は形状特性
を変化させて振動計算を行い、前記振動測定工程で測定
された振動数及び振幅に対する前記有限要素の材料特性
及び/又は形状特性の感度を調べる感度計算工程と、計
算された感度に基づいて、前記許容誤差内に収まるま
で、有限要素の材料特性及び/又は形状特性を変更して
収束計算を行うモデル最適化工程とを備えることを特徴
とする請求項1に記載の光学機器の振動解析方法。
2. The numerical model creating step includes: a tolerance specifying step of designating a tolerance of an optical path; an element selecting step of selecting a finite element in the simulation model; A sensitivity calculation step of performing a vibration calculation by changing the shape characteristics, and examining the sensitivity of the material characteristics and / or the shape characteristics of the finite element to the frequency and amplitude measured in the vibration measurement step; The optical apparatus according to claim 1, further comprising: a model optimization step of performing a convergence calculation by changing a material property and / or a shape property of a finite element until the finite element falls within the allowable error. Vibration analysis method.
【請求項3】 前記曲面フィッティング工程では、前記
形状関数を空間変数x,y,zの多項式で表現し、前記
シミュレーションモデルによって離散的に計算された変
形状態から最小2乗法によって前記多項式の各係数を求
めて曲面フィッティングを行うことを特徴とする請求項
1に記載の光学機器の振動解析方法。
3. In the curved surface fitting step, the shape function is represented by a polynomial of spatial variables x, y, z, and each coefficient of the polynomial is calculated by a least square method from a deformation state discretely calculated by the simulation model. 2. The method for analyzing vibration of an optical device according to claim 1, wherein a curved surface fitting is performed by calculating the following equation.
【請求項4】 前記曲面フィッティング工程は、前記多
項式による曲面の近似精度の許容値を指定する許容値指
定工程と、前記多項式の次数を指定する次数指定工程
と、前記多項式によって計算された光学部品の変形量と
測定された変形量とを比較し、その差が前記許容値内で
あれば前記次数指定工程において定義された多項式を採
用し、許容値外であれば許容値内に収まるまで前記多項
式の次数を増加させて計算する収束計算工程とを備える
ことを特徴とする請求項3に記載の光学機器の振動解析
方法。
4. The curved surface fitting step includes a tolerance value designation step of designating a tolerance value of approximation accuracy of the curved surface by the polynomial, a degree designation step of designating an order of the polynomial, and an optical component calculated by the polynomial. The deformation amount is compared with the measured deformation amount, and if the difference is within the allowable value, the polynomial defined in the degree designation step is adopted. 4. The method according to claim 3, further comprising a convergence calculation step of calculating by increasing the degree of the polynomial.
【請求項5】 光学機器の実稼動状況下における、該光
学機器を構成する光学部品の振動が前記光学機器の性能
に与える影響を解析するための光学機器の振動解析装置
であって、 前記光学部品の振動数及び振幅を実際に測定する振動測
定手段と、 測定された前記光学部品の振動数及び振幅の情報に基づ
いて有限要素法のシミュレーションモデルを作成する数
値モデル作成手段と、 該数値モデル作成手段により生成されたシミュレーショ
ンモデルを用いて振動計算を行う振動計算手段と、 該振動計算手段により計算された前記光学部品の振動に
よる変形状態を、3次元空間に対する一意の形状表現が
可能な形状関数で近似する曲面フィッティング手段と、 前記形状関数に基づいて前記光学機器の性能に関係する
特性の計算を行う光学計算手段と、 を具備することを特徴とする光学機器の振動解析装置。
5. An apparatus for analyzing vibration of an optical device for analyzing the influence of the vibration of an optical component constituting the optical device on the performance of the optical device in an actual operation state of the optical device, wherein: Vibration measurement means for actually measuring the frequency and amplitude of the component; numerical model creation means for creating a simulation model of the finite element method based on the measured information on the frequency and amplitude of the optical component; Vibration calculating means for calculating vibration using the simulation model generated by the generating means; and a shape capable of expressing a deformation state of the optical component due to vibration calculated by the vibration calculating means in a three-dimensional space. Curved surface fitting means approximating with a function, and optical calculation means for calculating characteristics related to the performance of the optical device based on the shape function , Vibration analyzer of an optical instrument, characterized by comprising.
【請求項6】 前記数値モデル作成手段は、光学経路の
許容誤差を指定する許容誤差指定手段と、前記シミュレ
ーションモデル中の有限要素を選択する要素選択手段
と、選択された有限要素の材料特性及び/又は形状特性
を変化させて振動計算を行い、前記振動測定手段により
測定された振動数及び振幅に対する前記有限要素の材料
特性及び/又は形状特性の感度を調べる感度計算手段
と、計算された感度に基づいて、前記許容誤差内に収ま
るまで、有限要素の材料特性及び/又は形状特性を変更
して収束計算を行うモデル最適化手段とを備えることを
特徴とする請求項5に記載の光学機器の振動解析装置。
6. The numerical model creating means includes a tolerance designating means for designating a tolerance of an optical path, an element selecting means for selecting a finite element in the simulation model, a material property of the selected finite element, And / or sensitivity calculation means for performing a vibration calculation by changing the shape characteristic and examining the sensitivity of the material property and / or shape property of the finite element to the frequency and amplitude measured by the vibration measurement means; 6. An optical apparatus according to claim 5, further comprising: a model optimizing unit that changes a material property and / or a shape property of a finite element and performs a convergence calculation until the finite element falls within the allowable error. Vibration analysis equipment.
【請求項7】 前記曲面フィッティング手段は、前記形
状関数を空間変数x,y,zの多項式で表現し、前記シ
ミュレーションモデルによって離散的に計算された変形
状態から最小2乗法によって前記多項式の各係数を求め
て曲面フィッティングを行うことを特徴とする請求項5
に記載の光学機器の振動解析装置。
7. The surface fitting means expresses the shape function by a polynomial of spatial variables x, y, z, and calculates each coefficient of the polynomial by a least square method from a deformation state discretely calculated by the simulation model. 6. A curved surface fitting is performed by calculating
A vibration analysis device for an optical device according to item 1.
【請求項8】 前記曲面フィッティング手段は、前記多
項式による曲面の近似精度の許容値を指定する許容値指
定手段と、前記多項式の次数を指定する次数指定手段
と、前記多項式によって計算された光学部品の変形量と
測定された変形量とを比較し、その差が前記許容値内で
あれば前記次数指定手段により定義された多項式を採用
し、許容値外であれば許容値内に収まるまで前記多項式
の次数を増加させて計算する収束計算手段とを備えるこ
とを特徴とする請求項7に記載の光学機器の振動解析装
置。
8. The surface fitting unit includes: a tolerance designating unit that designates a tolerance of approximation accuracy of a curved surface by the polynomial; an order designating unit that designates an order of the polynomial; and an optical component calculated by the polynomial. The deformation amount is compared with the measured deformation amount, and if the difference is within the allowable value, the polynomial defined by the degree designation means is adopted. The apparatus according to claim 7, further comprising: a convergence calculation unit configured to calculate by increasing a degree of the polynomial.
JP2001071091A 2001-03-13 2001-03-13 Method and apparatus for analyzing vibration of optical equipment Expired - Fee Related JP3595775B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001071091A JP3595775B2 (en) 2001-03-13 2001-03-13 Method and apparatus for analyzing vibration of optical equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001071091A JP3595775B2 (en) 2001-03-13 2001-03-13 Method and apparatus for analyzing vibration of optical equipment

Publications (2)

Publication Number Publication Date
JP2002267960A true JP2002267960A (en) 2002-09-18
JP3595775B2 JP3595775B2 (en) 2004-12-02

Family

ID=18928863

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001071091A Expired - Fee Related JP3595775B2 (en) 2001-03-13 2001-03-13 Method and apparatus for analyzing vibration of optical equipment

Country Status (1)

Country Link
JP (1) JP3595775B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115017754A (en) * 2022-05-09 2022-09-06 南京航空航天大学 Finite element model correction method considering manufacturing error

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102015217200A1 (en) * 2015-09-09 2017-03-09 Sauer Gmbh Method and device for determining a vibration amplitude of a tool
CN106847054A (en) * 2016-11-25 2017-06-13 安徽省春谷3D打印智能装备产业技术研究院有限公司 Laser scanning controls Teaching Module
CN109357747B (en) * 2018-10-25 2019-11-12 南京航空航天大学 A kind of identification of online train and speed estimation method based on fiber-optic vibration signal

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115017754A (en) * 2022-05-09 2022-09-06 南京航空航天大学 Finite element model correction method considering manufacturing error

Also Published As

Publication number Publication date
JP3595775B2 (en) 2004-12-02

Similar Documents

Publication Publication Date Title
JP5064864B2 (en) Optical deflection apparatus, image forming apparatus, and driving method of optical deflection apparatus
JP6821343B2 (en) Image forming device
KR20060053235A (en) Information processing apparatus and its control method
JP2018066731A (en) Determining position of movable part of coordinate measuring machine
JP3595775B2 (en) Method and apparatus for analyzing vibration of optical equipment
JPH07324928A (en) Method for measuring coordinate of surface of workpiece
JP3199047B2 (en) Optical scanning device
Rothberg et al. A universal framework for modelling measured velocity in laser vibrometry with applications
JP2008299297A (en) Rocking body apparatus and method of driving and controlling oscillation system of rocking body apparatus
JP2009098032A (en) Oscillating body apparatus, light deflector, and resonance frequency sensing method
JP5009840B2 (en) Optical device design support method and recording medium
JP2024507061A (en) Testing and manufacturing method of mainspring for watches
Vega et al. On the precision of optical imaging to study free surface dynamics at high frame rates
Yu et al. Design sensitivity analysis of acoustical damping and its application to design of musical bells
Eichenberger et al. Introduction to multipath Doppler vibrometry (MDV) for validating complex models accurately and without contact
US6546357B2 (en) Estimation of the configuration of an optical element for an optical writing device
Pesaresi et al. Automated measurement grid generation for scanning laser doppler vibrometers
JP6590634B2 (en) Shape measuring method and shape measuring apparatus
Anac et al. Model validation and performance prediction in the design of micro systems
JPH11119136A (en) Design aiding method of optical instrument, device therefor and recording medium recording program therefor
JP2001142917A (en) Cad device and method for optical pate simulation
JPWO2013084565A1 (en) Data correction apparatus and data correction program
JP4170088B2 (en) Reading position evaluation system and evaluation method
WO2022091458A1 (en) Information processing method, information processing device, and information processing system
Anac et al. Model validation methodology for design of micro systems

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040507

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040528

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040722

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040820

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040906

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070910

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080910

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090910

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090910

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100910

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100910

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110910

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110910

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120910

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120910

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130910

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees