(Luận văn thạc sĩ) xây dựng mô hình kết hợp phần tử rời rạc và phần tử chất lỏng nhằm mô phỏng ứng xử của đất dính trong thí nghiệm 3 trục

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN VĂN BÌNH XÂY DỰNG MƠ HÌNH KẾT HỢP PHẦN TỬ RỜI RẠC VÀ PHẦN TỬ CHẤT LỎNG NHẰM MƠ PHỎNG ỨNG XỬ CỦA ÐẤT DÍNH TRONG THÍ NGHIỆM TRỤC NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CÔNG NGHIỆP - 60580208 S K C0 Tp Hồ Chí Minh, tháng 4/2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN VĂN BÌNH XÂY DỰNG MƠ HÌNH KẾT HỢP PHẦN TỬ RỜI RẠC VÀ PHẦN TỬ CHẤT LỎNG NHẰM MƠ PHỎNG ỨNG XỬ CỦA ĐẤT DÍNH TRONG THÍ NGHIỆM TRỤC NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP - 60580208 Tp Hồ Chí Minh, tháng 04/2016 GVHD: TS TRẦN VĂN TIẾNG CHƯƠNG 1: TỔNG QUAN 1.1 Tính cấp thiết đề tài Ngày nay, phát triển xã hội dân số, nhu cầu nhà văn phòng ngày tăng cao, cao ốc với hệ tầng hầm sâu không ngừng mọc lên Việc thi công hệ tầng hầm cơng trình thường xãy cố trượt lở thành hố đào, ổn định thành hố, bục đáy hố đào gây ảnh hưởng đến cơng trình xây dựng cơng trình lân cận Thiệt hại cố gây không nhỏ người vật chất Bên cạnh đó, cố sạt lở đất Việt Nam vào mùa mưa xãy với tần suất ngày nhiều Dọc theo đường Trường Sơn liên tiếp xảy vụ sạt lở vào mùa mưa làm gián đoạn giao thông Gần vụ sạt lở Quốc lộ (Hồ Bình) làm người chết, sạt lở đường lên núi Cấm (An Giang) làm người chết tắc đường nhiều ngày, Lào Cai làm 15 người chết Chiều ngày 16/10/2013 xã Tiên Cảnh, huyện Tiên Phước, tỉnh Quảng Nam xãy vụ sạt lở làm cha thiệt mạng, chiều tối ngày 22/11, vụ sạt lở núi nghiêm trọng xảy địa bàn huyện Sơn Tây (Quảng Ngãi) vùi lấp quan hành … Hình Sạt lở Sơn Tây (Quảng Ngãi) http://www.baomoi.com/sat-lo-nui-dat-da-vui-lap-3-co-quan-hanh-chinh-o-quangngai/c/12492271.epi HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG Sạt lở đất Lạng Sơn làm công nhân thiệt mạng xảy ngày 17/9/2014 Hình Sạt lở đất Lạng Sơn http://vnexpress.net/tin-tuc/thoi-su/sat-lo-dat-o-lang-son-lam-6-nguoi-tu-vong3080660.html Các cố thi cơng móng cơng trình hay các cố sạt lở đất vào mùa mưa ảnh hưởng nghiêm trọng đến tài sản tính mạng người cơng tác quản lý xây dựng, quản lý rủi ro thiên tai Do đó, đề tài đề xuất xây dựng mơ hình ứng xủa đất dính bão hịa nước phương pháp phần tử rời rạc (DEM) phần tử chất lỏng có khả ứng dụng để tính tốn mơ nhằm dự đốn cố xãy từ đưa biện pháp phịng ngừa hợp lý, có biện pháp ứng phó kịp thời 1.2.Tình hình ngun cứu nước Nghiên cứu thực nghiệm ứng xử vật liệu đất, đá tác dụng nhiều loại tải trọng khác thực nhiều giới (Sulem et al 2005, Tien et al 1988,…) Tuy nhiên đất loại vật liệu đa dạng, phức tạp, tính chất thay đổi theo điều kiện địa lý Hơn nữa, vấn đề xảy với đất đa dạng phá hủy trượt lở, phá hủy đẩy nổi, hóa lỏng, vấn đề sức HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG chịu tải, vấn đề lún, vấn đề tương tác đất cơng trình … Do nghiên cứu thực nghiệm ứng xử vật liệu đất đá q trình khơng ngừng Về mặt tính tốn mơ số, tác giả giới chủ yếu thực mô ứng xử đất sử dụng phương pháp phần tử hữu hạn Phương pháp khơng thích hợp để mơ phá hủy đất Với phát triển công nghệ máy tính, phương pháp phần tử rời rạc đời ngày phát triển, ứng dụng rộng rãi tính tốn mơ số ứng xử đất đá (M Hakuno 1995, R.R.O Bonilla 2004, R Shafipour et al 2008, Lui et al 2003, Zhang et al 2006, Huang et al 2008, Alonso et al 2008, Mustafa Alsaleh et al 2010,…) Tuy nhiên nghiên cứu mơ với mơ hình 2D tập trung vào số toán ứng xử đất, hay mô với mô hình 3D áp dụng cho loại đất có độ bão hịa nước 30% (L Scholtes 2008) Ở Việt Nam, phương pháp phần tử rời rạc phương pháp mẽ, vài tác giả tiếp cận với phương pháp (Nguyễn Tiến Cường et al 2012, Trần et al 2016) Các tác giả nghiên cứu chủ yếu cho bê tông sử dụng luật ứng xử có sẵng để mơ tốn thực tế Mơ ứng xử khơng nước đất dính bão hịa nước sử mơ hình rời rạc 3D mơ hình chất lỏng vấn đề mà luận văn đặt để giải HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG 1.3 Nhiệm vụ giới hạn đề tài - Xây dựng mơ hình ứng xử đất dính phương pháp phần tử rời rạc - Kết hợp với mơ hình chất lỏng cho nước lổ rỗng có sẵng - Số hóa mơ hình rời rạc ngơn ngữ C++, đưa mơ hình vào chương trình tính mã nguồn mở Yade - Mơ số - Đánh giá kết 1.4 Phương pháp nghiên cứu - Mơ số ứng xử đất đính thí nghiệm nén trục phần tử chất lỏng từ so sánh với kết thực nghiệm nhằm đánh giá tính xác mơ hình số HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG CHƯƠNG 2: CƠ SỞ LÝ THUYẾT 2.1 Phương pháp phần tử rời rạc Phương pháp phần tử rời rạc (DEM) nghiên cứu phát triển ứng dụng rộng rãi Trong đó, có nhiều dạng phương pháp phần tử rời rạc ứng dụng mô ứng xử vật liệu đất, bê tông… Trong nghiên cứu này, phương pháp phần tử rời rạc cổ điển đề xuất Cundall & Strack (1979) sử dụng Với DEM, vật liệu rời rạc hóa phần tử hình học dạng đĩa, dạng hình cầu, dạng đa diện, …Các phần tử độc lập với biến dạng tuyệt đối cứng Trong nghiên cứu này, mẫu vật liệu mô tập hợp phần tử rời rạc dạng hình cầu khơng biến dạng chịu tác dụng tải trọng, thuật toán phương pháp DEM thực theo bước sau : Xác định tương tác phần tử rời rạc dựa vào thơng số bán kính tương tác γ Sử dụng “luật tương tác” (luật ứng xử cục bộ) để tính lực tương tác hai phần tử rời rạc Lực gồm hai thành phần, lực pháp tuyến lực tiếp tuyến Tính tổng lực tương tác tác dụng lên phần tử rời rạc Định luật Newton áp dụng để tính gia tốc cho phần tử rời rạc, gia tốc tích phân theo thời gian để tìm vị trí cho phần tử Q trình tuần hồn mô kết thúc HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG Với biến dạng nhỏ, vật liệu địa kỹ thuật có tính ma sát – dính thể ứng xử đàn hồi tuyến tính Để mơ ứng xử này, lực tương tác phần tử rời rạc cần lực đàn hồi tuyến tính giảm thời gian tính tốn mơ Do nghiên cứu này, tác giả đề xuất mơ hình ứng xử tuyến tính cho đất ma sát – dính 2.2 Mơ hình ứng xử đất dính Trong mơ hình đề xuất, vector lực F, đại diện cho tác động phần tử rời rạc (DE) a lên b, không kể đến phần tử tiếp xúc, mà kể đến khoảng cách chúng nhỏ so với bán kính tương tác, điều chỉnh hệ số  , hai DE gọi tương tác khi:  (R a  R b )  Dab (1) Trong Dab khoảng cách từ trọng tâm DE a DE b, Ra Rb bán kính DE a DE b tương ứng   Đây khác biệt quan trọng so với phương pháp rời rạc cổ điển xem xét tương tác có tiếp xúc thật (  = 1) Trong mô hình này, véc tơ lực tương tác F phân tích thành véc tơ lực pháp tuyến Fn véc tơ lực tiếp tuyến Fs Hai lực tính từ chuyển vị pháp tuyến tương đối gia số chuyển vị tiếp tuyến tương ứng, thông qua hệ số độ cứng Kn (Tran et al 2011) theo phương pháp tuyến Ks (Tran et al 2011) theo phương tiếp tuyến (Hart et al 1988) Ra Rb  K n  E Ra  Rb    K  K n  s HVTH: NGUYỄN VĂN BÌNH (2) Trang GVHD: TS TRẦN VĂN TIẾNG Trong E mơ đun biến dạng thực đất xem thơng số đầu vào mơ hình, α hệ số không thứ nguyên xác định dựa theo nghiên cứu Donzé et al (Donzé 1995) Hình Lực tương tác phần tử rời rạc Lực tương tác pháp tuyến tính thông qua cập nhật luật ứng xử cho tương tác phần tử (hình 2) chia thành hai giai đoạn, chịu nén chịu kéo Trong giai đoạn chịu nén, ứng xử đất tuyến tính lực tương tác pháp tuyến tính bởi, Fn  K n  Deq - Dab  (3) Fn lực tương tác pháp tuyến, Deq, Dab tương ứng khoảng cách ban đầu hai phần tử a b Do vật liệu đất nghiên cứu đất dính nên tồn lực dính hạt đất Trong mơ rời rạc, lực dính hạt đất giúp cho tương tác phần tử rời rạc chịu lực kéo (hình 2) Lực tương tác pháp tuyến giai đoạn chịu kéo tính với công thức (2) Tuy nhiên, tương tác bị phá hủy khoảng cách hai phần tử rời HVTH: NGUYỄN VĂN BÌNH Trang GVHD: TS TRẦN VĂN TIẾNG rạc vượt giá trị khoảng cách giới hạn kéo (Dab > Drupture); đó, lực tương tác không tồn Véc tơ lực tiếp tuyến Fs tính cập nhật định hướng nó, định hướng phụ thuộc vào hướng xuyên tâm hai phần tử rời rạc tương tác với nhau, cộng thêm gia số ΔFs (Hart et al 1988), gia số định nghĩa bởi: Fs  K s Us (4) Với ΔUs gia số véc tơ chuyển vị tiếp tuyến (thể thay đổi vị trí điểm tương tác hai phần tử sau bước thời gian Δt) Hình Lực tương tác pháp tuyến hai phần tử rời rạc Để mô ứng xử trạng thái tới hạn đất, tiêu chuẩn Mohr – Coulomb sử dụng (hình 3) Trị cực đại lực tương tác Fn,max định nghĩa hàm độ bền kéo T (độ bền kéo phụ thuộc vào góc ma sát lực dính) Lực tiếp tuyến cực đại Fs,max đặt trưng lực pháp tuyến Fn, lực dính C=T, góc ma sát tương tác Φc góc nội ma sát Φi (xem Tran et al 2011) Do vậy, lực pháp tuyến cực đại định nghĩa sau: Fn,max  TAint HVTH: NGUYỄN VĂN BÌNH (5) Trang GVHD: TS TRẦN VĂN TIẾNG ncb->engines.begin();vector::iterator itLast = ncb->engines.end();for ( ;itFirst!=itLast; ++itFirst ){ if ( ( *itFirst )->getClassName() == "TriaxialCompressionEngine" ){ cout checkpoint("Triangulating"); UpdateVolumes ( ncb ); timingDeltas>checkpoint("Update_Volumes"); ///Compute flow and and forces here flow->GaussSeidel ( ); timingDeltas->checkpoint("Gauss-Seidel"); // flow->MGPost(flow->T[flow>currentTes].Triangulation()); flow->tess_based_force=tess_based_force; flow->Compute_Forces ( ); timingDeltas>checkpoint("Compute_Forces"); ///End Compute flow and forces CGT::Finite_vertices_iterator vertices_end = flow->T[flow>currentTes].Triangulation().finite_vertices_end (); Vector3r f;int id; for ( CGT::Finite_vertices_iterator V_it = flow->T[flow>currentTes].Triangulation().finite_vertices_begin (); V_it != vertices_end; V_it++ ) HVTH: NGUYỄN VĂN BÌNH Trang 79 GVHD: TS TRẦN VĂN TIẾNG { id = V_it->info().id(); for ( int y=0;yinfo().forces ) [y]; ncb->forces.addForce ( id, f ); //ncb->forces.addTorque(id,t); } timingDeltas->checkpoint("Applying Forces"); Real time = Omega::instance().getSimulationTime(); int j = Omega::instance().getCurrentIteration(); // int j = Omega::instance().getSimulationTime(); char file [50]; string consol = flow->key+"%d_Consol"; const char* keyconsol = consol.c_str(); sprintf (file, keyconsol, j); char *g = file; timingDeltas->checkpoint("Writing cons_files"); MaxPressure = flow>PermeameterCurve(flow->T[flow->currentTes].Triangulation(), g, time); if ( Omega::instance().getCurrentIteration() % PermuteInterval == ) { Update_Triangulation = true; } timingDeltas->checkpoint("Storing Max Pressure"); // } } } } void FlowEngine::Oedometer_Boundary_Conditions() { flow->boundary ( flow->y_min_id ).flowCondition=0; flow->boundary ( flow->y_max_id ).flowCondition=0; flow->boundary ( flow->y_min_id ).value=0; flow->boundary ( flow->y_max_id ).value=0; triaxialCompressionEngine->wall_left_activated=0; triaxialCompressionEngine->wall_right_activated=0; triaxialCompressionEngine->wall_front_activated=0; triaxialCompressionEngine->wall_back_activated=0; triaxialCompressionEngine->wall_top_activated=1; triaxialCompressionEngine->wall_bottom_activated=1; HVTH: NGUYỄN VĂN BÌNH Trang 80 GVHD: TS TRẦN VĂN TIẾNG triaxialCompressionEngine>sigma_iso=(triaxialCompressionEngine->sigma_iso)*loadFactor; } void FlowEngine::Build_Triangulation ( Scene* ncb, double P_zero ) { if (first) { flow->currentTes=0; flow->Vtotalissimo=0; flow->Vsolid_tot=0; flow>Vporale=0; flow->Ssolid_tot=0; flow->SLIP_ON_LATERALS=slip_boundary; flow->key = triaxialCompressionEngine->Key; flow->TOLERANCE=Tolerance; } else { if (compute_K) {K = flow->Sample_Permeability ( flow->T[flow->currentTes].Triangulation(), flow->x_min, flow>x_max, flow->y_min, flow->y_max, flow->z_min, flow->z_max, flow->key );} flow->currentTes=!flow->currentTes; } currentTes=flow->currentTes; flow->x_min = 1000.0, flow->x_max = -10000.0, flow->y_min = 1000.0, flow->y_max = -10000.0, flow->z_min = 1000.0, flow>z_max = -10000.0; AddBoundary ( ncb ); Triangulate ( ncb ); cout k_factor = permeability_factor; flow->K_exp = K_exp; //flow->Compute_Permeability (); flow->Compute_Permeability_Darcy(); if (first) { CGT::Finite_cells_iterator cell_end = flow->T[flow>currentTes].Triangulation().finite_cells_end(); for ( CGT::Finite_cells_iterator cell = flow>T[flow->currentTes].Triangulation().finite_cells_begin(); cell != cell_end; cell++ ){cell->info().dv() = 0; cell->info().p() = 0;} if (compute_K) { K = flow->Sample_Permeability ( flow->T[flow->currentTes].Triangulation(), flow->x_min, flow>x_max, flow->y_min, flow->y_max, flow->z_min, flow->z_max, flow->key );} Oedometer_Boundary_Conditions(); first=!first; HVTH: NGUYỄN VĂN BÌNH Trang 81 GVHD: TS TRẦN VĂN TIẾNG } else { flow->Interpolate ( flow->T[!flow->currentTes], flow->T[flow->currentTes] ); Update_Triangulation=!Update_Triangulation; } flow->Initialize_pressures( P_zero ); Initialize_volumes ( ncb ); flow->DisplayStatistics (); } void FlowEngine::AddBoundary ( Scene* ncb ) { cout bodies ) { if ( !b ) continue; if ( b->shape->getClassIndex() == Bx_Index ) { Box* w = YADE_CAST ( b->shape.get() ); // const body_id_t& id = b->getId(); Real center [3], Extent[3]; for ( int h=0;hstate>pos[h]; Extent[h] = w->extents[h];} // ); flow->AddBoundingPlanes ( center, Extent, id flow->x_min = ( flow->x_min, center[0]- wall_thickness); flow->x_max center[0]+wall_thickness); flow->y_min wall_thickness); flow->y_max center[1]+wall_thickness); flow->z_min wall_thickness); flow->z_max center[2]+wall_thickness); } } = max ( flow->x_max, = ( flow->y_min, center[1]= max ( flow->y_max, = ( flow->z_min, center[2]= max ( flow->z_max, flow->AddBoundingPlanes(); } void FlowEngine::Triangulate ( Scene* ncb ) { shared_ptr sph ( new Sphere ); HVTH: NGUYỄN VĂN BÌNH Trang 82 GVHD: TS TRẦN VĂN TIẾNG int Sph_Index = sph->getClassIndexStatic(); int contator = 0; FOREACH ( const shared_ptr& b, *ncb->bodies ) { if ( !b ) continue; if ( b->shape->getClassIndex() == Sph_Index ) { Sphere* s=YADE_CAST ( b->shape.get() ); const body_id_t& id = b->getId(); Real rad = s->radius; Real x = b->state->pos[0]; Real y = b->state->pos[1]; Real z = b->state->pos[2]; flow->T[flow->currentTes].insert(x, y, z, rad, id); contator+=1; } } double SectionArea = ( flow->x_max - flow->x_min ) * ( flow->z_max-flow->z_min ); // cout cout cout cout cout cout cout cout info().dv() = ( Volume_cell_single_fictious ( cell,ncb )-cell->info().volume() ) /deltaT; cell->info().volume() = Volume_cell_single_fictious ( cell,ncb ); }break; case ( ) : { cell->info().dv() = ( Volume_cell ( cell,ncb )-cell->info().volume() ) /deltaT; cell->info().volume() = Volume_cell ( cell,ncb ); }break; } } } Real FlowEngine::Volume_cell_single_fictious ( CGT::Cell_handle cell, Scene* ncb ) { Real V[3][3]; HVTH: NGUYỄN VĂN BÌNH Trang 84 GVHD: TS TRẦN VĂN TIẾNG int b=0; int w=0; Real Wall_point[3]; for ( int y=0;yvertex ( y )->info().isFictious ) ) { const shared_ptr& sph = Body::byId ( cell->vertex ( y )->info().id(), ncb ); for ( int g=0;gstate>pos[g]; w++; } else { b = cell->vertex ( y )->info().id(); const shared_ptr& wll = Body::byId ( b , ncb ); for ( int i=0;iboundaries[b].p[i]; Wall_point[flow->boundaries[b].coordinate] = wll->state->pos[flow->boundaries[b].coordinate]-wall_thickness; } } double v1[3], v2[3]; for ( int g=0;gboundaries[b].normal ) * ( 0.33333333333* ( V[0][flow->boundaries[b].coordinate]+ V[1][flow->boundaries[b].coordinate]+ V[2][flow>boundaries[b].coordinate] ) - Wall_point[flow>boundaries[b].coordinate] ); // cout info().id(); const shared_ptr& wll = Body::byId ( b[j] , ncb ); for ( int i=0;iboundaries[b[j]].p[i]; Wall_point[j][flow>boundaries[b[j]].coordinate] = wll->state->pos[flow>boundaries[b[j]].coordinate] - wall_thickness; j++; } else if ( first_sph ) { const shared_ptr& sph1 = Body::byId ( cell->vertex ( g )->info().id(), ncb ); for ( int k=0;kstate->pos[k]; first_sph=false;} } else { const shared_ptr& sph2 = Body::byId ( cell->vertex ( g )->info().id(), ncb ); for ( int k=0;kstate->pos[k]; } } } AS[flow->boundaries[b[0]].coordinate]=BS[flow>boundaries[b[0]].coordinate] = Wall_point[0][flow>boundaries[b[0]].coordinate]; AT[flow->boundaries[b[1]].coordinate]=BT[flow>boundaries[b[1]].coordinate] = Wall_point[1][flow>boundaries[b[1]].coordinate]; for ( int h=0;hboundaries[b[0]].coordinate]- CS[flow>boundaries[b[0]].coordinate]; Real Volume = ( CGAL::cross_product ( v1,v2 ) *flow>boundaries[b[0]].normal ) *h; HVTH: NGUYỄN VĂN BÌNH Trang 86 GVHD: TS TRẦN VĂN TIẾNG // cout info().id(); const shared_ptr& wll = Body::byId ( b[j] , ncb ); for ( int i=0;iboundaries[b[j]].p[i]; Wall_point[j][flow>boundaries[b[j]].coordinate] = wll->state->pos[flow>boundaries[b[j]].coordinate]-wall_thickness; j++; } else { const shared_ptr& sph = Body::byId ( cell->vertex ( g )->info().id(), ncb ); for ( int k=0;kstate->pos[k];} } } AS[flow->boundaries[b[0]].coordinate]= AT[flow>boundaries[b[0]].coordinate]= AW[flow>boundaries[b[0]].coordinate]= Wall_point[0][flow>boundaries[b[0]].coordinate]; AT[flow->boundaries[b[1]].coordinate]= Wall_point[1][flow>boundaries[b[1]].coordinate]; AW[flow->boundaries[b[2]].coordinate]= Wall_point[2][flow>boundaries[b[2]].coordinate]; CGT::Vecteur v1 ( AS[0]-AT[0],AS[1]-AT[1],AS[2]-AT[2] ); CGT::Vecteur v2 ( AS[0]-AW[0],AS[1]-AW[1],AS[2]-AW[2] ); CGT::Vecteur h ( AT[0] - A[0], AT[1] - A[1], AT[2] - A[2] ); Real Volume = ( CGAL::cross_product ( v1,v2 ) ) * h; HVTH: NGUYỄN VĂN BÌNH Trang 87 GVHD: TS TRẦN VĂN TIẾNG // cout vel[1]; // // dz[m] = b->state->vel[2]; HVTH: NGUYỄN VĂN BÌNH Trang 88 GVHD: TS TRẦN VĂN TIẾNG // // // // ( v[m] ) [0] = b->state>pos[0]; // // ( v[m] ) [1] = b->state>pos[1]; // // ( v[m] ) [2] = b->state>pos[2]; // // // // } // // // // CGT::Vecteur v1 ( ( v[1] ) [0]v[0] ) [0], ( v[1] ) [1]- ( v[0] ) [1], ( v[1] ) [2]- ( v[0] [2] ); // // CGT::Vecteur v2 ( ( v[2] ) [0]v[1] ) [0], ( v[2] ) [1]- ( v[1] ) [1], ( v[2] ) [2]- ( v[1] [2] ); // // // // // // CGT::Vecteur V = 0.33333333333*CGT::Vecteur ( dx[0]+dx[1]+dx[2], dy[0]+dy[1]+dy[2], dz[0]+dz[1]+dz[2] ); // // CGT::Vecteur S = CGAL::cross_product ( v1,v2 ) /2.f; // // // // CGT::Somme ( grad_u, V, S ); // // } // // cell->info().dv() = grad_u.Trace(); // } // else // { ( ) ( ) // if ( triple_fictious ) // { // double deltaT = 1; // cell->info().dv() = (Volume_cell_double_fictious(cell,ncb)-cell>info().volume())/deltaT; // cell->info().volume() = Volume_cell_double_fictious(cell,ncb); /* int id_real_local=0,id_real_global=0,V_fict=0; double pos[3], surface=0; for ( int g=0;gvertex ( g )->info().isFictious ) { id_real_local=g; id_real_global=cell->vertex ( g )->info().id(); } } const shared_ptr& sph = Body::byId ( id_real_global, ncb ); for ( int i=0;istate->pos[i]; for ( int j=0;jvertex ( j )->info().isFictious ) { CGT::Boundary b = flow->boundaries[cell->vertex ( j )->info().id() ]; const shared_ptr& wall = Body::byId ( cell->vertex ( j )->info().id(), ncb ); surface = flow->surface_external_triple_fictious ( pos, cell, b ); Real Vs = sph->state->vel[b.coordinate]; Real Vw = wall->state->vel[b.coordinate]; Real Vrel = Vs - Vw; cell->info().dv() += Vrel*surface; } }*/ // } // if ( double_fictious ) // { // double deltaT = 1; // cell->info().dv() = (Volume_cell_double_fictious(cell,ncb)-cell>info().volume())/deltaT; // cell->info().volume() = Volume_cell_double_fictious(cell,ncb); // double A[3], AS[3], AT[3]; // double B[3], BS[3], BT[3]; // bool first=true, first_boundary=true; // Vector3r Vel_A, Vel_B, Vel_W1, Vel_W2; // // CGT::Boundary b1, b2; // // for ( int g=0;gvertex ( g )>info().isFictious && first ) // { // const shared_ptr& sph1 = Body::byId // ( cell->vertex ( g )->info().id(), ncb ); // for ( int y=0;ystate>pos[y]; AS[y]=A[y]; AT[y]=A[y];} // Vel_A = sph1->state->vel; // // first = false; // } // else if ( !cell->vertex ( g )>info().isFictious ) // { // const shared_ptr& sph2 = Body::byId // HVTH: NGUYỄN VĂN BÌNH Trang 90 GVHD: TS TRẦN VĂN TIẾNG ( cell->vertex ( g )->info().id(), ncb ); // for ( int y=0;yvertex ( g )->point() ) [y]; BS[y]=B[y]; BT[y]=B[y];} // // Vel_B = sph2->state->vel; // } // else if ( first_boundary ) // { // b1 = flow->boundaries[cell>vertex ( g )->info().id() ]; // const shared_ptr& wll1 = Body::byId // ( cell->vertex ( g )->info().id() , ncb ); // // Vel_W1=wll1->state->vel; // first_boundary=false; // } // else // { // b2 = flow->boundaries[cell>vertex ( g )->info().id() ]; // const shared_ptr& wll2 = Body::byId // ( cell->vertex ( g )->info().id() , ncb ); // // Vel_W2=wll2->state->vel; // } // } // // AS[b1.coordinate]=BS[b1.coordinate]=b1.p[b1.coordinate]; // AT[b2.coordinate]=BT[b2.coordinate]=b2.p[b2.coordinate]; // // double Vmoy[3]; // // for ( int y=0;yinfo().dv() += CGT::Vecteur ( Vmoy[0],Vmoy[1],Vmoy[2] ) *Surface; // // for ( int y=0;yinfo().dv() += CGT::Vecteur ( Vmoy[0],Vmoy[1],Vmoy[2] ) *Surface; // // for ( int y=0;yinfo().dv() += CGT::Vecteur ( Vmoy[0],Vmoy[1],Vmoy[2] ) *Surface; // // for ( int y=0;yinfo().dv() += CGT::Vecteur ( Vmoy[0],Vmoy[1],Vmoy[2] ) *Surface; // } // if ( single_fictious ) // { // double deltaT = 1; // cell->info().dv() = (Volume_cell_single_fictious(cell,ncb)-cell>info().volume())/deltaT; // cell->info().volume() = Volume_cell_single_fictious(cell,ncb); // } // } HVTH: NGUYỄN VĂN BÌNH Trang 92 S K L 0 ... DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN VĂN BÌNH XÂY DỰNG MƠ HÌNH KẾT HỢP PHẦN TỬ RỜI RẠC VÀ PHẦN TỬ CHẤT LỎNG NHẰM MƠ PHỎNG ỨNG XỬ CỦA ĐẤT DÍNH... rời rạc tương tác với Hình Mơ men lăn phần tử tương tác 2 .3 Mơ hình kết hợp chất lỏng - phần tử rời rạc Trong nghiên cứu này, mơ hình dịng chất lỏng (Tran et al 2016) áp dụng để kết hợp với mơ hình. .. VĂN TIẾNG 1 .3 Nhiệm vụ giới hạn đề tài - Xây dựng mơ hình ứng xử đất dính phương pháp phần tử rời rạc - Kết hợp với mơ hình chất lỏng cho nước lổ rỗng có sẵng - Số hóa mơ hình rời rạc ngơn ngữ