OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 6 — Mar. 17, 2008
  • pp: 3955–3960
« Show journal navigation

Stress distribution and induced birefringence analysis for pressure vector sensing based on single mode fibers

Zhengyong Li, Chongqing Wu, Hui Dong, P. Shum, C. Y. Tian, and S. Zhao  »View Author Affiliations


Optics Express, Vol. 16, Issue 6, pp. 3955-3960 (2008)
http://dx.doi.org/10.1364/OE.16.003955


View Full Text Article

Acrobat PDF (638 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We propose a theoretical approach to analyze the pressure stress distribution in single mode fibers (SMFs) and achieve the analytical expression of stress function, from which we obtain the stress components with their patterns in the core and compute their induced birefringence. Then we perform a pressure vector sensing based on ~2 km SMF. Using Mueller matrix method we measure the birefringence vectors which are employed to compute the pressure magnitudes and their orientation. When rotating the pressure around the fiber, the corresponding birefringence vector rotates around a circle with double speed. Statistics show the average deviation of calculated pressure-magnitude to practical value is ~0.17 N and it is ~0.85° for orientation.

© 2008 Optical Society of America

Introduction

Polarization dependent effects are becoming a major topic for modern optical communication and fiber sensing [1

1. D. Chowdhury and D. Wilcox, “Comparison between optical fiber birefringence induced by stress anisotropy and geometric deformation,” IEEE J. Sel. Topics Quantum Electron. 6, 227–232 (2000). [CrossRef]

]. The fiber birefringence is one of the basic polarization effects which contribute to the well-known PMD, meanwhile it is an important physical parameter for many fiber-based sensing applications such as pressure sensors, temperature sensors, and etc.. Several methods have been developed to analyze the stress-induced birefringence especially the intrinsic birefringence in single mode fibers (SMFs), for example, a vector perturbation model proposed in [1

1. D. Chowdhury and D. Wilcox, “Comparison between optical fiber birefringence induced by stress anisotropy and geometric deformation,” IEEE J. Sel. Topics Quantum Electron. 6, 227–232 (2000). [CrossRef]

], a modified coupled-mode theory described in [2

2. J. Sakai and T. Kimura, “Birefringence and polarization characteristics of single mode optical fibers under electric deformations,” IEEE J. Quantum Electron. , 17, 1041–1051 (1981). [CrossRef]

], an optical tomographic technique applied in [3

3. Y. Park, U. C. Paek, and D.Y. Kim, “Determination of stress-induced intrinsic birefringence in a single-mode fiber by measurement of the two-dimensional stress profile,” Opt. Lett. , 27, 1291–1293 (2002). [CrossRef]

], an analogy between the pressure-induced strain and the thermal strain employed in [4

4. K.S. Chang, “Pressure-induced birefringence in a coated highly birefringent optical fiber,” J. Lightwave Technol. , 8, 1850–1855 (1990). [CrossRef]

], and etc.. Among them the stress distribution is the key issue in calculating the birefringence. In Ref. 3

3. Y. Park, U. C. Paek, and D.Y. Kim, “Determination of stress-induced intrinsic birefringence in a single-mode fiber by measurement of the two-dimensional stress profile,” Opt. Lett. , 27, 1291–1293 (2002). [CrossRef]

and Ref. 5

5. S. M. Pietralunga, M. Ferrario, M. Tacca, and M. Martinelli, “Local Birefringence in Unidirectionally Spun Fibers,” J. Lightwave Technol. , 24, 4030–4038 (2006). [CrossRef]

the 2-D stress profile was measured with an optical tomographic technique while Ref. 6

6. K. Saitoh, M. Koshiba, and Y. Tsuji, “Stress analysis method for elastically anisotropic material based optical waveguides and its application to strain-induced optical waveguides,” J. Lightwave Technol. , 17, 255–259 (1999). [CrossRef]

introduced the finite element method for stress analysis, which are general methods for analysis of stress distribution in optical waveguides. To investigate the stress distribution of a round fiber under external forces, we think there should be a simple method when taking account of the circular symmetry of the fiber. In this paper, we propose a convenient approach to analyze the stress distribution in SMFs, by which we easily obtain all the stress components and the stress patterns in the core. Then we compute the stress-induced birefringence and perform a pressure vector sensing based on ~2 km SMF with the validity analyzed.

1. Pressure stress distribution in SMFs

Consider a round SMF with the core and cladding under the lateral pressure P as shown in Fig. 1, whose right part is the cross section of the fiber in the coordinate system, where a and R are the radii of the core and cladding. Since the residual stresses are relatively very small to the pressure stress in the SMF, we omit their effects on stress distribution.

Fig. 1. (a) The schematic graph and (b) enlarged cross section of SMF under lateral pressure.

In the case that pressed length (5 mm in our experiment) is much larger than the size (0.125 mm) of fiber cross-section, the stress distribution in pressed fiber can be treated as a plane stress problem which means it only relates to the cross section while independent of the fiber length. Thus the stress function on A(r, θ) and the stress tensor components in the fiber are described by [7

7. S. Timoshenko and J. N. Goodier, Theory of elasticity (McGraw-Hill, 1970), Chap. 4.

]

(2r2+1rr+1r22θ2)2A(r,θ)=0,
(1)
σr=1rAr+1r22Aθ2,σθ=2Ar2,τrθ=r(1rAθ),σz=τrz=τzθ=0.
(2)

The integral formula solution of biharmonic equation (1) for a circle plane with radius R is [8

8. D.H. Yu, Mathematical theory of natural boundary element method (Science Press, Beijing, 1993), Chap. 3.

]

A(r,θ)=02π[Rrcos(θφ)2πRXμ(φ)R2r24πRv(φ)]Xdφ,
X=(R2r2)[R2+r22Rrcos(θφ)],
(3)

in which μ(φ)=AQ = ΣM⃗i, v(φ)=(An)Q=Fi . Here Q is an arbitrary point on the circle boundary, n is the normal of the integration boundary, φ is the integration variable, while F⃗i and M⃗i are respectively a force element existing in the integration interval ( BQ in Fig. 1b) and its moment with respect to point Q [7

7. S. Timoshenko and J. N. Goodier, Theory of elasticity (McGraw-Hill, 1970), Chap. 4.

]. By simple calculation we get

μ(φ)={PRsinφ,0φπ0πφ2π,v(φ)={Psinφ,0φπ0πφ2π.
(4)

Inserting expressions (4) into Eq. (3), and compute the integration we obtain

A(r,θ)=P4πR[2πRrsinθ+4Rrsinθtan1(2RrsinθR2r2)+2(R2r2)],(0r<R).
(5)

Since most optical field just exists in the core, we focus on its stress distribution which mainly determines the birefringence of the fiber under press. Considering = r/R≪1(ra) we take it into account till the second power ( 2) and neglect the higher order terms, then substituting Eq. (5) into Eq. (2) we get the non-zero stress components

σr=σ(8r˜2cos2θ3),σθ=σ[4(cos2θr˜2)3],τ=2σ(1r˜2)sin(2θ),
σ=P(πR)1[14r˜2cos(2θ)],(r~=rR).
(6)

Transforming these components to the rectangular coordinates shown in Fig. 1(b) we have

σx=σ{3+2r˜2[1+3cos(2θ)]},σy=σ[1+4r˜2cosθcos(3θ)],
τxy=4σr˜2sin(2θ)[(r˜21)cos(2θ)+r˜2].
(7)

Given parameters 2R = 125 µm, a = 4.25 µm, P = 200 N/m which means a force of 1 N acts on a 5-mm-length bare fiber, by Eq. (7) we obtain the stress patterns in the fiber core which are plotted in Fig. 2. Comparing the three parts we find the normal stresses σx and σy are nearly 100 times more than shearing stress τxy and the latter just fluctuates around zero, so the normal stresses bring main contributions to the birefringence in SMFs. In addition, the magnitude of σx is about three times as high as σy, therefore the normal stress along the direction of pressure P is the greatest stress component in the fiber.

Fig. 2. Stress distributions (a): σx, (b): σy, and (c): τxy in the core of the fiber under lateral pressure.

2. Stress-induced birefringence vector with relation to lateral pressure

Based on above stress tensor components in Eq. (6), we can get the stress-induced perturbation in relative permittivity tensor as [3

3. Y. Park, U. C. Paek, and D.Y. Kim, “Determination of stress-induced intrinsic birefringence in a single-mode fiber by measurement of the two-dimensional stress profile,” Opt. Lett. , 27, 1291–1293 (2002). [CrossRef]

]

Δε˜=(ΔεrΔεΔεzrΔεΔεθΔεθzΔεzrΔεθzΔεz),
(8)

where Δεr = 2n 1(C 1 σr+C 2 σθ), Δεθ = 2n 1(C 1 σθ+C 2 σr), Δεz = 2n 1 C 2(σr+σθ), Δε = 2n 1(C 1-C 2)τ, and Δεθz = Δεzr = 0. The parameters C1 and C2 are photoelastic constants for pure silica.

Since the elastic modulus of the fiber is very large, its circular symmetry can not be broken by lateral pressure. And the stress-induced perturbation of relative permittivity is much smaller than the difference of the core and cladding refractive indices, so we can assume the eigen modes (two orthogonal solutions of the wave equation) do not vary and the perturbed electromagnetic fields are represented in terms of their linear combination. Then the coupled equations relating two eigen modes can be described as [1

1. D. Chowdhury and D. Wilcox, “Comparison between optical fiber birefringence induced by stress anisotropy and geometric deformation,” IEEE J. Sel. Topics Quantum Electron. 6, 227–232 (2000). [CrossRef]

]

ddz(A1A2)=j(α11α12α21α22)(A1A2),αmn=ωε0sEm*·Δε˜Endssẑ·(Em*×Hm+Em×Hm*)ds,
(9)

where A 1 and A 2 are the field amplitudes for two eigen modes which are only dependent on the fiber axis coordinate z.

Our fiber parameters are given as n1 = 1.449, Δ=n1-n2 = 0.005, v = 2.072, u = 1.552, and w = 1.373. By solving the simultaneous equations (6), (8) and (9) with the benefit of the computational method described in [2

2. J. Sakai and T. Kimura, “Birefringence and polarization characteristics of single mode optical fibers under electric deformations,” IEEE J. Quantum Electron. , 17, 1041–1051 (1981). [CrossRef]

] we obtain α 12=α 21 = 0 and the theoretical birefringence

B=α11α2217CP2[11.166(aR)2].
(10)

In Eq. (10) the stress-optic coefficient C = C 1-C 2 and it is 3.184 × 10−12 m2/N for fused silica [9

9. A. J. Barlow and D. N. Payne, “The stress-optic effect in optical fibers,” IEEE J. Quantum Electron. , 19, 834–839 (1983). [CrossRef]

]. Inserting above values of R, a and C into Eq. (10) we get the birefringence B = 0.2779P rad/m for wavelength λ = 1550nm.

In the case that pressure P squeezes the fiber at different orientation it should be written as a vector P⃗ = PP̂. Correspondingly, we note the birefringence in the vector format B⃗ = BB̂. Concerning the relationship between vectors P⃗ and B⃗, we find when P⃗ rotates around the fiber for angle θ, B⃗ will rotates around a circle for double θ on the Poincaré sphere, which we will prove later experimentally. Since the birefringence vector can be measured by Mueller matrix method (MMM), the orientation angle θ will be determined by

θ=cos1(B̂·B̂0)2,
(11)

in which 0 corresponds to the pressing azimuth at 0°. While the magnitude of the pressure is

P=3.5984B.
(12)

3. Experiment and analysis for pressure vector sensing

Based on Eqs. (11) and (12) we can realize the pressure vector sensing. Our experimental setup is shown schematically in Fig. 3, where the right part is the pressure vector generator (PVG) made of a force gauge and a rotatable PZT-based squeezing device. We use a computer-controlled measure system to obtain the Mueller matrix M of the fiber under different pressure.

Fig. 3. Experimental setup for pressure vector measurement based on Mueller matrix method.

MB=MMD1,MD=m11(1DTDmD),D=DD̂=DDD=m111(m12,m13,m14)TmD=1D2diag[1,1,1]+(11D2)D̂DT
(14)

here mij denotes the component of matrix M. Then the birefringence vector is calculated by

B=BB̂=B(s1,s2,s3)T,B=B=cos1[Tr(MB)21],
s1=(m34Bm43B)(2sinB),s2=(m42Bm24B)(2sinB),s3=(m23Bm32B)(2 sinB),
(15)

where Tr(MB) and mBij(i, j = 1,2,3,4) are the trace and components of matrix MB.

To get the initial value 0, recalling the coordinate system in Fig. 1(b), squeezing the fiber with P = 1 N as shown in the figure while measuring the corresponding birefringence we obtain B 0 = 0.2839 and its unit vector 0=[-0.6218,0.4995,-0.6045]T. Then we change the magnitude of pressure and measure the birefringence while recording the practical pressure values to analyze the error, using Eq. (12) we get the pressure and plot them in Fig. 4(a). By statistics we find the average error between calculated and practical pressure is ~0.17 N.

Furthermore, we rotate the device in Fig. 3 to press the fiber at different orientation determined by geometrical method. Then we calculate birefringence vectors by Eq. (15) and their pressure azimuthes by Eq. (11). Some results are shown in Fig. 4(b), which demonstrates when pressure rotating around the fiber the corresponding birefringence vector rotates around a circle with double speed.

Fig. 4. (a) Measurement and calculation of pressure magnitude, (b) birefringence vectors under different pressure orientations.

To know the validity of our analysis for pressure vector orientation, we compare the calculated and practical (geometrical) azimuthes. Statistical analysis shows the average deviation is ~0.85°. Our experimental results are partly shown in Fig. 5.

Fig. 5. Comparison for calculated and practical azimuthes of lateral pressureure.

4. Conclusions

In summary, based on stress function equation we propose a theoretical method to analyze the pressure stress distribution on the SMF profile, and obtain the analytical expression of the stress function with the stress patterns in the core. The stress-induced birefringence is calculated and applied to realize the pressure vector sensing. We demonstrate the birefringence vector rotates around a circle with double speed as the pressure rotating around the fiber. Experimental results confirm our analysis and the average errors of pressure magnitude and direction are ~0.17 N and ~0.85°. Our conclusions are applicable to high-speed optical communications and to fiber-based pressure sensing.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (grant 60577020 and 60672004) and Open Fund of Key Laboratory of Optical Communication and Lightwave Technologies (BUPT), Ministry of Education, P.R.C.

References and links

1.

D. Chowdhury and D. Wilcox, “Comparison between optical fiber birefringence induced by stress anisotropy and geometric deformation,” IEEE J. Sel. Topics Quantum Electron. 6, 227–232 (2000). [CrossRef]

2.

J. Sakai and T. Kimura, “Birefringence and polarization characteristics of single mode optical fibers under electric deformations,” IEEE J. Quantum Electron. , 17, 1041–1051 (1981). [CrossRef]

3.

Y. Park, U. C. Paek, and D.Y. Kim, “Determination of stress-induced intrinsic birefringence in a single-mode fiber by measurement of the two-dimensional stress profile,” Opt. Lett. , 27, 1291–1293 (2002). [CrossRef]

4.

K.S. Chang, “Pressure-induced birefringence in a coated highly birefringent optical fiber,” J. Lightwave Technol. , 8, 1850–1855 (1990). [CrossRef]

5.

S. M. Pietralunga, M. Ferrario, M. Tacca, and M. Martinelli, “Local Birefringence in Unidirectionally Spun Fibers,” J. Lightwave Technol. , 24, 4030–4038 (2006). [CrossRef]

6.

K. Saitoh, M. Koshiba, and Y. Tsuji, “Stress analysis method for elastically anisotropic material based optical waveguides and its application to strain-induced optical waveguides,” J. Lightwave Technol. , 17, 255–259 (1999). [CrossRef]

7.

S. Timoshenko and J. N. Goodier, Theory of elasticity (McGraw-Hill, 1970), Chap. 4.

8.

D.H. Yu, Mathematical theory of natural boundary element method (Science Press, Beijing, 1993), Chap. 3.

9.

A. J. Barlow and D. N. Payne, “The stress-optic effect in optical fibers,” IEEE J. Quantum Electron. , 19, 834–839 (1983). [CrossRef]

10.

S. Y. Lu and R. A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A 13, 1106–1113 (1996). [CrossRef]

OCIS Codes
(060.2310) Fiber optics and optical communications : Fiber optics
(060.2330) Fiber optics and optical communications : Fiber optics communications
(060.2370) Fiber optics and optical communications : Fiber optics sensors
(260.1440) Physical optics : Birefringence

ToC Category:
Fiber Optics and Optical Communications

History
Original Manuscript: January 29, 2008
Revised Manuscript: March 5, 2008
Manuscript Accepted: March 5, 2008
Published: March 10, 2008

Citation
Zhengyong Li, Chongqing Wu, Hui Dong, P. Shum, C. Y. Tian, and S. Zhao, "Stress distribution and induced birefringence analysis for pressure vector sensing based on single mode fibers," Opt. Express 16, 3955-3960 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-6-3955


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. D. Chowdhury and D. Wilcox, "Comparison between optical fiber birefringence induced by stress anisotropy and geometric deformation," IEEE J. Sel. Topics Quantum Electron. 6, 227-232 (2000). [CrossRef]
  2. J. Sakai and T. Kimura, "Birefringence and polarization characteristics of single mode optical fibers under electric deformations," IEEE J. Quantum Electron. 17, 1041-1051 (1981). [CrossRef]
  3. Y. Park, U. C. Paek, and D.Y. Kim, "Determination of stress-induced intrinsic birefringence in a single-mode fiber by measurement of the two-dimensional stress profile," Opt. Lett. 27, 1291-1293 (2002). [CrossRef]
  4. K.S. Chang, "Pressure-induced birefringence in a coated highly birefringent optical fiber," J. Lightwave Technol. 8, 1850-1855 (1990). [CrossRef]
  5. S. M. Pietralunga, M. Ferrario, M. Tacca, and M. Martinelli, "Local Birefringence in Unidirectionally Spun Fibers," J. Lightwave Technol. 24, 4030-4038 (2006). [CrossRef]
  6. K. Saitoh, M. Koshiba, and Y. Tsuji, "Stress analysis method for elastically anisotropic material based optical waveguides and its application to strain-induced optical waveguides," J. Lightwave Technol. 17, 255-259 (1999). [CrossRef]
  7. S. Timoshenko and J. N. Goodier, Theory of elasticity (McGraw-Hill, 1970), Chap. 4.
  8. D.H. Yu, Mathematical theory of natural boundary element method (Science Press, Beijing, 1993), Chap. 3.
  9. A. J. Barlow and D. N. Payne, "The stress-optic effect in optical fibers," IEEE J. Quantum Electron. 19, 834-839 (1983). [CrossRef]
  10. S. Y. Lu and R. A. Chipman, "Interpretation of Mueller matrices based on polar decomposition," J. Opt. Soc. Am. A 13, 1106-1113 (1996). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited