3.2 MATLAB Functions Used 吧 if p==1 y =1/(Vf/G12f Vm/Gm); elseif p ==2 y =1/((Vf/G12f EtaPrime*Vm/Gm)/(Vf EtaPrime*Vm)); elseif p ==3 y Gm*((Gm +G12f)-Vf*(Gm -G12f))/((Gm G12f)+ Vf*(Gm G12f)) end function y Alpha1(Vf,Eif,Em,Alphaif,Alpham) lpha1 This function returns the coefficient of thermal % expansion in the longitudinal direction. Its input are five values: Vf -fiber volume fraction 名 E1f longitudinal Young's modulus of the fiber Em - Young's modulus of the matrix % Alphalf -coefficient of thermal expansion in the % 1-direction for the fiber % Alpham coefficient of thermal expansion for the matrix Vm 1-Vf; y (Vf+E1f+Alpha1f Vm*Em*Alpham)/(E1f*Vf Em*Vm); function y Alpha2(Vf,Alpha2f,Alpham,E1,E1f,Em,NU1f,NUm, Alphaif,p) %Alpha2 This function returns the coefficient of thermal % expansion in the transverse direction. % Its input are ten values: Vf fiber volume fraction % Alpha2f coefficient of thermal expansion in the % 2-direction for the fiber % Alpham -coefficient of thermal expansion for the matrix % E1 -longitudinal Young's modulus of the lamina E1f -longitudianl Young's modulus of the fiber 名 Em -Young's modulus of the matrix % NU1f -Poisson's ratio of the fiber NUm -Poisson's ratio of the matrix % Alphaif -coefficient of thermal expansion in the 名 1-direction p parameter used to determine which equation to use % p =1 -use equation (3.8) p =2 -use equation (3.7) Use the value zero for any argument not needed in % the calculation Vm 1 Vf; 1fp=1 y Vf*Alpha2f Vm*Alpham; elseif p ==2 y =(Alpha2f -(Em/E1)*NU1f*(Alpham -Alphalf)+Vm)*Vf+ (Alpham (Eif/E1)*NUm*(Alpham Alphaif)*Vf)*Vm; end3.2 MATLAB Functions Used 31 if p == 1 y = 1/(Vf/G12f + Vm/Gm); elseif p == 2 y = 1/((Vf/G12f + EtaPrime*Vm/Gm)/(Vf + EtaPrime*Vm)); elseif p == 3 y = Gm*((Gm + G12f) - Vf*(Gm - G12f))/((Gm + G12f) + Vf*(Gm - G12f)); end function y = Alpha1(Vf,E1f,Em,Alpha1f,Alpham) %Alpha1 This function returns the coefficient of thermal % expansion in the longitudinal direction. % Its input are five values: % Vf - fiber volume fraction % E1f - longitudinal Young’s modulus of the fiber % Em - Young’s modulus of the matrix % Alpha1f - coefficient of thermal expansion in the % 1-direction for the fiber % Alpham - coefficient of thermal expansion for the matrix Vm = 1 - Vf; y = (Vf*E1f*Alpha1f + Vm*Em*Alpham)/(E1f*Vf + Em*Vm); function y = Alpha2(Vf,Alpha2f,Alpham,E1,E1f,Em,NU1f,NUm, Alpha1f,p) %Alpha2 This function returns the coefficient of thermal % expansion in the transverse direction. % Its input are ten values: % Vf - fiber volume fraction % Alpha2f - coefficient of thermal expansion in the % 2-direction for the fiber % Alpham - coefficient of thermal expansion for the matrix % E1 - longitudinal Young’s modulus of the lamina % E1f - longitudianl Young’s modulus of the fiber % Em - Young’s modulus of the matrix % NU1f - Poisson’s ratio of the fiber % NUm - Poisson’s ratio of the matrix % Alpha1f - coefficient of thermal expansion in the % 1-direction % p - parameter used to determine which equation to use % p = 1 - use equation (3.8) % p = 2 - use equation (3.7) % Use the value zero for any argument not needed in % the calculation Vm = 1 - Vf; if p == 1 y = Vf*Alpha2f + Vm*Alpham; elseif p == 2 y = (Alpha2f - (Em/E1)*NU1f*(Alpham - Alpha1f)*Vm)*Vf + (Alpham + (E1f/E1)*NUm*(Alpham - Alpha1f)*Vf)*Vm; end