Phân tích chuyển vị lớn khung phẳng đàn dẻo sử dụng công thức lagrange cập nhật