PhotoElasticCorrectionΒΆ

Type:section
Appearance:multiple

The photo-elastic correction relates the stresses within a body to a change of the relative permittivity. To formulate this correction term it is necessary to bring the stress tensor \TField{\sigma}, given by a symmetric 3 \times 3 - matrix, into diagonal form using a coordinate rotation \pvec{x}'=\TField{R}^{\mathrm{T}} \pvec{x}:

\begin{eqnarray*}
\TField{\sigma} & = & \TField{R} \cdot \left [ \begin{array}{ccc}
             \sigma_{x'x'} & 0 & 0 \\
             0 & \sigma_{y'y'} & 0 \\
             0 & 0 & \sigma_{z'z'}
             \end{array} \right] \cdot \TField{R}^{\mathrm{T}}
\end{eqnarray*}

The columns of the rotation matrix \TField{R} are the three principle axes of the stress tensor and \sigma_{x'x'}, \sigma_{y'y'} and \sigma_{z'z'} are the stresses into the directions of the principle axes.

Within the x'y'z' - coordinate system the photo-elastic correction is in diagonal form and is given as follows:

\begin{eqnarray*}
\Delta \TField{\varepsilon}^{(\mathrm{rel})}_{x'x'} & = & -C_1 \sigma_{x'x'}-C_2\sigma_{y'y'}-C_2\sigma_{z'z'} \\
\Delta \TField{\varepsilon}^{(\mathrm{rel})}_{y'y'} & = & -C_2 \sigma_{x'x'}-C_1\sigma_{y'y'}-C_2\sigma_{z'z'} \\
\Delta \TField{\varepsilon}^{(\mathrm{rel})}_{z'z'} & = & -C_2 \sigma_{x'x'}-C_2\sigma_{y'y'}-C_1\sigma_{z'z'}
\end{eqnarray*}

C_1 and C_2 are called photo-elastic constants. The first coefficient C_1 relates the stress in the direction of a principle axis to a permittivity change in this direction, whereas C_2 gives rise to a permittivity correction perpendicular to the considered principle axis.

Transforming to the original xyz' - coordinate system we get that

\begin{eqnarray*}
\Delta \TField{\varepsilon}^{(\mathrm{rel})} & = & \TField{R} \cdot \left [ \begin{array}{ccc}
             \Delta \TField{\varepsilon}^{(\mathrm{rel})}_{x'x'} & 0 & 0 \\
             0 & \Delta \TField{\varepsilon}^{(\mathrm{rel})}_{y'y'} & 0 \\
             0 & 0 & \Delta \TField{\varepsilon}^{(\mathrm{rel})}_{z'z'}
             \end{array} \right] \cdot \TField{R}^{\mathrm{T}}.
\end{eqnarray*}

The full relative permittivity is then given by

\begin{eqnarray*}
\TField{\varepsilon}^{(\mathrm{rel})} & =  \TField{\varepsilon}^{(\mathrm{rel})}_{0} + \Delta \TField{\varepsilon}^{(\mathrm{rel})},
\end{eqnarray*}

where \TField{\varepsilon}^{(\mathrm{rel})}_{0} is the relative permittivity of the stress free body.

For too large stresses the linear correction term may loose its validity. Especially, extremely large stresses occurs near corners of the geometry, where the stress may exhibit field singularities. To prevent obscure relative permittivity definitions one may restrict the absolute value of the correction term by a threshold value \delta_{\max},

\begin{eqnarray*}
\Delta \TField{\varepsilon}^{(\mathrm{rel})}_{ij} & =  & \min \left ( \max \left ( \Delta \TField{\varepsilon}^{(\mathrm{rel})}_{ij}, -\delta_{\max} \right ), \delta_{\max} \right ),
\end{eqnarray*}

The JCM - syntax looks like this:

# define the relative permittivity with photo-elastic correction
RelPermittivity {
  # define relative permittivity of stress free body
  Constant = ...
  ...
  PhotoElasticCorrection {
     C1 = ... # set first photo-elastic coefficient
     C2 = ... # set second photo-elastic coefficient
     CutOff .. # set threshold here
  }
}

Warning

Often, the photo-elastic coefficient is defined as a refractive index - stress relation.

In a linear sense the refractive index based definition is equivalent to the definition in JCMsuite: In the x'y'z' - coordinate system the refractive index change be given by

\begin{eqnarray*}
\Delta n_{x'x'} = -\widetilde{C}_1 \sigma_{x'x'}-\widetilde{C}_2 \sigma_{y'y'}-\widetilde{C}_2 \sigma_{z'z'},
\end{eqnarray*}

and corresponding expression for \Delta n_{x'x'} and \Delta n_{y'y'}.

With the refractive index n_0 of the stress-free body this results in the following definition of the relative permittivity:

\begin{eqnarray*}
\TField{\varepsilon}^{(\mathrm{rel})}_{x'x'} & = & (n_{0}+\Delta n_{x'x'})^2 =
(n_{0}-\widetilde{C}_1 \sigma_{x'x'}-\widetilde{C}_2 \sigma_{y'y'}
-\widetilde{C}_2 \sigma_{z'z'})^2 \\
\TField{\varepsilon}^{(\mathrm{rel})}_{x'x'} & = & n_0^2-2 n_0 \widetilde{C}_1 \sigma_{x'x'}-2 n_0
\widetilde{C}_2 \sigma_{y'y'} - 2 \widetilde{C}_2 \sigma_{z'z'}+\mathcal{O} (|\sigma|^2) \\
{ } & \sim & \TField{\varepsilon}^{(\mathrm{rel})}_{0}-C_1  \sigma_{x'x'}-C_2  \sigma_{y'y'} -C_2  \sigma_{z'z'}
\end{eqnarray*}

The approximation in the last line is due to linearization. The last line is precisely the definition of the photo-elastic relation in JCMsuite when setting C_1=2n_0\widetilde{C}_1 and C_2=2n_0\widetilde{C}_2.