CN101562409B - Piezoelectric structure damping control object compensation method - Google Patents

Piezoelectric structure damping control object compensation method Download PDF

Info

Publication number
CN101562409B
CN101562409B CN2009100840041A CN200910084004A CN101562409B CN 101562409 B CN101562409 B CN 101562409B CN 2009100840041 A CN2009100840041 A CN 2009100840041A CN 200910084004 A CN200910084004 A CN 200910084004A CN 101562409 B CN101562409 B CN 101562409B
Authority
CN
China
Prior art keywords
piezoelectric
strain
moving plate
epsiv
measured sheet
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.)
Expired - Fee Related
Application number
CN2009100840041A
Other languages
Chinese (zh)
Other versions
CN101562409A (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN2009100840041A priority Critical patent/CN101562409B/en
Publication of CN101562409A publication Critical patent/CN101562409A/en
Application granted granted Critical
Publication of CN101562409B publication Critical patent/CN101562409B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position Or Direction (AREA)
  • Micromachines (AREA)

Abstract

The invention discloses a piezoelectric structure damping control object compensation method, which comprises the following: step 1, arranging a substrate; step 2, arranging an actuating piece; step 3, arranging a sensing vane; and step 4, adjusting gain g<c>. The damping control method provided by the invention overcomes the defects in the prior art, eliminates the feedback of partial excitation strain of the piezoelectric actuating piece so as to make a mathematic model of a piezoelectric structure return to a theoretical model, has skillful design and simple operation, and can be widely applied in the technical field of structural vibration active control.

Description

Piezoelectric structure damping control object compensation method
(1) technical field
The present invention relates to a kind of piezoelectric structure damping control method, relate in particular to a kind of piezoelectric structure damping control object compensation method, belong to the active control in structural vibration technical field.
(2) background technology
What the piezoelectric structure damping control of putting down in writing in numerous documents at present, was adopted is typical control model---the closed-loop system that negative velocity feedback constitutes.Explain typical control model principle and control method of the present invention with piezoelectric beam as typical structure among the present invention, but this method also is applicable to various other structures.
Fig. 1 is the typical control model sketch map that numerous documents adopt, and is visible from scheming:
1, the lower surface at beam 1 posts a piezoelectric actuating sheet 2, the upper surface same position post same size piezo-electric measurement sheet 3 and after connect operational amplifier 5 and constitute strain rate sensors with resistance 4; This actuator and transducer " coordination " configuration can guarantee to be connected into
2, the stability of the system after the negative velocity feedback.
2, piezoelectric strain rate transducer output u sThrough gain g sBack and external excitation voltage u eSubtract each other, through power amplifier gain g aBack excitation piezoelectric actuating sheet constitutes closed loop.
The shortcoming and defect that the Typical control method is brought has:
Damping ratio can be along with loop gain g in theory sg aAnd increase.But experimental result is really not so satisfactory.Along with g sg aIncrease, system but trends towards instability, produces high frequency or the very self-excited vibration of low frequency, and before this modal damping also do not reach expectation than increment requirement as 0.1.
The strain open loop frequency response that observation is surveyed, finding has significantly different with theoretical model:
1. phase-frequency characteristic: behind the pi/2 of resonance region whereabouts, but return to soon and approach zero.
2. amplitude-frequency characteristic: after crossing resonance region is not to decay down with the 40db/oct slope, but returns to an approximate constant always.
(3) summary of the invention
The purpose of this invention is to provide a kind of piezoelectric structure damping control object compensation method, it has overcome deficiency of the prior art, is ingenious, the simple to operate damping control method of a kind of design.
Existing that inventive principle is as described below:
1. foreword
Changing structural damping, always is one of basic assignment of vibration control, the topmost means of vibration suppression especially---increase damping.From traditional strigil and the Passive Control means such as distributed damping material of centralized damper to recent two decades, each means all are to utilize the energy dissipation behavior relevant with the basal body structure deformation velocity of accessory structure to reach the target that increases structural damping.Typical case's representative of classical ACTIVE CONTROL damping is the closed-loop system that constitutes with negative velocity feedback; In single-mode system; Extremely successful in theory and practice, to many-degrees of freedom system, the puzzlement of overflowing and controlling problems such as overflowing can appear observing; Possibly cause controlling deleterious, even loss of stability.
Another central issue of negative velocity feedback damping control is the equipment that what right sensors and actuator is arranged and match.For the large-scale flexible structure, more need large quantities of distributed transducers and actuator.The development of piezoelectric; People have just in time been adapted to for a long time to the requirement of this transducer and actuator hope; Thereby over nearly 20 years, control about negative velocity feedback structure active damping with piezoelectric patches actuator and transducer---the research mushroom development of so-called " electronic damping " is got up.
At first, on the basis of 2 joint piezoelectric structure modelings,, find that the result can not show a candle to the so good of expection so that the local excitation of piezoelectric cantilever---the damping of sensing negative velocity feedback is controlled to be example and launches to discuss.Particularly piezoelectric actuator---the Frequency Response Analysis and the computer simulation experiment of transducer are probed into its reason to experimental phenomena, are " partial excitation strain " interference that this special piezoelectricity " strain excitation---strain-responsive " of piezoelectric structure brings.This has caused our self-examination to 2 joint piezoelectric structure modelings, has proposed qualitatively but enough corrections of practicality at 3.2 joints.Then, 3.3 the joint, the partial excitation strain compensation scheme has been proposed---let the piezoelectric structure Mathematical Modeling revert to 2 the joint theoretical model in go, the damping that the result controls in experiment is greatly improved.
2. the transfer function model of piezoelectric structure
On a structure, paste some piezoelectric patches, utilize the vibrational state of their direct piezoelectric effect measurement structure, call " piezoelectric transducer " to these piezoelectric patches and auxiliary equipment thereof; Paste some piezoelectric patches in addition again and use piezoelectric excitation, its reversed piezoelcetric effect will evoke structural vibration, calls " piezoelectric actuator " to these piezoelectric patches and auxiliary equipment thereof.The mechanical-electric coupling structure that basal body structure, piezoelectric actuator and piezoelectric transducer lump together is called " piezoelectric structure ".
2.1 four types of piezoelectric equations of piezoelectric ceramic piece
The machine of piezoelectrics---electric constitutive relation also is that piezoelectric equations is described the coupled relation between its mechanical quantity (stress T and strain S) and electrical quantities (electric field strength E and electric displacement D).This joint will be derived four types of piezoelectric equations of the piezoelectric ceramic piece of using in theoretical and the experimental study (belonging to 6mm point group crystal).
If piezoelectric sheet plane internal coordinate is x, y, normal direction also is that polarised direction is z.To be main with the symbol of habitually practising in the mechanical analysis herein, but analyze in the argumentation that piezoelectric patches is the master at Ben Jie and later on relevant chapters and sections, be the original appearance that keeps the piezoelectrics analysis, analyzes the general symbol in the document with continuing to use piezoelectrics, and the two table of comparisons is following:
Axle Stress Strain
From the mechanics angle x?y?z σ xyzyzzxxy ε xyzyzzxxy
From the piezoelectricity angle 1?2?3 T 1?T 2?T 3?T 4?T 5?T 6 S 1?S 2?S 3?S 4?S 5?S 6
Piezoelectric ceramic has 2 kinds of piezoelectric effect equations the most basic as 6mm point group crystal; One of which, electrical short (E 1=E 2=E 3=0) direct piezoelectric effect under the condition
D 1 D 2 D 3 = 0 0 0 0 d 15 0 0 0 0 d 25 0 0 d 31 d 32 d 33 0 0 0 T 1 T 2 T 3 T 4 T 5 T 6 , ( E 1 = E 2 = E 3 = 0 ) - - - ( 2 - 1 )
Its two, the free (T of mechanics i=0, the i=1-6) inverse piezoelectric effect under the condition
S 1 S 2 S 3 S 4 S 5 S 6 = 0 0 d 31 0 0 d 32 0 0 d 33 0 d 25 0 d 15 0 0 0 0 0 E 1 E 2 E 3 , ( T i = 0 , i = 1 - 6 ) - - - ( 2 - 2 )
E wherein iAnd D iBe respectively inside field intensity and the pole-face electric displacement of piezoelectrics, d along the i axle IjPiezoelectric constant for relevant i and j axle has
d p=d 31=d 32 (2-3)
Wherein subscript p indicates " piezoelectric patches ", and this paper will continue to use this subscript later always.
2 dimension piezoelectric patches for we are concerned about have
T 3=0 and E 1=E 2=0 (2-4)
And we also only are concerned about D 3, S 1And S 2, equation (2-1) (2-2) becomes
D 3=d 31T 1+d 32T 2=d p(T 1+T 2)(E 3=0) (2-5)
S 1 S 2 = d 31 d 32 E 3 = d p d p E 3 , ( T 1 = T 2 = 0 ) - - - ( 2 - 6 )
At the free (T of mechanics 1=T 2=0) under the condition, piezoelectric sheet also has the dielectric equation as a common electric capacity
D 3 = &epsiv; p T E 3 , ( T 1 = T 2 = 0 ) - - - ( 2 - 7 )
ε wherein p TBe the free (T of mechanics 1=T 2=0, by subscript T sign) dielectric coefficient under the condition
As elastomeric piezoelectric sheet, electrical short (E is arranged 3=0) elastic constitutive relation under the condition
S 1 S 2 = 1 E p 1 - &mu; p - &mu; p 1 T 1 T 2 , ( E 3 = 0 ) - - - ( 2 - 8 )
E wherein pAnd μ pBe respectively the modulus of elasticity and the Poisson's ratio of piezoelectric patches.
In conjunction with equation (2-5)---(2-8), obtain the machine of piezoelectric patches---the electric coupling constitutive equation
S 1 S 2 D 3 = 1 E p - &mu; p E p d p - &mu; p E p 1 E p d p d p d p &epsiv; p T T 1 T 2 E 3 - - - ( 2 - 9 )
It has reflected 6 amount (4 mechanical quantity T of piezoelectric patches all sidedly 1, T 2, S 1, S 2With 2 electrical quantities D 3And E 3) between coupled relation, both comprised direct piezoelectric effect, comprised inverse piezoelectric effect again.
In different occasions, with different independents variable and dependent variable with more convenient.For this reason, call first kind piezoelectric equations to equation (2-9); Therefrom, can also lead with its other of equal value three types of piezoelectric equations.
Second type of piezoelectric equations
T 1 T 2 D 3 = E p 1 - &mu; p 2 &mu; p E p 1 - &mu; p 2 - e p &mu; p E p 1 - &mu; p 2 E p 1 - &mu; p 2 - e p e p e p &epsiv; p S S 1 S 2 E 3 - - - ( 2 - 10 )
Piezoelectric modulus wherein
e p = E p d p 1 - &mu; p - - - ( 2 - 11 )
And
&epsiv; p S = &epsiv; p T - 2 e p d p = &epsiv; p T - 2 E p d p 2 1 - &mu; p - - - ( 2 - 12 )
For piezoelectric patches at clamping (S 1=S 2=0, with subscript S sign) dielectric constant under the condition.
The 3rd type of piezoelectric equations
S 1 S 2 E 3 = 1 E p - d p 2 &epsiv; p T - &mu; p E p - d p 2 &epsiv; p T g p - &mu; p E p - d p 2 &epsiv; p T 1 E p - d p 2 &epsiv; p T g p - g p - g p 1 &epsiv; p T T 1 T 2 D 3 - - - ( 2 - 13 )
Piezoelectric modulus wherein
g p = d p &epsiv; p T - - - ( 2 - 14 )
The 4th type of piezoelectric equations
T 1 T 2 E 3 = E p 1 - &mu; p 2 + e p 2 &epsiv; p S &mu; p E p 1 - &mu; p 2 + e p 2 &epsiv; p S - h p &mu; p E p 1 - &mu; p 2 + e p 2 &epsiv; p S E p 1 - &mu; p 2 + e p 2 &epsiv; p S - h p - h p - h p 1 &epsiv; p S S 1 S 2 D 3 - - - ( 2 - 15 )
Piezoelectric modulus wherein
h p = e p &epsiv; p S = E p d p &epsiv; p S ( 1 - &mu; p ) - - - ( 2 - 16 )
Piezoelectric patches is usually used in the one-dimensional stress structure like beam and so on, T 2≡ 0, and corresponding four types of piezoelectric equations become
S 1 D 3 = 1 E p d p d p &epsiv; p T T 1 E 3 (first kind T 2≡ 0) (2-17)
T 1 D 3 = E p - e p e p &epsiv; p S S 1 E 3 (second type of T 2≡ 0) (2-18)
S 1 E 3 = 1 E p &epsiv; p S &epsiv; p T g p - g p 1 &epsiv; p T T 1 D 3 (the 3rd type of T 2≡ 0) (2-19)
T 1 E 3 1 E p &epsiv; p T &epsiv; p S - h p - h p 1 &epsiv; p S S 1 D 3 (the 4th type of T 2≡ 0) (2-20)
Wherein
e p=E pd p (2-21)
g p = d p / &epsiv; p T - - - ( 2 - 22 )
h p = e p / &epsiv; p S = E p d p &epsiv; p S - - - ( 2 - 23 )
&epsiv; p S = &epsiv; p T - e p d p = &epsiv; p T - E p d p 2 - - - ( 2 - 24 )
2.2 two kinds of piezoelectric ceramic piece strain transducers
See Fig. 3 c, establish piezoelectric patches two interpolars and constitute the loop through resistance R.By the 4th type of piezoelectric equations (2-15), field intensity does in the piezoelectric patches
E 3 ( t ) = - h p ( S 1 ( t ) + S 2 ( t ) ) + D 3 ( t ) / &epsiv; p S - - - ( 2 - 25 )
Output voltage does
u p ( t ) = &delta; p E 3 ( t ) = - &delta; p h p ( S 1 ( t ) + S 2 ( t ) ) + ( &delta; p / &epsiv; p S ) D 3 ( t ) - - - ( 2 - 26 )
δ wherein pBe piezoelectric patches thickness.The electric current that the unit are piezoelectric patches provides does
i ( t ) = - D &CenterDot; 3 ( t ) = u p ( t ) R = - &delta; p h p R ( S 1 ( t ) + S 2 ( t ) ) + &delta; p &epsiv; p S R D 3 ( t ) - - - ( 2 - 27 )
Therefrom can be about electric displacement D 3(t) differential equation of first order
R &epsiv; p S &delta; p D &CenterDot; 3 ( t ) + D 3 ( t ) = &epsiv; p S h p ( S 1 ( t ) + S 2 ( t ) ) - - - ( 2 - 28 )
Separating in frequency domain is (in this article, a variable will adopt prosign at time domain t in frequency domain ω or the Laplace territory s)
D 3 ( j&omega; ) = h p &epsiv; p S 1 + j&omega; R &epsiv; p S &delta; p ( S 1 ( j&omega; ) + S 2 ( j&omega; ) ) - - - ( 2 - 29 )
In generation, returned (2-27), obtains in frequency domain
i ( j&omega; ) = - j&omega; D 3 ( j&omega; ) = - j&omega; h p &epsiv; p S 1 + j&omega; R &epsiv; p S &delta; p ( S 1 ( j&omega; ) + S 2 ( j&omega; ) ) - - - ( 2 - 30 )
Consider two kinds of special cases.One of which, R=0 has (noting equation 2-16)
i(jω)=-jωe p(S 1(jω)+S 2(jω)) (2-31)
i ( t ) = - e p ( S &CenterDot; 1 ( t ) + S &CenterDot; 2 ( t ) )
The electric current that piezoelectric patches provides is directly proportional with its strain rate sum along 1 and 2; Its two, R is very big, so that in frequency domain, have
&omega; R &epsiv; p S &delta; p > > 1 - - - ( 2 - 32 )
Equation (2-30) is approximately
i ( j&omega; ) &ap; - &delta; p h p R ( S 1 ( j&omega; ) + S 2 ( j&omega; ) ) - - - ( 2 - 33 )
i ( t ) &ap; - &delta; p h p R ( S 1 ( t ) + S 2 ( t ) )
The electric current that piezoelectric patches provides is directly proportional with its normal strain sum along 1 and 2
Equation (2-30) or (2-31) with (2-33) constituted the basis of piezoelectric patches as structural strain rate or strain transducer.
If paste piezoelectric patches on the structure " ideal "---the local strain facies that so-called " the desirable stickup " is meant strain that piezoelectric patches pastes the face place and structure is together; It is very thin to establish piezoelectric patches again, can ignore the variation of its strain along thickness, like this, and the strain S in the piezoelectric patches plane domain 1And S 2Respectively with the local strain stress of structure xAnd ε yConsistent
S 1(x,y,t)=ε x(x,y,t) S 2(x,y,t)=ε y(x,y,t) (2-34)
2.2.1 piezoelectric strain rate transducer
Directly receive piezoelectric patches the negative input end of a linear operational amplifier and (see Fig. 3 a).Because negative input end is " virtual earth ", current potential is approximately zero, so the approximate establishment of equation (2-31), and the amplifier output voltage does
u s ( t ) = - R f &Integral; &Integral; &Omega; S i ( x , y , t ) dxdy = R f e p &Integral; &Integral; &Omega; S ( &epsiv; &CenterDot; x ( x , y , t ) + &epsiv; &CenterDot; y ( x , y , t ) ) dxdy - - - ( 2 - 35 )
R wherein fBe amplifier feedback resistance, Ω sThe structural region that spreads all over for piezoelectric patches; It shows: piezoelectric patches---the amplifier combination is as " strain rate sensor " of structure; Output voltage is proportional to the integration of the local normal strain rate sum of structure; When enough hour of patch area, will be tending towards directly being proportional to structure when place normal strain rate sum.
2.2.2 piezoelectric strain transducer
The negative input end (see Fig. 3 b) of piezoelectric patches through a big resistance R access amplifier, the piezoelectric patches situation will be as Fig. 3 c, and under condition (2-32), the amplifier output voltage does
u s ( t ) = - R f &Integral; &Integral; &Omega; s i ( x , y , t ) dxdy
= &delta; p h p R f R &Integral; &Integral; &Omega; s ( &epsiv; x ( x , y , t ) + &epsiv; y ( x , y , t ) ) dxdy , ( &omega;R &epsiv; p S &delta; p > > 1 ) - - - ( 2 - 36 )
It shows: piezoelectric patches---amplifier combination is as " strain transducer " of structure, and output voltage is proportional to the integration of the local normal strain sum of structure, when enough hour of patch area, will be tending towards directly being proportional to when place normal strain sum.
Because the input impedance of amplifier positive input terminal is very big, therefore directly receive piezoelectric patches the positive input terminal (seeing Fig. 3 b) of amplifier, also will constitute strain transducer.
For one-dimensional stress structure (T 2(x, t)=σ y(x, t)=0), output is respectively with (2-36) strain transducer corresponding to the strain rate of equation (2-35)
u s ( t ) &ap; b p R f e p &Integral; &Omega; s &epsiv; &CenterDot; x ( x , t ) dx - - - ( 2 - 37 )
u s ( t ) &ap; b p &delta; p h p R f R &Integral; &Omega; s &epsiv; x ( x , t ) dx - - - ( 2 - 38 )
B wherein pBe the piezoelectric patches width.Note piezoelectric constant e wherein pAnd h pConfirm by equation (21-24).
2.3 the transfer function model of piezoelectric beam
We will paste some piezoelectric patches on structure to be controlled, utilize its direct piezoelectric effect to make structural strain or strain rate sensor (" piezoelectric transducer "); Simultaneously, paste some piezoelectric patches again in addition and receive voltage drive, utilize its inverse piezoelectric effect as actuator (" piezoelectric actuator ").In this article, call " piezoelectric structure " to this structure that has piezoelectric transducer and actuator simultaneously, like " piezoelectric beam ", " piezoelectric board " etc.
Beam because of its mechanical model is simple, has the analytic solutions of succinct explicit physical meaning as the most basic a kind of continuous structure, usually selected basic research object as the development of novel vibrating control technology; In addition, have a large amount of flexible structures also can be classified as the category of beam really in the actual engineering, therefore, the modeling of piezoelectric beam and vibration control become one of research topic that people pay close attention to the most, development is also comparatively ripe.But be noted that the present invention discusses with beam as an example, the method for this invention is applicable to other various structures, and just the theoretical modeling of structure has difference.
2.3.1 the piezoelectric actuator analysis on the beam
Consult Fig. 4, length and width and the height of establishing beam are respectively L, and b and h correspond respectively to x, y and z axle.Be located at the latitude of emulsion (x of beam 1, x 2) respectively to post a width be b to interior upper and lower surfaces p, thickness is δ pPiezoelectric ceramic piece, polarised direction is all along positive z axle, simultaneously extrinsic motivated voltage u a(t).
If beam is the Bernoulli-Euller beam.Because beam is one-dimensional stress structure (σ y(t) ≡ 0 for x, z), so piezoelectric patches is also located one-dimensional stress state (T 2(x, t) ≡ 0), according to second type of piezoelectric equations (2-18), upper surface piezoelectric patches section stress does
T 1=E pS 1-e pE 3(e p=E pd p) (2-39)
Field intensity in the piezoelectric patches wherein
E 3=u ap (2-40)
If (the consulting condition 2-34) that piezoelectric patches is pasted at Liang Shangshi " ideal ", the upper surface strain of beam does
ε x0=ε x| z=h/2=S 1 (2-41)
It with (2-40) together for returning (2-39), have
T 1 = E p &epsiv; x 0 - E p d p &delta; p u a - - - ( 2 - 42 )
Can know that by symmetry lower surface piezoelectric patches section stress is-T 1Thereby piezoelectric patches to the excitation moment of flexure of beam does
M = - ( T 1 b p &delta; p ) h = - E p b p &delta; p h &epsiv; x 0 + h b p E p d p &delta; p u a - - - ( 2 - 43 )
Above-mentioned excitation moment is converted the cross section moment of flexure of beam
M = 2 EI &epsiv; x 0 h , ( I = 1 12 b h 3 ) - - - ( 2 - 44 )
Wherein E and I are respectively the modulus of elasticity and the cross section moments of inertia of beam.Couplet is solved an equation (2-43) and (2-44) is had
M ( t ) = h b p E p d p 1 + 6 E p &delta; p b p Ehb u a ( t ) - - - ( 2 - 45 )
Equation (2-45) shows
(1) because u a(t) do not become, so the cross section of piezoelectric actuator excitation moment M (t) do not become with x yet, be the beam latitude of emulsion (x that spreads all at piezoelectric patches with coordinate x 1, x 2) interior even moment; Perhaps, also can be regarded as x at the piezoelectric patches two ends 1And x 2The place has applied a pair of reverse moment M (t) (see figure 4).
(2) fillip of piezoelectric actuator size is removed its piezoelectric constant of direct direct ratio d pDepend on that also section rigidity between piezoelectric patches and beam is than (Rigidity Matching) outward.Elastic modulus E pBigger piezoelectric ceramic piece is much littler that piezoelectric membrane (PVDF) has bigger driving force than modulus of elasticity, this just we to select piezoelectric ceramic piece for use be one of fundamental cause of actuator.
2.3.2 the infinitesimal piezoelectric actuator is to the transfer function of transducer
According to the modal theory of beam, the amount of deflection of beam expands into by its natural mode of vibration
w ( x , t ) = &Sigma; r = 1 &infin; &Phi; r ( x ) q r ( t ) - - - ( 2 - 46 )
Φ wherein r(x) and q r(t) be respectively r rank natural mode of vibration and corresponding modal coordinate, its equation of motion in Laplace territory (s) does
q r ( s ) = f r ( s ) m r s 2 + k r - - - ( 2 - 47 )
M wherein r, k r, f r(s) be respectively r rank modal mass, rigidity and mode generalized force.
Impose on x=x aAnd x=x a+ dx aThe mode generalized force that a pair of opposing torque M (t) at place produces
f r ( x a , t ) = M ( t ) ( &Phi; r x ( x a + d x a ) - &Phi; r x ( x a ) ) = M ( t ) &Phi; r xx ( x a ) d x a - - - ( 2 - 48 )
Wherein subscript x representative is to the x differentiate; In generation, returned (2-47) and (2-46), and the response of beam in the Laplace territory does
w ( x , x a , s ) = &Sigma; r = 1 &infin; &Phi; r ( x ) f r ( x a , t ) m r s 2 + k r = M ( s ) d x a &Sigma; r = 1 &infin; &Phi; r ( x ) &Phi; r xx ( x a ) m r s 2 + k r - - - ( 2 - 49 )
Be located at x=x at present aThe a pair of wide b that is is posted at the place p, long is dx aPiezoelectric actuator, according to the analysis of 2.3.1 joint,, just can obtain the response of beam to piezoelectric actuator as long as equation (2-45) substitution (2-49)
w ( x , x a , s ) = h E p d p 1 + 6 E p &delta; p b p Ehb ( b p d x a ) u a ( s ) &Sigma; r = 1 &infin; &Phi; r ( x ) &Phi; r xx ( x a ) m r s 2 + k r - - - ( 2 - 50 )
Accordingly, the surface strain of beam does
&epsiv; x 0 ( x , x a , s ) = h 2 w xx ( x , x a , s ) - - - ( 2 - 51 )
Be located at x=x at present sA wide b of being is posted at the place Ps, thick is δ Ps, long is dx s(δ wherein Ps, h PsDeng the subscript s in variable sign " piezoelectric transducer ", with the respective amount difference of piezoelectric actuator) the infinitesimal piezoelectric patches and after connect amplifier and constitute piezoelectric strain transducer (Fig. 3 b), according to equation (2-38), (2-51) with sensor output voltage (2-50) is arranged
u s ( s ) = &delta; ps h ps R f R &epsiv; x 0 ( x s , x a , s ) b ps d x s
= &alpha; u a ( s ) &Sigma; r = 1 &infin; &Phi; r xx ( x s ) &Phi; r xx ( x a ) m r s 2 + k r ( b ps d x s ) ( b p d x a ) - - - ( 2 - 52 )
Wherein
&alpha; = h &delta; ps h ps R f 2 R h E p d p 1 + 6 E p b p h p Ehb - - - ( 2 - 53 )
Hence one can see that, from unit are (b pDx a=1) the driving voltage u of piezoelectric actuator aTo unit are (b PsDx s=1) the output voltage u of piezoelectric transducer sTransfer function do
H ( x s , x a , s ) = u s ( s ) u a ( s ) = &alpha; &Sigma; r = 1 &infin; &Phi; r xx ( x s ) &Phi; r xx ( x a ) m r s 2 + k r - - - ( 2 - 54 )
In like manner, if usefulness be that (Fig. 3 a) will have piezoelectric strain rate transducer
H ( x s , x a , s ) = u s ( s ) u a ( s ) = &alpha; v &Sigma; r = 1 &infin; &Phi; r xx ( x s ) &Phi; r xx ( x a ) m r s 2 + k r - - - ( 2 - 55 )
Wherein
&alpha; v = - h R f e ps 2 h E p d p 1 + 6 E p b p &delta; p Ehb - - - ( 2 - 56 )
2.3.3 the finite size piezoelectric actuator is to the transfer function of transducer
Establish piezoelectric actuating sheet and vane all is finite size at present, parameter d p, b p, d Ps, b PsDeng not becoming with beam coordinate x.According to equation (2-54), the transfer function from finite size piezoelectric patches actuator to finite size piezoelectric patches strain transducer does
H ( &Omega; s , &Omega; a , s ) = u s ( s ) u a ( s ) = &alpha; b ps b p &Integral; &Omega; s d x s &Integral; &Omega; a H ( x s , x a , s ) d x a
Figure G2009100840041D00108
Wherein
Figure G2009100840041D00109
Ω s, Ω aRepresent the zone of piezoelectric sensing sheet and the beam that spreads all over as moving plate respectively.If usefulness is strain rate sensor, have accordingly
Figure G2009100840041D00111
3 piezoelectric structure damping controls and partial excitation strain compensation
3.1 the local negative velocity feedback damping control of piezoelectric beam experimental exploring
Fig. 5 is the typical control model that numerous documents adopt, and has also provided detailed diagram in the background technology part, will do quantitative explanation here, so that better understand the argumentation in the background technology part.
Modeling according to 2 joint beams has open-loop transfer function equation (equation 2-59)
Figure G2009100840041D00112
Figure G2009100840041D00113
square journey (2-58) wherein.A is the beam zone that piezoelectric actuating sheet and vane spread all over.Being that the upper and lower surfaces of beam respectively pastes one and makes moving plate when noting 2 joint modelings, is only to paste a slice at lower surface here, but from the overall bending response of beam, is equivalent to the alpha of equation (3-1) vReduce by half.
Press the control block diagram on Fig. 5 the right, can be in the hope of doing to closed loop transfer function,
H &OverBar; v ( s ) = u s ( s ) u e ( s ) = g a H v ( s ) 1 + H v ( s ) g a g s - - - ( 3 - 2 )
If each rank mode natural frequency of beam disperses the natural frequency Ω on the r rank rContiguous, the response of all the other each rank mode can be ignored, thereby has
Generation is returned (3-2) to be had
Figure G2009100840041D00116
This shows that r rank mode has increased damping coefficient
Figure G2009100840041D00117
Or damping ratio
Along with loop gain g sg aAnd increase, it seems it is very desirable.
But experimental result is really not so satisfactory.Along with g sg aIncrease, system but trends towards instability, produces high frequency or the very self-excited vibration of low frequency, and before this modal damping also do not reach expectation than increment requirement as 0.1.
It and theoretical model are found in the strain open loop frequency response (Fig. 6 b) that observation is surveyed
H(jω)=H v(jω)/jω (3-7)
(see that Fig. 6 a) has a remarkable difference:
(1) phase-frequency characteristic: behind the pi/2 of resonance region whereabouts, but return to soon and approach zero.
(2) amplitude-frequency characteristic: after crossing resonance region is not to decay down with the 40db/oct slope, but returns to an approximate constant always.
To piezoelectric actuator particularly the modeling of transducer analyze again and computer simulation experiment; Show that this difference derives from the interference of actuator partial excitation strain to the strain of piezo-electric measurement sheet, and this of retraining just that the closed-loop system feedback oscillator can not improve important reasons very.
3.2 self-examination and correction to the piezoelectric structure modeling
When 2 joint piezoelectric structure modelings, we have continued to use traditional excitation in fact---and a kind of acquiescence of sensing pattern: actuator encourages the local deformation that only produces the structural strain's of structure and ignore application of force place; Transducer is only experienced the structural strain's of structure and is ignored the influence that the coordination actuator encourages the local deformation that produces.For traditional power or basic excitation---accelerometer response and so on, this hypothesis are enough reasonably really,---piezoelectric sensing sheet this " strain excitation---strain sensing now at the piezoelectric actuating sheet " new situation under, just possibly become problem.
With the piezoelectric beam is example, piezoelectric actuating sheet strain stress aIn comprise two parts
ε a=-ε mta (3-8)
ε wherein mBe the overall bending strain response of beam, and ε TaIt is the partial excitation strain component that produces as driving source because of it.The strain stress of piezoelectric sensing sheet sAlso comprise corresponding two parts
ε s=ε mts (3-9)
ε wherein TsBe that it is to piezoelectric actuating sheet partial excitation strain ε TaThe response that comes through structure-borne.Like this, from voltage for piezoelectric actuation u aTo strain transducer ε sTransfer function do
H &epsiv; ( s ) = &epsiv; s ( s ) u a ( s ) = &epsiv; m ( s ) u a ( s ) + &epsiv; ts ( s ) u a ( s ) = &epsiv; m ( s ) u a ( s ) + &epsiv; ts ( s ) &epsiv; ta ( s ) &epsiv; ta ( s ) u a ( s ) - - - ( 3 - 10 )
First wherein is equivalent to save the transfer function that obtains in the modelings 2; Factor ε in the second portion Ta(s)/u a(s) depend on the drive characteristic of piezoelectric actuating sheet, can be considered constant; And another factor ε Ts(s)/ε Ta(s) then depend on the propagation characteristic of stress wave in structure, relevant with vane with the relative position of making moving plate, when two coordination configurations, because the two at a distance of recently, will reach maximum.
Like this, to the piezoelectricity modeling of 2 joints, at least at vane from making moving plate when very near, must introduce one and represent the correction term of making moving plate partial excitation strain transmission characteristic in structure.For example equation (3-1) should change into
Figure G2009100840041D00122
To this correction term H t(s) quantitative analysis will be a problem that is worth research, when having a plurality of piezoelectric patches to encourage, can seem more complicated.It was gratifying that it is not difficult to discern through experiment parameter and estimates and go to eliminate or " compensation " through some simple strategies, go that this will be 3.3 to save the problem that will study thereby the piezoelectric structure Mathematical Modeling is revert to again in " routine ".
We have introduced normal correction term H in computer simulation experiment t(s)=Const, when it when 0 increases, can see that frequency response curve will (Fig. 6 a) turns to " piezoelectricity pattern " (Fig. 6 b) gradually from the pattern of routine.In experiment, we add survey from making the moving plate frequency response curve u of the piezoelectric sensing sheet in (" strange land ") nearby St(j ω)/u a(j ω)/(j ω) be illustrated among Fig. 6 c, and it approaches conventional frequency response pattern very much, and (Fig. 6 a); This is because make the cause (Saint Venant's principle) that the propagation of moving plate partial excitation strain decays rapidly with distance.
At 2 joints, we post piezoelectric patches with two sides up and down and constitute piezoelectric actuator together, as if with present analysis have only the different of a piezoelectric patches, actually this is not so.In fact, moving plate excitation strain component, ε are made in the direct impression of piezoelectric sensing sheet at this moment Ts(s)/ε Ta(s)=1, thus make the response of partial excitation strain become bigger.
3.3 the compensation of piezoelectric actuator partial excitation strain
Note
Figure G2009100840041D00131
Open loop frequency response (3-11) becomes
H v ( s ) = u s ( s ) u a ( s ) = b ps b p &alpha; v s ( H m ( s ) + H t ( s ) ) - - - ( 3 - 13 )
The control block diagram of Fig. 5 becomes Fig. 7 (not having Δ H (s) branch road for the time being).The closed loop frequency response does
H &OverBar; v ( s ) = u s ( s ) u e ( s ) &ap; g a b ps b p &alpha; v s ( H m ( s ) + H t ( s ) ) 1 + g s g a b ps b p &alpha; v s ( H m ( s ) + H t ( s ) ) - - - ( 3 - 14 )
Corresponding to the approximate equation in the arrowband (3-3) and (3-4)
Figure G2009100840041D00134
Figure G2009100840041D00135
Become quite complicated, no longer include succinct damping incremental solution (3-5) or (3-6).
In fact, just because of H t(s) objective reality after constituting closed-loop system, has been damaged the controlling performance of flexural vibrations.Manage to eliminate H t(s) influence becomes a major issue that improves the piezoelectric structure control of quality.
See Fig. 7, in theory, do not have too big difficulty: increase a bypass correction link
ΔH(s)=-H t(s) (3-17)
That's all.In other words, the compensation of this link will revert to the transfer function model of piezoelectric structure among the result that 2 joints derive and go.
Comprehensive the above; The present invention realizes through following technical scheme: with piezoelectric beam on moving plate two piezoelectric patches identical of doing with vane be attached on a material little substrate identical with the beam matrix and constituted " piezoelectricity compensating plate " with thickness; It do receive same excitation as moving plate on moving plate and the beam; And connect amplifier behind the vane, its arrange also with beam on transducer the same.
Because adopt technique scheme, damping control method provided by the invention has been eliminated the feedback of piezoelectric actuating sheet partial excitation strain,, the Mathematical Modeling of piezoelectric structure goes thereby being revert in the theoretical model again.
See shown in Figure 2ly, a kind of piezoelectric structure damping control object compensation method is that as an exampleBSEMGVR takeN-PSVSEMOBJ is bright with the beam, and these method concrete steps are following:
Step 1: matrix is set.Get a little substrate 8, the material of this substrate is all identical with beam 1 with thickness.
Step 2: be provided with and make moving plate.With beam 1 on the moving plate 2 identical moving plates 9 of doing of doing be attached on the little substrate 8, make moving plate 9 and receive same excitation with the moving plate 2 of doing on the beam 1.
Step 3: vane is set.With beam 1 on measured sheet 3 identical measured sheet 10 be attached to the another side of little substrate 8, with original do moving plate 9 over against.Connect operational amplifier 12 and resistance 11 after the measured sheet 10; The same with measured sheet 3 resistance 4 at the back on the beam 1 with operational amplifier 5; And measured sheet 10, the annexation of operational amplifier 12 and resistance 11 also with measured sheet 3, resistance 4 is the same with the annexation of operational amplifier 5.
Step 4: regulate gain g cBecause chip area is very little; In the non-high frequency mode frequency band of beam; The bending strain of substrate can be ignored, and has only the tension and compression strain, thus it can do the propagation of the local excitation of moving plate tension and compression strain by the direct modeling beam from doing the moving plate strain to the propagation of vane to vane; In other words, the tension and compression strain transfer function ε of the two Ts(s)/ε Ta(s) identical or proportional.Like this, as long as regulate gain g cJust can reach " compensation " condition (3-17)
&Delta;H ( s , g c ) = &Delta; u s ( s , g c ) u a ( s ) = - H t ( s ) - - - ( 3 - 18 )
g cShould approach-1.
The advantage and the effect of this method are: it has eliminated the feedback of piezoelectric actuating sheet partial excitation strain, goes thereby the Mathematical Modeling of piezoelectric structure is revert in the theoretical model again.So it is ingenious, the simple to operate damping control method of a kind of design.
(4) description of drawings
The local negative velocity feedback damping control of Fig. 1 piezoelectric beam sketch map
Damping control sketch map under the compensation of Fig. 2 partial excitation strain
Fig. 3 piezoelectric patches transducer sketch map (a strain rate sensor b strain transducer c basic circuit)
Piezoelectric actuator is analyzed sketch map on Fig. 4 beam
The local negative velocity feedback damping control of Fig. 5 piezoelectric beam sketch map
Fig. 6 strain open loop frequency response sketch map (a theoretical model b actual measurement---local c actual measurement---strange land)
The influence of Fig. 7 partial excitation strain and compensation sketch map
Wherein symbol description is following among the figure:
1 beam, 2 make moving plate, 3 measured sheet, 4 resistance R _ f, 5 operational amplifiers, 6 gain g s, 7 gain g a,
8 little substrates, 9 make moving plate, 10 measured sheet, 11 resistance, 12 operational amplifiers, 13 gain g c
(5) embodiment
See shown in Figure 2ly, the present invention is a kind of piezoelectric structure damping control object compensation method, is that as an exampleBSEMGVR takeN-PSVSEMOBJ is bright with the beam, and these method concrete steps are following:
Step 1: matrix is set.Get a little substrate 8, the material of this substrate is all identical with beam 1 with thickness.
Step 2: be provided with and make moving plate.With beam 1 on the moving plate 2 identical moving plates 9 of doing of doing be attached on the little substrate 8, make moving plate 9 and receive same excitation with the moving plate 2 of doing on the beam.
Step 3: vane is set.With beam 1 on measured sheet 3 identical measured sheet 10 be attached to the another side of little substrate 8, with original do moving plate 9 over against.Connect operational amplifier 12 and resistance 11 after the measured sheet 10; The same with measured sheet 3 resistance 4 at the back on the beam 1 with operational amplifier 5; And measured sheet 10, the annexation of operational amplifier 12 and resistance 11 also with measured sheet 3, resistance 4 is the same with the annexation of operational amplifier 5.
Step 4: regulate gain g cBecause chip area is very little; In the non-high frequency mode frequency band of beam; The bending strain of substrate can be ignored, and has only the tension and compression strain, thus it can do the propagation of the local excitation of moving plate tension and compression strain by the direct modeling beam from doing the moving plate strain to the propagation of vane to vane; In other words, the tension and compression strain transfer function ε of the two Ts(s)/ε Ta(s) identical or proportional.Like this, as long as regulate gain g cJust can reach " compensation " condition (3-17)
&Delta;H ( s , g c ) = &Delta; u s ( s , g c ) u a ( s ) = - H t ( s ) - - - ( 3 - 18 )
g cShould approach-1.

Claims (1)

1. piezoelectric structure damping control object compensation method, it is characterized in that: these method concrete steps are following:
Step 1: matrix is set; Get a little substrate (8), the material of this substrate is all identical with beam (1) with thickness;
Step 2: be provided with and make moving plate; With beam (1) on first do identical second the making moving plate (9) and be attached on the little substrate (8) of moving plate (2), second makes first on moving plate (9) and the beam makes moving plate (2) and receives same excitation;
Step 3: vane is set; With beam (1) on identical second measured sheet (10) of first measured sheet (3) be attached to the another side of little substrate (8), with original second do moving plate (9) over against; Connect operational amplifier (12) and resistance (11) after second measured sheet (10); The same with first measured sheet (3) resistance (4) at the back on the beam (1) with operational amplifier (5); And second measured sheet (10); The annexation of operational amplifier (12) and resistance (11) also with first measured sheet (3), resistance (4) is the same with the annexation of operational amplifier (5);
Step 4: regulate gain g cBecause chip area is very little; In the non-high frequency mode frequency band of beam; The bending strain of substrate can be ignored; Have only the tension and compression strain, thus it can make the propagation of the local excitation of moving plate tension and compression strain to vane, i.e. the tension and compression strain transfer function ε of the two by the direct modeling beam from doing the moving plate strain to the propagation of vane Ts(s)/ε Ta(s) identical, promptly proportional, like this, as long as regulate gain g cJust can reach " compensation " condition,
&Delta;H ( s , g c ) = &Delta;u s ( s , g c ) u a ( s ) = - H t ( s )
g cShould approach-1.
CN2009100840041A 2009-05-12 2009-05-12 Piezoelectric structure damping control object compensation method Expired - Fee Related CN101562409B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100840041A CN101562409B (en) 2009-05-12 2009-05-12 Piezoelectric structure damping control object compensation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100840041A CN101562409B (en) 2009-05-12 2009-05-12 Piezoelectric structure damping control object compensation method

Publications (2)

Publication Number Publication Date
CN101562409A CN101562409A (en) 2009-10-21
CN101562409B true CN101562409B (en) 2012-06-20

Family

ID=41221061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100840041A Expired - Fee Related CN101562409B (en) 2009-05-12 2009-05-12 Piezoelectric structure damping control object compensation method

Country Status (1)

Country Link
CN (1) CN101562409B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2498076A1 (en) 2011-03-11 2012-09-12 Hexagon Technology Center GmbH Wear-Monitoring of a Gearbox in a Power Station
EP2505956A1 (en) * 2011-03-29 2012-10-03 Hexagon Technology Center GmbH Coordinate measuring machine
CN113078848B (en) * 2021-05-20 2022-05-17 江苏科技大学 Vibrating piezoelectric power generation device and piezoelectric method thereof

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030209953A1 (en) * 2002-05-08 2003-11-13 Pohang University Of Science And Technology Foundation Multi-mode vibration damping device and method using negative capacitance shunt circuits
CN101017989A (en) * 2007-02-28 2007-08-15 西南科技大学 Self-supply micro radio sensing network node based on the piezoelectric vibration power generation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030209953A1 (en) * 2002-05-08 2003-11-13 Pohang University Of Science And Technology Foundation Multi-mode vibration damping device and method using negative capacitance shunt circuits
CN101017989A (en) * 2007-02-28 2007-08-15 西南科技大学 Self-supply micro radio sensing network node based on the piezoelectric vibration power generation

Also Published As

Publication number Publication date
CN101562409A (en) 2009-10-21

Similar Documents

Publication Publication Date Title
Cho et al. Five-port equivalent electric circuit of piezoelectric bimorph beam
Damanpack et al. Active control of geometrically non-linear transient response of sandwich beams with a flexible core using piezoelectric patches
CN101562409B (en) Piezoelectric structure damping control object compensation method
Jin et al. The effect of different shapes of cantilever beam in piezoelectric energy harvesters on their electrical output
Dunsch et al. Unified mechanical approach to piezoelectric bender modeling
Hosseini et al. Nonlinear-forced vibrations of piezoelectrically actuated viscoelastic cantilevers
ur Rahman et al. An experimental study on dynamic analysis and active vibration control of smart laminated plates
Shivashankar et al. Design, modeling and testing of d33-mode surface-bondable multilayer piezoelectric actuator
Narita et al. Bending and polarization switching of piezoelectric laminated actuators under electromechanical loading
Parashar et al. Nonlinear shear-induced flexural vibrations of piezoceramic actuators: experiments and modeling
CN101561685B (en) Piezoelectric structure damping control electronic compensation method
Wereley et al. Demonstration of uniform cantilevered beam bending vibration using a pair of piezoelectric actuators
Cheng et al. An impedance approach for vibration response synthesis using multiple PZT actuators
Hennings et al. Numerical simulation of wave propagation using spectral finite elements
CN201159847Y (en) Time-sharing piezo-electricity self-perception actuator
Jalili et al. Modeling a multilayer piezo-electric transducer by equivalent electro-mechanical admittance matrix
Zehetner et al. Control of torsional vibrations in piezolaminated rods
CN101789711A (en) Construction method for quasi-independent modal actuator
Balmès et al. Validation and verification of FE models of piezo based SHM systems
CN115394384A (en) Piezoelectric mixed functional gradient sandwich hyperbolic plate dynamic response analysis method
Proulx et al. Dynamic analysis of piezoceramic actuation effects on plate vibrations
CN101800485A (en) Construction method of quasi-independent modal sensor
Heidary et al. Dynamic analysis of distributed piezothermoelastic composite plate using first-order shear deformation theory
CN201656808U (en) Quasi-independent modal actuator
Satir et al. Transmit optimization of CMUTs in non-collapse mode using a transient array model

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120620

Termination date: 20140512