LỰC NÚT TƯƠNG ĐƯƠNG BÊN TRONG CỦA PHẦN TỬ ỨNG SUẤT PHẲNG TRONG BÀI TOÁN PHẦN TỬ HỮU HẠN PHI TUYẾN ThS. NGUYỄN ĐẠI VIÊN Sở Xây dựng Thừa Thiên Huế 1. Đặt vấn đề Trong thực tế tính toán kết cấu, chúng ta thường phải phân tích bài toán phi tuyến. Bài toán này đưa về giải phương trình chứa các số hạng phi tuyến đối với ẩn số. Nói chung, không thể giải một cách chính xác dưới dạng đóng những phương trình phi tuyến mà phải dùng các thuật toán đúng đắn, trong đó tiêu chuẩn hội tụ là vấn đề cần quan tâm. Trong phương pháp phần tử hữu hạn (PTHH), tải trọng tác dụng lên hệ được thay thế gần đúng bằng một hệ lực đặt tại các nút của phần tử (gọi là lực nút tương đương bên ngoài). Tương tự, các thành phần ứng suất của hệ có thể được thay thế gần đúng bằng một hệ lực đặt tại các nút của phần tử (tạm gọi là lực nút tương đương bên trong). Theo nguyên tắc cân bằng của hệ thì hai hệ lực này phải cân bằng nhau. Tuy nhiên, trong các bước lặp của bài toán phân tích phi tuyến, sự chênh lệch giữa hai hệ lực này luôn tồn tại và có xu hướng giảm dần khi số vòng lặp tăng lên. Khi độ chênh lệch này nhỏ hơn một giá trị quy định, ta nói bài toán hội tụ. Bài báo nhằm xác định lực nút tương đương bên trong của phần tử chữ nhật trong trạng thái ứng suất phẳng. Đây là dạng bài toán có phạm vi ứng dụng tương đối rộng rãi trong ngành xây dựng, chế tạo máy bay, đóng tàu 2. Cơ sở lý thuyết a. Phương trình cân bằng của hệ PTHH ở vòng lặp i có thể viết: 0 ii FR (1) Trong đó: R i và F i là véctơ lực nút tương đương bên ngoài và bên trong của hệ ở vòng lặp i . Giả sử lời giải tại vòng lặp i đã biết, cần xác định lời giải tại vòng lặp ( i +1). Viết lại (1) cho vòng lặp ( i +1): 0 11   ii FR (2) Vì lời giải tại vòng lặp i đã biết nên ta có: F F F ii   1 (3) Với F là số gia lực nút tương đương bên trong. Mặt khác, véctơ F có thể được viết: UKF i . (4) Trong đó: K i là ma trận độ cứng của hệ ở vòng lặp thứ i và U là véctơ gia số chuyển vị nút. Thay (3), (4) vào (2): iii FRUK  1 . (5) Từ đó, tính được chuyển vị tại vòng lặp i +1: UUU ii  1 (6) Với U là nghiệm của (5). Vòng lặp trên đây được dùng cho nhiều phương pháp khác nhau. Trong đó, phương pháp Newton- Raphson (xem hình 1) được dùng phổ biến vì có độ hội tụ cao. Hình 1. Phương pháp lặp Newton - Raphson Chuyển vị Tải trọng U i U i+1 F (0) i F i+1 R i R i+1 U (2) U (1) K (0) i K (1) i R i+1 – F (1) i R i+1 – F (0) i Nghiệm của (5) là chuyển vị gây ra bởi hiệu số giữa lực nút tương đương bên ngoài ở vòng lặp đang xét và lực nút tương đương bên trong ở vòng lặp trước . Do vậy, việc xác định lực nút tương đương bên trong tại mỗi vòng lặp là vấn đề cần thiết đối với việc giải bài toán PTHH phi tuyến theo phương pháp lặp. b. Xét phần tử tứ giác 4 nút định vị trong hệ tọa độ Oxy , có trạng thái ứng suất phẳng. Để tiện tính toán, ta sử dụng phần tử tham chiếu định vị trong hệ toạ độ O  . Trong hệ này, các thành phần toạ độ  và  của những điểm trên phần tử tham chiếu mang các giá trị trong đoạn [-1, 1] (xem hình 2a và 2b). Mối quan hệ giữa tọa độ một điểm M(x,y) của phần tử thực trong hệ Oxy và điểm tương ứng M’(  ,  ) của phần tử tham chiếu trong hệ O  :              4 1 4 1 .),( .),( i ii i i i yNy xNx   (7) Với N i là các hàm dạng của phần tử. c. Lực nút tương đương bên trong của phần tử theo [1]:         v e T e dvBF  (8) Trong đó: [ B ] là ma trận đạo hàm; { F e } gồm các thành phần lực nút tương đương theo phương x, y của các nút. Vấn đề đặt ra là làm thế nào để tính được (8) khi kết quả của bài toán PTHH không phải là các hàm ứng suất mà là các trị số ứng suất rời rạc tại các điểm nút:              xy yy xx e     (9) d. Thay vi phân thể tích dv bởi ( e dx dy ), với e là chiều dày phần tử; (8) trở thành:         s e T e dydxeBF  (10) Với quan hệ (7), đổi biến số trong tính tích phân kép với:                       ddJBF yy xx T e ),( ),( (11) Với J là ma trận Jacobi. Trường hợp phần tử chữ nhật:                 1 1 1 1 1 1    ddJBF T e (12) Tích phân (12) có thể tính gần đúng bằng phương pháp Gauss:         1 1 1 1 1 1 ),( ),( r i r j jiji fwwddf  (13) Với: 1 Hình 2a. Phần tử thực -1 -1   l k y x o i j   o 1 Hình 2b. Phần tử tham chiếu (  i ,  j ) - Tọa độ các điểm Gauss nằm trên diện tích của phần tử tham chiếu, có giá trị thuộc đoạn [- 1,1]; w i , w j - Trọng số ứng với điểm Gauss có tọa độ (  i ,  j ); r - Số điểm Gauss nằm trên một phương  hoặc  . Ta chọn r = 2 thì tọa độ điểm Gauss là  0,577350269189626 và trọng số tương ứng là 1. Như vậy, ta có thể tính được lực nút tương đương bên trong của phần tử từ ứng suất tại các nút theo (13). Ghép các lực nút tương đương bên trong của các phần tử { F e }, ta sẽ có ma trận lực nút tương đương bên trong của hệ { F }. 3. Ví dụ Xét vách cứng chịu tải trọng ngang, liên kết ngàm ở chân có kích thước như hình 3. Vật liệu có môđun đàn hồi E = 2,65.10 9 KG/m 2 ; hệ số poisson là 0,3. - Chia hệ thành 4 phần tử và đánh số nút và phần tử như hình 3. - Ma trận lực nút tương đương bên ngoài từ tải trọng đã cho:     yxyxyxyxyxyxyxyxyx T rrrrrrrrrrrrrrrrrrR 998877665544332211      000001000000000200000000010000 T R - Lập trình và tính toán theo tích phân (13), ta có ma trận lực nút tương đương bên trong:     yxyxyxyxyxyxyxyxyx T ffffffffffffffffffF 998877665544332211      00000100000000020000396211614075615414037812318  T F Nhận xét: - Tại các nút không có liên kết (4, 5, 6, 7, 8, 9): { R } và { F } cân bằng nhau; - Lực nút tương đương bên ngoài r 1x thực ra không gây ra nội lực trong hệ. Do vậy, cộng giá trị r 1x vào f 1x của { F }, ta có:     00000100000000020000396211614075615414037822318 *  T F Lúc này, các giá trị f ix , f iy ( i = 1 ÷ 3) chính là các phản lực của hệ tại các nút liên kết. 4. Kết luận - Lực nút tương đương bên trong có thể tính từ các giá trị ứng suất rời rạc tại các điểm nút của hệ bằng phương pháp gần đúng Gauss. - Trong quá trình giải bài toán đàn hồi phi tuyến theo phương pháp lặp, giữa ma trận lực nút tương đương bên ngoài ở vòng lặp hiện tại luôn có chênh lệch với ma trận lực nút tương đương bên trong ở vòng lặp trước. Giá trị chênh lệch này được dùng làm tiêu chuẩn hội tụ để giải bài toán phi tuyến theo phương pháp lặp. - Khi bài toán hội tụ, ta nhận được ma trận lực nút tương đương bên trong gồm các thành phần: + Tại các nút không có liên kết: ma trận lực nút tương đương bên trong (do nội lực) là cân bằng với ma trận lực nút tương đương bên ngoài (do ngoại lực). + Tại các nút có liên kết: các số hạng của ma trận lực nút tương đương bên trong chính là phản lực của hệ. 10 m q = 20000 KG/m Hình 3. Ví dụ 0,2m 5m 5m 10 m 1 2 3 4 5 6 7 8 9 1 2 3 4 TÀI LIỆU THAM KHẢO 1. M. Y. H. BANGASH. Concrete and conrete structures: Numerical modelling and Application. Elsevier science publisher Ltd. 1989. 2. HỒ ANH TUẤN - TRẦN BÌNH. Phương pháp phần tử hữu hạn . NXB Khoa học kỹ thuật, Hà Nội, 1978. . lực nút tương đương bên trong của phần tử từ ứng suất tại các nút theo (13). Ghép các lực nút tương đương bên trong của các phần tử { F e }, ta sẽ có ma trận lực nút tương đương bên trong của. nói bài toán hội tụ. Bài báo nhằm xác định lực nút tương đương bên trong của phần tử chữ nhật trong trạng thái ứng suất phẳng. Đây là dạng bài toán có phạm vi ứng dụng tương đối rộng rãi trong. LỰC NÚT TƯƠNG ĐƯƠNG BÊN TRONG CỦA PHẦN TỬ ỨNG SUẤT PHẲNG TRONG BÀI TOÁN PHẦN TỬ HỮU HẠN PHI TUYẾN ThS. NGUYỄN ĐẠI VIÊN Sở Xây dựng Thừa Thiên Huế 1. Đặt vấn đề Trong thực tế tính toán