<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Vijay,<div class=""><br class=""></div><div class="">Thank you very much for your email. </div><div class=""><br class=""></div><div class="">We are trying to track the deformed configuration and write it into .geo file for GMSH to read and to generate the new mesh.</div><div class=""><br class=""></div><div class="">Now it can track some very simple geometry models based on the ‘Spline’ command in GMSH. </div><div class=""><br class=""></div><div class="">We have not yet deal with the complex geometries now (like with internal boundaries, curves), but we are trying to make the code extensible so it is re-useable for the future work.</div><div class=""><br class=""></div><div class="">Another difficulty in my head is automatically applying the new boundary conditions after deformation. Because the boundary conditions applied usually (in FEAP) based on the node indexes or node coordinates, and both of these 2 terms changed on the new mesh. For some simple boundary conditions, this is not very difficult.</div><div class=""><br class=""></div><div class="">Now the step we have not yet done anything is the transformation of variables between meshes. With that step implemented, the re-meshing algorithm could be completed. We will try to see if we can write a .vtu reader and writer to deal with this problem.</div><div class=""><br class=""></div><div class="">Thank you again and let's keep in touch!</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Jie</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">
<br class=""><div><blockquote type="cite" class=""><div class="">On Jun 2, 2016, at 10:28 PM, Vijay S. Mahadevan <<a href="mailto:vijay.m@gmail.com" class="">vijay.m@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><p dir="ltr" class="">Hi Jie, </p><p dir="ltr" class="">Thanks for detailing your workflow so that we understand your requirements better. I've considered adding support for vtu before but wasn't sure if we could do it without adding dependency on vtk. Our configuration already includes vtk support however, for our Paraview visualization plugin. So this may simplify things a bit. If you have the time to write the text/binary VTU I/O implementations (as an add-on to VTK I/O), we would be more than happy to take a PR submission for it. We can also help you along the way to get the implementation refined and tested. </p><p dir="ltr" class="">Regarding your workflow, we are attempting something similar with a spectral element code right now with conforming quad/hex local refinements, in combination with mesh optimization. This would be the task where you are using Gmsh in your step 2.</p><p dir="ltr" class="">I would be interested to know how your workflow is progressing and what difficulties, if any, you encounter in the mesh regeneration stages, especially with propagating deformations back to the geometry model. Perhaps with simpler domains, this may not prove to be a bottleneck.</p><p dir="ltr" class="">Anyway, do let us know if you have more questions. </p><p dir="ltr" class="">Best,<br class="">
Vijay</p>
<div class="gmail_quote">On Jun 2, 2016 5:24 PM, "Jie Wu" <<a href="mailto:jie.voo@gmail.com" class="">jie.voo@gmail.com</a>> wrote:<br type="attribution" class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Hi Iulian,<div class=""><br class=""></div><div class="">Thank you very much for your instruction. After talking about this with my colleagues, we think about writing a .vtu reader for MOAB.</div><div class=""><br class=""></div><div class="">And the overall workflow looks like this:</div><div class=""><ol class=""><li class="">Use python to call Gmsh to generate finite element mesh;</li><li class="">Prepare input files for FEAP (Finite Element Analysis Program);</li><li class="">While load steps not yet finished, Do:</li><ol class=""><li class="">Use python to call FEAP to run the calculation, and generate .vtu file based on the calculation results;</li><li class="">Generate a new mesh, with Gmsh, based on the deformed configuration;</li><li class="">Use the MOAB .vtu reader to read .vtu file and pass the information to mbcoupler;</li><li class="">Transfer data from old to new mesh by mbcoupler, and write the data on the new mesh into a .vtu file.</li><li class="">Prepare input files for FEAP based on the new .vtu file;</li></ol></ol><div class="">Run simulation using a mesh generated with GMSH and write vtu files to export the solution of the simulation.</div><div class="">When the mesh is too distorted, stop the simulation, generate a new mesh with GMSH and transfer the data from the last vtu file to the new mesh.</div><div class="">mbcoupler would be used to perform the transfer of the data from the last vtu to the new mesh.</div><div class=""><br class=""></div><div class="">Any suggestion about this work is greatly appreciated!</div></div><div class=""><br class=""></div><div class="">Best,</div><div class="">Jie</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">
<br class=""><div class=""><blockquote type="cite" class=""><div class="">On Jun 1, 2016, at 7:44 PM, Grindeanu, Iulian R. <<a href="mailto:iulian@mcs.anl.gov" target="_blank" class="">iulian@mcs.anl.gov</a>> wrote:</div><br class=""><div class=""><div style="font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;direction:ltr;font-family:Tahoma;font-size:10pt" class="">Hi<br class="">vtk can be unstructured too;<br class="">the issue is that the vtk file converted from your vtu file using VisIt cannot be read directly by MOAB.<span class=""> </span><br class="">It would make more sense to make a vtu reader in MOAB, or make a h5m writer in VisIt.<br class="">Actually, I am working on that (h5m writer) , but it is in a branch now, it is not part of visit yet.<br class=""><br class="">I am not sure what is the best solution to your problem now.<span class=""> </span><br class="">How do you plan to use the coupler?<span class=""> </span><br class=""><br class="">Iulian<br class=""><br class=""><div style="font-family:'Times New Roman';font-size:16px" class=""><hr class=""><div style="direction:ltr" class=""><font face="Tahoma" size="2" class=""><b class="">From:</b><span class=""> </span>Jie Wu [<a href="mailto:jie.voo@gmail.com" target="_blank" class="">jie.voo@gmail.com</a>]<br class=""><b class="">Sent:</b><span class=""> </span>Wednesday, June 01, 2016 5:07 PM<br class=""><b class="">To:</b><span class=""> </span>Grindeanu, Iulian R.<br class=""><b class="">Cc:</b><span class=""> </span><a href="mailto:moab-dev@mcs.anl.gov" target="_blank" class="">moab-dev@mcs.anl.gov</a><br class=""><b class="">Subject:</b><span class=""> </span>Re: [MOAB-dev] How can I prepare a .h5m file base on my own calculation<br class=""></font><br class=""></div><div class=""></div><div class=""><div style="word-wrap:break-word" class="">Hi Iulian,<div class=""><br class=""></div><div class="">I will try to convert from .vtu to .vtk by the tools you mentioned.</div><div class=""><br class=""></div><div class="">But as I know, .vtk file is usually used for structured mesh, while .vtu for unstructured mesh. I was just wondering if the conversion from .vtu to .vtk make sense?</div><div class=""><br class=""></div><div class="">Here is a small vtu file I used to store mesh and solution information. </div><div class=""><br class=""></div><div class="">On the nodes, we have displacements and a scaler damage parameter. We also have some other solution variables on the Integration nodes.</div><div class=""><br class=""></div><div class="">Thank you so much for your help.</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Jie</div><div class=""><br class=""></div><div class=""></div></div><div style="word-wrap:break-word" class=""><div class=""></div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Jun 1, 2016, at 5:31 PM, Grindeanu, Iulian R. <<a href="mailto:iulian@mcs.anl.gov" target="_blank" class="">iulian@mcs.anl.gov</a>> wrote:</div><br class=""><div class=""><div style="font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;direction:ltr;font-family:Tahoma;font-size:10pt" class="">Hi,<br class="">You can try paraview to convert from vtu format to vtk;<br class="">Moab vtk reader is not very robust, in the sense that there are vtk files it cannot read.<br class="">You can also try VisIt to convert from vtu to vtk<br class=""><br class="">Do you have a small vtu file that you can share? I can try to look at it, to see how it goes.<br class="">What kind of data do you have at nodes?<br class=""><br class="">Iulian<br class=""><br class=""><div style="font-family:'Times New Roman';font-size:16px" class=""><hr class=""><div style="direction:ltr" class=""><font face="Tahoma" size="2" class=""><b class="">From:</b><span class=""> </span>Jie Wu [<a href="mailto:jie.voo@gmail.com" target="_blank" class="">jie.voo@gmail.com</a>]<br class=""><b class="">Sent:</b><span class=""> </span>Wednesday, June 01, 2016 2:42 PM<br class=""><b class="">To:</b><span class=""> </span>Grindeanu, Iulian R.<br class=""><b class="">Cc:</b><span class=""> </span><a href="mailto:moab-dev@mcs.anl.gov" target="_blank" class="">moab-dev@mcs.anl.gov</a><br class=""><b class="">Subject:</b><span class=""> </span>Re: [MOAB-dev] How can I prepare a .h5m file base on my own calculation<br class=""></font><br class=""></div><div class=""></div><div class="">Hi Iulian,<div class=""><br class=""></div><div class="">Thank you very much for your help. The format we use to store mesh and solution information is .vtu file. </div><div class=""><br class=""></div><div class="">I tried to run function ./mbconvert. It can successfully transfer from .h5m to .vtk, and vice versa. </div><div class=""><br class=""></div><div class="">But I cannot convert either these 2 format to/from .vtu file. </div><div class=""><br class=""></div><div class="">An example error message is like following. </div><div class=""><br class=""></div><div class="">Is there any solution for this problem? Any suggestion is greatly appreciated!</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Jie</div><div class=""><br class=""></div><div class=""><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">dyn-160-39-11-225:tools jiewu$ ./mbconvert in.vtu out.h5m</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: in.vtu: File is not HDF5!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: set_up_read() line 472 in src/io/ReadHDF5.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: NULL file handle.!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: is_error() line 133 in src/io/ReadHDF5.hpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: ReadNCDF:: problem opening Netcdf/Exodus II file in.vtu!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: load_file() line 303 in src/io/ReadNCDF.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: Trouble opening file in.vtu!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: load_file() line 68 in src/io/ReadNC.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: Expected Keyword!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: load_file() line 178 in src/io/ReadABAQUS.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">0x0 hahhahahah</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: This doesn't appear to be a .cub file!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: load_file() line 310 in src/io/Tqdcfr.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: Parsing error at line 1: expected one of { $NOD $MeshFormat } got "<?xml"!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: match_token() line 378 in src/io/FileTokenizer.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: in.vtu: No such file or directory!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: ascii_read_triangles() line 201 in src/io/ReadSTL.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: in.vtu: not a binary STL file!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: binary_read_triangles() line 312 in src/io/ReadSTL.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: in.vtu.node: cannot read file!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: open_file() line 59 in src/io/ReadTetGen.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: in.vtu: Trouble reading vertices!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: read_vertices() line 129 in src/io/ReadTemplate.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: --------------------- Error Message ------------------------------------</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: Failed to load file after trying all possible readers!</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: serial_load_file() line 635 in src/Core.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">[0]MOAB ERROR: load_file() line 526 in src/Core.cpp</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">Failed to load "in.vtu".</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">Error code: MB_INDEX_OUT_OF_RANGE (1)</span></div><div style="margin:0px;font-size:14px;line-height:normal;font-family:'Andale Mono';color:rgb(254,253,201);background-color:rgb(0,0,0)" class=""><span class="">Error message: Failed to load file after trying all possible readers</span></div></div><div class=""><span class=""><br class=""></span></div><div class=""><br class=""></div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Jun 1, 2016, at 12:26 PM, Grindeanu, Iulian R. <<a href="mailto:iulian@mcs.anl.gov" target="_blank" class="">iulian@mcs.anl.gov</a>> wrote:</div><br class=""><div class=""><div style="font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;direction:ltr;font-family:Tahoma;font-size:10pt" class="">Hello,<br class="">Coupler.cpp depends on a moab instance, so even if you do not use h5m format, if you want to use the coupler, your model has to reside inside a moab instance.<span class=""> </span><br class="">Do you care about parallelism or not? Then h5m format is required, otherwise it will be hard to employ parallelism<br class=""><br class="">there are several converters available, to h5m format, you can do tools/mbconvert -l to see what is available based on your configuration<br class=""><br class="">What format do you use?<br class=""><br class="">Iulian<br class=""><div style="font-family:'Times New Roman';font-size:16px" class=""><hr class=""><div style="direction:ltr" class=""><font face="Tahoma" size="2" class=""><b class="">From:</b><span class=""> </span><a href="mailto:moab-dev-bounces@mcs.anl.gov" target="_blank" class="">moab-dev-bounces@mcs.anl.gov</a><span class=""> </span>[<a href="mailto:moab-dev-bounces@mcs.anl.gov" target="_blank" class="">moab-dev-bounces@mcs.anl.gov</a>] on behalf of Jie Wu [<a href="mailto:jie.voo@gmail.com" target="_blank" class="">jie.voo@gmail.com</a>]<br class=""><b class="">Sent:</b><span class=""> </span>Wednesday, June 01, 2016 9:17 AM<br class=""><b class="">To:</b><span class=""> </span><a href="mailto:moab-dev@mcs.anl.gov" target="_blank" class="">moab-dev@mcs.anl.gov</a><br class=""><b class="">Subject:</b><span class=""> </span>[MOAB-dev] How can I prepare a .h5m file base on my own calculation<br class=""></font><br class=""></div><div class=""></div><div class="">Hi All,<div class=""><br class=""></div><div class="">In order to transfer variables from a Source mesh to a Target mesh, I am studying the use of function<span class=""> </span><i class="">mbcoupler_test.cpp</i><span class=""> </span>these days. </div><div class=""><br class=""></div><div class="">The stuff confusing me is the .h5m file. From the MOAB user guide, I know that .h5m is theMOAB’s native file format. </div><div class=""><br class=""></div><div class="">But I don’t know how to prepare a .h5m file based on my own calculation results (which usually contains node coordinates, node connectivities, and nodal variables). As a result, I cannot employ a function like<span class=""> </span><i class="">./mbcoupler_test</i><span class=""> </span>to help me to transfer the mesh variables.</div><div class=""><br class=""></div><div class="">Any instruction is deeply appreciated!</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Jie</div><div class=""><br class=""></div><br class=""><div class=""><div class=""><div class="">Jie Wu, Ph.D. candidate</div><div class="">Department of Geotechnical Engineering, Tongji University</div><div class="">Visiting student at Civil Engineering & Engineering Mechanics, Columbia University <br class="">Add: 500 West 120th Street 641 S.W. Mudd New York, NY 10027<br class="">Ph: <a href="tel:%2B1%28917%29822-0942" value="+19178220942" target="_blank" class="">+1(917)822-0942</a></div></div></div></div></div></div></div></blockquote></div></div></div></div></div></div></blockquote></div></div></div></div></div></div></div></blockquote></div><br class=""></div></div></blockquote></div>
</div></blockquote></div><br class=""></div></body></html>