You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
61 lines
2.1 KiB
61 lines
2.1 KiB
2 years ago
|
import numpy as np
|
||
|
|
|
||
|
|
|
||
|
|
# The indices of the full stiffness matrix of (orthorhombic) interest
|
||
|
|
voigt_notation = [(0, 0), (1, 1), (2, 2), (1, 2), (0, 2), (0, 1)]
|
||
|
|
|
||
|
|
|
||
|
|
def full_3x3_to_voigt_6_index(i, j):
|
||
|
|
if i == j:
|
||
|
|
return i
|
||
|
|
return 6 - i - j
|
||
|
|
|
||
|
|
|
||
|
|
def voigt_6_to_full_3x3_strain(strain_vector):
|
||
|
|
"""
|
||
|
|
Form a 3x3 strain matrix from a 6 component vector in Voigt notation
|
||
|
|
"""
|
||
|
|
e1, e2, e3, e4, e5, e6 = np.transpose(strain_vector)
|
||
|
|
return np.transpose([[1.0 + e1, 0.5 * e6, 0.5 * e5],
|
||
|
|
[0.5 * e6, 1.0 + e2, 0.5 * e4],
|
||
|
|
[0.5 * e5, 0.5 * e4, 1.0 + e3]])
|
||
|
|
|
||
|
|
|
||
|
|
def voigt_6_to_full_3x3_stress(stress_vector):
|
||
|
|
"""
|
||
|
|
Form a 3x3 stress matrix from a 6 component vector in Voigt notation
|
||
|
|
"""
|
||
|
|
s1, s2, s3, s4, s5, s6 = np.transpose(stress_vector)
|
||
|
|
return np.transpose([[s1, s6, s5],
|
||
|
|
[s6, s2, s4],
|
||
|
|
[s5, s4, s3]])
|
||
|
|
|
||
|
|
|
||
|
|
def full_3x3_to_voigt_6_strain(strain_matrix):
|
||
|
|
"""
|
||
|
|
Form a 6 component strain vector in Voigt notation from a 3x3 matrix
|
||
|
|
"""
|
||
|
|
strain_matrix = np.asarray(strain_matrix)
|
||
|
|
return np.transpose([strain_matrix[..., 0, 0] - 1.0,
|
||
|
|
strain_matrix[..., 1, 1] - 1.0,
|
||
|
|
strain_matrix[..., 2, 2] - 1.0,
|
||
|
|
strain_matrix[..., 1, 2] + strain_matrix[..., 2, 1],
|
||
|
|
strain_matrix[..., 0, 2] + strain_matrix[..., 2, 0],
|
||
|
|
strain_matrix[..., 0, 1] + strain_matrix[..., 1, 0]])
|
||
|
|
|
||
|
|
|
||
|
|
def full_3x3_to_voigt_6_stress(stress_matrix):
|
||
|
|
"""
|
||
|
|
Form a 6 component stress vector in Voigt notation from a 3x3 matrix
|
||
|
|
"""
|
||
|
|
stress_matrix = np.asarray(stress_matrix)
|
||
|
|
return np.transpose([stress_matrix[..., 0, 0],
|
||
|
|
stress_matrix[..., 1, 1],
|
||
|
|
stress_matrix[..., 2, 2],
|
||
|
|
(stress_matrix[..., 1, 2] +
|
||
|
|
stress_matrix[..., 1, 2]) / 2,
|
||
|
|
(stress_matrix[..., 0, 2] +
|
||
|
|
stress_matrix[..., 0, 2]) / 2,
|
||
|
|
(stress_matrix[..., 0, 1] +
|
||
|
|
stress_matrix[..., 0, 1]) / 2])
|