WO2018165999A1 - Fiber reinforced composite material parameter identifying method based on laser nondestructive scanning, and device - Google Patents

Fiber reinforced composite material parameter identifying method based on laser nondestructive scanning, and device Download PDF

Info

Publication number
WO2018165999A1
WO2018165999A1 PCT/CN2017/078513 CN2017078513W WO2018165999A1 WO 2018165999 A1 WO2018165999 A1 WO 2018165999A1 CN 2017078513 W CN2017078513 W CN 2017078513W WO 2018165999 A1 WO2018165999 A1 WO 2018165999A1
Authority
WO
WIPO (PCT)
Prior art keywords
vibration
laser
platform
excitation
composite beam
Prior art date
Application number
PCT/CN2017/078513
Other languages
French (fr)
Chinese (zh)
Inventor
李晖
刘�东
高宁
王雨波
薛鹏程
闻邦椿
Original Assignee
东北大学
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 东北大学 filed Critical 东北大学
Priority to JP2019548923A priority Critical patent/JP7109798B2/en
Publication of WO2018165999A1 publication Critical patent/WO2018165999A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • G01N29/045Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/36Detecting the response signal, e.g. electronic circuits specially adapted therefor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0231Composite or layered materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters

Definitions

  • the invention relates to the field of machinery, in particular to a fiber reinforced composite material parameter identification method and device.
  • Structural composites are widely used in aviation, aerospace, marine, sports equipment due to their high specific strength, high specific modulus, designability, thermal stability, and high bearing capacity and light weight. , electrical equipment, medicine, weapons industry and chemical industry.
  • fiber reinforced composite thin-walled members such as composite blades, composite whole leaf discs and composite cylindrical shells, often work in harsh environments such as high-speed rotation, high temperature, corrosive gas erosion, etc.
  • the resulting vibration fatigue and vibration failure problems are becoming more and more prominent, making it impossible to perform the functions that people have previously envisaged.
  • Patent CN201610166104.9 invented a new type of composite material parameter identification instrument. By changing the structure of the identification instrument, the identified parameters are more accurate. However, the principle is still to measure the composite material by traditional methods such as pressing, shearing and pulling. The parameters do not fundamentally solve the waste, there are certain security risks and other issues.
  • the detection method of fiber reinforced composites needs to be improved, and the application of vibration and laser non-destructive scanning technology can solve the related problems well.
  • the present invention provides a fiber reinforced composite structural material parameter tester based on laser non-destructive scanning and a working method thereof.
  • the specific technical solutions are as follows:
  • a method for identifying a fiber reinforced composite material based on laser non-destructive scanning includes the following steps:
  • Step 1 Open the vacuum casing, adjust the tightening screw, and install the test piece of the composite beam to be tested in the reference position of the clamping mechanism in sequence;
  • Step 2 Gradually adjust the pressure of the clamping mechanism by adjusting the tightening screw, according to the pressure output by the pressure sensor After the judgment, it is judged that the appropriate clamping force has been reached, and the adjustment is stopped;
  • Step 3 Adjust the lifting platform to a suitable position, tighten the vacuum casing; exhaust the gas inside the casing with an air extractor to create a vacuum environment;
  • Step 4 The position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device to be in the free end position of the composite beam test piece to be tested; then, the excitation device is turned on, and the sine sweep is performed in a large frequency range.
  • the frequency excitation test monitors the time domain waveform data of the frequency sweep excitation signal according to the acceleration sensor on the movable clamp body, and obtains the spectrum of the vibration response signal of the free end position of the composite beam test piece to be tested by the laser vibration meter, and passes the half power
  • the bandwidth method identifies the frequency corresponding to each peak and the frequency of the half power point, thereby obtaining the first three natural frequencies and damping ratios of the composite beam specimen to be tested;
  • Step 5 Adjust the frequency of the excitation device to the first-order natural frequency, and the excitation composite beam test piece reaches the first-order resonance state, and the excitation amplitude corresponding to the first-order resonance state is determined by the acceleration sensor on the vibration platform; meanwhile, the laser is turned on.
  • the control switch of the scanning vibration measuring device moves the position of the laser spot emitted by the laser vibrometer by the laser scanning vibration measuring device, and realizes the scanning test of the test piece of the composite beam from the cantilever end position to the free end position, and obtains the composite beam test to be tested.
  • the amplitude of the vibration response signal at each scanning point position, and then the curve of the vibration amplitude of the composite beam specimen with the length of the first-order resonance state is plotted as a first-order resonance amplitude curve.
  • the horizontal axis of the first-order resonance amplitude curve is the length, and the vertical axis is the vibration amplitude;
  • Step 6 adjusting the frequency of the excitation device to the second-order natural frequency and the third-order natural frequency, and illuminating the steps used in step 5 to obtain a second-order resonance amplitude curve and a third-order resonance amplitude curve;
  • Step 7 According to the parameter identification method of the fiber reinforced composite structural material, input the first three natural frequencies, the damping ratio and the resonance amplitude curve, and set the allowable value of the recognition error to 10% to obtain the longitudinal elastic modulus E 1 of the fiber and the transverse elasticity of the fiber.
  • the modulus E 2 , the shear modulus G 12 , the Poisson's ratio is ⁇ 21 and the fiber longitudinal loss factor ⁇ 11 , the fiber transverse loss factor ⁇ 22 and the shear loss factor ⁇ 12 .
  • the invention further discloses a device for realizing the method, which is mainly composed of a double working platform, a lifting platform, a clamping mechanism, an excitation device, a laser scanning vibration measuring system and a vacuum device;
  • the double working platform comprises a first working platform And the second working platform, the two platforms are supported by the "work" word supporting steel, and the bottom is supported by the base;
  • the clamping mechanism is used for fixing and monitoring the composite beam test piece to be tested, including the reference platform and the upper pressing block thereof, The pressure block is driven by bolts to fix the composite beam test piece to be tested;
  • two circular hole structures are arranged on the reference platform for placing the circular pressure sensing gasket of the pressure sensor, and quantitatively adjusting according to the value of the pressure sensor Acting on the pressure of the tested composite material, the quantitative characterization test of the constraint condition is realized;
  • the excitation device consists of two parts: the vibration exciter and the vibration platform, the function of which is to generate the excitation force and transmit the vibration to the test piece of the
  • the laser scanning vibration measuring system comprises a laser scanning vibrometer, a lead screw, a coupling and an electric motor; the platform carrying the laser scanning vibrometer is penetrated by a set of lead screws, and the lead screw is connected with the coupling to form a transmission mechanism; Under the power of the motor, the laser scanning vibrometer can realize the horizontal movement of the laser vibration measuring device through the transmission mechanism; the position of the laser spot emitted by the laser vibrometer can realize the test specimen of the composite beam from the cantilever at a certain scanning speed.
  • the scanning test from the end position to the free end position obtains the amplitude of the vibration response signal of the composite beam test piece at each scanning measuring point position, and then draws a corresponding resonance amplitude curve; the laser scanning vibration measuring system is arranged on the lifting platform
  • the adjustment of the position of the laser scanning vibrometer in the vertical direction can be realized, and the lifting frame of the lifting platform is fixed on the second working platform.
  • a vacuum device is arranged outside the tester, which is composed of a casing and a baffle; the baffle is tightly connected with the casing to prevent air from entering; and the baffle has a circular hole structure connected to the air pump to extract the gas inside the device, thereby creating Vacuum environment reduces experimental error.
  • An advantage of the present invention is that the present invention provides a fiber reinforced composite structural material parameter tester based on laser non-destructive scanning.
  • the parameter test method combining vibration transmission and laser displacement sensing technology was designed, which broke the traditional measurement mode by stretching and other physical methods, and achieved the goal of measuring the mechanical properties of the material to be tested in all directions at one time.
  • the earth reduces the workload of related work;
  • the instrument does not have any impact on the material during the whole process of testing, saves raw materials, reduces economic losses, and is more environmentally friendly; in addition, the external vacuum housing design makes The test process is carried out under vacuum, avoiding the influence of air damping, and the measurement result is more accurate.
  • the work is easy to assemble and disassemble, and the quality of each component is not large, and can be loaded into the pull box for better portability. .
  • the parameters of the composite materials can be measured at different locations according to actual conditions, which provides great convenience for production measurement.
  • FIG. 1 is a front elevational view of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention
  • FIG. 2 is a top plan view of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention
  • FIG. 3 is a schematic structural diagram of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention.
  • Figure 4 is a theoretical model of the fiber reinforced composite beam specimen under the basic excitation of Figure 1;
  • Figure 5 is a theoretical and experimental obtained first-order resonance amplitude curve of the composite beam and its error upper and lower limits
  • the device of the present invention mainly comprises a double working platform, a lifting platform, a clamping mechanism, an excitation device, a laser scanning vibration measuring system and a vacuum device;
  • the compound working platform comprises a first working platform 4 and a second working platform 2, wherein the two platforms are connected by the "work" word supporting steel 3, and the lower part is supported by the base 1;
  • the clamping mechanism 8 is used for fixing and monitoring the composite beam test piece to be tested 15, comprising a reference platform 17 and an upper pressing block 18, the pressing block is driven by a bolt to press the composite beam test piece to be tested;
  • the upper surface of the reference platform 17 is provided with two circular hole structures for placing a circular pressure sensing gasket of the pressure sensor, and quantitatively adjusting the pressure acting on the composite material according to the value of the pressure sensor to realize quantitative characterization of the constraint condition. test.
  • the excitation device is composed of two parts, the vibration exciter 16 and the vibration platform 5, and the function is to generate an excitation force and transmit vibration to the test piece of the composite beam to be tested; the vibration exciter is connected with the signal source to generate vibration of a certain frequency, and the vibration passes through the vibration.
  • the platform 5 and the clamping mechanism 8 are transmitted to the beam test piece 15 to be tested; in order to realize the monitoring of the excitation vibration state, an acceleration sensor 9 is mounted on the upper portion of the clamping mechanism 8 for measuring the magnitude of the vibration amplitude.
  • the vibration exciter 16 is fixed to the first working platform 4 by bolt connection for generating an excitation force and applying it to the vibration platform 5; under the excitation platform 5, there are four shock-absorbing screws uniformly distributed at the four corners, The shock spring 6 can be inserted into the shock absorbing screw to be connected with the second working platform 2, and only the displacement in the direction of the spring axis can be generated between the two.
  • the laser scanning vibration measuring system includes a laser scanning vibrometer 11, a lead screw 10, a coupling, and an electric motor.
  • the platform carrying the laser scanning vibrometer is penetrated by a set of lead screws, and the lead screw is connected with the coupling to form a transmission mechanism; under the power of the motor, the laser scanning vibrometer can realize the horizontal direction of the laser vibration measuring device through the transmission mechanism. exercise.
  • the laser spot position of the laser vibrometer is used to scan the test specimen from the cantilever end position to the free end position at a certain scanning speed, and obtain the position of the composite beam test piece to be tested at each scanning point.
  • the vibration responds to the amplitude of the signal, which in turn plots the corresponding resonance amplitude curve.
  • the laser scanning vibration measuring system is disposed on the lifting platform 14 to adjust the position of the laser scanning vibrometer in the vertical direction.
  • the lifting platform of the lifting platform 14 is fixed on the second working platform 2, and the two cross-hinged Supporting beams, one fixed on the working platform 2, the other one is connected with the sliding block 13, and the sliding block is connected with the motor through the screw; under the driving of the motor, the screw drives the sliding block 13 to slide on the working platform to drive the lifting The platform moves up and down;
  • a vacuum device is provided on the outside of the tester, which is composed of a casing and a baffle.
  • the baffle can be tightly connected to the casing to prevent air from entering; the baffle has a circular hole structure that can be connected to the air pump to extract gas from the device, create a vacuum environment, and reduce experimental errors.
  • the fiber reinforced composite structural material parameter identification method is:
  • the composite beam specimen is composed of n layers of fibers with orthogonal anisotropy characteristics and matrix materials; assuming that the layers are firmly bonded, there is no slip between the layers, no relative displacement, solid can be ignored
  • the effect of the interlayer coupling effect first, the midplane is used as the reference plane, and the xoy coordinate system is established; the angle between the fiber direction and the x-axis direction of the global coordinate system is ⁇ , the plate length is a, the plate width is b, and the plate thickness is For h, each layer is located between the lower surface h k-1 of the z coordinate axis and the higher surface h k , and the thickness of each layer is the same; 1 in the figure represents the longitudinal direction of the fiber, 2 represents the transverse direction of the fiber, and 3 represents the vertical direction. 1-2 plane direction;
  • i a virtual unit
  • the elastic modulus of the composite material is expressed as follows
  • the Poisson's ratio of the directional strain is ⁇ 21 ;
  • z represents the displacement in the z-axis direction
  • u, v, w represents the displacement of any point in the plate
  • u 0 , v 0 , w 0 represents the surface displacement of the plate
  • h is the thickness of the composite beam test piece
  • t represents time
  • ⁇ x represents the positive strain of the point in the x direction
  • ⁇ y represents the positive strain of the point in the y direction
  • ⁇ xy indicates the shear strain of the point on the xy plane
  • the bending curvature and the twist rate of the middle surface of the thin plate can be expressed as
  • k x represents the flexural curvature of the point in the x direction
  • k y represents the flexural curvature of the point in the y direction
  • k xy represents the distortion rate of the point on the xy plane
  • ⁇ x z ⁇ x
  • ⁇ y z ⁇ y
  • ⁇ xy z ⁇ xy
  • 1 indicates the longitudinal direction of the fiber
  • 2 indicates the transverse direction of the fiber
  • 6 indicates the direction of the vertical plate surface
  • Q* indicates the elastic modulus
  • the stress-strain relationship of the k-th layer in the global coordinate system is calculated by the stress-strain rotation formula.
  • k represents the k-th layer of the composite beam test piece
  • ⁇ k represents the angle between the fiber direction of the k-th layer plate and the x-axis of the global coordinate system
  • M x represents the bending moment in the x-axis direction
  • M y represents the bending moment in the y-axis direction
  • M xy represents the torque on the xy plane
  • D* represents the bending stiffness coefficient
  • the basic excitation of the composite beam specimen is equivalent to the uniform inertial force external load.
  • the kinetic energy of the vibration of the thin plate can be expressed by the following formula
  • is the density of the thin plate and h is the thickness of the thin plate
  • the strain energy stored in the bending of the thin plate is expressed by the following formula
  • the uniform inertial force of the thin plate is the uniform inertial force of the thin plate.
  • W ij ( ⁇ , ⁇ ) is a vibration mode function, and has the following form
  • M ⁇ N non-homogeneous linear algebraic equations can be obtained. For the convenience of solving, write it as a matrix form.
  • the natural frequency and mode shape of the composite beam specimen can be obtained; further, it is assumed that the fiber reinforced composite beam specimen is subjected to the basic excitation load and the vibration response ⁇ (x, y) is solved under the fundamental excitation. , t) expression; taking into account the experimental test to obtain the absolute vibration response of the composite beam specimen, that is, including the sum of its own vibration response and the base excitation displacement; therefore, the vibration response of the composite beam specimen under the foundation excitation ⁇ (x, y, t) is expressed as
  • Equation (26) gives the expression of the vibration response of the fiber reinforced composite beam specimen under the basic excitation.
  • the composite can be calculated. The vibration response of the beam test piece at any point;
  • R is the modal order
  • ⁇ f i is the theoretically calculated ith order natural frequency and the experimentally obtained ith order natural frequency difference, The i-th order natural frequency obtained for the experimental test;
  • the material parameters are iterated in a manner of arrangement and combination.
  • the least squares relative error function e fre takes a minimum value
  • the material parameters E 1 , E 2 , G 12 , ⁇ 12 are obtained by preliminary calculation;
  • ⁇ r is the modal damping ratio obtained by the experiment.
  • the strain energies U 1 , U 2 and U 12 of the fiber longitudinal direction, the fiber transverse direction and the shear direction are respectively
  • the modal loss factor and the loss factor in all directions of the fiber have the following relationship
  • the first three natural frequencies and damping ratios of the composite beam specimens are obtained through experimental tests, and the composite beam specimens are excited to reach the resonance state by the above-mentioned natural frequencies, and then the first three-order resonance amplitudes of the composite beams are obtained by laser non-destructive scanning experiments.
  • the identification process includes:
  • the position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device to be at the free end position of the test piece of the composite beam to be tested; then, the excitation device is turned on, and the sine sweep excitation is performed in a large frequency range.
  • Test monitor the time domain waveform data of the sweep excitation signal according to the acceleration sensor on the vibration platform, and obtain the spectrum of the vibration response signal of the free end position of the test specimen to be tested by the laser vibration measuring instrument, and pass the half power bandwidth
  • the method identifies the frequency corresponding to each peak and the frequency of the half power point, thereby obtaining the first three natural frequencies and damping ratios of the composite beam specimen to be tested;
  • the material parameter vector is constructed with the mean value of the material parameters provided by the manufacturer, and the material parameters are iterated in a way of arrangement and combination.
  • the relative error function of the i-th natural frequency obtained by theoretical calculation and the experimentally obtained i-th natural frequency is obtained.
  • the longitudinal elastic modulus E 1 , the transverse elastic modulus E 2 , the shear modulus G 12 and the Poisson's ratio ⁇ 12 of the fiber can be preliminarily calculated.
  • the experiment is obtained.
  • the results of the third-order modal loss factor of the composite beam specimen were calculated, and the longitudinal loss factor ⁇ 1 , the transverse loss factor ⁇ 2 and the shear loss factor ⁇ 12 were calculated.
  • the first three-order resonance amplitude curve of the composite beam is obtained by accurate test, and the frequency of the excitation device is adjusted to the first-order natural frequency, and the test beam of the composite beam is excited to reach the first-order resonance state, and the vibration is passed.
  • the acceleration sensor on the platform determines the excitation amplitude corresponding to the first-order resonance state; at the same time, the control switch of the laser scanning vibration measuring device is turned on, and the position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device at a certain scanning speed.
  • step size constructs the material parameter vector and iterates the parameters in a permutation and combination manner.
  • a set of resonance amplitude curves corresponding to the first three natural frequencies can be obtained by theoretical calculation; then, a theoretical calculation is obtained.
  • the material parameters used at the time are the final material parameters obtained by the identification.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention discloses a fiber reinforced composite material parameter identifying method based on laser nondestructive scanning. The method comprises the following steps: imputing first three order natural frequency, damping ratio and resonance amplitude curve, and setting an identification permissible error value to be 10%, so that fiber longitudinal elasticity modulus E 1, fiber transverse elasticity modulus E 2, shear modulus G 12, Poisson's ratio ν 21, fiber longitudinal loss factor η 22, fiber transverse loss factor η 22, and shear loss factor η 12 are obtained. Also provided is a device for implementing the method. The method combines vibration transmission with laser displacement sensing technology to achieve a one-time measurement of the mechanical properties of the test material in all directions, reducing the associated workload.

Description

基于激光无损扫描的纤维增强复合材料参数辨识方法及装置Fiber reinforced composite material parameter identification method and device based on laser non-destructive scanning 技术领域Technical field
本发明涉及机械领域,具体说是一种纤维增强复合材料参数辨识方法及装置。The invention relates to the field of machinery, in particular to a fiber reinforced composite material parameter identification method and device.
背景技术Background technique
结构复合材料因其结构的比强度高、比模量高、材料具有可设计性、热稳定性好,而且还具有承载力大、重量轻等优点,广泛应用于航空、航天、船舶、体育器械、电气设备、医学、兵器工业与化学工业等领域。随着现代工业水平的不断提高,很多纤维增强复合薄壁构件,例如复合材料叶片、复合材料整体叶盘以及复合圆柱壳,其经常工作在高速旋转、高温、腐蚀气体冲蚀等恶劣环境中,由此带来的振动疲劳、振动失效问题越来越突出,使其不能正常履行人们对其预先设想的功能。Structural composites are widely used in aviation, aerospace, marine, sports equipment due to their high specific strength, high specific modulus, designability, thermal stability, and high bearing capacity and light weight. , electrical equipment, medicine, weapons industry and chemical industry. With the continuous improvement of the modern industrial level, many fiber reinforced composite thin-walled members, such as composite blades, composite whole leaf discs and composite cylindrical shells, often work in harsh environments such as high-speed rotation, high temperature, corrosive gas erosion, etc. The resulting vibration fatigue and vibration failure problems are becoming more and more prominent, making it impossible to perform the functions that people have previously envisaged.
工业生产上应用纤维增强复合材料时时,常常要求对这类材料的各项材料性能有着清楚的认识和了解。然而,由于复合材料的各向异性特点,在检测其纤维纵向弹性模量,纤维横向弹性模量,剪切模量,泊松比以及纤维纵向损耗因子,纤维横向损耗因子和剪切损耗因子等材料参数时,目前工业上大多通过疲劳测试机进行测量。该方法在诸多方面上有其局限性:成本方面:需要破坏较多的材料样板,在人力物力方面造成大量浪费;测量手段方面:高强度的拉伸对仪器的稳定性有特殊要求,同时测量的准确度难以保证,测量结果存在较大误差;安全方面:测量过程中材料会断裂,在安全上有较大隐患。When fiber-reinforced composite materials are used in industrial production, it is often required to have a clear understanding and understanding of the properties of various materials of such materials. However, due to the anisotropic characteristics of the composite, the fiber longitudinal elastic modulus, fiber transverse elastic modulus, shear modulus, Poisson's ratio and fiber longitudinal loss factor, fiber transverse loss factor and shear loss factor are measured. In the case of material parameters, most of the current industrial measurements are made by fatigue testers. This method has its limitations in many aspects: cost: it needs to destroy more material samples, causing a lot of waste in manpower and material resources; measuring means: high-strength stretching has special requirements on the stability of the instrument, and measurement The accuracy is difficult to guarantee, and the measurement results have large errors. In terms of safety: the material will break during the measurement process, and there is a big hidden danger in safety.
目前,对复合材料无损检测这一方面的研究并不是很深入,但也有部分相关的研究有一定的参考价值。专利CN201610166104.9发明了一种新型的复合材料参数识别仪,通过改变识别仪的结构从而使得识别出的参数更为精确,然而其原理仍是通过压、剪、拉等传统方式测得复合材料的参数,并没有从根本上解决浪费,有一定安全隐患等问题。对纤维增强复合材料的检测方法亟需改善,应用振动及激光无损扫描技术能很好地解决相关问题。At present, the research on non-destructive testing of composite materials is not very deep, but some related researches have certain reference value. Patent CN201610166104.9 invented a new type of composite material parameter identification instrument. By changing the structure of the identification instrument, the identified parameters are more accurate. However, the principle is still to measure the composite material by traditional methods such as pressing, shearing and pulling. The parameters do not fundamentally solve the waste, there are certain security risks and other issues. The detection method of fiber reinforced composites needs to be improved, and the application of vibration and laser non-destructive scanning technology can solve the related problems well.
发明内容Summary of the invention
针对现有技术存在的问题,本发明提供一种基于激光无损扫描的纤维增强复合结构材料参数测试仪及其工作方法。具体技术方案如下:In view of the problems existing in the prior art, the present invention provides a fiber reinforced composite structural material parameter tester based on laser non-destructive scanning and a working method thereof. The specific technical solutions are as follows:
基于激光无损扫描的纤维增强复合材料参数辨识方法,包括以下步骤:A method for identifying a fiber reinforced composite material based on laser non-destructive scanning includes the following steps:
步骤1:打开真空外壳,调节旋紧螺杆,将待测复合梁试件依次安装于夹持机构的基准位置;Step 1: Open the vacuum casing, adjust the tightening screw, and install the test piece of the composite beam to be tested in the reference position of the clamping mechanism in sequence;
步骤2:通过调节旋紧螺杆,逐渐调紧夹持机构的压力,根据压力传感器输出的压力 示数,判断已到达适宜的夹持力度后,停止调节;Step 2: Gradually adjust the pressure of the clamping mechanism by adjusting the tightening screw, according to the pressure output by the pressure sensor After the judgment, it is judged that the appropriate clamping force has been reached, and the adjustment is stopped;
步骤3:调节升降平台到合适位置,拧紧真空外壳;用抽气机将外壳内的气体排出,制造真空环境;Step 3: Adjust the lifting platform to a suitable position, tighten the vacuum casing; exhaust the gas inside the casing with an air extractor to create a vacuum environment;
步骤4:通过激光扫描测振装置移动激光测振仪发出的激光点位置,使其处于待测复合梁试件的自由端位置;然后,开启激振装置,在较大频率范围内进行正弦扫频激振测试,根据活动钳身上加速度传感器监测扫频激振信号的时域波形数据,并通过激光测振仪获得待测复合梁试件的自由端位置的振动响应信号的频谱,通过半功率带宽法辨识每个峰值对应的频率以及半功率点的频率,进而获得待测复合梁试件的前三阶固有频率和阻尼比;Step 4: The position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device to be in the free end position of the composite beam test piece to be tested; then, the excitation device is turned on, and the sine sweep is performed in a large frequency range. The frequency excitation test monitors the time domain waveform data of the frequency sweep excitation signal according to the acceleration sensor on the movable clamp body, and obtains the spectrum of the vibration response signal of the free end position of the composite beam test piece to be tested by the laser vibration meter, and passes the half power The bandwidth method identifies the frequency corresponding to each peak and the frequency of the half power point, thereby obtaining the first three natural frequencies and damping ratios of the composite beam specimen to be tested;
步骤5:调节激振装置的频率为第一阶固有频率,激发测复合梁试件达到第一阶共振状态,通过振动平台上加速度传感器确定第一阶共振状态对应的激励幅度;同时,开启激光扫描测振装置的控制开关,通过激光扫描测振装置移动激光测振仪发出的激光点位置,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出第一阶共振状态下复合梁试件振动幅值随与其长度的变化曲线,称之为第一阶共振振幅曲线,所述第一阶共振振幅曲线横轴为长度,纵轴为振动幅值;Step 5: Adjust the frequency of the excitation device to the first-order natural frequency, and the excitation composite beam test piece reaches the first-order resonance state, and the excitation amplitude corresponding to the first-order resonance state is determined by the acceleration sensor on the vibration platform; meanwhile, the laser is turned on. The control switch of the scanning vibration measuring device moves the position of the laser spot emitted by the laser vibrometer by the laser scanning vibration measuring device, and realizes the scanning test of the test piece of the composite beam from the cantilever end position to the free end position, and obtains the composite beam test to be tested. The amplitude of the vibration response signal at each scanning point position, and then the curve of the vibration amplitude of the composite beam specimen with the length of the first-order resonance state is plotted as a first-order resonance amplitude curve. The horizontal axis of the first-order resonance amplitude curve is the length, and the vertical axis is the vibration amplitude;
步骤6:调节激振装置的频率为第二阶固有频率和第三阶固有频率,安照步骤5采用的步骤,获得第二阶共振振幅曲线和第三阶共振振幅曲线;Step 6: adjusting the frequency of the excitation device to the second-order natural frequency and the third-order natural frequency, and illuminating the steps used in step 5 to obtain a second-order resonance amplitude curve and a third-order resonance amplitude curve;
步骤7:根据纤维增强复合结构材料参数识别方法,输入前三阶固有频率、阻尼比和共振振幅曲线,设置识别误差允许值为10%,便可获得纤维纵向弹性模量E1,纤维横向弹性模量E2,剪切模量G12,泊松比为ν21以及纤维纵向损耗因子η11,纤维横向损耗因子η22和剪切损耗因子η12Step 7: According to the parameter identification method of the fiber reinforced composite structural material, input the first three natural frequencies, the damping ratio and the resonance amplitude curve, and set the allowable value of the recognition error to 10% to obtain the longitudinal elastic modulus E 1 of the fiber and the transverse elasticity of the fiber. The modulus E 2 , the shear modulus G 12 , the Poisson's ratio is ν 21 and the fiber longitudinal loss factor η 11 , the fiber transverse loss factor η 22 and the shear loss factor η 12 .
本发明进一步公开了一种实现所述方法的装置,主要是由复式工作平台、升降平台、夹持机构、激振装置、激光扫描测振系统及真空装置组成;复式工作平台包括第一工作平台和第二工作平台,两平台间由“工”字支撑钢连接,下方由底座支撑;夹持机构用于固定及监测待测复合梁试件,包括基准平台及其上部的压块,所述压块由螺栓驱动下压,将待测复合梁试件固定;基准平台上设置了两个圆孔结构,用于放置压力传感器的圆形压力传感垫片,根据压力传感器的数值来定量调节作用于被测复合材料的压力,实现约束条件的量化表征测试;激振装置由激振器和振动平台两部分组成,其作用是产生激励力并向被测复合梁试件传递振动;激振器与信号源连接,产生一定频率的振动,振动通过振动平台 和夹持机构传递给待测梁试件;为实现对激励振动状态的监测,在夹持机构上部装有加速度传感器,用于测量振动幅度的大小;激振器通过螺栓连接固定于第一工作平台上,用以产生激励力,并将其施加到振动平台上;激振平台下方有四个均匀分布在四角的减震螺杆,减震弹簧可套进减震螺杆内,实现与第二工作平台相连,两者之间只能产生弹簧轴线方向上的位移;The invention further discloses a device for realizing the method, which is mainly composed of a double working platform, a lifting platform, a clamping mechanism, an excitation device, a laser scanning vibration measuring system and a vacuum device; the double working platform comprises a first working platform And the second working platform, the two platforms are supported by the "work" word supporting steel, and the bottom is supported by the base; the clamping mechanism is used for fixing and monitoring the composite beam test piece to be tested, including the reference platform and the upper pressing block thereof, The pressure block is driven by bolts to fix the composite beam test piece to be tested; two circular hole structures are arranged on the reference platform for placing the circular pressure sensing gasket of the pressure sensor, and quantitatively adjusting according to the value of the pressure sensor Acting on the pressure of the tested composite material, the quantitative characterization test of the constraint condition is realized; the excitation device consists of two parts: the vibration exciter and the vibration platform, the function of which is to generate the excitation force and transmit the vibration to the test piece of the tested composite beam; Connected to the signal source to generate a certain frequency of vibration, vibration through the vibration platform And the clamping mechanism is transmitted to the test piece of the beam to be tested; in order to realize the monitoring of the state of the excitation vibration, an acceleration sensor is arranged on the upper part of the clamping mechanism for measuring the magnitude of the vibration amplitude; the vibration exciter is fixed to the first work by bolt connection On the platform, the excitation force is generated and applied to the vibration platform; under the excitation platform, there are four shock-absorbing screws evenly distributed at the four corners, and the shock-absorbing spring can be inserted into the shock-absorbing screw to realize the second work. The platforms are connected, and only the displacement in the direction of the spring axis can be generated between the two;
激光扫描测振系统包括激光扫描测振仪、丝杠、联轴器以及电动机;承载激光扫描测振仪的平台被一组丝杠贯穿,丝杠与联轴器相连接,构成传动机构;在电动机的动力下,激光扫描测振仪能够通过传动机构实现激光测振装置水平方向的运动;激光测振仪发出的激光点位置,在一定的扫描速度下,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出相应的共振振幅曲线;激光扫描测振系统设置在升降平台上可实现激光扫描测振仪在竖直方向上位置的调整,所述升起平台的升降架固定在第二工作平台上。The laser scanning vibration measuring system comprises a laser scanning vibrometer, a lead screw, a coupling and an electric motor; the platform carrying the laser scanning vibrometer is penetrated by a set of lead screws, and the lead screw is connected with the coupling to form a transmission mechanism; Under the power of the motor, the laser scanning vibrometer can realize the horizontal movement of the laser vibration measuring device through the transmission mechanism; the position of the laser spot emitted by the laser vibrometer can realize the test specimen of the composite beam from the cantilever at a certain scanning speed. The scanning test from the end position to the free end position obtains the amplitude of the vibration response signal of the composite beam test piece at each scanning measuring point position, and then draws a corresponding resonance amplitude curve; the laser scanning vibration measuring system is arranged on the lifting platform The adjustment of the position of the laser scanning vibrometer in the vertical direction can be realized, and the lifting frame of the lifting platform is fixed on the second working platform.
在测试仪的外部设有真空装置,由壳体和挡板组成;挡板与壳体紧密连接,防止空气进入;在挡板上有圆孔结构,与气泵相连,从而抽出装置内气体,创造真空环境,减小实验误差。A vacuum device is arranged outside the tester, which is composed of a casing and a baffle; the baffle is tightly connected with the casing to prevent air from entering; and the baffle has a circular hole structure connected to the air pump to extract the gas inside the device, thereby creating Vacuum environment reduces experimental error.
本发明的优点是:本发明提供了一种基于激光无损扫描的纤维增强复合结构材料参数测试仪。首先设计了振动传递与激光位移传感技术结合的参数测试方法,打破了传统通过拉伸等物理方法进行测量的模式,并实现了一次性测量待测材料在各个方向的力学性能的目标,极大地减少了相关工作的工作量;其次,仪器在测试的整个过程不会对材料产生任何影响,节约了原材料,减少了经济损失,更加绿色环保;另外,外部加装真空壳体的设计,使得测试过程是在真空下进行的,避免了空气阻尼的影响,测量结果更准确;最后,本作品便于组装拆卸,而且各部件的质量均不大,可装入拉箱,具有较好的便携性。可根据实际情况在不同地点测量复合材料的参数,为生产测量提供了极大的便利。An advantage of the present invention is that the present invention provides a fiber reinforced composite structural material parameter tester based on laser non-destructive scanning. Firstly, the parameter test method combining vibration transmission and laser displacement sensing technology was designed, which broke the traditional measurement mode by stretching and other physical methods, and achieved the goal of measuring the mechanical properties of the material to be tested in all directions at one time. The earth reduces the workload of related work; secondly, the instrument does not have any impact on the material during the whole process of testing, saves raw materials, reduces economic losses, and is more environmentally friendly; in addition, the external vacuum housing design makes The test process is carried out under vacuum, avoiding the influence of air damping, and the measurement result is more accurate. Finally, the work is easy to assemble and disassemble, and the quality of each component is not large, and can be loaded into the pull box for better portability. . The parameters of the composite materials can be measured at different locations according to actual conditions, which provides great convenience for production measurement.
附图说明DRAWINGS
图1是本发明具体实施方式的纤维增强复合结构材料参数测试仪的正视图;1 is a front elevational view of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention;
图2是本发明具体实施方式的纤维增强复合结构材料参数测试仪的俯视图;2 is a top plan view of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention;
图3是本发明具体实施方式的纤维增强复合结构材料参数测试仪的结构简;3 is a schematic structural diagram of a fiber reinforced composite structural material parameter tester according to an embodiment of the present invention;
图4是图1基础激励下纤维增强复合梁试件的理论模型;Figure 4 is a theoretical model of the fiber reinforced composite beam specimen under the basic excitation of Figure 1;
图5是理论和实验获得复合梁第一阶共振振幅曲线及其误差上下限图;Figure 5 is a theoretical and experimental obtained first-order resonance amplitude curve of the composite beam and its error upper and lower limits;
图中:1-底座,2-第二工作平台,3-“工”字支撑钢,4-第一工作平台,5-振动平台, 6-减震弹簧,7-挡板,8-夹持机构,9-加速度传感器,10-丝杠,11-激光扫描测振仪,12-壳体,13-滑块,14-升降平台,15-待测梁试件,16-激振器,17-联轴器,18-电动机,19-可伸缩平台,20-排气孔,21-基准平台。In the picture: 1-base, 2-second working platform, 3-"gong" word support steel, 4-first working platform, 5-vibration platform, 6-damper spring, 7-baffle, 8-clamping mechanism, 9-acceleration sensor, 10-screw, 11-laser scanning vibrometer, 12-shell, 13-slider, 14-lift platform, 15-beam specimen to be tested, 16-vibrator, 17-coupling, 18-motor, 19-retractable platform, 20-venting, 21-reference platform.
具体实施方式detailed description
下面结合附图具体说明本发明。The invention will be specifically described below with reference to the accompanying drawings.
1.测试仪硬件结构设计1. Tester hardware structure design
如图1-图5所示,本发明的装置主要是由复式工作平台、升降平台、夹持机构、激振装置、激光扫描测振系统及真空装置组成;As shown in FIG. 1 to FIG. 5, the device of the present invention mainly comprises a double working platform, a lifting platform, a clamping mechanism, an excitation device, a laser scanning vibration measuring system and a vacuum device;
复式工作平台包括第一工作平台4和第二工作平台2,两平台间由“工”字支撑钢3连接,下方由底座1支撑;夹持机构8用于固定及监测待测复合梁试件15,包括基准平台17及其上部的压块18,所述压块由螺栓驱动下压,将待测复合梁试件固定;The compound working platform comprises a first working platform 4 and a second working platform 2, wherein the two platforms are connected by the "work" word supporting steel 3, and the lower part is supported by the base 1; the clamping mechanism 8 is used for fixing and monitoring the composite beam test piece to be tested 15, comprising a reference platform 17 and an upper pressing block 18, the pressing block is driven by a bolt to press the composite beam test piece to be tested;
基准平台17上表面设置了两个圆孔结构,用于放置压力传感器的圆形压力传感垫片,根据压力传感器的数值来定量调节作用于被测复合材料的压力,实现约束条件的量化表征测试。The upper surface of the reference platform 17 is provided with two circular hole structures for placing a circular pressure sensing gasket of the pressure sensor, and quantitatively adjusting the pressure acting on the composite material according to the value of the pressure sensor to realize quantitative characterization of the constraint condition. test.
激振装置由激振器16和振动平台5两部分组成,其作用是产生激励力并向被测复合梁试件传递振动;激振器与信号源连接,产生一定频率的振动,振动通过振动平台5和夹持机构8传递给待测梁试件15;为实现对激励振动状态的监测,在夹持机构8上部装有加速度传感器9,用于测量振动幅度的大小。激振器16通过螺栓连接固定于第一工作平台4上,用以产生激励力,并将其施加到振动平台5上;激振平台5下方有四个均匀分布在四角的减震螺杆,减震弹簧6可套进减震螺杆内,实现与第二工作平台2相连,两者之间只能产生弹簧轴线方向上的位移。The excitation device is composed of two parts, the vibration exciter 16 and the vibration platform 5, and the function is to generate an excitation force and transmit vibration to the test piece of the composite beam to be tested; the vibration exciter is connected with the signal source to generate vibration of a certain frequency, and the vibration passes through the vibration. The platform 5 and the clamping mechanism 8 are transmitted to the beam test piece 15 to be tested; in order to realize the monitoring of the excitation vibration state, an acceleration sensor 9 is mounted on the upper portion of the clamping mechanism 8 for measuring the magnitude of the vibration amplitude. The vibration exciter 16 is fixed to the first working platform 4 by bolt connection for generating an excitation force and applying it to the vibration platform 5; under the excitation platform 5, there are four shock-absorbing screws uniformly distributed at the four corners, The shock spring 6 can be inserted into the shock absorbing screw to be connected with the second working platform 2, and only the displacement in the direction of the spring axis can be generated between the two.
激光扫描测振系统包括激光扫描测振仪11、丝杠10、联轴器以及电动机。承载激光扫描测振仪的平台被一组丝杠贯穿,丝杠与联轴器相连接,构成传动机构;在电动机的动力下,激光扫描测振仪能够通过传动机构实现激光测振装置水平方向的运动。激光测振仪发出的激光点位置,在一定的扫描速度下,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出相应的共振振幅曲线。激光扫描测振系统设置在升降平台14上可实现激光扫描测振仪在竖直方向上位置的调整,所述升起平台14的升降架固定在第二工作平台2上,两根交叉铰接的支撑梁,一根固定在工作平台2上,另外一根与滑块13轴接,滑块通过丝杠与电机连接;在电机驱动下,丝杠带动滑块13在工作平台上滑动进而带动升降平台上下移动; The laser scanning vibration measuring system includes a laser scanning vibrometer 11, a lead screw 10, a coupling, and an electric motor. The platform carrying the laser scanning vibrometer is penetrated by a set of lead screws, and the lead screw is connected with the coupling to form a transmission mechanism; under the power of the motor, the laser scanning vibrometer can realize the horizontal direction of the laser vibration measuring device through the transmission mechanism. exercise. The laser spot position of the laser vibrometer is used to scan the test specimen from the cantilever end position to the free end position at a certain scanning speed, and obtain the position of the composite beam test piece to be tested at each scanning point. The vibration responds to the amplitude of the signal, which in turn plots the corresponding resonance amplitude curve. The laser scanning vibration measuring system is disposed on the lifting platform 14 to adjust the position of the laser scanning vibrometer in the vertical direction. The lifting platform of the lifting platform 14 is fixed on the second working platform 2, and the two cross-hinged Supporting beams, one fixed on the working platform 2, the other one is connected with the sliding block 13, and the sliding block is connected with the motor through the screw; under the driving of the motor, the screw drives the sliding block 13 to slide on the working platform to drive the lifting The platform moves up and down;
在测试仪的外部设有真空装置,由壳体和挡板组成。挡板可与壳体紧密连接,防止空气进入;在挡板上有圆孔结构,可与气泵相连,从而抽出装置内气体,创造真空环境,减小实验误差。A vacuum device is provided on the outside of the tester, which is composed of a casing and a baffle. The baffle can be tightly connected to the casing to prevent air from entering; the baffle has a circular hole structure that can be connected to the air pump to extract gas from the device, create a vacuum environment, and reduce experimental errors.
2.所述纤维增强复合结构材料参数识别方法是:2. The fiber reinforced composite structural material parameter identification method is:
2.1纤维增强复合梁试件固有特性和基础激励下振动响应求解;2.1 The inherent characteristics of the fiber reinforced composite beam specimen and the vibration response under the fundamental excitation;
复合梁试件是由n层具有正交各向异性特点的纤维和基体材料组合而成的;假设各层之间是牢固粘结的,层间无滑移,无相对位移,固可以不考虑层间耦合效应的影响;首先,将其中面作为参考平面,并建立xoy坐标系;假设纤维方向与整体坐标系x轴方向的夹角为θ,板长为a,板宽为b,板厚为h,每一层位于z坐标轴较低表面hk-1和较高表面hk之间,每层的厚度均相同;图中的1代表纤维纵向,2代表纤维横向,3代表垂直于1-2平面的方向;The composite beam specimen is composed of n layers of fibers with orthogonal anisotropy characteristics and matrix materials; assuming that the layers are firmly bonded, there is no slip between the layers, no relative displacement, solid can be ignored The effect of the interlayer coupling effect; first, the midplane is used as the reference plane, and the xoy coordinate system is established; the angle between the fiber direction and the x-axis direction of the global coordinate system is θ, the plate length is a, the plate width is b, and the plate thickness is For h, each layer is located between the lower surface h k-1 of the z coordinate axis and the higher surface h k , and the thickness of each layer is the same; 1 in the figure represents the longitudinal direction of the fiber, 2 represents the transverse direction of the fiber, and 3 represents the vertical direction. 1-2 plane direction;
假设复合梁试件受到基础激励载荷的影响,且该基础激励的运动表达式为y(t)=Yeiωt(1)It is assumed that the composite beam specimen is affected by the base excitation load, and the motion expression of the base excitation is y(t)=Ye iωt (1)
i表示虚单位i represents a virtual unit
t表示时间t indicates time
其中,Y为激励幅值,ω为激励频率;Where Y is the excitation amplitude and ω is the excitation frequency;
考虑纤维方向的影响,将复合材料的弹性模量表示成如下形式Considering the influence of the fiber direction, the elastic modulus of the composite material is expressed as follows
Figure PCTCN2017078513-appb-000001
Figure PCTCN2017078513-appb-000001
其中,
Figure PCTCN2017078513-appb-000002
分别代表平行纤维方向和垂直纤维方向的复弹性模量,
Figure PCTCN2017078513-appb-000003
代表1-2平面内的复剪切模量,E′1、E′2和G′12分别复弹性模量
Figure PCTCN2017078513-appb-000004
和复剪切模量
Figure PCTCN2017078513-appb-000005
的实部;并且1方向作用应力引起1、2方向应变的泊松比为ν12,2方向作用应力引起1、2方向应变的泊松比为ν21
among them,
Figure PCTCN2017078513-appb-000002
Relating to the parallel fiber direction and the direction of the vertical fiber, respectively,
Figure PCTCN2017078513-appb-000003
Representing the complex shear modulus in plane 1-2, the complex elastic modulus of E' 1 , E' 2 and G' 12 respectively
Figure PCTCN2017078513-appb-000004
And complex shear modulus
Figure PCTCN2017078513-appb-000005
The real part of the stress; and the Poisson's ratio of the strain in the direction of 1 and 2 is ν 12 , and the Poisson's ratio of the strain in the direction of 2 and 2 is ν 21 ;
方向应变的泊松比为ν21The Poisson's ratio of the directional strain is ν 21 ;
基于经典层合板理论,可将纤维增强复合梁试件的位移场写为如下形式Based on the classical laminate theory, the displacement field of the fiber reinforced composite beam specimen can be written as follows
Figure PCTCN2017078513-appb-000006
Figure PCTCN2017078513-appb-000006
Figure PCTCN2017078513-appb-000007
Figure PCTCN2017078513-appb-000007
w(x,y,z,t)=w0(x,y,t) w(x,y,z,t)=w 0 (x,y,t)
其中,z表示在z轴方向上的位移;u,v,w代表板内任意一点的位移;u0,v0,w0代表板中面位移;h为复合梁试件的厚度;t表示时间;Where z represents the displacement in the z-axis direction; u, v, w represents the displacement of any point in the plate; u 0 , v 0 , w 0 represents the surface displacement of the plate; h is the thickness of the composite beam test piece; t represents time;
由于是对称层合板,面内振动和横向振动不存在耦合,因此只需考虑薄板的横向振动,即忽略中面位移u0和v0;根据经典层合板理论的假设可知,正应变εz和剪应变γyz、γxz都为0,即εz=γyz=γxz=0,由应变和位移的关系,板内任意一点的应变可以表示为Since it is a symmetrical laminate, there is no coupling between in-plane vibration and lateral vibration. Therefore, it is only necessary to consider the lateral vibration of the thin plate, that is, ignore the mid-surface displacements u 0 and v 0 ; according to the assumption of the classical laminate theory, the positive strain ε z and The shear strains γ yz and γ xz are all 0, that is, ε z = γ yz = γ xz =0. From the relationship between strain and displacement, the strain at any point in the plate can be expressed as
Figure PCTCN2017078513-appb-000008
Figure PCTCN2017078513-appb-000008
Figure PCTCN2017078513-appb-000009
Figure PCTCN2017078513-appb-000009
Figure PCTCN2017078513-appb-000010
Figure PCTCN2017078513-appb-000010
εx表示该点在x方向上的正应变ε x represents the positive strain of the point in the x direction
εy表示该点在y方向上的正应变ε y represents the positive strain of the point in the y direction
εxy表示该点在xy平面上的剪应变ε xy indicates the shear strain of the point on the xy plane
薄板中面弯曲挠曲率和扭曲率可表示为The bending curvature and the twist rate of the middle surface of the thin plate can be expressed as
Figure PCTCN2017078513-appb-000011
Figure PCTCN2017078513-appb-000011
kx表示该点在x方向上的挠曲率k x represents the flexural curvature of the point in the x direction
ky表示该点在y方向上的挠曲率k y represents the flexural curvature of the point in the y direction
kxy表示该点在xy平面上的扭曲率k xy represents the distortion rate of the point on the xy plane
which is
εx=zκx,εy=zκy,γxy=zκxy ε x = zκ x , ε y = zκ y , γ xy = zκ xy
对于正交各向异性材料,材料主轴方向的应力-应变关系为For orthotropic materials, the stress-strain relationship in the direction of the major axis of the material is
Figure PCTCN2017078513-appb-000012
Figure PCTCN2017078513-appb-000012
1表示纤维纵向,2表示纤维横向,6表示垂直板面方向,Q*表示弹性系数 1 indicates the longitudinal direction of the fiber, 2 indicates the transverse direction of the fiber, 6 indicates the direction of the vertical plate surface, and Q* indicates the elastic modulus.
其中,among them,
Figure PCTCN2017078513-appb-000013
Figure PCTCN2017078513-appb-000013
Figure PCTCN2017078513-appb-000014
Figure PCTCN2017078513-appb-000014
当材料主轴方向与整体坐标系之间有一定夹角θ时,用应力-应变转轴公式计算得到第k层板在整体坐标系下的应力-应变关系如下When there is a certain angle θ between the main axis direction of the material and the global coordinate system, the stress-strain relationship of the k-th layer in the global coordinate system is calculated by the stress-strain rotation formula.
Figure PCTCN2017078513-appb-000015
Figure PCTCN2017078513-appb-000015
其中,among them,
Figure PCTCN2017078513-appb-000016
Figure PCTCN2017078513-appb-000016
Figure PCTCN2017078513-appb-000017
Figure PCTCN2017078513-appb-000017
Figure PCTCN2017078513-appb-000018
Figure PCTCN2017078513-appb-000018
Figure PCTCN2017078513-appb-000019
Figure PCTCN2017078513-appb-000019
Figure PCTCN2017078513-appb-000020
Figure PCTCN2017078513-appb-000020
Figure PCTCN2017078513-appb-000021
Figure PCTCN2017078513-appb-000021
其中,k表示复合梁试件的第k层,θk表示第k层板的纤维方向与整体坐标系x轴的夹角;Where k represents the k-th layer of the composite beam test piece, and θ k represents the angle between the fiber direction of the k-th layer plate and the x-axis of the global coordinate system;
薄板所受弯矩和扭矩为The bending moment and torque of the thin plate are
Figure PCTCN2017078513-appb-000022
Figure PCTCN2017078513-appb-000022
Mx表示x轴方向上的弯矩M x represents the bending moment in the x-axis direction
My表示y轴方向上的弯矩M y represents the bending moment in the y-axis direction
Mxy表示xy平面上的扭矩M xy represents the torque on the xy plane
D*表示弯曲刚度系数 D* represents the bending stiffness coefficient
其中,among them,
Figure PCTCN2017078513-appb-000023
Figure PCTCN2017078513-appb-000023
为了便于理论分析建模,将复合梁试件所受基础激励等效为均布惯性力外载In order to facilitate the theoretical analysis and modeling, the basic excitation of the composite beam specimen is equivalent to the uniform inertial force external load.
Figure PCTCN2017078513-appb-000024
Figure PCTCN2017078513-appb-000024
Y表示基础激励的位移幅值Y represents the displacement amplitude of the base excitation
则薄板振动的动能可以用下式表示The kinetic energy of the vibration of the thin plate can be expressed by the following formula
Figure PCTCN2017078513-appb-000025
Figure PCTCN2017078513-appb-000025
其中,ρ为薄板的密度,h为薄板的厚度;Where ρ is the density of the thin plate and h is the thickness of the thin plate;
则薄板弯曲储存的应变能用下式表示The strain energy stored in the bending of the thin plate is expressed by the following formula
Figure PCTCN2017078513-appb-000026
Figure PCTCN2017078513-appb-000026
薄板所受均布惯性力做功为The uniform inertial force of the thin plate is
Wq=∫∫Rq(t)w0dxdy          (13)W q= ∫∫ R q(t)w 0 dxdy (13)
假设薄板横向振动的振动位移可以表示为It is assumed that the vibration displacement of the lateral vibration of the thin plate can be expressed as
w0(x,y,t)=eiωtW(ξ,η)        (14)w 0 (x,y,t)=e iωt W(ξ,η) (14)
其中,ω为薄板振动的圆频率,与激励频率相同,Wij(ξ,η)为振型函数,有如下形式Where ω is the circular frequency of the vibration of the thin plate, and the excitation frequency is the same, W ij (ξ, η) is a vibration mode function, and has the following form
Figure PCTCN2017078513-appb-000027
Figure PCTCN2017078513-appb-000027
其中,aij为待定系数,pi(ξ)(i=1,...,M)和qj(η)(j=1,...,N)为一系列的正交多项式;Where a ij is a pending coefficient, p i (ξ) (i=1, . . . , M) and q j (η) (j=1, . . . , N) are a series of orthogonal polynomials;
通过对满足边界条件的多项式函数进行正交化处理来获得一系列的正交多项式Obtain a series of orthogonal polynomials by orthogonalizing the polynomial functions that satisfy the boundary conditions
P1(ξ)=χ(ξ),P1(η)=κ(η)P 1 (ξ)=χ(ξ), P 1 (η)=κ(η)
P2(ζ)=(ζ-B2)P1(ζ)P 2 (ζ)=(ζ-B 2 )P 1 (ζ)
Pk(ζ)=(ζ-Bk)Pk-1(ζ)-CkPk-2(ζ)P k (ζ)=(ζ-B k )P k-1 (ζ)-C k P k-2 (ζ)
ζ=ξ,η,k>2        (16)ζ=ξ,η,k>2 (16)
其中,Bk和Ck为系数函数,其表达式分别为 Where B k and C k are coefficient functions, and their expressions are respectively
Figure PCTCN2017078513-appb-000028
Figure PCTCN2017078513-appb-000028
Figure PCTCN2017078513-appb-000029
Figure PCTCN2017078513-appb-000029
其中,W(ζ)是权函数,通常取W(ζ)=1;而χ(ξ)和κ(η)是满足固支、简支、自由等边界条件的多项式函数,且具有如下的形式Where W(ζ) is a weight function, usually taking W(ζ)=1; and χ(ξ) and κ(η) are polynomial functions satisfying the boundary conditions of solid support, simply supported, free, etc., and have the following form
χ(ξ)=ξp(1-ξ)q,κ(η)=ηr(1-η)s χ(ξ)=ξ p (1-ξ) q , κ(η)=η r (1-η) s
ξ=x/a,η=y/b           (18)ξ=x/a, η=y/b (18)
由于是悬臂边界条件,所以取p=2,r=0,q=0,s=0;将公式(15)代入公式(11)、(12)和(13)中,可以得到复合梁试件振动的最大动能Tmax、最大应变能Umax和均布惯性力做功最大值Wqmax分别为Since it is a cantilever boundary condition, take p=2, r=0, q=0, s=0; substituting formula (15) into equations (11), (12) and (13), the composite beam specimen can be obtained. The maximum kinetic energy T max , the maximum strain energy U max and the uniform inertial force work maximum value W qmax of the vibration are respectively
Figure PCTCN2017078513-appb-000030
Figure PCTCN2017078513-appb-000030
Figure PCTCN2017078513-appb-000031
Figure PCTCN2017078513-appb-000031
Wqmax=ρhYω2∫∫RWdxdy           (21)W qmax =ρhYω 2 ∫∫ R Wdxdy (21)
定义拉格朗日能量函数L的表达式为The expression defining the Lagrangian energy function L is
L=Tmax+Wqmax-Umax        (22)L=T max +W qmax -U max (22)
通过使能量函数L对待定系数aij的偏导数等于零,即By making the energy function L the partial derivative of the coefficient a ij to be equal to zero, ie
Figure PCTCN2017078513-appb-000032
Figure PCTCN2017078513-appb-000032
可以获得M×N个非齐次线性代数方程,为了求解方便,将其写为矩阵形式为M×N non-homogeneous linear algebraic equations can be obtained. For the convenience of solving, write it as a matrix form.
(K+iC-ω2M)a=F      (24)(K+iC-ω 2 M)a=F (24)
其中,K、C和M分别为结构刚度矩阵、材料阻尼矩阵和结构质量矩阵,广义位移向量a=(a11,a12,…aij)T,F为激振力向量;Where K, C and M are structural stiffness matrix, material damping matrix and structural mass matrix, respectively, generalized displacement vector a = (a 11 , a 12 , ... a ij ) T , F is the excitation force vector;
对于复合梁试件的自由振动问题,只需令材料阻尼矩阵C和激振力向量F为零,即 For the free vibration problem of composite beam specimens, it is only necessary to make the material damping matrix C and the excitation force vector F zero, ie
(K-ω2M)a=0        (25)(K-ω 2 M)a=0 (25)
由式(25)即可求得复合梁试件的固有频率和模态振型;进一步,假设纤维增强复合梁试件受到基础激励载荷作用,并求解基础激励下其振动响应λ(x,y,t)的表达式;考虑到实验测试获得的是复合梁试件绝对振动响应,即包括其自身的振动响应与基础激励位移之和;因此,可将基础激励下复合梁试件的振动响应λ(x,y,t)表示为From the formula (25), the natural frequency and mode shape of the composite beam specimen can be obtained; further, it is assumed that the fiber reinforced composite beam specimen is subjected to the basic excitation load and the vibration response λ (x, y) is solved under the fundamental excitation. , t) expression; taking into account the experimental test to obtain the absolute vibration response of the composite beam specimen, that is, including the sum of its own vibration response and the base excitation displacement; therefore, the vibration response of the composite beam specimen under the foundation excitation λ(x, y, t) is expressed as
λ(x,y,t)=y(t)+w0(x,y,t)          (26)λ(x,y,t)=y(t)+w 0 (x,y,t) (26)
式(26)即给出了基础激励下纤维增强复合梁试件振动响应的求解表达式,在明确基础激励表达式(1)和薄板振动响应表达式(14)的情况下,可以计算获得复合梁试件任意一点的振动响应;Equation (26) gives the expression of the vibration response of the fiber reinforced composite beam specimen under the basic excitation. In the case of clarifying the basic excitation expression (1) and the thin plate vibration response expression (14), the composite can be calculated. The vibration response of the beam test piece at any point;
2.2基于前三阶固有频率和阻尼测试结果初步计算获得复合材料参数2.2 Based on the first three natural frequencies and damping test results, preliminary calculation of composite parameters
首先,以厂商提供的材料参数均值
Figure PCTCN2017078513-appb-000033
为中心,考虑Rerr=10%~20%误差,给出材料参数的取值范围如下
First, the average value of the material parameters provided by the manufacturer
Figure PCTCN2017078513-appb-000033
For the center, consider R err = 10% ~ 20% error, give the range of material parameters as follows
Figure PCTCN2017078513-appb-000034
Figure PCTCN2017078513-appb-000034
Figure PCTCN2017078513-appb-000035
Figure PCTCN2017078513-appb-000035
Figure PCTCN2017078513-appb-000036
Figure PCTCN2017078513-appb-000036
Figure PCTCN2017078513-appb-000037
Figure PCTCN2017078513-appb-000037
在各材料参数取值范围内选取合适的步长构造材料参数向量E1,E2,G1212,具体表达式为Select the appropriate step size material parameter vector E 1 , E 2 , G 12 , ν 12 within the range of values of each material parameter. The specific expression is
Figure PCTCN2017078513-appb-000038
Figure PCTCN2017078513-appb-000038
Figure PCTCN2017078513-appb-000039
Figure PCTCN2017078513-appb-000039
Figure PCTCN2017078513-appb-000040
Figure PCTCN2017078513-appb-000040
Figure PCTCN2017078513-appb-000041
Figure PCTCN2017078513-appb-000041
在分别获取理论与实验固有频率的基础上,基于最小二乘法构造频率相对误差函数efreBased on the theoretical and experimental natural frequencies, respectively, the frequency relative error function e fre is constructed based on the least squares method:
Figure PCTCN2017078513-appb-000042
Figure PCTCN2017078513-appb-000042
其中,R为模态阶数,Δfi为理论计算获得的第i阶固有频率与实验获得的第i阶固有频率差值,
Figure PCTCN2017078513-appb-000043
为实验测试获得的第i阶固有频率;
Where R is the modal order, Δf i is the theoretically calculated ith order natural frequency and the experimentally obtained ith order natural frequency difference,
Figure PCTCN2017078513-appb-000043
The i-th order natural frequency obtained for the experimental test;
以排列组合的方式对材料参数进行迭代,当最小二乘相对误差函数efre取得最小值时,即初步计算获得材料参数E1,E2,G1212The material parameters are iterated in a manner of arrangement and combination. When the least squares relative error function e fre takes a minimum value, the material parameters E 1 , E 2 , G 12 , ν 12 are obtained by preliminary calculation;
然后,由阻尼比和损耗因子的关系可以获得各阶模态损耗因子ηrThen, the relationship between the damping ratio and the loss factor can be used to obtain the modal loss factor η r
ηr=2ζr        (42)η r =2ζ r (42)
其中,ζr为实验获得的模态阻尼比;Where ζ r is the modal damping ratio obtained by the experiment;
根据模态应变能法,纤维纵向、纤维横向和剪切方向的应变能U1、U2和U12分别为According to the modal strain energy method, the strain energies U 1 , U 2 and U 12 of the fiber longitudinal direction, the fiber transverse direction and the shear direction are respectively
Figure PCTCN2017078513-appb-000044
Figure PCTCN2017078513-appb-000044
Figure PCTCN2017078513-appb-000045
Figure PCTCN2017078513-appb-000045
Figure PCTCN2017078513-appb-000046
Figure PCTCN2017078513-appb-000046
模态损耗因子和纤维各个方向的损耗因子有如下关系The modal loss factor and the loss factor in all directions of the fiber have the following relationship
Figure PCTCN2017078513-appb-000047
Figure PCTCN2017078513-appb-000047
其中,U为复合梁试件总的应变能;Where U is the total strain energy of the composite beam test piece;
由式(44)可知,只要通过实验获得复合梁试件前3阶模态阻尼比结果,就可以初步确定纤维纵向、纤维横向和剪切方向的损耗因子η1、η2和η12It can be known from equation (44) that as long as the results of the third-order modal damping ratio of the composite beam specimen are obtained experimentally, the loss factors η 1 , η 2 and η 12 of the fiber longitudinal direction, the fiber transverse direction and the shear direction can be initially determined;
如此,便可通过前三阶固有频率和阻尼测试结果,初步计算获得E1,E2,G1212η1、η2和η12等7个材料参数;In this way, seven material parameters such as E 1 , E 2 , G 12 , ν 12 η 1 , η 2 and η 12 can be obtained by preliminary calculation of the first three natural frequencies and damping test results;
2.3基于激光无损扫描实验数据准确辨识复合材料参数2.3 Based on laser non-destructive scanning experimental data to accurately identify composite parameters
首先,通过实验测试获得复合梁试件的前三阶固有频率和阻尼比,并以上述固有频率激发复合梁试件达到共振状态,然后通过激光无损扫描实验分别获得复合梁的前三阶共振振幅曲线;然后,在初步确定纤维纵向弹性模量E1、纤维横向弹性模量E2、剪切模量G12、泊松比ν12、纤维纵向损耗因子η1、纤维横向损耗因子η2和剪切损耗因子η12的基础上,考虑更小的误差范围(例如10%),以更小的步长构造材料参数向量,并以排列组合的方式对参数进行迭代,每迭代一次,便可通过理论计算获得一组前三阶固有频率对应的共振振幅曲线;Firstly, the first three natural frequencies and damping ratios of the composite beam specimens are obtained through experimental tests, and the composite beam specimens are excited to reach the resonance state by the above-mentioned natural frequencies, and then the first three-order resonance amplitudes of the composite beams are obtained by laser non-destructive scanning experiments. Curve; then, preliminary determination of fiber longitudinal elastic modulus E 1 , fiber transverse elastic modulus E 2 , shear modulus G 12 , Poisson's ratio ν 12 , fiber longitudinal loss factor η 1 , fiber transverse loss factor η 2 and Based on the shear loss factor η 12 , consider a smaller error range (for example, 10%), construct the material parameter vector in a smaller step size, and iterate the parameters in a permutation and combination manner. A set of resonance amplitude curves corresponding to a set of first three natural frequencies is obtained by theoretical calculation;
最后,对比某次理论计算获得的前三阶固有频率对应的共振振幅曲线与测试获得的前 三阶共振振幅曲线的偏差;以对比第一阶共振振幅曲线为例,当理论计算获得的曲线位于误差上下限曲线范围内时,即可认为上述材料参数是准确的,此时所采用的材料参数即为辨识获得的最终材料参数。Finally, compare the resonance amplitude curve corresponding to the first three natural frequencies obtained by a theoretical calculation with the test obtained before The deviation of the third-order resonance amplitude curve; taking the first-order resonance amplitude curve as an example, when the theoretically calculated curve is within the range of the upper and lower error curves, the material parameters are considered to be accurate, and the material used at this time is used. The parameter is the final material parameter obtained by the identification.
辨识流程包括:The identification process includes:
(1)对纤维增强复合梁试件的固有特性和振动响应进行求解(1) Solving the inherent characteristics and vibration response of fiber reinforced composite beam specimens
首先,基于经典层合板理论建立了纤维增强复合梁试件的理论模型,并将其材料参数表示成复模量的形式;然后,基于Ritz能量法求解了复合梁试件的固有特性以及基础激励下复合梁试件的振动响应;Firstly, based on the classical laminated plate theory, the theoretical model of the fiber reinforced composite beam specimen is established, and its material parameters are expressed in the form of complex modulus. Then, the inherent characteristics and basic excitation of the composite beam specimen are solved based on the Ritz energy method. The vibration response of the lower composite beam specimen;
(2)测试获得复合梁试件的前三阶固有频率和阻尼比(2) Test the first three natural frequencies and damping ratios of the composite beam specimens
通过激光扫描测振装置移动激光测振仪发出的激光点位置,使其处于待测复合梁试件的自由端位置;然后,开启激振装置,在较大频率范围内进行正弦扫频激振测试,根据振动平台上加速度传感器监测扫频激振信号的时域波形数据,并通过激光测振仪件测试获得待测复合梁试件的自由端位置的振动响应信号的频谱,通过半功率带宽法辨识每个峰值对应的频率以及半功率点的频率,进而获得待测复合梁试件的前三阶固有频率和阻尼比;The position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device to be at the free end position of the test piece of the composite beam to be tested; then, the excitation device is turned on, and the sine sweep excitation is performed in a large frequency range. Test, monitor the time domain waveform data of the sweep excitation signal according to the acceleration sensor on the vibration platform, and obtain the spectrum of the vibration response signal of the free end position of the test specimen to be tested by the laser vibration measuring instrument, and pass the half power bandwidth The method identifies the frequency corresponding to each peak and the frequency of the half power point, thereby obtaining the first three natural frequencies and damping ratios of the composite beam specimen to be tested;
(3)初步计算复合材料参数(3) Preliminary calculation of composite parameters
首先以厂商提供的材料参数均值为中心,构造材料参数向量,以排列组合的方式对材料参数进行迭代,当理论计算获得的第i阶固有频率与实验获得的第i阶固有频率的相对误差函数取最小值时,便可初步计算获得纤维纵向弹性模量E1、纤维横向弹性模量E2、剪切模量G12和泊松比ν12;然后,基于模态应变能法,通过实验获得复合梁试件前3阶模态损耗因子结果,初步计算纤维纵向损耗因子η1、纤维横向损耗因子η2和剪切损耗因子η12Firstly, the material parameter vector is constructed with the mean value of the material parameters provided by the manufacturer, and the material parameters are iterated in a way of arrangement and combination. The relative error function of the i-th natural frequency obtained by theoretical calculation and the experimentally obtained i-th natural frequency is obtained. When taking the minimum value, the longitudinal elastic modulus E 1 , the transverse elastic modulus E 2 , the shear modulus G 12 and the Poisson's ratio ν 12 of the fiber can be preliminarily calculated. Then, based on the modal strain energy method, the experiment is obtained. The results of the third-order modal loss factor of the composite beam specimen were calculated, and the longitudinal loss factor η 1 , the transverse loss factor η 2 and the shear loss factor η 12 were calculated.
(4)通过激光无损扫描实验,准确测试获得复合梁的前三阶共振振幅曲线,调节激振装置的频率为第一阶固有频率,激发测复合梁试件达到第一阶共振状态,通过振动平台上加速度传感器确定第一阶共振状态对应的激励幅度;同时,开启激光扫描测振装置的控制开关,通过激光扫描测振装置移动激光测振仪发出的激光点位置,在一定的扫描速度下,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出第一阶共振振幅曲线;重复该步骤,可依次获得第二阶共振振幅曲线和第三阶共振振幅曲线;(4) Through the laser non-destructive scanning experiment, the first three-order resonance amplitude curve of the composite beam is obtained by accurate test, and the frequency of the excitation device is adjusted to the first-order natural frequency, and the test beam of the composite beam is excited to reach the first-order resonance state, and the vibration is passed. The acceleration sensor on the platform determines the excitation amplitude corresponding to the first-order resonance state; at the same time, the control switch of the laser scanning vibration measuring device is turned on, and the position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device at a certain scanning speed. To realize the scanning test of the test piece of the composite beam from the cantilever end position to the free end position, obtain the amplitude of the vibration response signal of the test piece of the composite beam to be tested at each scanning point, and then draw the first-order resonance amplitude curve. Repeating this step, the second-order resonance amplitude curve and the third-order resonance amplitude curve are sequentially obtained;
(5)精确辨识复合材料参数(5) Accurate identification of composite parameters
以步骤3中初步获得的材料参数为中心,考虑更小的误差范围(例如10%),以更小 的步长构造材料参数向量,并以排列组合的方式对参数进行迭代,每迭代一次,便可通过理论计算获得一组前三阶固有频率对应的共振振幅曲线;然后,对比某次理论计算获得的前三阶固有频率对应的共振振幅曲线与测试获得的前三阶共振振幅曲线的偏差;当理论计算获得的曲线位于误差上下限曲线范围内时,即可认为上述材料参数是准确的,此时所采用的材料参数即为辨识获得的最终材料参数。 Centering on the material parameters initially obtained in step 3, consider a smaller error range (eg 10%) to be smaller The step size constructs the material parameter vector and iterates the parameters in a permutation and combination manner. Each iteration, a set of resonance amplitude curves corresponding to the first three natural frequencies can be obtained by theoretical calculation; then, a theoretical calculation is obtained. The deviation between the resonance amplitude curve corresponding to the first three natural frequencies and the first three-order resonance amplitude curve obtained by the test; when the theoretically calculated curve is within the range of the upper and lower limits of the error, the above material parameters are considered to be accurate. The material parameters used at the time are the final material parameters obtained by the identification.

Claims (3)

  1. 基于激光无损扫描的纤维增强复合材料参数辨识方法,其特征在于包括以下步骤:A fiber reinforced composite material parameter identification method based on laser non-destructive scanning, characterized in that the method comprises the following steps:
    步骤1:打开真空外壳,调节旋紧螺杆,将待测复合梁试件依次安装于夹持机构的基准位置;Step 1: Open the vacuum casing, adjust the tightening screw, and install the test piece of the composite beam to be tested in the reference position of the clamping mechanism in sequence;
    步骤2:通过调节旋紧螺杆,逐渐调紧夹持机构的压力,根据压力传感器输出的压力示数,判断已到达适宜的夹持力度后,停止调节;Step 2: gradually adjust the pressure of the clamping mechanism by adjusting the tightening screw, and according to the pressure indication outputted by the pressure sensor, judge that the appropriate clamping force has been reached, and then stop the adjustment;
    步骤3:调节升降平台到合适位置,拧紧真空外壳;用抽气机将外壳内的气体排出,制造真空环境;Step 3: Adjust the lifting platform to a suitable position, tighten the vacuum casing; exhaust the gas inside the casing with an air extractor to create a vacuum environment;
    步骤4:通过激光扫描测振装置移动激光测振仪发出的激光点位置,使其处于待测复合梁试件的自由端位置;然后,开启激振装置,在较大频率范围内进行正弦扫频激振测试,根据活动钳身上加速度传感器监测扫频激振信号的时域波形数据,并通过激光测振仪获得待测复合梁试件的自由端位置的振动响应信号的频谱,通过半功率带宽法辨识每个峰值对应的频率以及半功率点的频率,进而获得待测复合梁试件的前三阶固有频率和阻尼比;Step 4: The position of the laser spot emitted by the laser vibrometer is moved by the laser scanning vibration measuring device to be in the free end position of the composite beam test piece to be tested; then, the excitation device is turned on, and the sine sweep is performed in a large frequency range. The frequency excitation test monitors the time domain waveform data of the frequency sweep excitation signal according to the acceleration sensor on the movable clamp body, and obtains the spectrum of the vibration response signal of the free end position of the composite beam test piece to be tested by the laser vibration meter, and passes the half power The bandwidth method identifies the frequency corresponding to each peak and the frequency of the half power point, thereby obtaining the first three natural frequencies and damping ratios of the composite beam specimen to be tested;
    步骤5:调节激振装置的频率为第一阶固有频率,激发测复合梁试件达到第一阶共振状态,通过振动平台上加速度传感器确定第一阶共振状态对应的激励幅度;同时,开启激光扫描测振装置的控制开关,通过激光扫描测振装置移动激光测振仪发出的激光点位置,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出第一阶共振状态下复合梁试件振动幅值随与其长度的变化曲线,称之为第一阶共振振幅曲线,所述第一阶共振振幅曲线横轴为长度,纵轴为振动幅值;Step 5: Adjust the frequency of the excitation device to the first-order natural frequency, and the excitation composite beam test piece reaches the first-order resonance state, and the excitation amplitude corresponding to the first-order resonance state is determined by the acceleration sensor on the vibration platform; meanwhile, the laser is turned on. The control switch of the scanning vibration measuring device moves the position of the laser spot emitted by the laser vibrometer by the laser scanning vibration measuring device, and realizes the scanning test of the test piece of the composite beam from the cantilever end position to the free end position, and obtains the composite beam test to be tested. The amplitude of the vibration response signal at each scanning point position, and then the curve of the vibration amplitude of the composite beam specimen with the length of the first-order resonance state is plotted as a first-order resonance amplitude curve. The horizontal axis of the first-order resonance amplitude curve is the length, and the vertical axis is the vibration amplitude;
    步骤6:调节激振装置的频率为第二阶固有频率和第三阶固有频率,安照步骤5采用的步骤,获得第二阶共振振幅曲线和第三阶共振振幅曲线;Step 6: adjusting the frequency of the excitation device to the second-order natural frequency and the third-order natural frequency, and illuminating the steps used in step 5 to obtain a second-order resonance amplitude curve and a third-order resonance amplitude curve;
    步骤7:根据纤维增强复合结构材料参数识别方法,输入前三阶固有频率、阻尼比和共振振幅曲线,设置识别误差允许值为10%,便可获得纤维纵向弹性模量E1,纤维横向弹性模量E2,剪切模量G12,泊松比为ν21以及纤维纵向损耗因子η11,纤维横向损耗因子η22和剪切损耗因子η12Step 7: According to the parameter identification method of the fiber reinforced composite structural material, input the first three natural frequencies, the damping ratio and the resonance amplitude curve, and set the allowable value of the recognition error to 10% to obtain the longitudinal elastic modulus E 1 of the fiber and the transverse elasticity of the fiber. The modulus E 2 , the shear modulus G 12 , the Poisson's ratio is ν 21 and the fiber longitudinal loss factor η 11 , the fiber transverse loss factor η 22 and the shear loss factor η 12 .
  2. 一种实现权利要求1所述方法的装置,其特征在于:主要是由复式工作平台、升降平台、夹持机构、激振装置、激光扫描测振系统及真空装置组成;复式工作平台包括第一工作平台和第二工作平台,两平台间由“工”字支撑钢连接,下方由底座支撑;夹持机构 用于固定及监测待测复合梁试件,包括基准平台及其上部的压块,所述压块由螺栓驱动下压,将待测复合梁试件固定;基准平台上设置了两个圆孔结构,用于放置压力传感器的圆形压力传感垫片,根据压力传感器的数值来定量调节作用于被测复合材料的压力,实现约束条件的量化表征测试;激振装置由激振器和振动平台两部分组成,其作用是产生激励力并向被测复合梁试件传递振动;激振器与信号源连接,产生一定频率的振动,振动通过振动平台和夹持机构传递给待测梁试件;为实现对激励振动状态的监测,在夹持机构上部装有加速度传感器,用于测量振动幅度的大小;激振器通过螺栓连接固定于第一工作平台上,用以产生激励力,并将其施加到振动平台上;激振平台下方有四个均匀分布在四角的减震螺杆,减震弹簧可套进减震螺杆内,实现与第二工作平台相连,两者之间只能产生弹簧轴线方向上的位移;An apparatus for implementing the method of claim 1 is characterized in that: mainly consists of a duplex working platform, a lifting platform, a clamping mechanism, an excitation device, a laser scanning vibration measuring system and a vacuum device; the duplex working platform includes the first Working platform and second working platform, the two platforms are supported by the word "work" and the bottom is supported by the base; the clamping mechanism It is used for fixing and monitoring the test piece of the composite beam to be tested, including the reference platform and the upper part of the pressure block, the pressure block is driven by the bolt, and the test piece of the composite beam to be tested is fixed; two circular holes are arranged on the reference platform. Structure, a circular pressure sensing gasket for placing a pressure sensor, quantitatively adjusting the pressure acting on the composite material to be tested according to the value of the pressure sensor, and realizing a quantitative characterization test of the constraint condition; the excitation device is composed of a vibration exciter and vibration The platform is composed of two parts, the function of which is to generate excitation force and transmit vibration to the test piece of the composite beam to be tested; the vibration exciter is connected with the signal source to generate vibration of a certain frequency, and the vibration is transmitted to the beam to be tested through the vibration platform and the clamping mechanism. In order to realize the monitoring of the excitation vibration state, an acceleration sensor is arranged on the upper part of the clamping mechanism for measuring the magnitude of the vibration amplitude; the vibration exciter is fixed on the first working platform by bolt connection for generating the excitation force, and Apply it to the vibration platform; there are four shock-absorbing screws evenly distributed at the four corners under the excitation platform, and the shock-absorbing spring can be inserted into the shock-absorbing screw to realize Working platform is connected between the spring can only be displaced in the axial direction;
    激光扫描测振系统包括激光扫描测振仪、丝杠、联轴器以及电动机;承载激光扫描测振仪的平台被一组丝杠贯穿,丝杠与联轴器相连接,构成传动机构;在电动机的动力下,激光扫描测振仪能够通过传动机构实现激光测振装置水平方向的运动;激光测振仪发出的激光点位置,在一定的扫描速度下,实现对待测复合梁试件从悬臂端位置到自由端位置的扫描测试,获得待测复合梁试件在每个扫描测点位置的振动响应信号幅值,进而绘制出相应的共振振幅曲线;激光扫描测振系统设置在升降平台上可实现激光扫描测振仪在竖直方向上位置的调整,所述升起平台的升降架固定在第二工作平台上。The laser scanning vibration measuring system comprises a laser scanning vibrometer, a lead screw, a coupling and an electric motor; the platform carrying the laser scanning vibrometer is penetrated by a set of lead screws, and the lead screw is connected with the coupling to form a transmission mechanism; Under the power of the motor, the laser scanning vibrometer can realize the horizontal movement of the laser vibration measuring device through the transmission mechanism; the position of the laser spot emitted by the laser vibrometer can realize the test specimen of the composite beam from the cantilever at a certain scanning speed. The scanning test from the end position to the free end position obtains the amplitude of the vibration response signal of the composite beam test piece at each scanning measuring point position, and then draws a corresponding resonance amplitude curve; the laser scanning vibration measuring system is arranged on the lifting platform The adjustment of the position of the laser scanning vibrometer in the vertical direction can be realized, and the lifting frame of the lifting platform is fixed on the second working platform.
  3. 根据权利要求1所述的装置,其特征在于:在测试仪的外部设有真空装置,由壳体和挡板组成;挡板与壳体紧密连接,防止空气进入;在挡板上有圆孔结构,与气泵相连,从而抽出装置内气体,创造真空环境,减小实验误差。 The device according to claim 1, wherein a vacuum device is disposed outside the tester, and is composed of a casing and a baffle; the baffle is tightly connected to the casing to prevent air from entering; and the baffle has a circular hole. The structure is connected to the air pump to extract the gas inside the device, creating a vacuum environment and reducing experimental errors.
PCT/CN2017/078513 2017-03-16 2017-03-29 Fiber reinforced composite material parameter identifying method based on laser nondestructive scanning, and device WO2018165999A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019548923A JP7109798B2 (en) 2017-03-16 2017-03-29 Method and Apparatus for Nondestructive Laser Scanning Fiber Reinforced Composite Parameter Identification

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201710155498.2A CN106950280B (en) 2017-03-16 2017-03-16 Fibre reinforced composites parameter identification method based on the lossless scanning of laser
CN201710155498.2 2017-03-16

Publications (1)

Publication Number Publication Date
WO2018165999A1 true WO2018165999A1 (en) 2018-09-20

Family

ID=59472594

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/078513 WO2018165999A1 (en) 2017-03-16 2017-03-29 Fiber reinforced composite material parameter identifying method based on laser nondestructive scanning, and device

Country Status (3)

Country Link
JP (1) JP7109798B2 (en)
CN (1) CN106950280B (en)
WO (1) WO2018165999A1 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109270169A (en) * 2018-10-19 2019-01-25 中国工程物理研究院总体工程研究所 A kind of fibre reinforced composites elastic parameter rapid assay methods and system
CN110008501A (en) * 2019-01-28 2019-07-12 华北水利水电大学 The modeling method and application of the free vibration model of non-linear rectangular section convexity spring
CN110849973A (en) * 2019-12-03 2020-02-28 上海海事大学 High-frequency vibration system and method for nondestructive testing of micro-cracks on surface layer of small-size component
CN111639457A (en) * 2020-06-23 2020-09-08 广州电力机车有限公司 Design method and test method of electric wheel test tool of mining dump truck
CN112100895A (en) * 2020-10-28 2020-12-18 北京机电工程研究所 Layered model correction method and system for bolt connection combined structure
CN112818288A (en) * 2019-11-18 2021-05-18 中车时代电动汽车股份有限公司 Parameter identification method and system of suspension vibration isolation foundation for road simulation test machine
CN113221061A (en) * 2021-04-29 2021-08-06 河南科技大学 Parabolic-substrate resonant curved bulldozer blade and method for setting soil-contacting curved surface
CN115111991A (en) * 2022-06-01 2022-09-27 华中科技大学 Method and equipment for precisely measuring dynamic fit clearance in narrow space

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107515247B (en) * 2017-07-24 2019-11-22 河海大学 A kind of composite panel part splitting method for detecting position
CN107818222A (en) * 2017-11-01 2018-03-20 东北大学 Heat is shaken fiber composite plate nonlinear kinetics parameter test method and system under environment
CN109000872B (en) * 2018-05-31 2019-10-25 东北大学 Unmanned plane wing damage detection device and method under the dust storm of laser scanning-heat is shaken
CN108760020A (en) * 2018-06-01 2018-11-06 东北大学 Fiber composite thin plate Non-Linear Vibration characterization test method based on vibration measurement with laser
CN109115594B (en) * 2018-09-03 2020-10-23 东北大学 Device and method for testing mechanical property of fiber reinforced composite material during degradation
CN109798973B (en) * 2019-03-14 2020-12-25 浙江润久机械科技有限公司 Method for testing natural frequency of non-contact ultrasonic transducer
CN110133103B (en) * 2019-05-29 2021-04-16 东北大学 Fiber reinforced composite material parameter identification method based on plane pulse sound wave excitation
CN110377940B (en) * 2019-06-10 2020-11-13 大连理工大学 Method for obtaining milling vibration amplitude of CFRP workpiece
CN110231405B (en) * 2019-06-19 2021-12-17 中汽研(天津)汽车工程研究院有限公司 Method for testing elastic modulus and damping loss factor of viscoelastic material
CN112172190B (en) * 2020-09-08 2022-02-01 武汉大学 Nondestructive testing guidance composite material laying process on-line adjustment system and method
CN112611659A (en) * 2020-11-30 2021-04-06 潍柴动力股份有限公司 Method and device for testing material parameters
KR102520394B1 (en) * 2021-06-11 2023-04-11 부경대학교 산학협력단 Apparatus and method for analyzing dynamic mode change of anisotropic materials

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0735668A (en) * 1993-07-16 1995-02-07 Toshiba Tungaloy Co Ltd Ultrasonic vibration type fatigue tester
US7798000B1 (en) * 2005-10-28 2010-09-21 Trustees Of Boston University Non-destructive imaging, characterization or measurement of thin items using laser-generated lamb waves
US20120166135A1 (en) * 2009-08-11 2012-06-28 Mathias Fink Method for Determining a Physical Parameter, Imaging Method, and Device for Implementing Said Method
CN103713048A (en) * 2013-12-10 2014-04-09 同济大学 Ultrasonic field non-contact visualization method for nondestructive inspection and device thereof
CN104483218A (en) * 2014-04-25 2015-04-01 华侨大学 Material fatigue characteristic testing method based on supersonic vibration under microscale
US20150177195A1 (en) * 2012-06-06 2015-06-25 Nec Corporation Structure analyzing device and a structure analyzing method
CN104897268A (en) * 2015-06-08 2015-09-09 东北大学 Laser-scanning-based apparatus and method for modal shape testing of high-grade numerical control machine tool
CN105510435A (en) * 2015-11-27 2016-04-20 航天晨光股份有限公司 Method for detecting defects of metal bellows based on laser vibration

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0785047B2 (en) * 1987-04-06 1995-09-13 大阪府 Measuring method and measuring device for loss coefficient, dynamic elastic coefficient, dynamic shear elastic coefficient and dynamic Poisson's ratio of rod or plate
US5086775A (en) * 1990-11-02 1992-02-11 University Of Rochester Method and apparatus for using Doppler modulation parameters for estimation of vibration amplitude
JPH05302880A (en) * 1992-04-28 1993-11-16 Ishikawajima Harima Heavy Ind Co Ltd Apparatus for fatigue test
JPH063237A (en) * 1992-06-24 1994-01-11 Fujitsu Ltd Measuring device for young's modulus and measuring method for young's modulus
US5533399A (en) 1992-09-30 1996-07-09 Wayne State University Method and apparatus for non-destructive measurement of elastic properties of structural materials
JP3785785B2 (en) * 1998-02-04 2006-06-14 ソニー株式会社 Material property measuring device
JP3374973B2 (en) 2000-11-10 2003-02-10 独立行政法人産業技術総合研究所 Measuring method and measuring device for output pressure of constant section working body
DE102006013770A1 (en) 2006-03-24 2007-09-27 Occlutech Gmbh Occlusion instrument and method for its production
CN102538950A (en) * 2010-12-23 2012-07-04 沈阳黎明航空发动机(集团)有限责任公司 Acoustic testing method for inherent frequency of engine part
JP2013203788A (en) 2012-03-27 2013-10-07 Teijin Ltd Vibration-damping structure
CN103528782B (en) * 2013-10-23 2015-11-04 东北大学 Based on thin-walled workpiece vibration-testing apparatus and the method for piezoelectric actuator
CN110867469B (en) 2014-03-17 2023-12-29 三星显示有限公司 Organic EL display device
JP6413596B2 (en) * 2014-10-10 2018-10-31 横河電機株式会社 Resonance frequency measurement system, resonance frequency measurement method
CN105628496A (en) * 2016-03-22 2016-06-01 北京航空航天大学 Multifunctional shearing test fixture used for composite material parameter identification

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0735668A (en) * 1993-07-16 1995-02-07 Toshiba Tungaloy Co Ltd Ultrasonic vibration type fatigue tester
US7798000B1 (en) * 2005-10-28 2010-09-21 Trustees Of Boston University Non-destructive imaging, characterization or measurement of thin items using laser-generated lamb waves
US20120166135A1 (en) * 2009-08-11 2012-06-28 Mathias Fink Method for Determining a Physical Parameter, Imaging Method, and Device for Implementing Said Method
US20150177195A1 (en) * 2012-06-06 2015-06-25 Nec Corporation Structure analyzing device and a structure analyzing method
CN103713048A (en) * 2013-12-10 2014-04-09 同济大学 Ultrasonic field non-contact visualization method for nondestructive inspection and device thereof
CN104483218A (en) * 2014-04-25 2015-04-01 华侨大学 Material fatigue characteristic testing method based on supersonic vibration under microscale
CN104897268A (en) * 2015-06-08 2015-09-09 东北大学 Laser-scanning-based apparatus and method for modal shape testing of high-grade numerical control machine tool
CN105510435A (en) * 2015-11-27 2016-04-20 航天晨光股份有限公司 Method for detecting defects of metal bellows based on laser vibration

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109270169A (en) * 2018-10-19 2019-01-25 中国工程物理研究院总体工程研究所 A kind of fibre reinforced composites elastic parameter rapid assay methods and system
CN110008501B (en) * 2019-01-28 2023-08-04 华北水利水电大学 Modeling method and application of free vibration model of convex spring with nonlinear rectangular section
CN110008501A (en) * 2019-01-28 2019-07-12 华北水利水电大学 The modeling method and application of the free vibration model of non-linear rectangular section convexity spring
CN112818288B (en) * 2019-11-18 2024-04-02 中车时代电动汽车股份有限公司 Parameter identification method and system for suspension vibration isolation foundation for road simulation testing machine
CN112818288A (en) * 2019-11-18 2021-05-18 中车时代电动汽车股份有限公司 Parameter identification method and system of suspension vibration isolation foundation for road simulation test machine
CN110849973A (en) * 2019-12-03 2020-02-28 上海海事大学 High-frequency vibration system and method for nondestructive testing of micro-cracks on surface layer of small-size component
CN111639457A (en) * 2020-06-23 2020-09-08 广州电力机车有限公司 Design method and test method of electric wheel test tool of mining dump truck
CN111639457B (en) * 2020-06-23 2023-06-09 广州电力机车有限公司 Design method and test method of electric wheel test fixture of mining dump truck
CN112100895B (en) * 2020-10-28 2023-09-26 北京机电工程研究所 Hierarchical model correction method and system for bolt connection combined structure
CN112100895A (en) * 2020-10-28 2020-12-18 北京机电工程研究所 Layered model correction method and system for bolt connection combined structure
CN113221061A (en) * 2021-04-29 2021-08-06 河南科技大学 Parabolic-substrate resonant curved bulldozer blade and method for setting soil-contacting curved surface
CN113221061B (en) * 2021-04-29 2024-05-14 河南科技大学 Parabolic substrate resonant curved surface bulldozer blade and soil contact curved surface setting method
CN115111991A (en) * 2022-06-01 2022-09-27 华中科技大学 Method and equipment for precisely measuring dynamic fit clearance in narrow space

Also Published As

Publication number Publication date
JP2020510837A (en) 2020-04-09
CN106950280A (en) 2017-07-14
CN106950280B (en) 2019-04-12
JP7109798B2 (en) 2022-08-01

Similar Documents

Publication Publication Date Title
WO2018165999A1 (en) Fiber reinforced composite material parameter identifying method based on laser nondestructive scanning, and device
WO2020047734A1 (en) Device and method for testing mechanical properties of fiber reinforced composite material during degradation
CN106442115B (en) Ultrasonic experimental device for rock joint under complex stress and control system thereof
Zapoměl et al. Identification of material damping of a carbon composite bar and study of its effect on attenuation of its transient lateral vibrations
Balkan et al. Nonlinear dynamic behavior of viscoelastic sandwich composite plates under non-uniform blast load: Theory and experiment
Wesolowski et al. Air damping influence on dynamic parameters of laminated composite plates
Vanwalleghem et al. External damping losses in measuring the vibration damping properties in lightly damped specimens using transient time-domain methods
Biscontin et al. Vibrations of steel-concrete composite beams
CN110133103A (en) Fibre reinforced composites parameter identification method based on planar pulse sound wave exciting
Körük Quantification and minimization of sensor effects on modal parameters of lightweight structures
Dey et al. Modal characteristics of cracked thin walled unsymmetrical cross-sectional steel beams curved in plan
Ashcroft et al. Mixed-mode crack growth in bonded composite joints under standard and impact-fatigue loading
Wang et al. Structural-acoustic modeling and analysis of carbon/glass fiber hybrid composite laminates
Xu et al. Application of split Hopkinson tension bar technique to the study of dynamic fracture properties of materials
Hammami et al. Vibration behavior of composite material with two overlapping delaminations
Toygar et al. Vibration and buckling analysis of a curved sandwich composite beam with FEM
Taheri et al. Computational Modelling of Delamination and Disbond in Adhesively Bonded Joints and the Relevant Damage Detection Approaches
Li et al. Structural Design and Characteristic Research of Product Vibration Test Fixture Based on High Frequency Response Characteristics
Song et al. Comparisons of caustics and strain gauge methods for determining mixed mode fracture toughness of PMMA
Zhang et al. A modal displacement unevenness coefficient method for multi-point matrix damage detection in composite laminates
Lopez-Lopez et al. Novel method for determination the dynamic elastic modulus of composite wind turbine blades
Saber et al. Experimental Stress Analysis for Woven Carbon, Glass, and Kevlar Laminate Shells
Bisagni et al. Strain energy method to model composite structures Damping
Katunin et al. Investigation of dynamic behaviour of laminated composite plates under cyclic loading
Dubinskii et al. Study of composite skin-stringer fatigue behavior under broadband random vibration

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17900299

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019548923

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17900299

Country of ref document: EP

Kind code of ref document: A1