Elastic #

The elastic calculation can give you

1. effective elastic stiffness of the system
2. stress distribution
3. strain distribution

Solved Equation #

$\nabla \cdot \sigma = 0$, where $\sigma = C \epsilon T$

• $\sigma$ The stress tensor
• C The elastic stiffness tensor
• $\epsilon$ The strain tensor

Input options #

The following options are needed for the elastic calculation.

• input->system The choice of calculation system must be elastic
• input->external->elastic/type Either stress or strain, this will determine what external constraints are expected.
• input->external->elastic/stress External applied stress constraint, if you want to calculate the stress/strain distribution within the composite
  • tensor11$\sigma_{11}$
  • tensor12$\sigma_{12}$
  • tensor13$\sigma_{13}$
  • tensor22$\sigma_{22}$
  • tensor23$\sigma_{23}$
  • tensor33$\sigma_{33}$
• input->external->elastic/strain External applied strain constraint, if you want to calculate the stress/strain distribution within the composite
  • tensor11$\epsilon_{11}$
  • tensor12$\epsilon_{12}$
  • tensor13$\epsilon_{13}$
  • tensor22$\epsilon_{22}$
  • tensor23$\epsilon_{23}$
  • tensor33$\epsilon_{33}$
• input->solver->ref A property tensor for the reference phase called stiffness is needed. Tensor definition is explained [here]

<system>
<type>elastic</type>
    <external>
        <elastic>
            <type>strain</type> 
            <strain>
                <tensor11>0.0</tensor11>
                <tensor22>0.0</tensor22>
                <tensor33>0.0</tensor33>
                <tensor23>0.0</tensor23>
                <tensor13>0.001</tensor13>
                <tensor12>0.0</tensor12>
            </strain>
        </elastic>
    </external>
    <solver>
        <ref>
            <tensor>
                <pointGroup>Oh</pointGroup>
                <name>stiffness</name>
                <rank>4</rank>
                <component>
                    <value>1000</value>
                    <index>1111</index>
                </component>
                <component>
                    <value>1200</value>
                    <index>1122</index>
                </component>
                <component>
                    <value>1100</value>
                    <index>1212</index>
                </component>
                <component>
                    <value>1100</value>
                    <index>1221</index>
                </component>
            </tensor>
        </ref>
    </solver>
</system>

Output files #

There are a total of 4 possible output files for the dielectric calculation

1. effElasticStiffness.dat
2. avElasticVariables.dat
3. stress.00000000.dat
4. strain.00000000.dat