Hi! Great application, it works great for estimating the effective properties of a reinforced material.

I wanted to ask of it would be possible to estimate the micro strains and micro stresses in the matrix and inclusion, after going through the code I don’t see a mention of that.

This can be done through
S_matrix=B_matrix x S_global
where
B_matrix=C_matrix(I + v(Eshelby- I)Localization_tensor)^-1 x C_matrix^-1

but this is only true for a uni directional case. As soon as we add an orientation tensor to this where a1!=a2 or a1!=a3 this does not work anymore as the Eshelby tensor and hence the localization tensor are both defined in the local coordinate system and the stress, Sigma is defined in the global system.

Maybe I’m being stupid about this, but what should be rotated and in which way? Maybe this is already solved and I have just overseen it in the code.

Hi, I came across your code when trying to find something I could test my own Mori Tanaka calculation code against. My code still doesn't seem to give the right results but meanwhile I can use your code :). However, when l looked at your calculation, it seemed the be a bit different than in the literature

According to eq. (14) in Benveniste (1987, https://doi.org/10.1016/0167-6636(87)90005-6), the composite stiffness tensor is given by

However, in your code also the term L_1 (or C0) is multiplied by the [ ... ]^-1 factor. In practice there seems to be hardly any difference in the results but I'd be interested to know if there is a reason for the different form?