Nghiên cứu lập trình phần tử hữu hạn trên ngôn ngữ Matlab

Chia sẻ: Phạm Đức Linh | Ngày: | Loại File: PDF | Số trang:7

Thêm vào BST

Báo xấu

256
lượt xem 71
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

The Finite element method is one of strong numerical methods that applied in technology problem widely and it is the essence of famous finite element program for example: SAP, ANSYS,… MATLAB is a very strong matrix analysis language and it is interested to apply for solving professional problem by researchers, engineers. This thesis research to apply CALFEM tool to set up finite element program on MATLAB determining the force adjustment of stayed-cable bridge.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Nghiên cứu lập trình phần tử hữu hạn trên ngôn ngữ Matlab

Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 THE RESEARCHING ABOUT PROGAMMING FINTE ELEMENT ON MATLAB LANGUAGE, APPLY TO FORCE ADJUSTMENT OF CABLE STAYED BRIDGE Đ B GVHD: Đ B đư . ABSTRACT The Finite element method is one of strong numerical methods that applied in technology problem widely and it is the essence of famous finite element program for example: SAP, ANSYS,… MATLAB is a very strong matrix analysis language and it is interested to apply for solving professional problem by researchers, engineers. This thesis research to apply CALFEM tool to set up finite element program on MATLAB determining the force adjustment of stayed-cable bridge. 1. được ứng dụng ,. . Nhiều phần mềm ứng dụng t đã ra đời dựa trên cơ sở phương pháp phần tử hữu hạn như: SAP,MIDAS, ANSYS.... Matlab là một công cụ phần mềm của MathWork ) có những lợi thế trong kỹ thuật lập trình đáp ứng được những hết sức đa dạng: từ các lĩnh vực kỹ thuật chuyên ngành như điện, điện tử, hống kê, … , MATLAB . . (Cable Stayed Bridge Force Adjustment). 2. Giới thiệu ”Computer Aided Learning of the FiniteElement Method”, 317
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 : hàm về ma trận, hàm về vật liệu, hàm về phần tử, môi tr hàm hệ thống (hàm phân tích tĩnh, phân tích động) 3. Cơ sở tính toán của phương pháp PTHH Phương trình cân bằng kết cấu chịu tải trọng ngoài theo phương pháp PTHH [4]: M.U’’(t) + C. U’(t) + K.U (t) = F(t) (1) Trong đó: , ma trận cản của kết cấu. U’’(t), U’ (t), U(t), F(t): Véc tơ gia tốc, vận tốc, chuyển vị nút và véc tơ tải trọng thay đổi theo thời gian. Các ma trận độ cứng, khối lượng, ma trận cản đều là ma trận vuông đối xứng, chúng được lắp ghép từ các ma trận tương ứng của từng phần tử trong kết cấu. – Trường hợp phân tích tĩnh (Static Analysis): F(t)= F Phương trình (1) trở thành: K. U = F (2) Giải hệ phương trình (2) tìm tất cả các thành phần chuyển vị tại các nút, sau đó tính nội lực ứng suất cho từng phần tử. – Trường hợp phân tích tần số dao động riêng (Eigen value Annalysis): Khi tải trọng ngoài bằng zero, bỏ qua lực cản của môi trường lúc đó kết cấu dao dộng điều hòa chuyển vị của hệ có dạng: U=U. sin( t) và U’’ = -U. 2 . sin( t) (3) 2 -M.U. . sin( t) + K. U. sin( t) = {0} 2. (K - M). U = {0} (4) Giải phương trình (4) bằng phương pháp SUBSPACE sẽ cho các giá trị riêng và véc tơ riêng từ đó tính được các tần số riêng (eigen frequencies) và dạng dao động riêng (mode shape) tương ứng. Trong . 4. 4.1. : [ Y ].{X}+{Y0}+{YII}+{YC}=0 [5] Trong đó: [ Y ]: ma trận ảnh hưởng độ võng, phần tử yij là độ võng tại nút thứ i do lực căng dây tại nút thứ j có giá trị bằng đơn vị gây ra (tương ứng với sơ đồ của hệ ở thời điểm căng dây tại nút thứ j). 318
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 ) {Y0}:vectơ độ võng của hệ xuất phát (A). {YC ). {YII}: Vectơ độ võng do tĩnh tải phần 2 gây ra trong hệ ở trạng thái hoàn chỉnh 4.2. : 145 + 293 + 145m : DC = 632.73 kN/m n II: DW = 41.36 kN/m = 1000 kN, M = 500 kN.m 4.3. Mô hình phần tử hữu hạn và các yếu tố được xem xét 31 24 975 6 8 10 12 14 16 18 20 22 39 37 35 33 31 29 27 25 23 21 19 17 15 13 11 24 26 28 30 32 34 36 38 40 1: PTHH 3 bậc tự d ). u5 u4 a. b. u6 u4 u3 (x2,y2) (x2,y2) y y E, A, I E,A x x u2 u2 u3 u1 u1 x x (x1,y1) (x1,y1) nh 2: ; S 3: 319
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 % ----- Geometry ------------------------------------ BÀÕT ÂÁÖU - %%BEGIN GEOMETRY Geometry=gnew('geometry'); %%END GEOMETRY NHÁÛP SÄÚ LIÃÛU ÂÁÖU VAÌO % ----- Materials ----------------------------------- - ÂÀÛC TRÆNG HÇNH HOÜC, VÁÛT LIÃÛU, TAÍI TROÜNG % Mtrl=mat(Mtrl,pos,ep,D); %%BEGIN MATERIALIDX Mtrlidx=gnew('materialidx'); MÄ HÇNH HOÏA CÁÖU QUA TÆÌNG GIAI ÂOAÛN: %%END MATERIALIDX LÁÛP CAÏC MA TRÁÛN PHÁÖN TÆÍ, NUÏT gmdraw2(Geometry); % ----- System matrices ----------------------------- - ndomains=size(Geometry.surfaces,2); Gdof=cell(ndomains,1); LÁÛP CAÏC MA TRÁÛN ÂÄÜ CÆÏNG CHO TÆÌNG Dd=cell(ndomains,1); PHÁÖN TÆÌ VAÌ LÀÕP GHEÏP VAÌO MA TRÁÛN ÂÄÜ Es=cell(ndomains,1); CÆÏNG CAÍ HÃÛ TRONG TÆÌNG GIAI ÂOAÛN Et=cell(ndomains,1); % ----- Show element mesh --------------------------- - figure; LÁÛP CAÏC MA TRÁÛN LÆÛC VAÌ MA TRÁÛN ÂIÃÖU ddraw2(Domains); KIÃÛN BIÃN TRONG TÆÌNG GIAI ÂOAÛN zoom on; % ----- Create Gdof -------------------------------- -- [Gdof,Domains,dofcount]=gtopo(Geometry,Domains); GIAÍI CAÏC PHÆÅNG TRÇNH MA TRÁÛN ÂÃØ XAÏC % ----- Assemble stiffness matrix ------------------- ÂËNH CHUYÃØN VË TAÛI NUÏT - K=zeros(dofcount); (TAÛI LIÃN KÃÚT GIÆÎA DÁY VÀNG VÅÏI DÁÖM) f=zeros(dofcount,1); for i=1:ndomains [K,f]=assem(Gdof{i},K,Domains(i).K,f,Domains(i).f); XUÁÚT RA CAÏC CHUYÃØN VË TRONG TÆÌNG end; GIAI ÂOAÛN RÄÖI GHEÏP VAÌO MA TRÁÛN CHUYÃØN VË CAÍ CÁÖU. % ----- Load vector % f=ldcurve(f,Geometry,Domains,icurve,qx,qy); % f=ldcurve(f,Geometry,Domains,icurve,q); GIAÍI PHÆÅNG TRÇNH: Ydv*X + Ytt = 0 TA COÏ ÂÆÅÜC KÃÚT QUAÍ LÆÛC CÀNG % ----- Boundary conditions ------------------------- TRONG CAÏP X - bc=[]; % ----- Solve equation system ----------------------- TÊNH LAÛI LAÛI ÂÄÜ VOÎNG TOAÌN CÁÖU SAU - ÂIÃÖU CHÈNH . XUÁÚT KÃÚT QUAÍ TÊNH TOAÏN [a,Q]=solveq(K,f,bc); THEO YÃU CÁÖU % ----- Element forces ------------------------------ -for i=1:ndomains Dd{i}=gextract(Domains,i,a); KÃÚT THUÏC [Es{i},Et{i}]=planqss(Domains,i,Mtrl(Mtrlidx(i).Idx), Dd{i}); end; % ----- Post processing ----------------------------- 3: CBFA - 4.4. như: Ma Edof ep, bc,... Kei: Kei như sau: BEAM 320
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 Kei như sau: fei tr . Kei fei K F . K.U = F . . . 4.5. L : Dây văng X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 (kN) 3919 4071.5 3265.5 3721 3138 4017.5 3107.5 4447 3134.5 4861 Dây văng X11 X12 X13 X14 X15 X16 X17 X18 X19 X20 (kN) 3222.5 5172 3378 5353 3591 5426 3838 5434.5 4093.5 5420.5 Dây văng X21 X22 X23 X24 X25 X26 X27 X28 X29 X30 (kN) 4341 5410.5 4572 5418 4784 5446.5 4978 5494 5156 5556.5 Dây văng X31 X32 X33 X34 X35 X36 X37 X38 X39 X40 (kN) 5319.5 5630 5471 5711 5611.5 5796.5 5800 5862 7029.5 7450 1: ): X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 - -0.01 0.009 -0.05 0.038 -0.17 0.121 -0.40 0.291 0.488 (mm) 0.720 321
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 X11 X12 X13 X14 X15 X16 X17 X18 X19 X20 -0.94 0.503 -0.77 0.052 -0.12 -1.037 0.959 -2.735 2.267 -4.864 (mm) X21 X22 X23 X24 X25 X26 X27 X28 X29 X30 3.564 -7.21 4.653 -9.581 5.408 -11.86 5.764 -13.97 5.703 -15.86 (mm) X31 X32 X33 X34 X35 X36 X37 X38 X39 X40 5.246 -17.49 4.445 -18.85 3.381 -19.91 2.167 -20.65 6.737 -20.82 (mm) 2: . 4.6. 7.01 CBFA MIDAS v7.01 2. (a) (b) (c) 4: 2 trong MIDAS v7.01; 1kN trong MIDAS v7.01 3: CBFA MIDAS X1 X2 X3 X4 X1 X2 X3 X4 -2.9322 -2.9322 -7.0567 -7.0567 -3.2193 -3.2193 -7.0638 -7.0638 3 0.0018 -0.0012 0.0033 -0.0019 0.0018 -0.0013 0.0035 -0.0018 p X4 -0.0012 0.0018 -0.0019 0.0033 -0.0013 0.0018 -0.0018 0.0035 3: 7.01 (%) 3 4 X1 -8.9181 0 -7.6923 X2 -8.9181 -7.6923 0 X3 -0.1005 -5.7143 5.55556 X4 -0.1005 5.55556 -5.7143 4: 322
Tuyển tập Báo cáo Hội nghị Sinh viên Nghiên cứu Khoa học lần thứ 7 Đại học Đà Nẵng năm 2010 5. Kết luận và kiến nghị au: . to MIDAS (sai s 10%). ,.. . : x . ,... [1] , 2007 [2] , 2005 [3] , , Nxb. K , 2008 [4] Lund university, C A L F E M a finite element toolbox to MATLAB Version 3.3, February 1999 [5] 323