材料力学之动力学分析算法：模态分析：材料力学基础理论.Tex.header.docx

M = [1, 0; 0, 1]; %

K = [2, -1; -1, 2]; %

[V, D] = eig(K, M);

omega = sqrt(diag(D)); %

phi = V; %

alpha = 0.1; %

beta = 0.05; %

C = alpha*M + beta*K; %

epsilon = np.array([[0.001, 0.0005, 0.0],

[0.0005, 0.002, 0.0],

[0.0, 0.0, 0.0]])

stress = E / (1 - nu**2) * (epsilon + nu * np.transpose(epsilon) - (nu/(1-2*nu)) * np.trace(epsilon)

* np.eye(3))

strain_energy_density = 0.5 * np.sum(stress * epsilon)