Soil-structure interaction with the contact algorithm
[email protected] Thu, 9 Jul 2015 23:18:33 -0500
| Newsgroups | gmane.comp.mathematics.tochnog.user |
|---|---|
| Message-ID | <[email protected]> |
--===============0217984259697615771== Content-Type: multipart/alternative; boundary="Apple-Mail=_776EBA42-3F47-4A89-BF24-A3BD56C532FF" --Apple-Mail=_776EBA42-3F47-4A89-BF24-A3BD56C532FF Content-Transfer-Encoding: quoted-printable Content-Type: text/plain; charset=utf-8 Hi Tochnogers, I=E2=80=99m trying to simulate a soil-structure interaction. It consists = of a thin-wall pipe embedded within a soil block, which has prescribed = displacement at its boundaries. Obviously, as the soil moves downwards, = and due to its lower stiffness, should =E2=80=98flow=E2=80=99 around the = pipe, but nonetheless inducing some stresses to it, which eventually = will end up in deformation of the pipe (way smaller than the soil = surrounding it). My initial idea was to use the contact algorithm because interfaces = haven=C2=B4t been implemented in tochnog yet. However, now I=E2=80=99m = not sure if it=E2=80=99s gonna work. I modified the file given in http://tochnog.sourceforge.net/ho_mech1.dat = <http://tochnog.sourceforge.net/ho_mech1.dat>, which is a rubber ball = being pressed downwards. I substituted the contact boundary by a block = of the same material and set the contact with the =E2=80=98node_boundary=E2= =80=99 option. It worked and the ball continues being deformed, but the = block below was=C2=B4t affected too much (some deformation was seen). = Therefore, I decided to modify the stiffness of the material of the ball = to be 100 times larger than the rubber, which now stays beneath the = ball. I was expecting a larger deformation of the block below due to its = smaller stiffness (100 lower than the ball pressing it). However, it = wasn=C2=B4t the case and the ball continues deforming more than the = block, which actually presents some deformations upwards!! Does anyone have an idea of how can I have this =E2=80=98two-bodies=E2=80=99= interaction. I=E2=80=99m attaching the modified data file if someone = might find some errors. Thank you very much in advance for your help. cheers, Ferney echo -yes number_of_space_dimensions 2 materi_velocity materi_displacement materi_stress materi_strain_elasti materi_strain_total end_initia node 0 0. -0.2 node 1 0. -0.1 node 2 0. 0.0 node 3 .1 0.0 node 4 .2 0.0 node 5 .1414 -0.14142 node 6 .0805 -0.06667 node 7 0 -.21 node 8 .3 -.21 node 9 .3 -.4 node 10 0 -.4 element 0 -quad4 0 5 1 6 element 1 -quad4 5 4 6 3 element 2 -quad4 1 6 2 3 element 3 -quad4 10 9 7 8 element_group 0 1 element_group 1 1 element_group 2 1 element_group 3 0 group_type 0 -materi group_materi_memory 0 -total group_materi_elasti_young 0 2.82 group_materi_elasti_poisson 0 0.49967 group_type 1 -materi group_materi_memory 1 -total group_materi_elasti_young 1 282 group_materi_elasti_poisson 1 0.249967 (left edge) geometry_line 0 0.0 -0.2 0.0 0.0 1.e-3 (upper edge) geometry_line 1 0.0 0.0 0.2 0.0 1.e-3 (contact line) geometry_line 2 0.0 -0.2 1.0 -0.2 1.e-2 geometry_line 3 0 -.4 1 -.4 0.001 geometry_line 4 0 -.21 0 -.4 0.001 bounda_unknown 0 -geometry_line 0 -velx bounda_time 0 0. 0. 100. 0. bounda_unknown 1 -geometry_line 1 -vely bounda_time 1 0.0 -0.15 100.0 -0.15 bounda_unknown 2 -geometry_line 3 -vely bounda_time 2 0.0 -0.0 100.0 -0.0 bounda_unknown 3 -geometry_line 4 -velx bounda_time 3 0. 0. 100. 0. node_boundary -ra 7 8 -ra -yes contact_penalty_velocity 1.e5 contact_geometry 1 -geometry_line 2 post_node 1 -node_rhside -sum -geometry_line 1 control_mesh_refine_globally 7 -h_refinement control_mesh_refine_globally_geometry 7 -geometry_set 0 control_mesh_refine_globally 8 -h_refinement control_mesh_refine_globally_geometry 8 -geometry_set 0 control_timestep 10 0.01 1. control_timestep_iterations 10 10 control_print_gid 10 -yes target_item 1 -post_node_result 1 -vely=20 target_value 1 1.41727 0.02 end_data --Apple-Mail=_776EBA42-3F47-4A89-BF24-A3BD56C532FF Content-Transfer-Encoding: quoted-printable Content-Type: text/html; charset=utf-8 <html><head><meta http-equiv=3D"Content-Type" content=3D"text/html = charset=3Dutf-8"></head><body style=3D"word-wrap: break-word; = -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" = class=3D"">Hi Tochnogers,<div class=3D""><br class=3D""></div><div = class=3D"">I=E2=80=99m trying to simulate a soil-structure interaction. = It consists of a thin-wall pipe embedded within a soil block, which has = prescribed displacement at its boundaries. Obviously, as the soil moves = downwards, and due to its lower stiffness, should =E2=80=98flow=E2=80=99 = around the pipe, but nonetheless inducing some stresses to it, which = eventually will end up in deformation of the pipe (way smaller than the = soil surrounding it).</div><div class=3D""><br class=3D""></div><div = class=3D"">My initial idea was to use the contact algorithm because = interfaces haven=C2=B4t been implemented in tochnog yet. However, now = I=E2=80=99m not sure if it=E2=80=99s gonna work.</div><div class=3D""><br = class=3D""></div><div class=3D"">I modified the file given in <a = href=3D"http://tochnog.sourceforge.net/ho_mech1.dat" = class=3D"">http://tochnog.sourceforge.net/ho_mech1.dat</a>, which = is a rubber ball being pressed downwards. I substituted the contact = boundary by a block of the same material and set the contact with the = =E2=80=98node_boundary=E2=80=99 option. It worked and the ball continues = being deformed, but the block below was=C2=B4t affected too much (some = deformation was seen). Therefore, I decided to modify the stiffness of = the material of the ball to be 100 times larger than the rubber, which = now stays beneath the ball.</div><div class=3D""><br class=3D""></div><div= class=3D"">I was expecting a larger deformation of the block below due = to its smaller stiffness (100 lower than the ball pressing it). However, = it wasn=C2=B4t the case and the ball continues deforming more than = the block, which actually presents some deformations upwards!!</div><div = class=3D""><br class=3D""></div><div class=3D"">Does anyone have an idea = of how can I have this =E2=80=98two-bodies=E2=80=99 interaction. I=E2=80=99= m attaching the modified data file if someone might find some = errors.</div><div class=3D""><br class=3D""></div><div class=3D""><br = class=3D""></div><div class=3D""><div class=3D"">Thank you very much in = advance for your help.</div><div class=3D""><br class=3D""></div><div = class=3D"">cheers,</div><div class=3D""><br class=3D""></div><div = class=3D"">Ferney</div></div><div class=3D""><br class=3D""></div><div = class=3D""><br class=3D""></div><div class=3D""><br class=3D""></div><div = class=3D""><div class=3D"">echo -yes</div><div = class=3D"">number_of_space_dimensions 2</div><div = class=3D"">materi_velocity</div><div = class=3D"">materi_displacement</div><div = class=3D"">materi_stress</div><div = class=3D"">materi_strain_elasti</div><div = class=3D"">materi_strain_total</div><div class=3D"">end_initia</div><div = class=3D""><br class=3D""></div><div class=3D"">node 0 0. = -0.2</div><div class=3D"">node 1 0. -0.1</div><div = class=3D"">node 2 0. 0.0</div><div class=3D"">node = 3 .1 0.0</div><div class=3D"">node 4 .2 = 0.0</div><div class=3D"">node 5 .1414 = -0.14142</div><div class=3D"">node 6 .0805 = -0.06667</div><div class=3D"">node 7 0 -.21</div><div class=3D"">node 8 = .3 -.21</div><div class=3D"">node 9 .3 -.4</div><div class=3D"">node 10 = 0 -.4</div><div class=3D""><br class=3D""></div><div class=3D"">element = 0 -quad4 0 5 1 6</div><div class=3D"">element 1 -quad4 5 4 6 3</div><div = class=3D"">element 2 -quad4 1 6 2 3</div><div class=3D"">element 3 = -quad4 10 9 7 8</div><div class=3D""><br class=3D""></div><div = class=3D"">element_group 0 1</div><div class=3D"">element_group 1 = 1</div><div class=3D"">element_group 2 1</div><div = class=3D"">element_group 3 0</div><div class=3D""><br = class=3D""></div><div class=3D"">group_type 0 -materi</div><div = class=3D"">group_materi_memory 0 -total</div><div = class=3D"">group_materi_elasti_young 0 2.82</div><div = class=3D"">group_materi_elasti_poisson 0 0.49967</div><div = class=3D""><br class=3D""></div><div class=3D"">group_type 1 = -materi</div><div class=3D"">group_materi_memory 1 = -total</div><div class=3D"">group_materi_elasti_young 1 = 282</div><div class=3D"">group_materi_elasti_poisson 1 = 0.249967</div><div class=3D""><br class=3D""></div><div = class=3D"">(left edge)</div><div class=3D"">geometry_line 0 0.0 = -0.2 0.0 0.0 1.e-3</div><div class=3D""><br class=3D""></div><div = class=3D"">(upper edge)</div><div class=3D"">geometry_line 1 0.0 = 0.0 0.2 0.0 1.e-3</div><div class=3D""><br = class=3D""></div><div class=3D"">(contact line)</div><div = class=3D"">geometry_line 2 0.0 -0.2 1.0 -0.2 1.e-2</div><div = class=3D""><br class=3D""></div><div class=3D""><br class=3D""></div><div = class=3D"">geometry_line 3 0 -.4 1 -.4 0.001</div><div = class=3D"">geometry_line 4 0 -.21 0 -.4 0.001</div><div class=3D""><br = class=3D""></div><div class=3D"">bounda_unknown 0 -geometry_line 0 = -velx</div><div class=3D"">bounda_time 0 0. 0. 100. = 0.</div><div class=3D"">bounda_unknown 1 -geometry_line 1 = -vely</div><div class=3D"">bounda_time 1 0.0 -0.15 = 100.0 -0.15</div><div class=3D"">bounda_unknown 2 -geometry_line 3 = -vely</div><div class=3D"">bounda_time 2 0.0 -0.0 = 100.0 -0.0</div><div class=3D"">bounda_unknown 3 -geometry_line 4 = -velx</div><div class=3D"">bounda_time 3 0. 0. 100. = 0.</div><div class=3D""><br class=3D""></div><div class=3D"">node_boundary= -ra 7 8 -ra -yes</div><div class=3D""><br class=3D""></div><div = class=3D"">contact_penalty_velocity 1.e5</div><div = class=3D"">contact_geometry 1 -geometry_line 2</div><div = class=3D""><br class=3D""></div><div class=3D"">post_node 1 -node_rhside = -sum -geometry_line 1</div><div class=3D""><br class=3D""></div><div = class=3D"">control_mesh_refine_globally = 7 -h_refinement</div><div = class=3D"">control_mesh_refine_globally_geometry 7 = -geometry_set 0</div><div class=3D"">control_mesh_refine_globally = 8 -h_refinement</div><div = class=3D"">control_mesh_refine_globally_geometry 8 = -geometry_set 0</div><div class=3D""><br class=3D""></div><div = class=3D"">control_timestep = 10 0.01 1.</div><div = class=3D"">control_timestep_iterations = 10 10</div><div class=3D""><br class=3D""></div><div = class=3D""><br class=3D""></div><div class=3D"">control_print_gid<span = class=3D"Apple-tab-span" style=3D"white-space:pre"> = </span>10 -yes</div><div class=3D""><br class=3D""></div><div = class=3D"">target_item 1 -post_node_result 1 -vely </div><div = class=3D"">target_value 1 1.41727 0.02</div><div = class=3D"">end_data</div></div><div class=3D""><br class=3D""></div><div = class=3D""><br class=3D""></div><div class=3D""><br = class=3D""></div></body></html>= --Apple-Mail=_776EBA42-3F47-4A89-BF24-A3BD56C532FF-- --===============0217984259697615771== Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit Content-Disposition: inline ------------------------------------------------------------------------------ Don't Limit Your Business. Reach for the Cloud. GigeNET's Cloud Solutions provide you with the tools and support that you need to offload your IT needs and focus on growing your business. Configured For All Businesses. Start Your Cloud Today. https://www.gigenetcloud.com/ --===============0217984259697615771== Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit Content-Disposition: inline _______________________________________________ Tochnog-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/tochnog-users --===============0217984259697615771==--