三向应力状态求剪应力

Question

请问大家，单元处于三向应力状态，已知单元所受的三个正应力，三个主应力，能否求出单元所受的三个剪应力呢，如何求解？

Answer 1

阿里嘎多学长整理AIGC生成，因移动端显示问题导致当前答案未能完全显示，请使用PC端查看更加详细的解答过程

三向应力状态求剪应力

在三向应力状态下，单元受到三个正应力$\sigma_1$, $\sigma_2$, $\sigma_3$和三个主应力$\varepsilon_1$, $\varepsilon_2$, $\varepsilon_3$。剪应力可以通过使用Mohr's circle法或者Voigt notation来求解。

使用Mohr's circle法，首先计算主应力矩阵：

sigma_1 = 10 MPa
sigma_2 = 5 MPa
sigma_3 = 3 MPa

然后，绘制Mohr's circle：

import matplotlib.pyplot as plt

sigma1 = 10
sigma2 = 5
sigma3 = 3

r = ((sigma1 - sigma2) ** 2 + (sigma2 - sigma3) ** 2 + (sigma3 - sigma1) ** 2) / 2
theta = 0.5 * np.arctan2(2 * sigma2 * (sigma1 - sigma3), (sigma1 + sigma3 - 2 * sigma2))

circle = plt.Circle((sigma1 + sigma2 + sigma3) / 3, r, fill=False)

plt.gca().add_artist(circle)

plt.plot([sigma1, sigma1], [0, r], 'k-')
plt.plot([sigma2, sigma2], [0, r], 'k-')
plt.plot([sigma3, sigma3], [0, r], 'k-')

plt.show()

接着，根据Mohr's circle，计算剪应力：

tau = r * np.sin(theta)

使用Voigt notation，可以直接计算剪应力：

 tau = np.sqrt((sigma_1 - sigma_2) ** 2 + (sigma_2 - sigma_3) ** 2 + (sigma_3 - sigma_1) ** 2)

通过上述方法，计算出的剪应力为：

tau_x = sqrt((sigma_1 - sigma_2) ** 2 + (sigma_2 - sigma_3) ** 2) = sqrt((10 - 5) ** 2 + (5 - 3) ** 2) = 5.49 MPa
tau_y = sqrt((sigma_2 - sigma_3) ** 2 + (sigma_3 - sigma_1) ** 2) = sqrt((5 - 3) ** 2 + (3 - 10) ** 2) = 7.07 MPa
tau_z = sqrt((sigma_3 - sigma_1) ** 2 + (sigma_1 - sigma_2) ** 2) = sqrt((3 - 10) ** 2 + (10 - 5) ** 2) = 8.61 MPa

注意：上述方法假设单元所受的三维应力状态是无关联的，即三轴之间的夹角为零。如果存在关联，需要使用更复杂的方法来求解。