Giáo trình ứng dụng ansys vào bài toán kỹ thuật phần 1

LỜI NÓI ĐÀU Phương pháp phần tử hữu hạn FEM từ lâu đã trở thành một trong những phương pháp cơ bản và chính yếu cho việc tính toán và mô phỏng các trạng thái vật lý xảy ra trong các hệ t

Trang 1

ĐẠI HỌC QUỐC GIA TP HO CHi MINH

Trang 2

LỜI NÓI ĐÀU

Phương pháp phần tử hữu hạn (FEM) từ lâu đã trở thành một trong những phương pháp cơ bản và chính yếu cho việc tính toán và mô phỏng

các trạng thái vật lý xảy ra trong các hệ thống kỹ thuật Đồng thời trong

công nghiệp, trong các công trình nghiên cứu tại các trường đại học và

ngay cả trong nghiên cứu quân sự ở nhiều nước trên thế giới, phương

pháp phần tử hữu hạn đã và được xem như là một công cụ hiệu quả nhất Với việc ra đời của phương pháp phần tử hữu hạn, nhiều phần mềm tính toán (dựa trên phương pháp phần tử hữu hạn) cũng đã ra đời sau đó Trong đó phần mềm ANSYS là một trong những phần mềm lớn và rất mạnh trong việc tính toán ANSYS có thê giải được rât nhiều bài toán

trong nhiều lĩnh vực như cơ học kết cấu, cơ sinh học, cơ, nhiệt, truyền

nhiệt, cơ lưu chất Nhằm giúp ích cho việc học tập, thiết kế, tính toán,

mô phỏng và có thể là tối ưu hóa các bài toán kỹ thuật cũng như hiểu TỐ hon vé phan mềm ANSYS nên chúng tôi đã biên soạn quyên sách “Ứng dung Ansys vào bài toán kỹ thuật” Quyền sách này sẽ giúp cho bạn đọc

có thể dùng phần mềm ANSYS giải được các bài toán cơ bản trong nhiều lĩnh vực của kỹ thuật và từ đó có thể nghiên cứu và giải quyết được những bài toán lớn hơn Trong quá trình biên soạn không thể tránh khỏi những lỗi và sai sót, chúng tôi mong nhận được nhiều ý kiến đóng góp

của các bạn đọc

Mọi ý kiến đóng góp xin vui lòng liên hệ: Phan Đức Huynh, Bộ môn Cơ Học — Khoa Xây dựng & Cơ học Ứng dụng — Trường Đại học

Sư phạm Kỹ thuật TP Hồ Chí Minh

01 Võ Văn Ngân, Phường Linh Trung, Quận Thủ Đức, TP Hồ Chí Minh

Email: huynhpd@hcmute.edu.vn

Các tác giả

Trang 3

MỤC LỤC

CHUONG 1: SƠ LƯỢC VẺ PHƯƠNG PHÁP PHÀN TỬ

HỮU HẠN ‹cserrtitekreeEttrrtirrrrirocrc 9 _ L Các bước giải bài toán FEM

CHUONG 4: PHAN TICH DAO DONG

I Tìm đáp ứng chuyển vị của hệ dao động

Trang 4

CHUONG 5: UNG DUNG ANSYS VAO BAI TOÁN UNG

I Nguyên lý ST Venant

IL Anh hưởng của lỗ tròn đối với sự phân bố ứng suất và

II Tập trung ứng suất do lỗ elip

IV Ứng suất cho bài toán tiếp xúc hai chiều

CHUONG 6: UNG DUNG ANSYS VÀO BÀI TOAN CO

HOC PHA HUY

bố đều,

II Phân tích phá hủy của

CHUONG 7: UNG DUNG ANSYS VÀO BÀI TOÁN LƯU

II Giải bài toán dòng chảy trong kênh có van bướm 2 17

II Phân tích dòng khí xả đi qua một bộ giảm thanh (muffler) 227

IV Phân tích dòng chảy của gió khi đi qua dãy núi Rocky 238

TRUYEN NEUUBID sscsctssdsnjcsccsestnstonsnsobbcsossessocsbettltoets

I Đặc điểm chung của bài toán truyền HHIỆ sc0622016066i66x

IV Bài toán dẫn nhiệt ba c|

Phân tích cơ nhiệt của một thiết bị điện : -+

VI, ` Mỗi hắn gilp THÔI có cách g 2h ga À g2 ghaghÌ HH danh

Trang 5

CHUONG 9: ỨNG DỤNG ANSYS VÀO BÀI TOAN TIE

XUC GIU'A CAC CHI TIET MAY

I Lap phép có độ dôi của chốt,

1 Tiếp xúc giữa trụ tròn và hai khải vật thể -.sc5:ssccsse 299

II Tiếp xúc giữa bánh xe & đường ray co 300) CHƯƠNG 10: ỨNG DỤNG ANSYS VÀO BÀI TOÁN TÓI

ƯU HÓA

I Tối ưu hóa bài toán dầm 2D

Il Tối ưu hóa bài toán thanh - dầm 2D

CHƯƠNG 11: ỨNG DỤNG ANSYS VÀO BÀI TOÁN

TƯƠNG TÁC LƯU CHÁT-KÉT CÁU 357

I Luu chat chảy qua một kênh có vật cản sac 357

II Lưu chất chảy qua một kênh thẳng đứng có vật cản 370

CHƯƠNG 12: SỬ DỤNG LỆNH TRONG ANSYS 385

I Những lệnh ANSYS căn bản 225-cccccsccrrxee 385

390

+394 400

Il Bài toán dầm cong

II Bài toán biến dạng phẳng

IV Bài toán lưu chất

Trang 6

mô tả và biên khảo sát rât phức tạp, do đó việc tìm lời giải chính xác cho bài toán gặp rất nhiều khó khăn, mà đôi lúc không thể giải được Để khắc phục những khó khăn này, nhiều phương pháp sô đã được ra đời Trong những phương pháp số đó, nôi bật và thường dùng nhất là phương pháp phan tir hữu hạn (finite element method) Ý tưởng của phương pháp phần

tử hữu hạn là chia vật thể ra thành một tập hữu hạn các miền con (gọi là phần tử) liền nhau nhưng không liên kết hoàn toàn với nhau trên khắp từng mặt biên của chúng Trường chuyển vị, biến dạng, ứng suất được xác định trên từng phần tử Dạng phân tử có thể là thanh, dầm, giàn, tắm, vỏ, khối

Các phần tử nói kết với nhau qua các nút

I CÁC BƯỚC GIẢI BÀI TOÁN FEM

1 Rời rạc hoá miền khảo sát (tạo lưới cho phần tử): chia vật Tra thành nhiều phan | tử sao cho tinh chat vật lí của mỗi phần tử không thay đôi Nếu vật có biên dạng phức tạp thì chia các phần

tử ở gân biên sao cho thật nhuyễn

Xây dựng các phương trình phần tử

Lắp ghép các phương trình phần tử và tìm ra ma trận độ cứng Khử các điều kiện biên

Giải hệ phương trình toàn cục

am

Tính toán các kết quả

I VÍ DỤ 1

Xét một hệ giàn như hình 1.1 Tat cả các thanh đều làm bằng gỗ

với môđun đàn hồi Young E=l 9x10° (Ncm? ) và diện tích mặt cắt

ngang là 8 cm” Tính ứng suất, chuyên vị trong mỗi phần tử thanh giàn

và phản lực tại hai gối

Trang 7

Hình 1.1 a) Mô hình bài toán, b) Mối liên hệ gữa tọa độ tòan cục và hệ

tọa độ địa phương

Giải bài toán này sẽ trải qua ba giai đoạn: tiền xử lý, giải, hậu xử lý

1) Giai đoạn tiền xử lý

a) Nút và phần tử

Xét mỗi thanh giàn là một phần tử, điểm nối chúng lại với nhau là một nút Vì vậy, hệ giàn đã cho có thể chỉ thành 5 nút và 6 phần tử như biểu diễn ở hình 1.1 Sau đây là bảng liệt kê nút và phần tử của hệ giàn:

Phần tử (1), @), (4) và (6) có cùng chiều dài, diện tích mặt cắt ngang, môđun đàn hồi, nên hệ số độ cứng tương đương cho mỗi phần tử

10

Trang 8

[KJ?=k sin Øcos Ø sin? Ø —sinØcosØ -—sin°Ø

—sinØcosØ -—sin°Ø sin ØcosØ sin? Ø

=4/22xig5| Sin0cos0 sin?0 —sin0cos0 —sin?0

—cos?0 —sin0cos0 cos°0 sỉn 0cos0

—sin0cos0 —sin?0 sin0cos0 sin?0 ” Phần tử thứ (1)

Trang 10

cos* (90) sin (90) cos (90) sin (90)cos (90) sin? (90)

—cos? (90) —sin (90)cos (90)

—sin (90)cos(90) —sin? (90)

(1.10)

13

Trang 11

Địa phương: [K]f Poem

Trang 12

—cos? (135) —sin (135)cos(135)

—sin (135) cos (135) —sin (135) (1.13) cos? (135) sin (135)cos(135) sin (135) cos (135) sin? (135)

Trang 14

e) Áp đặt điều kiện biên và tải trọng

Điều kiện biên: nút 1 và 3 bị kẹp chặt nên ¡y= Uy= Uax= U›y= 0

Tải trọng: ngoại lực tác dụng lên nút 4 và 5 theo phương thang

dimg: F4y = —500 N va Fsy=—500 N

Ủy

(1.20)

Từ những điều kiện biên thì phương trình độ cứng toàn cục có thể

khử đi hàng thứ nhất, thứ hai, thứ năm và thứ sáu; cùng với những cột tương ứng Từ đó phương trình độ cứng toàn cục chỉ còn ma trận 6x6 :

17

Trang 15

2) Giai đoạn giải

Giải phương trình (1.22) tìm được:

ie = 0.0035 cm U,, =-0.0102 cm fe = 0.00118 cm Uy =-0.0114cem — |U,, =-0.0195 cm ie = 0.00240 cm (1.22)

Ma trận chuyển vị toàn cục có thể viết là:

U,„|_ |-0.0035 U„ | |-0.0102

U,,| |000118 U„| |-00114

Us, | |0.00240 Us| |-0.0195 3) Giai đoạn hậu xử lý

Sau khi tìm được chuyển vị của các nút thì có thể ứng dụng để tìm

Trang 16

Ry

Ry Ryy

0 0 _| 0 0

149 -1.49

0 -422 -149 -149 -149 -1.49

0

0 -149 1.49 3.71 -1.49 4.22

Sử dụng các phép tính ma trận cho phương trình trên sẽ tìm được các phản lực:

Trang 17

b) Nội lực và ứng suất pháp

Các nội lực phần tử f, va Sic Va bằng nhau và ngược chiều:

dụ =k(u, —w„)

Fig = (typ tu)

Trong công thức (1.26), uix Va uj là chuyển vị tại hai đầu nút của

phân tử ở hệ địa phương

{0}=[T]'{U}

HH, = —sinØ cosØở 0 - 0 ||Uy„ (129)

Huy 0 0 —sinØ cosØ || Uy

_ Sau khi tính được nội lực, có thể tính được ứng suất pháp cho mỗi

phân tử:

AE

_f k{, —w„) LL (4 -u,) _ :[* ~ity ) (130)

Ví dụ, tính nội lực và ứng suất pháp cho phần tử thứ (5) Những dữ

liệu đã có của phần tử thứ (5): Ø=45°, U;y =-0.00355 cm,

Uy =-0.01026 cm, Us, = 0.0024 cm va U,, =—0.0195 cm Dau tiên,

tinh chuyén vj dia phuong cua nút 2 và 5:

20

Trang 18

tạ, cos4Š5 sin45 0 0 —0.00355

Us, 0 0 -—sin45 cos45 || ~0.01950

thực hiện phép tính ma trận ta tìm được ; =-0.00976m và

us, =—0.01209 m, thé vao công thức (1.26) và (1.30) sẽ có được nội lực

và ứng suất pháp cho phầmtử thứ (5) lần lược là 696N và §7 N/mỲ

Tương tự tìm được nội lực và ứng suât pháp cho các phân tử còn lại

Il Vi DU 2

Sử dụng phan tử tam giác tuyến tính để xác định chuyên vị của các nút và ứng suât trong các phần tử của một tắm mỏng được ràng buộc chuyển vị và chịu tác dụng của các lực bề mặt như hình 1 6 Ngoài ra tắm mỏng còn đặt trong môi trường nhiệt độ thay đổi A7 = 10°C Tam mong

có bê dày 0.5cm, môđun đàn hồi Young E=15x10°N/cm?, hé s6 Possion v=0.25 và hệ số giãn nở nhiệt = 6x10” re

600N/cm?

Hình L6

Chia tắm mỏng thành hai phần tử tam giác như biểu diễn ở hình

1.7 Tọa độ toàn cục của các nút được cho trong bảng I.2

21

Trang 19

c Ma trận thuộc tính vat ligu (material property matrix)

Giả sử bài toán thuộc bài toán ứng suất phẳng

Trang 20

& Ma trận hình dạng phần tir (element shape matrix)

yh =? -yW=y,-ya 3» [x eaP aa) =x,-2 =

Trang 21

e Ma tran d@ ciing phan tit (element stiffness matrix)

ø Vector lực nhiệt của mỗi phần tử (thermal load vector)

Bién dang ban dau (initial strain) do sự thay đổi từ nhiệt độ (AT = 10°C, a = 6x10"):

Vector lực nhiệt của mỗi phần tử:

Trang 22

€ 1 7

{p£}=:A°[B'"]Ï[D]{e'}=(-900 -300 900 -300 0 600}N aan {p?}=:^®[B° Ï'Ip]feF}={o -300 900 300 -900 O}N

h Biến đổi ma trận độ cứng, vector lực và vector lực nhiệt về dạng toàn cục

Trang 23

1 Cộng các ma trận độ cứng, vector lực và vector lực nhiệt

Từ đó phương trình hệ thống toan cuc (global system of equation) như sau:

[K]{U} ={P}

15 15 -69 3 0 0 =6 -l2] (Uy) [300v10-900)

l5 35 3 -19 0 0 -1§ -l6| |Uy ~300

+69 3 81 -15 -6 -I8 -6 30| |Uay 900 (1.55) 1ữ|-3 -19 -15 41 -12 -16 30 -l6| |Ưy| | 600

Trang 25

Chương 2

TONG QUAN VE PHAN MEM ANSYS

ANSYS là một phần mềm mạnh được phát triển và dp dung rong rãi trên toàn thế giới, bằng cách sử dụng phương pháp phần tử hữu hạn

(FEM - finite element method) trong tính toán Những ưu điểm nổi bậc

của ANSYS như sau:

ø - Giảm chỉ phí sản xuất vì có thể tính toán thử nghiệm

I KHOI DONG ANSYS

Như các phần mềm khác thì khởi động ANSYS bằng các cách sau:

Trên thanh Taskbar chon: Start / Programs / ANSYS / ANSYS Product Launcher

s* Click đúp chuột vao biéu tuong “ANSYS Product Launcher” trén

desktop (nêu có)

Sau đó màn hình sẽ xuất hién ctra sé “Launcher”, ké tiếp chọn mục

“File manegement” Trong cửa số mới chọn đường dẫn cho ƒïie cần lưu

gitt (working directory) va dat tén cho bai toan (job name) Sau khi chon xong chon “Run” thi chuong trinh ANSYS sẽ được khởi động

II GIAO DIEN ANSYS

Sau khi khởi động chương trình ANSYS thì màn hình sẽ hiện giao

diện sau:

29

Trang 26

Hình 2.1 Giao diện ANSYS

[A]: MENU TIEN [CH (Utility Menu)

List (liệt kê)

Plot (hiển thị dạng đồ hoạ)

PlotCtrls (lựa chọn trước khi lựa chọn dạng đồ hoạ)

WorkPlane (mat phẳng làm việc)

Parameters (quản lý các thông số)

Macro (tạo các dạng chuẩn)

MenuCtrls (chọn lựa cách hiển thị các Menu)

Help (trợ giúp)

[B]: NHAP DU’ LIEU BANG LENH (Input)

Sử dụng cho hướng giao tiếp bằng lệnh

[C]: THANH CÔNG CỤ (Toolbar)

“Thường xuyên sử dụng cho hướng giao tiếp bằng chuột

30

Trang 27

[D]: MENU CHÍNH (Main menu)

Thường xuyên sử dụng trong ANSYS cho hướng giao tiếp bằng chuột

[E]: VÙNG HIẾN THỊ (Graphics)

Hiển thị các kết quả bằng đồ hoạ

Phóng to thu nhỏ đối tượng trên Graphics: cũng như các phần mềm khác thì ANSYS cũng cung cập các công cụ đề có thể di chuyển khung nhìn theo ý muốn Có thê phóng to thu nhỏ hoặc xoay khung nhìn bằng thanh công cụ bên tay phải màn hình hoặc có thể thực hiện theo hộp thoại

băng cách thực hiện theo đường dẫn:

COMMAND: Utility Menu— PlotCtrls— Pan Zoom Rotate Xuất hiện hộp thoại Pan-Zoom-Rorare như biêu diễn ở hình 2.2

TT Dynanie Mode ˆ [F]: Hiển thị mô

Oz] =] o>}

Hinh 2.2

epee

Trang 28

II LƯU TRỮ VÀ HÒI PHỤC

COMMAND: Uriliy Menu > File — Save as— dat tén file, chon

thư mục —>ØK (dữ liệu lưu với phân mở rộng 1a *.db)

|ANSYS Multiphysics U@l Save DataBase

fle Select Uist Bot Plo Save Database lo Lx_—]

Clear & Start New m= con vd lke

List Files of Type:

Database Files ("db") _v fan

Hinh 2.3

Muốn nở lại bài toán đã tính:

COMMAND: Utility Menu File > Resume from Sau d6 chon file cần mở, sau d6 bam OK

Resume Database

Resume Database Fiom Directories

c:\documents and settings\t

List Files of Type: Dives

I Ignore parameters on db fle I Suppress auto plotting

Hinh 2.4

Iv CAU TRUC BAI TOAN TRONG ANSYS

Khi tính toán bất kỳ một bài toán ANSYS nào thì cũng qua các công đoạn sau:

1 Tạo mô hình tính - Preprocessing: gồm các bước:

() Dinh nghĩa điểm chính/đường/diện tíchhể tích

32

Trang 29

(keypoints/lines/areas/volumes),

(ii) Dinh nghia kiéu phan tir (element type) va tinh chất vật

ligu/hinh hoc (material/geometric)

(ii) Tạo lưới đường/diện tich/thé tích (lines/areas/volumes) Mức độ chỉ tiết phụ thuộc số lượng chiều của bài toán, ví dụ 1D, 2D, bất đối xứng, 3D

2 Tinh todn - Solution: đặt tải trọng (lực điểm hoặc áp suất ), các

ràng buộc (độ dịch chuyên, góc xoay) và tính toán

3 Xử lý kết quả - Posiprocessing: các bước bổ sung và xem kết quả Trong phân này có thể bao gồm:

Giải bài toán cho ở ví dụ về giàn đã được cho ở chương 1 bằng ANSYS

Xét một hệ giàn như hình 2.5 Tất cả các thanh đều làm bằng gỗ ai

médun dan hồi E=1.9x10° (Niem* ) va dign tich mat cat ngang 1a 8 cm’, Tinh ứng suất, chuyển vị trong mỗi phần tử thanh giàn và phản lực tại hai gối

Trang 30

ANSYS đưa ra hai kiểu phần tử để phân tích bài toán về giàn là: LINKI và LINKS Ò đây đề giải bài toán trên ta dùng phan tir LINK1

Phần tử LINKI là phần tử đối xứng trục, có hai nút @ và j), có hai

bậc tự do (Uy và Uy ), như hình 2.6 Để tìm hiểu thêm về tham số của

phần tử này, ta có thể vao “Help” cla ANSYS

° Nhấp vào Ađđ làm xuất hién cira s6 “Library of Element Types”

để chọn các phần tử cần sử dụng

Nhấp chọn Link rồi chọn 2D spar 1, sau đó chon Apply

e — Cửa số “Elemenr Types” xuất hiện dạng phần từ thứ nhất là LINK

¢ Bam Close

2 Nhập diện tích mặt cắt ngang cho từng phần tử giàn

COMMAND: Main Menu— Preprocessor — Real Constans >

Add/ Edit/Delete

e Nhấp vào Add lam xuất hiện cửa số “Element Type for Real Constans”, nhdp chọn phần tử LINK1 roi bam OK

e Cửa số “Real Constans for LINKI" xuất hiện, nhập thông số: AREA = § (lưu ý rằng trong ANSYS không nhập đơn vị mà thay vào đó ANSYS ngâm hiểu đơn vị đã được đồng nhất)

Trang 31

3 Nhập tính chất vật liệu

COMMAND: Main Menu— Preprocessor —> Material Props > Material Models

° Nhấp d6i vao Structural, Linear, Elastic, va Isotropic lam xuất

hiện hộp thoai “Linear Isotropic Properties forMaterials Number

hiện hình mũi tên hướng lên (Ì )

35

Trang 32

Y Element (1): nhấp node 1 va node 2—>Apply

Element (2): nhấp node 2 va node 3— Apply

Element (3): nhap node 3 va node 4—> Apply

v

Y Element (4): nhap node 2 va node 4— Apply

Vv Element (5): nhap node 2 va node 5—> Apply

v Element (6): nhap node 4 va node 5— OK

5 Áp đặt điều kiện biên và tải trọng

a) Điều kiện biên về chuyển vị

COMMAND: Main Menu — Solhdion —> Define Loads => Apply

— Structural + Displacement — On Nodes

cũng xuất hiện hình mũi tên hướng lên (Ï )

«_ Nhấp chọn node 1 và node 3, rồi bắm ØK, cửa số “Apply U ROT

on Nodes” khac xuat hién

¢ Chon ALL DOF trong 6 Lab2 và nhập giá trị “0” vào ô

Displacement value

36

Trang 33

bj Lực tập trung — -

COMMAND: Main Menu — Solution > Define Loads + Apply

— Structural — Force/Moment — On Nodes

cũng xuất hiện hình mũi tên hướng lên ( Ï )

® - Nhấp chọn nøde 4 và node 5 rồi bam OK, thi cita s6 “Apply F/M

øn Nodes” khác xuất hiện

¢ Chon FY trong 6 4b và nhập giá trị “-500” vào ô VALUE, rồi nhap OK

ELEMENTS ELEM NUM

Hinh 2.8

6 Giải bài toán

COMMAND: Main Menu —> Solution —> Solve —> Cuzrent LS

7 Xử lý kết quả

COMMAND: Main Menu — General Postproc > Plot Results

— Deformed shape — chon muc Def + undeformed

37

Trang 34

[\ Piot Deformed Shape

[PLDISP] Plot Deformed Shape

KUND Items to be plotted

COMMAND: Mai Menu —> General Postproc —> Plot Resulis

— Contour Plot — Nodal Solution

Deformed shape with undeformed edge trong 6 Undisplaced shape

key R6i bam OK, thi trudng chuyén vi duge cho ở hình 2.10

nút

38

Trang 35

PRNSOL Command

File

PRINT U NODAL SOLUTION PER NODE

wwe POST1 NODAL DEGREE OF FREEDOM LISTING »*+~

LOAD STEP= 1 SUBSTEP=

TIME= 1.8888 LOAD case 8

THE Ree DEGREE OF FREEDOM RESULTS ARE IN GLOBAL

So sánh với kết quả tính ở chương 1

U;y =-0.0035 m U,, =0.00118 m U,, = 0.00240 m

Uy =-0.0102 m Uy =-0.0114m Us, =-0.0195 m

Ta thấy rằng kết quả giải bằng ANSYS phù hợp với kết quả tính ở

chương 1

b) Phản lực

COMMAND: Utility Menu—> List Results — Reaction solution

Kết quả phản lực giải bằng ANSYS cũng phù hợp với kết quả tính

ở chương I

[ÏPRRSOL Command

PRINT REACTION SOLUTIONS PER NODE

suns POST1 TOTAL REACTION SOLUTION LISTING

xe

LOAD STEP= 1 SUBSTEP= 1

THE FOLLOWING X.¥.Z SOLUTIONS ARE IN GLOBAL COORDINATES

Trang 36

e) Trường ứng suất và nội lực

Vào “Help” để xem đầu ra của phần tử LIVKI

Output Quantity Name ETADLE and ESOL Command Input

Hinh 2.13 Lấy nội lực:

COMMAND: Main Menu — General Postproc — Element Table

— Define Table

Xuất hiện cửa s6 “Element Table Data’

Nhấp vào Ađ¿ làm xuất hiện cửa sé “Define Additional

Element Table Items”

Nhập “AXFORCE” vào ô Lab

Trong phần Irem, Comp kéo thanh cuộn xuống dé chọn mục

By sequence num

Trong 6 ké bén chon SMISC

Trong ô ở dưới thêm số “1” vao nhu sau: “SMISC, 1”, ri bam OK

Lấy ứng suất:

COMMAND: Main Menu — General Postproc — Element

Table — Define Table

40

Xuất hiện cửa s6 “Element Table Data”.

Trang 37

¢ Nhap vào Add làm xuất hiện cửa số “Define Addirional Element Table Items”

e Nhap “AXSTRESS” vao 6 Lab

e Trong phần /em, Comp kéo thanh cuộn xuống để chọn mục

PRINT ELEMENT TABLE ITEMS PER ELEMENT

.~wxxe POST1 ELEMENT TABLE LISTING xxx»

STAT CURRENT CURRENT ELEM AXFORCE AXSTRESS

2 1414.2 176.278

3 SI

MINIMUM UALUES ELEM 1 1

MAXIMUM UALUES ELEM 2 2 VALUE 1414.2 176.78

Hinh 2.14 Xem két qua:

COMMAND: Main Menu — General Postproc — List Results

— Elem Table Data

¢ Chon AXFORCE va AXSTRESS trong 6 Lab1-9

e _ Rồi bấm ÓK Thì xuất hiện bang PRETAB nhv hinh 2.14

So sánh với kết quả nội lực và ứng suất trên phần tử thứ (5):ở chương | lan luge 1A 696 N va 87 N/m? Kết quả khá phù hợp với kết quả tính ở chương 1

41

Trang 38

8 Giải bằng lệnh trong ANSYS

Trang 39

Giải bài toán ứng suất phẳng được cho ở ví dụ 2 của chương 1

Sử dụng phần tử tam giác tuyến tính và tứ giác tuyến tính để xác

định chuyển vị của các nút và ứng suất trong các phần tử của một tắm

mỏng được ràng buộc chuyền vị và chịu tác dụng của các lực bể mặt như

hình 2.15 Ngoài ra tim mỏng còn đặt trong môi trường nhiệt độ thay đổi

tir AT = 10°C đến nhiệt độ chuẩn Tắm mỏng có bề dày 0.5cm, médun

đàn hồi Young E = 15x10” N/cmỶ , hệ số Possion w= 0.25, hệ số giãn nở

Trang 40

1 Nhập tính chất vật liệu

COMMAND: Main Menu — Preprocessor + Material Props > Material Models

Cửa số “Define Material Model Behavior’ sé xuất hiện

Nhấp đôi vao Structural, Linear, Elastic va Isotropic làm xuat hién hộp thoai “Linear Isotropic Properties forMaterials Number 1, Nhập vào các giá trị EX = 15e6 và PRXY = 0.25

Trén cita s6 “Define Material Model Behavior”

Nhấp đôi vào Structural, Thermal Expansion, Secant Coefficient

va Isotropic lam xuất hiện hộp thoai “Linear Isotropic Properties forMaterials Number 1”

« _ Nhập vào các gid tri APLX = 6e-6, rồi bắm OK

2 Chọn kiểu phần tử

COMMAND: Main Menu — Preprocessor — Element Type > Add/ Edit/Delete

Nhấp vào Add lam xuat hign cia sé “Library of Element Types”

dé chọn các phân tử cân sử dung

Nhap chon phan tir Structural Mass — Solid

Chon kidu phan tir Quad 4node 42 rdi bam OK

Cia sé “Element Types” xuất hiện dạng phần từ PLANE42

Nhấp chon muc Option dé lam xuất hiện cửa số “PLANE42

element type options”

Chon điều kiện Plane strs w/thk trong hộp” thọai “Element

behavior” dé dinh trang thái cho phần tử (ứng suât phăng)

Định dạng
Số trang	179
Dung lượng	10,49 MB