CN113514178B - Photoelastic stress measurement system and method based on three-wavelength illumination light source - Google Patents
Photoelastic stress measurement system and method based on three-wavelength illumination light source Download PDFInfo
- Publication number
- CN113514178B CN113514178B CN202110410966.2A CN202110410966A CN113514178B CN 113514178 B CN113514178 B CN 113514178B CN 202110410966 A CN202110410966 A CN 202110410966A CN 113514178 B CN113514178 B CN 113514178B
- Authority
- CN
- China
- Prior art keywords
- wavelength
- light source
- arccos
- illumination
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005286 illumination Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005259 measurement Methods 0.000 title claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 13
- 230000003287 optical effect Effects 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 claims description 16
- 230000003247 decreasing effect Effects 0.000 claims description 15
- 238000005516 engineering process Methods 0.000 claims description 8
- 230000009191 jumping Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 238000000691 measurement method Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 12
- 238000004088 simulation Methods 0.000 description 11
- 230000001678 irradiating effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
- G01L1/241—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet by photoelastic stress analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a photoelastic stress measurement system and method based on a three-wavelength illumination light source, and belongs to the technical field of photoelastic measurement. The optical elasticityThe stress measurement system comprises an image acquisition module for respectively acquiring wavelength lambda1Wavelength lambda of2And wavelength lambda3Eliminating the fringe intensity distribution of background light under the illumination of a light source; the main stress difference calculation module is used for acquiring the image stripe information obtained by the image acquisition module and solving the main stress difference of the sample to be detected; and the main stress direction calculation module is used for acquiring the image stripe information obtained by the image acquisition module and solving the main stress direction of the sample to be tested. The invention has convenient use and high measuring speed; the error in the high-level stripe processing is greatly reduced, and the measuring range of photoelastic stress is expanded; the method can process the curve segment containing errors, accurately measure the main stress direction and the main stress difference of the sample, obviously eliminate the errors generated by the dual-wavelength measuring method, and has the advantage of high measuring precision.
Description
Technical Field
The invention relates to a photoelastic stress measurement system and method based on a three-wavelength illumination light source, and belongs to the technical field of photoelastic measurement.
Background
The photoelastic measurement technology is a stress measurement technology based on the birefringence characteristics of a light-transmitting material, and can intuitively reflect the stress distribution condition in the material through a fringe pattern. In addition to the conventional phase shift method, the dual-wavelength method for performing series analysis on photoelastic fringes provides a new idea for quantitative analysis of the magnitude and direction of the internal stress of the material, but still has certain limitations. Errors caused by frequency difference of illumination light are ignored in the calculation process of the dual-wavelength method, when the measured fringe level is low, the error between the result and the actual result is small, but for some samples with large stress variation, the calculation errors are also cumulatively increased along with the increase of the measurement level, so that the measurement accuracy of the samples cannot meet the requirement.
FIG. 1 shows a photoelastic measuring optical path of a two-wavelength method (Chen T Y. digital determination of photoelastic birefringence using two wav)elengths[J]Experimental Mechanics 1997,37(3): 232-. Chen proposes to place 5 x 5 cross-shaped points within one stripe width of the detected stripe pattern for accurate stripe search. The fringe order error of each point under the irradiation of two kinds of wavelength light is calculated respectively, the accurate fringe order of the central point is determined firstly by searching the minimum error value, and then the accurate fringe order of other points is determined by tree search by taking the minimal error value as the seed point. Thus, the value of Nf equals t (σ)1-σ2) The value (sigma) of the difference in principal stresses at each point can be determined1-σ2) Wherein N is the accurate fringe order, f is the fringe value under the corresponding wavelength, and t is the sample thickness. Although the method does not need special treatment of phase unwrapping, the method is only suitable for judging the fringe level in a certain level, and when the fringe density is higher, a larger error is generated (Hanyong, Zhangong Sheng, Ro\2815638]Engineering mechanics, 2008,25(002):62-65.), thus resulting in limited applicability.
Disclosure of Invention
[ problem ] to
The invention aims to solve the problems that: the existing dual-wavelength photoelasticity measuring method is only suitable for judging the fringe order in a certain order, and when the fringe density is higher, a larger error is generated, so that the application range is limited.
[ solution ]
The invention provides a photoelastic stress measurement system and method based on a three-wavelength illumination light source, which are applicable to the technical field of photoelastic measurement.
The photoelastic stress measuring method based on the three-wavelength illumination light source is used for measuring by using a photoelastic stress measuring system based on the three-wavelength illumination light source, and the photoelastic stress measuring system based on the three-wavelength illumination light source comprises the following steps:
the image acquisition module comprises a three-wavelength surface light source, wherein a polarizer, an analyzer, an imaging lens and a photoelectric detector are sequentially arranged on a propagation path of emergent light of the three-wavelength surface light source and are used for respectively acquiring wavelength lambda1Wavelength lambda of2And wavelength lambda3The light source illuminates the sample to be detected to eliminate the stripe intensity distribution of the background light;
the main stress difference calculation module is used for acquiring the image stripe information obtained by the image acquisition module and solving the main stress difference of the sample to be detected;
the main stress direction calculation module is used for acquiring the image stripe information obtained by the image acquisition module and solving the main stress direction of the sample to be tested;
the photoelastic stress measuring method based on the three-wavelength illumination light source comprises the following steps:
the method comprises the following steps: making the polarizer orthogonal to the analyzer at wavelength lambda1Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained1(ii) a At a wavelength λ2Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained2(ii) a At a wavelength λ3Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained3(ii) a The included angle between the polarizer and the analyzer is 45 DEG at wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained4At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained5(ii) a The polarizer and the analyzer are orthogonal and synchronously rotated by 45 DEG at the wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained6At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained7；
Step two: acquiring image stripe information obtained in the first step, and solving a main stress difference sigma (x, y) of a sample to be detected by using Fourier transform, low-pass filtering and phase unwrapping techniques;
step three: and (4) acquiring the image stripe information obtained in the first step, and solving the main stress direction theta of the sample to be tested by using a phase unwrapping technique.
In one embodiment of the present invention, the three-wavelength surface light source is controlled by a switch to emit the wavelength λ at different timings1Wavelength lambda of2And wavelength lambda3Light of three frequencies.
In an embodiment of the present invention, the photodetector is configured to convert optical information transmitted from the imaging lens into an electrical signal, and collect a fringe pattern generated by the three-wavelength surface light source irradiating the sample to be measured.
In one embodiment of the present invention, the wavelengths λ of the three illumination light sources1Wavelength lambda of2Wavelength lambda of3The following relationship is satisfied:
λ1＞λ2＞λ3。
in an embodiment of the present invention, the second step of solving the main stress difference σ (x, y) of the sample to be measured by using fourier transform, low-pass filtering and phase unwrapping techniques specifically includes the following steps:
step1, under an illumination light source with three wavelengths, the normalized fringe intensity distribution after eliminating the influence of background light is respectively as follows:
wherein c is the stress optical coefficient of the material and is a known constant, and t is the thickness of the sample to be measured and is a measurable constant;
step2. for Img respectively1、Img2、Img3Using the formula of multiple angles to obtain I1(x,y)、I2(x,y)、I3(x,y)：
Step3. mixing I2(x, y) are each independently of I1(x, y) and I3(x, y) are compared at the same stripe level to give arccos [ I2(x,y)]In an increasing or decreasing range, arccos [ I ]2(x,y)]The increasing and decreasing turning parts are preliminarily distinguished, and two images { arccos [ I ] which have jumping discontinuities and still contain partial turning parts are obtained respectively2(x,y)]}21And { arccos [ I ]2(x,y)]}23；
Step4, { arccos [ I ] obtained by Step32(x,y)]}21And { arccos [ I ]2(x,y)]}23Calculating the error curve interval, and assigning the value in the error interval as 0 to obtain the function { arccos [ I [ ]2(x,y)]}230；
Step5. pair { arccos [ I ]2(x,y)]}230Fourier transform is carried out to obtain the spectrum distribution of the spectrum, a pass band is selected for low-pass filtering, and the filtering result is converted into space domain information by inverse Fourier transform and is marked as { arccos [ I ]2(x,y)]}230′；
Step6. by judging { arccos [ I ]2(x,y)]}230The positive and negative values of the segment value of the median error curve construct a reversal judgment matrix, which is marked as L (x, y), and the L (x, y) and { arccos [ I ]2(x,y)]}23Multiplication to make { arccos [ I ]2(x,y)]}23Corrected error of (2) to obtain an image { arccos [ I ] having only jumping discontinuities2(x,y)]}23′；
Step7. pair { arccos]I2(x,y)]}23' phase unwrapping processing is performed by adopting a phase unwrapping algorithm, and continuous stripe distribution can be spliced and recorded as { arccos [ I2(x,y)]}23FI.e., [ arccos (1-2 Img)2)]23F。
In an embodiment of the present invention, the principal stress difference σ (x, y) of the sample to be tested satisfies the following relation:
in an embodiment of the present invention, the step three of solving the principal stress direction θ of the sample to be tested by using the phase unwrapping technique includes the following steps:
wherein, I20、I30Respectively represent the wavelength lambda2And wavelength lambda3The intensity of the illumination when illuminated.
In an embodiment of the present invention, the principal stress direction θ of the sample to be measured satisfies the following relation:
and (5) obtaining the actual stress direction of the sample to be detected by adopting a phase unwrapping algorithm on the theta.
Advantageous effects
(1) According to the photoelastic stress measurement system based on the three-wavelength illumination light source, the illumination light sources with three frequencies are used for respectively irradiating a sample to be measured to obtain a fringe pattern containing sample stress information, so that the stress magnitude and direction can be measured, and the photoelastic stress measurement system based on the three-wavelength illumination light source has the advantages of convenience in use and high measurement speed.
(2) The photoelastic stress measuring system based on the three-wavelength illumination light source disclosed by the invention can measure a large-caliber sample under the condition of not obviously increasing the weight of an instrument, greatly reduces the error in advanced stripe processing, and expands the measurement range of photoelastic stress.
(3) When the dual-wavelength method is used for solving the main stress difference, an error interval is generated along with the increase of the fringe order, and the error interval is increased along with the increase of the fringe order (see c in fig. 5), so that the principle of measuring the stress by the dual-wavelength method cannot be corrected. However, the photoelastic stress measuring method based on the three-wavelength illumination light source disclosed by the invention is realized by additionally introducing a wavelength lambda3The illumination light source can process the curve segment containing errors by applying the modern digital image processing technology, accurately measure the main stress direction and the main stress difference of the sample, obviously eliminate the errors generated by the dual-wavelength measuring method (see c in figure 6), and has the advantage of high measuring precision.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic structural diagram of a dual-wavelength stress measuring instrument in the prior art of photoelasticity measurement;
FIG. 2 is a schematic structural diagram of a three-wavelength illumination-based photoelastic stress measurement system according to the present invention;
FIG. 3 is a schematic structural diagram of an image acquisition module according to the present invention;
FIG. 4 is a flowchart of a method for measuring photoelastic stress based on a three-wavelength illumination source according to the present invention;
FIG. 5 is a diagram illustrating a simulation result of stress calculation by a dual-wavelength method in the prior art photoelastic measurement;
FIG. 6 is a schematic diagram of a simulation result of stress calculation using a three-wavelength illumination light source-based photoelastic stress measurement system and method of the present invention;
FIG. 7 shows a schematic view of the present invention2(x, y) are each independently of I1(x, y) and I3(x, y) comparison at the same stripe level, { arccos [ I ] with jumping discontinuities and still containing partial breakover2(x,y)]}21And { arccos [ I ]2(x,y)]}23A simulation result schematic diagram;
FIG. 8 is a diagram illustrating simulation results of screening error curve segment positions according to the present invention;
FIG. 9 is a diagram illustrating simulation results of processing the error curve segment filtered in FIG. 8 using Fourier transform, inverse Fourier transform and low pass filtering techniques in accordance with the present invention;
FIG. 10 is a diagram illustrating simulation results of the symbol decision matrix L (x, y) according to the present invention;
FIG. 11 shows { arccos [ I ] according to the present invention2(x,y)]}23The increasing and decreasing turning parts are further differentiated to obtain { arccos [ I ] with jumping discontinuities only2(x,y)]}23' simulation result schematic diagram;
FIG. 12 is a graph of the present invention vs { arccos [ I ] in FIG. 112(x,y)]}23' A schematic diagram of a simulation result of true phase distribution obtained by unwrapping processing is made.
In the figure, 1, a dual-wavelength laser light source; 2. a polarizer; 3. a first quarter wave plate; 4. a collimating lens; 5. a sample to be tested; 6. an imaging lens; 7. a second quarter wave plate; 8. an analyzer; 9. a photodetector; 10. a three-wavelength surface light source; 11. an image acquisition module; 12. a main stress difference calculation module; 13. and a principal stress direction calculation module.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Example 1:
as shown in fig. 2, the present embodiment provides a photoelastic stress measurement system based on a three-wavelength illumination light source, which includes an image acquisition module 11, configured to acquire a wavelength λ when a relative direction of deflection axes of a polarizer 2 and an analyzer 8 is changed1、λ2And λ3Eliminating the fringe intensity distribution of background light under the illumination of a light source; the main stress difference calculating module 12 is used for solving the main stress difference of the image stripe information obtained by the image obtaining module 11 by using modern digital image processing technologies such as Fourier transform, low-pass filtering, phase unwrapping and the like; and the main stress direction calculation module 13 is configured to collect the image stripe information obtained by the image obtaining module 11, and solve the main stress direction by using a phase unwrapping technique.
As shown in fig. 3, the image acquisition module 11 includes a three-wavelength surface light source 10, and the three-wavelength surface light source 10 is controlled by a switch to emit a wavelength λ at different times1Wavelength lambda of2And wavelength lambda3Light of three frequencies; a polarizer 2, an analyzer 8, an imaging lens 6 and a photoelectric detector 9 are sequentially arranged on a propagation path of the light emitted by the three-wavelength surface light source 10; a sample 5 to be tested is placed between the polarizer 2 and the analyzer 8; the photoelectric detector 9 converts the light information transmitted from the imaging lens 6 into an electric signal, and collects a fringe pattern generated by irradiating the sample 5 to be measured by the three-wavelength surface light source 10.
Example 2:
as shown in fig. 4, the present embodiment discloses a method for measuring photoelastic stress based on a three-wavelength illumination light source, where the method for measuring photoelastic stress based on a three-wavelength illumination light source uses the system for measuring photoelastic stress based on a three-wavelength illumination light source disclosed in embodiment 1 to perform measurement, and includes the following steps:
step (ii) ofFirstly, the method comprises the following steps: making the polarizer 2 orthogonal to the analyzer 8 at the wavelength lambda1Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained1(ii) a At a wavelength λ2Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained2(ii) a At a wavelength λ3Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained3(ii) a The included angle between the polarizer 2 and the analyzer 8 is 45 DEG at the wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained4At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained5(ii) a The polarizer 2 and the analyzer 8 are orthogonal and synchronously rotated by 45 DEG at the wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained6At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained7；
Step two: acquiring image stripe information obtained in the first step, and solving a main stress difference sigma (x, y) of the sample 5 to be detected by using Fourier transform, low-pass filtering and phase unwrapping techniques;
step three: and (4) acquiring the image stripe information obtained in the first step, and solving the main stress direction theta of the sample 5 to be tested by using a phase unwrapping technique.
In the second step, the main stress difference sigma (x, y) of the sample 5 to be measured is solved by using Fourier transform, low-pass filtering and phase unwrapping technologies, and the method specifically comprises the following steps:
step1, under an illumination light source with three wavelengths, the normalized fringe intensity distribution after eliminating the influence of background light is respectively as follows:
wherein c is the stress optical coefficient of the material and is a known constant, and t is the thickness of the sample 5 to be measured and is a measurable constant;
step2. for Img respectively1、Img2、Img3Using the formula of multiple angles to obtain I1(x,y)、I2(x,y)、I3(x,y)：
The normalized stripe intensity distribution Img after the background light influence is eliminated, which is obtained in the step one1The square of the sine function with the principal stress difference σ (x, y) as the argument cannot be applied directly to Img1The real main stress difference distribution is obtained through root-opening processing, so a multiple angle formula is used to obtain a cosine function with the main stress difference sigma (x, y) as an independent variable: i is1(x, y) in this case, the two-wavelength method by introducing Img2To obtain I2(x, y) to determine I1(x, y) and I2(x, y) relative size at the same stripe order, to obtain arccos [ I2(x,y)]The main stress difference is obtained by phase unwrapping technique, but as the fringe order increases, I1(x, y) and I2The relative size of (x, y) does not sufficiently guarantee arccos [ I ]2(x,y)]Is monotonically increasing or monotonically decreasing, as shown by a in FIG. 5, at [280, 310]Within the interval, I2(x,y)>I1(x, y) but arccos [ I2(x,y)]Is not monotonically decreasing, inThe area marked in the rectangular box is monotonically increased (see b in fig. 5), so that the area which does not need to be inverted originally is also inverted (see c in fig. 5), and a curve segment which is concave downwards, namely an error interval is generated, the error interval is increased along with the increase of the fringe order, and the principle of measuring stress by the dual-wavelength method cannot be corrected. In order to solve the error and accurately obtain the main stress difference, the invention applies two wavelengths which are lambda respectively in a dual-wavelength method1And λ2On the basis of the illumination light source, an additional introduction of a wavelength of lambda3And let the wavelengths λ of the three illumination light sources1、λ2、λ3Successively smaller, i.e. λ1>λ2>λ3To 1, pair1(x,y)、I2(x,y)、I3(x, y) processing. The simulation result of stress calculation by the photoelastic stress measurement method based on the three-wavelength illumination light source provided by the invention is shown in fig. 6, wherein a in fig. 6 is I1(x,y)、I2(x,y)、I3(x, y) simulation diagram, b in FIG. 6 is arccos [ I ]1(x,y)]、arccos[I2(x,y)]、arccos[I3(x,y)]C in fig. 6 is the phase distribution simulation result of the unwrapped phase distribution which can reflect the main stress difference distribution, and it can be found that c in fig. 6 does not have the concave curve segment, i.e. the error interval, similar to c in fig. 5. By adopting the photoelastic stress measurement method based on the three-wavelength illumination light source, inevitable errors in a dual-wavelength method are obviously eliminated. Further, continue to I1(x,y)、I2(x,y)、I3(x, y) processing;
step3. mixing I2(x, y) are each independently of I1(x, y) and I3(x, y) the comparisons are made at the same stripe order:
for I2(x, y) and I1(x, y), the same stripe order needs to be satisfied
Namely, it is
If I1(x,y)>I2(x, y) then I2(x, y) can be regarded as [0, π/2 ]]Internal decreasing function, its inverse arccos [ I2(x,y)]Can be regarded as [0, pi/2 ]]Is increased function of
If I1(x,y)<I2(x, y) then I2(x, y) can be regarded as [ pi/2, pi]Increasing function of inner, its inverse function arccos [ I2(x,y)]Can be regarded as [ pi/2, pi]Is a decreasing function of
After recording treatment
The fringe image is { arccos [ I ]2(x,y)]}21As shown in fig. 7 a;
for I2(x, y) and I3(x, y), the same stripe order needs to be satisfied
Namely, it is
If I2(x,y)>I3(x, y) then I2(x, y) can be regarded as [0, π/2 ]]Internal decreasing function, its inverse arccos [ I2(x,y)]Can be regarded as [0, pi/2 ]]An increasing function of the inside, an
If I2(x,y)<I3(x, y) then I2(x, y) can be regarded as [ pi/2, pi]Increasing function of inner, its inverse function arccos [ I2(x,y)]Can be regarded as [ pi/2, pi]Within a decreasing function, an
After recording treatment
The fringe image is { arccos [ I ]2(x,y)]}23Such asB in fig. 7;
thus by separately judging I2(x, y) and I1(x, y) and I2(x, y) and I3(x, y) the relationship of the sizes of the same stripe order, arccos [ I ] can be obtained2(x,y)]Coarse range of increase and decrease, thereby decreasing arccos [ I ]2(x,y)]The increasing and decreasing turning parts are preliminarily distinguished to respectively obtain two images { arccos [ I ] with jumping discontinuities and still containing partial turning parts2(x,y)]}21And { arccos [ I ]2(x,y)]}23(see FIG. 7);
step4, { arccos [ I ] obtained by Step32(x,y)]}21And { arccos [ I ]2(x,y)]}23Calculating the error curve interval, and assigning the value in the error interval as 0 to obtain the function { arccos [ I [ ]2(x,y)]}230；
Due to three wavelengths lambda1、λ2、λ3With a certain difference, I2(x, y) are each independently of I1(x,y)、I3(x, y) obtained by comparison with respect to I2The increasing (x, y) part of the curve contains decreasing tiny curve segments, and with respect to I2The part of the curve with decreasing (x, y) comprises the tiny curve segments with increasing (x, y) and simultaneously satisfies I corresponding to the same stripe level1(x,y)<I2(x, y) and I2(x,y)>I3The curve segment of (x, y) and the same stripe order satisfy I simultaneously1(x,y)>I2(x, y) and I2(x,y)<I3(x, y), which is unavoidable, and these errors are larger as the fringe order increases;
let { arccos [ I ]2(x,y)]}21The error curve segment of the ith-level stripe is [ Ei]21，{arccos[I2(x,y)]}23The error curve segment of the ith-level stripe is [ Ei]23I is 1,2,3 …, by
Can determine the position of these error curve segments: { arccos [ I2(x,y)]}21And { arccos [ I ]2(x,y)]}23In satisfy
Is error-free and satisfies
Is erroneous, so that [ Ei ] can be selected]21And [ Ei]23I is 1,2,3 …, where only [ Ei ] needs to be screened out]21Or [ Ei]23. To screen [ Ei]23For example, the value in the error interval is set to 0, and the obtained fringe image is recorded as { arccos [ I [ ]2(x,y)]}230The result is shown in FIG. 8, where the positions of a series of line segments parallel to the x-axis represent the error curve segment [ Ei ]]23I is the position of 1,2,3 …;
step5. pair { arccos [ I ]2(x,y)]}230Fourier transform is carried out to obtain the spectrum distribution of the spectrum, a proper passband is selected for low-pass filtering, and the filtering result is converted into space domain information by inverse Fourier transform and is marked as { arccos [ I ]2(x,y)]}230′；
For { arccos [ I ]2(x,y)]}230Fourier transform is carried out to obtain a in figure 9, a proper passband is selected according to a measuring object, low-pass filtering is carried out on a in figure 9, and inverse Fourier transform is carried out on the filtering result to obtain a space domain image which is marked as { arccos [ I ]2(x,y)]}230' as shown in b of fig. 9, wherein all the descending curve segments are processed error curve segments, and the straight line parallel to the x-axis is the reference line in the figure. The intersection point of the ith descending curve segment b in FIG. 9 and the reference line is the error interval [ Ei]23(b in fig. 9 marks the middle of the partial error interval with an arrow). From b in FIG. 9, error interval [ Ei ]]23The curve value of the part on the left side of the middle point is larger than zero, the part needs to be turned, the curve value of the part on the right side of the middle point is smaller than zero, and the part does not need to be turned;
step6. by judging { arccos [ I ]2(x,y)]}230The positive and negative values of the segment value of the middle error curve construct a turnover judgment matrix, which is marked as L (x, y),for judging { arccos [ I ]2(x,y)]}23Whether the curve in each interval needs to be inverted.
B error interval [ Ei ] in FIG. 9]23The curve value of the left part of the midpoint is larger than zero, the part originally needs to be overturned, but cannot be overturned when the dual-wavelength method is adopted for judgment, so that overturning correction is needed, and an error interval [ Ei ]]23The value of the curve in the portion to the right of the midpoint is less than zero and the portion does not need to be flipped. If { arccos [ I ]2(x,y)]}230If the value of the middle error curve segment is positive, let the value of the corresponding position in the L (x, y) matrix be-1, and the remaining position values be 1, then L (x, y) is obtained as shown in fig. 10. Mixing L (x, y) with { arccos [ I ]2(x,y)]}23Multiplying by the error interval [ Ei]23The left part of the midpoint is flipped and the right part remains unchanged, so that { arccos [ I2(x,y)]}23Corrected error of (2) to obtain an image { arccos [ I ] having only jumping discontinuities2(x,y)]}23', as shown in FIG. 11;
step7. pair { arccos [ I ]2(x,y)]}23' phase unwrapping processing is performed by adopting a phase unwrapping algorithm, and continuous stripe distribution can be spliced and recorded as { arccos [ I2(x,y)]}23FI.e., [ arccos (1-2 Img)2)]23F(see fig. 12), the magnitude of the primary stress difference σ (x, y) of the sample 5 to be measured satisfies the following relation:
thirdly, solving the main stress direction theta of the sample 5 to be tested by using a phase unwrapping technology, wherein the method comprises the following specific steps:
wherein, I20、I30Respectively represent the wavelength lambda2And wavelength lambda3The intensity of the illumination when illuminated. Therefore, it is not only easy to use
Therefore, it is
And (5) processing the theta by adopting a phase unwrapping algorithm to obtain the actual stress direction of the sample to be measured.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (4)
1. A photoelastic stress measurement method based on a three-wavelength illumination light source is characterized in that a photoelastic stress measurement system based on a three-wavelength illumination light source is applied for measurement, and the photoelastic stress measurement system based on the three-wavelength illumination light source comprises the following steps:
the image acquisition module (11) comprises a three-wavelength surface light source (10), wherein a polarizer (2), an analyzer (8), an imaging lens (6) and a photoelectric detector (9) are sequentially arranged on a propagation path of emergent light of the three-wavelength surface light source (10) and used for respectively acquiring wavelength lambda1Wavelength lambda of2And wavelength lambda3The light source illuminates the sample (5) to be detected to eliminate the stripe intensity distribution of the background light;
the main stress difference calculation module (12) is used for acquiring the image stripe intensity distribution obtained by the image acquisition module (11) and solving the main stress difference of the sample (5) to be detected;
the main stress direction calculation module (13) is used for acquiring the image stripe intensity distribution obtained by the image acquisition module (11) and solving the main stress direction of the sample (5) to be detected;
the photoelastic stress measuring method based on the three-wavelength illumination light source comprises the following steps:
the method comprises the following steps: the polarizer (2) is orthogonal to the analyzer (8) at a wavelength lambda1Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained1(ii) a At a wavelength λ2Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained2(ii) a At a wavelength λ3Under the illumination of the light source, the normalized fringe intensity distribution Img after eliminating the influence of background light is obtained3(ii) a The included angle between the polarizer (2) and the analyzer (8) is 45 DEG at the wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained4At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained5(ii) a The polarizer (2) and the analyzer (8) are orthogonal and synchronously rotated for 45 DEG at the wavelength lambda2Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained6At a wavelength λ3Under the illumination of the light source, the fringe intensity distribution Img after eliminating the influence of the background light is obtained7；
Step two: acquiring the intensity distribution of the image stripes obtained in the first step, and solving the main stress difference sigma (x, y) of the sample (5) to be detected by using Fourier transform, low-pass filtering and phase unwrapping techniques;
step three: acquiring the intensity distribution of the image stripes obtained in the first step, and solving the main stress direction theta of the sample (5) to be detected by using a phase unwrapping technique;
solving the main stress difference sigma (x, y) of the sample (5) to be detected by applying Fourier transform, low-pass filtering and phase unwrapping technology in the second step specifically comprises the following steps:
step1, under an illumination light source with three wavelengths, the normalized fringe intensity distribution after eliminating the influence of background light is respectively as follows:
wherein c is the stress optical coefficient of the material and is a known constant, and t is the thickness of the sample (5) to be measured and is a measurable constant;
step2. for Img respectively1、Img2、Img3Using the formula of multiple angles to obtain I1(x,y)、I2(x,y)、I3(x,y)：
Step3. mixing I2(x, y) are each independently of I1(x, y) and I3(x, y) are compared at the same stripe level to give arccos [ I2(x,y)]In an increasing or decreasing range, arccos [ I ]2(x,y)]The increasing and decreasing turning parts are preliminarily distinguished, and two images { arccos [ I ] which have jumping discontinuities and still contain partial turning parts are obtained respectively2(x,y)]}21And { arccos [ I ]2(x,y)]}23；
Step4, { arccos [ I ] obtained by Step32(x,y)]}21And { arccos [ I ]2(x,y)]}23Calculating the error curve interval, and assigning the value in the error interval as 0 to obtain the function { arccos [ I [ ]2(x,y)]}230；
Step5. pair { arccos [ I ]2(x,y)]}230Fourier transform is carried out to obtain the spectrum distribution of the spectrum, a pass band is selected for low-pass filtering, and the filtering result is converted into space domain information by inverse Fourier transform and is marked as { arccos [ I ]2(x,y)]}230′；
Step6. by judging { arccos [ I ]2(x,y)]}230The positive and negative values of the segment value of the median error curve construct a reversal judgment matrix, which is marked as L (x, y), and the L (x, y) and { arccos [ I ]2(x,y)]}23Multiplication to make { arccos [ I ]2(x,y)]}23Corrected error of (2) to obtain an image { arccos [ I ] having only jumping discontinuities2(x,y)]}23′；
Step7. pair { arccos [ I ]2(x,y)]}23' phase unwrapping processing is performed by adopting a phase unwrapping algorithm, and stripe intensity distribution can be spliced and recorded as { arccos [ I ]2(x,y)]}23FI.e., [ arccos (1-2 Img)2)]23F；
The main stress difference sigma (x, y) of the sample (5) to be tested satisfies the following relational expression:
the third step of solving the main stress direction theta of the sample (5) to be tested by applying the phase unwrapping technology comprises the following processes:
wherein, I20、I30Respectively represent the wavelength lambda2And wavelength lambda3The intensity of the illumination light during illumination;
the main stress direction theta of the sample (5) to be tested satisfies the following relational expression:
and (5) obtaining the main stress direction of the sample (5) to be detected by adopting a phase unwrapping algorithm on the theta.
2. The photoelastic stress measurement method based on three-wavelength illumination source of claim 1, wherein the photoelastic stress measurement system based on three-wavelength illumination source is characterized in that the three-wavelength surface light source (10) is controlled by a switch to emit the wavelength λ at different times1Wavelength lambda of2And wavelength lambda3Light of three frequencies.
3. The photoelastic stress measurement method based on three-wavelength illumination light source according to claim 1 or 2, characterized in that in the photoelastic stress measurement system based on three-wavelength illumination light source, the photodetector (9) is used to convert the light information transmitted from the imaging lens (6) into an electrical signal and collect the fringe intensity distribution generated by the illumination of the sample (5) to be measured by the three-wavelength surface light source (10).
4. The photoelastic stress measurement method based on three-wavelength illumination light source according to claim 2, wherein the wavelength λ of the three-wavelength surface light source (10)1Wavelength lambda of2Wavelength lambda of3The following relationship is satisfied:
λ1>λ2>λ3。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110410966.2A CN113514178B (en) | 2021-04-16 | 2021-04-16 | Photoelastic stress measurement system and method based on three-wavelength illumination light source |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110410966.2A CN113514178B (en) | 2021-04-16 | 2021-04-16 | Photoelastic stress measurement system and method based on three-wavelength illumination light source |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113514178A CN113514178A (en) | 2021-10-19 |
| CN113514178B true CN113514178B (en) | 2022-04-15 |
Family
ID=78062547
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110410966.2A Active CN113514178B (en) | 2021-04-16 | 2021-04-16 | Photoelastic stress measurement system and method based on three-wavelength illumination light source |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113514178B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114061803B (en) * | 2021-11-29 | 2022-09-27 | 江南大学 | Circular polarization photoelastic stress measuring system and measuring method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103499405A (en) * | 2013-10-15 | 2014-01-08 | 上海理工大学 | Device and method for detecting residual stress of transparent plastic product |
| CN103759831A (en) * | 2014-01-03 | 2014-04-30 | 南京邮电大学 | Spectral measurement device and spectral measurement method based on elasto-optical effect |
| US9372127B1 (en) * | 2015-05-04 | 2016-06-21 | Eastman Kodak Company | Printing system controlled using photoelastic measurement device |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN85100284B (en) * | 1985-04-01 | 1985-09-10 | 北京大学 | Automatic measuring method and apparatus |
| JP3746287B2 (en) * | 2004-01-15 | 2006-02-15 | 学校法人東京電機大学 | Stress measuring method and apparatus |
| FR2867561B1 (en) * | 2004-03-11 | 2007-02-02 | Commissariat Energie Atomique | DISTRIBUTED MEASUREMENT SYSTEM OF THE CURVES OF A STRUCTURE |
| CN100590425C (en) * | 2006-05-30 | 2010-02-17 | 长沙科创计算机***集成有限公司 | Method for detecting stress of transparency flat plate or model material |
| JP2010117326A (en) * | 2008-11-14 | 2010-05-27 | Harison Toshiba Lighting Corp | Device and method for measuring stress |
| US9134842B2 (en) * | 2012-10-04 | 2015-09-15 | Corning Incorporated | Pressure sensing touch systems and methods |
| CN105675186B (en) * | 2016-01-26 | 2018-09-18 | 中国科学院声学研究所 | Method for measuring stress based on dynamic photoelasticity sexual system |
| CN106273094A (en) * | 2016-10-17 | 2017-01-04 | 中国航空动力机械研究所 | Epoxy resin pattern manufacture method for photoelastic test |
| CN108562385A (en) * | 2018-03-08 | 2018-09-21 | 北京理工大学 | A kind of microoperation Real-time force feedback device and method based on photoelastic effect |
| CN108760112B (en) * | 2018-05-25 | 2020-06-30 | 中国科学院上海光学精密机械研究所 | Stress measuring device and method based on diffraction overlapping iterative algorithm |
| CN109000829B (en) * | 2018-10-09 | 2020-03-24 | 江南大学 | Capacitive pressure distribution monitoring material and method |
| CN111879708B (en) * | 2019-06-03 | 2021-07-27 | 江南大学 | Device and method for measuring laser damage three-dimensional structure |
-
2021
- 2021-04-16 CN CN202110410966.2A patent/CN113514178B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103499405A (en) * | 2013-10-15 | 2014-01-08 | 上海理工大学 | Device and method for detecting residual stress of transparent plastic product |
| CN103759831A (en) * | 2014-01-03 | 2014-04-30 | 南京邮电大学 | Spectral measurement device and spectral measurement method based on elasto-optical effect |
| US9372127B1 (en) * | 2015-05-04 | 2016-06-21 | Eastman Kodak Company | Printing system controlled using photoelastic measurement device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113514178A (en) | 2021-10-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8654350B2 (en) | Inspecting method and inspecting apparatus for substrate surface | |
| US7495762B2 (en) | High-density channels detecting device | |
| US7566517B1 (en) | Feature printability optimization by optical tool | |
| WO2009093453A1 (en) | Analysis device and analysis method | |
| KR940016660A (en) | Thin Film Thickness Measurement Apparatus and Method | |
| CN109470154B (en) | Film thickness initial value measuring method suitable for spectrum ellipsometer | |
| EP0796424A1 (en) | A method of obtaining information | |
| CN113514178B (en) | Photoelastic stress measurement system and method based on three-wavelength illumination light source | |
| CN113538381A (en) | Method and system for rapidly detecting Mueller matrix of sample in weak light field | |
| CN115356270A (en) | Method for noise suppression and decomposition in mixture Mueller offset measurement | |
| CN113008426B (en) | Double-frequency quantitative photoelastic measuring instrument and measuring method | |
| CN113447126B (en) | Spatial modulation polarization detection system and design method | |
| JP7031664B2 (en) | Object detection device, object detection method, and program | |
| CN117128877A (en) | Film thickness detection method, computer and system | |
| CN110118754B (en) | Rapid measurement method for ultrathin film optical constant | |
| CN114061803B (en) | Circular polarization photoelastic stress measuring system and measuring method | |
| CN109115695B (en) | Method for extracting optical constants and Euler angles of anisotropic body materials | |
| JPH0781840B2 (en) | Optical film thickness measuring device | |
| CN113008427B (en) | Liquid crystal photoelastic stress measuring instrument and measuring method | |
| CN107643135B (en) | Stress analysis method based on time domain phase unwrapping | |
| US20230058064A1 (en) | Film thickness measuring device and film thickness measuring method | |
| JPH07294220A (en) | Method and apparatus for detecting film thickness of multilayered thin film | |
| KR102176199B1 (en) | Ellipsometer | |
| JP2000111472A (en) | Device for measuring birefringence and device for measuring orientation of film | |
| US20120050710A1 (en) | Method of measuring mark position and measuring apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |