0% found this document useful (0 votes)
141 views5 pages

Laminate Strength Analysis with MATLAB

This document analyzes the strength of two laminate configurations - (0@2t/25@t) and (0@t/25@1.5t) - under uniaxial loading using MATLAB. It calculates shear strain, applies maximum stress and first ply failure theories, and determines that configuration (0@2t/25@t) has higher shear strain and lower first ply failure load than configuration (0@t/25@1.5t).

Uploaded by

Vijay Kishore
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
141 views5 pages

Laminate Strength Analysis with MATLAB

This document analyzes the strength of two laminate configurations - (0@2t/25@t) and (0@t/25@1.5t) - under uniaxial loading using MATLAB. It calculates shear strain, applies maximum stress and first ply failure theories, and determines that configuration (0@2t/25@t) has higher shear strain and lower first ply failure load than configuration (0@t/25@1.5t).

Uploaded by

Vijay Kishore
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

ASSIGNMENT 2:

LAMINATE STRENGTH ANALYSIS

Submitted by
AE12M004- DHARAMENDR KUMAR BHARADWAJ
AE12M005- JADHAV KRISHNA BALU
AE12M019- S VIJAY KISHORE
AE12M022- D VISHNURAM

MATLAB CODE FOR DETERMINATION OF SHEAR STRAIN(


and
ACCORDING TO MAXIMUM STRESS THEORY AND FIRST PLY FAILURE CRITERIA
clear all;
clc;
format long eng;
% MATERIAL STRENGTH PARAMETERS in Pa
XT=1.364*10^9;XC=1.364*10^9;YT=0.0448*10^9;YC=0.0448*10^9;
% MATERIAL PROPERTIES IN Pa
EL=200*10^9;ET=20*10^9;GL=10*10^9;NU=0.3;PU=0.03;
Q=[(EL/(1-PU*NU)) EL*PU/(1-PU*NU) 0;ET*NU/(1-PU*NU) ET/(1-PU*NU) 0;0 0 GL];
% PLY POSITIONS AND ORIENTATIONS
theta=[0 25 25 0];
c=input('enter 1 for (0@2t/25@t) or enter 2 for (0@t/25@1.5t)= ');
switch(c)
case 1
% For 0 deg at 2t and 25 deg ply at t
t=1*10^(-3);
h0=-t/2;
h1=h0+2*t/6;
h2=h1+t/6;
h3=h2+t/6;
h4=h3+2*t/6;
otherwise
t=1*10^(-3);
% For 0 deg ply at t and 25 deg ply at 1.5t
h0=-t/2;
h1=h0+t/5;
h2=h1+1.5*t/5;
h3=h2+1.5*t/5;
h4=h3+t/5;
end
h=[h0 h1 h2 h3 h4];
% REDUCED STIFFNESS MATRIX QBAR=QQ
for PLY=1:4
c=cosd(theta(PLY));s=sind(theta(PLY));
% TRANSFORMATION MATRIX
T(:,:,PLY)=[c.^2 s.^2 2.*s.*c;
s.^2 c.^2 -2.*s.*c;
-c.*s c.*s c.^2-s.^2];
R=[1 0 0;0 1 0;0 0 2];
QQ(:,:,PLY)=T(:,:,PLY)^-1*Q*R*T(:,:,PLY)*R^-1;
End

% COMPUTATION OF EXTENSIONAL STIFFNESS MATRIX


A=zeros(3);
for i=1:4
A=A+QQ(:,:,i)*(h(i+1)-h(i));
end
NX=1;NY=0;NXY=1;
N=[NX NY NXY];
% MID-PLANE STRAINS
MS=A^-1*N';
disp('GAMMAXY_N in mm/N is');
GAMMAXY_N=MS(3)*1000
% STRESS DETERMINATION IN GLOBAL REFERENCE AXES SYSTEM
for PLY=1:4
SIGMA(:,:,PLY)=QQ(:,:,PLY)*MS;
end
% STRESSS DETERMINATION IN MATERIAL AXES SYSYTEM
for PLY=1:4
MSIGMA(:,:,PLY)=T(:,:,PLY)*SIGMA(:,:,PLY);
End
% APPLYING FAILURE THEORY
for PLY=1:2
j=1;
i=1;
if MSIGMA(i,j,PLY)<0
SR1(:,:,PLY)=abs(MSIGMA(i,j,PLY))/XC;
else
SR1(:,:,PLY)=abs(MSIGMA(i,j,PLY))/XT;
end
if MSIGMA(i+1,j,PLY)<0
SR2(:,:,PLY)=abs(MSIGMA(i+1,j,PLY))/YC;
else
SR2(:,:,PLY)=abs(MSIGMA(i+1,j,PLY))/YT;
end
end
i=1;j=1;
SR=[SR1(i,j,1) SR2(i,j,1) SR1(i,j,2) SR2(i,j,2)];
FP=1/max(SR);
for i=1:4
if max(SR)==SR(i)
location=i;
switch(location)
case 1
disp('Failure occurs in longitudinal direction of 0 deg ply')

case 2
disp('failure occurs in transverse direction of 0 deg ply')
case 3
disp('Failure occurs in longitudinal direction of 25 deg ply')
otherwise
disp('Failure occurs in transverse direction of 25 deg ply')
end
end
end
disp('The first ply failure load in N/mm is');
fprintf('%f\n,%f',FP/1000);

PROBLEM
Case1: (0@2t/25@t)s
Case2: (0@t/25@1.5t)s
Loading:
Total laminate thickness = 1mm
Material strength Parameters in Pa
=1.364*10^9; =1.364*10^9; =0.0448*10^9; ==0.0448*10^9;
Material Properties in Pa
=200*10^9 Pa, =20*10^9 Pa,
=10*10^9Pa , =0.3
RESULTS
CONFIGURATION
(0@2t/25@t)s
(0@t/25@1.5t)s
The

(mm/N)
5.7024e-05
4.7051e-05

(N/mm)
86.9222
85.4407

and shear strain per unit load are higher for (0@2t/25@t)s configuration.

Common questions

Powered by AI

Ply orientations affect the distribution of strains through the laminate thickness, impacting shear strain per unit load. The configuration (0@2t/25@t)s exhibits higher shear strain per unit load compared to (0@t/25@1.5t)s, indicating that an increased thickness of plies oriented at zero degrees influences shear deformation .

Shear strain (GAMMAXY_N) is calculated by determining mid-plane strains from the inverse extensional stiffness matrix applied to load N, and converting the third component to millimeters per Newton . Shear strain is important as it provides insights into the deformation of the laminate under in-plane shear loading, influencing the laminate's design and potential failure points.

The transformation matrix is used to convert stresses and strains from the global reference axis to the material axes of each ply, accommodating differences in orientation between plies . This conversion is necessary because the material properties and failure criteria are defined within the material axes, requiring accurate comparison against material limits to predict failure accurately.

Material properties like moduli of elasticity and strength define how a laminate can withstand loads, while ply layup configurations determine the stress and strain distribution within the laminate. A suitable combination maximizes strength and delay failure by aligning high-strength materials in load-bearing orientations, optimizing laminate endurance against specified operational conditions .

Material strength parameters, such as longitudinal tensile and compressive strengths (1.364x10^9 Pa) and transverse tensile and compressive strengths (0.0448x10^9 Pa), significantly affect computational results by setting thresholds for permissible stress before failure . These parameters ensure that calculations align with realistic material behavior, thereby influencing laminate strength and performance evaluations.

The orientation of plies significantly affects the failure characteristics of the laminate. For instance, in the case of (0@2t/25@t)s configuration, the failure commonly occurs in the longitudinal direction of the 0-degree ply due to its orientation that allows greater stress in fiber direction . Conversely, different orientations lead to different stress distributions and hence different failure points.

First ply failure occurs at different loads for different configurations, such as 86.9222 N/mm for (0@2t/25@t)s and 85.4407 N/mm for (0@t/25@1.5t)s . These observations imply that ply orientation affects the initial failure point, influencing design decisions to optimize strength and delay failure. Such insights help in selecting ply layups and orientations that suit specific loading conditions.

The extensional stiffness matrix (A) is calculated by summing the contributions of each ply's reduced stiffness matrix (QBAR) across the laminate thickness . This matrix is significant because it represents the overall stiffness of the laminate when subjected to in-plane forces. It helps predict how the laminate will deform under applied loads, essential for structural analysis and design.

The material properties of the laminate are: Longitudinal modulus (EL) of 200 GPa, transverse modulus (ET) of 20 GPa, shear modulus (GL) of 10 GPa, and Poisson's ratio (NU) of 0.3 . These properties influence the laminate's performance by determining its stiffness, strength, and ability to undergo deformation under stress. The higher longitudinal modulus increases stiffness along the fiber direction, while the lower transverse modulus and shear modulus control deformation in off-axis directions.

Applying failure theories is crucial to predict the load at which a laminate will fail, facilitating design safety and reliability. The analysis employed the maximum stress theory and first ply failure criteria to predict failure points by comparing calculated stresses against tensile and compressive material strengths in different directions .

You might also like