Summary of the invention
The purpose of this invention is to provide a kind of method that combines quantitatively analyzing material interface properties with Non-Destructive Testing with finite element analogy.
Behind the compound substance plus load that bonds together, because different matrix binding material and reinforcing material modulus, tensile strength, and the difference of the bonding strength between them, to produce inhomogenous deformation displacement around the interface, this variation and material character, adhesive property have direct relation.
As everyone knows, the laser speckle interferometry measuring technique is a kind ofly rough surface to be can't harm the technology of full-filled optical measurements (Robert K. ell husband etc. writes, " speckle measuring: principle and application ", China Machine Press, 58-80), owing to the speckle interference measuring technique has that structural approach is simple, the noncontact measurement of full field, light requirement is learned smooth surface, strong interference immunity can be used for outstanding advantages such as operating mode scene, be widely used in the measurement of each side such as displacement, strain.When body surface deformed, the displacement of the derivative of its displacement and displacement and space speckle had definite relation, and this determines that relation is to infer the foundation of body surface distortion according to measured speckle movement.Produce speckle field and corresponding phase size reduction mode by the twin-beam speckle interference, as follows:
In the formula, ψ is the various phases at random of mixing rough surface,
Be the one dimension phase change that thing light is produced.Because the laser speckle sequence chart that the modulation of sine function obtains will present with strip form.
In the formula,
For relative displacement changes, λ is an optical wavelength,
Be the displacement increment of loading direction, n is a fringe order, and size is out of shape in how much directly reflecting of striped.
Thus, the inventor uses the variable that lossless detection method records a certain loading procedure of material, and find with finite element method and to change identical correct emulation mode therewith, thereby obtain every physical parameter of exosyndrome material performance, comprise friction factor, stress, strain and counter-force etc. between fiber and matrix.
The method of using Non-Destructive Testing to combine quantitatively analyzing material interface properties with finite element analogy provided by the invention shown in Fig. 1 process flow diagram, in turn includes the following steps:
1) composite sample to be measured is held on the omnipotent puller system, loading is loaded in described testing sample, omnipotent puller system stretch described testing sample and real time record load and stretching displacement, laser speckle lossless detection method or holograph detection method write down the speckle bar graph or the hologram of described composite sample deformation to be measured simultaneously;
2) drawing process of composite sample to be measured application finite element software simulation steps 1), modulus, yield strength and Poisson ratio are set, adjust the coupling friction factor, obtain the different loads choppy in step, the distortion parameter tabulation in described different loads step is derived, the deflection of same coordinate point subtracts each other, be converted into the two-dimensional matrix value, draw the identical potentiometric map of fringe spacing value size with the speckle bar graph, obtain the emulation mode consistent, finish the quantitative test of described composite sample interface performance to be measured with experimental result.
In this method, lossless detection method is laser speckle lossless detection method or holograph detection method.The loading mode of load is that mechanics loads, heat loads or sound loads.The speed that loads load is 0.005mm/min~0.05mm/min, and required loading load is for can produce the load of the little distortion of 0.5~50 μ m, and this load can not damage material interface.
In the step 1), the real time record mode is that sequence is subtracted each other or phase size reduction mode uniformly-spaced; Accordingly, step 2) in, the method that described same position point displacement subtracts each other is subtracted each other or is uniformly-spaced subtracted each other for its corresponding sequence.In addition, for obtaining speckle image more clearly, in the time of can keeping deformation state to stablize, in reference light or thing light, introduce three steps, four steps or five step phase-shift phases by the piezoelectric micromotor displacement driver, change the phase contraposition phase of two relevant wavelength, and then find the solution and obtain the insensitive clear speckle stripe pattern of background, contrast and noise.Three steps, four steps or five step phase shift phase methods are to carry out according to the method in the following document: Self-marking phase-stepping electronic specklepattern interferometry (ESPI) for determining a phase map with least residues, M.J.Huang, Bo-Son Yun, Optics ﹠amp; Laser Technology, 2009,39 (1), 136-148.
Step 2) in, the tabulation of the distortion parameter in the different loads of derivation step is XYZ Work List at random, i.e. X, and the Y row correspond respectively to X, and Y position coordinate value, Z be listed ass corresponding to being out of shape big or small.The deflection of same coordinate position point is subtracted each other, obtain subtracting each other or uniformly-spaced subtract each other the distortion variable with the corresponding sequence of described step 1).For further obtaining potentiometric map, needing at random, Work List is converted into the regular grid matrix.It is to carry out according to the programming of the Kriging algorithm in the following document that X Y Z Work List at random is converted into grid matrix: Statistics and data analysis ingeology, Davis, John C., John Wiley; Sons, Inc.second edition, pp.383.
Composite sample to be measured is made up of matrix material and reinforcing material; Wherein, matrix material is thermoplastic resin or thermoset resin, be selected from polyolefin, polystyrene, polyamide, polyimide, polyester, polyoxymethylene, polycarbonate, poly-(methyl) acrylate, polyphenylene sulfide, polyphenylene oxide, polyethersulfone, polyetheretherketone, epoxy resin, any one in phenolics, polyurethane, silicones, silicon rubber, butadiene-styrene rubber, neoprene, butyl rubber, butadiene-acrylonitrile rubber and the EP rubbers or several arbitrarily potpourri that mixes with arbitrary proportion.Reinforcing material is carbon fiber, glass fibre, boron fibre, ceramic fiber, basalt fibre, carbon fibre precursor or aromatic polyamide fibre.
The length of above-mentioned composite sample to be measured is 50~200mm, and width is 5~50mm, and thickness is 0.5~5mm; The buried depth on described reinforcing material distance composite sample observation to be measured surface does not exceed sample surfaces, is 0.25~2.5mm; The embedment length of reinforcing material accounts for 10%~100% of composite sample total length to be measured; The embedding form of reinforcing material in composite sample to be measured is that fiber list bundle is embedding or lamella is embedding; Fiber list beam diameter is 0.3~0.7mm; Lamella length is 50~200mm, and width is 5~50mm, and thickness does not exceed described composite sample to be measured surface, is 0.2~1mm.
The method of using Non-Destructive Testing to combine quantitatively analyzing material interface properties with finite element analogy provided by the invention combines the high precision of Non-Destructive Testing, real-time advantage with the comprehensive advantage of finite element analysis, have the advantage of two kinds of methods concurrently.This method adopts the change in displacement under the laser speckle detection system whole audience record compound substance loaded condition, omnipotent puller system real time record pulling force size of while, therefore, not only can directly the be enhanced stressing conditions of body of this method, and can the real-time monitored fiber and matrix material carefully see destruction.This method has been owing to adopted the laser speckle optical detecting method, and is highly sensitive, can carry out harmless quantitative to adhesive property and detect, and need not to destroy fully compound interface.This method can be made quantitative detecting analysis to material adhesive, provides reliable foundation for whether compound substance is fit to actual the use.This method quantitative values highly sensitive, that record is reliable, can carry out harmless quantitative to interface performance and detect, and need not to destroy adhesive linkage.
Embodiment
The invention will be further described below in conjunction with specific embodiment, but the present invention is not limited to following examples.
The method of using Non-Destructive Testing to combine quantitatively analyzing material interface properties with finite element analogy provided by the invention, shown in Figure 1 as flow process, concrete implementation step is as follows:
1) white paint is sprayed on the testing sample surface, be held on after the drying on the omnipotent puller system;
2) with steady weighing apparatus speed load sample, omnipotent puller system real time record stretching displacement and load write down the speckle bar graph simultaneously;
Fig. 2 is the laser speckle pick-up unit synoptic diagram under the composite material interface performance loaded condition.Wherein, 1 is laser instrument, and 2 is charge-coupled image sensor (Charge Coupled Device), and 3 are charger (can be load device such as omnipotent puller system or thermal source), and 4 is computing machine, and 5 is the piezoelectric micromotor displacement driver.
3) use finite element software modeling interface drawing process, set modulus, yield strength and Poisson ratio etc., adjust the coupling friction factor, obtain the different loads choppy in step, simultaneously shift value is derived, application computer programs is subtracted each other and is obtained the potentiometric map consistent with the fringe spacing value;
4) adjust the Finite Element Method intrinsic parameters and find correct emulation mode with speckle bar graph experimental result, and obtain reflecting each physical parameter of interface performance, comprise friction factor, stress, strain or counter-force etc. between fiber and matrix.
The detection at embodiment 1, single bundle carbon-fibre reinforced epoxy resin composite material interface
Matrix material is that (the bisphenol A type epoxy resin model is E-51 to epoxy resin, hardening agent low-molecular-weight polyamide resin model is 651, Wuxi City A Erzi chemical industry company limited), strengthening body is single diameter 7 μ m, 1000 of every bundles, the carbon fibre tow of length 150mm (Dalian Xingke Carbon Fiber Co., Ltd).
This detection method is carried out according to following steps:
1) pours epoxy resin and hardening agent into 80 * 20 * 2mm by 10: 3 quality after than mixing
3The tetrafluoro mould in, a branch of carbon fiber is inserted wherein 1mm thickness place.The follow procedure curing that slowly heats up in baking oven: 60 ℃, 3 hours, 70 ℃, 3 hours, 80 ℃, 3 hours, 105 ℃, 3 hours, 120 ℃ kept 3 hours, and obtained golden yellow translucent batten until curing.Sample is taken out from mould, and white paint is sprayed on the surface, treats to be held on the omnipotent puller system after the white paint drying.With the steady weighing apparatus speed of 0.01mm/min this sample that stretches, omnipotent puller system real time record stretching displacement and load, CCD record speckle bar graph simultaneously.
Dotted line is depicted as the load-displacement curve of the used epoxy resin-carbon fibre composite of this embodiment among Fig. 3.Because epoxy resin matter is hard, modulus is big, and tight with the carbon fiber bonding, and load reaches 200N, and load-displacement curve is smoothly stable, does not still have tangible slippage.Laser speckle shown in Fig. 4 (a) uniformly-spaced subtracts each other figure (be 120s interval time) and reflects this phenomenon equally, and incremental deformation reaches unanimity in the constant speed drawing process.
2) use finite element Ansys10.0 software simulation interface tension process, it is 2.8E9 that the epoxy resin elastic modulus is set, and Poisson ratio is 0.44, and the fibre reinforcement modulus is 4E11, and Poisson ratio is 0.3, uses dichotomy adjustment coupling friction factor (initial range R
0=0.5, R
e=1.0), obtain the operation result in different loads step under arbitrary coupling friction co-efficient value; The distortion parameter tabulation in different loads step is derived, the tabulation of deriving is XYZ Work List at random, be X, the Y row correspond respectively to X, the Y position coordinate value, Z is listed as corresponding to the distortion size, and the deflection of same coordinate position point is subtracted each other, and obtains and the corresponding distortion variable list that uniformly-spaced subtracts each other of step 1).For further obtaining potentiometric map, Work List is converted into the regular grid matrix at random.It is to carry out according to the programming of the Kriging algorithm in the following document that X Y Z Work List at random is converted into grid matrix: Statistics and data analysis in geology, Davis, John C., John Wiley; Sons, Inc.second edition, pp.383.
After being converted into the two-dimensional matrix value, draw and the identical potentiometric map of fringe spacing value size, consistent with analog result until experimental result shown in Fig. 5 (a), binary search finishes, and the friction factor of this moment is 0.88.Fig. 6 (a) has further showed under the correct mode of finite element analysis, along the whole audience displacement potentiometric map of draw direction.Because fiber combine firmly with epoxy matrix, and matrix has very big modulus equally, and the same appearance of fiber is born very big load than strain greatly.Simultaneously, can also obtain more useful indicators from software analysis, as displacement deformation, tension, VonMise stress etc.
Embodiment 2, single bundle carbon fiber strengthen the detection at silicon rubber composite material interface
In the compound substance of selecting for use, matrix material is silicon rubber (GN512, the dimethyl siloxane catalysis addition of ethenyl blocking, Hong Kong advanced technology Industrial Co., Ltd), strengthening body is single diameter 7 μ m, 1000 of every bundles, the carbon fiber of the carbon fibre tow of length 150mm (Dalian Xingke Carbon Fiber Co., Ltd).
This detection method is carried out according to following steps:
1) silicon gel 5g is introduced into 80 * 20 * 2mm
3The tetrafluoro mould in, a branch of carbon fiber inserts wherein 1mm thickness place.The follow procedure sulfuration that slowly heats up in baking oven: 60 ℃, 3 hours, 70 ℃, 3 hours, 80 ℃, 3 hours, 105 ℃, 3 hours, 120 ℃ kept 3 hours, and obtained colourless translucent batten until sulfuration.Sample is taken out from mould, and white paint is sprayed on the surface, is held on the omnipotent puller system (3) after waiting to coat with lacquer drying.With the steady weighing apparatus speed stretching sample of 0.01mm/min, omnipotent puller system record stretching displacement and load, CCD real time record speckle bar graph simultaneously.
Solid line is depicted as the load-displacement curve of the used silicon rubber-carbon fibre composite of this embodiment among Fig. 3, because silicon rubber matter is soft, modulus is little, and not firm with the carbon fiber bonding, reaches 4N and promptly begins slippage, and power loads silicon rubber easy deformation down, overcome friction energy when the distortion aggravation acquires a certain degree, relative displacement will take place and reach a new relaxation balance state in fiber and silicon rubber, continue to stretch, deformed begins a new cycle again, thereby the cycle occurs repeatedly.Laser speckle shown in Fig. 4 (b) uniformly-spaced subtracts each other figure (be 120s interval time) and reflects this phenomenon equally, and the silicon rubber distortion is violent, and the cycle occurs repeatedly.
2) use finite element Ansys10.0 software simulation interface tension process, it is 7E6 that the silicon rubber elastic modulus is set, and Poisson ratio is 0.48, and the fiber reinforcement elastic modulus is 7E11, and Poisson ratio is 0.3, uses dichotomy adjustment coupling friction factor (initial range R
0=0.5, R
e=1.0), obtain the different loads choppy in step, simultaneously shift value is derived, same coordinate point deformation amount is subtracted each other, according to data being converted to the two-dimensional matrix value with embodiment 1 identical method, drawing and the identical potentiometric map of fringe spacing value size, shown in Fig. 5 (b), is 0.65 until the final experimental result friction factor consistent with analog result.Fig. 6 (b) further showed under the correct mode of finite element analysis, along the whole audience displacement potentiometric map of draw direction, because a little less than the silicon rubber matrix combines with fiber, and matrix modulus is little much than fiber, and deformation is more to be born by matrix.Further can obtain more useful indicators from software analysis, as displacement deformation, tension, VonMise stress etc.
The detection of embodiment 3, carbon fibre slice reinforced epoxy composite material interface
Matrix material is that (the bisphenol A type epoxy resin model is E-51 to epoxy resin, hardening agent low-molecular-weight polyamide resin model is 651, Wuxi City A Erzi chemical industry company limited), strengthening body is single diameter 7 μ m, 1000 of every bundles, the carbon fiber of length 150mm (Dalian Xingke Carbon Fiber Co., Ltd).
This detection method is carried out according to following steps:
1) after being mixed by the mass ratio of 10:3, epoxy resin and hardening agent pour 80 * 20 * 2mm into
3The tetrafluoro mould in, complete fiber is smooth with wound form, brush the thin layer epoxy resin prepolymer and be solidified into 100 * 20 * 0.2mm
3Behind the thin slice, 1mm thickness place carefully is placed in one.The follow procedure curing that slowly heats up in baking oven: 60 ℃, 3 hours, 70 ℃, 3 hours, 80 ℃, 3 hours, 105 ℃, 3 hours, 120 ℃ kept 3 hours, and obtained golden yellow translucent batten until curing.Sample is taken out from mould, and white paint is sprayed on the surface, is held on the omnipotent puller system after waiting to coat with lacquer drying.With the steady weighing apparatus speed stretching sample of 0.01mm/min, omnipotent puller system real time record stretching displacement and load, CCD record speckle bar graph simultaneously.
Because epoxy resin matter is hard, modulus is big, and with the carbon cloth bonding closely, load reaches that the curve of load does not still have tangible slippage in the 200N.Laser speckle uniformly-spaced subtracts each other figure (be 120s interval time) and obtains the parallel vertical striped, and the fringe spacing vary stable.
2) use finite element Ansys10.0 software simulation interface tension process, it is 2.8E9 that the epoxy resin elastic modulus is set, and Poisson ratio is 0.44, and the fibre reinforcement modulus is 4E11, and Poisson ratio is 0.3, uses dichotomy adjustment coupling friction factor (initial range R
0=0.5, R
e=1.0), obtain the different loads choppy in step, same coordinate point deformation amount is subtracted each other, according to data being converted to the two-dimensional matrix value with embodiment 1 identical method, drawing and the identical potentiometric map of fringe spacing value size, is 0.85 until the final experimental result friction factor consistent with analog result.The results shows, because fiber combine firmly with epoxy matrix, and matrix has very big modulus equally, and the same appearance of fiber is born very big load than strain greatly.Further, can also obtain more useful indicators from software analysis, as displacement deformation, tension, VonMise stress etc.
Embodiment 4, carbon fibre slice strengthen the detection at silicon rubber composite material interface
In the compound substance of selecting for use, matrix material is silicon rubber (GN512, the dimethyl siloxane catalysis addition of ethenyl blocking, Hong Kong advanced technology Industrial Co., Ltd), strengthening body is single diameter 7 μ m, a branch of 1000, the carbon fiber of the carbon fibre tow of length 150mm (Dalian Xingke Carbon Fiber Co., Ltd).
This detection method is carried out according to following steps:
1) 5g silicon gel is introduced into 80 * 20 * 2mm
3The tetrafluoro mould in, complete fiber is smooth with wound form, brushing the precuring of thin-layer silicon rubber resin becomes 100 * 20 * 0.2mm
3Behind the thin slice, 1mm thickness place carefully is placed in one.The follow procedure sulfuration that slowly heats up in baking oven: 60 ℃, 3 hours, 70 ℃, 3 hours, 80 ℃, 3 hours, 105 ℃, 3 hours, 120 ℃ kept 3 hours, and obtained colourless translucent batten until sulfuration.Sample is taken out from mould, and white paint is sprayed on the surface, is held on the omnipotent puller system (3) after waiting to coat with lacquer drying.With the steady weighing apparatus speed stretching sample of 0.01mm/min, omnipotent puller system record stretching displacement and load, CCD real time record speckle bar graph simultaneously.
Because silicon rubber matter is soft, modulus is little, and not firm with the carbon fiber bonding, reaches 4N and promptly begins slippage.The laser speckle constant duration subtracts each other figure (be 120s interval time) and obtains the parallel vertical striped, and fringe spacing changes violent, overcomes deformation energy when distortion acquires a certain degree, and fiber and silicon rubber generation relative displacement reach a new equilibrium state.Continue to stretch, deformed begins a new cycle, thereby spacing compactness-sparse-compact alternately deformation tendency occurs.
2) use finite element Ansys10.0 software simulation interface tension process, it is 7E6 that the silicon rubber modulus is set, and Poisson ratio is 0.48, and the fibre reinforcement modulus is 4E11, and Poisson ratio is 0.3, uses dichotomy adjustment coupling friction factor (initial range R
0=0.5, R
e=1.0), obtain the different loads choppy in step, simultaneously shift value is derived, same coordinate point deformation amount is subtracted each other, according to data being converted to the two-dimensional matrix value with embodiment 1 identical method, draw and the consistent potentiometric map of fringe spacing value, the friction factor when final experimental result is consistent with analog result is 0.64.The results shows, because a little less than the silicon rubber matrix combines with fiber, and matrix modulus is little much than fiber, and deformation is more to be born by matrix; Further can obtain more useful indicators from software analysis, as displacement deformation, tension, VonMise stress etc.