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

3D theo hướng tiếp cận giải tích. Trong bài báo thực hiện sự thiết lập các phương trình ba chiều trong lý thuyết đàn hồi phi tuyến thành các phương trình phi tuyến hai chiều đối với vỏ trụ bằng cách sử dụng phương pháp biến phân và phân tích trường chuyển vị thành chuỗi hàm đa thức theo chiều dày vỏ. Trên cơ sở các phương trình nhận được đã đưa ra hệ phương trình cân bằng theo trường chuyển vị và các điều kiện biên tương ứng đối với trường hợp phân tích chuyển vị thành đa thức bậc ba. Từ đó khảo sát các trạng thái ứng suất và biến dạng của vỏ trụ composite lớp với các tham số khác nhau trong điều kiện liên kết ngàm hai đầu và đưa ra phạm vi ứng dụng cho từng trường hợp. Từ khóa: vỏ trụ composite lớp, trạng thái ứng suất biến dạng, biến dạng trượt bậc cao, chuyển vị. 1. ĐẶT VẤN ĐỀ Trong tính toán kết cấu thành mỏng dạng tấm vỏ, dựa trên các giả thiết của Kirchhoff– Love, nhiều tác giả đã xây dựng các lý thuyết tấm vỏ tuyến tính hai chiều như Timoshenko và Woinowsky-Krieger [1], Flügge [3], Các lý thuyết nêu trên là các phương án khác nhau của lý thuyết cổ điển (CST). Nói chung, lý thuyết tấm vỏ cổ điển không cho kết quả phù hợp khi áp dụng để tính toán tấm, vỏ dày hoặc làm từ vật liệu composite lớp. Để khắc phục những nhược điểm của lý thuyết cổ điển, Reissner đề xuất lý thuyết biến dạng trượt bậc nhất (FSDT) [4]. Tuy nhiên, cũng giống như CST, trong FDST đòi hỏi phải đưa thêm vào hệ số hiệu chỉnh cắt để hiệu chỉnh kết quả tính toán ứng suất cắt. Nhiều nghiên cứu về tấm\vỏ composite lớp có sử dụng lý thuyết biến dạng cắt bậc cao, như mô hình lý thuyết bậc 3 của Ready [5], Aydogdu [6], Asadi [7], Kant [8],... Mô hình lý thuyết biến dạng cắt bậc cao đã khắc phục được nhược điểm của lý thuyết cổ điển CST và bậc nhất FSTD là không cần đưa vào hệ số hiệu chỉnh cắt, trong đó ứng suất cắt phân bố theo chiều dày là hàm phi tuyến. Tuy nhiên, cần phải lưu ý rằng, các giả thiết mà các tác giả sử dụng đã không thỏa mãn điều kiện biên tại bề mặt trên và dưới của vỏ. Do đó, trong tính toán cần phải thực hiện quá trình xác định lại ứng suất và biến dạng của vỏ tại những vị trí đặc biệt của vỏ. Trong bài báo này, nghiên cứu được thực hiện trên cơ sở mô hình lý thuyết biến dạng cắt bậc cao có tính đến ảnh hưởng của biến dạng và ứng suất pháp ngang trình bày trong [9-10]. Trường chuyển vị trong trường hợp này thỏa mãn các điều kiện tương thích về năng lượng được đề xuất bởi Vasilev và Lurie [12]. Nguyên lý biến phân Lagrange được sử dụng để thiết lập các phương trình cân bằng và điều kiện biên tương ứng. Phân tích trạng thái ứng suất của vỏ trụ composite cross-ply được thực hiện dựa trên cơ sở phép biến đổi Laplace. Các thành phần ứng suất ngang được xác đinh dựa trên các phương trình của lý thuyết đàn hồi 3D nên trường chuyển vị, trường biến dạng và ứng suất được xác định hoàn toàn. Điều này, khẳng đinh qua các kết quả tính toán số nhận được phù hợp với kết quả tính toán theo lý thuyết đàn hồi 3D. 2. THIẾT LẬP MÔ HÌNH VÀ XÂY DỰNG HỆ PHƯƠNG TRÌNH TÍNH Nghiên cứu khoa học công nghệ Tạp chí Nghiên cứu KH&CN quân sự, Số 65, 02 - 2020 187 2.1. Thiết lập các phương trình cơ bản 2.1.1. Mô hình chuyển vị Xét vỏ trụ composite lớp trong hệ tọa độ trụ .O z Kết cấu vỏ trụ composite gồm n lớp, có tổng chiều dày bằng ,h mỗi lớp là vật liệu composite cốt sợi đồng phương với các thông số hình học như hình 1. x z q R L uw v x q 31 z 2 b h 2 h 2 1 2 i n Z1 Zi Z0 Z q Hình 1. Vỏ trụ composite lớp và các hệ trục tọa độ. Chuyển vị của vỏ trụ trong hệ tọa độ trụ O z có dạng:             1 0 0 0 , , , , , , , , , , , ! ! ! K K Kk k k k k k k k k z z z u z u v z v w z w k k k                      (1) Trong đó:  , ,u z  ,  , ,v z  ,  , ,w z  là các thành phần chuyển vị ba chiều của điểm  , ,P z  ở khoảng cách z tính từ mặt trung bình (z=0) theo các trục tọa độ. 0u , 0 ,v 0w là các thành phần chuyển vị 2D của điểm  , ,0P   theo các trục tọa độ. 1,u 1v là góc quay so với mặt trung bình theo các trục ,  tương ứng. Các thành phần khác trong công thức Error! Reference source not found. là thành phần chuyển vị 2D bậc cao trong phân tích chuỗi Taylor. 2.1.2. Biến dạng và ứng suất vỏ composite Quan hệ tuyến tính giữa các thành phần biến dạng và chuyển vị: 1 1 1 1 , , , 1 1 1 , , .z z z u v v u w R R z R R z w v v w u w R z z R z R z R z                                                     (2) Khi vỏ trụ làm việc trong giới hạn đàn hồi thì mối liên hệ giữa biến dạng và ứng suất tuân theo định luật Hooke đối với lớp vật liệu thứ k trong hệ tọa độ địa phương gắn với từng lớp có dạng:        ,k k kC     (3) Trong đó:   ( ) 11 22 33 12 23 13, , , , , Tk k           ( ) 11 22 33 12 23 13, , , , , Tk k        Phương trình (3) là mối quan hệ giữa trường ứng suất và trường biến dạng của lớp k Cơ kỹ thuật & Kỹ thuật cơ khí động lực N. T. Thanh, , T. V. Hùng, “Tính toán vỏ trụ composite lớp hướng tiếp cận giải tích.” 188 trong hệ tọa độ cục bộ của lớp vật liệu. Nhưng hệ tọa độ cục bộ này không trùng với hệ tọa độ chung nên cần thực hiện chuyển trường chuyển vị và ứng suất sang hệ tọa độ chung. Mối quan hệ giữa trường ứng suất và chuyển vị trong hệ tọa độ chung theo công thức:      ( ) ( ) ,kk kQ      (4) ở đây:    ( ) ( ) ( ) ( ) ( ) ( ) ( ), , , , , T k k k k k k k z z z           là trường ứng suất của vỏ trong hệ tọa độ chung;    ( ) ( ) ( ) ( ) ( ) ( ) ( ), , , , , T k k k k k k k z z z           là trường chuyển vị trong hệ tọa độ chung; ma trận độ cứng ( )kQ   được cho bởi: ( ) ( ) ( ) ( ) ; tk k k kQ T C T               ma trận độ cứng  kC    và ma trận chuyển hệ tọa độ ( ) ,kT   xác định theo [13]. 2.1.3. Phương trình cân bằng và các điều kiện biên Biến phân thế năng biến dạng đàn hồi được xác định theo công thức sau:  zz zz z z z z V U dV                           ( ) ( 1) ( ) ( ) ( ) 2 ( ) ( ) ( ) 1 1 i i h k k kn z z k k k i z z z zh z R dz d d R                                                (5) Xét vỏ trụ composite lớp chịu tác dụng của tải trọng phân bố hướng tâm ở bề mặt bên ngoài q+ và tải trọng phân phối hướng tâm trên bề mặt bên trong q-. Khi đó biến phân của công ngoại lực tác dụng lên vỏ được xác định theo công thức sau: 21 1 2 2 h h A q w q w R d d R R                               (6) ở đây, 2 0 1 2 , 2 8 h h w w w w    2 0 1 2 2 8 h h w w w w    . Để xây dựng phương trình cân bằng ta sử dụng nguyên lý biến phân Lagrange: ( ) 0U A    (7) Thay (5) và (6) vào (7), thông qua biến đổi thu được hệ phương trình cân bằng theo chuyển vị như sau: 0, 0, N N NN Q                    0 0, 0, Q M MQ N Rp RQ                       10, 0,z M SM S RQ M RQ Rp                        (8) * * * * *0, 0, N N N N RS RS Q                          * * ** * * 2 0, 0,z Q M MQ N RS Rp RQ                        khi K=3, Nghiên cứu khoa học công nghệ Tạp chí Nghiên cứu KH&CN quân sự, Số 65, 02 - 2020 189 * * * *2 0, M M RQ S              khi K=3. Trong hệ phương trình (8), sử dụng các ký hiệu nội lực suy rộng như sau: ( ) ( 1) 2 3 * * ( ) 1 1 1 , 2 6 k k hn k k h z z z N M N M z dz R                        ( ) ( 1) 2 3 * * ( ) 1 1 , 2 6 k k hn k k h z z N M N M z dz                       ( ) ( 1) ( ) 1 1 1 , k k hn k z z z k h z Q S z dz R             (9) ( ) ( 1) 2 3 * * ( ) 1 1 , 2 6 k k hn k k h z z N M N M z dz                   ( ) ( 1) 2 3 * * ( ) 1 1 , 2 6 k k hn k z k h z z Q S Q S z dz                     1 1 , 0,1,2 2 22 ! 2 ! ii i i i hh h h p q q i R Ri i                  2.1.4. Phương pháp giải tích Xét vỏ trụ với điều kiện liên kết khác nhau: Tựa đơn Tại: 1 2, :   * * * *0, 0,N M N M N M N M               0,jw  với i=0,..K; j=0,..K-1; hoặc 0, 0, 0,i jN M v w     với i=0,..K; j=0,..K-1; (10) Ngàm chặt Tại: 1 2, :   0, 0, 0,i i ju v w   với i=0,..K; j=0,..K-1 (11) Tự do Tại: 1 2, :   * * * * * 0, 0, 0; N M N M N M N M N N M                       (12) Đặt chuyển vị và tải trọng dưới dạng khai triển chuỗi Fourie dạng: Cơ kỹ thuật & Kỹ thuật cơ khí động lực N. T. Thanh, , T. V. Hùng, “Tính toán vỏ trụ composite lớp hướng tiếp cận giải tích.” 190                                                 1 2 0 1 1 2 0 1 1 2 0 1 1 2 0 1 , cos sin , , cos sin , 0..3 , sin cos , 0..3 , cos sin , 0..2. m m m k km km k m k km km k m l lm lm l m q Q m Q m Q u U m U m U k v V m V m V k w W m W m W l                                                                       (13) Hệ phương trình cân bằng (8) cho ta hệ 3 2K  phương trình vi phân đạo hàm riêng đối với 3 2K  chuyển vị , , ,k k ku v w và có bậc bằng  2 3 2K  . Giải hệ phương trình trên, ta nhận được nghiêm tổng quát có chứa các hằng số tích phân, các hằng số này tìm được từ các điều kiện liên kết tại biên theo (10), (11) hoặc (12). Thay biểu thức nghiệm tìm được đối với chuyển vị vào các công thức (4) cho ta biến dạng của vỏ. Để tìm các ứng suất , , ,     sử dụng công thức định luật Hooke, các ứng suất cắt tìm được bằng cách tích phân phương trình cân bằng của lý thuyết đàn hồi ba chiều: /2 1 1 , z z h z dz R z R                          2 1 1 , 2( ) /2 zR z z dzz R RR z h                            (14) /2 1 / 2 1 . / 2 z z z z h z R h dz q R z R R h                               2.2. Kiểm tra mô hình và chương trình tính toán Xét vỏ trụ composite lớp có chiều dài tương đối L/R=4, độ dày tương đối S=R/H thay đổi, các lớp composite có chiều dày bằng nhau và có thông số vật liệu như sau: E1= 25E2, G23= 0.2E2, G13= G12 =0.5E2, hệ số Poisson 12=0.25; vỏ chịu tải hình sin tác dụng trên thành trong vỏ dạng 0 sin cosmm m R q Q n L   ; trong điều kiện liên kết tựa đơn hai đầu. Trong bảng 1 trình bày kết quả tính toán độ võng không thứ nguyên w tại vị trí giữa vỏ đối với các độ dày tương đối khác nhau và các kiểu phân bố lớp vật liệu khác nhau nhận được trong công trình này và so sánh với kết quả của Salvatore Brischetto [11]. Chuyển vị và ứng suất không thứ nguyên được sử dụng theo công thức sau:    1 4 2 0 00 0 10 . 1 , , , , , . z z z z E w w Q S QQ HS Q S                  (15) Bảng 1. Độ võng w không thứ nguyên theo tại vị trí giữa của vỏ composite lớp chịu tải trọng hình sin bên trong ( / ; 1; 4/ 4;L mR S R h n    ). S [90] [90/0] [90/0/90] [90/0/90/0/90]2 TT* B* [11] TT B* [11] TT* B* [11] TT* B* [11] 4 2.7723 2.7830 7.3555 6.1000 3.8021 4.0090 3.9979 4.2060 Nghiên cứu khoa học công nghệ Tạp chí Nghiên cứu KH&CN quân sự, Số 65, 02 - 2020 191 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 *Ghi chú: TT*: kết quả tính toán trong bài báo; B*: Brischetto Nhận xét: Các kết quả tính toán của Salvatore Brischetto [11] là lời giải sử dụng lý thuyết đàn hồi 3D. Phân tích, so sánh kết quả nhận được cho thấy, mô hình toán học và chương trình tính toán đảm bảo độ tin cậy. 3. KHẢO SÁT ẢNH HƯỞNG CỦA CÁC THAM SỐ 3.1. Nghiên cứu ảnh hưởng của số lớp Trong mục này, để đánh giá ảnh hưởng của các tham số khác nhau, ta nghiên cứu vỏ trụ composite lớp, có độ dày h = 0.1, vật liệu E1= 25E2, G23= 0.2E2, G13= G12 =0.5E2, 12=0.25, chịu tải cơ phân bố đều thành trong vỏ; vỏ trụ chịu điều kiên liên kết ngàm hai đầu như hình 2. Hình 2. Mô hình vỏ trụ chịu tải phân bố đối xứng thành trong. Khảo sát ảnh hưởng của số lớp vật liệu (n) đến chuyển vị và ứng suất của các vỏ. Kết quả tính toán cho chuyển vị và ứng suất không thứ nguyên của vỏ trụ có cùng thông số kết cấu, vật liệu nhưng khác nhau về số lớp vật liệu được thể hiện trong bảng 2, theo công thức (15). Bảng 2. Chuyển vị và ứng suất không thứ nguyên của vỏ trụ composite chịu cơ phân bố đều thành trong; vỏ chịu liên kết ngàm 2 đầu với thông số (L/R=4, S=30). n (Số lớp) w , 2 2         0  , 2 2         / 2h  , 2 2         / 2h z , 2 2         / 4h z , 2          / 3h [900/0] 0.0208 -0.00345 0.00375 0.02513 0.63830 0.320384 0.589979 -0.48437 -0.53962 [900/0]2 0.0209 -0.00333 0.00366 0.02480 0.64237 0.345408 0.613430 -0.74105 -0.83077 [900/0]3 0.0210 -0.00329 0.00363 0.02463 0.64379 0.353199 0.625293 -0.64788 -0.84186 [900/0]4 0.0210 -0.00327 0.00361 0.02464 0.64452 0.359199 0.362529 -0.59192 -1.0014 [900/0]5 0.0210 -0.00326 0.00360 0.02461 0.64496 0.363554 0.628858 -0.58412 -0.91854 192 võng c vào s 3.2 trạng thái ứng suất biến dạng của vỏ E2 L/R hợp tải trọng sau: ở đây, có diễn ở h Nh - Trong trư - Đ . Ảnh h Kh =E =4; Nghiên c Trư Khi đó k N. T. Thanh ận xét: ủa vỏ trụ composite lớp c ộ v ố lư ảo sát ảnh h 3=10.3Gpa, G23=6.2GPa S ờng hợp 1: ình 3, hình 4.a và b a1) b1). õng c ợng lớp v ư =[10..50]; 0 Ứ Ứng ởng của một số dạng tải trọng c ứu chuyển vị v =L/R ết quả tính toán chuyển vị v ng su , ờng hợp phân bố lớp đối xứng th ủa vỏ trụ ảnh h , w su à cách s ưởng của các dạng tải trọng (không kể tới tác động của nhiệt độ) đến h= Vỏ chịu tải trọng phân bố theo quy luật tuyến tính theo công thức: Q10   ất ất , T. V. Hùng 0.1; =1.5 0 10 .E w Q hS   Q 1 3 tại biên tại bi ắp xếp các lớp có sự thay đổi chuyển vị v Q0=10 à tr 0; ảng 3, theo công thức sau: ; , , , , “ ùng d ưởng nhiều khi tổng số độ d , ạng thái ứng suất của vỏ trụ composite lớp trong các tr Q11       ên Tính toán v G 6), v ( )q Q Q =     ạng v 13= G ỏ hai đầu ng   Q0. 0 0 11 10 z zz à t 12 0   à ứng suất không thứ nguy ỏ trụ composite lớp ỷ lệ chiều d ơ h trụ composte lớp, vật liệu =6.9  ì s ọc GPa àm.     ố lớp ảnh h    a2). b2). Cơ k , , , , ỹ thuật & Kỹ thuật c ày. ày thay đ 12 0Q Ứng suất Ứng suất =0.22, z zz ưởng không nhiều đến độ hư ổi. Nghĩa l à     ớng tiếp cận giải tích. ứng suất. 13=0.49 ên c   ủa vỏ đ tại tại ơ khí đ (E ),     à ph 1=131GPa [90 0 ộng lực ụ thuộc 0/0/90 ược biểu / 2 0 / 2 ường (16 (17 ” , 0], ) ) Nghiên cứu khoa học công nghệ Tạp chí Nghiên cứu KH&CN quân sự, Số 65, 02 - 2020 193 c). Ứng suất zz tại 0 / 2  d). Ứng suất z tại biên 0  Hình 3. Sự biến đổi của ứng suất không thứ nguyên theo chiều dày vỏ. a). Vỏ chịu tác dụng tải trọng theo quy luật tuyến tính. b). Vỏ chịu tác dụng tải trọng theo quy luật hàm tuần hoàn. Hình 4. Chuyển vị không thứ nguyên của vỏ composite có độ dày khác nhau chịu các dạng tải trọng áp phân bố thành trong theo các quy luật tương ứng. Bảng 3. Chuyển vị và ứng suất không thứ nguyên của vỏ composite có độ dày khác nhau chịu tác dụng của tải trọng tuyến tính . , 2 2 (0) w          , 2 2 / 2h            , 2 2 / 2h            , 2 2 / 2 zz h            , 2 / 2h             , 2 / 2h             , 2 / 3 z h             , 2 / 4 zz h           Q1S10 0.011062 0.15254 12.809 0.20713 -3.5430 -1.1658 -41.238 -2.259 -0.2094 14.926 0.72486 4.02295 1.32376 -7.8935 5.6420 Q1S20 0.045186 0.34988 27.252 0.22248 -7.2532 -2.3867 -75.641 -2.8194 -0.0159 29.181 0.74300 7.87738 2.5921 -13.678 11.402 Q1S30 0.102293 0.55368 41.663 0.22805 -10.812 -3.5578 -106.89 -2.2536 0.18676 43.526 0.74883 11.584 3.8118 -18.551 17.494 Q1S40 0.182383 0.75918 56.064 0.23093 -14.288 -4.7016 -136.55 -1.1302 0.39177 57.895 0.75171 15.210 5.0050 -22.967 23.742 Q1S50 0.285456 0.96541 70.462 0.23268 -17.710 -5.8277 -165.21 0.3386 Cơ kỹ thuật & Kỹ thuật cơ khí động lực N. T. Thanh, , T. V. Hùng, “Tính toán vỏ trụ composite lớp hướng tiếp cận giải tích.” 194 0.59772 72.272 0.75342 18.784 6.1809 -27.10 30.084 Nhận xét: - 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 (khoảng >5h) cũng có dạng tuyến tính tương tự. Khi đó, 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 vỏ. Trường hợp 2: Vỏ chịu tải trọng phân bố theo quy luật tuần hoàn dạng: 21 20 0 2 ( ) sin ,q Q Q           (18) với (0=L/R, (Q20=1.5Q0; Q21= Q0), khi đó kết quả tính toán chuyển vị và ứng suất không thứ nguyên của vỏ được biểu thị ở hình 4.b và bảng 4, theo công thức (17). Bảng 4. Chuyển vị và ứng suất không thứ nguyên của vỏ composite có độ dày khác nhau chịu tác dụng tải trọng tuần hoàn. , 2 2 (0) w          , 2 2 / 2h            , 2 2 / 2h            , 2 2 / 2 zz h            , 2 / 2h             , 2 / 2h             , 2 / 3 z h             , 2 / 4 zz h           Q2S10 0.01106 0.1525 12.8089 0.2071 -7.3601 -2.4219 -84.953 -4.8816 -0.2094 14.9259 0.7249 8.3279 2.7403 -16.218 10.9063 Q2S20 0.04519 0.3499 27.2519 0.2225 -14.993 -4.9336 -155.10 -6.2156 -0.0160 29.1814 0.7430 16.098 5.2970 -27.995 22.6252 Q2S30 0.10229 0.5537 41.6629 0.2281 -22.290 -7.3344 -218.64 -5.3269 0.1868 43.5267 0.7488 23.536 7.7444 -37.905 34.9352 Q2S40 0.18238 0.7592 56.0646 0.2309 -29.403 -9.675 -278.86 -3.3377 0.3918 57.8951 0.7517 30.795 10.133 -46.886 47.522 Q2S50 0.28546 0.9654 70.4625 0.2327 -36.399 -11.977 -337.01 -0.6676 0.5977 72.273 0.7534 37.939 12.484 -55.287 60.273 Nhận xét: - 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 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ỏ nhắm giảm hiện tượng phá hủy do hiệu ứng tại các vị trí này. 4. KẾT LUẬN Trên cơ sở tính toán lý thuyết và phân tích số trình bày trong công trình này, có thể rút ra những kết luận chủ yếu sau: 1. Mô hình được trình bày trong công trình có tính tới ảnh hưởng của ứng suất, biến dạng pháp ngang (bị bỏ qua trong hầu kết các lý thuyết bậc nhất và lý thuyết biến dạng trượt bậc cao), do đó mở rộng phạm vi ứng dụng của mô hình tính toán đang nghiên cứu. Nghiên cứu khoa học công nghệ Tạp chí Nghiên cứu KH&CN quân sự, Số 65, 02 - 2020 195 2. Thực hiện khảo sát ảnh hưởng của số lớp và các một số quy luật phân bố tải trọng cơ học tác dụng lên các vỏ có độ dày khác nhau đến chuyển vị, ứng suất của vỏ trụ composite lớp với điều kiện liên kết ngàm hai đầu. 3. Hiện tượng tập trung ứng suất tại biên và các thiết diện có bước nhảy về tải tác dụng là đáng kể. Điều này có ý nghĩa quan trọng trong việc thiết kế tối ưu các kết cấu chịu tải trọng trong thực tiễn. TÀI LIỆU THAM KHẢO [1]. Timoshenko.S, Woinowsky-Krieger S (1959), “Theory of plates and shells”, McGrawHill. [2]. L. Gol'Denveizer, Th. Von Kármán, H. L. Dryden (2014), “Theory of elastic thin shells”. Pergamon Press; 1961 [3]. Flügge.W (1960), “Stresses in shells”, Springer-Verlag. [4]. Reissner E (1975), “On transverse bending of plates, including the effect of transverse shear deformation”, Int J Solids Struct, 11:569–73. [5]. J.N. Reddy (2004), “Mechanics of laminated composite plate and shell theory and analysis”, second edition, CRC Press. [6]. Aydogdu M (2009). “A new shear deformation theory for laminated composite plates”. Composite Structures, 89, pp. 94–101. [7]. Asadi.E, Wang.W, Qatu.MS (2012), “Static and vibration analyses of thick deep laminated cylindrical shells using 3D and various shear deformation theories”, Compos Struct, 94:494–500 [8]. Kant.T, Swaminathan.K (2002), “Analytical solutions for the static analysis of laminated composite and sandwich plates based on a higher order refined theory”, Compos Struct, 56, pp. 329–44. [9]. Firsanov.V.V, Doan.T.N (2015), “Investigation of the statics and free vibrations of cylindrical shells on the basis of a nonclassical theory”, Composites Mechanics Computations Applications: An International Journal, Vol. 6, Issue 2, pp. 135-166. [10]. N.T Thanh, T.N Đoàn, N.V Hà, T.V Hùng và H.M 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 Cơ học Vật rắn biến dạng toàn quốc lần thứ XIV, tr 578-585. [11]. Brischetto.S (2017), “Exact three-dimensional static analysis of single- and multi- layered plates and shells”, Composite Part B, pp.1 - 39. [12]. Vasiliev.V.V. and Lurie.S.A (1990), “To a problem of improvement of the theory of shallow shells”, Izv. RAN, Mekh-Tverd-Tela, no. 6, pp.139–146. ABSTRACT CALCULATING COMPOSITE GRADES BASED ON THE THEORY OF SLIDING HIGH QUASI-3D ACCORDING TO ANALYSIS APPROACH The paper presents the problem of calculating the stress state of the composite cylindrical casing deformation based on the Quasi-3D high-grade shear deformation theory according to the analytical approach. In the paper, the establishment of three-dimensional equations in nonlinear elastic theory into two- dimensional nonlinear equations for the cylindrical method using the method of transforming and analyzing displacement fields into multi-function chains awake according to the shell thickness. Based on the received equations, the equilibrium Cơ kỹ thuật & Kỹ thuật cơ khí động lực N. T. Thanh, , T. V. Hùng, “Tính toán vỏ trụ composite lớp hướng tiếp cận giải tích.” 196 equations and the corresponding boundary conditions for the case of displacement polynomial analysis are given. Thereby survey the stress and deformation states of composite cylindrical shell with different parameters and give the application scope for each case. Keywords: Composite layer cylinder shell; Strain stress state; Variational method; High-order shear deformation theory; Displacement; Equilibrium equation system. Nhận bài ngày 26 tháng 08 năm 2019 Hoàn thiện ngày 13 tháng 12 năm 2019 Chấp nhận đăng ngày 17 tháng 02 năm 2020 Địa chỉ: 1Viện Khoa học và Công nghệ quân sự; 2Học viện Kỹ thuật Quân sự. *Email: thanhvtv2010@gmail.com.