5/27/2018 Matlab Code

1/3

% definitions.......................................% Vf = volume of the fiber% Vm = volume of the matrix% segmaf = fiber stress% segmam = matrix stress% segma1 = lamina stress in longitudinal direction% segma2 = matrix stress in transverse direction% epslonf = fiber strain% epslonm = matrix strain% epslon1 = lamina strain in longitudinal direction% epslon2 = lamina strain in transverse direction% Ef = fiber young modulus% Em = matrix young modulus% E1 = lamina young modulus direction1% E2 = lamina young modulus direction2% alphaf = thermal expansion coefficient of fiber

% alpham = thermal expansion coefficient of matrix% alpha1 = thermal expansion coefficient of lamina in direction1% alpha2 = thermal expansion coefficient of lamina in direction2% etaf = fiber poisson ratio

% etam = matrix poisson ratio% eta12 = inplane poisson ratio% Gf = fiber shear modulus% Gm = matrix shear modulus% G12 = inplane shear modulusclear;clc;close all;Vm = 0.45;Vf = 1-Vm;

%fiber.......................segmaf = 12;epslonf = 27;Ef = 72.345;

etaf = 0.22;Gf = 3.45;alphaf = 5.4;

%matrix......................segmam = 12;epslonm = 27;Em = 3.4;etam = 0.38;Gm = 75.9;alpham = 30;

% in direction (1)...........E1 = (Vf*Ef)+(Vm*Em)segma1 = (Vf*segmaf)+(Vm*segmam)epslon1 = segma1/E1alpha1 = ((Vf*Ef*alphaf)+(Vm*Em*alpham))/E1

% in direction (2)...........E2 = (Em*Ef)/((Vf*Em)+(Vm*Ef))epslon2 = ((segmaf*Vf)/Ef)+((segmam*Vm)/Em)

5/27/2018 Matlab Code

2/3

segma2 = E2*epslon2alpha2 = (Vf*alphaf)+(Vm*alpham)

eta12 =(etam*Vm)+(etaf*Vf)G12 = (Gm*Gf)/((Vf*Gm)+(Vm*Gf))

:utputsO

5/27/2018 Matlab Code

3/3