# Vibration Analysis of an Optical Fiber Coupler Using the Differential Quadrature Method

Authors: Ming-Hung Hsu
Journal: Mathematical Problems in Engineering
 (2010)

---

## Abstract

This study elucidates the dynamic characteristics of an optical fiber coupler via the differential quadrature simulations. A novel modeling scheme is suitable for developing an optical fiber coupler application. Exactly how the locations of bonding points, string tension, and spring stiffness of the rubber pad affect the dynamic behavior of the optical fiber coupler are investigated.


---

## Body

## 1. Introduction

An optical fiber coupler is an optical device with several input fibers and several output fibers. The optical fibers in couplers may be under shock and impact. Cheng and Zu [1] and Sun [2] studied  vibration of an optical fiber coupler subjected to a half-sine shock. Malomed and Tasgal [3]  analyzed the dynamics of small internal vibrations in a two-component gap soliton. They found three oscillation modes, which are composed of dilation-contraction of each component’s width, and a relative translation of the two components. Brown et al. [4] performed vibration tests on commercial grade fiber optic connectors and splices. Huang et al. [5] presented optical coupling loss and vibration characterization for packaging of 2×2 MEMS vertical torsion mirror switches.  Thomes et al.  [6] presented vibration performance of current fiber optic connector.In this study the idea of differential quadrature formulation is extended to an optical fiber coupler. During the last decade, the differential quadrature approach applied to engineering and science problems has attracted considerable attention [7–31]. Liew et al. [7–16] applied the differential quadrature method to Mindlin plates on Winkler foundations and developed an application of the differential quadrature method to thick symmetric cross-ply laminates with first-order shear flexibility. Liew et al. [7–16] also employed the generalized differential quadrature method for buckling analysis and examined static and free vibration of beams and rectangular and annular plates. Sherbourne and Pandey [17] investigated the buckling behavior of beams and composite plates using the differential quadrature method. Mirfakhraei and Redekop [18] evaluated the buckling of circular cylindrical shells using the differential quadrature method. Tomasiello [19] applied the differential quadrature method to evaluate initial-boundary-value problems. Moradi and Taheri [20] conducted buckling analysis of general laminated composite beams using the differential quadrature method. De Rosa and Franciosi [21] solved the dynamic problem of circular arches using the differential quadrature method. Sun and Zhu [22] investigated incompressible viscous flow problems using the differential quadrature method. Via the differential quadrature method, Tanaka and Chen [23] solved transient elastodynamic problems. Chen and Zhong [24] observed that the differential quadrature and differential cubature methods, due to their global domain characteristics, are more efficient in solving nonlinear problems than conventional numerical schemes, such as the finite element and finite difference methods. Civan [25] solved multivariable mathematical models using the differential quadrature and differential cubature methods. Hua and Lam [26] identified the frequency characteristics of a thin rotating cylindrical shell using the differential quadrature method. Wang et al. [27–30] employed new versions of the differential quadrature and differential quadrature element methods to analyze anisotropic rectangular plates, frame structures, nonuniform beams, circular annular plates, and isotropic skew plates. The dynamic behavior of an optical fiber coupler is elucidated using the differential quadrature method in this work. Few studies have conducted vibration analysis of an optical fiber coupler using the differential quadrature method.


## 2. Differential Quadrature Method

Solutions to numerous complex beam problems have been efficiently acquired using fast computers and various numerical schemes, including the Galerkin technique, finite element method, boundary element method, and Rayleigh-Ritz method. In this study, the differential quadrature scheme is employed to generate discrete eigenvalue problems for an optical fiber coupler. The basic concept of the differential quadrature method is that the derivative of a function at a given point can be approximated as a weighted linear sum of functional values at all sample points in the domain of that variable. The partial differential equation is then reduced to a set of algebraic equations. For a function,f(x), the differential quadrature approximation for the mth-order derivative at the ith sample point is given by
(2.1)[dmf(x1)dxmdmf(x2)dxm⋮dmf(xN)dxm]≅[Dij(m)][f(x1)f(x2)⋮f(xN)]fori,j=1,2,…,N,
where f(xi) is the value of the function at sample point xi,Dij(m) is the weighted coefficient of the mth-order differentiation attached to these functional values, N is the number of sample points, and xi is the location of the ith sample point in the domain. The most convenient technique is to distribute sample points uniformly [31]. A Lagrangian interpolation polynomial is utilized to eliminate possible adverse conditions when determining the weighted coefficients Dij(m) [31], which are as follows:(2.2)f(x)=M(x)(x-xi)M1(xi)fori=1,2,…,N,
where
(2.3)M(x)=∏j=1N(x-xj),M1(xi)=∏j=1,j≠iN(xi-xj)fori=1,2,…,N.
Inputting (2.2) into (2.1) yields
(2.4)Dij(1)=M1(xi)(xi-xj)M1(xj)fori,j=1,2,…,N,i≠j,Dii(1)=-∑j=1,j≠iNDij(1)fori=1,2,…,N.
The coefficients of the weighted matrix can be acquired using (2.4). For the mth-order derivative, the weighted coefficients can be obtained using the following recurrence relation equations:
(2.5)Dij(m)=∑k=1NDik(1)Dkj(m-1)fori,j=1,2,…,N.
The selection of sample points always is very important to solution accuracy when using the differential quadrature approach. For a beam problem, the most convenient technique is to choose equally spaced sample points [31]. Unequally spaced sample points, such as Chebyshev-Gauss-Lobatto sample points, have been utilized by a number of studies. With the Chebyshev-Gauss-Lobatto distribution, the sample points of an optical fiber coupler are distributed as(2.6)xi=12(1-cos(i-1)πN-1)fori=1,2,…,N.


## 3. Vibration of an Optical Fiber Coupler with a Continuous Elastic Support

Figure1 shows a sectional view of the optical fiber coupler with a continuous elastic support. The fibers, the substrate, and the rubber pad are on the inside of the steel tube. The optical fibers are placed onto the substrate. A beam represents the substrate. The optical fibers model as a string. The rubber pad is placed between the substrate and steel tube. The linear model assumes that fiber tension is constant. The equation of motion for the optical fiber coupler is derived as [1, 2](3.1)P∂2u1∂x2-ρ1A1∂2u1∂t2=0forx:0tol1,P∂2u2∂x2-ρ1A1∂2u2∂t2=0forx:l1tol1+l2,P∂2u3∂x2-ρ1A1∂2u3∂t2=0forx:l1+l2tol1+l2+l3,E2I2∂4v1∂x4+kfv1+ρ2A2∂2v1∂t2=0forx:0tol1,E2I2∂4v2∂x4+kfv2+ρ2A2∂2v2∂t2=0forx:l1tol1+l2,E2I2∂4v3∂x4+kfv3+ρ2A2∂2v3∂t2=0forx:l1+l2tol1+l2+l3,
where u1, u2, and u3 are displacements of the fibers; v1, v2, and v3 are displacements of the substrate; P is string tension; t is time; kf is the constant determined by the material constants of the silicon rubber pad; A1 is cross-section area of fibers; A2 is cross-section area of substrate; ρ1 is density of the fiber material; ρ2 is density of the substrate material; I2 is the second moment of the cross-sectional area A2; E2 is Young’s modulus of the substrate. The optical fibers and substrate are bonded at four points. The boundary conditions of the optical fiber coupler are(3.2)u1(0,t)-v1(0,t)=0,u1(l1,t)-v1(l1,t)=0,u2(l1,t)-v2(l1,t)=0,u2(l1+l2,t)-v2(l1+l2,t)=0,u3(l1+l2,t)-v3(l1+l2,t)=0,u3(l1+l2+l3,t)-v3(l1+l2+l3,t)=0,E2I2∂2v1(0,t)∂x2=0,E2I2∂3v1(0,t)∂x3=0,v1(l1,t)-v2(l1,t)=0,∂v1(l1,t)∂x-∂v2(l1,t)∂x=0,E2I2∂2v1(l1,t)∂x2-E2I2∂2v2(l1,t)∂x2=0,E2I2∂3v1(l1,t)∂x3-E2I2∂3v2(l1,t)∂x3=0,v2(l1+l2,t)-v3(l1+l2,t)=0,∂v2(l1+l2,t)∂x-∂v3(l1+l2,t)∂x=0,E2I2∂2v2(l1+l2,t)∂x2-E2I2∂2v3(l1+l2,t)∂x2=0,E2I2∂3v2(l1+l2,t)∂x3-E2I2∂3v3(l1+l2,t)∂x3=0,E2I2∂2v3(l1+l2+l3,t)∂x2=0,E2I2∂3v3(l1+l2+l3,t)∂x3=0.
To obtain frequencies, a harmonic movement of the optical fiber coupler is assumed as
(3.3)us(x,t)=u̅s(x)cos(ωt)fors=1,2,3,vs(x,t)=v̅s(x)cos(ωt)fors=1,2,3,
where u̅1(x),u̅2(x),u̅3(x),v̅1(x),v̅2(x), and v̅3(x) are vibrational modes, and ω is the natural frequency of the optical fiber coupler. Substituting (3.3) into (3.1) yields(3.4)Pd2u̅sdx2=-ω2ρ1A1u̅sfors=1,2,3,E2I2d4v̅sdx4+kfv̅s=ω2ρ2A2v̅sfors=1,2,3.
The boundary conditions of the optical fiber coupler are rewritten as
(3.5)u̅1(0)-v̅1(0)=0,u̅1(l1)-v̅1(l1)=0,u̅2(l1)-v̅2(l1)=0,u̅2(l1+l2)-v̅2(l1+l2)=0,u̅3(l1+l2)-v̅3(l1+l2)=0,u̅3(l1+l2+l3)-v̅3(l1+l2+l3)=0,E2I2d2v̅1(0)dx2=0,E2I2d3v̅1(0)dx3=0,v̅1(l1)-v̅2(l1)=0,dv1(l1)dx-dv2(l1)dx=0,E2I2d2v̅1(l1)dx2-E2I2d2v̅2(l1)dx2=0,E2I2d3v̅1(l1)dx3-E2I2d3v̅2(l1)dx3=0,v̅2(l1+l2)-v̅3(l1+l2)=0,dv̅2(l1+l2)dx-dv̅3(l1+l2)dx=0,E2I2d2v̅2(l1+l2)dx2-E2I2d2v̅3(l1+l2)dx2=0,E2I2d3v̅2(l1+l2)dx3-E2I2d3v̅3(l1+l2)dx3=0,E2I2d2v̅3(l1+l2+l3)dx2=0,E2I2d3v̅3(l1+l2+l3)dx3=0.
The equations of motion of the optical fiber coupler can be rearranged in the differential quadrature method formula by substituting (2.1) into (3.4) and (3.5). The equations of motion of the optical fiber coupler are derived as(3.6)∑j=1NPDi,j(2)ls2u̅s,j=-ω2ρ1A1u̅s,ifori=1,2,…,N,s=1,2,3,∑j=1NE2I2Di,j(4)ls4v̅s,j+kfv̅s,i=ω2ρ2A2v̅s,ifori=1,2,…,N,s=1,2,3.
Using the differential quadrature method, the boundary conditions of the optical fiber coupler can be rearranged into the matrix form as
(3.7)u̅1,1-v̅1,1=0,u̅1,N-v̅1,N=0,u̅2,1-v̅2,1=0,u̅2,N-v̅2,N=0,u̅3,1-v̅3,1=0,u̅3,N-v̅3,N=0,∑j=1NE2I2D1,j(2)l12v̅1,j=0,∑j=1NE2I2D1,j(3)l13v̅1,j=0,v̅1,N-v̅2,1=0,∑j=1NDN,j(1)l1v̅1,j-∑j=1ND1,j(1)l2v̅2,j=0,∑j=1NE2I2DN,j(2)l12v̅1,j-∑j=1NE2I2D1,j(2)l22v̅2,j=0,∑j=1NE2I2DN,j(3)l13v̅1,j-∑j=1NE2I2D1,j(3)l23v̅2,j=0,v̅2,N-v̅3,1=0,∑j=1NDN,j(1)l2v̅2,j-∑j=1ND1,j(1)l3v̅3,j=0,∑j=1NE2I2DN,j(2)l22v̅2,j-∑j=1NE2I2D1,j(2)l32v̅3,j=0,∑j=1NE2I2DN,j(3)l23v̅2,j-∑j=1NE2I2D1,j(3)l33v̅3,j=0,∑j=1NE2I2DN,j(2)l32v̅3,j=0,∑j=1NE2I2DN,j(3)l33v̅3,j=0.Figure 1
Sectional view of the optical fiber coupler with an elastic foundation [1, 2].


## 4. The Optical Fiber Coupler with Two Spring Supports

Figure2 shows a sectional view of the optical fiber coupler with two rubber pads at each end of the coupler. The rubber pads are placed between the substrate and steel tube. The rubber pads model as the two spring supports. The equations of motion for the optical fiber coupler are [1, 2]
(4.1)P∂2u1∂x2-ρ1A1∂2u1∂t2=0forx:0tol1,P∂2u2∂x2-ρ1A1∂2u2∂t2=0forx:l1tol1+l2,P∂2u3∂x2-ρ1A1∂2u3∂t2=0forx:l1+l2tol1+l2+l3,E2I2∂4v1∂x4+ρ2A2∂2v1∂t2=0forx:0tol1,E2I2∂4v2∂x4+ρ2A2∂2v2∂t2=0forx:l1tol1+l2,E2I2∂4v3∂x4+ρ2A2∂2v3∂t2=0forx:l1+l2tol1+l2+l3.
The boundary conditions of the optical fiber coupler are
(4.2)u1(0,t)-v1(0,t)=0,u1(l1,t)-v1(l1,t)=0,u2(l1,t)-v2(l1,t)=0,u2(l1+l2,t)-v2(l1+l2,t)=0,u3(l1+l2,t)-v3(l1+l2,t)=0,u3(l1+l2+l3,t)-v3(l1+l2+l3,t)=0,E2I2∂2v1(0,t)∂x2=0,E2I2∂3v1(0,t)∂x3=-kspringv1(0,t),v1(l1,t)-v2(l1,t)=0,∂v1(l1,t)∂x-∂v2(l1,t)∂x=0,E2I2∂2v1(l1,t)∂x2-E2I2∂2v2(l1,t)∂x2=0,E2I2∂3v1(l1,t)∂x3-E2I2∂3v2(l1,t)∂x3=0,v2(l1+l2,t)-v3(l1+l2,t)=0,∂v2(l1+l2,t)∂x-∂v3(l1+l2,t)∂x=0,E2I2∂2v2(l1+l2,t)∂x2-E2I2∂2v3(l1+l2,t)∂x2=0,E2I2∂3v2(l1+l2,t)∂x3-E2I2∂3v3(l1+l2,t)∂x3=0,E2I2∂2v3(l1+l2+l3,t)∂x2=0,E2I2∂3v3(l1+l2+l3,t)∂x3=kspringv3(l1+l2+l3,t),
where kspring is the constant determined by material constants of the spring support. Substituting (3.3) into (4.1) yields(4.3)Pd2u̅sdx2=-ω2ρ1A1u̅sfors=1,2,3,E2I2d4v̅sdx4=ω2ρ2A2v̅sfors=1,2,3.
The boundary conditions of the optical fiber coupler are rewritten as
(4.4)u̅1(0)-v̅1(0)=0,u̅1(l1)-v̅1(l1)=0,u̅2(l1)-v̅2(l1)=0,u̅2(l1+l2)-v̅2(l1+l2)=0,u̅3(l1+l2)-v̅3(l1+l2)=0,u̅3(l1+l2+l3)-v̅3(l1+l2+l3)=0,E2I2d2v̅1(0)dx2=0,E2I2d3v̅1(0)dx3=-kspringv̅1(0),v̅1(l1)-v̅2(l1)=0,dv1(l1)dx-dv2(l1)dx=0,E2I2d2v̅1(l1)dx2-E2I2d2v̅2(l1)dx2=0,E2I2d3v̅1(l1)dx3-E2I2d3v̅2(l1)dx3=0,v̅2(l1+l2)-v̅3(l1+l2)=0,dv̅2(l1+l2)dx-dv̅3(l1+l2)dx=0,E2I2d2v̅2(l1+l2)dx2-E2I2d2v̅3(l1+l2)dx2=0,E2I2d3v̅2(l1+l2)dx3-E2I2d3v̅3(l1+l2)dx3=0,E2I2d2v̅3(l1+l2+l3)dx2=0,E2I2d3v̅3(l1+l2+l3)dx3=kspringv̅3(l1+l2+l3).
The equations of motion of the optical fiber coupler can be rearranged in the differential quadrature method formula by substituting (2.1) into (4.3). The equations of motion of the optical fiber coupler then become(4.5)∑j=1NPDi,j(2)ls2u̅s,j=-ω2ρ1A1u̅s,ifori=1,2,…,N,s=1,2,3,∑j=1NE2I2Di,j(4)ls4v̅s,j=ω2ρ2A2v̅s,ifori=1,2,…,N,s=1,2,3.
Using the differential quadrature method, the boundary conditions of the optical fiber coupler can be rearranged into the matrix form as
(4.6)u̅1,1-v̅1,1=0,u̅1,N-v̅1,N=0,u̅2,1-v̅2,1=0,u̅2,N-v̅2,N=0,u̅3,1-v̅3,1=0,u̅3,N-v̅3,N=0,∑j=1NE2I2D1,j(2)l12v̅1,j=0,∑j=1NE2I2D1,j(3)l13v̅1,j=-kspringv̅1,1,v̅1,N-v̅2,1=0,∑j=1NDN,j(1)l1v̅1,j-∑j=1ND1,j(1)l2v̅2,j=0,∑j=1NE2I2DN,j(2)l12v̅1,j-∑j=1NE2I2D1,j(2)l22v̅2,j=0,∑j=1NE2I2DN,j(3)l13v̅1,j-∑j=1NE2I2D1,j(3)l23v̅2,j=0,v̅2,N-v̅3,1=0,∑j=1NDN,j(1)l2v̅2,j-∑j=1ND1,j(1)l3v̅3,j=0,∑j=1NE2I2DN,j(2)l22v̅2,j-∑j=1NE2I2D1,j(2)l32v̅3,j=0,∑j=1NE2I2DN,j(3)l23v̅2,j-∑j=1NE2I2D1,j(3)l33v̅3,j=0,∑j=1NE2I2DN,j(2)l32v̅3,j=0,∑j=1NE2I2DN,j(3)l33v̅3,j=kspringv̅3,N.Figure 2
Sectional view of the optical fiber coupler with two spring supports [1, 2].


## 5. Nonlinear Dynamic Analysis of an Optical Fiber Coupler with  a Continuous Elastic Support

The nonlinear model assumes that fiber tension varies. The coupler is subjected to a sine shock motion. The equations of motion for the optical fiber coupler with a continuous elastic support are [1, 2](5.1)P∂2u1∂x2+kstring(∫0l11+(∂u1∂x)2dx+∫l1l1+l21+(∂u2∂x)2dx+∫l1+l2l1+l2+l31+(∂u3∂x)2dx-l1-l2-l3)-ρ1A1∂2u1∂t2=ρ1A1asin(Ωt)forx:0tol1,P∂2u2∂x2+kstring(∫0l11+(∂u1∂x)2dx+∫l1l1+l21+(∂u2∂x)2dx+∫l1+l2l1+l2+l31+(∂u3∂x)2dx-l1-l2-l3)-ρ1A1∂2u2∂t2=ρ1A1asin(Ωt)forx:l1tol1+l2,P∂2u3∂x2+kstring(∫0l11+(∂u1∂x)2dx+∫l1l1+l21+(∂u2∂x)2dx+∫l1+l2l1+l2+l31+(∂u3∂x)2dx-l1-l2-l3)-ρ1A1∂2u3∂t2=ρ1A1asin(Ωt)forx:l1+l2tol1+l2+l3,(5.2)E2I2∂4v1∂x4+kfv1+ρ2A2∂2v1∂t2=-ρ2A2asin(Ωt)forx:0tol1,E2I2∂4v2∂x4+kfv2+ρ2A2∂2v2∂t2=-ρ2A2asin(Ωt)forx:l1tol1+l2,E2I2∂4v3∂x4+kfv3+ρ2A2∂2v3∂t2=-ρ2A2asin(Ωt)forx:l1+l2tol1+l2+l3,
where kstring is the elastic coefficient of the string, a is the acceleration of shock motion, and Ω is the circular frequency of shock motion. With the following approximation equation
(5.3)1+(∂us∂x)2≈1+12(∂us∂x)2fors=1,2,3(5.1) can be rewritten as [1, 2]
(5.4)P∂2u1∂x2+kstring2(∫0l1(∂u1∂x)2dx+∫l1l1+l2(∂u2∂x)2dx+∫l1+l2l1+l2+l3(∂u3∂x)2dx)-ρ1A1∂2u1∂t2=ρ1A1asin(Ωt)forx:0tol1,P∂2u2∂x2+kstring2(∫0l1(∂u1∂x)2dx+∫l1l1+l2(∂u2∂x)2dx+∫l1+l2l1+l2+l3(∂u3∂x)2dx)-ρ1A1∂2u2∂t2=ρ1A1asin(Ωt)forx:l1tol1+l2,P∂2u3∂x2+kstring2(∫0l1(∂u1∂x)2dx+∫l1l1+l2(∂u2∂x)2dx+∫l1+l2l1+l2+l3(∂u3∂x)2dx)-ρ1A1∂2u3∂t2=ρ1A1asin(Ωt)forx:l1+l2tol1+l2+l3.
The boundary conditions of the optical fiber coupler are
(5.5)u1(0,t)-v1(0,t)=0,u1(l1,t)-v1(l1,t)=0,u2(l1,t)-v2(l1,t)=0,u2(l1+l2,t)-v2(l1+l2,t)=0,u3(l1+l2,t)-v3(l1+l2,t)=0,u3(l1+l2+l3,t)-v3(l1+l2+l3,t)=0,E2I2∂2v1(0,t)∂x2=0,E2I2∂3v1(0,t)∂x3=0,v1(l1,t)-v2(l1,t)=0,∂v1(l1,t)∂x-∂v2(l1,t)∂x=0,E2I2∂2v1(l1,t)∂x2-E2I2∂2v2(l1,t)∂x2=0,E2I2∂3v1(l1,t)∂x3-E2I2∂3v2(l1,t)∂x3=0,v2(l1+l2,t)-v3(l1+l2,t)=0,∂v2(l1+l2,t)∂x-∂v3(l1+l2,t)∂x=0,E2I2∂2v2(l1+l2,t)∂x2-E2I2∂2v3(l1+l2,t)∂x2=0,E2I2∂3v2(l1+l2,t)∂x3-E2I2∂3v3(l1+l2,t)∂x3=0,E2I2∂2v3(l1+l2+l3,t)∂x2=0,E2I2∂3v3(l1+l2+l3,t)∂x3=0.
The equation of motion of the optical fiber coupler can be rearranged in the differential quadrature method formula by substituting (2.1) into (5.2) and (5.4). The equations of motion of the optical fiber coupler are
(5.6)∑j=1NPDi,j(2)ls2us,j+kstring2(∫0l1(∂u1∂x)2dx+∫l1l1+l2(∂u2∂x)2dx+∫l1+l2l1+l2+l3(∂u3∂x)2dx)-ρ1A1∂2us,i∂t2=ρ1A1asin(Ωt)fori=1,2,…,N,s=1,2,3,∑j=1NE2I2Di,j(4)ls4vs,j+kfvs,i+ρ2A2∂2vs,i∂t2=-ρ2A2asin(Ωt)fori=1,2,…,N,s=1,2,3.
Using the differential quadrature method, the boundary conditions of the optical fiber coupler can be rearranged into the matrix form as
(5.7)u1,1-v1,1=0,u1,N-v1,N=0,u2,1-v2,1=0,u2,N-v2,N=0,u3,1-v3,1=0,u3,N-v3,N=0,∑j=1NE2I2D1,j(2)l12v1,j=0,∑j=1NE2I2D1,j(3)l13v1,j=0,v1,N-v2,1=0,∑j=1NDN,j(1)l1v1,j-∑j=1ND1,j(1)l2v2,j=0,∑j=1NE2I2DN,j(2)l12v1,j-∑j=1NE2I2D1,j(2)l22v2,j=0,∑j=1NE2I2DN,j(3)l13v1,j-∑j=1NE2I2D1,j(3)l23v2,j=0,v2,N-v3,1=0,∑j=1NDN,j(1)l2v2,j-∑j=1ND1,j(1)l3v3,j=0,∑j=1NE2I2DN,j(2)l22v2,j-∑j=1NE2I2D1,j(2)l32v̅3,j=0,∑j=1NE2I2DN,j(3)l23v2,j-∑j=1NE2I2D1,j(3)l33v3,j=0,∑j=1NE2I2DN,j(2)l32v3,j=0,∑j=1NE2I2DN,j(3)l33v3,j=0.


## 6. Numerical Results and Discussion

Figure3presents the natural frequencies of the optical fiber coupler for various values of P. The material and geometric parameters of the optical fiber coupler are kf=50000N/m2, A1=3.1×10-8m2, A2=6.61×10-6m2, ρ1=2.2×103kg/m3, ρ2=2.2×103kg/m3, I2=4.34×10-12m4, l1=0.1333m, l2=0.1333m, l3=0.1333m, and E2=7.24×1010N/m2 [1, 2]. The first and second natural frequencies of the optical fiber coupler are robust to the string tension, P. The third and fourth natural frequencies of the optical fiber coupler increase as the string tension, P, increases. Figure 4 shows the natural frequencies of the optical fiber coupler for various values of kf. The material and geometric parameters of the optical fiber coupler are P=0.01N, A1=3.1×10-8m2, A2=6.61×10-6m2, ρ1=2.2×103kg/m3, ρ2=2.2×103kg/m3, I2=4.34×10-12m4, l1=0.1333m, l2=0.1333m, l3=0.1333m, and E2=7.24×1010N/m2 [1, 2]. The first and third natural frequencies of the optical fiber coupler increase as the rubber pad stiffness increases. The rubber pad stiffness does not significantly affect the second and fourth natural frequencies of the optical fiber coupler. Figure 5 lists the natural frequencies of the optical fiber coupler with bonding points at various locations. The material and geometric parameters of the optical fiber coupler are P=0.01N, kf=50000N/m2, A1=3.1×10-8m2, A2=6.61×10-6m2, ρ1=2.2×103kg/m3, ρ2=2.2×103kg/m3, I2=4.34×10-12m4, l1+l2+l3=0.4m, and E2=7.24×1010N/m2 [1, 2]. The fourth natural frequency of the optical fiber coupler generally increases rapidly as lengths l1 and l3 increase. The locations of bonding points markedly impact the second and third natural frequencies of the optical fiber coupler. Figure 6 plots the natural frequencies of the optical fiber coupler for various values of kspring. The material and geometric parameters of the optical fiber coupler are P=0.01N, A1=3.1×10-8m2, A2=6.61×10-6m2, ρ1=2.2×103kg/m3, ρ2=2.2×103kg/m3, I2=4.34×10-12m4, l1=0.1333m, l2=0.1333m, l3=0.1333m, and E2=7.24×1010N/m2 [1, 2]. The spring constant, kspring, does not affect the first, third and fourth natural frequencies of the optical fiber coupler. Notably, the spring constant, kspring, increases the second natural frequencies of the optical fiber coupler. Figures 7 and 8 show the displacements of the center of the fibers and the substrate under a shock, respectively. The material and geometric parameters of the optical fiber coupler are P=0.01N, A1=3.1×10-8m2, A2=6.61×10-6m2, ρ1=2.2×103kg/m3, ρ2=2.2×103kg/m3, I2=4.34×10-12m4, l1=0.1333m, l2=0.1333m, l3=0.1333m, E2=7.24×1010N/m2, and kstring=5000N/m [1, 2]. The fibers and substrate stiffen when the foundation stiffness, kf, is large. The differential quadrature method is effective in treating this problem.Figure 3
Natural frequencies of the optical fiber coupler for various values ofP.Figure 4
Natural frequencies of the optical fiber coupler for various values ofkf.Figure 5
Natural frequencies of the optical fiber coupler with bonding points at various locations.Figure 6
Natural frequencies of the optical fiber coupler for various values ofkspring.Figure 7
Displacements at the fiber center for various values ofkf.Figure 8
Displacements at the substrate center for various values ofkf.


## 7. Conclusions

This study demonstrates the value of the differential quadrature method for vibration analysis of an optical fiber coupler. The effects of string tensionP, bonding locations, surrounding medium, spring constant kspring, and rubber pad stiffness kf on the natural frequencies of the optical fiber coupler are discussed. The effect of stiffness of the silicon rubber pad during vibrations is significant and should be incorporated into the designs of the optical fiber couplers.



---

*Source: 325656-2010-05-17.xml*