弹性力学材料模型：弹塑性材料：塑性理论与屈服准则.docx

import numpy as np

#

stress_tensor = np.array([[100, 50, 0],

[50, 100, 0],

[0, 0, 0]])

#

stress_dev = stress_tensor - np.mean(np.diag(stress_tensor)) * np.eye(3)

#

von_mises_stress = np.sqrt(3/2 * np.dot(stress_dev.flat, stress_dev.flat))

print(f"von Mises Stress: {von_mises_stress}")