1

我想知道如何绘制 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 求解它,但我还没有尝试过这种方法。

4

0 回答 0