CN107515247B - A kind of composite panel part splitting method for detecting position - Google Patents
A kind of composite panel part splitting method for detecting position Download PDFInfo
- Publication number
- CN107515247B CN107515247B CN201710607919.0A CN201710607919A CN107515247B CN 107515247 B CN107515247 B CN 107515247B CN 201710607919 A CN201710607919 A CN 201710607919A CN 107515247 B CN107515247 B CN 107515247B
- Authority
- CN
- China
- Prior art keywords
- shearing strain
- composite panel
- partial derivative
- splitting
- precision
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/045—Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0231—Composite or layered materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/26—Scanned objects
- G01N2291/263—Surfaces
- G01N2291/2632—Surfaces flat
Abstract
The invention discloses a kind of composite panel part splitting method for detecting position, step includes: the vibration deformation for measuring composite panel;The shearing strain formula that vibration deformation substitutes into plate is obtained into shearing strain;Second order mixing partial derivative operator is applied to shearing strain and obtains shearing strain second-order mixed partial derivative;Shearing strain second-order mixed partial derivative and more precision Spline convolution of functions obtain more precision shearing strain second-order mixed partial derivatives;According to whether unusual peak value occur judges composite panel part splitting in more precision shearing strain second-order mixed partial derivatives, and splitting position is determined by peak position.The present invention can be used for the vibration deformation of COMPOSITE MATERIALS plate, and then judges whether there is local splitting and position splitting;This method, as damage criterion, can not only eliminate the influence of ambient noise using more precision shearing strain second-order mixed partial derivatives of explicit physical meaning, moreover it is possible to significantly reinforcement damage characteristic.
Description
Technical field
The invention discloses a kind of damage position detection methods of composite panel, and in particular to plate structural damage detection skill
Art field.
Background technique
Composite panel is also easy to produce small local splitting under percussion, these small local splittings can be held
Continuous extension causes structure that Delamination Buckling even overall collapse destruction occurs.Therefore, in time, effectively diagnosis composite panel is small
Local splitting is impacted to be of great significance.Currently, the local lossless detection method based on ultrasonic wave, electromagnetic field, vortex etc.
It has been widely used in the detection of composite panel part splitting.But these methods generally require the substantially position of precognition damage
It sets, and splitting position must be easily accessible, this more difficult satisfaction in practical applications.And the vibration based on composite panel becomes
Shape forms more precision shearing strain second-order mixed partial derivatives with clear physical significance, and for detecting composite panel part
The research of splitting has not been reported.
Local splitting detection method based on ultrasonic wave, electromagnetic field, vortex etc. for test environment in noise very
Sensitivity is tested in environment at the scene, and noise contribution is likely to cover true local splitting information, influences detection effect.
Summary of the invention
The technical problems to be solved by the present invention are: in view of the drawbacks of the prior art, proposing with clear physical significance
The more precision shearing strain mixing local derviation detecting factors of composite panel part splitting, reach elimination Environmental Noise Influence, together
The purpose of Shi Jiaqiang damage information feature.
The present invention uses following technical scheme to solve above-mentioned technical problem:
The present invention proposes a kind of composite panel part splitting method for detecting position, and specific steps include:
Step 1: measuring the vibration deformation of composite panel, the shearing strain formula that vibration deformation substitutes into plate is obtained cutting and is answered
Become;
Step 2: second order mixing partial derivative operator, which is applied to shearing strain, obtains shearing strain second-order mixed partial derivative;
Step 3: shearing strain second-order mixed partial derivative and more precision Spline convolution of functions are obtained more precision shearing strains two
Rank mixed partial derivative;
Step 4: according to whether unusual peak value occur judges composite material in more precision shearing strain second-order mixed partial derivatives
Plate part splitting, and splitting position is determined by peak position.
As a further preferred embodiment of the present invention, the specific calculation of the step 1 is as follows:
(1) the vibration deformation w (x, y) of composite panel is obtained;
(2) w (x, y) shearing strain formula for substituting into plate is obtained into shearing strain γxy(x, y):
Wherein, x, y, z be respectively the length of plate, width, thickness direction coordinate.
As a further preferred embodiment of the present invention, the specific calculation of the step 2 is as follows:
Second order mixing partial derivative operator is applied to γxyOn (x, y), shearing strain second-order mixed partial derivative is obtained
As a further preferred embodiment of the present invention, the specific calculation of the step 3 is as follows:
It willWith more precision Spline function Spu,v,s(x, y) convolution:
Spu,v,s(x, y)=s-1Sp(x-u)/s,(y-v)/s)
Wherein,Convolution algorithm is represented, s is scale parameter, and u, v are respectively the translation parameters in the direction x, y,For
More precision shearing strain second-order mixed partial derivatives.
As a further preferred embodiment of the present invention, the acquisition modes of the vibration deformation of composite panel are specific in step 1
Are as follows:
Modality vibration exciter is used to apply sinusoidal excitation in the back side geometric center position of composite panel as power driver, together
When with laser sensor measure the positive vibration of composite panel, obtain the vibration deformation of composite panel.
As a further preferred embodiment of the present invention, the laser sensor uses PSV-400 laser scanning vialog, institute
It states power driver and uses 4890 modality vibration exciters.
The invention adopts the above technical scheme compared with prior art, has following technical effect that
Detection method proposed by the invention can be used for the vibration deformation of analysis plates structure, and then judge whether there is part
Splitting simultaneously positions splitting.The invention proposes a kind of more precision shearing strain mixing of composite panel part splitting
Local derviation detecting factor, compared to conventional local lossless detection method, this method has specific physical significance, eliminates environment and makes an uproar
Sound shadow sound, the effect for reinforcing damage characteristic are obvious.Meanwhile the present invention can advanced sensors be mating makes with laser scanning vialog etc.
With efficient application is in noise circumstance lower plate structure damage check.
Detailed description of the invention
Fig. 1 is method flow schematic diagram of the invention.
Fig. 2 is the 6th first order mode schematic diagram of plate in the present invention.
Fig. 3 is shearing strain schematic diagram of the plate under the 6th first order mode in the present invention.
Fig. 4 is shearing strain second-order mixed partial derivative schematic diagram of the plate under the 6th first order mode in the present invention.
Fig. 5 is more precision shearing strain second-order mixed partial derivative schematic diagrames of the plate under the 6th first order mode in the present invention.
Specific embodiment
Technical solution of the present invention is described in further detail with reference to the accompanying drawing:
Those skilled in the art can understand that unless otherwise defined, all terms used herein (including skill
Art term and scientific term) there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.
As shown in Figure 1, composite panel part splitting method for detecting position proposed by the present invention, including following step
It is rapid:
1. obtaining the vibration deformation w (x, y) of composite panel;
2. w (x, y) to be substituted into the shearing strain formula of plate, shearing strain γ is obtainedxy(x, y):
Wherein, x, y, z be respectively the length of plate, width, thickness direction coordinate.
3. second order mixing partial derivative operator is applied to γxyOn (x, y), shearing strain second-order mixed partial derivative is obtained
4.With more precision Spline function Spu,v,s(x, y) convolution:
Spu,v,s(x, y)=s-1Sp((x-u)/s,(y-v)/s) (3)
Wherein,Convolution algorithm is represented, s is scale parameter, and u, v are respectively the translation parameters in the direction x, y,For
More precision shearing strain second-order mixed partial derivatives.
5. basisIn whether there is unusual peak value to judge composite panel part splitting, and by peak value
Position determines splitting position.
The working principle of the invention is specifically described further below as follows:
Composite panel part splitting can make shearing strain mutate, and accordingly, can identify partial layer by shearing strain
Between remove.However, slight, local splitting can only cause the shearing strain for being not easy to observe to be mutated.In order to sufficiently disclose shearing strain
Mutation, is applied to γ for second order mixing partial derivative operatorxyOn (x, y), shearing strain second-order mixed partial derivative is obtainedExisting
During field measures, hardened structure vibration deformation w (x, y) inevitably contains ambient noise.Shearing strain second order mixing local derviation
Although number is sensitive to damage, also inevitably influence of noise can be amplified.At the same time, current domestic and international damage check
Used in laser scanning vialog, as advanced vibration signal sensor, play high-precision, high resolution advantage it is same
When, it can also aggravate this noise jamming effect.
Based on this, the present invention is by shearing strain second-order mixed partial derivativeWith more precision Spline functions
Convolution obtains more precision shearing strain second-order mixed partial derivativesPass through the ruler of more precision shearing strain second-order mixed partial derivatives
Characteristic is spent, noise signal is filtered out in large scale, retains damage signal.According toIn whether occur unusual peak value come
Judge composite panel part splitting, and splitting position is determined by peak position.
Lift a specific example below to illustrate detection method of the invention:
The long 1000mm of GFRP composite panel used in embodiment, wide 1000mm, thick 1.5mm, apart from two adjacent side 300mm of plate
At 550mm, there are a prefabricated long 50mm, the local splitting of wide 50mm.
Firstly, modality vibration exciter is used to apply sinusoidal excitation in back geometric center position as power driver, use simultaneously
The positive vibration of laser scanning vialog measurement plate, obtains the vibration deformation of plate.6th first order mode of plate is as shown in Fig. 2, wherein
Horizontal axis x indicates that the coordinate long along plate, longitudinal axis y indicate that the coordinate wide along plate, vertical pivot w indicate vibration deformation.
Secondly, acquiring the shearing strain γ of plate with formula (1)xy(x, y), as shown in figure 3, wherein vertical pivot γxyExpression, which is cut, answers
Become.
Then, second order mixing partial derivative operator is applied to by γ according to formula (2)xyOn (x, y), the mixing of shearing strain second order is obtained
Partial derivativeAs shown in figure 4, wherein vertical pivotIndicate shearing strain second-order mixed partial derivative.Figure 4, it is seen that by
In the influence of noise, the noise signal that true damage characteristic has been amplified is covered, it is difficult to identification damage.
Then, according to formula (4)With more precision Spline function Spu,v,s(x, y) convolution, scale s take 0.35,
Obtain more precision shearing strain second-order mixed partial derivativesAs shown in Figure 5.It can be found that there is a unusual peak from Fig. 5,
Show there is damage.Peak value corresponding position is the splitting position detected, is consistent with practical splitting position.
It thereby confirms, compared to conventional local lossless detection method, this method not only has specific physical significance,
Environmental Noise Influence can also be eliminated, effectively reinforce damage characteristic, so that quickly and accurately judgement damage occurs and detects layer
Between remove position.
In the present embodiment: laser sensor is the PSV-400 laser scanning vialog of Polytec company, Germany production;Power
Driver is 4890 modality vibration exciters of Denmark B&K company production.
In conclusion the invention proposes a kind of composite panel part splitting method for detecting position, based on mostly smart
Shearing strain mixing local derviation detecting factor is spent, compared to conventional local lossless detection method, this method is anticipated with specific physics
Justice can not only eliminate the influence of ambient noise, moreover it is possible to reinforce the information characteristics of local splitting.Quickly and accurately judge
Local interlayer separates existing and determines splitting position.
The above examples only illustrate the technical idea of the present invention, and this does not limit the scope of protection of the present invention, all
According to the technical idea provided by the invention, any changes made on the basis of the technical scheme each falls within the scope of the present invention
Within.
Claims (4)
1. a kind of composite panel part splitting method for detecting position, which is characterized in that specific steps include:
Step 1: the vibration deformation of measurement composite panel, obtains shearing strain for the shearing strain formula that vibration deformation substitutes into plate;
Step 2: second order mixing partial derivative operator, which is applied to shearing strain, obtains shearing strain second-order mixed partial derivative;
It is mixed Step 3: shearing strain second-order mixed partial derivative is obtained more precision shearing strain second orders with more precision Spline convolution of functions
Close partial derivative;
Step 4: in more precision shearing strain second-order mixed partial derivatives, according to whether unusual peak value occur judges composite panel
Local splitting, and splitting position is determined by peak position;
The specific calculation of the step 1 is as follows:
(1) the vibration deformation w (x, y) of composite panel is obtained;
(2) w (x, y) shearing strain formula for substituting into plate is obtained into shearing strain γxy(x, y):
Wherein, x, y, z be respectively the length of plate, width, thickness direction coordinate;
The specific calculation of the step 2 is as follows:
Second order mixing partial derivative operator is applied to γxyOn (x, y), shearing strain second-order mixed partial derivative is obtained
The specific calculation of the step 3 is as follows:
With more precision Spline function Spu,v,s(x, y) convolution:
Spu,v,s(x, y)=s-1Sp((x-u)/s,(y-v)/s)
Wherein,Convolution algorithm is represented, s is scale parameter, and u, v are respectively the translation parameters in the direction x, y,It is mostly smart
Spend shearing strain second-order mixed partial derivative.
2. a kind of composite panel part splitting method for detecting position as described in claim 1, which is characterized in that step
The acquisition modes of the vibration deformation of composite panel in one specifically:
Apply sinusoidal excitation in the back side geometric center position of composite panel with vibration excitor, while being measured again with laser sensor
The positive vibration of condensation material plate, obtains the vibration deformation of composite panel.
3. a kind of composite panel part splitting method for detecting position as claimed in claim 2, which is characterized in that described
Vibration excitor is 4890 modality vibration exciters.
4. a kind of composite panel part splitting method for detecting position as claimed in claim 2, which is characterized in that described
Laser sensor is PSV-400 laser scanning vialog.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710607919.0A CN107515247B (en) | 2017-07-24 | 2017-07-24 | A kind of composite panel part splitting method for detecting position |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710607919.0A CN107515247B (en) | 2017-07-24 | 2017-07-24 | A kind of composite panel part splitting method for detecting position |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107515247A CN107515247A (en) | 2017-12-26 |
CN107515247B true CN107515247B (en) | 2019-11-22 |
Family
ID=60722475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710607919.0A Active CN107515247B (en) | 2017-07-24 | 2017-07-24 | A kind of composite panel part splitting method for detecting position |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107515247B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988119A (en) * | 2019-10-31 | 2020-04-10 | 河海大学 | Method for detecting layered damage of composite laminated plate by measuring equivalent pseudo load through laser |
CN115876883B (en) * | 2022-12-29 | 2024-03-29 | 南京航空航天大学 | Method and system for detecting layered damage position of composite material laminated plate |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101918787A (en) * | 2007-12-17 | 2010-12-15 | 明产株式会社 | Seat sag evaluation method and device |
CN103575807A (en) * | 2013-10-24 | 2014-02-12 | 河海大学 | Method for detecting structural damage of beam based on Teager energy operator-wavelet transformation curvature mode |
CN104887266A (en) * | 2015-05-29 | 2015-09-09 | 西安交通大学 | Method for small-area three-dimensional passive cavitation imaging and three-dimensional composite imaging based on area array |
CN106950280A (en) * | 2017-03-16 | 2017-07-14 | 东北大学 | Fibre reinforced composites parameter identification method and device based on the lossless scanning of laser |
-
2017
- 2017-07-24 CN CN201710607919.0A patent/CN107515247B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101918787A (en) * | 2007-12-17 | 2010-12-15 | 明产株式会社 | Seat sag evaluation method and device |
CN103575807A (en) * | 2013-10-24 | 2014-02-12 | 河海大学 | Method for detecting structural damage of beam based on Teager energy operator-wavelet transformation curvature mode |
CN104887266A (en) * | 2015-05-29 | 2015-09-09 | 西安交通大学 | Method for small-area three-dimensional passive cavitation imaging and three-dimensional composite imaging based on area array |
CN106950280A (en) * | 2017-03-16 | 2017-07-14 | 东北大学 | Fibre reinforced composites parameter identification method and device based on the lossless scanning of laser |
Non-Patent Citations (3)
Title |
---|
Structural damage identification using damping: a compendium of uses and features;Maosen Cao 等;《Smart Mater. Struct》;20170307;第26卷;第1-14页 * |
Two-Dimensional Modal Curvature for Damage Detection in Plates;Maosen Cao 等;《Transactions of Nanjing University of Aeronautics and Astronautic》;20150630;第32卷(第3期);第255-260页 * |
多向层合板中的脱层与超声检测成像;杨学永 等;《南京航空航天大学学报》;20091231;第41卷;第67-71页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107515247A (en) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102043016B (en) | Lamb wave-based autonomous damage identification imaging method | |
La et al. | Autonomous robotic system for bridge deck data collection and analysis | |
Kharkovsky et al. | Dual-polarized near-field microwave reflectometer for noninvasive inspection of carbon fiber reinforced polymer-strengthened structures | |
CN107515247B (en) | A kind of composite panel part splitting method for detecting position | |
CN100573181C (en) | Utilize the method for ferromagnetic materials surface stray magnetic field signal monitoring fatigue damage | |
CN103063124B (en) | A kind of detection method of austenitic stainless steel amount of plastic deformation | |
CN107917957B (en) | Damage detection method for plate-shaped structure | |
WO2021098197A1 (en) | Apparatus for testing rebar in concrete, and method therefor | |
CN104597083B (en) | Ultrasonic guided wave and electromechanical impedance-based mobile damage detection method | |
CN111597653B (en) | Bridge structure defect dynamic detection and identification method for bridge crane | |
CN112630270A (en) | Method for detecting state of reinforcing steel bar in concrete structure | |
CN104729435A (en) | Sheet thickness online measurement system based on laser ultrasound | |
CN103884777A (en) | Plate structure damage detection method based on wavelet transformation-fractal analysis | |
CN104535656A (en) | All-optical non-contact type composite material plate layer crack damage detection system and method | |
CN104849348A (en) | Beam structure damage detection method based on singular value decomposition-wavelet transform | |
CN112683207A (en) | Method and device for detecting thickness of concrete protective layer and distance between reinforcing steel bars | |
CN110988119A (en) | Method for detecting layered damage of composite laminated plate by measuring equivalent pseudo load through laser | |
CN104181237A (en) | Structural member flaw detection monitoring temperature compensating method and system thereof | |
US20150008912A1 (en) | Method and device for detecting mechanical changes in a component by means of a magnetoelastic sensor | |
CN106198368A (en) | Inside concrete steel bar corrosion method for detecting position | |
US7690260B2 (en) | Method and system having ultrasonic sensor movable by translation device for ultrasonic profiling of weld samples | |
CN103743818B (en) | Based on the damage diagnosis method of the energy flow chart of ripple and the damage system of enforcement the method | |
CN113984881A (en) | Reinforced concrete corrosion detection method, detection device and use method | |
CN102426271A (en) | Kelvin probe measurement method for bonding strength of metal alloy connecting piece interface | |
JP2018132426A (en) | Reinforcement diameter of ferroconcrete, measuring device of covering, and reinforcement arrangement direction measuring method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |