Header Page of 166 Tính toán composite cốt hạt có tính đến truyền nhiệt Nghiêm Thị Thu Hà 20/11/2011 Footer Page of 166 Header Page of 166 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA TOÁN - CƠ - TIN HỌC Nghiêm Thị Thu Hà Tính toán composite cốt hạt có tính đến truyền nhiệt LUẬN VĂN THẠC SĨ KHOA HỌC Ngành: Cơ học vật thể rắn Người hướng dẫn: PGS TSKH Nguyễn Đình Đức Hà Nội - 2011 Footer Page of 166 Header Page of 166 Lời cảm ơn Đầu tiên, em xin gửi lời cảm ơn chân thành sâu sắc tới PGS TSKH Nguyễn Đình Đức Thầy tận tình bảo, hướng dẫn em kể từ em làm khóa luận tốt nghiệp đại học (2008) liên tục em hoàn thành luận văn (2011) Em vô cám ơn thầy cô giáo khoa Toán - Cơ - Tin học, trường ĐH Khoa học Tự nhiên - ĐH Quốc gia Hà Nội tận tình dạy bảo, tạo điều kiện giúp đỡ chúng em suốt năm học qua Nhân đây, em xin gửi lời cảm ơn chân thành, sâu sắc tới gia đình bạn bè động viên, khích lệ em năm học qua Em xin chân thành cảm ơn! Cuối cùng, em xin chúc thầy cô lời chúc sức khỏe, công tác tốt Chúc bạn mạnh khỏe, thành công sống Hà Nội, tháng 11 năm 2011 Học viên Nghiêm Thị Thu Hà Footer Page of 166 Header Page of 166 Mục lục Lời cảm ơn Lời mở đầu Chương Các hệ thức 1.1 Phương trình truyền nhiệt 1.2 Liên hệ ứng suất - chuyển vị 11 1.3 Lực dãn, lực tiếp, mômen uốn mômen xoắn 13 1.4 Phương trình xác định uốn 14 1.5 Điều kiện biên 17 Chương Uốn composite mỏng có truyền nhiệt dừng 19 2.1 Modun đàn hồi hệ số dãn nở nhiệt composite cốt hạt 19 2.2 Sự phân bố nhiệt độ 20 2.3 Uốn composite mỏng có truyền nhiệt dừng 21 2.3.1 Mặt không biến dạng 21 2.3.2 Biểu thức nghiệm xác định uốn 23 2.4 Tính toán số 27 2.5 Kết luận 29 Chương Uốn composite mỏng có truyền nhiệt không dừng 30 3.1 Uốn composite mỏng có truyền nhiệt không dừng 30 3.1.1 Sự phân bố nhiệt độ 30 3.1.2 Biểu thức nghiệm xác định uốn 36 3.2 Tính toán số 36 3.3 Kết luận 41 Footer Page of 166 Header Page of 166 Kết luận chung 43 Những kết nghiên cứu luận văn công bố 45 Tài liệu tham khảo 46 Phụ lục 48 Phụ lục 1: Sự phân bố nhiệt độ có truyền nhiệt dừng 48 Phụ lục 2: Độ võng có truyền nhiệt dừng 50 Phụ lục 3: Giải phương trình siêu việt phương pháp chia đôi 54 Phụ lục 4: Sự phân bố nhiệt độ có truyền nhiệt không dừng 55 Phụ lục 5: Độ uốn t = 1200s 58 Phụ lục 6: Độ uốn điểm 61 Footer Page of 166 Header Page of 166 Lời mở đầu Vật liệu composite vật liệu chế tạo tổng hợp từ hai hay nhiều vật liệu thành phần khác nhau, nhằm tạo vật liệu có tính ưu việt hẳn vật liệu thành phần ban đầu, vật liệu làm việc riêng rẽ Vì vậy, có nhiều tính ưu việt trội nhẹ, bền, đáp ứng đòi hỏi khắt khe kĩ thuật công nghệ đại Và nhờ ưu điểm bật mà chúng ngày ứng dụng rộng rãi ngành công nghiệp đại ngành chế tạo máy, hàng không, vũ trụ, tên lửa, xây dựng, ô tô, chế tạo tàu thuyền, đời sống Ví dụ composite ứng dụng làm bảng biển, pano ngành quảng cáo, trang trí nội thất, ngoại thất công trình xây dựng, ốp mặt nhà, làm trần nhà, mái vòm, hay ốp nội thất cho ô tô, tàu thuyền, Trong năm gần đây, ứng xử tác dụng tải nhiệt nhiều tác giả nghiên cứu Shariyat M [14] nghiên cứu giải tích uốn nhiệt nhiều lớp composite hình chữ nhật có tính chất vật liệu biến đổi với nhiệt độ tăng nhiệt độ sử dụng lý thuyết lớp lớn, xác định nhiệt độ uốn, từ nghiên cứu ảnh hưởng tham số tính chất hình học học composite vào nhiệt độ uốn Shiau, Kuo Chen [15] sử dụng phương pháp phần tử hữu hạn nghiên cứu chi tiết ứng xử uốn nhiệt composite nhiều lớp Wu Lanhe [10] dựa lý thuyết biến dạng trượt cấp suy phương trình cân ổn định dày vừa phải hình chữ nhật tựa lề làm từ FGM ảnh hưởng hai loại tải nhiệt tăng nhiệt gradient nhiệt thông qua bề dày tấm, suy nhiệt độ uốn, thảo luận ảnh hưởng tỉ số hướng, dày tương đối số gradient trượt ngang vào nhiệt độ uốn Trong [11, 13], tác giả trình bày giải tích uốn nhiệt chức hình chữ nhật [11], tác giả nghiên cứu tác dụng nhiệt riêng mặt phẳng tăng nhiệt thông qua bề dày tấm, đánh giá ảnh hưởng tính không đồng vật liệu, tỉ Footer Page of 166 Header Page of 166 số hướng khoảng nhiệt vào nhiệt độ uốn tới hạn, [13], với lý thuyết cổ điển suy phương trình cân bằng, ổn định, tương thích FGM không hoàn hảo tác dụng ba loại tải nhiệt tăng nhiệt đều, tăng nhiệt phi tuyến thông qua bề dày tấm, tăng nhiệt dọc trục, thu nghiệm hoàn toàn cho biến đổi nhiệt độ uốn tới hạn Trong luận văn, tác giả nghiên cứu độ võng composite hình chữ nhật có độn hạt hình cầu tựa lề cạnh chịu ảnh hưởng trình truyền nhiệt dừng không dừng Tác giả thu biểu thức nghiệm giải tích uốn có truyền nhiệt dừng không dừng Trên sở nghiệm giải tích tìm được, tác giả tính toán số để nghiên cứu ứng xử uốn làm từ vật liệu composite PVC cốt hạt Titan, qua làm rõ vai trò hạt Hiện nay, Vật liệu composite polyme độn hạt Titan ứng dụng rộng rãi Việt Nam giới Ở Việt Nam, composite polyme hạt Titan ứng dụng rộng rãi công nghiệp đóng tàu, ống dẫn dầu khí, hóa chất gần chíp sinh học sử dụng vật liệu phát quang OLED Các hạt Titan có vai trò cải thiện đáng kể tính lý vật liệu Lưu ý toán truyền nhiệt không dừng cho ống kỹ thuật composite độn hạt Titan nghiên cứu [3] Luận văn gồm: Chương 1: Các hệ thức Chương 2: Uốn composite mỏng có truyền nhiệt dừng Chương 3: Uốn composite mỏng có truyền nhiệt không dừng Kết luận chung Mặc dù cố gắng trình bày vấn đề cách mạch lạc cô đọng chắn luận văn tránh khỏi thiếu sót Vì vậy, tác giả mong nhận nhận xét, đánh giá góp ý quý thầy cô bạn để luận văn hoàn thiện Footer Page of 166 Header Page of 166 Chương Các hệ thức 1.1 Phương trình truyền nhiệt Tính truyền nhiệt môi trường đàn hồi đẳng hướng tuân theo định luật truyền nhiệt Fourier [1]: c j = −kT, j (1.1) đó, c j ( j = 1, 2, 3) thành phần vectơ dòng nhiệt, k hệ số truyền nhiệt môi trường, phải dương để bảo toàn tốc độ sản entropi dương Quá trình nhiệt đàn hồi trình thuận nghịch, nên phương trình lượng có dạng: du = σi j d εi j + dq, ρ (1.2) định luật thứ hai nhiệt động lực học có dạng dq = T ds, (1.3) tốc độ dòng nhiệt đơn vị khối lượng môi trường dq = − c j, j dt ρ (1.4) Kết hợp (1.2) (1.3), ta thu được: du = σi j d εi j + T ds, ρ (1.5) Đưa vào hàm lượng tự Helmholz f (εi j , T ) xác định hệ thức f = u −sT với biến đổi trạng thái vô nhỏ môi trường, d f vi phân toàn phần: d f = du − sdT − T ds, Footer Page of 166 (1.6) Header Page of 166 Thay (1.5) vào (1.6) ta có: df = σi j d εi j − sdT, ρ (1.7) ∂f ∂f dT + ∂ εi j ∂ T (1.8) mặt khác: df = So sánh hai hệ thức d f , suy σi j = ρ ∂f ∂f , s=− ∂ εi j ∂T Gọi F = ρ f , S = ρ s hàm lượng tự entropi đơn vị thể tích, hệ thức tên viết dạng σi j = ∂F ∂F , S=− ∂ εi j ∂T (1.9) Định luật nhiệt đàn hồi tuyến tính có dạng [1]: εij = 1+ν ν σij − σkk δij + α ∆T δij , E E (1.10) đó, ∆T = T − T0 , (1.11) với T0 nhiệt độ tuyệt đối trạng thái tự nhiên từ đây, ta biểu thị ngược lại ứng suất qua biến dạng: σij = λ εkk δij + 2µεij − (3λ + 2µ ) α ∆T δij (1.12) Từ hệ thức đầu (1.9) (1.12) ta tính biểu thức hàm lượng tự F= λ (εkk )2 + µεij εij − (3λ + 2µ ) α ∆T εkk + F0 , F0 hàm T Thay F vào hệ thức thứ hai (1.9) ta tính entropi S = (3λ + 2µ ) αεkk − dF0 dT Đặt i = j = k (1.10) ta εkk = Footer Page of 166 σkk + 3α ∆T , (3λ + 2µ ) (1.13) Header Page 10 of 166 đem thay vào (1.13), kết nhận biểu thức khác entropi S = ασkk + (3λ + 2µ ) α ∆T − dF0 dT (1.14) Nhờ biểu thức (1.12) (1.13) entropi tính tỉ nhiệt Cv biến dạng không đổi tỉ nhiệt Cp ứng suất không đổi Kết hợp (1.3) (1.4) ta ds − c j, j = T , ρ dt mặt khác ∂ s d εi j ∂ s dT ds = + , dt ∂ εi j dt ∂ T dt suy −c j, j = T ∂ S d εi j ∂ S dT + ∂ εi j dt ∂ T dt (1.15) từ suy biến dạng không đổi d εi j = T ∂∂St xác định tỉ nhiệt Cv Vậy, ∂ S (εi j , T ) d F0 = −T Cv = T ∂T dT tương tự Cp = T ∂ S (σi j , T ) d F0 = (3λ + 2µ ) α T − T ∂T dT Xem Cv ,Cp , α λ , µ số vật liệu không phụ thuộc nhiệt dộ, từ hệ thức Cv tìm biểu thức F0 T T F0 = T0 T0 Cv dT dT T Đặt kết vào (1.13) ta nhận biểu thức entropi S = (3λ + 2µ ) αεkk +Cv ln T T0 (1.16) Dùng biểu thức c j theo (1.1) S, ta đưa phương trình (1.15) dạng kT, j j = Cv Footer Page 10 of 166 ∂ εkk ∂T + (3λ + 2µ ) α T ∂t ∂t Header Page 51 of 166 xi = 0.15; k = (1 − xi) ∗ k1 + xi ∗ k2; ro = (1 − xi) ∗ ro1 + xi ∗ ro2; beta1 = 60; beta2 = 40; gama1 = beta1/k; gama2 = beta2/k; syms x y z; T = 330; T = 300; T = 293; T = −T + (T ∗ gama1 ∗ (1 + gama2 ∗ 0.01)+ T ∗ gama2 ∗ (1 + gama1 ∗ 0.01)+ gama1 ∗ gama2 ∗ (T − T 1) ∗ z)/ (gama1 + gama2 + gama1 ∗ gama2 ∗ 0.02); d = [−0.01 : 0.001 : 0.01]; f or i = : length(d) T 015(i) = subs(T, z, d(i)); end (tương tự với ξ = 0.2, ξ = 0.3) plot(d, T 01, k− ) hold on plot(d, T 015, k − − ) plot(d, T 02, k − ) plot(d, T 03, k ) Footer Page 51 of 166 49 Header Page 52 of 166 Phụ lục 2: Độ uốn có truyền nhiệt dừng clear all Em = 3e9; poism = 0.2; anpham = 8e − 5; km = 0.16; Ec = 100e9; poisc = 0.34; anphac = 4.8e − 6; kc = 22.1; Km = Em/(3 ∗ (1 − ∗ poism)); Gm = Em/(2 ∗ (1 + poism)); Kc = Ec/(3 ∗ (1 − ∗ poisc)); Gc = Ec/(2 ∗ (1 + poisc)); xi = 0.1; k = (1 − xi) ∗ km + xi ∗ kc; beta1 = 60; beta2 = 40; gama1 = beta1/k; gama2 = beta2/k; syms x y z; T = 330; T = 300; T = 293; T = −T + (T ∗ gama1 ∗ (1 + gama2 ∗ 0.01)+ T ∗ gama2 ∗ (1 + gama1 ∗ 0.01)+ gama1 ∗ gama2 ∗ (T − T 1) ∗ z)/ (gama1 + gama2 + gama1 ∗ gama2 ∗ 0.02); T momen = T ∗ z; NT = int(T, −0.01, 0.01); MT = int(T momen, −0.01, 0.01); K = Km + (Kc − Km) ∗ xi/(1 + (Kc − Km)/(Km + ∗ Gm/3)); G = Gm − 15 ∗ (1 − poism) ∗ (Gm − Gc) ∗ xi/ (7 − ∗ poism + (8 − 10 ∗ poism) ∗ Gc/Gm); Footer Page 52 of 166 50 Header Page 53 of 166 anpha = anpham + (anphac − anpham) ∗ Kc ∗ (3 ∗ Km + ∗ Gm) ∗ xi/ (Km ∗ (3 ∗ Kc + ∗ Gm) + ∗ (Kc − Km) ∗ Gm ∗ xi); S1 = 0; f or m = : : S = 0; f or n = : : w(m, n) = 864 ∗ K ∗ anpha ∗ MT /(pi2 ∗ m ∗ n ∗ 0.023 ∗ (3 ∗ K + G))/ (m2 ∗ pi2 /2.252 + n2 ∗ pi2 /1.52 − 54 ∗ K ∗ anpha ∗ NT / (0.023 ∗ (3 ∗ K + G))); w1(m, n) = w(m, n) ∗ sin(pi ∗ m ∗ x/2.25) ∗ sin(pi ∗ n ∗ y/1.5); S = S + w1(m, n); end S1 = S1 + S; end a = [0 : 0.045 : 2.25]; b = [0 : 0.03 : 1.5]; f or i = : length(a) wx01(i) = subs(S1, x, a(i)); wy01(i) = subs(wx01(i), y, b(i)); end wy01 Footer Page 53 of 166 51 Header Page 54 of 166 xi = 0.15; k = (1 − xi) ∗ km + xi ∗ kc; beta1 = 60; beta2 = 40; gama1 = beta1/k; gama2 = beta2/k; syms x y z; T = 330; T = 300; T = 293; T = −T + (T ∗ gama1 ∗ (1 + gama2 ∗ 0.01)+ T ∗ gama2 ∗ (1 + gama1 ∗ 0.01)+ +gama1 ∗ gama2 ∗ (T − T 1) ∗ z) /(gama1 + gama2 + gama1 ∗ gama2 ∗ 0.02); T momen = T ∗ z; NT = int(T, −0.01, 0.01); MT = int(T momen, −0.01, 0.01); K = Km + (Kc − Km) ∗ xi/(1 + (Kc − Km)/(Km + ∗ Gm/3)); G = Gm − 15 ∗ (1 − poism) ∗ (Gm − Gc) ∗ xi/ (7 − ∗ poism + (8 − 10 ∗ poism) ∗ Gc/Gm); anpha = anpham + (anphac − anpham) ∗ Kc ∗ (3 ∗ Km + ∗ Gm) ∗ xi/ (Km ∗ (3 ∗ Kc + ∗ Gm) + ∗ (Kc − Km) ∗ Gm ∗ xi); S1 = 0; f or m = : : S = 0; Footer Page 54 of 166 52 Header Page 55 of 166 f or n = : : w(m, n) = 864 ∗ K ∗ anpha ∗ MT /(pi2 ∗ m ∗ n ∗ 0.023 ∗ (3 ∗ K + G))/ (m2 ∗ pi2 /2.252 + n2 ∗ pi2 /1.52 − 54 ∗ K ∗ anpha ∗ NT / (0.023 ∗ (3 ∗ K + G))); w1(m, n) = w(m, n) ∗ sin(pi ∗ m ∗ x/2.25) ∗ sin(pi ∗ n ∗ y/1.5); S = S + w1(m, n); end S1 = S1 + S; end a = [0 : 0.045 : 2.25]; b = [0 : 0.03 : 1.5]; f or i = : length(a) wx015(i) = subs(S1, x, a(i)); wy015(i) = subs(wx015(i), y, b(i)); end wy015 (tương tự với ξ = 0.2, ξ = 0.3) plot(a, wy01, k− ) hold on plot(a, wy015, k − − ) plot(a, wy02, k − ) plot(a, wy03, k ) Footer Page 55 of 166 53 Header Page 56 of 166 Phụ lục 3: Giải phương trình siêu việt phương pháp chia đôi f unction v = nghiemppchiadoi( f , x0, eps) j = 1; f or i = : (length(x0) − 1) a = x0(i); b = x0(i + 1); f a = subs( f , a); f b = subs( f , b); if fa∗ fb < while (b − a) > eps x = (b + a)/2; f x = subs( f , x); f a = subs( f , a); if fx∗ fa > a = x; else b = x; end end x1 = (a + b)/2; f x1 = subs( f , x1); i f abs( f x1) < eps v( j) = x1; j = j + 1; end end end Footer Page 56 of 166 54 Header Page 57 of 166 Phụ lục 4: Sự phân bố nhiệt độ có truyền nhiệt không dừng clear all Em = 3e9; poism = 0.2; anpham = 8e − 5; km = 0.16;Cm = 900; rom = 1380; Ec = 100e9; poisc = 0.34; anphac = 4.8e − 6; kc = 22.1;Cc = 523; roc = 4500; Km = Em/(3 ∗ (1 − ∗ poism)); Gm = Em/(2 ∗ (1 + poism)); Kc = Ec/(3 ∗ (1 − ∗ poisc)); Gc = Ec/(2 ∗ (1 + poisc)); xi = 0.1; k = (1 − xi) ∗ km + xi ∗ kc; ro = (1 − xi) ∗ rom + xi ∗ roc; C = (Cm ∗ rom ∗ (1 − xi) +Cc ∗ roc ∗ xi)/(rom ∗ (1 − xi) + roc ∗ xi); beta1 = 60; beta2 = 40; gama1 = beta1 ∗ 0.02/k; gama2 = beta2 ∗ 0.02/k; syms t u z x0 = [0 : 0.01 : 25]; f = tan(u) − (gama1 + gama2) ∗ u/(u2 − gama1 ∗ gama2); muy = nghiemppchiadoi( f , x0, 0.00001) muy = 3.3907 6.4154 9.5140 12.6336 15.7619 18.8945 22.0297 T = 330; T = 300; T = 293; A0 = gama1 ∗ (T − T 0) ∗ (1 + gama2 ∗ (1/2 − z/0.02))+ gama2 ∗ (T − T 0) ∗ (1 + gama1 ∗ (1/2 + z/0.02)); B0 = gama1 + gama2 + gama1 ∗ gama2; asao = k/(C ∗ ro ∗ 0.022); T muy = 0; Footer Page 57 of 166 55 Header Page 58 of 166 f or i = : length(muy) Amuy = −2 ∗ gama1 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 − z/0.02))+ gama2 ∗ sin(muy(i) ∗ (1/2 − z/0.02)))− ∗ gama2 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 + z/0.02))+ gama1 ∗ sin(muy(i) ∗ (1/2 + z/0.02))); Bmuy = ((1 + gama1 + gama2) ∗ muy(i)2 + gama1 ∗ gama2) ∗ sin(muy(i))+ (muy(i)2 − gama1 ∗ gama2) ∗ muy(i) ∗ cos(muy(i)); deltaT muy = Amuy ∗ exp(−muy(i)2 ∗ asao ∗ t)/Bmuy; T muy = T muy + deltaT muy; end deltaT = A0/B0 + T muy; deltaT = subs(deltaT,t, 0); d = [−0.01 : 0.001 : 0.01]; f or i = : length(d) deltaT 0z(i) = subs(deltaT 0, z, d(i)); end plot(deltaT 0z, d, b− , linewidth , 1.2) deltaT 10 = subs(deltaT,t, 10); f or i = : length(d) deltaT 10z(i) = subs(deltaT 10, z, d(i)); end hold on plot(deltaT 10z, d, r − − , linewidth , 1.2) Footer Page 58 of 166 56 Header Page 59 of 166 deltaT 60 = subs(deltaT,t, 60); f or i = : length(d) deltaT 60z(i) = subs(deltaT 60, z, d(i)); end plot(deltaT 60z , d, g − , linewidth , 1.2) deltaT 120 = subs(deltaT,t, 120); f or i = : length(d) deltaT 120z(i) = subs(deltaT 120, z, d(i)); end plot(deltaT 120z, d, k , linewidth , 1.2) deltaT 600 = subs(deltaT,t, 600); f or i = : length(d) deltaT 600z(i) = subs(deltaT 600, z, d(i)); end plot(deltaT 600z, d, k∗ , linewidth , 1.2) deltaT 1200 = subs(deltaT,t, 1200); plot(deltaT 1200z, d, k− , linewidth , 1.2); (tương tự với ξ = 0.15, ξ = 0.2, ξ = 0.3) %ξ = 0.15 muy = 3.3159 6.3740 9.4859 12.6123 15.7448 18.8802 22.0175 %ξ = 0.2 muy = 0.6499 3.2756 6.3524 9.4712 12.6013 15.7359 18.8729 22.0111 %ξ = 0.3 muy = 0.5373 3.2333 6.3300 9.4561 12.5899 15.7268 18.8653 22.0046 Footer Page 59 of 166 57 Header Page 60 of 166 Phụ lục 5: Độ uốn t = 1200s clear all Em = 3e9; poism = 0.2; anpham = 8e − 5; km = 0.16;Cm = 900; rom = 1380; Ec = 100e9; poisc = 0.34; anphac = 4.8e − 6; kc = 22.1;Cc = 523; roc = 4500; Km = Em/(3 ∗ (1 − ∗ poism)); Gm = Em/(2 ∗ (1 + poism)); Kc = Ec/(3 ∗ (1 − ∗ poisc)); Gc = Ec/(2 ∗ (1 + poisc)); xi = 0.1; k = (1 − xi) ∗ km + xi ∗ kc; ro = (1 − xi) ∗ rom + xi ∗ roc; C = (Cm ∗ rom ∗ (1 − xi) +Cc ∗ roc ∗ xi)/(rom ∗ (1 − xi) + roc ∗ xi); beta1 = 60; beta2 = 40; gama1 = beta1 ∗ 0.02/k; gama2 = beta2 ∗ 0.02/k; syms t u z x y x0 = [0 : 0.01 : 25]; f = tan(u) − (gama1 + gama2) ∗ u/(u2 − gama1 ∗ gama2); muy = nghiemppchiadoi( f , x0, 0.00001); T = 330; T = 300; T = 293; A0 = gama1 ∗ (T − T 0) ∗ (1 + gama2 ∗ (1/2 − z/0.02))+ gama2 ∗ (T − T 0) ∗ (1 + gama1 ∗ (1/2 + z/0.02)); B0 = gama1 + gama2 + gama1 ∗ gama2; asao = k/(C ∗ ro ∗ 0.022); T muy = 0; Footer Page 60 of 166 58 Header Page 61 of 166 f or i = : length(muy) Amuy = −2 ∗ gama1 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 − z/0.02))+ gama2 ∗ sin(muy(i) ∗ (1/2 − z/0.02)))− ∗ gama2 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 + z/0.02))+ gama1 ∗ sin(muy(i) ∗ (1/2 + z/0.02))); Bmuy = ((1 + gama1 + gama2) ∗ muy(i)2 + gama1 ∗ gama2) ∗ sin(muy(i))+ (muy(i)2 − gama1 ∗ gama2) ∗ muy(i) ∗ cos(muy(i)); deltaT muy = Amuy ∗ exp(−muy(i)2 ∗ asao ∗ t)/Bmuy; T muy = T muy + deltaT muy; end deltaT = A0/B0 + T muy; deltaT t = subs(deltaT,t, 1200); deltaT tmomen = z ∗ deltaT t; N0sao = int(deltaT t, −0.01, 0.01); M0sao = int(deltaT tmomen, −0.01, 0.01); K = Km + (Kc − Km) ∗ xi/(1 + (Kc − Km)/(Km + ∗ Gm/3)); G = Gm − 15 ∗ (1 − poism) ∗ (Gm − Gc) ∗ xi/ (7 − ∗ poism + (8 − 10 ∗ poism) ∗ Gc/Gm); anpha = anpham + (anphac − anpham) ∗ Kc ∗ (3 ∗ Km + ∗ Gm) ∗ xi/ (Km ∗ (3 ∗ Kc + ∗ Gm) + ∗ (Kc − Km) ∗ Gm ∗ xi); S1 = 0; f or m = : : S = 0; Footer Page 61 of 166 59 Header Page 62 of 166 f or n = : : w(m, n) = 864 ∗ K ∗ anpha ∗ M0sao/(pi2 ∗ m ∗ n ∗ 0.023 ∗ (3 ∗ K + G))/ (m2 ∗ pi2 /2.252 + n2 ∗ pi2 /1.52 − 54 ∗ K ∗ anpha ∗ N0sao/ (0.023 ∗ (3 ∗ K + G))); w1(m, n) = w(m, n) ∗ sin(pi ∗ m ∗ x/2.25) ∗ sin(pi ∗ n ∗ y/1.5); S = S + w1(m, n); end S1 = S1 + S; end a = [0 : 0.045 : 2.25]; b = [0 : 0.03 : 1.5]; f or i = : length(a) wx01(i) = subs(S1, x, a(i)); wy01(i) = subs(wx01(i), y, b(i)); end wy01 (tương tự với ξ = 0.15ξ = 0.2, ξ = 0.3) plot(a, wy01, k− ) hold on plot(a, wy015, k ) plot(a, wy02, k − − ) plot(a, wy03, k − ) Footer Page 62 of 166 60 Header Page 63 of 166 Phụ lục 6: Độ uốn điểm clear all Em = 3e9; poism = 0.2; anpham = 8e − 5; km = 0.16;Cm = 900; rom = 1380; Ec = 100e9; poisc = 0.34; anphac = 4.8e − 6; kc = 22.1;Cc = 523; roc = 4500; Km = Em/(3 ∗ (1 − ∗ poism)); Gm = Em/(2 ∗ (1 + poism)); Kc = Ec/(3 ∗ (1 − ∗ poisc)); Gc = Ec/(2 ∗ (1 + poisc)); xi = 0.2; k = (1 − xi) ∗ km + xi ∗ kc; ro = (1 − xi) ∗ rom + xi ∗ roc; C = (Cm ∗ rom ∗ (1 − xi) +Cc ∗ roc ∗ xi)/(rom ∗ (1 − xi) + roc ∗ xi); beta1 = 60; beta2 = 40; gama1 = beta1 ∗ 0.02/k; gama2 = beta2 ∗ 0.02/k; syms t u z x y x0 = [0 : 0.01 : 25]; f = tan(u) − (gama1 + gama2) ∗ u/(u2 − gama1 ∗ gama2); muy = nghiemppchiadoi( f , x0, 0.00001); T = 330; T = 300; T = 293; A0 = gama1 ∗ (T − T 0) ∗ (1 + gama2 ∗ (1/2 − z/0.02))+ gama2 ∗ (T − T 0) ∗ (1 + gama1 ∗ (1/2 + z/0.02)); B0 = gama1 + gama2 + gama1 ∗ gama2; asao = k/(C ∗ ro ∗ 0.022); T muy = 0; Footer Page 63 of 166 61 Header Page 64 of 166 f or i = : length(muy) Amuy = −2 ∗ gama1 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 − z/0.02))+ gama2 ∗ sin(muy(i) ∗ (1/2 − z/0.02)))− ∗ gama2 ∗ (T − T 0) ∗ (muy(i) ∗ cos(muy(i) ∗ (1/2 + z/0.02))+ gama1 ∗ sin(muy(i) ∗ (1/2 + z/0.02))); Bmuy = ((1 + gama1 + gama2) ∗ muy(i)2 + gama1 ∗ gama2) ∗ sin(muy(i))+ (muy(i)2 − gama1 ∗ gama2) ∗ muy(i) ∗ cos(muy(i)); deltaT muy = Amuy ∗ exp(−muy(i)2 ∗ asao ∗ t)/Bmuy; T muy = T muy + deltaT muy; end deltaT = A0/B0 + T muy; deltaT momen = z ∗ deltaT ; N0sao = int(deltaT, −0.01, 0.01); M0sao = int(deltaT momen, −0.01, 0.01); K = Km + (Kc − Km) ∗ xi/(1 + (Kc − Km)/(Km + ∗ Gm/3)); G = Gm − 15 ∗ (1 − poism) ∗ (Gm − Gc) ∗ xi/ (7 − ∗ poism + (8 − 10 ∗ poism) ∗ Gc/Gm); anpha = anpham + (anphac − anpham) ∗ Kc ∗ (3 ∗ Km + ∗ Gm) ∗ xi/ (Km ∗ (3 ∗ Kc + ∗ Gm) + ∗ (Kc − Km) ∗ Gm ∗ xi); S1 = 0; f or m = : : S = 0; Footer Page 64 of 166 62 Header Page 65 of 166 f or n = : : w(m, n) = 864 ∗ K ∗ anpha ∗ M0sao/(pi2 ∗ m ∗ n ∗ 0.023 ∗ (3 ∗ K + G))/ (m2 ∗ pi2 /2.252 + n2 ∗ pi2 /1.52 − 54 ∗ K ∗ anpha ∗ N0sao/ (0.023 ∗ (3 ∗ K + G))); w1(m, n) = w(m, n) ∗ sin(pi ∗ m ∗ x/2.25) ∗ sin(pi ∗ n ∗ y/1.5); S = S + w1(m, n); end S1 = S1 + S; end wx02 = subs(S1, x, 1.125); wy02 = subs(wx01, y, 0.75); d = [0 : 10 : 1800]; f or i = : length(d) wt02(i) = subs(wy,t, d(i)); end (tương tự với ξ = 0.3) plot(d, wt02, k − − ) hold on plot(d, wt03, k − ) Footer Page 65 of 166 63 ... Uốn composite mỏng có truyền nhiệt dừng 19 2.1 Modun đàn hồi hệ số dãn nở nhiệt composite cốt hạt 19 2.2 Sự phân bố nhiệt độ 20 2.3 Uốn composite mỏng có truyền nhiệt. .. 1: Sự phân bố nhiệt độ có truyền nhiệt dừng 48 Phụ lục 2: Độ võng có truyền nhiệt dừng 50 Phụ lục 3: Giải phương trình siêu việt phương pháp chia đôi 54 Phụ lục 4: Sự phân bố nhiệt độ có truyền. .. Header Page 21 of 166 Chương Uốn composite mỏng có truyền nhiệt dừng 2.1 Modun đàn hồi hệ số dãn nở nhiệt composite cốt hạt Giả thiết composite có độn hạt hình cầu có độ dài hai cạnh a, b, chiều