弹性力学数值方法：有限差分法(FDM)：二维弹性问题的有限差分法.docx

for iteration in range(max_iterations):

u_old = u.copy()

for i in range(1, nx-1):

for j in range(1, ny-1):

u[i, j] = (u[i+1, j] + u[i-1, j] + u[i, j+1] + u[i, j-1] - h**2 * f(X[i, j], Y[i, j])) / 4

if np.linalg.norm(u - u_old) < tolerance:

break

#