材料力学之应变分析算法：大变形应变分析：大变形分析中的几何非线性.docx

:param coordinates:

:return:

"""

n = displacement.shape[0]

strain = np.zeros((n, 3, 3))

grad_u = np.gradient(displacement, coordinates)

for i in range(n):

strain[i] = 0.5 * (grad_u[i] + grad_u[i].T + np.dot(grad_u[i], grad_u[i].T))