class representing the elastic stiffness tensor
The stiffness tensor C is the rank 4 tensor relating stress and strain by sigma = C : epsilon. It is usually given as a 6x6 Voigt matrix, in GPa and without the factor of two of the double convention. The constructor warns if the matrix is not symmetric or not positive definite.
Syntax
C = stiffnessTensor(M,cs)
C = stiffnessTensor(M,cs,'unit','GPa')Input
| M | 6x6 Voigt matrix or 3x3x3x3 array |
| cs | crystal symmetry |
Output
| C | stiffnessTensor |
Options
| unit | physical unit of the entries, GPa by default |
Class Properties
| M | the tensor coefficients |
| rank | always 4 |
| CS | symmetry the coefficients refer to |
Example
cs = crystalSymmetry('m-3m');
C = stiffnessTensor(diag([230 230 230 117 117 117]),cs)C = stiffnessTensor (crystal)
unit: GPa
rank: 4 (3 x 3 x 3 x 3)
tensor in Voigt matrix representation:
230 0 0 0 0 0
0 230 0 0 0 0
0 0 230 0 0 0
0 0 0 117 0 0
0 0 0 0 117 0
0 0 0 0 0 117