write a general program which computes the constitutive matrix,

(1 ) write a general program which computes the constitutive matrix,
Aij= (A11,A12, A16, A22, A26, A66)
Bij= (B11, B12, B16, B22, B26, B66)
Dij= (D11, D12, D16, D22, D26, D66)

Base on input parameter,
Given N= total number of layers
Given fiber orientations (theta1,theta2, theta3, …,thetaN-1, ThetaN)
Given E1/E2= 10,
given poisson12 = 0.25 Therefore poisson21 = poison 12 (E2/E1)

(1a) we all know for each layer, we have
Q matrix for the kth layer = (Q11,Q12, Q16, Q22,Q26,Q66)
We define Q non-dim = (1/(E2)) (Q11,Q12,Q16,Q22,Q26,Q66)

For example, use equation 2.66, Q11 /E2 = (E1/E2)/(1-poisson12*poisson 21)
Computer Q non-dim
The above Q is assume we are referring to theta= 0

We now specialize the above to arbitrary value of the fiber angle, see equation 2.85, page 77
Q11bar (1/E2) = (1/E2) [ Q11 cos4 + 2(Q12+2Q66) sin2cos2 + Q22 sin4]

A11 non-dim = A11 /(E2t) A12 non-dim = A12 / (E2t) A16 non-dim= A16/ (E2t)
A22 non-dim = A22 /(E2t) A26 non-dim = A26 / (E2t) A66 non-dim= A66/ (E2t)

Then use equation 4.24 , page 198, to compute
A11/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q11bar(function of theta) . (zk-z k-1)]
A12/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q12bar(function of theta) . (zk-z k-1)]

A16/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q16bar(function of theta) . (zk-z k-1)]

A22/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q22bar(function of theta) . (zk-z k-1)]

A26/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q26bar(function of theta) . (zk-z k-1)]

A66/(E2t) = (1/E2t) [ summation from k=1 to k=N for Q66bar(function of theta) . (zk-z k-1)]

Remember t is total thickness of the entire laminate
heis the thickness of each individual layer, therefore, (N) ( he) = t

Similarly, compute Bij non-dim, and computer Dij non-dim where

B11 non-dim = B11 /(E2t2) B12 non-dim = B12 / (E2t2) B16 non-dim= B16/ (E2t2)
B22 non-dim = B22 /(E2t2) B26 non-dim = B26 / (E2t2) B66 non-dim= B66/ (E2t2)

D11 non-dim = D11 /(E2t3) D12 non-dim = B12 / (E2t3) D16 non-dim= D16/ (E2t3)
D22 non-dim = D22 /(E2t3) D26 non-dim = D26 / (E2t3) D66 non-dim= D66/ (E2t3)

2. do the specific example of layups

2a (0,0,0,0,90,90,90,90,0,0,0,0,90,90,90,90)
Question, is the above the same as (0,90,0,90) if so, why
(90,90,90,90,0,0,0,0,90,90,90,90,0,0,0,0)
Question, is the above the same as (90,0,90,0) if so, why

2b (0,45,90,-45,0,45,90,-45)symm this means total of 16 layers

2c for case of total of only four layers, computer by hand, and also by computer
And make sure your hand calculations agrees with your computer
calculations for Aijnon, Bijnon, Dij non

(22, -17,-43, -36) total of four layers
(32,-45,-23,11) total of four layers
(0,0,0,0) total of four layer

2d do the special case of isotropic orthotropic layers
(0,0,0,0) but using E1/E2= 1.0

3. read chapter 16,William D. Callister book “Materials Science and Engineering”
Do problems, 16.2, 16.6, 16.9, 16.10