elasticipy 3.0.0__tar.gz → 4.0.0__tar.gz

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Elasticipy_vs_pymatgen.py

@@ -1,7 +1,5 @@
-from docutils.utils.math.latex2mathml import layout_styles
-
-from Elasticipy.StressStrainTensors import StrainTensor
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.stress_strain import StrainTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 from matplotlib import pyplot as plt
 import matplotlib as mpl
 mpl.use('Qt5Agg')   # Ensure interactive plot

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Elate_vs_Elasticipy.py

@@ -1,9 +1,9 @@
 import numpy as np
 from elastic import Elastic
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 import matplotlib as mpl
 mpl.use('Qt5Agg')   # Ensure interactive plot
-from Elasticipy.SphericalFunction import sph2cart, _plot3D
+from Elasticipy.spherical_function import sph2cart, _plot3D
 import time
 import matplotlib.pyplot as plt
 

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Example_Stiffness_tensor.py

@@ -1,4 +1,4 @@
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 from scipy.spatial.transform import Rotation
 import matplotlib as mpl
 mpl.use('Qt5Agg')   # Ensure interactive plot

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Example_StressStrain_arrays.py

@@ -2,7 +2,7 @@
 # Imports and simple example of stress
 # ======================================================
 import numpy as np
-from Elasticipy.StressStrainTensors import StressTensor, StrainTensor
+from Elasticipy.tensors.stress_strain import StressTensor, StrainTensor
 
 
 stress = StressTensor([[0, 1, 0],
@@ -23,7 +23,7 @@ print(strain.volumetric_strain())
 # ======================================================
 # Linear elasticity
 # ======================================================
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 
 C = StiffnessTensor.fromCrystalSymmetry(symmetry='cubic', phase_name='ferrite',
                                         C11=274, C12=175, C44=89)

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Example_WaveVelocity.py

@@ -1,7 +1,6 @@
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 import matplotlib as mpl
 mpl.use('Qt5Agg')   # Ensure interactive plot
-from matplotlib import pyplot as plt
 
 C = StiffnessTensor.fromCrystalSymmetry(symmetry='orthorombic', phase_name='forsterite',
                                         C11=320, C12=68.2, C13=71.6,

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/Multiple_phases.py

@@ -1,4 +1,4 @@
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 import numpy as np
 from matplotlib import pyplot as plt
 

elasticipy-4.0.0/Examples/SelfConsistent.py

@@ -0,0 +1,89 @@
+import numpy as np
+from Elasticipy.tensors.fourth_order import FourthOrderTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
+from scipy.integrate import trapezoid
+from scipy.spatial.transform import Rotation
+
+I = FourthOrderTensor.identity()
+
+def gamma(C_macro_local, phi, theta, a1, a2, a3):
+    s1 = np.sin(theta)*np.cos(phi) / a1
+    s2 = np.sin(theta)*np.sin(phi) / a2
+    s3 = np.cos(theta) / a3
+    s = np.array([s1, s2, s3])
+    D = np.einsum('lmnp,pqr,lqr->qrmn', C_macro_local.full_tensor(), s, s)
+    Dinv = np.linalg.inv(D)
+    return np.einsum('qrjk,iqr,lqr->qrijkl', Dinv, s, s) # The symmetrization is made afterward (see below)
+
+def polarization_tensor(C_macro_local, a1, a2, a3, n_phi, n_theta):
+    theta = np.linspace(0, np.pi, n_theta)
+    phi = np.linspace(0, 2 * np.pi, n_phi)
+    phi_grid, theta_grid = np.meshgrid(phi, theta, indexing='ij')
+    g = gamma(C_macro_local, phi_grid, theta_grid, a1, a2, a3)
+    gsin = (g.T*np.sin(theta_grid.T)).T
+    a = trapezoid(gsin, phi, axis=0)
+    b= trapezoid(a, theta, axis=0)/(4*np.pi)
+    return FourthOrderTensor(b, force_minor_symmetry=True) # Symmetrization here (see above)
+
+def localization_tensor(C_macro_local, C_incl, n_phi, n_theta, a1, a2, a3):
+    E = polarization_tensor(C_macro_local, a1, a2, a3, n_phi, n_theta)
+    delta = FourthOrderTensor(C_incl.matrix - C_macro_local.matrix)
+    Ainv = E.ddot(delta) + I
+    return Ainv.inv()
+
+def Kroner_Eshelby(Ci, g, method='strain', max_iter=5, atol=1e-3, rtol=1e-3, display=False, n_phi=100, n_theta=100, particle_size=None):
+    Ci_rotated = (Ci * g)
+    C_macro = Ci_rotated.Hill_average()
+    eigen_stiff = C_macro.eig_stiffnesses
+    keep_on = True
+    k = 0
+    message = 'Maximum number of iterations is reached'
+    m = len(g)
+    A_local = FourthOrderTensor.zeros(m)
+    if particle_size is None:
+        a1 = a2 = a3 = 1
+    else:
+        a1, a2, a3 = particle_size
+    while keep_on:
+        eigen_stiff_old = eigen_stiff
+        C_macro_local = C_macro * (g.inv())
+        for i in range(m):
+            A_local[i] = localization_tensor(C_macro_local[i], Ci, n_phi, n_theta, a1, a2, a3)
+        A = A_local * g
+        Q = Ci_rotated.ddot(A)
+        if method=='stress':
+            CiAi_mean = Q.mean()
+            C_macro = StiffnessTensor.from_Kelvin(CiAi_mean.matrix, force_symmetries=True)
+            err = A.mean() - FourthOrderTensor.identity()
+        else:
+            B = Q.ddot(C_macro.inv())
+            R = Ci_rotated.inv().ddot(B)
+            R_mean = R.mean()
+            C_macro = StiffnessTensor.from_Kelvin(R_mean.inv().matrix, force_symmetries=True)
+            err = B.mean() - FourthOrderTensor.identity()
+
+        # Stopping criteria
+        eigen_stiff = C_macro.eig_stiffnesses
+        abs_change = np.abs(eigen_stiff - eigen_stiff_old)
+        rel_change = np.max(abs_change / eigen_stiff_old)
+        max_abs_change = np.max(abs_change)
+        k += 1
+        if  max_abs_change < atol:
+            keep_on = False
+            message = 'Absolute change is below threshold value'
+        if rel_change < rtol:
+            keep_on = False
+            message = 'Relative change is below threshold value'
+        if k == max_iter:
+            keep_on = False
+        if display:
+            err = np.max(np.abs(err.matrix))
+            print('Iter #{}: abs. change={:0.5f}; rel. change={:0.5f}; error={:0.5f}'.format(k, max_abs_change, rel_change,err))
+    return C_macro, message
+
+Cstrip = StiffnessTensor.transverse_isotropic(Ex= 10.2, Ez=146.8, nu_zx=0.274, nu_yx=0.355, Gxz=7)
+Cstrip = Cstrip * Rotation.from_euler('Y', 90, degrees=True)
+#orientations = Rotation.from_euler('Z', np.linspace(0, 180, 10, endpoint=False), degrees=True)
+orientations = Rotation.random(100, random_state=1234)
+Cstrip = StiffnessTensor.cubic(C11=186, C12=134, C44=77)
+C_stress, reason = Kroner_Eshelby(Cstrip, orientations, display=True, method='strain')

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/essai_plasticity.py

@@ -1,6 +1,6 @@
-from Elasticipy.Plasticity import JohnsonCook
-from Elasticipy.StressStrainTensors import StressTensor, StrainTensor
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.plasticity import JohnsonCook
+from Elasticipy.tensors.stress_strain import StressTensor, StrainTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 import numpy as np
 from matplotlib import pyplot as plt
 import matplotlib as mpl

{elasticipy-3.0.0 → elasticipy-4.0.0}/Examples/example_readwrite.py

@@ -1,4 +1,4 @@
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+from Elasticipy.tensors.elasticity import StiffnessTensor
 
 C = StiffnessTensor.from_MP("mp-1048")
 C.save_to_txt("TiNi.txt")

{elasticipy-3.0.0 → elasticipy-4.0.0}/JOSS/paper.md

@@ -142,7 +142,9 @@ with ease. It already implements thermal expansion.
 \autoref{fig:Young}.a) and b) were rendered with the following syntax:
 
 ````python
-from Elasticipy.FourthOrderTensor import StiffnessTensor
+
+from Elasticipy.tensors.elasticity import StiffnessTensor
+
 C = StiffnessTensor.cubic(C11=186, C12=134, C44=77)
 E = C.Young_modulus
 E.plot3D(n_phi=500, n_theta=500)
@@ -181,10 +183,11 @@ Efforts have been made to provide out-of-the-box simple syntaxes for common oper
 will create a tensor array corresponding to evenly-spaced strain along $[1,0,0]$ axis:
 
 ````python
-from Elasticipy.StressStrainTensors import StrainTensor
-m = 1000                     # length of tensor array
-mag = np.linspace(0, 0.1, m) # Strain magnitude
-strain = StrainTensor.tensile([1,0,0], mag)
+from Elasticipy.tensors.stress_strain import StrainTensor
+
+m = 1000  # length of tensor array
+mag = np.linspace(0, 0.1, m)  # Strain magnitude
+strain = StrainTensor.tensile([1, 0, 0], mag)
 ````
 
 Given the stiffness tensor ``C`` (see above), one can compute the corresponding stress array with:

{elasticipy-3.0.0/src/elasticipy.egg-info → elasticipy-4.0.0}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: elasticipy
-Version: 3.0.0
+Version: 4.0.0
 Summary: A Python library for elasticity tensor computations
 Author-email: Dorian Depriester <dorian.dep@gmail.com>
 License: MIT
@@ -27,8 +27,8 @@ Requires-Dist: pytest-cov; extra == "dev"
 Requires-Dist: pymatgen; extra == "dev"
 Requires-Dist: orix; extra == "dev"
 Requires-Dist: mp_api; extra == "dev"
+Dynamic: license-file
 
-# Elasticipy 
 [![PyPI - Version](https://img.shields.io/pypi/v/Elasticipy?link=https%3A%2F%2Fpypi.org%2Fproject%2FElasticipy%2F)](https://pypi.org/project/elasticipy/)
 [![PyPI - Downloads](https://img.shields.io/pypi/dm/Elasticipy?link=https%3A%2F%2Fpypi.org%2Fproject%2FElasticipy%2F)](https://pypistats.org/packages/elasticipy)
 [![PyPI - License](https://img.shields.io/pypi/l/Elasticipy)](https://github.com/DorianDepriester/Elasticipy/blob/main/LICENSE)
@@ -38,9 +38,9 @@ Requires-Dist: mp_api; extra == "dev"
 ![PyPI - Python Version](https://img.shields.io/pypi/pyversions/Elasticipy)
 [![status](https://joss.theoj.org/papers/8cce91b782f17f52e9ee30916cd86ad5/status.svg)](https://joss.theoj.org/papers/8cce91b782f17f52e9ee30916cd86ad5)
 
+# ![Elasticipy](docs/source/logo/logo_text.svg)
 
-
-A python toolkit to manipulate strain and strain tensors, and other linear elasticity-related tensors (e.g. stiffness). 
+A python toolkit to manipulate stress and strain tensors, and other linear elasticity-related tensors (e.g. stiffness). 
 This package also provides a collection of easy-to-use and very fast tools to work on stress and strain tensors.
 
 ## Main features

{elasticipy-3.0.0 → elasticipy-4.0.0}/README.md

@@ -1,4 +1,3 @@
-# Elasticipy 
 [![PyPI - Version](https://img.shields.io/pypi/v/Elasticipy?link=https%3A%2F%2Fpypi.org%2Fproject%2FElasticipy%2F)](https://pypi.org/project/elasticipy/)
 [![PyPI - Downloads](https://img.shields.io/pypi/dm/Elasticipy?link=https%3A%2F%2Fpypi.org%2Fproject%2FElasticipy%2F)](https://pypistats.org/packages/elasticipy)
 [![PyPI - License](https://img.shields.io/pypi/l/Elasticipy)](https://github.com/DorianDepriester/Elasticipy/blob/main/LICENSE)
@@ -8,9 +7,9 @@
 ![PyPI - Python Version](https://img.shields.io/pypi/pyversions/Elasticipy)
 [![status](https://joss.theoj.org/papers/8cce91b782f17f52e9ee30916cd86ad5/status.svg)](https://joss.theoj.org/papers/8cce91b782f17f52e9ee30916cd86ad5)
 
+# ![Elasticipy](docs/source/logo/logo_text.svg)
 
-
-A python toolkit to manipulate strain and strain tensors, and other linear elasticity-related tensors (e.g. stiffness). 
+A python toolkit to manipulate stress and strain tensors, and other linear elasticity-related tensors (e.g. stiffness). 
 This package also provides a collection of easy-to-use and very fast tools to work on stress and strain tensors.
 
 ## Main features

elasticipy-4.0.0/docs/source/Elasticipy.plasticity.rst

@@ -0,0 +1,7 @@
+Elasticipy.plasticity module
+============================
+
+.. automodule:: Elasticipy.plasticity
+   :members:
+   :undoc-members:
+   :show-inheritance:

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Elasticipy.rst

@@ -7,11 +7,13 @@ Submodules
 .. toctree::
    :maxdepth: 4
 
-   Elasticipy.FourthOrderTensor
-   Elasticipy.SecondOrderTensor
-   Elasticipy.SphericalFunction
-   Elasticipy.StressStrainTensors
-   Elasticipy.Plasticity
+   Elasticipy.tensors.fourth_order
+   Elasticipy.tensors.elasticity
+   Elasticipy.tensors.second_order
+   Elasticipy.tensors.stress_strain
+   Elasticipy.spherical_function
+   Elasticipy.plasticity
+
 
 Module contents
 ---------------

elasticipy-4.0.0/docs/source/Elasticipy.spherical_function.rst

@@ -0,0 +1,7 @@
+Elasticipy.spherical_function module
+====================================
+
+.. automodule:: Elasticipy.spherical_function
+   :members:
+   :undoc-members:
+   :show-inheritance:

elasticipy-4.0.0/docs/source/Elasticipy.tensors.elasticity.rst

@@ -0,0 +1,7 @@
+Elasticipy.tensors.elasticity module
+====================================
+
+.. automodule:: Elasticipy.tensors.elasticity
+   :members:
+   :undoc-members:
+   :show-inheritance:

elasticipy-4.0.0/docs/source/Elasticipy.tensors.fourth_order.rst

@@ -0,0 +1,7 @@
+Elasticipy.tensors.fourth_order module
+======================================
+
+.. automodule:: Elasticipy.tensors.fourth_order
+   :members:
+   :undoc-members:
+   :show-inheritance:

elasticipy-4.0.0/docs/source/Elasticipy.tensors.second_order.rst

@@ -0,0 +1,7 @@
+Elasticipy.tensors.second_order module
+======================================
+
+.. automodule:: Elasticipy.tensors.second_order
+   :members:
+   :undoc-members:
+   :show-inheritance:

elasticipy-4.0.0/docs/source/Elasticipy.tensors.stress_strain.rst

@@ -0,0 +1,7 @@
+Elasticipy.tensors.stress_strain module
+=======================================
+
+.. automodule:: Elasticipy.tensors.stress_strain
+   :members:
+   :undoc-members:
+   :show-inheritance:

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Plasticity.rst

@@ -39,17 +39,44 @@ First, let us create the model:
 
 .. doctest::
 
-    >>> from Elasticipy.Plasticity import JohnsonCook
+    >>> from Elasticipy.plasticity import JohnsonCook
     >>> JC = JohnsonCook(A=363, B=792.7122, n=0.5756)
 
 The parameters are taken from [1]_. As we will also take elastic behaviour into account, we also need:
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.isotropic(E=210000, nu=0.27)
 
 Now, let say that we want to investigate the material's response for the tensile stress ranging from 0 to 725 MPa:
 
-    >>> from Elasticipy.StressStrainTensors import StressTensor, StrainTensor
+    >>> from Elasticipy.tensors.stress_strain import StressTensor, StrainTensor
+    >>> import numpy as np
+    >>> n_step = 100
+    >>> stress_mag = np.linspace(0, 725, n_step)
+    >>> stress = StressTensor.tensile([1,0,0], stress_mag)
+
+At least, we can directly compute the elastic strain for each step:
+
+    >>> elastic_strain = C.inv() * stress
+
+So now, the plastic strain can be computed using an iterative approach:
+
+    >>> plastic_strain = StrainTensor.zeros(n_step)
+    >>> for i in range(1, n_step):
+    ...       strain_increment = JC.compute_strain_increment(stress[i])
+    ...       plastic_strain[i] = plastic_strain[i-1] + strain_increment
+
+That's all. Finally, let us plot the applied stress as a function of the overall elongation:
+
+    >>> from matplotlib import pyplot as plt
+    >>> elong = elastic_strain.C[0,0]+plastic_strain.C[0,0]
+    >>> fig, ax = plt.subplots()
+    >>> ax.plot(elong, stress_mag, label='Stress-controlled') # doctest: +SKIP
+    >>> ax.set_xlabel(r'$\varepsilon_{xx}$') # doctest: +SKIP
+    >>> ax.set_ylabel('Tensile stress (MPa)') # doctest: +SKIP
+
+Now, let say that we want to investigate the material's response for the tensile stress ranging from 0 to 725 MPa:
+
     >>> import numpy as np
     >>> n_step = 100
     >>> stress_mag = np.linspace(0, 725, n_step)

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_AveragingMethods.rst

@@ -17,7 +17,7 @@ As an example, let's consider the stiffness tensor for monoclinic TiNi:
 
 .. doctest::
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.monoclinic(phase_name='TiNi',
     ...                                C11=231, C12=127, C13=104,
     ...                                C22=240, C23=131, C33=175,
@@ -36,7 +36,7 @@ behaviour. In this case, the aforementioned averages can be estimated as follows
     Let's see how ``C_Hill`` looks like:
 
     >>> print(C_Hill)
-    Stiffness tensor (in Voigt notation):
+    Stiffness tensor (in Voigt mapping):
     [[200.92610744 119.59583533 119.59583533   0.           0.
         0.        ]
      [119.59583533 200.92610744 119.59583533   0.           0.

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_MultiplePhases.rst

@@ -42,7 +42,7 @@ Let's start by creating the two stiffness tensors in Python:
 
 .. doctest::
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C_austenite = StiffnessTensor.cubic(C11=204, C12=137, C44=126)
     >>> C_ferrite = StiffnessTensor.cubic(C11=242, C12=146, C44=116)
 

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_ReadWriteFiles.rst

@@ -12,7 +12,7 @@ To save a given stiffness tensor to a plain text file, just use the ``save_to_tx
 
 .. doctest::
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.cubic(C11=186, C12=134, C44=77, phase_name='Cu')
     >>> C.save_to_txt('Stiffness_Cu.txt')
 
@@ -25,13 +25,14 @@ To read entries from an existing text file, just run ``StiffnessTensor.from_txt_
 Using the example above:
 
     >>> StiffnessTensor.from_txt_file('Stiffness_Cu.txt')
-    Stiffness tensor (in Voigt notation) for Cu:
+    Stiffness tensor (in Voigt mapping):
     [[186. 134. 134.   0.   0.   0.]
      [134. 186. 134.   0.   0.   0.]
      [134. 134. 186.   0.   0.   0.]
      [  0.   0.   0.  77.   0.   0.]
      [  0.   0.   0.   0.  77.   0.]
      [  0.   0.   0.   0.   0.  77.]]
+    Phase: Cu
     Symmetry: cubic
 
 

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_StiffnessTensor.rst

@@ -12,7 +12,7 @@ First, create a stiffness tensor with a given symmetry (let say, monoclinic):
 
 .. doctest::
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.monoclinic(phase_name='TiNi',
     ...                                C11=231, C12=127, C13=104,
     ...                                C22=240, C23=131, C33=175,

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_StressStrain.rst

@@ -11,7 +11,7 @@ equivalent stresses:
 
 .. doctest::
 
-    >>> from Elasticipy.StressStrainTensors import StressTensor, StrainTensor
+    >>> from Elasticipy.tensors.stress_strain import StressTensor, StrainTensor
             >>> stress = StressTensor.shear([1, 0, 0], [0, 1, 0], 1.0) # Unit XY shear stress
             >>> print(stress.vonMises(), stress.Tresca())
             1.7320508075688772 2.0
@@ -32,10 +32,10 @@ Linear elasticity
 This section is dedicated to linear elasticity, hence introducing the fourth-order stiffness tensor.
 As an example, create a stiffness tensor corresponding to steel:
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.isotropic(E=210e3, nu=0.28)
     >>> print(C)
-    Stiffness tensor (in Voigt notation):
+    Stiffness tensor (in Voigt mapping):
     [[268465.90909091 104403.40909091 104403.40909091      0.
            0.              0.        ]
      [104403.40909091 268465.90909091 104403.40909091      0.
@@ -64,7 +64,7 @@ Conversely, one can compute the compliance tensor:
 
     >>> S = C.inv()
     >>> print(S)
-    Compliance tensor (in Voigt notation):
+    Compliance tensor (in Voigt mapping):
     [[ 4.76190476e-06 -1.33333333e-06 -1.33333333e-06  0.00000000e+00
        0.00000000e+00  0.00000000e+00]
      [-1.33333333e-06  4.76190476e-06 -1.33333333e-06  0.00000000e+00
@@ -97,13 +97,14 @@ As an illustration, we consider the anisotropic behaviour of ferrite:
     >>> C = StiffnessTensor.fromCrystalSymmetry(symmetry='cubic', phase_name='ferrite',
     ...                                         C11=274, C12=175, C44=89)
     >>> print(C)
-    Stiffness tensor (in Voigt notation) for ferrite:
+    Stiffness tensor (in Voigt mapping):
     [[274. 175. 175.   0.   0.   0.]
      [175. 274. 175.   0.   0.   0.]
      [175. 175. 274.   0.   0.   0.]
      [  0.   0.   0.  89.   0.   0.]
      [  0.   0.   0.   0.  89.   0.]
      [  0.   0.   0.   0.   0.  89.]]
+    Phase: ferrite
     Symmetry: cubic
 
 Let's start by creating an array of 10 stresses:

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_ThermalExpansion.rst

@@ -8,7 +8,7 @@ At first, let's consider steel as an example, assuming an isotropic behaviour:
 
 .. doctest::
 
-    >>> from Elasticipy.ThermalExpansion import ThermalExpansionTensor as ThEx
+    >>> from Elasticipy.tensors.thermal_expansion import ThermalExpansionTensor as ThEx
     >>> alpha = ThEx.isotropic(11e-6)
     >>> print(alpha)
     Thermal expansion tensor

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/Tutorials/Tutorial_wave-velocities.rst

@@ -11,18 +11,19 @@ We thus start by defining the stiffness tensor for forsterite:
 
 .. doctest::
 
-    >>> from Elasticipy.FourthOrderTensor import StiffnessTensor
+    >>> from Elasticipy.tensors.elasticity import StiffnessTensor
     >>> C = StiffnessTensor.orthorhombic(phase_name='forsterite',
     ...                                  C11=320, C12=68.2, C13=71.6, C22=196.5, C23=76.8,
     ...                                  C33=233.5, C44=64, C55=77, C66=78.7)
     >>> print(C)
-    Stiffness tensor (in Voigt notation) for forsterite:
+    Stiffness tensor (in Voigt mapping):
     [[320.   68.2  71.6   0.    0.    0. ]
      [ 68.2 196.5  76.8   0.    0.    0. ]
      [ 71.6  76.8 233.5   0.    0.    0. ]
      [  0.    0.    0.   64.    0.    0. ]
      [  0.    0.    0.    0.   77.    0. ]
      [  0.    0.    0.    0.    0.   78.7]]
+    Phase: forsterite
     Symmetry: orthorhombic
 
 And define the mass density:

{elasticipy-3.0.0 → elasticipy-4.0.0}/docs/source/conf.py

@@ -40,6 +40,8 @@ language = 'english'
 
 html_theme = "sphinx_rtd_theme"
 html_static_path = ['../_static']
+html_logo = "logo/logo_text_whitebg.png"
+html_favicon = 'logo/favicon.png'
 numpydoc_class_members_toctree = False
 autoclass_content = 'both'
 

elasticipy-4.0.0/docs/source/logo/logo.svg

@@ -0,0 +1,100 @@
+<?xml version="1.0" encoding="UTF-8" standalone="no"?>
+<!-- Created with Inkscape (http://www.inkscape.org/) -->
+
+<svg
+   width="35mm"
+   height="35mm"
+   viewBox="0 0 35 35"
+   version="1.1"
+   id="svg1"
+   sodipodi:docname="logo.svg"
+   inkscape:export-filename="favicon.png"
+   inkscape:export-xdpi="23.222858"
+   inkscape:export-ydpi="23.222858"
+   inkscape:version="1.4 (86a8ad7, 2024-10-11)"
+   xmlns:inkscape="http://www.inkscape.org/namespaces/inkscape"
+   xmlns:sodipodi="http://sodipodi.sourceforge.net/DTD/sodipodi-0.dtd"
+   xmlns:xlink="http://www.w3.org/1999/xlink"
+   xmlns="http://www.w3.org/2000/svg"
+   xmlns:svg="http://www.w3.org/2000/svg">
+  <sodipodi:namedview
+     id="namedview1"
+     pagecolor="#ffffff"
+     bordercolor="#000000"
+     borderopacity="0.25"
+     inkscape:showpageshadow="2"
+     inkscape:pageopacity="0.0"
+     inkscape:pagecheckerboard="0"
+     inkscape:deskcolor="#d1d1d1"
+     inkscape:document-units="mm"
+     inkscape:zoom="4.1071317"
+     inkscape:cx="50.887095"
+     inkscape:cy="57.339286"
+     inkscape:window-width="1920"
+     inkscape:window-height="1017"
+     inkscape:window-x="1912"
+     inkscape:window-y="-8"
+     inkscape:window-maximized="1"
+     inkscape:current-layer="layer1"
+     inkscape:export-bgcolor="#ffffff00" />
+  <defs
+     id="defs1">
+    <linearGradient
+       inkscape:collect="always"
+       xlink:href="#linearGradient4671"
+       id="linearGradient11"
+       gradientUnits="userSpaceOnUse"
+       gradientTransform="matrix(0.26458333,0,0,-0.26458333,63.222824,140.3324)"
+       x1="347.23236"
+       y1="-16.795702"
+       x2="318.25162"
+       y2="-62.776524"
+       spreadMethod="pad" />
+    <linearGradient
+       id="linearGradient4671">
+      <stop
+         style="stop-color:#ffd43b;stop-opacity:1"
+         offset="0"
+         id="stop4673" />
+      <stop
+         style="stop-color:#ffe873;stop-opacity:1"
+         offset="1"
+         id="stop4675" />
+    </linearGradient>
+    <linearGradient
+       inkscape:collect="always"
+       xlink:href="#linearGradient4689"
+       id="linearGradient9"
+       gradientUnits="userSpaceOnUse"
+       gradientTransform="matrix(0.26458333,0,0,0.26458333,80.574734,138.09756)"
+       x1="181.30275"
+       y1="-17.131027"
+       x2="256.79184"
+       y2="52.408573" />
+    <linearGradient
+       id="linearGradient4689">
+      <stop
+         style="stop-color:#5a9fd4;stop-opacity:1"
+         offset="0"
+         id="stop4691" />
+      <stop
+         style="stop-color:#306998;stop-opacity:1"
+         offset="1"
+         id="stop4693" />
+    </linearGradient>
+  </defs>
+  <g
+     inkscape:label="Layer 1"
+     inkscape:groupmode="layer"
+     id="layer1"
+     transform="translate(-126.73542,-133.35)">
+    <path
+       id="text1-4-5-2-54-3"
+       style="font-size:192px;text-align:start;writing-mode:lr-tb;direction:ltr;text-anchor:start;display:inline;fill:url(#linearGradient11);fill-opacity:1;stroke-width:0.264583"
+       d="m 150.74998,141.06271 c -0.96996,0 -1.83852,0.007 -2.66237,0.0186 0.52425,0.50605 0.97673,1.07142 1.28261,1.74873 0.50892,1.12689 0.82992,2.28417 0.82992,3.53829 v 1.39991 c 2.66736,0.0162 5.07681,0.0196 8.15455,0.0196 0,-0.87942 0.0155,-1.76812 0.0155,-2.4474 0,-2.51527 -0.32529,-4.27777 -7.62021,-4.27777 z m -11.61583,2.07016 c -0.74411,1.15235 -0.76426,2.93071 -0.76946,5.67458 v 0.0248 c 0,1.14102 0.14056,2.1084 0.42168,2.90215 0.11054,0.30092 0.24218,0.57875 0.39378,0.83406 1.45247,-0.19775 2.71407,-2.31732 2.71404,-4.72736 3e-5,-2.43536 -1.28837,-4.55112 -2.76004,-4.70824 z m 16.1029,0.14314 a 1.3449186,1.3449186 0 0 1 1.34513,1.34514 1.3449186,1.3449186 0 0 1 -1.34513,1.34513 1.3449186,1.3449186 0 0 1 -1.34462,-1.34513 1.3449186,1.3449186 0 0 1 1.34462,-1.34514 z m -5.05448,5.9893 c -0.0858,2.29191 -0.60955,4.25462 -1.45572,5.93038 -0.90878,1.81756 -2.28539,2.88556 -4.03025,3.84887 -1.74485,0.96331 -3.67874,1.21646 -6.13244,1.21646 -0.067,0 -0.12977,-0.004 -0.19585,-0.005 -3.2e-4,0.0383 -0.004,0.0729 -0.004,0.11162 v 0.0248 c 0.0134,7.04478 0.0291,7.74475 12.38529,7.74475 7.29492,0 7.62021,-1.7625 7.62021,-4.27777 0,-0.67928 -0.0155,-1.56798 -0.0155,-2.4474 -4.40154,0 -7.05527,-0.01 -11.76001,0.046 -0.74191,0.01 -0.77069,-0.8692 -0.01,-0.8692 4.27995,0 5.96966,-1.55757 5.96966,-6.01823 0,-2.77849 -0.67428,-4.46657 -2.37195,-5.30562 z" />
+    <path
+       id="text1-1-5-1"
+       style="font-weight:bold;font-size:192px;font-family:Bahnschrift;-inkscape-font-specification:'Bahnschrift Bold';display:inline;fill:url(#linearGradient9);stroke-width:0.705556"
+       d="m 143.93283,133.35 c -2.87272,0 -8.44558,0.27355 -10.72338,1.48363 -1.5875,0.84336 -2.80293,2.06706 -3.64629,3.6711 -0.82682,1.60402 -1.24023,3.53053 -1.24023,5.77949 v 4.19199 c 0,2.36472 0.41341,4.39043 1.24023,6.07715 0.84336,1.67018 2.05879,2.94349 3.64629,3.81992 1.5875,0.87644 3.49746,1.31465 5.72988,1.31465 2.23242,0 4.14239,-0.43821 5.72989,-1.31465 1.5875,-0.87643 2.79466,-2.14147 3.62148,-3.79511 0.84336,-1.67019 1.26504,-3.67937 1.26504,-6.02754 v -2.03399 c 0,-1.14101 -0.23151,-2.22415 -0.69453,-3.24941 -0.46302,-1.02526 -1.14102,-1.87689 -2.03399,-2.55488 -0.29302,-0.22791 -0.60567,-0.42439 -0.93431,-0.59532 0,0 6.299,0.0465 8.87181,0 10e-4,-0.72917 0.0207,-1.55216 0.0207,-2.13475 0,-3.69368 -0.73316,-4.63228 -10.85257,-4.63228 z m 7.72666,2.2779 a 1.3383084,1.3383084 0 0 1 1.33841,1.33842 1.3383084,1.3383084 0 0 1 -1.33841,1.3379 1.3383084,1.3383084 0 0 1 -1.33791,-1.3379 1.3383084,1.3383084 0 0 1 1.33791,-1.33842 z m -12.67365,6.68641 c 1.97408,3e-4 3.57433,2.43687 3.57446,5.44257 6e-5,3.00589 -1.60025,5.44278 -3.57446,5.44308 -1.97441,1.4e-4 -3.57503,-2.43688 -3.57497,-5.44308 1.3e-4,-3.00601 1.60069,-5.44271 3.57497,-5.44257 z" />
+  </g>
+</svg>