Nghiên cứu chống sạt lở bờ sông Đồng Tháp, dùng phương pháp phân tích trực tiếp

TR B ăGIÁOăD CăVẨă ẨOăT O NGă IăH CăC UăLONG TR NHăV NăD N NGHIểNăC UăCH NGăS TăL ăB ăSỌNGă DỐNGăPH NGăTHÁP,ă NGăPHÁPăPHỂNăTÍCHăTR CăTI P CHUYÊN NGÀNH: K ăTHU TăXỂYăD NG CỌNGăTRỊNHăDỂNăD NGăVẨăCỌNGăNGHI P Mà NGÀNH: 60 58 02 08 LU NăV NăTH CăS K ăTHU TăXỂYăD NGăCỌNGăTRỊNHă DỂNăD NGăVẨăCỌNGăNGHI P NG IăH NGăD NăKHOAăH C:ă PGS.TSăLểăV NăC NH V nhăLong,ăthángă11ăn mă2016 L IăCAMă OAN Tơiăxinăcamăđoan,ăđơyălƠăcơngătrìnhăkhoaăh cădoăchínhătơiăth căhi n.ăCácăk tăqu ăvƠă s ăli uătrongălu năv nălƠătrungăth c,ăch aăđ căcôngăb ătrongăb tăkìăcơngătrìnhăkhoaă h cănƠoăkhác.ăTơiăxinăhoƠnătoƠnăch uătráchănhi măv ătínhăxácăth căvƠătínhăngună b năc aălu năv n Tácăgi ălu năv n Tr nhăV năD n L IăC Mă N Tr c h t, xin bày t lòng bi tă nă chơnă thƠnhă đ i v iă υGS.TSă Lêă V nă C nhăđƣăt nătìnhăh ng d n, cho nhi u ch d n khoa h c có giá tr vƠăth ng xuyên đ ng viên, t oăđi u ki n thu n l i,ăgiúpăđ tơi su t q trình h c t păc ngănh ă nghiên c u hoàn thành lu n v nănƠy.ăQuaăquáătrìnhăti p xúc, h c t p nghiên c u Th y, đƣă giúpă tôiă nơngă caoă ngh l c tính sáng t o cơng vi c, ngồi ki n th c khoa h c, tơi h c Th y tính c n cù, kiên nh n gi n d Tôi xin g i l i c mă năđ n th y PGS.TS Nguy năHoƠiăS n,ăng n n t ng v b môn Ph n t h u h n, th y PGS.TS D iăđƣăgiúp cho ng H ng Th măđƣăti p l a cho s say mê v c ăh căđ t, n n móng cơng trình, th y TS Nguy năS ăHùngăđƣătruy n d y nh ng ki n th c v c ăh căđ tănơngăcaoăđ góp ph n hồn thành lu năv nănƠy Tác gi xin trân tr ngă cámă nă cácă cánă b Khoa Qu n lý khoa h c-Sauă h c HTQT, quý th y cô đƣăgi ng d y L p Cao h c xây d ngăKhóaăIăTr h c C u Long,ălƣnhăđ o Chi c c Th y l iă ngă i i ng Tháp b năđ ng nghi p, đƣă t oăđi u ki n thu n l i,ăđ ng viên, h p tác giúpăđ tơi q trình nghiên c u hoàn thành lu n v nănƠy hoƠnăthƠnhăđ c lu năv năc a tác gi đ c s đ ng viên, ng h , s chia s k p th i t giaăđìnhătrongănh ngălúcăkhóăkh nănh t, tơi xin bày t lòng bi t năvƠăchiaăs nh ngăthƠnhăcơngăcóăđ c c a b năthơnăđ năgiaăđình V nh Long, ngày 10 tháng 11 n m 2016 H c viên Tr nhăV năD n M CăL C L I C Mă σ M C L C DANH M C CÁC KÝ HI U, CH VI T T T DANH M C HÌNH NH M U GI I THI U 1.1.ă t v năđ 1.2 T ng quan v l nhăv c nghiên c u 1.3 Tính c p thi t c aăđ tài M C TIÊU, N IăDUσGăVẨăυH σGăυHÁυăσGHIểσăC U 11 2.1 M c tiêu c aăđ tài 11 2.2 N i dung nghiên c u 11 2.3ăυh CH ngăphápălu năvƠăph ngăphápănghiênăc u 11 σGă1.ăT NG QUAN LÝ THUY T 13 1.1 D τăLụăT NG VÀ TIÊU CHU N PHÁ H YăCHτă T 13 1.1.1 Gi i h năđƠnăh i hàm ch y 13 1.1.2 Lu t ch y d o k t h p 14 1.1.3 Hàm ch y d o Morh-Coulomb 15 1.2 LÝ THUY T PHÂN TÍCH GI I H N 17 1.2.1ă nh lý c năd 1.2.2ă nh lý c n 20 i 19 1.3.ăυH σGăυHÁυăS VẨăCÁCHăTệσHăσ σGăL NG TIÊU TÁN D O CHO PH N T 22 1.3.1ăυh ngăphápăph n t h u h n: 22 1.3.2ăσ ngăl ng tiêu tán d o ph n t s d ngăph ngăphápăph n t h u h n 24 1.3.3.ăυh ngăphápăph n t h u h nătr năc nh ( ES-FEM) 25 1.3.4ăσ ngăl ng tiêu tán d o ph n t s d ngăph ngăphápăph n t h u h n tr năd a c nh 26 1.3.5 M t s nh n xét v 2ăph ngăphápăs FEM-T3 ES-FEM-T3 27 CH σGă2.ăCH σGăăTRỊσHăσịσăVẨăCÁCHăTHI T L P BÀI TỐN T Iă Uă HĨA T L I GI I C N TRÊN 28 2.1ă σHăσGH A 28 2.2 CÁC D NG HÌNH NĨN 29 2.3 S D σGăCH σGăTRỊσHăNĨN CHO BÀI TỐN BI N D NG PH NG 29 2.4 THI T L P BÀI TOÁN T Iă U 30 2.4.1 Thi t l p toán t iă uăkhiădùngăph ngăphápăph n t h u h n s d ng ph n t tam giác ba nút: 30 2.4.2ăυh CH ngăphápăph n t h u h n d a c nh (ES-FEM-T3) 31 σGă3.ăVệăD TÍNH TỐN 33 3.1ă Că I Mă 3.1.1 L păđ t ph : A CH T CƠNG TRÌNH 33 t ph cát pha, s n s i, g ch đá, l n r th c v t 34 3.1.2 L păđ t s 1: Sét màu xám nâu, nâu đ , tr ng thái d o ch y đ n d o c ng 34 3.1.3.Th u kính 1: Cát l n b i màu xám đen, tr ng thái r t r i 35 3.1.4 L păđ t s 2: B i l n sét màu xám nâu đen, tr ng thái ch y đ n d o ch y 35 3.1.5 Th u kính 2: Sét màu xám xanh, tr ng thái d o c ng đ n n a c ng 36 3.1.6 L păđ t s 3: Cát l n b i màu xám nâu xanh, tr ng thái r i đ n ch t v a 36 3.1.7 L păđ t s 4: Cát pha b i màu xám h ng, tr ng thái d o 37 3.2.ăMÔăHỊσHăTệσHăTτÁσăVẨăL I PH N T H U H N: 37 3.3 K T QU TÍNH TỐN: 38 CH σGă4.ăυHỂσăTệCHă σă NH MÁI D C S D NG ES-FEM VÀ CH σGăTRỊσHăHỊσHăσịσ 44 4.1 PHÂN TÍCH S 4.1.1ăS ăl 4.1.2ă σă căph NH MÁI D C 44 ngăphápătínhătốn 44 t v năđ 44 4.1.3 Bài toán t iă u 45 4.1.4 Mơ hình tính tốn 46 4.2ă σă σHăKHÔσGăTHτÁTăσ CăKHIăMịσGăCÔσGăTRỊσHă T LÊN MÁI D C 46 4.2.1ăS ăl CH căph ngăphápătínhătốn 46 σGă5.ăK T LU N VÀ KI N NGH 53 CODE THAMăKH O 56 DANHăM CăCÁCăKụăHI U,ăCH ăVI TăT T ăă Góc n i ma sát c L c dính ES-FEM υh ngăphápăph n t h u h nătr năd a c nh (edgebased smoothed finite element method) FEM Ph n t h u h n (finite element method) PTHH Ph n t h u h n SCOCP Ch x, y, z Các ng su t pháp xy, xz, yz Các ng su t ti p  Góc giãn n ngătrìnhănónăb c hai Bi n d ngăt cu ngăđ i Gócămaăsátătrongăđi u ki n có c k t-khơngăthốtăn k Mi nătr n B B r ng móng L Kho ng cách t mépăTaluyăđ n móng Ns H s năđ nh mái d c c DANHăM CăHỊNHă NH Hình 1.1 Quan h gi a ng su t bi n d ng c a v t li u ng x đƠn-d oălỦăt ng 14 Hình S minh h a hình h c c a lu t ch y d o k t h p 15 Hình ng x th t c aăđ t ng x đƠnăd oălỦăt Hình 1.4 Mơ hình Morh s c ch ng c tăthoátăn ng .16 c c aăđ t 16 Hình 1.5 υh ngăc aăvect ăgiaăs bi n d ng d o h tr c    choăhaiătr ng h p:ăa)ăđ tăkhơngăthốtăn căvƠăb)ăđ tăthốtăn c .17 Hình 1.6 Nghi m c a l i gi i c n c năd i cho tốn phân tích gi i h n 18 Hìnhă1.7ăS ăđ phân tích gi i h n .18 Hìnhă1.ă8ăă i u ki n biên l c chuy n v 19 Hình 1.9 T aăđ ph n t tam giác 22 Hình 1.10 Mi nătr nă  k d a c nh 25 Hình 2.1 Khơng gian hình nón ầầầầầầầầầầầầầầầầầầầầ.28 Hình 2.2 Di n tích mi nătr năk d a c nh 31 Hình 3.1 M t c tăngangăđ a ch t d c b sơng S Th ngầầầầầầầầầầ.33 Hình M t c tăngangăđ a ch t b sông Ti n 34 Hình 3 M t c tăđi n hình mái d c b sông MC3 37 Hình L i ph n t h u h n: MC1 38 Hình C ăc uătr t c a mái d c theo m t c t MC3 .39 Hình T căđ bi n d ng c a mái d c theo m t c t MC3 .39 Hình σ ngăl Hình C ăc uătr ng tiêu tán d o c a mái d c theo m t c t MC3 39 t c a mái d c theo m t c t MC3 theo thông s l păđ t .40 Hình T căđ bi n d ng c a mái d c theo m t c t MC3 theo thông s l păđ t 40 Hình 3.10 σ ngăl ng tiêu tán d o c a mái d c c a MC3 theo thơng s l păđ t 41 Hình 11 Góc mái d c 500: h s năđ nh c a mái d c Ns = 1.5685 41 Hình 12 Góc mái d c 700: h s năđ nh c a mái d c Ns = 1.2082 42 Hình 13 Góc mái d c 800: h s năđ nh c a mái d c Ns = 1.0572 42 Hình 14 Góc mái d c 900: h s năđ nh c a mái d c Ns = 0.9127 42 Hình s ăđ mái d cầầầầầầầầầầầầầầầầầầầầầầầ.45 Hình Phân tích năđ nh mái d c:ăs ăđ hình h căvƠăcáchăchiaăl i 46 Hình Mơ hình hình h căbƠiătốnămóng,ăcơngătrìnhăđ t mái d c .48 Hình 4 năđ nh mái d c: L/B = 0, Ns = 11.7234 .48 Hình năđ nh mái d c: L/B = 1, Ns = 18.5593 .49 Hình năđ nh mái d c: L/B = 2, Ns = 22.3178 .49 Hình năđ nh mái d c: L/B = 3, Ns = 26.8507 .50 Hình năđ nh mái d c: L/B = 4, Ns = 31.4405 .50 Hình năđ nh mái d c: L/B = 5, Ns = 35.7359 .51 Hình 10 năđ nh mái d c: L/B = 6, Ns = 39.5911 51 Hình 11 năđ nh mái d c: L/B = 7, Ns = 41.8449 52 M ă U 1.ăGI IăTHI U 1.1.ă tăv năđ Hi n nay, tình hình s t l trênă đ a bàn t nhă ng b ng Sông C u Long nói chung t nhă ngăThápănóiăriêngăđangădi n h t s c ph c t p làm t n hao đ n tính m ng tài s n c aănhơnădơn.ă ng Tháp t nhăđ u ngu n sông C u Long c a ch y qua Vi t Nam, có Sơng Ti n Sơng H u ch y qua có nhi u thu n l i 02ăconăSôngănƠyăđemăl iănh : + Là n giao thông huy t m ch n i li n t nh khu v c g n k t v iăn c b n Campuchia, +ă σ iă cungă c pă n c ng t, phù sa, tôm cá, l n nh t cho dân sinh, nông nghi p, th y l i, công nghi p, lâm nghi p, xây d ng cho tồn T nh + Là năthốtăl ăchínhăc a tồn b Khu v că ng b ng sông C u Long, + K t h p v i h th ng sơng ngòi, kênh r ch góp ph n b o t n, năđ nh, cân b ngămôiătr ngăsinhătháiăchoă ng Tháp t nh lân c n + Là ngu n cung c p v t li u xây d ng (cát san l p cát xây d ng cơng trình) +ăCóăỦăngh aăquanătr ngăđ i v i an ninh qu c phòng c a T nh nói riêng c a c vùng Nam B nói chung Tuy nhiên, bên c nhăđóăti m năc ngăăkhơngăítănh ng hi m nguy s t l b sơng gây nên T n mă2005ăậ 2014ăcóă84ăđo n b sơng thu căđ a ph n c a 42 xã b s t l , t ng chi u dài s t l 703 Km, di n tích s t l 283 ha, t ng s h dân c n ph i di d i 18.854 h ,ăđƣădiăd i 6.449 h l i 12.405 h ph i ti p t c di d i, t ng thi t h i 227,18 t đ ng căđi măđ a ch t cơng trình d c theo n sơng Ti n, sông H u nguyên nhân ch y u gây hi năt ng m t năđ nh c a b sông Ti n, sông H u (trên góc đ đ a ch t) là:  ng b ng sông C uăLongăăđ c c u t o b iăcácăđ t tr m tích phù sa tr có ngu n g c bi n sông - bi n h n h p v i b dày l n  Các l păđ t t o nên n năđ tăđ uălƠăđ t y u v i tr ng thái v t lý kém, tính ch t c ăh c kém, mang nh ngăđ cătr ngăc a nh ng lo iăđ tăđ c bi t V i góc ma sát = 4029', C= 0,071 kG/cm2(đ sâu: 2-6m); = 6028', C= 0,086 kG/cm2 (đ sâu: 7-14m); = 22041', C= 0,06 kG/cm2 (đ sâu: 16-30m)  σ c ng m m t nh ng y u t quan tr ngăđiătheoăgơyăraăcácăhi n t ng xói ng m, xúc bi năvƠăđ c bi t cát ch y Hi nătr ngăv ăs tăl ăt iăT ngăkhoăx ngăd uăTr năQu căTo n ậ ngăTháp 1.2ăT ngăquanăv ăl nhăv cănghiênăc u σg i đ u tiên nghiên c u v s c ch u t i c a móng nơng n nă đ ng nh t Prandlt [12], nghiên c u Prandlt đƣ đ a c c u m tă tr t xác c a n năđ t đ ng nh t t xác đ nh đ c giá tr t i tr ng gi i h n Sau đó, Terzaghi [16] c ng đƣ đ a công th c xác đ nh s c ch u t i gi i h n d a mơ hình m t tr t tam giác Tuy nhiên, công th c g p nhi u h n ch áp d ng cho n n đ t khác,ăđ c bi tăđ i v i n n đ t th c t v i nhi u l p đ t khác Vì v y, nhi u nghiên c u phân tích gi i h n d a đ nh lý c n c năd i k t h p v i ph ng pháp ph n t h u h n đƣ đ c phát tri n.ă u tiên nghiên c u c a Lysmer [9], tác gi đƣ dùng ph ng pháp PTHH t i u n tính đ xác đ nh c năd i c a t i tr ng gi i h n.ăTr ng ng su t đ c x p x b i PTHH b t liên t c, ngh a m t t a đ măđ c gán nút c a ph n t khác nhau, khác v i PTHH truy n th ng t a đ t ng ng v i ch nút Sau đó, nhi u nghiên c u khác đƣ đ c phát tri n,ăđ c bi t nhóm nghiên c u c a GS Sloan [5-8;14,15] Trong [14],ătr ng chuy n v b t liên t c đ c s đ tính tốn c n c a t i gi i h n, cách r i r c ph n t gi ng nh [9] Trong nghiên c u này, tiêu chu n Morh-Coloumbăđ c n tính hóa thành nhi u đo n th ng thu t tốn t i u n tính đ c s d ng đ thu đ cătr ng chuy n v b t liên t c.ă gi m kh i l ng tính tốn t ngă đ xác, d ng phi n c a tiêu chu n Morh-Coloumbăđ c dùng [1,6,5]ăđ thuăđ c c năd i c n c a s c ch u t i c a n n đ t Trong phân tích gi i h n, tr ng ng su t chuy n v tr ng thái gi i h n có tính c c b - ng su t bi n d ng d o ch t p trung t i m t Cao độ bờ H Cao độ đáy sông D Cao độ Hình 4.1 s ăđ ămáiăd c 4.1.3ăBƠiătoánăt iă u Bài toán t iă uăchoăb i l i gi i c n    min( D(&)  Wext0 ) Ràng bu c: (4.1) Wext   u  u0 tr ên biên Thành ph năgơyăraătr t mái d c tr ngăl ng b n thân c aăđ t Do khơng có thành ph n l c m t tác d ng nên: Wext0  Bài tốn t iă uăhóaătr thành:    D(&) Ràng bu c: (4.2) Wext   u  u0 tr ên biên Bài toán t iă uăhóaăkhiădùngăph ngăphápăăES-FEM T3 tr thành ned       Ai cosi ti i 1 Ràng bu c:  u  u0 tr ên biên  Wext   2 ti   xx   yy    xy2 i  1, 2, , ned  ned: t ng s c nh mơ hình phân tích gi i h n (4.3) 4.1.4 Mơ hình tính tốn Trongă tr đ ng h pă đ nă gi n nh t,ă đ tă đ că xemă lƠă đ ng nh t thông s đ t c xem xét tốn phân tích gi i h n    ; c  Hình υhơnătíchă năđ nhămáiăd c:ăs ăđ ăhìnhăh căvƠăcáchăchiaăl Trên Hình 4.2 th hi năl 3ănútăvƠă1711ănútăđ i ph n t h u h nădùngăđ tính tốn v i 3311 ph n t tam giác c t oătheoăkíchăth c ph n t khác dùng ph n m m GMSH 4.2 ă Nă NHăKHỌNGăTHOÁTăN LểNăMÁIăD C 4.2.1 S ăl căph i CăKHIăMịNGăCỌNGăTRỊNHă Tă ng pháp tính tốn υhơnătíchăđ năđ nh c a mái d c có v t th ,ămóngăcơngătrìnhăđ t lên mái d c r t phúc t p, vìăc ăch tr t s k t h p c a y u t g m t i lên móng c tr ngăl ng c aăđ t Nhi uăph ngăphápăđ căđ xu tăđ gi i quy t v năđ bao g m:ăđ ngătr t th ng c a Sokolovski, cân b ng gi i h n c a Meyerhof , lý thuy t m t ch y d o c aă Buhană vƠă Garnier,ầ Tuy nhiên, hai th p niên qua S W Sloan vƠă cácă đ ng nghi pă đƣă ápă d ng lý thuy t phân tích gi i h n gi i quy t r t nhi u toánătrongăđ a k thu t s d ng ph n t h u h n ph n t b t liên t c Trong ph n lý thuy t phân tích gi i h n s d ng l i gi i c n dùng ESFEMT3ăvƠăch ngătrìnhănónăs đ c ti n hành kh o sát cho v năđ nƠyăcóăxétăđ n s nhăh ng c a v tríămóngăđ t mái d c Th c tr ng hi n b sơng có nhi u cơng trình cơng c ng nhà đ c xây d ng, v i nhi u k t c u ph c t p, v tríăđ t khơng theo m t tr t t nh t đ nh Vì v yăđ đ m b o vi căđánhăgiáăanătoƠnăchoăcácăcơngătrìnhănƠy,ăbài tốn phân tích năđ nh mái d căcóăxétăđ n t i tr ng c a móng cơng trình c năđ s đ c quan tâm c trình bày ph n V năđ ch u t i c a móng c ngăđ t g n mái d c cho b i Hình 4.3 Móng có b r ngăBăđ t n năđ ng nh t v i góc mái d c  v i chi u cao H kho ng cách t mépămóngăđ năđ nh mái d c L chu năđƠnăậ d oălỦăt t n n gi đ nhălƠăđ ng nh t ng x theo tiêu ng Morh-Coulomb T i gi i h n c a móng c ng khơng nh ng b nhăh ng b i góc mái d c  kho ngăcáchăđ t móng L so v iăđ nh mái d c Kh n ngăch u t i gi i h năc ngăph thu c vào tr ngăl ng riêng  c aăđ t, nhăh ngăđ n năđ nh t ng th c a mái d c i u khơng gi ng v iă móngă đ t m tăđ t ph ng, mà kh n ngăch u t i c c h năđ c l p v i  Kh n ngăch u t i c c h n c aăbƠiătoánăđ th đ c phát bi uăd c xem xét đơyăcóă i d ng:  L c p q H ,   f , , u , B  B B B B Trongăđó: - p:ălƠăápăl căgi iăh nătrungăbìnhătácăđ ngălênămóngăvƠă - q:ălƠăt iăphơnăb ă Theoăđóăkh n ngăch u t iăđ cu/B,ătrongăđóă nhăh c trình bày theo d ng không th nguyên , L/B, ng c aăđ thơ c a móng t i phân b q/B s đ tích riêng l H/B = t t c phân tích, đ đ s phá h y x y mái d c c phân chân L B Móng c ng H q Cu, , u =0 Hình 4.3 Mơăhìnhăhìnhăh căbƠiătốnămóng, cơng trình đ tătrênămáiăd c Trong th c t , m t c t b sơng b xóiămònăvƠăthayăđ i theo th iăgian,ănênăđ an toƠnăchúngătaăxemăxétătr ng h p b t l i nh t mái d c th ngăđ ng (b qua ph n mái d c,ăngh aălƠăt ngăđ an toàn cho tính tốn) Các thơng s đ a ch tăc ngăs l y theoătr ng h p b t l i nh t L p đ a ch t s 2, đ c thi t l p tốn phân tích gi i h n cu/B = 5, q/B = 0, kho ngăcáchăđ t móng so v iăđ nh mái d c L/B = ÷ K t qu tính tốn v c ăc uătr hi n Hình 4.4 đ n Hình 4.11 t h s năđ nh mái d c đ Hình 4.4 năđ nhămáiăd c:ăL/B = 0, Ns = 11.7234 c th Hình năđ nhămáiăd c:ăL/B = 1, Ns = 18.5593 Hình năđ nhămáiăd c:ăL/B = 2, Ns = 22.3178 Hình năđ nhămáiăd c:ăL/B = 3, Ns = 26.8507 Hình năđ nhămáiăd c:ăL/B = 4, Ns = 31.4405 Hình năđ nhămáiăd c:ăL/B = 5, Ns = 35.7359 Hình 10 năđ nhămáiăd c:ăL/B = 6, Ns = 39.5911 Hình 11 năđ nhămáiăd c:ăL/B = 7, Ns = 41.8449 Nh n xét: T k t qu ta th y: - KhiămóngăcơngătrìnhăcƠngăxaăđ nh mái d c (L/B l n) h s năđ nh mái d căcƠngăt ng,ăngh aălƠăvi c xây d ng cơng trình xa mép b sơng an tồn - Khi L/B = 7, móng cơng trình khơng nhăh ngăđ n mái d c n a, mà ch tácăđ ng c c b nh ăbƠiătốnămóngănơngătrênăn n thiên nhiên CH NGă5.ăK TăLU NăVẨăKI NăNGH 5.1 K T LU N Trong lu năv nănƠy,ăthu t tốn phân tích gi i h n d aătrênăph ngăphápăph n t h u h nă tr nă ES-FEM k thu t t iă uă nónă b că haiă đ că đ xu t b i Le [3] đ c phát tri n cho mơ hình Morh-Coloumb m r ng áp d ng cho toán n đ nh mái d căt ngă ng v i s li uăđ a ch t c a b sông Ti năđiăngangăquaăđ a ph n t nhă ng Tháp H s năđ nh mái d căvƠăc ăc uătr tăt ngă ng v i mái d c th c c a b sông Ti năđ căxácăđ nh khuy n cáo cho công tác thi t k ch ng s t l trênăđ a bàn t nhă ng Tháp υh ngăphápăđƣăđ c áp d ngătínhătốnăchoăcácătr ng h p b t l i nh t v đ a ch t góc mái d c.ăTrênăc ăs đóăđ aăraănh ng giá tr tham chi uăchoăcácăđ năv liên quan công tác th măđ nh ki m tra, giám sát cơng trình b kè cơng trình khác v ch ng s t l Trong lu năv n,ăs nhăh ng c a móng cơng trình xây d ng b sơng c ngăđ c phân tích tính tốn K t qu tính tốn cho th y r ng, cơng trình xây d ng cƠngă xaă đ nh mái d c h s nă đ nh mái d că cƠngă t ng,ă đ ng th i kho ng cách t đ nh mái d că đ n mép móng cơng trình l nă h nă 7ă l nă kíchă th c móng (L>7B) móng cơng trình khơng gây nhăh ngăđ n mái d c mà ch tácăđ ng c c b nh ătr ng h p móng nơng n n thiên nhiên 5.2 KI N NGH Tuyănhiên,ăđ t ngăhi u qu tính tốn ng d ng th c ti n m t s v năđ sauăđơyăc năđ c quan tâm: đ tăđ c giá giá tr xác c a h s năđ nh mái d c s ph n t c a l i tính tốn c năt ngălên,ătuyănhiênăs d năđ n chi phí tính tốn s t ngătheo.ă kh c ph c v năđ k thu t tái t oăl i t đ ng (mesh adaptivity) c năđ c xemăxét,ăngh aălƠăl i t i vùng g năđ ngătr t c năđ c t đ ng chia nh đ gi m sai s tính tốn - C n xem xét tính tốn cho n năkhơngăđ ng nh t - Trongăph ngăphápăph n t h u h nătr n,ăbi n d ng c nh c a hai ph n t lân c năđ călƠmătr nădoăđóăd năđ năkhóăkh nătrongăvi c khai báo thông s đ a ch t c a ph n t n m biên ti p giáp c a l păđ t Trong nghiên c u ti p c n x lý v năđ đ có th mơăhìnhăđúngăbƠiătốnăcóănhi u l păđ tăcóăđ cătr ngăkhácă Vi c nghiên c u s t l b sông v năđ r t ph c t p, c n k t h p y u t nh :ă đ a ch t, th yăv n, áp l c ch đ ng, áp l c b đ ng,ăl uăl ng dòng ch y, áp l c dòng ch y, hi năt ng bán nh t tri u,ầđ c bi t tình hình khai thác cát sông hi n nay./ TẨIăLI UăTHAMăKH O [1] A V Lyamin and S W Sloan Upper bound limit analysis using linear finite elements and non-linear programming International Journal for Numerical and Analytical Methods in Geomechanics, 26:181-216, 2002 [2] C V Le, H Askes, and M Gilbert A locking-free stabilized kinematic EFG model for plane strain limit analysis Computers and Structures, 106-107:1ậ8, 2012 [3] C.V Le, Nguyen-Xuan H, Askes H, Rabczuk T, Nguyen-Thoi T Computation of limit load using edge-based smoothed finite element method and secondorder cone programming International Journal of Computational Methods 2013; 10(1):1340004 [4] C.V Le, H Nguyen-Xuan, H Askes, S Bordas, T Rabczuk, and H NguyenVinh A cell-based smoothed finite element method for kinematic limit analysis International Journal for Numerical Methods in Engineering, 83:1651-1674, 2010 [5] Lyamin, A V & Sloan, S W 2002 Upper bound limit analysis using linear finite elements and non-linear programming International Journal for Numerical and Analytical Methods in Geomechanics 26, 181ậ216 [6] Lyamin, A V & Sloan, S W 2002 Lower bound limit analysis using nonlinear programming International Journal for Numerical Methods in Engineering 55, 573ậ611 [7] Lyamin, A V & Sloan, S.W 2003 Mesh generation for lower bound limit analysis Advances in Engineering Software 34, 321ậ338 [8] Lyamin, A V., Sloan, S W., Krabbenhoft, K & Hjiaj, M 2005 Lower bound limit analysis with adaptive remeshing International Journal for Numerical Methods in Engineering 63, 1961ậ1974 [9] Lysmer, J 1970 Limit analysis of plane problems in soil mechanics Journal of the Soil Mechanics and Foundations Division, ASCE 96, 1311ậ1334 [10] Makrodimopoulos A, Martin CM Lower bound limit analysis of cohesivefrictional materials using second-order cone programming Int J Numer Methods Eng 2006; 66:604ậ34 [11] Makrodimopoulos A, Martin CM Upper bound limit analysis using simplex strain elements and second-order cone programming International Journal for Numerical and Analytical Methods in Geomechanics 2006; 31:835ậ865 [12] Prandtl L Ueber die haerte plastischer koerper Nachrichtex der Akademie der Wissenschaften in Gottingen II Mathematisch-Physikalische Klasse II 1920; 12:74ậ85 [13] Smith C, Gilbert M 2007 Application of discontinuity layout optimization to plane plasticity problems Proc R Soc A 463, 2461ậ2484 [14] S W Sloan and P W Kleeman Upper bound limit analysis using discontinuous velocity fields Computer Methods in Applied Mechanics and Engineering, 127:293-314, 1995 [15] S W Sloan, Randolph MF Numerical prediction of collapse loads using finite element methods Int J Numer Anal Methods Geomech 1982;6:47ậ76 [16] Terzaghi, K., Theoretical Soil Mechanics, John Wiley and Sons, New York (1943) CODE THAMăKH O % -% Slope stability analysis using egde-based finite elements % Developed by Canh V Le, International Univeristy - VNU HCMC % clear all, format short global node element; global edges area_cell enrich_node ng; % cohesion force c = 0.088; phi = 6.9176*pi/180; % friction angle % weight of soil rho = 1.686; q = 0; % surcharge load % read mesh file -fid1 = fopen('mc60node.m','r'); count = 0; while tline = fgetl(fid1); if isnumeric(tline) break else [C2]= sscanf(tline,'%f %f %f %f') ; if size(C2,1)>0 count=count+1; node(count,1:2)=[C2(2) C2(3)]; end end end fclose(fid1); fid2=fopen('mc60element.m','r'); count=0; while tline = fgetl(fid2); if isnumeric(tline)%~ischar(tline) break else %[C2]= sscanf(tline,'%f %f %f %f ') ; [C2]= sscanf(tline,'%f %f %f %f %f %f %f %f') ; if size(C2,1)>0 count=count+1; %element(count,1:3)=[C2(2) C2(3) C2(4)]; element(count,1:3)=[C2(6) C2(7) C2(8)]; end end end fclose(fid2); % end read file ndof = 2; % number of displacement dofs per node % total sampling node number nnode = size(node,1); nel = size(element,1); % number of element sdof = nnode*ndof; % total system dofs for mechanical displacement Index = cell(nel,1); for i = 1:nel Index{i} = element(i,:); end % mesh plot patch('faces',element,'vertices',node,'facecolor','none'); axis equal, hold on % find nodes on boundary tol = 1e-9; x_max = max(node(:,1)); bc_left = find(abs(node(:,1))