Question

我想知道如何绘制3D屈服面，如本页http://en.wikipedia.org/wiki/Von_Mises_yield_criterion所示。我一直试图将冯米塞斯屈服面作为一个开始，但根本没有成功。我想我根本不知道如何开始。如果有人能给我看一个例如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解决它，但我还没有尝试过这种方法。

使用matplotlib绘制3D屈服面

0 个答案: