Phân tích tiếp xúc (phần 1)¶
Phân tích này sử dụng dữ liệu trong tutorial/09_contact_hertz.
Đối tượng phân tích¶
Đây là bài toán tiếp xúc Hertz. Hình dạng của đối tượng phân tích được thể hiện trong Hình 4.9.1 và dữ liệu lưới được thể hiện trong Hình 4.9.2.
| Hạng mục | Nội dung | Ghi chú | Tham chiếu |
|---|---|---|---|
| Loại phân tích | Phân tích tĩnh phi tuyến (đàn dẻo, tiếp xúc) | !SOLUTION,TYPE=NLSTATIC !CONTACT | |
| Số nút | 408 | ||
| Số phần tử | 168 | ||
| Loại phần tử | Phần tử lục diện bậc nhất 8 nút | !ELEMENT,TYPE=361 | |
| Tên vật liệu | MAT1 | !MATERIAL,NAME=MAT1 | |
| Tính chất vật liệu | ELASTIC | !ELASTIC | |
| Điều kiện biên | Ràng buộc, chuyển vị chỉ định | ||
| Bộ giải ma trận | CG/SSOR | !SOLVER,METHOD=CG,PRECOND=1 |
Hình 4.9.1: Hình dạng đối tượng phân tích
Hình 4.9.2: Dữ liệu lưới của đối tượng phân tích
Nội dung phân tích¶
Thực hiện phân tích tiếp xúc bằng phương pháp nhân tử Lagrange tăng cường, trong đó áp đặt chuyển vị chỉ định theo hướng nén lên bề mặt trên của mô hình một phần tư đĩa. Dữ liệu điều khiển phân tích được trình bày dưới đây.
Dữ liệu điều khiển phân tích cgs3.cnt¶
# Tệp điều khiển FISTR
## Điều khiển phân tích
!VERSION
3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## Điều khiển bộ giải
### Điều kiện biên
!BOUNDARY, GRPID=1
ALL, 3, 3, 0.0
BOTTOM, 2, 2, 0.0
CENTER, 1, 1, 0.0
UPPER, 2, 2, -0.306
!CONTACT_ALGO, TYPE=ALAGRANGE
!CONTACT, GRPID=1
CP1, 0.0
### Bước
!STEP, SUBSTEPS=5, CONVERG=1.0e-5
BOUNDARY, 1
CONTACT, 1
### Vật liệu
!MATERIAL, NAME=MAT1
!ELASTIC
1100.0, 0.0
### Thiết lập bộ giải
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES
1000, 1
1.0e-10, 1.0, 0.0
## Điều khiển hậu xử lý
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
Quy trình phân tích¶
Thực thi lệnh chạy FrontISTR fistr1.
(Chạy với 4 luồng.)
Kết quả phân tích¶
Đối với kết quả phân tích ở bước con thứ 5, REVOCAP_PrePost được dùng để tạo hình biến dạng có thêm đường đồng mức chuyển vị theo hướng y, như thể hiện trong Hình 4.9.3. Ngoài ra, một phần tệp nhật ký kết quả phân tích được trình bày dưới đây dưới dạng dữ liệu số của kết quả phân tích.
Hình 4.9.3: Kết quả phân tích biến dạng và chuyển vị theo hướng y
Nhật ký kết quả phân tích 0.log¶
fstr_setup: OK
#### Result step= 0
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 0.0000E+00 1 0.0000E+00 1
//U2 0.0000E+00 1 0.0000E+00 1
//U3 0.0000E+00 1 0.0000E+00 1
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Global Summary @Node :Max/IdMax/Min/IdMin####
//U1 0.0000E+00 1 0.0000E+00 1
//U2 0.0000E+00 1 0.0000E+00 1
//U3 0.0000E+00 1 0.0000E+00 1
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
...
#### Result step= 5
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 1.3079E-02 31 -4.1442E-02 54
//U2 4.9597E-03 1008 -3.0604E-01 69
//U3 0.0000E+00 1 0.0000E+00 1
//E11 7.1171E-03 1049 -4.7627E-02 1
//E22 2.0644E-03 1012 -5.3428E-02 50
//E33 9.3301E-03 50 -4.0103E-03 54
//E12 4.3231E-02 1000 -2.7666E-02 54
//E23 1.3028E-12 1046 -1.2515E-12 2
//E31 3.5660E-13 48 -2.1028E-12 51
//S11 7.9455E+00 1049 -4.9987E+01 1
//S22 2.2819E+00 1012 -5.3282E+01 50
//S33 1.1359E+01 50 -4.4937E+00 54
//S12 2.0802E+01 1000 -1.4668E+01 54
//S23 6.7682E-10 1046 -6.3666E-10 2
//S31 2.0235E-10 48 -1.1497E-09 51
//SMS 6.0944E+01 1000 1.0198E-01 69
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.6311E-03 1024 -3.5443E-02 8
//E22 1.7091E-03 1031 -5.4262E-02 8
//E33 1.5903E-04 7 -1.6631E-05 101
//E12 1.8590E-02 1009 -2.6073E-02 14
//E23 1.5204E-13 1000 -4.1521E-13 16
//E31 1.3902E-13 5 -8.2391E-13 7
//S11 5.2516E+00 1024 -3.8099E+01 8
//S22 1.8872E+00 1031 -5.5262E+01 8
//S33 2.2824E-01 7 -1.8263E-02 101
//S12 9.8516E+00 1009 -1.3352E+01 14
//S23 7.8767E-11 1000 -2.1075E-10 16
//S31 7.7553E-11 5 -4.4780E-10 7
//SMS 4.9470E+01 8 2.4786E-01 1007
##### Global Summary @Node :Max/IdMax/Min/IdMin####
//U1 1.3079E-02 31 -4.1442E-02 54
//U2 4.9597E-03 1008 -3.0604E-01 69
//U3 0.0000E+00 1 0.0000E+00 1
//E11 7.1171E-03 1049 -4.7627E-02 1
//E22 2.0644E-03 1012 -5.3428E-02 50
//E33 9.3301E-03 50 -4.0103E-03 54
//E12 4.3231E-02 1000 -2.7666E-02 54
//E23 1.3028E-12 1046 -1.2515E-12 2
//E31 3.5660E-13 48 -2.1028E-12 51
//S11 7.9455E+00 1049 -4.9987E+01 1
//S22 2.2819E+00 1012 -5.3282E+01 50
//S33 1.1359E+01 50 -4.4937E+00 54
//S12 2.0802E+01 1000 -1.4668E+01 54
//S23 6.7682E-10 1046 -6.3666E-10 2
//S31 2.0235E-10 48 -1.1497E-09 51
//SMS 6.0944E+01 1000 1.0198E-01 69
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.6311E-03 1024 -3.5443E-02 8
//E22 1.7091E-03 1031 -5.4262E-02 8
//E33 1.5903E-04 7 -1.6631E-05 101
//E12 1.8590E-02 1009 -2.6073E-02 14
//E23 1.5204E-13 1000 -4.1521E-13 16
//E31 1.3902E-13 5 -8.2391E-13 7
//S11 5.2516E+00 1024 -3.8099E+01 8
//S22 1.8872E+00 1031 -5.5262E+01 8
//S33 2.2824E-01 7 -1.8263E-02 101
//S12 9.8516E+00 1009 -1.3352E+01 14
//S23 7.8767E-11 1000 -2.1075E-10 16
//S31 7.7553E-11 5 -4.4780E-10 7
//SMS 4.9470E+01 8 2.4786E-01 1007