Nghiên cứu trạng thái ứng suất của vỏ trụ composite lớp dưới tác dụng của áp suất trong và nhiệt trên cơ sở lý thuyết biến dạng..

Nội dung tài liệu: Nghiên cứu trạng thái ứng suất của vỏ trụ composite lớp dưới tác dụng của áp suất trong và nhiệt trên cơ sở lý thuyết biến dạng..

1. BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG VIỆN KHOA HỌC VÀ CÔNG NGHỆ QUÂN SỰ NGUYỄN TRƯ ỜNG THANH NGHIÊN CỨU TRẠNG THÁI ỨNG SUẤT CỦA VỎ TRỤ COMPOSITE LỚP DƢỚI TÁC DỤNG CỦA ÁP SU ẤT TRONG VÀ NHIỆT TRÊN CƠ SỞ LÝ THUYẾT BIẾN DẠNG TRƢỢT BẬC CAO QUASI-3D Ngành: Cơ kỹ thuật Mã số: 9 52 01 01 TÓM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT Hà Nội, 2023
2. CÔNG TRÌNH ĐƢỢC HOÀN THÀNH TẠI VIỆN KHOA HỌC VÀ CÔNG NGH Ệ QUÂN SỰ-BỘ QUỐC PHÒNG Người hướng dẫn khoa học: 1. PGS.TS Trần Ngọc Đoàn 2. TS Phan Văn Chương Phản biện 1: GS. TS Hoàng Xuân Lượng Học viện Kỹ thuật quân sự Phản biện 2: PGS. TS Nguyễn Mạnh Cường Đại học Bách khoa Hà Nội Phản biện 3: PGS. TS Trần Ngọc Thanh Viện Khoa học và Công nghệ quân sự Luận án được bảo vệ tại hội đồng đánh giá luận án tiến sĩ cấp Viện Khoa học công nghệ quân sự họp tại Viện Khoa học và Công nghệ quân sự vào hồi giờ phút, ngày tháng năm 20 Có thể tìm hiểu luận án tại: - Thư viện Viện Khoa học và Công nghệ quân sự; - Thư viện Quốc gia Việt Nam.
3. 1 MỞ ĐẦU 1. Tính cấp thiết của đề tài Nghiên cứu vỏ trụ composite có tính đến ảnh hưởng của biến dạng trượt bậc cao cho phép đánh giá chính xác hơn trạng thái ứng suất của vỏ, đặc biệt là dưới tác dụng của tải trọng phức tạp. Từ đó, đánh giá chính xác hơn trạng thái ứng suất tại những vùng chuyển tiếp của kết cấu, vùng biên, vùng chịu tải trọng tập trung, phục vụ cho quá trình thiết kế, chế tạo các sản phẩm công nghệ cao, đáp ứng các yêu cầu khắt khe về khối lượng và độ bền như trong lĩnh vực kỹ thuật tên lửa, động cơ phản lực và lĩnh vực hàng không vũ trụ. 2. Mục tiêu nghiên cứu Mục tiêu nghiên cứu của luận án là nghiên cứu trạng thái ứng suất của vỏ trụ composite lớp chịu tác dụng cuả tải trọng cơ-nhiệt khác nhau trên cơ sở lý thuyết biến dạng trượt bậc cao kiểu Quasi-3D; thông qua việc phân tích ứng suất và đánh giá ảnh hưởng của các tham số khác nhau như điều kiện biên, tham số vật liệu, hình học, tải trọng, lên trạng thái ứng suất vỏ, rút ra các kết luận có giá trị khoa học và thực tiễn 3. Nội dung nghiên cứu - Nghiên cứu xây dựng mô hình bài toán, xây dựng các quan hệ ứng xử cơ học và thành lập hệ phương trình xác định trạng thái ứng suất-biến dạng cho vỏ composite lớp chịu tác dụng tải trọng cơ-nhiệt. - Nghiên cứu xây dựng thuật toán giải hệ phương trình bằng phương pháp giải tích và lập chương trình tính toán trạng thái ứng suất-biến dạng của vỏ trụ composite lớp dưới tác dụng tải trọng cơ- nhiệt với các điều kiện biên trên cơ sở lý thuyết biến dạng trượt bậc cao. Các thành phần ứng suất cắt được chính xác hóa bằng phương trình cân bằng của lý thuyết đàn hồi ba chiều. - Khảo sát ảnh hưởng của các thông số kết cấu, tải trọng, vật liệu, điều kiện biên, đến trạng thái ứng suất-biến dạng của vỏ. Trên cơ sở đó đề xuất các khuyến nghị có ý nghĩa khoa học phục vụ trong thiết kế, chế tạo, khai thác, sửa chữa đối với kết cấu loại này trong thực tiễn và trong kỹ thuật. 4. Đối tƣợng và phạm vi nghiên cứu - Đối tượng nghiên cứu + Về kế cấu: Vỏ trụ composite lớp với các điều kiện biên; + Về tải trọng: Vỏ trụ chịu áp suất trong phân bố đối xứng với điều kiện nhiệt độ ổn định theo thời gian.
4. 2 - Phạm vi nghiên cứu: Xác định trạng thái ứng suất, biến dạng cho vỏ composite lớp chịu tác dụng của áp suất và nhiệt độ trên cơ sở biến dạng trượt bậc cao Quasi-3D với điều kiện vỏ làm việc trong giới hạn đàn hồi tuyến tính. 5. Phƣơng pháp nghiên cứu Sử dụng phương pháp giải tích để giải bài toán biên, kết hợp tính toán, phân tích số và so sánh với các kết quả đã công bố trên các tạp chí chuyên ngành uy tín trong nước và trên thế giới. 6. Ý nghĩa khoa học và thực tiễn của luận án - Ý nghĩa khoa học: Xây dựng được phương pháp giải tích xác định trạng thái ứng suất vỏ trụ composite lớp dưới tác dụng của áp suất trong và nhiệt trên cơ sở lý thuyết biến dạng trượt bậc cao với các điều kiện biên khác nhau trong điều kiện tải đối xứng trục. - Ý nghĩa thực tiễn: Kết quả nghiên cứu của luận án đã giải thích được vùng ứng suất mạnh tại biên bằng cơ sở toán. Kết quả nghiên cứu này có giá trị cao trong ứng dựng thực tế thiết kế, chế tạo các dạng kết cấu vỏ. 7. Bố cục của luận án Toàn bộ luận án gồm: phần mở đầu, 4 chương, kết luận chung, danh mục 07 công trình nghiên cứu khoa học được công bố và 97 tài liệu tham khảo. Chƣơng 1. TỔNG QUAN VỀ PHƢƠNG PHÁP TÍNH TOÁN KẾT CẤU VỎ COMPOSITE LỚP 1.1. Tổng quan về vật liệu composite lớp Trình báy khái quát về chung về vật liệu composite lớp, ưu nhược điểm và một số ứng dụng của kết cấu vỏ sử dụng vật liệu composite. 1.2. Tổng quan về lý thuyết tính toán kết cấu vỏ composite Phân tích một số nghiên cứu trong và ngoài nước đã công bố về kết cấu vỏ composite trên cơ sở các mô hình lý thyết vỏ như (lý thuyết vỏ cổ điển, lý thuyết trượt bậc nhất, lý thuyết trượt bậc cao và lý thuyết đàn hồi 3D). 1.3. Tổng quan tình hình nghiên cứu kết cấu vỏ composite Nghiên cứu tổng quan các tình hình nghiên cứu kết cấu vỏ composite từ các công bố trong và ngoài nươc. 1.4. Những vấn đề cần tiếp tục nghiên cứu Từ phân tích, tổng hợp trên rút ra các vấn đề cần tiếp tục nghiên cứu, phát triển của luận án tập trung giải quyết như sau:
5. 3 - Nghiên cứu trạng thái ứng suất-biến dạng của vỏ composite lớp bằng phương pháp giải tích với các điều kiện biên khác nhau dựa trên lý thuyết biến dạng trượt bậc cao kiểu Quasi-3D; - Nghiên cứu trạng thái ứng suất tại vùng ứng suất mạnh đối với các kết cấu vỏ composite lớp khi kết cấu dưới tác dụng của tải cơ, nhiệt và cơ-nhiệt đồng thời; - Nghiên cứu đáp ứng tĩnh, động và ổn định kết cấu tấm, vỏ composite theo lý thuyết bậc cao phi tuyến; - Tối ưu hóa kết cấu để giảm hiện tượng tập trung ứng suất đối với các kết cấu composite. 1.5. Kết luận chƣơng 1 Nghiên cứu tổng quan về kết cấu làm bằng vật liệu composite lớp đã và đang được ứng dụng rộng rãi trong thực tế. Nhiều nghiên cứu về trạng thái ứng suất biến dạng của vỏ trụ composite lớp của các nhà nghiên cứu trong và ngoài nước sử dụng các lý thuyết khác nhau như: lý thuyết cổ điển, lý thuyết bậc nhất, lý thuyết bậc cao với các điều kiện biên khác nhau. Trên cơ sở sử dụng đa dạng các phương pháp nghiên cứu (phương pháp phần tử hữu hạn, phương pháp giải tích hoặc là cả thực nghiệm kiểm chứng). Tuy nhiên, trong các nghiên cứu vẫn chỉ tập trung nhiều vào lý thuyết bậc nhất hoặc bậc cao với điều kiện biên biên tựa đơn. Mặc dù gần đây đã có các nghiên cứu cho kết cấu vỏ trụ composite với điều kiện biên bất kỳ nhưng hạn chế chủ yếu là sử dụng phương pháp phần tử hữu hạn hoặc phương pháp số khác. Phân tích vùng ứng suất mạnh đối với vỏ trụ composite lớp với điều kiện biên khác nhau trên cơ sở lý thuyết biến dạng trượt bậc cao kiểu quasi-3D theo hướng tiếp cận giải tích còn chưa được quan tâm nghiên cứu đầy đủ, đặc biệt là trong trường hợp vỏ chịu tác dụng tải trọng cơ-nhiệt đồng thời. Luận án sẽ đi nghiên cứu xây dựng mô hình, thiết lập hệ phương trình vi phân mô tả trạng thái ứng suất vỏ trụ composite lớp từ đó chỉ ra vùng ứng suất mạnh và khảo sát các yếu tố ảnh hưởng đến hiệu ứng này. Chƣơng 2. NGHIÊN CỨU TÍNH TOÁN TRẠNG THÁI ỨNG SUẤT CỦA VỎ TRỤ COMPOSITE LỚP 2.1. Mô hình bài toán và các giả thiết Kết cấu vỏ trụ composite có n lớp, mỗi lớp là vật liệu composite cốt sợi đồng phương, có các thông số hình học, như hình 2.4.
6. 4 Z q w z u n x h i 2 v Zi q Z1 R h 2 Z0 q(x) 2 1 T0 T1 L Hình 2.4. Kết cấu vỏ trụ composite lớp dưới tác dụng của tải trọng 2.2. Xây dựng mô hình tính toán vỏ trụ composite lớp Trường chuyển vị của vỏ trụ theo [76], được tính theo công thức: u u x,,;,,;,, q z v x q z w x q z KKKzi z i 1 z i (2.1) ui x,;,;, q v i x q w i x q i!!!  i  i i 0 i 0 i 0 Mối quan hệ tuyến tính giữa các thành phần của trường biến dạng và trường chuyển vị của vỏ trụ trong hệ tọa độ trụ OZxq :     (2.2) k0 kT k Ở đây, lần lượt là biến dạng của vỏ do tác dụng của tải cơ học và tải nhiệt độ, được thể hiện như sau. t 0 x,,,,,  q z  xq  q z  x z  k k (2.3)  TTT, , ,0,0,0t T kk1 2 3 với: 1u 1  v 1  v 1  u x ;;;  q w  xq Rx R z  q R  x R z  q (2.4) 1w  v v 1  w  u 1  w  ;;.   qxzR z qx  z R z z R   z z R  z Quan hệ ứng suất-biến dạng lớp k trong hệ tọa độ lớp vật liệu: ()k  ()()()()kkkk x x QQQQ11 12 13 14 00 ()k  ()()()()kkkk  q QQQQ21 22 23 24 00 q ()k ()()()()kkkk  z QQQQ 00  z  31 32 33 34 , ()k ()()()()kkkk  (2.10) xq QQQQ41 42 43 44 00 xq ()k 0000QQ()()kk  x z 55 56 x z ()()kk ()k 0000QQ 65 66  q z q z Khi đó hệ phương trình đạo hàm riêng cân bằng của kết cấu đối
7. 5 với mỗi mô hình K khác nhau. Ở đây với K=3 được thiết lập như sau: NNNN(0)  (0)  (0)  (0) x qx 0, xq x  N( 0 ) 0, x  q  x  q q z ( 1 ) ( 1 ) ( 1 ) ( 1 ) NNNx  qx( 0 )  xq Nq ( 0 ) RNxqzz 0, RN 0, x  q  x  q NNN( 2 )  ( 2 )  ( 2 ) N ( 2 ) x qx RN( 1 ) 0, xq q RN ( 1 )  N ( 2 ) 0, x  qxz  x  q q z q z NNN( 3 )  ( 3 )  ( 3 ) N ( 3 ) x qx RN( 2 ) 0, xq q RN ( 2 ) 2N ( 3 ) 0, (2.15b) x  qxz  x  q q z q z ( 0 ) ( 0 ) Nx z Nq z ( 0 ) Nq Rp0 0, xq N( 1 ) N( 1 ) x z q z  N( 1 ) RN ( 0 ) Rp 0, xqq z z 1 N( 2 ) N ( 2 ) x z q z N( 2 ) RN ( 1 ) Rp 0, xqq z2 Các thành phần lực màng và mômen uốn được xác định như sau: h NL ( k ) z N(0)(1) N N (K)  (k) 1 1 z z K dz, x x x  x k1 R h( k 1 ) h NL ( k ) z N(0)(1) N N (K)  (k) 1 1 z z K dz, q q q  q k1 R h( k 1 ) h NL ( k ) z N(0)(1) N N (K2)  (k) 1 1 z z K2 dz, z z z  z k1 R h( k 1 ) h NL ( k ) N(0)(1) N N (K)  (k) 1 z z K dz, qx qx qx  qx (2.16) k1 h( k 1 ) h NL ( k ) z N(0)(1) N N (K)  (k) 1 1 z z K dz, xq xq xq  xq k1 R h( k 1 ) h NL ( k ) z N(0)(1) N N (K1)  (k) 1 1 z z K1 dz, xz x z x z  x z k1 R h( k 1 ) h NL ( k ) N(0)(1) N N (K)  (k) 1 z z K dz, qz q z q z  q z k1 h( k 1 )
8. 6 h NL ( k ) QQ11 xq 12 z NNNN(0)(1)(2)(3) T 1 1zzzdz, 23 TTTTx x x x  k1 QQ R h( k 1 ) 13 z 14 xq h NL ( k ) QQ NNNN(0)(1)(2)(3) T21xq 22 1zzzdz, 23 TTTTq q q q  k1 QQ h( k 1 ) 23 z 24 xq h NL ( k ) QQ31 xq 32 z NNNN(0)(1)(2)(3) T 1 1zzzdz, 23 Tz Tz Tz Tz  k1 QQ R h( k 1 ) 33 z 34 xq h NL ( k ) QQ41 xq 42 z NNNN(0)(1)(2)(3) T 1 1zzzdz, 23 TTTTxq xq xq xq  k1 QQ R h( k 1 ) 43 z 44 xq h NL ()k QQ N(0)(1)(2)(3) N N N T41xq 42 1, z z 23 z dz TTTTqx qx qx qx  k 1 QQ h(k 1) 43z 44 xq i i h h h h pq1i ii q1 ,i0,1,1  2R 2 i! 2R 2 i! Các điều kiện biên thường gặp được xác định như sau: - Đối với biên ngàm chặt: Tại xx 0, 0 : ui v i w j 0;i 1,K;j 1,K 1 (2.17) Khi vỏ dạng panel thì tại qq 0, 0 : u v w 0;i 1,K;j 1,K 1 i i j - Đối với biên gối tựa: Tại xx 0, 0 : (i) (i) (i) (i) Nx NT x N q N T q 0,w j 0,i 1,K;j 1,K1 ( i ) ( i ) (2.18) hoặc Nxx NT 0;v i w j 0;i 1,K;j 1,K1 ( i ) ( i ) Tại qq 0, 0 : Nxq N 0,wj 0,i 1,K;j 1,K1 ( i ) ( i ) hoặc Nxq N 0,wj 0,i 1,K;j 1,K1 - Đối với biên tự do: Tại xx 0, :N( i ) N ( i ) N ( j ) 0,i 1,K;j 1,K12 : 0zx xq x (2.19) ()()()i i j Tại qq 0,0 :NNNx q q z 0, i 1,K; j 1,K 1 2.3. Xây dựng phƣơng pháp giải tích tính toán trạng thái ứng suất vỏ Hệ phương trình cân bằng theo chuyển vị như sau:
9. 7 KK 2  2  2 Hl H l H l u H l v  1i 1 i ,1122 1 i ,22 i 2 i ,12 i ii 00 xq xq K 1 l   H3ii ,1 w 0, l 1, , K 1 , i 0 x KK2  2  2 Hm u H m H m H m v 1i ,12 i 2 i 2 i ,1122 2 i ,22 i (2.20) ii 00xq xq K 1 m   H3ii ,2 w 0, m K 2 , , 2 K 2 , i 0 q KKK 1 22   Hn u H n H n H n w H n v 1,1i i  3 i 3,11 i22 3,22 i i  2,2 i i i 0xq i 0 xq i 0 n n n n HqHqHTHTn , 2 K 3, ,3 K 2, qqT0 0 Ti i Các phương trình (2.20) cho phép xác định các thành phần của trường chuyển vị u ui,,, v i w i  của vỏ trụ. Từ đó tìm được các biến dạng x,,,,  q zz  xq  x và  q z . Tiếp theo, sử dụng ba biểu thức quan hệ biến dạng và ứng suất của Hooke để xác định các ứng suất mặt xq, và xq . Các thành phần ứng suất cắt còn lại được tìm dựa trên phương trình cân bằng của lý thuyết đàn hồi 3D như công thức sau: z 1 z x xq x z 1, dz R z R xq  h/2 z 2 R zz q xq  1 1 dz , (2.21) q z 2 RRqx ()Rz h/2 z 1 zx z q z R h / 2 z 1. q dz q R z R xq  R h /2 h/2 Trong nghiên cứu vỏ trụ tròn kín, điều kiện biên đối với biên qq 0, sẽ được thay thế bằng điều kiện tuần hoàn. Sử dụng khai 0 triển các hàm tải trọng và chuyển vị được cho bởi các công thức sau: (1) (2) q x, q Q0 x  Qmm x cos m q Q x sin m q m 1 (1) (2) To T o0 x  T om x cos m q T om x sin m q m 1 (1) (2) Tin T in0 x  T inm x cos m q T inm x sin m q (2.22) m 1
10. 8 (1) (2) ui x, q U i0 x  U im x cos m q U im x sin m q m 1 (1) (2) vi x, q V i0 x  V im x sin m q V im x cos m q m 1 (1) (2) wj x, q W j0 x  W jm x cos m q W jm x sin m q m 1 Rút ra được hệ phương trình vi phân thường cấp hai cho các hàm UVWi0,, i 0 j 0 biến dạng: K2 K 1 l ld l dWj0 H1i H 1i,11 U i0 H 3i,1 0,l1, ,K1, dx 2 dx i 0 i 0 KdU K 1 d 2 HHHWni0 n n 1i,1 3i 3i,112 j0 (2.23) i 0dx i 0 dx n n n n HQ HQ H T H T,n 2K 3, ,3K 2 qq0 0 To o0 Tin in0 K d 2 Hm H m mH 2 m V 0,m K2, ,2K2.  2i 2i,112 2i,22 i0 i0 dx 1 1 1 2 2 2 Còn đối với dạng hàm UVWim x ,, im x im x và UVWim x ,, im x im x nhận được hệ phương trình vi phân thường cấp hai sau: KK dd2 Hl H l m 2 H l U m H l V  1i 1i,112 1i,22 im 2i,12 im i 0 dx i 0 dx K1 l d  H3i,1 W im 0, l 1, , K 1 , i0 dx KK2 mdd m m2 m m H1i ,12 U im H 2 i H 2 i ,11 m H 2 i ,22 V im dx dx 2 ii 00 (2.24) K 1 m m H3i ,2 W im 0, m K 2 , , 2 K 2 , i 0 KKKdd 1 2 Hn U H n H n m2 H n W m H n V 1,1i im  3 i 3,11 i2 3,22 i im  2,2 i im i 0dx i 0 dx i 0 n()()()() s n s n s n s HQHQHTHTnK , 2 3 , , 3 K 2 qqm m To om Tin inm Biến đổi Laplace từ hệ phương trình vi phân thường sang hàm ảnh, nhận được hệ phương trình đại số ẩn p: KKK l2 l l l  H1i pH 1 i ,11 UppHC i().()  1 i 1 i 0 p  HWp 3 i i i 0 i 0 i 0
11. 9 KK qq13 13 ll H0 LapQ 13 H 0 LapQ 13  H 111i C i H 330 i C i , l 1,2,3,4; i 00 KK m2 m m  H2i p H 2 i V i() p pH 1 i C 2 i 0 ii 00 K (2.30) qq23 23 m H0 Q 23 H 0 Q 23  H 1ii C 2 1, m 5,6,7,8; i 0 KKK j2 j m j qq33 33  H0 p H 1i W i().() p  pH 130 i C i  p H 1 U i p H 033033 Q H Q i 0 i 0 i 0 KK jj jjHTHTTo o00 Tin in H11 C 3i 1 H 1 i C 1 i 0 , j 9,10,11 ii 00 p Ở dạng ma trận, hệ phương trình (2.58) được biến đổi thành: A .[U ,U ,U ,U ,W ,W ,W ,W ]T = 0 ,   00 10 20 30 00 10 20 30  (2.34) trong đó, A là một ma trận hệ số của các phương trình (2.51). Khi đó Uij00, W , i 0,3, j 0,2, được xác định bởi các công thức: 8 8 8 8 U det A ,U det A ,U det A ,U det A , 00 1s 10  2s 20  3s 30  4s s 1 s 1 s 1 s 1 (2.36) 8 8 8 8 W det A ,W det A ,W det A ,W det A , 00 5s 10  6 s 20  7s 30  8s s 1 s 1 s 1 s 1 Chuyển vị của vỏ trụ composite lớp được xác định dư ới dạng tổng của các hàm phân số dạng: 0i 1i 2i2 3i3 03 2 K K p K p K p K UUUUU0 0 0 0 0 U0 p LapQ 33 p  ,  2222 i1 p p q p p q p  i i i i 0i 1i 2i 2 3i 3 2 K K p K p K p UUUUj j j j Uj p LapQ 33 p  , j 1,2,3, (2.42)  2222 i1 p p q p p q  i i i i 0i 1i 2i 2 3i 3 2 K K p K p K p WWWWl l l l Wl p LapQ 33 p  ,l 0,3.  2222 i1 p p q p p q i i i i  Các hệ số trong biểu thức (2.70) được xác định bằng cách so sánh các hệ số tại p trong các biểu thức tương ứng cho hàm ảnh của mỗi
12. 10 chuyển vị. 2.4. Xây dựng sơ đồ thuật toán chƣơng trình tính toán ĐẦU VÀO Nhập các thông số vỏ: thông số kết cấu (L, R, h); thông số vật liệu (E,G, v); thứ tự xếp lớp và số lớp; dạng liên kết; tải trọng q; Xác định các hệ số của hệ phương trình cân bằng Xác định hệ phương trình đối với hàm ảnh các biến số Xác định hàm ảnh của trường chuyển vị vỏ trụ Biến đổi Laplace ngược Biểu thức chuyển vị tổng quát tương ứng với điều kiện biên tại x L/R Xác định chuyển vị: u,v,w Xác định trường biến dạng x, q, z, xq, xz, qz Xác định ứng suất màng: x, q,xq Xác định ứng suất cắt: Z, xz,qz KẾT THÚC Xuất file dữ liệu: trường chuyển vị; trường ứng suất; Hình 2.10. Sơ đồ khối tính toán kết cấu vỏ composite lớp
13. 11 2.5. Kiểm chứng mô hình và chƣơng trình tính toán Bài toán 1: (theo bái toán Varadan-Bhaskar [88]) Bảng 2.2 Chuyển vị không thứ nguyên của vỏ composite lớp [900/0/900/0/900] [900] [900/0] [900/0/900] 2 Luận Luận án Luận án Luận án [88] [88] án [88] [88] S K=3 K=3 K=3 K=3 4 2.7723 2.7830 7.3555 6.1000 3.8021 4.0090 3.9979 4.2060 10 0.9172 0.9189 3.6490 3.3300 1.0971 1.2230 1.3437 1.3800 50 0.5384 0.5385 2.2516 2.2420 0.5436 0.5495 0.7609 0.7622 100 0.5169 0.5170 1.3682 1.3670 0.4703 0.4715 0.6259 0.6261 500 0.3060 0.3060 0.1005 0.1005 0.1027 0.1027 0.1006 0.1006 *Nhận xét: Sai số kết quả tính toán chuyển vị hướng kính không thứ nguyên vỏ trụ composite lớp của luận án với kết quả công bố của tác giả T.K.Varadan-K. Bhaskar [88] là dưới 10%. Như vậy, các kết quả tính toán này và tài liệu công bố là tin cậy và chấp nhận được. Bài toán 2: (theo bài toán Khare-Kant [54]) Bảng 2.3. Chuyển vị không thứ nguyên của các vỏ composite lớp Khare-Kant R/a CST FSDT HSDT Luận án [54] 5 1.1280 1.1248 1.1235 1.1261 1.1263 [0/900] 10 1.1447 1.1439 1.1421 1.1434 1.1443 50 1.1501 1.1501 1.1482 1.1493 1.1502 5 1.0247 1.0215 1.0216 - 1.0151 0 [0/90 ]5 10 1.0310 1.0302 1.0303 - 1.0225 50 1.0331 1.0330 1.0332 - 1.0250 *Nhận xét: - Trong bảng 2.3 thể hiện kết quả chuyển vị không thứ nguyên cho vỏ trụ composite có số lớp khác nhau dưới tác dụng của tải nhiệt. So sánh kết quả tính toán trong tài liệu này với kết quả tính theo HOST12 bởi Khare & Kant [54] là phù hợp.
14. 12 Từ đó, có thể kết luận kết quả tính toán của luận án so với kết quả công bố trong các tài liệu [88] và [54] là tương đồng. 2.6. Kết luận chƣơng 2 Trong Chương 2, luận án đã đạt được một số kết quả cụ thể sau: - Thiết lập được mô hình toán và quan hệ ứng xử cơ học cho vỏ trụ composite lớp; xây dựng thuật toán theo hướng giải tích và lập chương trình để tính toán kết cấu chịu tác dụng đồng thời tải trọng cơ-nhiệt. - Xây dựng chương trình tính toán trạng thái ứng suất vỏ trụ composite lớp trong môi trường Maple và được kiểm chứng bằng việc so sánh kết quả tính toán với các công bố, đảm bảo tin cậy. Chƣơng 3. NGHIÊN CỨU TRẠNG THÁI ỨNG SUẤT TẠI VÙNG ỨNG SUẤT MẠNH CỦA VỎ TRỤ COMPOSITE LỚP 3.1. Phân tích dạng nghiệm của hệ phƣơng trình cân bằng vỏ Đối với trường hợp vỏ làm từ vật liệu composite lớp dạng cross- ply hệ phương trình cân bằng có dạng (2.20) như đã trình bày trong Chương 2. Trạng thái ứng suất biến dạng của vỏ trụ composite phụ thuộc vào dạng nghiệm của đa thức đặc trưng hệ (2.20). Đa thức đặc trưng của hệ này cũng chính là đa thức đặc trưng của hệ phương trình đại số (2.34) với biến số p. Dạng nghiệm của đa thức đặc trưng này bao gồm các nhóm nghiệm tuân theo bất đẳng thức chặt. Mỗi nhóm nghiệm thể hiện một dạng của trạng thái ứng suất của vỏ. Thông thường, TTƯSBD được chia ra thành: TTƯSBD cơ bản và các hiệu ứng biên. Vì vậy, để phân tích TTƯSBD ta đi nghiên cứu các nghiệm đa thức đặc trưng của hệ phương trình (2.34) trong một số trường hợp cụ thể. Để làm ví dụ minh họa, dưới đây chỉ xét trong trường hợp vỏ chịu tác dụng của tải đối xứng trục, tức là khi xét hệ phương trình (2.23) đối với các giá trị cụ thể của K như sau: Đối với mô hình “K=2” Khi đó, hệ phương trình đặc trưng (2.34) có dạng: ( 2 ) t A .[U00 ,U 10 ,U 20 ,W 00 ,W 10 ] = 0 , (3.1a) Đa thức đặc trưng của hệ (3.1a) trở thành 4 2n 2 2n 1n p 1 H p 0. (3.2) n0 Ngoài nghiệm tầm thường, đa thức đặc trưng (3.2) còn có các trường hợp nghiệm sau:
15. 13 - Trường hợp 1: có một cặp nghiệm phức liên hợp ±p1±iq1 và hai cặp nghiệm thực ±p2,±p3 với bất đẳng thức chặt p1<<p2, p3. Khi đó, nghiệm tổng quát của hệ (3.1) được biểu diễn dưới dạng: pp11xx F11 a 1 sin q 1x a 2 cos q 1 x e a 3 sin q 1 x a 4 cos q 1 x e (3.3) pp22xx pp33xx a5 e a 6 e a 7 e a 8 e - Trường hợp 2: có hai cặp nghiệm phức liên hợp ±p1±iq1, ±p2±iq2 với bất đẳng thức chặt p1<<p2, q1<<q2. Trong trường hợp này, nghiệm tổng quát của hệ (3.1) được viết như sau: F a sinqx a cos q x epp11xx a sinq x a cos q x e 12 1 1 2 1 3 1 4 1 (3.4) a sinqx a cos q x epp22xx a sinq x a cos q x e 5 2 6 2 7 2 8 2 Đối với mô hình “K=3” Đa thức đặc trưng của hệ trở thành 6 n 3 2n 2n  1 H p 0. (3.5) n0 Ngoài nghiệm tầm thường, đa thức đặc trưng (3.5) còn có các trường hợp nghiệm sau: - Trường hợp 1: có một cặp nghiệm phức liên hợp ±p1±iq1 và ba cặp nghiệm thực ±p3, ±p4, ±p5 với bất đẳng thức chặt sau p1<<p2, p3, p4, p5. Ở đây, nghiệm tổng quát của hệ phương trình (3.1) được biểu diễn dưới dạng: F a sinqx a cosq x ep1x a sinq x a cosq x e p 1 x a e p 2 x 211121 3141 5 (3.6) ae p2x aep3x ae p 3 x ae p 4 x ae p 4 x ae p 5 x ae p 5 x 6 7 8 9 10 11 12 - Trường hợp 2: có hai cặp nghiệm phức liên hợp ±p1±iq1, ±p2±iq2 và một cặp nghiệm thực ±p3 với bất đẳng thức chặt p1<<p2, p3, q1<<q2. Khi đó, nghiệm tổng quát của hệ phương trình (3.1) được biểu diễn dưới dạng: pp11xx Faqaqe22 1sin 1x 2 cos 1 x aqaqe 3 sin 1 x 4 cos 1 x pp22xx a5sin q 2x a 6 cos q 2 x e a 7 sin q 2 x a 8 cos q 2 x e (3.7) a epp33xx a e a epp44xx a e 9 10 11 12 Đa thức đặc trưng chia ra làm hai nhóm: nghiệm nhỏ và nghiệm lớn. Nghiệm nhỏ đặc trưng cho TTƯSBD cơ bản của vỏ, TTƯSBD này có thể được xác định bởi các phương trình của lý thuyết vỏ Mindlin. Còn các nghiệm lớn đặc trưng cho các hiệu ứng biên của vỏ. Các hiệu ứng này chỉ tồn tại ở vùng trạng thái kết cấu vỏ bị biến hình (lân cận vùng biên, vùng chịu tác dụng của lực tập trung ), khi đi xa vùng trạng thái kết cấu vỏ bị biến hình thì hiệu ứng của
16. 14 chúng nhanh chóng bị triệt tiêu. 3.2. Ảnh hƣởng của một số tham số tới hiện tƣợng gia tăng ứng suất tại vùng ứng suất mạnh - Sự phù hợp kết quả tính toán giữa các lý thuyết khác nhau tại vị trí giữa vỏ có thể giải thích dựa vào giá trị nghiệm của đa thức đặc trưng (3.1). Ở vị trí xa biên trạng thái ứng suất được xác định chủ yếu bởi các nghiệm bé p11 iq của đa thức đặc trưng (3.1), còn các nghiệm lớn chỉ gây ảnh hưởng lớn tại những khu vực lân cận vùng trạng thái ứng suất biến hình. Giá trị các nghiệm nhỏ theo cả hai mô hình K=2 và K=3 khá gần nhau và xấp xỉ bằng nghiệm được xác định bởi lý thuyết FSDT. - Có sự sai khác lớn về trạng thái ứng suất khi tính toán với lý thuyết khác nhau. Ở đây, lý thuyết biến dạng trượt bậc cao đã chỉ ra có hiện tượng tập trung ứng suất tại biên. a1). Ứng suất x tại biên a2). Ứng suất x tại xx 0 2 b1). Ứng suất q tại biên b2). Ứng suất q tại xx 0 /2
17. 15   xx 2 c1). Ứng suất z tại biên c2). Ứng suất z tại 0   d1). Ứng suất x z tại biên d2). Ứng suất x z tại x 5/hR Hình 3.1. Sự biến đổi của ứng suất không thứ nguyên theo chiều dày vỏ theo các lý thuyết khác nhau L=4R, S=10 3.2. Ảnh hƣởng của một số tham số tới hiện tƣợng gia tăng ứng suất tại vùng ứng suất mạnh *Ảnh hưởng của tham số kết cấu Khi xét vỏ trụ composite lớp với góc xếp lớp và thứ tự xếp lớp khác nhau, nhận thấy: - Thứ tự và góc xếp lớp ảnh hưởng đến trạng thái ứng suất của vỏ trụ; - Tại những vị trí gần biên thì ứng suất cực đại của vỏ trụ thay đổi tùy thuộc vào thứ tự và góc xếp lớp.
18. 16 *Ảnh hưởng của điều kiện biên - Ở những vị trí xa biên thì điều kiện biên không ảnh hưởng nhiều đến chuyển vị và trạng thái ứng suất của vỏ trụ; - Tại những vị trí gần biên thì ứng suất cực đại của vỏ trụ thay đổi tùy thuộc vào điều kiện biên; - Khi vỏ trụ càng ngắn thì điều kiện biên càng ảnh hưởng nhiều tới chuyển vị lớn nhất tại điểm giữa vỏ trụ. Như vậy, kết cấu sẽ dễ bị phá hủy bởi hiện tượng tăng lên của một số thành phần ứng suất tại biên nên trong tính toán, thiết kế, chế tạo cần hạn chế sử dụng hoặc phải có biện pháp kỹ thuật khi sử dụng các kết cấu vỏ có điều kiện biên ngàm chặt. Nếu phải sử dụng kiểu liên kết này thì cần phải kiểm tra bền và kiểm tra các khuyết tật. Bảng 3.6. Ứng suất vỏ với điều kiện biên khác nhau của vỏ trụ composite lớp, (L/R=4, S=10) x q  x z  z Điều kiện q q q q x, x, x, x, biên 2 2 2 2 h /2 h /2 h /3 h /4 0.7720 0.3472 -2.6158 -1.3210 C-C -2.9086 -0.1117 3.7355 3.7998 0.0000 0.0000 -2.6184 0.7441 C-S 0.0000 0.0000 2.1405 1.5500 0.0059 -1.7978 -0.3876 0.5432 C-F 0.0137 -0.0824 -0.0023 0.0973 0.0000 0.0000 -2.6184 0.7441 S-S 0.0000 0.0000 2.1405 1.5500 *Ảnh hưởng của tải trọng - Khi vỏ chịu tác dụng tải trọng phân bố tuyến tính theo chiều dài thì chuyển vị của phần vỏ trụ xa hai biên cũng có dạng tuyến tính tương tự. Trạng thái ứng suất của vỏ trụ cũng tương ứng có thay đổi theo tải trọng. - Khi tải trọng giảm dần về biên thì hiện tượng tập trung ứng suất tại biên cũng giảm đáng kể điều này cho ta ý nghĩa thực tiễn khi thiết kế chế tạo phải có biện pháp làm giảm tải trọng tại vùng biên ngàm. 3.3. Kết luận chƣơng 3 Trong chương 3, luận án đã thực hiện nghiên cứu trạng thái ứng suất của vỏ trụ composite lớp tại vùng ứng suất mạnh và khảo sát- đánh giá ảnh hưởng của điều kiện biên, thông số kết cấu và một số
19. 17 dạng tải trọng đến trạng thái ứng suất tại khu vực này. Một số nhận xét chính được rút ra như sau: - Nghiệm của đa thức đặc trưng của hệ phương trình cân bằng vỏ được chia thành lớp nghiệm nhỏ và lớp nghiệm lớn: lớp nghiệm nhỏ đặc trưng cho trạng thái ứng suất tại vùng xa biên, còn lớp nghiệm lớn có ảnh hưởng lớn đến trạng thái ứng suất tại vùng biên. - Tại vùng biên trạng thái ứng suất chịu sự ảnh hưởng lớn của các hiệu ứng biên, tương ứng với các giá trị nghiệm lớn, vì vậy khi xác định trạng thái ứng suất tại vùng biên cần thiết phải sử dụng lý thuyết vỏ bậc cao. - Vùng ứng suất mạnh xảy ra rõ nét đối với khu vực biên ngàm. Phạm vi ảnh hưởng của hiện tượng này khá nhỏ với kích thước bằng chiều dày vỏ. Khi vỏ chịu tải cơ-nhiệt thì hiện tượng này xảy ra với cả biên ngàm và tựa. - Tại những khu vực xa biên, kết quả tính toán cho vỏ mỏng theo lý thuyết biến dạng trượt bậc cao và lý thuyết bậc nhất, lý thuyết cổ điển có sai lệch không đáng kể. Điều này có nghĩa là tại những vùng ở xa biên các tính toán theo lý thuyết đàn hồi 3D, lý thuyết bậc nhất, cổ điển và công trình này là tương đồng. Chƣơng 4. NGHIÊN CỨU ẢNH HƢỞNG CỦA THAM SỐ KẾT CẤU VÀ TẢI TRỌNG ĐẾN TRẠNG THÁI ỨNG SUẤT 4.1. Nghiên cứu ảnh hƣởng của tham số kết cấu *Ảnh hưởng của chiều dày vỏ (tỉ số R/h) - Độ dày của vỏ ảnh hưởng nhiều đến chuyển vị và ứng suất của vỏ; - Hiện tượng tập trung ứng suất tại biên xảy ra đối với cả vỏ mỏng và dày. Ứng suất tập trung này làm tăng khoảng 1,4 lần ứng suất tối đa; - Vỏ có độ dày càng lớn thì một số thành phần ứng suất tại vùng tập trung ứng suất càng tăng mạnh. *Ảnh hưởng của chiều dài vỏ (tỷ số L/R) - Khi chiều dài tương đối nhỏ L/R 2 (vỏ dài) càng tăng thì điều kiện biên biên ảnh hưởng không nhiều tới chuyển vị và ứng suất cực đại của vỏ. - Đặc biệt khi vỏ cực ngắn (L/R<1) thì chuyển vị và ứng suất của
20. 18 vỏ có sự khác biệt rõ rệt. Chuyển vị của vỏ giảm đi một nửa so với các vỏ cùng điều kiện khi chiều dài tương đối L/R=0,5. *Ảnh hưởng của số lớp, thứ tự xếp lớp - Vỏ trụ có cùng chiều dày thì số lớp có ảnh hưởng đáng kể đến chuyển vị và ứng suất. Số lớp càng ít thì chuyển vị tại vị trí giữa vỏ trụ càng nhiều hơn. - Số lớp càng nhiều thì ứng suất trượt  * là lớn nhất tại biên x z ngàm và có thể lớn hơn ứng suất chính. Như vậy, vỏ nhiều lớp dễ gây ra trạng thái bong lớp tại vị trí có bước nhảy về kết cấu hoặc lực. Nên khi tính toán thiết kế cần lưu ý lựa chọn số lớp cho phù hợp. * * a1) Ứng suất x tại vị trí biên a2) Ứng suất x tại xx 0 /2 * * b1) Ứng suất  tại vị trí biên b2) Ứng suất  tại xx /2 q q 0
21. 19 * * c1) Ứng suất  tại vị trí biên c2) Ứng suất  tại xx /2 z z 0 d1) Ứng suất  * tại vị trí biên d2) ) Ứng suất  * tại x 5/hR x z x z Hình 4.4. Ảnh hưởng số lớp đến ứng suất tại một số vị trí của vỏ * Ảnh hưởng của góc xếp lớp - Với vỏ trụ có cùng chiều dày và số lớp như nhau thì thứ tự xếp lớp có ảnh hưởng không nhiều đến chuyển vị. Tuy nhiên nó lại làm thay đổi rất nhiều tới trạng ứng suất. - Tại các vùng xa biên thì tuy trạng thái ứng suất của vỏ có thứ tự xếp lớp là khác nhau nhưng giá trị lớn nhất lại không có sự khác biệt. Điều này khẳng định rằng tại các vùng xa biên hoặc xa chỗ tập trung lực thì trạng thái ứng suất gần như không bị ảnh hưởng bởi trật tự xếp lớp. Như vậy, trật tự xếp lớp làm ảnh hưởng tới hiệu ứng tập trung ứng suất của vỏ trụ composite lớp, nên khi tính toán thiết kế cần phải kiểm tra sơ đồ xếp lớp hợp lý để tăng độ bền nhờ hạn chế hiện tượng phá hủy do hiện tượng tập trung ứng suất tại khu vực liên kết hoặc
22. 20 tải tập trung hoặc là kết cấu có bước nhảy. 4.2. Nghiên cứu ảnh hƣởng của một số dạng tải trọng cơ Khảo sát ảnh hưởng của các dạng tải trọng (không kể tới tác động của nhiệt độ) đến trạng thái ứng suất biến dạng của vỏ trụ composte lớp (vật liệu Gr.Ep-T, [900/0/900], L/R=4; S=[10 50]; h=0.1m), vỏ hai đầu ngàm, trong các trường hợp tải trọng cụ thể sau: - Tải trọng phân bố theo quy luật tuyến tính. + Khi vỏ chịu tác dụng tải trọng phân bố tuyến tính theo chiều dài thì chuyển vị của phần vỏ trụ xa hai biên cũng có dạng tuyến tính tương tự. Trạng thái ứng suất của vỏ trụ cũng tương ứng có thay đổi theo tải trọng. + Khi tải trọng giảm dần về biên thì hiện tượng tập trung ứng suất tại biên cũng giảm đáng kể. - Tải trọng phân bố theo quy luật tuần hoàn. + Khi vỏ chịu tác dụng tải trọng phân bố tuần hoàn theo chiều dài thì chuyển vị của phần vỏ trụ xa hai biên (khoảng >5h) cũng có dạng tuần hoàn tương tự hình 2.6, đồng thời khi đó, trạng thái ứng suất của vỏ trụ tương ứng cũng có thay đổi theo tải trọng. + Khi tải trọng giảm dần về biên thì hiện tượng tập trung ứng suất tại biên cũng giảm đáng kể điều này cho ta ý nghĩa thực tiễn khi thiết kế, chế tạo các kết cấu vỏ sao cho lực tác động nhỏ tại vùng biên như dạng tải trọng này. - Tải trọng phân bố theo quy luật hàm mũ Khi vỏ chịu tác dụng tải trọng phân bố theo quy luật hàm mũ dọc chiều dài thì chuyển vị của phần vỏ trụ xa hai biên (khoảng >5h) cũng có hình dạng tương tự hình 4.11. Trạng thái ứng suất tương ứng của vỏ cũng có thay đổi theo biên độ tải trọng tác dụng lên vỏ. - Tải trọng phân bố theo quy luật hàm đa thức - Khi vỏ chịu tác dụng tải trọng áp suất phân bố theo quy luật hàm đa thức dọc chiều dài vỏ thì chuyển vị của phần vỏ trụ xa hai biên cũng có hình dạng tương tự. - Ngoài trừ hai đầu biên, thì vị trí nào vỏ có biên độ lực lớn thì chuyển vị và ứng suất của vỏ cũng lớn hơn vùng khác. 4.3. Nghiên cứu ảnh hƣởng nhiệt độ Xét vỏ trụ composite lớp chịu tác dụng nhiệt độ với: thông số hình học (R=1; L/R=[0,5 4]; S=[10 20]); vật liệu Gr.Ep_AS, theo bảng 2.1; Xét vỏ trong các vỏ có điều kiện biên như sau: vỏ hai đầu
23. 21 ngàm, vỏ một đầu ngàm- một đầu tựa; vỏ một đầu ngàm - một đầu tự do và vỏ hai đầu tựa đơn. a1) Ứng suất không thứ nguyên a2) Ứng suất không thứ nguyên x tại vị trí biên xx 0 x tại xx 0 2 b1) Ứng suất không thứ nguyên b2) Ứng suất không thứ nguyên q tại vị trí biên xx 0 q tại xx 0 2 c1) Ứng suất không thứ nguyên c2) Ứng suất không thứ nguyên  z tại vị trí biên xx 0  z tại xx 0 2
24. 22 d1) Ứng suất không thứ nguyên x z d2) ) Ứng suất không thứ nguyên x z tại xx 0 2 tại vị trí biên xx 0 Hình 4.15. Biểu đồ phân bố ứng suất không thứ nguyên trong vỏ trụ composite 2 lớp [0/900] với LRS 0.5 , 20 *Nhận xét: Phân tích các kết quả thu được trong các trường hợp đối với vỏ chịu tác dụng của tải nhiệt độ cho thấy: - Đối với vỏ có chiều dài, điều kiện biên không ảnh hưởng nhiều đến chuyển vị và ứng suất tại vị trí xa biên. - Tại các cạnh liên kết của vỏ, ứng suất lớn nhất phụ thuộc vào điều kiện biên. Trong các trường hợp đã xét, ứng suất tại biên ngàm có giá trị lớn nhất, tại đó xuất hiện ứng suất tập trung. - Do ảnh hưởng của hiệu ứng biên, các thành phần ứng suất tăng lên nhiều so với ứng suất tại vị trí giữa của vỏ, trong đó các thành phần của ứng suất cắt tăng mạnh hơn. - Đối với vỏ ngắn, điều kiện biên có ảnh hưởng lớn đến chuyển vị và các thành phần ứng suất trong toàn bộ kết cấu. 4.4. Nghiên cứu TTƢS vỏ dƣới tác dụng của tải cơ-nhiệt - Nhiệt độ làm tăng ứng suất của kết cấu vỏ. - Nhiệt độ càng tăng hiệu ứng tập trung ứng suất tại biên. Trong * * đó, ứng pháp  z và suất tiếp x z tại lớp giáp lớp vỏ thành ngoài mạnh hơn vùng lớp thành trong vỏ. 4.5. Kết luận chƣơng 4 Từ việc khảo sát ảnh hưởng của các thông số kết cấu, vật liệu và tải trọng tới trạng thái ứng suất kết cấu vỏ, rút ra một số kết luận sau: - Đã thực hiện khảo sát ảnh hưởng của các điều kiện biên, độ dày tương đối (R/h), chiều dài tương đối (L/R) của vỏ đối với chuyển vị và ứng suất.
25. 23 - Tại vùng biên trạng thái ứng suất chịu sự ảnh hưởng lớn của các hiệu ứng biên, tương ứng với các giá trị nghiệm lớn, vì vậy khi xác định trạng thái ứng suất tại vùng biên cần thiết phải sử dụng lý thuyết vỏ bậc cao. - Khi chiều dài tương đối của vỏ lớn thì chuyển vị và ứng suất tại các vùng xa biên không ảnh hưởng nhiều bởi điều kiện biên . Điều này có nghĩa là tại những vùng ở xa biên các tính toán theo lý thuyết đàn hồi 3D, lý thuyết Mindlin và công trình này là tương đồng. Qua khảo sát ảnh hưởng thấy rằng các thông số kết cấu ảnh hưởng nhiều đến trạng thái ứng suất-biến dạng của vỏ composite lớp nên chế tạo, nghiệm thu cần đặc biệt chú ý: - Đối với vỏ nhiều lớp cần lưu ý khả năng vỏ bị phá hủy do hiện tượng tách lớp vật liệu vì có hiện tượng tập trung ứng suất trượt ngang khi số lớp của vỏ tăng lên, nên cần chú ý các biện pháp làm giảm hiệu ứng này như làm tăng độ dày tại vùng liên kết hoặc là chuyển điều kiện biên từ ngàm sang tựa hoặc tựa đàn hồi. - Đối với vỏ ít lớp thì cần tăng cứng tại các vùng xa biên để tăng khả năng chịu tải hay là tăng hiệu quả của toàn kết cấu vỏ. KẾT LUẬN I. Các kết quả chính đạt đƣợc Các nội dung chính của luận án đã thể hiện được: 1. Luận án đã nghiên cứu tổng quan về vật liệu composite, các đặc tính cơ học, công nghệ chế tạo, các ứng dụng của nó trong khoa học kỹ thuật và các mô hình tính toán thường được sử dụng trong nghiên cứu kết cấu vỏ trụ composite lớp. 2. Xây dựng được mô hình toán học tính toán vỏ trụ composite lớp chịu tác dụng đồng thời tải trọng cơ-nhiệt trên cơ sở lý thuyết biến dạng trượt bậc cao kiểu Quasi-3D. Độ tin cậy của mô hình toán học, phương pháp nghiên cứu và chương trình tính toán được kiểm chứng thông qua so sánh với các các kết quả nghiên cứu đã công bố trên các tạp chí quốc tế uy tín. 3. Luận án đã thực hiện nghiên cứu ảnh hưởng của các thông số hình học (, các thông số vật liệu, số lớp, trật tự xếp lớp, các điều kiện biên và dạng tải trọng đến trạng thái ứng suất- biến dạng của vỏ trụ composite lớp. Kết quả nghiên cứu nhận được có ý nghĩa rất quan trọng trong thực tiễn thiết kế, chế tạo, cũng như khai thác kết cấu vỏ trụ composite lớp. 4. Dựa vào phân tích nghiệm của đa thức đặc trưng, luận án đã chỉ ra hiện sự tồn tại của vùng ứng suất mạnh (khu vực biên ngàm chặt,
26. 24 khu vực chịu tải tập trung, ), tại đây trạng thái ứng suất có sự gia tăng mạnh. Sự suất hiện của vùng ứng suất mạnh có thể được giải thích do sự tồn tại của các hiệu ứng biên bổ sung tương ứng với các giá trị nghiệm lớn của đa thức đặc trưng. Đồng thời, trong luận án cũng chỉ ra giá trị tương đối lớn của ứng suất  z (thường bị bỏ qua trong các lý thuyết cổ điển, lý thuyết biến dạng trượt bậc nhất) tại những vùng trạng thái ứng suất bị suy biến là không thể bỏ qua. II. Những đóng góp mới của luận án 1. Dựa trên lý thuyết biến dạng trượt bậc cao kiểu quasi-3D, phương pháp tiếp cận giải tích đã xây dựng được mô hình và chương trình toán để xác định trạng thái ứng suất-biến dạng của vỏ trụ composite lớp chịu tải áp suất trong và nhiệt độ với các điều kiện biên khác nhau. 2. Thông qua phân tích hiệu ứng biến đổi xác định các vùng ứng suất mạnh trong kết cấu vỏ, đặc biệt là xác định được các thành phần ứng suất pháp và ứng suất tiếp theo chiều dày (z,,  yz  xz ) mà một số lý thuyết khác chưa đề cập tới. 3. Đánh giá ảnh hưởng của điều kiện biên, tải trọng và các thông số kết cấu đến trạng thái ứng suất-biến dạng của vỏ trụ composite lớp để làm cơ sở cho việc chọn và thiết kế kết cấu vỏ hợp lý. III. Những vấn đề cần tiếp tục nghiên cứu 1. Nghiên cứu kết cấu vỏ composite lớp chịu các loại tải trọng cơ nhiệt khi tính đến bài toán phi tuyến. 2. Tối ưu hóa kết cấu vỏ trụ composite lớp nhằm giảm sự gia tăng một số thành phần ứng suất tại vùng ứng suất mạnh khi chịu tác động đồng thời của nhiều loại tải trọng phức tạp khác nhau.
27. 25 DANH MỤC CÁC CÔNG TRÌNH KHOA HỌC ĐÃ CÔNG BỐ CT1. Trần Ngọc Đoàn, Phan Văn Chương, Nguyễn Trường Thanh, (2017) “Xây dựng mô hình tính toán vỏ trụ composite lướp theo lý thuyết biến dạng trượt bậc cao có tính tới ảnh hưởng của biến dạng pháp tuyến”, Hội nghị Cơ học toàn quốc lần thứ X, Tập 3, tr.260- 267. CT2. Nguyễn Trường Thanh, Trần Ngọc Đoàn, Nguyễn Việt Hà, Trần Văn Hùng, Hoàng Mai Phương, (2018), “Nghiên cứu trạng thái ứng suất, biến dạng vỏ trụ composite lớp với các dạng tải trọng khác nhau trên cơ sở lý thuyết biến dạng trượt bậc cao kiểu quasi-3D”, Hội nghị khoa học toàn quốc Cơ học vật rắn lần thứ XIV,tr 578-585. CT3. Tran Ngoc Doan, Do Van Thom, Nguyen Truong Thanh, Phan Van Chuong, Nguyen Chi Tho, Nguyen Tri Ta, Hoang Nam Nguyen (2019), “Analysis of stress concentration phenomenon of cylinder laminated shells using higher-order shear deformation Quasi-3D theory”. Composite Structures, 232 111526. CT4. Nguyễn Trường Thanh, Phan Văn Chương, Lê Song Tùng, Trần Ngọc Đoàn và Trần Văn Hùng (2020), “Tính toán vỏ trụ composite lớp trên cơ sở lý thuyết biến dạng trượt bậc cao Quasi-3D theo hướng tiếp cận giải tích”, Tạp chí nghiên cứu Khoa học và Công nghệ quân sự, Số 65/02-2020, tr.186-195. CT5. Tran Ngoc Doan, Nguyen Truong Thanh (2021), “Thermoelastic response and boundary effect of cross-ply laminated cylindrical shells based on a quasi-3D type higher-order shear deformation theory”, International Journal of Pressure Vessels and Piping, Vol 194, Part A, 12-2021, 104534. CT6. Nguyễn Trường Thanh, Trần Ngọc Đoàn, Phan Văn Chương (2021), “Phân tích đáp ứng cơ-nhiệt của vỏ trụ composite lớp sử dụng lý thuyết biến dạng trượt bậc cao kiểu quasi 3D”, Tạp chí Nghiên cứu Khoa học và Công nghệ quân sự, (Số Đặc san Viện Tên lửa, 10-2021), tr.105-114. CT7. Nguyễn Trường Thanh (2021), “Nghiên cứu trạng thái ứng suất nhiệt của vỏ trụ composite lớp trên cơ sở lý thuyết biến dạng trượt bậc cao theo hướng tiếp cận giải tích”, Tạp chí Nghiên cứu Khoa học và Công nghệ quân sự, (Số Đặc san HNKH, 11-2021), tr.85-93.