− Ứng xử phi tuyến ví dụ như độ võng lớn, biến dạng lớn, bài toán tiếp xúc, chảy dẻo, siêu đàn hồi, từ biến ...
• Phân tích động lực học:
− Bao gồm hiệu ứng khối lượng và giảm chấn.
− Phân tích điều hòa: xác định ứng xử của kết cấu khi tải trọng có dạng hình sin với biên độ và tần số xác định.
− Phân tích động lực học tức thời (Transient Dynamic Analysis): xác định ứng xử của kết cấu khi tải trọng thay đổi theo thời gian và có thể bao gồm cả ứng xử phi tuyến.
• Động lực học biến dạng lớn:
− Dùng để mô phỏng biến dạng rất lớn khi lực quán tính đóng vai trò quyết định.
− Dùng để mô phỏng các bài toán va chạm, phá huỷ, tạo hình nhanh,…
Hình 3.3 – Phân tích va chạm của một thí nghiệm đối với ô tô
Phân tích nhiệt
• Phân tích nhiệt được dùng để xác định trường phân bố nhiệt độ trong một vật thể. Các đại lượng đáng quan tâm khác bao gồm: lượng nhiệt mất đi hoặc tăng lên, gradient nhiệt, và dòng nhiệt.
• Tất cả 3 dạng truyền nhiệt cơ bản đều có thể được phân tích và mô phỏng: dẫn nhiệt, đối lưu, bức xạ.
• Trạng thái ổn định (Steady-State): Bỏ qua các ảnh hưởng phụ thuộc thời gian.
• Trạng thái tức thời hay chưa ổn định (Transient):
• Để xác định nhiệt độ và một số đại lượng khác như một hàm của thời gian.
Hình 3.4 – Phân bố nhiệt trong kết cấu
Phân tích điện từ
• Mô phỏng các thiết bị sử dụng nguồn điện 1 chiều, nguồn xoay chiều tần số thấp, các tín hiệu tức thời ngắn tần số thấp. Ví dụ: thiết bị khởi động từ (solenoid), các động cơ, máy biến thế.
• Các thông số đáng quan tâm bao gồm : mật độ thông lượng từ, cường độ từ trường, lực và mô men từ, trở kháng, độ tự cảm, dòng điện xoáy, công suất mất mát, và dòng rò.
• Mô phỏng các thiết bị truyền sóng điện từ.
• Ví dụ: các thiết bị thu vi sóng và sóng radio, dẫn sóng, thiết bị kết nối đồng trục.
• Các đại lượng đáng quan tâm gồm có: các thông số S, nhân tố Q, tổn thất đường về, tổn hao điện môi và tổn hao dẫn điện, và các trường điện và từ.
Hình 3.5 – Mật độ dòng chảy điện từ của van kiểm soát chất lỏng solenoid
• Tính toán trường điện khi kích thích bằng điện áp hoặc tích điện như thiết bị cao áp, các hệ vi cơ điện tử (MEMS), đường truyền.
• Các đại lượng điển hình là cường độ và điện dung của trường điện.
• Độ dẫn điện: để tính toán dòng điện trong dây dẫn khi áp đặt một điện áp
• Kết nối mạch: để kết nối mạch điện với các thiết bị điện từ.
• Các kiểu phân tích điện từ:
• Phân tích tĩnh: tính toán từ trường của dòng 1 chiều hoặc nam châm vĩnh cửu.
• Phân tích điều hòa: tính toán từ trường của dòng điện xoay chiều.
• Phân tích tức thời: được sử dụng với từ trường thay đổi theo thời gian. (adsbygoogle = window.adsbygoogle || []).push({});
Tính toán động lực học dòng chảy
• Xác định phân bố lưu lượng và nhiệt độ trong một dòng chảy.
• Mô phỏng dòng chảy tầng và dòng chảy rối, dòng nén được và dòng không nén được, và nhiều dòng chảy kết hợp.
• Ứng dụng cho hàng không vũ trụ, đóng gói điện tử, thiết kế ôtô.
• Các đại lượng đặc trưng đáng quan tâm là vận tốc, áp suất, nhiệt độ và các hệ số màng.
Hình 3.6 – Trường dòng chảy trong ống dẫn và phân bố áp suất của thùng trộn
Âm thanh
ANSYS Acoustics Structures kết hợp cơ học kết cấu của ANSYS với một thư viện phần tử âm học vô hạn và hữu hạn mạnh mẽ, nhanh chóng, và đáng tin cậy. Việc tích hợp công nghệ ACTRAN trong môi trường ANSYS Workbench giúp cho việc sử dụng được dễ dàng hơn rất nhiều. Điều này có nghĩa là các kỹ sư có thể tập trung vào phân tích kết quả thay vì tốn thời gian thiết lập mô hình.
Hình 3.7 – Đồ thị áp suất mức áp âm
Phân tích tương tác giữa các trường vật lí
Xem xét sự tương tác giữa hai hoặc nhiều trường khác nhau. Vì trên thực tế các trường đều phụ thuộc lẫn nhau, nên không thể giải quyết chúng một cách tách biệt, bởi vậy cần có một chương trình giải quyết đồng thời cả hai hiện tượng bằng cách kết hợp chúng.
• Phân tích nhiệt-ứng suất.
• Phân tích áp điện (điện và kết cấu)
• Âm thanh (dòng chảy và kết cấu)
• Phân tích nhiệt - điện
• Cảm ứng nhiệt (từ và nhiệt)
• Phân tích tĩnh điện - kết cấu
2. Giải bài toán cơ học kết bằng phần mềm Ansys
Trong hệ thống tính toán đa năng của Ansys, bài toán cơ kỹ thuật được giải quyết bằng phương pháp phần tử hữu hạn lấy chuyển vị làm gốc. Trong bài toán kết cấu (Structural), phần mềm Ansys dùng để giải các bài toántrường ứng suất – biến dạng, trường nhiệt cho các kết cấu. Giải các bài toán dạng tĩnh, dao động, cộng hưởng, bài toán ổn định, bài toán va đập, bài toán tiếp xúc. Các bài toán được giải cho các dạng phần tử kết cấu thanh, dầm, 2D và 3D, giải các bài toán với vật liệu đàn hồi: đàn hồi phi tuyến, đàn dẻo lý tưởng, dẻo nhớt, đàn nhớt…Ansys cung cấp trên 200 kiểu phần tử khác nhau. Mỗi kiểu phần tử tương ứng với một dạng bài toán. Khi chọn một phần tử, bộ lọc sẽ chọn các module tính toán phù hợp, và đưa ra các yêu cầu về việc nhập các tham số tương ứng để giải. Đồng thời việc chọn phần tử, Ansys yêu cầu chọn dạng bài toán riêng cho từng phần tử. Việc tính toán còn phụ thuộc vào dạng vật liệu. Mỗi bài toán cần đưa mô hình vật liệu, cần xác định rõ mô hình là vật liệu đàn hồi hay dẻo, là vật liệu tuyến tính hay phi tuyến tính, với mỗi vật liệu cần nhập đủ thông số vật lý của vật liệu. Ansys là phần mềm giải các bài toán bằng phương pháp số, chúng giải trên mô hình học thực. Vì vậy cần đưa vào mô hình học đúng. Ansys cho phép xây dựng các mô hình học 2D và 3D với các kích thước thực hình dáng đơn giản hóa hoặc mô hình như vật thật. Ansys có khả năng mô phỏng theo mô hình học với các điểm, đường, diện tích và mô hình phần tử hữu hạn với các nút và phần tử. Hai dạng mô hình được trao đổi và thống nhất với nhau để tính toán. Ansys là phần mềm giải bài toán bằng phương pháp phần tử hữu hạn, nên sau khi dựng mô hình hình học, Ansys cho phép chia lưới phần tử do người sử dụng hoặc tự động chia lưới. Số lượng nút và phần tử quyết định đến độ chính xác của bài toán, nên cần chia lưới càng nhỏ càng tốt. Nhưng việc chia lưới phụ thuộc năng lực của từng phần mềm.
Để giải một bài toán bằng phần mềm Ansys, cần đưa các điều kiện ban đầu và điều kiện biên cho mô hình hình học. Các ràng buộc, các nội lực hoặc ngoại lực (lực, chuyển vị, nhiệt độ, mật độ) được đưa vào tại từng nút, từng phần tử trong mô hình hình học.
Sau khi xác lập được các điều kiện bài toán, để giải chúng Ansys cho phép chọn các dạng bài toán. Khi giải các bài toán phi tuyến, đặt ra vấn đề là sự hội tụ của bài toán. Ansys cho phép xác lập các bước lặp để giải bài toán lặp với độ chính xác cao. Để theo dõi bước tính, Ansys cho biểu đồ các bước lặp và hội tụ. Các kết quả tính toán được ghi vào file dữ liệu. Việc xuất các dữ liệu được tính toán và lưu trữ, Ansys xử lý rất mạnh, cho phép xuất dữ liệu dưới dạng đồ thị, ảnh đồ, để có thể quan sát trường ứng suất và biến dạng, đồng thời cũng cho phép xuất kết quả dưới dạng bảng số.
2.1 Các bước phân tích của bài toán kết cấu bằng phần mềm Ansys
a) Pre-processing: tiền xử lý, bao gồm:
• Xây dựng mô hình (Model generation)
• Định nghĩa loại phân tử (Define element type)
• Định nghĩa vật liệu (Define material)
• Chia lưới – tạo mô hình phần tử hữu hạn (Meshing)
b) Solution: Giải bài toán, bao gồm
• Xác định loại phân tích (Analysis type)
• Áp đặt tải trọng và xác định điều kiện biên (Specify loads and boundary conditions) (adsbygoogle = window.adsbygoogle || []).push({});
• Giải bài toán (solver)
c) Post-processing: hậu xử lý, bao gồm:
• Phân tích hình ảnh hoặc đưa ra kết quả
• Kiểm tra kết quả
Ngoài 3 bước chính trên, quá trình phân tích bài toán trong ANSYS còn phải kể đến quá trình chuẩn bị (preferences) chính là quá trình định hướng cho bài tính. Trong quá trình này cần định hướng xem bài toán ta sắp giải dùng kiểu phân tích nào (kết cấu, nhiệt hay điện từ…), mô hình bài toán như thế nào (đối xứng trục hay đối xứng quay, mô hình 3 chiều đầy đủ…), dùng kiểu phần tử nào (Beam, Shell, Struss,…)
2.2 Hai phương pháp làm việc với Ansys
Ansys cung cấp cho người dùng hai phương pháp làm việc: dùng GUI (Graphical User Interface) với ưu điểm là trực quan, dễ thao tác và phương pháp dùng câu lệnh APDL (Ansys Parametric Design Language), đây là phương pháp làm việc mà đòi hỏi
người sử dụng phải hiểu rõ về cấu trúc của bài toán, cũng như danh sách các lệnh và cấu trúc của nó. Nhưng làm việc với ansys bằng phương pháp dùng câu lệnh có những ưu điểm sau:
• Người dùng có thể tạo ra Input file, thuận tiện cho việc thay đổi giá trị của các tham số (kết hợp với các lệnh vòng lặp, các lệnh có điều kiện), trong khi đó trong chế độ GUI thì khi thay đổi giá trị tham số chúng ta phải tiến hành phân tích lại bài toán từ đầu.
• Tạo ra các macro file và sử dụng chúng như một hàm tự định nghĩa
3. Ví dụ
Cho một hệ khung giàn có kết cấu như hình 4. Xác định chuyển vị tại các gối, phản lực và nội lực của hệ . Biết moodun đàn hồi E=200 GPa, diện tích mặt cắt ngang A= 3250 mm2 .
Hình 3.8 – Kết cấu khung giàn
Với bài toán kết cấu này, ta có thể tiến hành giải theo hai phương pháp đó là phương pháp GUI và phương pháp dùng câu lệnh của ngôn ngữ APDL.
Để giải bài toán, chúng ta cũng đi xây dựng các bước như đã nói ở phần 2.1 bao gồm:
• Xây dựng mô hình của kết cấu:
- Dùng lệnh tạo keypoint trong Ansys để tạo ra các điểm đặc biệt của kết cấu. - Dùng lệnh tạo đường thẳng (Line) để nối các keypoint đó thành hình dạng của
1 1 2 1 3 2 3 2 4 3 4 3 5 4 5 4 6 5 6 5 7 6 7 X Y Z
Bridge Truss Tutorial MAY 18 2012
00:02:57 LINES
TYPE NUM
Hình 3.9 – Mô hình của giàn khi chưa chia lưới
• Khai báo phần tử và thuộc tính của vật liệu
- Định nghĩa kiểu phần tử của kết cấu. Ở ví dụ này kiểu phần tử của kết cấu là kiểu thanh (LINK)
- Khai báo thuộc tính vật liệu: Khai báo hệ số mô đun đàn hồi và diện tích mặt cắt ngang của thanh
• Chia lưới phần tử
- Định kích thước phần tử - Chia lưới
1 1 2 1 3 2 3 2 4 3 4 3 5 4 5 4 6 5 6 5 7 6 7 X Y Z
Bridge Truss Tutorial MAY 18 2012
00:11:33 LINES (adsbygoogle = window.adsbygoogle || []).push({});
TYPE NUM
Hình 3.10 – Mô hình kết cấu khi đã được phần tử hóa
• Xác định điều kiện biên và đặt tải trọng
1 1 2 3 4 5 6 7 8 9 10 11 X Y Z
Bridge Truss Tutorial 1 1 2 1 3 2 3 2 4 3 4 3 5 4 5 4 6 5 6 5 7 6 7 X Y Z
Bridge Truss Tutorial 1 1 2 3 4 5 6 7 X Y Z
Bridge Truss Tutorial 1
X
Y Z
Bridge Truss Tutorial MAY 18 2012
00:16:03 E-L-K-N
U
F
Hình 3.11 – Mô hình của kết cấu khi đã đặt điều kiện biên và tải trọng
• Xuất kết quả
Hình 3.12 – Bảng kết quả phản lực tại các gối
1
X
Y Z
Bridge Truss Tutorial MAY 12 2012 21:25:03 DISPLACEMENT STEP=1 SUB =1 TIME=1 DMX =7.409
Hình 3.14 – Biểu đồ chuyển vị của các nút trong kết cấu
CHƯƠNG 4: TÍNH TOÁN KHẢ NĂNG PHÁ HỦY CỦA MỘT KẾT CẤU HAI VẬT LIỆU (BI – MATERIAL)
1. Phương pháp phân tích phá hủy[8]
Để tính toán khả năng phá hủy của vật liệu khi xuất hiện vết nứt, người ta tiến hành đi phân tích vết nứt để tính toán hệ số cường độ ứng suất K hay tỉ lệ giải phóng năng lượng G, J của vết nứt. Sau đó so sánh với các hệ số chống lại sự phá hủy tương ứng của vật liệu (Kc, Gc, Jc) để có thể dự đoán được sự phát triển của vết nứt.
Để tính toán được các giá trị đó, người ta thường dùng các phương pháp sau:
1.1 Phương pháp thực nghiệm
Phương pháp này được áp dụng để đo giới hạn phá hủy của các vật liệu. Hiện nay theo tiêu chuẩn ASTM E399, có 4 mô hình được dùng để đo giới hạn phá hủy
Hình 4.1 – Các mô hình thực nghiệm đo giới hạn phá hủy
Hình 4.2 – Mô hình thực nghiệm đo KIC
Ban đầu sẽ tiến hành đo lấy số liệu để xây dựng đường cong biểu diễn chuyển vị theo tải trọng. Chuyển vị được đo bằng cách sử dụng thiết bị Strain gages. Strain gagesđược gắn trên hai thanh kim loại mảnh. Để chắc chắn là đầu còn lại của thanh kim loại tự do xoay thì đầu này phải tiếp xúc với một cạnh sắc nhọn. Hai dụng cụ có cạnh sắc nhọn sẽ được gắn tại miệng vết nứt. Khi miệng vết nứt chuyển vị dẫn đến chuyển vị của thanh kim loại mảnh và điện thế đầu ra sẽ thay đổi do biến dạng thay đổi trong thiết bị strain gages. Qua các công thức trung gian, ta sẽ tính được chuyển vị. (adsbygoogle = window.adsbygoogle || []).push({});
Mẫu vật dùng để đo KIC phải tuân theo các quy định về kích thước như là chiều rộng W phải gấp 2 lần bề dày B. Người thí nghiệm sẽ tạo một vết nứt có sẵn trên mẫu vật và tỷ lệ chiều dài vết nứt / chiều rộng mẫu vật (a/W) phải nằm trong khoảng 0,45 và 0,55. Một điều vô cùng quan trọng là tất cả các kích thước đều phải rất lớn so với vùng chảy dẻo. Quá trình đo KIC thường xảy ra kết quả sai mà lỗi không phải do thao tác thí nghiệm. Đó là khi kích thước vùng chảy dẻo quá lớn thì sẽ không bao giờ đạt được giá trị đúng cho dù người thí nghiệm có kỹ thuật cao.
Hình 4.3 – Thiết bị dùng để đo KIC
1.2 Phương pháp tương quan chuyển vị (Displacement Correlation Methods)
Đây là một phương pháp áp dụng phần tử hữu hạn để tính hệ số K của tổ chức vết nứt. Nội dung cơ bản của phương pháp này là: chuyển vị tại một điểm của phần tử bên trong lưới được thế vào biểu thức giải tích tính chuyển vị gần đỉnh vết nứt, sau khi trừ đi chuyển vị của đỉnh vết nứt. Điểm được chọn nằm bên trên bề mặt nứt sao cho chuyển vị là lớn nhất.
Khi đó hệ số cường độ ứng suất K được tính theo công thức sau:
2 ( ) (2 2 ) yb ya I u u K r v µ π − = − (4.1) 2 ( ) (2 2 ) xb xa II u u K r v µ π − = − (4.2) 2 ( ) 2 zb za III u u K r µ π − = (4.3)
Hình 4.4 – Phương pháp tương quan chuyển vị
Ở đây b là điểm được chọn là điểm tương quan, a là đỉnh vết nứt, µ là mô đun đàn
Trong phương pháp này, độ chính xác phụ thuộc vào việc lựa chon điểm tương quan. Điểm tương quan được lựa chọn sao cho nằm trong vùng có sự ảnh hưởng lớn của hệ số K. Ta cũng có thể áp dụng phương pháp là lựa chọn một loại các điểm tiến