擅长:python、mysql、java
<p>我最终发现了问题所在。当试图打印CalcStress时,我没有包含参数。我现在做的这段代码给出了我想要的数组中的答案</p>
<pre><code>import numpy as np
Ixx = 5.6667 * 10**5
Iyy = 1.493 * 10**5
Ixy = - .8 * 10**5
x = 66.7
y = -50
Sy = np.array([400, 600, 800, 1000, 1200, -400, -800,-1000, -1200])
Sx = np.array([600, 600, 400, -400, -600, 600, 400, -400, 600])
Mx = -Sx * 1000
My = Sy * 1000
DEN = (Ixx * Iyy) - (Ixy**2)
def CalcStress(Ixx, Iyy, Ixy, Mx, My, x, y, DEN):
n = My.shape[0]
result = np.zeros_like(Mx)
for i in range(0, n):
result[i] = (((((My[i] * Ixx) - (Mx[i]*Ixy)) * x) + ((Mx[i] * Iyy) - (My[i] * Ixy)) * y)) / DEN
return result
print(CalcStress(Ixx, Iyy, Ixy, Mx, My, x, y, DEN))
</code></pre>
<p>结果=[189275356421502,-156,-334,-443,-502]</p>