Thành lập ma trận độ cứng xét đến hiệu ứng bậc hai của phần tử khung bằng hàm chuyển vị tổng quát được giải từ phương trình vi phân