Circular Hole Elasticity Solution
Open Hole With Far-Field Inplane Load
In[1]:=
Define a quasi-isotropic laminate.
The maximum stress at x=0, y=1, for a hole with radius of 1, and an applied far-field load of
= 1000,
= 0,
= 0.
The theoretical stress (or load) concentration factor for an isotropic laminate should be 3.
Here is a plot of x load as a function of distance from the edge of hole.
Polar Plot for an Anisotropic laminate
Next, consider an anisotropic laminate. Because the 45's are not balanced, the laminate is not orthotropic.
Apply a unit axial load.
Now plot the polar stress distribution around the hole.Note: The stress plot is skewed because of the anisotropy
Carpet Plot of Theoretical Stress Concentration Factor
The theoretical stress concentration factor
for a hole in a composite plate is a function of the percentage of layers in each direction. For the 0, 90, ±45 family of laminates, a convenient way to show all possible values of
is a carpet plot. A carpet plotting function is available in the Material Minds system. First, we define the layer elastic constants.
Carpet plot function requires a user defined function that takes the number of 0, 90 and ±45 layers, respectively, and returns a single value. The following function executes elasticHole for the above layerProperties, and the given number of layers. elasticHole is directly wrapped by the MMLoad function which evaluates the load at the location {0,1}, and takes the first value of the load vector. Because the input load is unity, the output of this new function is
The new user function can be passed to carpetPlot
Open Hole With Far-Field Bending Load
Define a quasi-isotropic laminate.
Solve for a 1 in. radius hole, with an applied far-field moment
of unity
Note that the moment concentration factor is not equal to the load concentration factor of 3.
All of the functions can compute stress anywhere in the interior of the plate.
Here is a plot of
as a function of distance from the edge of hole.
Distribution of θ moment around hole.
Bearing Forces
As discussed in the introduction, the package implements a common approximate technique for representing a bolt load as a known radial pressure distribution.
Consider an orthotropic laminate.
The following gives the stress at {x,y} = {0,1} for a 1 inch radius hole with a 1000 pound bolt load applied at an angle of 0 (default value of BoltLoadDirection).
Take a look at the peak tangential load caused by the bearing force
The average bearing load is
= 500 lb/in. As a rough rule, the maximum tangential load is approximately the same as the average bearing load. In this case 391 versus 500.
The follow shows the radial stress distribution, which is equal to the applied pressure.
The radial stress using the holePolarPlot function.
And the tangential stress distribution.
Failure Prediction
Consider an orthotropic laminate.
0.5 in diameter hole, with a 1000 lb/in far-field load
The maximum stress failure criterion is already defined in the laminate package.
where xt, yt, xc, yc, and s are the layer strengths for longitudinal tension, transverse tension, longitudinal compression, transverse compression, and shear, respectively. The second set of arguments [{sx, sy, sxy}] are the layer stresses in layer coordinates. This set of arguments will be supplied by the failure function.
The factor-of-safety implies that the loads would be proportionally scaled by the factor to be exactly at failure.
Created by Mathematica (March 7, 2004)