CN1840746A - Control system and method for crystal isodiametric growth - Google Patents

Control system and method for crystal isodiametric growth Download PDF

Info

Publication number
CN1840746A
CN1840746A CN 200510059325 CN200510059325A CN1840746A CN 1840746 A CN1840746 A CN 1840746A CN 200510059325 CN200510059325 CN 200510059325 CN 200510059325 A CN200510059325 A CN 200510059325A CN 1840746 A CN1840746 A CN 1840746A
Authority
CN
China
Prior art keywords
crystal
signal
control
controller
speed
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
CN 200510059325
Other languages
Chinese (zh)
Other versions
CN100371507C (en
Inventor
荀建华
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CNB2005100593258A priority Critical patent/CN100371507C/en
Publication of CN1840746A publication Critical patent/CN1840746A/en
Application granted granted Critical
Publication of CN100371507C publication Critical patent/CN100371507C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)
  • Crystals, And After-Treatments Of Crystals (AREA)

Abstract

The invention relates to a control system and a control method for crystal isodiametric growth. The system is provided with a crystal growth furnace, a crystal diameter detection part and a crucible temperature control part; a crystal pulling-up speed control part for collecting real-time pulling-up speed signals of the crystal and feeding back and adjusting the pulling-up speed of the crystal; there is also a common controller; the controller respectively receives detection signals from the crystal diameter detection part, the crystal pulling speed control part and the crucible temperature control part, and the controller calculates a crystal pulling speed regulation signal by using an algorithm combining Fuzzy control and PID control and feeds the crystal pulling speed regulation signal back to the crystal pulling speed control part to regulate the crystal pulling speed, and calculates a heating power regulation signal and feeds the heating power regulation signal back to the crucible temperature control part to regulate the temperature of the melt. The system achieves excellent control quality, greatly improves the diameter control precision of the crystal in the equal-diameter growth stage, and stably improves the one-time drawing yield of the large-diameter silicon single crystal rod.

Description

The Controlling System of isodiametric growth of crystal and method thereof
Technical field
The present invention relates to a kind of Controlling System and method thereof of isodiametric growth of crystal, relate in particular to the Controlling System and the method thereof of silicon single crystal isodiametric growth.
Background technology
Semiconductor material is the base mateiral of semi-conductor industry, and current semiconducter device more than 95% is made with silicon materials, and the unicircuit more than 99% is a silicon integrated circuit, and producing the needed raw material of unicircuit is silicon single crystal.In recent years, the whole world has strengthened the exploitation dynamics to this eco-friendly power source of sun power, the solar power station that with conversion of solar energy is civilian electric power energy rises just rapidly as an industry, simultaneously to solar energy highly effective transform novel material---the demand of silicon single crystal also is doubled and redoubled, but domestic producer that can the manufacture order crystal silicon is very few, and inefficiency of production far can not satisfy the demand of domestic market.
The technology of traditional manufacture order crystal silicon is when the silicon crystalline growth begins, and makes with the way that lifts.In the monocrystalline silicon production process, the crystalline pull rate can have influence on the crystalline quality.Simultaneously, the process of the crystal growth in the single crystal silicon growth furnace has non-linear and characteristics large time delay, uses conventional PID (proportion integration differentiation) control and is difficult to achieve effective control.Its reason is that whole process of production mainly is divided into the fusing of silicon raw material, seeding, shouldering, commentaries on classics shoulder, isodiametric growth and six successive stage of ending.In these six successive stage crystal growing process non-linear as having, the time become, the controlled member of large time delay, as isodiametric growth.Though present many control theories and technology are increasingly mature, the effect of using conventional PID controller control is still undesirable.And methods such as employing self-adaptive PID, computing time is long, and debugging cost is big, in real time application difficult.
With CZ method (Czochralski method, have in the Chinese article and be translated into Czochralski method) during pulling of silicon single crystal, the diameter control accuracy in the isometrical stage of single crystal silicon growth is relevant with the overall performance of equipment, as factors such as melt temperature, dynamic property in pull rate control, the crystal growing process.Therefore, particularly for the Controlling System of major diameter single crystal stove, traditional monotype is regulated and can not be reached ideal diameter control effect.
Diameter of silicon single crystal control method and equipment thereof are disclosed as Chinese patent application 91102922.2.In this application, the pulling single crystal measuring diameter value that Optical devices are measured with require diameter value to compare, to determine deviation, again the deviation that draws is carried out incomplete differential PID processing or smith's method processing, to calculate pull rate, again pull rate is added to crystal and lifts on the motor controller of equipment, thereby carry the diameter that crystal draws high body by the control of control pull rate.This application is the pull rate by the equipment that lifts in the control monocrystalline silicon growing process only, so that crystal diameter remains on the desired value, and has ignored in crystal growing process melt temperature to the influence of the diameter of crystal growth.Simultaneously, the dynamic property that incomplete differential PID handles or smith's method is handled of this application employing is relatively poor.
Chinese patent application 00805352.9 discloses the method and apparatus of control diameter of silicon crystal in the growth process.It is a kind of along with promoting the method and apparatus of the temperature of crystal bar by changing melt with the control boule diameter from melt that this application provides, and its crystal bar draws high from melt with the targeted rate according to predetermined velocity profile basically.Obviously this application is a controlled member with the temperature of melt only, organically combines and the temperature control of melt is drawn high speed control with crystal bar, and the dynamic property handled of PID is relatively poor simultaneously, so this application is more limited to the control ratio of crystal diameter growth.
Usually, the person of ordinary skill in the field only considers by the method for PID control melt temperature in pull rate control, the crystal growing process to be controlled, thereby cause the Controlling System of existing isodiametric growth of crystal and method thereof to be difficult to further improve in the crystal growing process to the control accuracy of diameter, be difficult to reach the production purpose of the isodiametric growth of crystal of expection.
Summary of the invention
The object of the present invention is to provide a kind of with Fuzzy (forecasting type is fuzzy) control and two kinds of algorithm combinations of PID control, improve the diameter control accuracy in isodiametric growth of crystal stage, thus the Controlling System and the method thereof of the isodiametric growth of crystal of the yield rate that stable raising large diameter silicon monocrystal rod once draws.
The technical scheme that realizes the object of the invention is: a kind of Controlling System of isodiametric growth of crystal has
Crystal growing furnace;
Detect the crystal diameter test section of crystal diameter variable signal in the interior crystal growing process of stove;
Gather melt temperature signal and the heater voltage current signal of well heater and the crucible temperature control part of the interior melt temperature of feedback control stove in the stove;
Gathering crystal draws high the crystal that speed signal and feedback adjustment crystal draw high speed in real time and draws high the speed control part;
It is characterized in that:
Also has a shared control unit; Described controller receives respectively from described crystal diameter test section, described crystal and draws high speed control part and described crucible temperature control part detection signal separately, described controller, use Fuzzy control and PID control bonded arithmetic calculation and go out crystal and draw high speed-regulating signal and feed back to described crystal and draw high the speed control part and draw province's speed to adjust, and calculate the heating power conditioning signal and feed back to described crucible temperature control part the temperature of melt is adjusted crystal.
In the technique scheme, described crucible temperature control partly comprises crucible temperature testing circuit, heater voltage current detection circuit and heating power regulating circuit; The output terminal of the crucible temperature detection signal of described crucible temperature testing circuit is electrically connected with the crucible temperature detection signal input terminus of described controller; The heater voltage current detection signal output terminal of described heater voltage current detection circuit is electrically connected with the heater voltage current detection signal input terminus of described controller; The heating power of described controller is regulated output terminal and is electrically connected with the heating power adjusting input terminus of described heating power regulating circuit.
In the technique scheme, described crystal draws high the speed controlling portion branch and comprises that crystal draws high tachometer circuit and crystal draws high alignment circuit; The crystal that described crystal draws high tachometer circuit draws high the tachometer signal output terminal and draws high the tachometer signal input terminus with the crystal of described controller and be electrically connected; The crystal of described controller draws high crystal that speed-regulating signal output terminal and described crystal draw high alignment circuit and draws high the speed-regulating signal input terminus and be electrically connected.
In the technique scheme, described crystal diameter test section comprises a crystal diameter testing circuit, and the crystal diameter detection signal output terminal of described crystal diameter testing circuit is electrically connected with the crystal diameter detection signal input terminus of controller.
In the technique scheme, controller is a programmable logic controller, and this programmable logic controller has simulating signal input block, CPU, simulating signal output unit, the I/O interface that connects by bus.
The Controlling System of isodiametric growth of crystal of the present invention, the method that adopts during its work has following steps:
1. measure and detect crystal diameter, gather the crystal diameter variable signal in the crystal growing process; Measure the detection crystal simultaneously and draw high speed, gather and draw high speed signal in real time; Measure the temperature that detects melt in the crucible, gather the melt temperature signal; And the heater voltage current signal of measuring detection heater heats electric current and voltage size; Described crystal diameter variable signal, described speed signal, described melt temperature signal and the described heater voltage current signal of drawing high in real time are sent on the shared control unit;
2. described controller adopts position model PID (proportional integral derivative) that the described crystal diameter variable signal and the described speed signal that draws high in real time that receive are handled, and obtains to draw high speed control signal; Simultaneously described controller adopts increment type PID that the melt temperature signal and the heater voltage current signal that receive are handled, and obtains the heating power conditioning signal;
3. the speed control signal that draws high that will draw after will handling is used for control and adjusts crystal and draw high speed; And the heating power conditioning signal that will draw after will handling is used for control adjustment heater power.
The step of technique scheme 2. in, controller drawing high the arithmetic expression that is adopted when alignment circuit is controlled growth in thickness by crystal is: SL=SP+F (Δ D)
In the formula, SL is that crystal draws high speed, and SP is that given crystal draws high speed, and F (Δ D) is diameter control PID output, and Δ D is a diameter deviation; The arithmetic expression to electrical signal that is adopted in the 3. middle temperature controlled processes of step is:
F ( ΔT ) = F [ Tset - Ti + Σ i = 0 n Tr * Δt ] + Fset
In the formula, F (Δ T) is temperature control PID output, and Fset is predetermined control PID output, and Tset is that temperature is set, and Ti is a crucible temperature detection signal,
Figure A20051005932500072
Compensation for warm school, Δ t is the time segment length;
Wherein:
Tr = Σ i = 0 n Trset * ΔL + F ( ΔV )
In the formula, Tr is warm school slope (μ V/min), and Trset is for being warm school slope preset value, and Δ L is that crystal draws high the segment length, and F (Δ V) is warm school control PID output, and Δ V is that crystal draws high real speed and given speed deviate.
The step of technique scheme 2. in, controller is in drawing high the process that alignment circuit controls growth in thickness by crystal, for suppressing high frequency interference and sampling action to spike interferential susceptibility, change desirable differential into restricted differential, deviation signal is through medium filtering; The formula of this process is as follows:
m n = 1 p { e n + 1 T I Σ i = 0 n e i ΔT + T D ΔT + T D K D ( e n - e n - 1 ) + T D K D ΔT + T D K D m n - 1 }
The step of technique scheme 3. in temperature control, for avoiding causing the vibration of power control circuit, realize manually control and the noiseless switching of control automatically; Carry out in the treating processes at increment type PID simultaneously, differential forward, promptly differential algorithm is not influenced by set(ting)value, only the deviation of measuring value signal is carried out differential algorithm; The formula of this process is as follows:
Δ m n = m n - m n - 1 = 1 p { ΔT T I e n + ( e n - e n - 1 ) + T D ΔT ( e ‾ n - 2 e ‾ n - 1 + e ‾ n - 2 ) }
In following formula, m nBe the n time output valve, e nBe the n time deviation signal, P is a proportional band, and Δ T is the sampling period, T IBe integration time constant, T DBe derivative time constant, K DBe the differential gain; Under the situation of large deviation, controller adopts the fuzzy processing mode from the pid parameter control of adjusting of forecasting type; The formula of this process is as follows:
m n = K P e n + K I Σ i = 0 n e i + K D e ‾ n
In the formula, K P, K IK DParameter for the Fuzzy sets theory; | E| is an absolute value of the bias, | EC| is the deviation change absolute value, according to different | E| and | EC| online self-tuning K P, K I, K DWhen | when E| is big, get bigger K PWith less K D, and K IBe 0; When | E| hour, get bigger K PAnd K I, less K DBy with K P, K I, K D, | E| and | the EC| obfuscation, determine the degree of membership of Fuzzy subclass, design P1D parameter tuning table, by finding corresponding correction factor in the table, adopt critical proportional band law to determine base value, thereby adjust pid parameter.
In the technique scheme, the user is provided with the process curve of crystal growth control process according to sampled data and by industrial computer, and control data, working order, input state and the output state of controller is set simultaneously; When industrial computer operation monitoring configuration software, the user carries out picture supervision and control by the touch of display screen or the keyboard of industrial computer to crystal growing process, and the user can adjust to further setting of process curve and controlled variable simultaneously.
The positively effect that the present invention has: (1) the present invention regulates with melt temperature and embeds the adjusting of crystal growth pulling rate, realizes the diameter control accuracy of expection.(2) the present invention combines diameter control and the control of warm school, realizes the isometrical automatic growth of crystal diameter.(3) adopt intelligent controller control Heating temperature and crystalline isodiametric growth, to improve the quality of products, reduce working strength, energy efficient has great meaning.(4) the present invention combines Fuzzy (bluring) control and two kinds of algorithms of PID (proportional integral derivative) control, has obtained all better controlled characteristics of dynamic property, steady-state behaviour.Consider PID and Fuzzy control algolithm relative merits separately, adopt Fuzzy control in the large deviation scope, the method for PID control in the little deviation range.Promptly, control the steady-state error of the system of elimination with PID, thereby make system reach good controlling performance with the good dynamic response of Fuzzy control assurance system.
Description of drawings
Fig. 1 is the circuit assembly block diagram of the Controlling System of isodiametric growth of crystal of the present invention.
Fig. 2 is the diameter control principle block diagram of controller of the present invention.
Fig. 3 is the temperature control principle block diagram of controller of the present invention.
Fig. 4 for the present invention on the basis of conventional PID control, use the Fuzzy sets theory, according to different | E| and | EC| online self-tuning K D, K I, K DThe time, the model block diagram of adjusting that is adopted.
Fig. 5 is the system architecture diagram of industrial computer 8 of the present invention.
Fig. 6 is the system architecture diagram of controller 7 of the present invention.
Embodiment
Referring to Fig. 1, the Controlling System of the isodiametric growth of crystal of present embodiment, have crystal growing furnace, crystal diameter testing circuit 1, crystal draw high tachometer circuit 2, crucible temperature testing circuit 3, heater voltage current detection circuit 4, crystal and draw high alignment circuit 5, heating power regulating circuit 6, controller 7, industrial computer 8 and display screen 9.
The crystal diameter detection signal output terminal of crystal diameter testing circuit 1 is electrically connected with the crystal diameter detection signal input terminus of controller 7, the crystal that crystal draws high tachometer circuit 2 draws high the tachometer signal output terminal and draws high the tachometer signal input terminus with the crystal of controller 7 and be electrically connected, the output terminal of crucible temperature testing circuit 3 crucible temperature detection signals is electrically connected with the crucible temperature detection signal input terminus of controller 7, the heater voltage current detection signal output terminal of heater voltage current detection circuit 4 is electrically connected with the heater voltage current detection signal input terminus of controller 7, the crystal of controller 7 draws high crystal that speed-regulating signal output terminal and crystal draw high alignment circuit 5 and draws high the speed-regulating signal input terminus and be electrically connected, and the heating power of controller 7 is regulated output terminal and is electrically connected with the heating power adjusting input terminus of heating power regulating circuit 6.
Industrial computer 8 is the industry control PC, and the screen of display screen 9 is tangible control; The two-way electrical connection of RS-232 serial data end of the RS-232 serial data end of controller 7 and industrial computer 8.
When controller 7 and industrial computer 8 communications, baud rate is 9600b/s, and the employing cycle is 1s, and controller 7 adopts interrupt mode when receiving data, sends and adopts the routine call launching means.The demonstration signal output part of industrial computer 8 is electrically connected with the demonstration signal input part of display screen 9, and the touch signal output terminal of display screen 9 is electrically connected with the touch signal input terminus of industrial computer 8.
Controller 7 is a programmable logic controller, and this programmable logic controller has: 8 simulating signal input blocks 71, its model are AD003; CPU72, model is Omron PLC C200HE/CPU42; 8 simulating signal output units 71, its model are DA003; 2 16 I/O input blocks 74, its model is 1D212; 2 16 I/O output units 75, its model is OC225; Simulating signal input block 71, CPU72, simulating signal output unit 71, I/O input block 74 and I/O output unit 75 all have EBI, and simulating signal input block 71, CPU72, simulating signal output unit 71, I/O input block 74 and I/O output unit 75 are connected electrically on the system bus mutually.
The method that adopts when the Controlling System of the isodiametric growth of crystal of present embodiment is worked has following steps: 1. measure and detect crystal diameter, gather the crystal diameter variable signal in the crystal growing process; Measure the detection crystal simultaneously and draw high speed, gather and draw high speed signal in real time; Measure the temperature that detects melt in the crucible, gather the melt temperature signal; And the heater voltage current signal of measuring detection heater heats electric current and voltage size; Described crystal diameter variable signal, described speed signal, described melt temperature signal and the described heater voltage current signal of drawing high in real time are sent on the shared control unit (7);
2. described controller (7) adopts position model PID (proportional integral derivative) that the described crystal diameter variable signal and the described speed signal that draws high in real time that receive are handled, and obtains to draw high speed control signal; Described controller (7) adopts increment type PID that the melt temperature signal and the heater voltage current signal that receive are handled simultaneously, obtains the heating power conditioning signal;
Step 2. in, controller 7 drawing high the arithmetic expression that is adopted when 5 pairs of growths in thickness of alignment circuit are controlled by crystal is: SL=SP+F (Δ D)
In the formula, SL is that crystal draws high speed, and SP is that given crystal draws high speed, and F (Δ D) is diameter control PID output, and Δ D is a diameter deviation;
3. the speed control signal that draws high that will draw after will handling is used for control and adjusts crystal and draw high speed; And the heating power conditioning signal that will draw after will handling is used for control adjustment heater power.
The arithmetic expression to electrical signal that is wherein adopted in the temperature controlled processes is:
F ( ΔT ) = F [ Tset - Ti + Σ i = 0 n Tr * Δt ] + Fset
In the formula, F (Δ T) is temperature control PID output, and Fset is predetermined control PID output, and Tset is that temperature is set, and Ti is a crucible temperature detection signal,
Figure A20051005932500102
Compensation for warm school, Δ t is the time segment length;
Wherein:
Tr = Σ i = 0 n Trset * ΔL + F ( ΔV )
In the formula, Tr is warm school slope μ V/min, and Trset is for being warm school slope preset value, and Δ L is that crystal draws high segment length mm, and F (Δ V) is warm school control PID output, and Δ V is that crystal draws high real speed and given speed deviate mm.
See Fig. 5 and Fig. 6, the user is provided with the process curve of crystal growth control process according to sampled data and by industrial computer 8, and control data, working order, input state and the output state of controller 7 is set simultaneously; When industrial computer 8 operation monitoring configuration softwares, the user carries out picture supervision and control by the touch of display screen 9 or the keyboard of industrial computer 8 to crystal growing process, and the user can adjust to further setting of process curve and controlled variable simultaneously.
Because the complicacy of controlled member is difficult to mathematical model and accurately expresses, and can only do qualitative analysis.From control principle, can determine the basic controlling strategy.As everyone knows, the dynamic property of PID control is relatively poor, but its integrating function can be eliminated static difference, and steady-state behaviour is improved; The dynamic property of Fuzzy control is better, but because quantification gradation limits, makes the stable state accuracy can not be satisfactory.Therefore, Fuzzy control and two kinds of algorithms of PID control are combined, just might obtain all better controlled characteristics of dynamic property, steady-state behaviour.Consider PID and Fuzzy control algolithm relative merits separately, adopt Fuzzy control in the large deviation scope, the method for PID control in the little deviation range.Promptly, control the steady-state error of the system of elimination with PID, thereby make system reach good control with the good dynamic response of Fuzzy control assurance system; Quality.
In step diameter control 2., adopt position model PID.For suppressing high frequency interference and sampling action to spike interferential susceptibility, change desirable differential into restricted differential, deviation signal is through medium filtering.Formula is as follows:
m n = 1 p { e n + 1 T I Σ i = 0 n e i ΔT + T D ΔT + T D K D ( e n - e n - 1 ) + T D K D ΔT + T D K D m n - 1 }
Step 3. in temperature control, for avoiding causing the vibration of power control circuit, realize manually and automatic noiseless switching, adopt increment type PID, and differential forward, promptly differential algorithm is not influenced by set(ting)value, only the deviation of measuring value signal is carried out differential algorithm.Formula is as follows:
Δ m n = m n - m n - 1 = 1 p { ΔT T I e n + ( e n - e n - 1 ) + T D ΔT ( e ‾ n - 2 e ‾ n - 1 + e ‾ n - 2 ) }
In the formula, m nBe the n time output valve; e nBe the n time deviation signal; P is a proportional band; Δ T is the sampling period; T IBe integration time constant; T DBe derivative time constant; K DBe the differential gain.
Under the situation of large deviation, adopt forecasting type Fuzzy from Adjustable PID Parameter Controller 7.Formula is as follows:
m n = K P e n + K I Σ i = 0 n e i + K D e ‾ n
On the basis of conventional PID control, use the Fuzzy sets theory, set up parameter K P, K I, K DSame absolute value of the bias | E| and deviation change absolute value | binary continuous function between EC| relation and based on the fuzzy rule of thumb rule, and according to different | E| and | EC| online self-tuning K P, K I, K DAdjust model as shown in Figure 4.
When | when E| is big, get bigger K PWith less K D, and K IBe 0; When | E| hour, get bigger K PAnd K I, less K DBy with K P, K I, K D, | E| and | the EC| obfuscation, determine the degree of membership of Fuzzy subclass, design P1D parameter tuning table.By finding corresponding correction factor in the table, adopt critical proportional band law to determine base value, thereby adjust pid parameter.
In the present embodiment, because it is that the ring of light signal during to crystal growth is measured indirectly by infrared thermometer that crystal diameter changes, through this ring of light signal is amplified, computer is sampled again afterwards, so this ring of light pickup electrode is vulnerable to external interference, introduce the velocity feedback of feed forward control for this reason, this ring of light signal is carried out dynamic compensation, suppressed external interference effectively.
In temperature control, introduce nonlinear Control, guarantee in the limit of error of predesignating, by crystal being drawn high the measurement of speed, the pulling rate setting curve that crystal is drawn high speed and process of growth compares, the crystalline growth temperature is regulated control, the crystal pull rate is changed by the technique initialization curve, reduced because PID regulates the hysteresis that produces.
See Fig. 2, in the present embodiment, when crystal diameter changed, controller 7 can be regulated crystal automatically and draw high speed, and crystal diameter is remained unchanged.
1, when diameter is controlled, the explanation of the input parameter of controller 7:
1. diameter signal D IBe variation, thereby measure the crystalline vary in diameter indirectly by the ring of light signal of measuring the solid-liquid interface in the crystal growing process.When crystal diameter becomes big or diminishes, diameter signal D IAlso become big or diminish.
2. diameter is set SV: when isodiametric growth of crystal, when controller 7 inputs were automatic, the diameter signal of getting this moment was as the diameter set(ting)value.
3. control speed output OP: be a bipolarity output parameter, the positive and negative variation of its numerical value, dimension is mm/min, to draw high the speed measure unit consistent with real crystal.
When crystal diameter becomes big, diameter signal D IBecome big, control fast output OP then and become big, crystal draws high speed SL change greatly then, and last crystal diameter diminishes.
When crystal diameter diminishes, diameter signal D IDiminish, control speed output OP diminishes then, and crystal draws high speed SL and diminishes then, and last crystal diameter becomes big.
No matter how diameter signal changes, controller 7 can both guarantee the variation of crystal diameter in certain scope, the effect and the purpose of Here it is controller 7.
2, when diameter is controlled, the inner parameter of controller 7:
1. proportional gain P n: (P n=0~1.99).P nShow that proportional action is strong more, P greatly more nMore for a short time, show that proportional action is more little.Work as P n=0 o'clock, proportional action was invalid.
2. integration constant I n: (I n=0~199).I nBig more, show that integral action is more little, I nMore little, show that integral action is big more.Work as I n>199 o'clock, integral action was invalid.
3. differential gain D n: (D n=0~9.9).D nBig more, show that the differential action is strong more, D nMore little, show that the differential action is weak more.
4. control cycle t (t=0.1~99.9 second).Control cycle be controller 7 every a control cycle time, once adjust controlling fast output valve.
5. control diameter parameter Md (Md=0.00~3.00).Change according to crystal diameter and crystal to be drawn high speed control, the value of E1 is chosen in E1=1.0.
6. differential slope Dk (Dk=00~99).The value of general E2 gets final product between 50~80.
7. the selection of controller 7 parameters.P in the present embodiment n=1.0, I n=30, D n=0.8.Md=1.20,Dk=80。
See Fig. 3, the effect of warm school control is in crystal growing process, and the control Heating temperature makes crystal draw high speed and changes by the setting of pulling rate setting curve SL ramp.When crystal draws high the speed deviations set(ting)value, change the size of temperature control output OP value, crystal diameter is changed, at this moment, the speed that controller 7 control crystal draw high makes crystal draw high speed SL near pulling rate set(ting)value SP.
When the control of warm school, the parameter declaration of controller 7:
1, I/O parameter
The I/O parameter comprises: crystal draws high speed SL, temperature control output OP, warm school speed t r, set pulling rate SP (being pulling rate setting curve SL ramp).
1. crystal draws high speed SL (0~9.99mm/min).It is the same that the displayed value allomeric of this parameter draws high the demonstration of speed unit, and the parameter dimension is mm/min.
2. temperature control is exported OP (0~4095).This parameter is an output digital parameters, and it exports numerical value 0~4095, and corresponding output voltage is (adjustable) about 0~12mv.
3. warm school speed t r, (99.9~99.9 μ V/MIN).Temperature school speed is meant the variable quantity of per minute temperature control output OP, works as t r<0 is cooling state, t r>0 is soaking condition, t r=0 is temperature constant state.
4. set pulling rate SP (SL ramp).Drop into when automatic at controller 7, the crystal in this is moment drawn high speed SL as setting pulling rate SP.The variation with crystal length changes by the parameter setting of pulling rate setting curve SL ramp to set pulling rate simultaneously.
2, when the control of warm school, the inner parameter of controller 7:
1. proportional gain P n(P n=0~5.00).P nBig more, proportional action is strong more, P oMore little, proportional action is weak more.
2. integration constant I n(I n=0~500).I nBig more, integral action is strong more, I nMore little, integral action is weak more.
3. control cycle: t (t=0.1~99.9 second).Control cycle shows every time t, and warm school speed is once adjusted.
4. the selection of controller 7 inner parameters.In the ordinary course of things, establish P earlier n=1.0, I n=20, t=20.0.
Principle and characteristics about PID control:
In engineering reality, the setter control law that is most widely used is ratio, integration, differential control, is called for short PID control, claims PID to regulate again.The PID controller existing so far nearly 70 years history of coming out, it is simple in structure with it, good stability, reliable operation, easy to adjust and become one of major technique of Industry Control.When the structure and parameter of controlled plant can not be on top of, or when can not get precise math model, when other technology of control theory was difficult to adopt, the structure and parameter of central controller must rely on experience and field adjustable to determine, at this moment the thereof using PID control techniques is the most convenient.Promptly understand a system and controlled plant when we are incomplete, in the time of maybe can not obtaining system parameter by effective measurement means, the most suitable PID control techniques of using.PID control also has PI and PD control in the reality.The PID controller is exactly the error according to system, and proportion of utilization, integration, difference gauge calculate that manipulated variable controls.
The output and the proportional relation of error originated from input signal of the controller of ratio P control.There is steady-state error in system's output when proportional control only.
In integration I control, the output of controller and error originated from input signal be integrated into proportional relation.To an automatic control system, if after entering stable state, there is steady-state error, then claim this Controlling System be steady-state error arranged or be called for short poor system arranged.In order to eliminate steady-state error, in controller, must introduce " integral ".Integral depends on the integration of time to error, and along with the increase of time, integral can increase.Like this, even if error is very little, integral also can strengthen along with the increase of time, and the output increase that it promotes controller further reduces steady-state error, up to equalling zero.Therefore, the PI controller of " proportional+integral " can make system not have steady-state error after entering stable state.
In differential D control, the differential of the output of controller and error originated from input signal is that the error change rate is proportional.Vibration even unstability may appear in automatic control system in overcoming the regulate process of error.Its reason is owing to have big inertia assembly link or hysteresis delay assembly is arranged, and has the effect that suppresses error, and its variation always lags behind error change.The way that solves is the variation " in advance " that makes the effect that suppresses error, promptly in error near zero the time, the effect that suppresses error just should be zero.In other words, it is not enough often only to introduce " ratio " item in controller, the effect of proportional only is the amplitude of fault in enlargement, and what need at present to increase is " differential term ", the trend that its can predicated error changes, like this, PD controller with " ratio+differential ", the effect of inhibition error Control is equalled zero, even be negative value, thereby avoid the serious overshoot of controlled volume.So to the controlled plant of big inertia or hysteresis is arranged, the PD controller can the kinetic characteristic of improvement system in regulate process.
The parameter tuning of PID controller is the core content of Control System Design.The parameter tuning of PID controller is a size of determining scale-up factor, integral time and the differential time of PID controller according to the characteristic of controlled process.The method that the PID controller parameter is adjusted is a lot, summarizes and to get up to have two big classes: the one, and the Theoretical Calculation method of adjusting.It mainly is the mathematical model according to system, determines controller parameter through Theoretical Calculation.The resulting data calculated of this method may not directly be used, and also must adjust and revise by engineering is actual.The 2nd, the engineering setting method, it mainly relies on engineering experience, directly in the test of Controlling System, carry out, and method simple, be easy to grasp, in engineering reality, be widely adopted.The engineering setting method of PID controller parameter mainly contains aritical ratio method, response curve method and damped method, and three kinds of methods respectively have its characteristics, and its common ground all is by test, according to the engineering experience formula controller parameter is adjusted then.But no matter adopt the resulting controller parameter of any method, all need in actual motion, adjust at last with perfect.General employing is the aritical ratio method now.It is as follows to utilize this method to carry out the step of adjusting of PID controller parameter: (1) at first preselected enough short sampling period allows system works; (2) additional proportion controlling unit only occurs threshold oscillation up to system to the step response of input, writes down at this moment ratio scale-up factor and threshold oscillation cycle; (3) under certain degree of control, calculate the parameter of PID controller by formula.

Claims (10)

1, a kind of Controlling System of isodiametric growth of crystal has
Crystal growing furnace;
Detect the crystal diameter test section of crystal diameter variable signal in the interior crystal growing process of stove;
Gather melt temperature signal and the heater voltage current signal of well heater and the crucible temperature control part of the interior melt temperature of feedback control stove in the stove;
Gathering crystal draws high the crystal that speed signal and feedback adjustment crystal draw high speed in real time and draws high the speed control part;
It is characterized in that:
Also has a shared control unit (7); Described controller (7) receives respectively from described crystal diameter test section, described crystal and draws high speed control part and described crucible temperature control part detection signal separately, described controller (7) is used Fuzzy control and PID control bonded arithmetic calculation and is gone out crystal and draw high speed-regulating signal and feed back to described crystal and draw high the speed control part and draw province's speed to adjust to crystal, and calculates the heating power conditioning signal and feed back to described crucible temperature control part the temperature of melt is adjusted.
2, the Controlling System of a kind of isodiametric growth of crystal according to claim 1 is characterized in that: described crucible temperature control partly comprises crucible temperature testing circuit (3), heater voltage current detection circuit (4) and heating power regulating circuit (6); The output terminal of the crucible temperature detection signal of described crucible temperature testing circuit (3) is electrically connected with the crucible temperature detection signal input terminus of described controller (7); The heater voltage current detection signal output terminal of described heater voltage current detection circuit (4) is electrically connected with the heater voltage current detection signal input terminus of described controller (7); The heating power of described controller (7) is regulated output terminal and is electrically connected with the heating power adjusting input terminus of described heating power regulating circuit (6).
3, the Controlling System of a kind of isodiametric growth of crystal according to claim 1 is characterized in that: described crystal draws high the speed controlling portion branch and comprises that crystal draws high tachometer circuit (2) and crystal draws high alignment circuit (5); The crystal that described crystal draws high tachometer circuit (2) draws high the tachometer signal output terminal and draws high the tachometer signal input terminus with the crystal of described controller (7) and be electrically connected; The crystal of described controller (7) draws high crystal that speed-regulating signal output terminal and described crystal draw high alignment circuit (5) and draws high the speed-regulating signal input terminus and be electrically connected.
4, the Controlling System of a kind of isodiametric growth of crystal according to claim 1, it is characterized in that: described crystal diameter test section comprises a crystal diameter testing circuit (1), and the crystal diameter detection signal output terminal of described crystal diameter testing circuit (1) is electrically connected with the crystal diameter detection signal input terminus of controller (7).
5, the Controlling System of a kind of isodiametric growth of crystal according to claim 1, it is characterized in that: controller (7) is a programmable logic controller, and this programmable logic controller has simulating signal input block (71), CPU (72), simulating signal output unit (71), the I/O interface that connects by bus.
6, a kind of control method that is used for the Controlling System of the described isodiametric growth of crystal of claim 1 to 5 has following steps:
1. measure and detect crystal diameter, gather the crystal diameter variable signal in the crystal growing process; Measure the detection crystal simultaneously and draw high speed, gather and draw high speed signal in real time; Measure the temperature that detects melt in the crucible, gather the melt temperature signal; And the heater voltage current signal of measuring detection heater heats electric current and voltage size; Described crystal diameter variable signal, described speed signal, described melt temperature signal and the described heater voltage current signal of drawing high in real time are sent on the shared control unit (7);
2. described controller (7) adopts position model PID (proportional integral derivative) that the described crystal diameter variable signal and the described speed signal that draws high in real time that receive are handled, and obtains to draw high speed control signal; Described controller (7) adopts increment type PID that the melt temperature signal and the heater voltage current signal that receive are handled simultaneously, obtains the heating power conditioning signal;
3. the speed control signal that draws high that will draw after will handling is used for control and adjusts crystal and draw high speed; And the heating power conditioning signal that will draw after will handling is used for control adjustment heater power.
7, method according to claim 6 is characterized in that: step 2. in, controller (7) is when drawing high alignment circuit (5) growth in thickness is controlled to crystalline by crystal, the arithmetic expression that is adopted is: SL=SP+F (Δ D)
In the formula, SL is that crystal draws high speed, and SP is that given crystal draws high speed, and F (Δ D) is diameter control PID output, and Δ D is a diameter deviation; The arithmetic expression to electrical signal that is adopted in the 3. middle temperature controlled processes of step is:
F ( ΔT ) = F [ Tset - Ti + Σ i = 0 n Tr * Δt ] + Fset
In the formula, F (Δ T) is temperature control PID output, and Fset is predetermined control PID output, and Tset is that temperature is set, and Ti is a crucible temperature detection signal, Compensation for warm school, Δ t is the time segment length;
Wherein: Tr = Σ i = 0 n Trset * ΔL + F ( ΔV )
In the formula, Tr is warm school slope (μ V/min), and Trset is for being warm school slope preset value, and Δ L is that crystal draws high segment length (mm), and F (Δ V) is warm school control PID output, and Δ V is that crystal draws high real speed and given speed deviate (mm).
The method that is adopted when 8, the Controlling System of a kind of isodiametric growth of crystal according to claim 6 is worked, it is characterized in that: step 2. in, controller (7) is in drawing high the process that alignment circuit (5) controls growth in thickness by crystal, adopt restricted differential, to deviation signal process medium filtering; The formula of this process is as follows:
m n = 1 p { e n + 1 T I Σ i = 0 n e i ΔT + T D ΔT + T D K D ( e n - e n - 1 ) + T D K D ΔT + T D K D m n - 1 }
9, method according to claim 6 is characterized in that; In step temperature control 2., carry out in the treating processes at increment type PID, the while differential forward, promptly differential algorithm is not influenced by set(ting)value, only the deviation of measuring value signal is carried out differential algorithm; The formula of this process is as follows:
Δm n = m n - m n - 1 = 1 p { ΔT T I e n + ( e n - e n - 1 ) + T D ΔT ( e ‾ n - 2 e ‾ n - 1 + e ‾ n - 2 ) }
In following formula, m nBe the n time output valve, e nBe the n time deviation signal, P is a proportional band, and Δ T is the sampling period, T IBe integration time constant, T DBe derivative time constant, K DBe the differential gain; Under the situation of large deviation, controller (7) adopts adjust the certainly processing mode of pid parameter control of forecasting type fuzzy (Fuzzy); The formula of this process is as follows:
m n = K P e n + K I Σ i = 0 n e i + K D e ‾ n
In the formula, K P, K IK DParameter for the Fuzzy sets theory; | E| is an absolute value of the bias, | EC| is the deviation change absolute value, according to different | E| and | EC| online self-tuning K P, K I, K DWhen | when E| is big, get bigger K PWith less K D, and K IBe 0; When | E| hour, get bigger K PAnd K I, less K DBy with K P, K I, K D, | E| and | the EC| obfuscation, determine the degree of membership of Fuzzy subclass, design the pid parameter table of adjusting, by finding corresponding correction factor in the table, the employing critical proportional band law is determined base value, thereby adjusts pid parameter.
10, method according to claim 6, it is characterized in that: the user is provided with the process curve of crystal growth control process according to sampled data to described controller input, and control data, working order, input state and the output state of controller (7) are set; When industrial computer (8) operation monitoring configuration software, user's output, input are carried out picture supervision and control to crystal growing process; According to processing requirement, the user can further set with controlled variable process curve and adjust.
CNB2005100593258A 2005-03-28 2005-03-28 Control system and method for crystal isodiametric growth Expired - Fee Related CN100371507C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005100593258A CN100371507C (en) 2005-03-28 2005-03-28 Control system and method for crystal isodiametric growth

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005100593258A CN100371507C (en) 2005-03-28 2005-03-28 Control system and method for crystal isodiametric growth

Publications (2)

Publication Number Publication Date
CN1840746A true CN1840746A (en) 2006-10-04
CN100371507C CN100371507C (en) 2008-02-27

Family

ID=37029930

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100593258A Expired - Fee Related CN100371507C (en) 2005-03-28 2005-03-28 Control system and method for crystal isodiametric growth

Country Status (1)

Country Link
CN (1) CN100371507C (en)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100560810C (en) * 2007-09-18 2009-11-18 山东大学 Integration type programmable crystal growth control system
CN102534757A (en) * 2012-01-06 2012-07-04 南京师范大学 Intelligent crystal growth control system based on industrial bus
CN102758250A (en) * 2012-07-20 2012-10-31 西安理工晶体科技有限公司 Automatic equal diameter controlling method for germanium single crystal straight-pull growth method
CN103261492A (en) * 2010-12-13 2013-08-21 Lg矽得荣株式会社 System of controlling diameter of single crystal ingot and single crystal ingot growing apparatus including the same
CN103834991A (en) * 2014-03-10 2014-06-04 马鞍山明鑫电气科技有限公司 Open-loop power self-control crystal growth control method without temperature signal treatment
CN104514030A (en) * 2013-09-29 2015-04-15 上海昀丰新能源科技有限公司 Crystal growth speed detection method, crystal growth speed control method and crystal growth speed control system
CN106245116A (en) * 2016-08-31 2016-12-21 天通银厦新材料有限公司 A kind of sapphire crystal growth equipment and technique thereof
CN106801250A (en) * 2017-01-13 2017-06-06 中山大学 The method of feedback transistor growth conditions, crystal growth control method and control system
CN107815729A (en) * 2016-09-12 2018-03-20 上海新昇半导体科技有限公司 A kind of single crystal growing furnace
CN109183141A (en) * 2018-10-29 2019-01-11 上海新昇半导体科技有限公司 A kind of crystal growth control method, device, system and computer storage medium
CN110004492A (en) * 2019-04-25 2019-07-12 苏州新美光纳米科技有限公司 Monitoring method and long crystal furnace in long crystal furnace
CN110528067A (en) * 2018-05-25 2019-12-03 隆基绿能科技股份有限公司 A kind of temprature control method of czochralski silicon monocrystal
CN111519244A (en) * 2020-05-28 2020-08-11 璨隆科技发展有限公司 Power supply control method of resistance heater applied to single crystal production equipment
CN112680789A (en) * 2019-12-13 2021-04-20 山东大学 HVPE temperature control system and method based on fuzzy control
CN113293433A (en) * 2020-02-21 2021-08-24 上海新昇半导体科技有限公司 Crystal bar growth control method and control system
CN113463184A (en) * 2021-06-29 2021-10-01 西门子(中国)有限公司 Method and device for controlling lifting speed in crystal growth process and readable medium
CN113584576A (en) * 2021-07-30 2021-11-02 上海众壹云计算科技有限公司 Method and apparatus for automatically controlling crystal diameter growth, electronic device, and storage medium
CN114318533A (en) * 2021-12-28 2022-04-12 安徽科瑞思创晶体材料有限责任公司 Intelligent control system for crystal growth

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI911857A (en) * 1990-04-27 1991-10-28 Nippon Kokan Kk FOERFARANDE OCH APPARAT FOER KONTROLL AV DIAMETERN HOS EN ENSKILD SILIKONKRISTALL.
US6171391B1 (en) * 1998-10-14 2001-01-09 Memc Electronic Materials, Inc. Method and system for controlling growth of a silicon crystal
US6776840B1 (en) * 1999-03-22 2004-08-17 Memc Electronic Materials, Inc. Method and apparatus for controlling diameter of a silicon crystal in a locked seed lift growth process
CN1396965A (en) * 2000-02-01 2003-02-12 Memc电子材料有限公司 Method for controlling growth of silicon crystal to minimize growth rate and diameter deviations
CN1393746A (en) * 2001-07-02 2003-01-29 李晓枫 Self-adaptive fuzzy controller for PID parameters

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100560810C (en) * 2007-09-18 2009-11-18 山东大学 Integration type programmable crystal growth control system
CN103261492A (en) * 2010-12-13 2013-08-21 Lg矽得荣株式会社 System of controlling diameter of single crystal ingot and single crystal ingot growing apparatus including the same
US9422637B2 (en) 2010-12-13 2016-08-23 Lg Siltron Inc. System of controlling diameter of single crystal ingot and single crystal ingot growing apparatus including the same
CN102534757A (en) * 2012-01-06 2012-07-04 南京师范大学 Intelligent crystal growth control system based on industrial bus
CN102758250A (en) * 2012-07-20 2012-10-31 西安理工晶体科技有限公司 Automatic equal diameter controlling method for germanium single crystal straight-pull growth method
CN104514030A (en) * 2013-09-29 2015-04-15 上海昀丰新能源科技有限公司 Crystal growth speed detection method, crystal growth speed control method and crystal growth speed control system
CN104514030B (en) * 2013-09-29 2017-01-04 内蒙古恒嘉晶体材料有限公司 Rate of crystalline growth detection method, control method and system
CN103834991A (en) * 2014-03-10 2014-06-04 马鞍山明鑫电气科技有限公司 Open-loop power self-control crystal growth control method without temperature signal treatment
CN106245116A (en) * 2016-08-31 2016-12-21 天通银厦新材料有限公司 A kind of sapphire crystal growth equipment and technique thereof
CN107815729A (en) * 2016-09-12 2018-03-20 上海新昇半导体科技有限公司 A kind of single crystal growing furnace
CN106801250A (en) * 2017-01-13 2017-06-06 中山大学 The method of feedback transistor growth conditions, crystal growth control method and control system
CN106801250B (en) * 2017-01-13 2019-05-07 中山大学 Method, crystal growth control method and the control system of feedback transistor growth conditions
CN110528067A (en) * 2018-05-25 2019-12-03 隆基绿能科技股份有限公司 A kind of temprature control method of czochralski silicon monocrystal
CN109183141A (en) * 2018-10-29 2019-01-11 上海新昇半导体科技有限公司 A kind of crystal growth control method, device, system and computer storage medium
CN110004492A (en) * 2019-04-25 2019-07-12 苏州新美光纳米科技有限公司 Monitoring method and long crystal furnace in long crystal furnace
CN112680789A (en) * 2019-12-13 2021-04-20 山东大学 HVPE temperature control system and method based on fuzzy control
CN113293433A (en) * 2020-02-21 2021-08-24 上海新昇半导体科技有限公司 Crystal bar growth control method and control system
CN113293433B (en) * 2020-02-21 2022-08-16 上海新昇半导体科技有限公司 Crystal bar growth control method and control system
CN111519244A (en) * 2020-05-28 2020-08-11 璨隆科技发展有限公司 Power supply control method of resistance heater applied to single crystal production equipment
CN111519244B (en) * 2020-05-28 2021-08-24 璨隆科技发展有限公司 Power supply control method of resistance heater applied to single crystal production equipment
CN113463184A (en) * 2021-06-29 2021-10-01 西门子(中国)有限公司 Method and device for controlling lifting speed in crystal growth process and readable medium
CN113584576A (en) * 2021-07-30 2021-11-02 上海众壹云计算科技有限公司 Method and apparatus for automatically controlling crystal diameter growth, electronic device, and storage medium
CN114318533A (en) * 2021-12-28 2022-04-12 安徽科瑞思创晶体材料有限责任公司 Intelligent control system for crystal growth

Also Published As

Publication number Publication date
CN100371507C (en) 2008-02-27

Similar Documents

Publication Publication Date Title
CN1840746A (en) Control system and method for crystal isodiametric growth
CN2825658Y (en) Control system for crystal isodiametric growth
CN105182740B (en) Raw material grinding autocontrol method
JP6325691B2 (en) PLC closed loop control method for sapphire single crystal growth
CN1344335A (en) Method and appts. for controlling diameter of silicon crystal in growth process
CN1212551C (en) Control equipment
CN1527170A (en) Sync control device
CN104753439B (en) A kind of PID intelligent speed-regulating methods of motor
CN104988577A (en) Sapphire automatic control system and control method
CN103834991A (en) Open-loop power self-control crystal growth control method without temperature signal treatment
CN102758250A (en) Automatic equal diameter controlling method for germanium single crystal straight-pull growth method
CN1082203C (en) Process controlling method and apparatus
CN114318533A (en) Intelligent control system for crystal growth
CN107037729A (en) A kind of design method based on RBF neural automatic disturbance rejection controller
CN1069134A (en) Program control system
CN1295632A (en) Open-loop method and system for controlling growth of semiconductor crystal
CN102351190A (en) Reducing furnace control method and device
CN115198351B (en) Heating power control method for automatically growing low-temperature sensitive artificial crystal material
CN208432894U (en) A kind of toughened glass production line control system
CN2913393Y (en) Control device for monocrystalline silicon furnace
CN106042327B (en) A kind of tire base width detection and transfer rate control method and detecting and controlling system
CN112301421B (en) Photosensitive rapid ending method suitable for large-size monocrystalline silicon
CN103558871A (en) Coating thickness automatic control system and method
CN114083086A (en) Control method and device for stabilizing electric arc additive manufacturing process
JP2001354490A (en) Method for producing single crystal

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20090206

Address after: Jiangsu city of Jintan Province Yao Tang Zhen Jin Hongwu Road No. 18

Patentee after: Changzhou EGing Photovoltaic Technology Co., Ltd.

Address before: No. 17 East Main Street, Tang Zhen Town, Jintan, Jiangsu

Patentee before: Xun Jianhua

ASS Succession or assignment of patent right

Owner name: CHANGZHOU YIJING OPTO-ELECTRICAL SCIENCE CO., LTD.

Free format text: FORMER OWNER: XUN JIANHUA

Effective date: 20090206

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080227

Termination date: 20150328

EXPY Termination of patent right or utility model