Bài giảng phần mềm ANSYS cho các lớp cơ bản: các thao tác cơ bản, tính toán kết cấu giàn, dầm, khung, bài toán phẳng........................................................................................................................
BÀI 1: GIỚI THIỆU PHẦN MỀM ANSYS PHẦN MỀM PHÂN TÍCH MƠ PHỎNG KẾT CẤU CƠNG TRÌNH DƯỚI TÁC DỤNG CỦA TẢI TRỌNG TĨNH ĐỘNG LỰC, PHI TUYẾN THEO PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Hà Nội, 02/2012 SIMQKE_GR VỊ TRÍ GiD AutoCad HyperMesh CATIA SeismoSpect SolidWorks SAP2000 GEO Sigma/w ABAQUS HyperGraph SAMCEF ADINA PREP SOLU HyperView POST ANSYS FLAC3D OriginPro PLAXIS Phase2 THAM KHẢO CHÍNH ứng dụng rộng rãi ngành khoa học kỹ thuật: -Xây dựng -Địa chất khoáng sản đủ khả tiến hành nghiên cứu: -Kết cấu -Nhiệt -Chất lỏng -Điện từ -Âm thanh… -Thủy lợi -Đường sắt -Đường -Công nghiệp quốc phòng -Hàng khơng vũ trụ -Chế tạo máy -Dầu khí -Cơng nghiệp nhẹ Có 70% trường Đại học Viện nghiên cứu giới sử dụng phần mềm ANSYS giảng dạy nghiên cứu -Điện tử -Điện gia dụng -Y học Khởi động ANSYS Khởi động nhanh Sử dụng mặc định thiết lập môi trường làm việc lần đầu Kiến nghị sử dụng phương thức Thoát khỏi ANSYS Phương pháp 1: GUI: Utility Menu > File > Exit Phương pháp 2: GUI: Nhấn nút COMMAND: /EXIT công cụ ngang Phương pháp 3: GUI: Nhấn biểu tượng góc phải hình Lựa chọn hình thức lưu trữ trước khỏi Giao diện phần mềm Menu tiện ích Cửa sổ nhập lệnh Thanh cơng cụ Menu Thanh trạng thái Cửa sổ hiển thị PHƯƠNG PHÁP THAO TÁC GUI Nhập lệnh trực tiếp -Thao tác đơn giản, trực quan phù hợp với số đông người quen sử dụng thao tác window -Với tốn khơng gian phức tạp, cần giải lặp, tốn nhiều thời gian công sức chí khơng thể thực Sử dụng công cụ -Tạo file dung lượng lớn Thao tác menu Menu Khối xử lý số liệu Khối xử lý kết Khối tính tốn Menu tiện ích Lựa chọn đối tượng Quản lý file Hiển thị đối tượng theo dạng bảng liệu Định nghĩa hệ tọa độ mặt phẳng làm việc Định nghĩa tham số tính tốn Hiển thị đối tượng cửa sổ hình Khống chế hiển thị PHƯƠNG PHÁP THAO TÁC APDL Soạn thảo câu lệnh Notepad với đuôi *.txt APDL (ANSYS Parametric Design Language) Ngơn ngữ tham số hóa thiết kế APDL sử dụng ngôn ngữ FORTRAN nên phải tuân thủ quy tắc ngơn ngữ lập trình -Với mơ hình phức tạp -Với mơ hình cần chỉnh sửa -Khi cần giải lặp -Tạo module thông dụng -Dung lượng file số liệu nhỏ, dễ dàng trao đổi -Câu lệnh nhiều, khó nhớ -Cần hồn chỉnh câu lệnh file thu kết -Không trực quan, dễ phát sinh chán nản Quản lý file Tên file Jobname.db Loại hình ANSYS, Inc Jobname.elem Text Document Jobname.emat EMAT File Jobname.err Error log Jobname.esav ESAV File Nội dung File kho liệu, ghi chép liệu PT, nút, tải trọng File định nghĩa PT, ghi chép liệu định nghĩa PT File ma trận phần tử, ghi chép liệu ma trận phần tử File thông tin lỗi với cảnh báo, ghi chép tất thông tin cảnh báo, lỗi tính tốn File bảo lưu liệu PT, bảo lưu liệu tính tốn PT Jobname.log Text Document File nhật ký, ghi chép tất lịch trình nhập lệnh theo phương thức bổ sung Jobname.out Text Document File xuất ANSYS, ghi chép tình hình thi hành lệnh Jobname.rfl RFL File File kết phân tích FLOTRAN Jobname.rmg RMG FIle File kết phân tích từ trường Jobname.rst RST File File kết phân tích trường kết cấu ngẫu hợp Jobname.rsh RSH File File kết phân tích trường nhiệt độ Jobname.snn Text Document Jobname.stat STAT File File bước tải trọng, ghi chép thông tin tải trọng bước tải trọng File trạng thái xử lý vận hành ANSYS Chỉ định tên file Phương pháp 1: Để tên file mặc định: Jobname = file Phương pháp 2: Người dùng định nghĩa tên file: Jobname = DapBTTL Phương pháp 3: GUI: Utility Menu > File > Change Jobname COMMAND: /FILNAME,DapBTTL Sử dụng phương pháp này, tất liệu sau thay đổi chuyển đến tên file làm việc Nhưng Jobname.log Jobname.err lưu tên file làm việc định phương pháp thứ hai Lưu file liệu Phương pháp 1: GUI: Utility Menu > File > Save as Jobname.db Utility Menu > File > Save as Phương pháp 2: Jobname.db tương đương với Save as Đọc file liệu Phương pháp 1: GUI: Utility Menu > File > Resume Jobname.db Utility Menu > File > Resume from Phương pháp 2: Jobname.db tương đương với Resume Nên tập thói quen: -Mỗi hạng mục phân tích nên cho vào thư mục riêng -Sử dụng tên file làm việc khác giải vấn đề mới, định nghĩa tên file làm việc cửa sổ khởi động ANSYS -Trong q trình phân tích, sau khoảng thời gian định nên lưu liệu -Trước tiến hành thao tác mà chưa biết hậu (ví dụ chia mạng lưới phần tử) có khả tạo thành ảnh hưởng lớn (ví dụ xóa bỏ) nên lưu liệu -Sau hồn thành q trình phân tích nên lưu giữ tốt file: file log (Jobname.log), file liệu (Jobname.db), file kết (Jobname.rst, Jobname.rsh ), file bước tải trọng (Jobname.s01, Jobname.s02 ), file xuất (Jobname.out) file môi trường vật lý (Jobname.ph1, Jobname.ph2 ) 10 Số liệu đầu phần tử LINK1 S(AXL) Ứng suất theo phương dọc trục phần tử EPEL(AXL) Biến dạng đàn hồi theo phương dọc trục phần tử EPTH(AXL) Biến dạng nhiệt phương dọc trục phần tử EPIN(AXL) Biến dạng ban đầu theo phương dọc trục phần tử EPPL(AXL) Biến dạng dẻo theo phương dọc trục phần tử EPCR(AXL) Biến dạng từ biến theo phương dọc trục phần tử EPSW(AXL) Biến dạng dãn nở theo phương dọc trục phần tử Phần tử LINK8 Phần tử chiều chịu kéo nén dọc trục thường dùng để mô kết cấu giàn không gian Số liệu đầu vào phần tử LINK8 Tên gọi LINK8 Điểm nút I, J Độ tự UX, UY, UZ Hằng số thực AREA (diện tích mặt cắt), ISTRN (biến dạng ban đầu) Đặc trưng vật liệu EX, ALPX, DENS, DAMP Tải trọng bề mặt Khơng có Tải trọng khối Nhiệt độ: T(I), T(J) Đặc tính Dẻo, từ biến, dãn nở, ứng suất cứng hóa, biến dạng lớn, phần tử sinh chết 32 Số liệu đầu phần tử LINK8 SAXL Ứng suất theo phương dọc trục phần tử EPELAXL Biến dạng đàn hồi theo phương dọc trục phần tử EPTHAXL Biến dạng nhiệt theo phương dọc trục phần tử EPINAXL Biến dạng ban đầu theo phương dọc trục phần tử EPPLAXL Biến dạng dẻo theo phương dọc trục phần tử EPCRAXL Biến dạng từ biến theo phương dọc trục phần tử EPSWAXL Biến dạng dãn nở theo phương dọc trục phần tử Một vài lưu ý tính tốn nội lực ứng suất giàn: * Định nghĩa đặc trưng hình học giàn: Prep > Real Constants > Add\Edit\Delete * Nội lực ứng suất giàn: General Postprocessor > Element Table > Define Table ETABLE, NAXL, SMISC,1 ETABLE, SAXL, LS, 33 * Hiển thị biểu đồ ứng suất: General Postprocessor > Plot Result > Contour Plot > Line Element Result * Xuất bảng liệt kê giá trị lực dọc ứng suất: General Postprocessor > List Result > Elem Table Data Ví dụ - Giàn phẳng hình thang Xác định chuyển vị, lực dọc ứng suất giàn giàn phẳng có sơ đồ tính tốn hình đây, bỏ qua trọng lượng thân giàn Tiết diện mặt cắt ngang giàn A=32cm2 Vật liệu thép CT3 có E=2.1x108kN/m2, =0.3 Sơ đồ tính tốn giàn 34 Ví dụ - Giàn khơng gian hình thang ngược Xác định chuyển vị, lực dọc, ứng suất phản lực liên kết giàn khơng gian có sơ đồ tính tốn hình Tiết diện mặt cắt ngang giàn A=0.0032m2 Vật liệu thép CT3 có E=2.1x108kN/m2, =0.3 Sơ đồ tính tốn BÀI 4: KẾT CẤU DẦM VÀ KHUNG Thanh chịu uốn gọi dầm Khung cấu tạo thẳng liên kết với nút hai đầu 35 Phần tử dầm (BEAM) Phần tử BEAM3 Phần tử chiều chịu kéo nén dọc trục chịu uốn thường dùng để mô dầm, khung phẳng kết cấu cơng trình Số liệu đầu vào phần tử BEAM3 Tên BEAM3 Điểm nút I, J Độ tự UX, UY, ROTZ Hằng số thực AREA, IZZ, HEIGHT,SHEARZ, ADDMAS KEYOPT(6) Lực mômen đầu ra: - Không xuất lực mômen - Xuất lực mômen hệ tọa độ phần tử KEYOPT(9) Xuất kết số điểm N nút I J (N=0, 1, 3, 5, 7, 9) KEYOPT(10) Dùng lệnh SFBEAM gán vị trí tải trọng với khoảng cách từ đầu I tới tải trọng phân bố - Giá trị chiều dài đoạn có đơn vị độ dài - Giá trị chiều dài đoạn có đơn vị tỷ số độ dài (từ đến 1.0) Số liệu đầu phần tử BEAM3 Tên gọi Định nghĩa SDIR Ứng suất lực dọc trục phần tử dầm SBYT Ứng suất uốn mép +Y phần tử dầm SBYB Ứng suất uốn mép -Y phần tử dầm SMAX Ứng suất lớn (ứng suất dọc trục + ứng suất uốn) SMIN Ứng suất nhỏ (ứng suất dọc trục - ứng suất uốn) EPELDIR Biến dạng đàn hồi dọc trục phần tử dầm EPELBYT Biến dạng uốn mép +Y phần tử dầm EPELBYB Biến dạng uốn mép -Y phần tử dầm EPTHDIR Biến dạng nhiệt dọc trục điểm cuối EPTHBYT Biến dạng nhiệt sinh uốn mép +Y phần tử dầm EPTHBYB Biến dạng nhiệt sinh uốn mép -Y phần tử dầm EPINAXL Biến dang dọc truc ban đầu phần tử MFOR(X,Y) Các lực thành phần theo phương X, Y hệ tọa độ phần tử MMOMZ Mômen trục Z hệ tọa độ phần tử 36 Phần tử BEAM54 Phần tử dầm chiều chịu kéo, nén dọc trục chịu uốn thường dùng để mơ kết cấu khung phẳng cơng trình đặt đàn hồi Số liệu đầu vào phần tử BEAM54 Số liệu vào BEAM54 Điểm nút I, J Độ tự UX, UY, ROTZ Hằng số thực AREA, IZ, HYT, HYB, AREA, AREAS, , EFS, ISTRN, ADDMAS Vật liệu EX, ALPX, DENS, GXY, DAMP KEYOPT(6) Lựa chọn ứng suất cứng hóa - Khi mở NLGEOM, sử dụng ma trận độ cứng hướng cắt 1- Khi mở NLGEOM SOLCONTROL, sử dụng ma trận độ cứng hướng cắt - Khi đóng SOLCONTROL, không sử dụng ma trận độ cứng hướng cắt KEYOPT(9) Lực mômen đầu ra: - Không xuất lực mômen - Xuất lực mômen hệ tọa độ phần tử KEYOPT(10) Dùng lệnh SFBEAM gán tải trọng bề mặt - Số liệu gán có đơn vị độ dài - Số liệu gán có đơn vị tỷ số độ dài Số liệu đầu phần tử BEAM54 Ký hiệu Định nghĩa SDIR Ứng suất lực dọc trục phần tử dầm SBYT Ứng suất uốn mép +Y phần tử dầm SBYB Ứng suất uốn mép -Y phần tử dầm SMAX Ứng suất lớn (ứng suất dọc trục + ứng suất uốn) SMIN Ứng suất nhỏ (ứng suất dọc trục - ứng suất uốn) SXY Lực cắt trung bình (phương Y) MFOR(X,Y) Các lực thành phần theo phương X, Y hệ tọa độ phần tử MMOMZ Mômen trục Z hệ tọa độ phần tử 37 Phần tử BEAM4 Phần tử dầm chiều chịu kéo nén dọc trục, chịu xoắn uốn thường dùng để mô kết cấu khung không gian Số liệu đầu vào phần tử BEAM4 Số liệu vào BEAM4 Điểm nút I, J, K (K nút định vị) Độ tự UX, UY, UZ, ROTX, ROTY, ROTZ Hằng số thực AREA, IXX, IYY, IZZ, HEIGHT,SHEARZ, ADDMAS KEYOPT(2) Lựa chọn ứng suất cứng hóa - Khi mở NLGEOM, sử dụng ma trận độ cứng hướng cắt 1- Khi mở NLGEOM SOLCONTROL, sử dụng ma trận độ cứng hướng cắt - Khi đóng SOLCONTROL, khơng sử dụng ma trận độ cứng hướng cắt KEYOPT(6) Lực mômen đầu ra: - Không xuất lực mômen - Xuất lực mômen HTĐ phần tử KEYOPT(6) Ma trận xoay tắt dần: - Không xét ma trận xoay tắt dần - Xét ma trận xoay tắt dần KEYOPT(9) Xuất kết với N điểm nút I J (N=0, 1, 3, 5, 7) KEYOPT(10) Dùng lệnh SFBEAM gán tải trọng bề mặt - Giá trị gán có đơn vị độ dài - Giá trị gán có đơn vị tỷ số độ dài Số liệu đầu phần tử BEAM4 Tên gọi Định nghĩa SDIR Ứng suất lực dọc trục phần tử dầm SBYT Ứng suất uốn mép +Y phần tử dầm SBYB Ứng suất uốn mép -Y phần tử dầm SBZT Ứng suất uốn mép +Z phần tử dầm SBZB Ứng suất uốn mép -Z phần tử dầm SMAX Ứng suất lớn (ứng suất dọc trục + ứng suất uốn) SMIN Ứng suất nhỏ (ứng suất dọc trục - ứng suất uốn) SXY Lực cắt trung bình (phương Y) MFOR(X,Y) Các lực thành phần theo phương X, Y hệ tọa độ phần tử MMOMZ Mômen trục Z hệ tọa độ phần tử 38 Phần tử BEAM44 Phần tử dầm chiều chịu kéo nén dọc trục, chịu xoắn uốn thường dùng để mơ kết cấu khung khơng gian cơng trình xây dựng Phần tử cho phép có đặc trưng hình học khơng đối xứng khác đầu cho phép đường trục có độ lệch với trọng tâm tiết diện, dùng phần tử dầm mặt cắt đối xứng phần tử BEAM4 Phần tử lò xo (COMBIN14) Phần tử lò xo COMBIN14 cho phép sử dụng toán chiều, chiều chiều Số liệu đầu vào phần tử BEAM4 Số liệu vào COMBIN14 KEYOPT (1) Loại hình tính tốn: – Tính tốn tuyến tính; – Tính tốn phi tuyến (u cầu CV2#0) KEYOPT(2) Lựa chọn độ tự I-D: - Lựa chọn sử dụng KEYOPT(3) 1- Phần tử lò xo-cản hướng dọc I-D (độ tự UX) - Phần tử lò xo-cản hướng dọc I-D (độ tự UY) - Phần tử lò xo-cản hướng dọc I-D (độ tự UZ) - Phần tử lò xo-cản hướng dọc I-D (độ tự ROTX) - Phần tử lò xo-cản hướng dọc I-D (độ tự ROTY) - Phần tử lò xo-cản hướng dọc I-D (độ tự ROTZ) - Độ tự chịu nén phần tử tự - Độ tự nhiệt độ phần tử tự Lựa chọn độ tự 2-D 3-D : -Lò xo-cản hướng dọc 3-D -Lò xo-cản dần xoay 3-D -Lò xo-cản dần hướng dọc 2-D (phần tử 2-D thiết phải mặt phẳng X-Y) 39 Một vài lưu ý tính tốn nội lực dầm khung: * Định nghĩa đặc trưng hình học cho phần tử Beam3: Preprocessor > Real Constants > Add\Edit\Delete > Real Constant Set Number 1, for Beam3 nhập số liệu sau: AREA=46.5x10-4 m2, Izz=7080x10-8m4, HEIGHT=0.30m > OK R,1, 46.5E-04,7080E-08,0.3 * Gán tải trọng phân bố hình tam giác: Solution > Define Loads > Pressure > On Beams > Xuất bảng Apply PRES on Beams > Chọn trực tiếp phần tử Beams nhịp thứ từ phần tử đến phần tử 14 > Chọn phần tử > Apply > Tại nút I nhập VAL1=0, nút J nhập VALJ=25 > Apply > Chọn phần tử > Apply > Nút I nhập VAL1=25, Nút J nhập VALJ=50 > Apply > ….> Chọn phần tử 14 > OK > Nút I nhập VAL1=25, nút J nhập VALJ=0 > OK * Định nghĩa số thực cho phần tử BEAM54: Preprocessor > Real Cosntants > Add\Edit\Delete > Real Constant Set Number 1, for BEAM54 nhập số liệu dầm đàn hồi hình vẽ bên Ngồi liên kết đứng lò xo chưa đủ để kết cấu khơng bị biến hình, cần đặt thêm liên kết ngang đặt vị trí dầm để khơng bị suy biến q trình giải thiếu lên kết * Định nghĩa đặc trưng hình học tiết diện dầm cột :Preprocessor > Beam > Common Sections > Xuất bảng Real Constant Set Number 1, for Beam 44 nhập số liệu sau: ID:1: IN40, W1= W2=0.155; W3=0.40, t1=t2=0.013, t3=0.008 ID:2: IN30, W1= W2=0.135, W3=0.30, t1=t2=0.0105, t2=0.0065 40 Ví dụ - Dầm liên tục hai nhịp Dầm liên tục nhịp có sơ đồ tính tốn hình vẽ đây, tiết diện chữ INo30 có A=132cm2, mơmen qn tính I=75450cm4, chiều cao dầm h=60cm Vật liêu thép CT3 có E=2.1×108kN/m2, =0.3 Xác định biểu đồ mômen uốn, lực cắt dầm Sơ đồ tính tốn dầm Ví dụ - Dầm đàn hồi Xác định biểu đồ mômen uốn lực cắt dầm.trên đàn hồi có sơ đồ tính tốn cho hình vẽ Tiết diện dầm có mặt cắt ngang hình chữ nhật bxh=1.0×0.8m, với A=0.8m2, Iz=0.04266m4 Cho biết P=300kN, M=90kNm, q=20kNm, Vật liệu bê tơng có mơđun đàn hồi Eb=1.94×107kN/m2, hệ số Poisson =0.2, hệ số ko=6x104kN/m3 Sơ đồ tính tốn dầm 41 Ví dụ - Khung tầng hai nhịp Xác định chuyển vị, nội lực ứng suất khung tầng hai nhịp có sơ đồ tính tốn cho hình vẽ Tiết diện dầm cột hình chữ I có đặc trưng hình học cho bảng Gối tựa phải (điểm 6) bị lún DY=-0.05m Vật liêu thép có mơđun đàn hồi E =2.4×108kN/m2, hệ số Poisson =0.3 BÀI 5: BÀI TOÁN PHẲNG Vật thể giới hạn mặt phẳng song song có chiều dày nhỏ so với hai cạnh gọi vật thể phẳng, lực tác dụng song song với mặt vật thể phân bố theo chiều dày gọi toán phẳng Trên mặt ứng suất theo phương pháp tuyến với mặt z=0 gọi tốn ứng suất phẳng, biến dạng z=0 gọi toán biến dạng phẳng 42 Phương thức xây dựng mơ hình hình học tốn phẳng * Kết hợp hình dạng định nghĩa trước ANSYS hình chữ nhật, tròn, tam giác hình đa giác * Xây dựng mơ hình bắt đầu tạo từ điểm, tiếp đến vẽ đường mặt Một vài lưu ý xây dựng mơ hình hình học tốn phẳng (1) Khi xây dựng mơ hình hình học theo hai phương thức cần sử dụng phép tốn Booleans (2) Khi xây dựng mơ hình cần phải tạo đường (Line) chung miền, dùng phép dán (Glue) dùng phép chồng (Overlap) phép tốn Booleans 43 (3) Phải gán thuộc tính vật liệu cho miền trước chia lưới phần tử (4) Có hai cách chia mạng lưới phần tử, chia tự theo kích thước chiều dài cạnh phần tử, hai chia theo mạng lưới đường biên tốn (5) Khi xây dựng mơ hình có đường biên cong coi tập hợp đoạn thẳng nhỏ, mã điểm nút cần đươc đánh số thứ tự tăng dần với gia số 1, để thuận tiện cho việc sử dụng phương thức lệnh (6) Khi thao tác chức ANSYS nên kết hợp dùng chuột lệnh có hiệu hơn, đặc biệt với tốn khơng gian Phần tử phẳng (PLANE) Phần tử PLANE42 Phần tử phẳng chiều phần tử ứng suất phẳng, phần tử biến dạng phẳng phần tử đối xứng trục Phần tử có điểm nút, nút có độ tự chuyển vị theo phương X Y 44 Phần tử PLANE183 Phần tử PLANE183 phần tử phẳng chiều phần tử ứng suất phẳng, phần tử biến dạng phẳng phần tử đối xứng trục Phần tử có điểm nút điểm nút, nút có độ tự chuyển vị theo phương X Y Hàm xấp xỉ chuyển vị phần tử PLANE183 có bậc cao hàm xấp xỉ chuyển vị phần tử PLANE42, toán với kích thước phần tử độ xác lời giải dùng phần tử PLANE183 cao dùng phần tử PLANE42 Ví dụ - Cống hộp Xác định ứng suất chuyển vị cống khoang đặt cứng, có kích thước chịu tải trọng hình vẽ Vật liệu bê tơng B20 có Eb=2.4×107kN/m2, b=0.2, b=24kN/m3 Cống đặt đá cứng Sơ đồ tính tốn cống 45 Ví dụ - Đập ngăn Xác định trạng thái ứng suất chuyển vị đập bê tơng trọng lực có kích thước chịu áp lực nước đến cao trình đỉnh đập hìnhvẽ Vật liệu bê tơng B150 có Eb=2.12x107kN/m2, b=0.2, b=24kN/m3, đá IIA có E2A=1.8×106kN/m2, 2A=0.17, IIB có E2B=8.3x106kN/m2, 2B=0.19, trọng lượng riêng nước n=10kN/m3 Mặt cắt ngang đập trọng lực 46 ... mềm ANSYS giảng dạy nghiên cứu -Điện tử -Điện gia dụng -Y học Khởi động ANSYS Khởi động nhanh Sử dụng mặc định thiết lập môi trường làm việc lần đầu Kiến nghị sử dụng phương thức Thoát khỏi ANSYS. .. SeismoSpect SolidWorks SAP2000 GEO Sigma/w ABAQUS HyperGraph SAMCEF ADINA PREP SOLU HyperView POST ANSYS FLAC3D OriginPro PLAXIS Phase2 THAM KHẢO CHÍNH ứng dụng rộng rãi ngành khoa học kỹ thuật:... hình Khống chế hiển thị PHƯƠNG PHÁP THAO TÁC APDL Soạn thảo câu lệnh Notepad với đuôi *.txt APDL (ANSYS Parametric Design Language) Ngơn ngữ tham số hóa thiết kế APDL sử dụng ngôn ngữ FORTRAN nên