我想知道如何绘制 3D 屈服面,如本页http://en.wikipedia.org/wiki/Von_Mises_yield_criterion所示。我一直在尝试绘制 von Mises 屈服面作为开始,但完全没有成功。我想我根本不知道如何开始。如果有人可以向我展示一个示例代码,例如 von Mises,那么我可以使用它来绘制其他一些屈服面,例如 Cam Clay 和 Modified Cam Clay。
最初我以为我可以做这样的事情,但它不起作用:
for i in range(0, MAXSTEP):
stressx = stressx + STEPSIZE
stressy = -STEPSIZE
for j in range(0, MAXSTEP):
stressy = stressy + STEPSIZE
stressz = -STEPSIZE
for k in range(0, MAXSTEP):
stressz = stressz + STEPSIZE
tI1 = calcI1(stressx, stressy, stressz)
tJ2 = calcJ2(stressx, stressy, stressz)
tyieldf = calcF(tI1, ycon)
tsqrtJ2 = np.sqrt(tJ2)
tdiff = abs(tsqrtJ2 - tyieldf)
if tdiff < TOL:
sigx[count] = stressx
sigy[count] = stressy
sigz[count] = stressy
I1[count] = tI1
sqrtJ2[count] = tsqrtJ2
diff[count] = tdiff
count = count + 1
X = np.outer(sigx[0:count], sigx[0:count])
Y = np.outer(sigy[0:count], sigy[0:count])
Z = np.outer(sigz[0:count], sigz[0:count])
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X, Y, Z)
ax.set_xlabel('X Label')
ax.set_ylabel('Y Label')
ax.set_zlabel('Z Label')
plt.show()
main()
我想我不知道如何生成数据。我有方程,但不知道如何使用它们来生成曲面图所需的数据。我想我需要生成一个 XY 网格,然后为 Z 求解它,但我还没有尝试过这种方法。