Stress with temperature expansion - Code Answer

Question 1

What line comes after ‘strain = [0.0015; 0.0011; 0.0003];’?

Answer 1

beta = [alpha * (tempr - temp0); alpha * (tempr - temp0; 0];

Question 2

What are the last 3 lines required to form the D-matrix? (Just look at the D-matrix given it makes sense)

Answer 2

dmatx(1, 2) = const * poiss;
dmatx(2, 1) = const * poiss;
dmatx(3, 3) = (1.0 - poiss) * const / 2.0;

Question 3

What ‘for loops’ cone after ‘% Compute stress vector (stress)’?

Answer 3

for i = 1 : 3
     stres(i) = 0.0;
     for = j : 3;
            stres(i) = stres(i) + dmatx(i, j) * (strain(j) - beta(j));
     end
end