计算欧拉角的晶界取向差，可先将欧拉角转换为旋转矩阵，再通过旋转矩阵计算晶界取向差。以下是具体步骤和Python代码示例： ### 方法 1. 将欧拉角转换为旋转矩阵。 2. 计算两个旋转矩阵的乘积，得到相对旋转矩阵。 3. 从相对旋转矩阵中提取欧拉角，即为晶界取向差。 ### Python代码示例 ```python import numpy as np from scipy.spatial.transform import Rotation as R def euler_to_rotation_matrix(euler_angles): """ 将欧拉角转换为旋转矩阵 :param euler_angles: 欧拉角 (α, β, γ)，单位为度 :return: 旋转矩阵 """ r = R.from_euler('zyx', euler_angles, degrees=True) return r.as_matrix() def calculate_grain_boundary_misorientation(euler1, euler2): """ 计算两个欧拉角表示的晶界取向差 :param euler1: 第一个欧拉角 (α, β, γ)，单位为度 :param euler2: 第二个欧拉角 (α, β, γ)，单位为度 :return: 晶界取向差的欧拉角，单位为度 """ # 将欧拉角转换为旋转矩阵 R1 = euler_to_rotation_matrix(euler1) R2 = euler_to_rotation_matrix(euler2) # 计算相对旋转矩阵 R_relative = np.dot(R2, np.linalg.inv(R1)) # 将相对旋转矩阵转换为欧拉角 r_relative = R.from_matrix(R_relative) misorientation = r_relative.as_euler('zyx', degrees=True) return misorientation # 示例使用 euler1 = [30, 45, 60] # 第一个欧拉角 euler2 = [15, 30, 45] # 第二个欧拉角 misorientation = calculate_grain_boundary_misorientation(euler1, euler2) print("晶界取向差的欧拉角 (α, β, γ):", misorientation) ``` ### 代码解释 1. `euler_to_rotation_matrix` 函数：将欧拉角转换为旋转矩阵，使用 `scipy.spatial.transform.Rotation` 类实现。 2. `calculate_grain_boundary_misorientation` 函数：计算两个欧拉角表示的晶界取向差，先将欧拉角转换为旋转矩阵，再计算相对旋转矩阵，最后将相对旋转矩阵转换为欧拉角。 3. 示例使用：定义两个欧拉角 `euler1` 和 `euler2`，调用 `calculate_grain_boundary_misorientation` 函数计算晶界取向差，并打印结果。 ###