Stress Free Edge Singularity (Stroh's Method): NPIB

Submit batch job to compute singularities at an interface seperating two
dissimilar anisotropic materials that terminates at a laminate stress free edge

Instructions for using this form:

    1. Enter appropriate numbers in the boxes,
    2. enter your email address,
    3. click on submit, and
    4. wait for results to be returned to your email. (this will take about 10-seconds)
  1. Web location of results: strohedge.html will be sent to your email address.
  2. If your email does not provide direct email links, copy and paste the web address: http://www.jwave.vt.edu/output/strohedge_"unique date-time"/strohedge.html from the email
    into your web browser Location: window and look at this file to see if the job ran correctly.
  3. If the job ran correctly you can now download files from your web browser to your computer for archiving (a 4.6Gbyte read/write optical disk is available for each student to use in the SMVC for archiving large simulation files).
  4. If you submit the form "as-is", a sample of results has been archived in the directory: http://www.jwave.vt.edu/output/ARCHIVE_Examples_SAVE/strohedge_12-18-2000-21:09:023:755/
    Select strohedge.html to see a summary of the results.

Below is a brief description of a method to calculate singuariltes that exists at an interface seperating two dissimilar anisotropic materials that terminates at a laminate stress free edge. A more complete description is give in Section 7. "Cracks Near Interfaces Between Dissimilar Anisotropic Materials"


Figure 1. Coordinates defined for stress free edge singularity Ref.[25]

It is important to note that singularities coincide with the coordinate system origin. Hence each interface singularity has a unique coordinate system. Since the objective here is to reproduce published results the same coordinate system established by Ting et al. will be used. For this module we have expanded on the solution introduced by Ting et al. and developed the numerical approach to solve for these singularities. The solution is organized into three parts: 1) transformation of fourth order stiffness tensor, Cijkl, for arbitrary rotation in the 1-3 plane, see eqn 35 Ref.[25], 2) calcuate eigenvalues, pk, and eigen vectors, ni and tij, from equations (9) and (10) Ref.[25], and 3) calcuate exponent, k (singularity) on radius, r, from the determinate of equations (30 and (31) of Ref.[25] using appropriate boundary conditions. From these boundary conditions 12 homogeneous equations are generated. Equation coefficients are functions of elastic properties and the exponent, k on the radius vector whose origin is located at the stress free edge where the singularity exists. With the elastic properties prescribed in the NPIB form the determinate of the coefficient matrix is only a function of the exponent, k. For comparison results from Ting et al. Ref.[25] are provided below in Table 1 using the Pipes-Pagano elastic property approximation which are the default elastic properties in the NPIB form.

E1=E2=+2.100E+06,    E1=+2.000E+07
G12=G13=G23=+0.850E+06
n21= n31= n32=+0.210E+00

Table 1. Stress singularity at the free edge interface of an angle-ply graphite/epoxy laminate
Exponent
of r-k
q'=0 q'=90 q'=-q
q=0 - 3.3388 x 10-2 -
q=15 1.3528 x 10-4 3.2814 x 10-2 6.4322 x 10-4
q=30 2.6286 x 10-3 2.8682 x 10-2 1.1658 x 10-2
q=45 9.6461 x 10-3 2.0575 x 10-2 2.5575 x 10-2
q=60 1.9866 x 10-2 1.0519 x 10-2 2.3346 x 10-2
q=75 2.9388 x 10-2 2.6785 x 10-3 8.9444 x 10-3
q=90 3.3388 x 10-2 - -