<div dir="ltr"><div>Hi Gaetano, <br></div><div><br></div><div>I'm glad it's working for you! <br></div><div>Yes, I think the same code should work, but the data labels will be incorrect for this new batch. </div><div>You can change the calls to gmsh.view.add( name ) to fix this. <br></div><div><br></div><div>Sincerely, <br></div><div>Max <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 19, 2019 at 10:02 AM Gaetano Camarda <<a href="mailto:gaetano_18@live.it">gaetano_18@live.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div>


<div class="gmail-m_5527032270517882208WordSection1">
<p class="MsoNormal">Hi Max,</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thank you again for your code, I’ve tried it on my unitcell and it work great, I’ll send you a picture just for letting you see what I was telling about.
</p>
<p class="MsoNormal">I would like to ask you another question, I have same information with deformation (strain) information. The matrix is same way as stress:</p>
<p class="MsoNormal">[ epsx, epsy epsz epsxy epsyz epsxz]</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I was wondering if you think I can use same code for doing this operation, what do you think?</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Sincerily,</p>
<p class="MsoNormal">Gaetano</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Inviato da <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank">
Posta</a> per Windows 10</p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5527032270517882208divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>Da:</b> Max Orok <<a href="mailto:morok@mevex.com" target="_blank">morok@mevex.com</a>><br>
<b>Inviato:</b> Saturday, February 16, 2019 5:35:43 PM<br>
<b>A:</b> Gaetano Camarda; <a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
<b>Oggetto:</b> Re: [Gmsh] Stress on Gauss Points</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>Hi Gaetano, <br>
</div>
<div><br>
</div>
<div>Here I modified the script to each of the stress values separately as scalars.
<br>
</div>
<div>The views will all be visible at once by default but can be toggled by clicking the "Post-processing" text in the sidebar and then cycled through.
<br>
</div>
<div><br>
</div>
<div>
<div>About that displacement question, I just meant I would defer to your judgment as to whether the displacement matched the stresses, etc.
<br>
</div>
<div><br>
</div>
<div>Sincerely, <br>
</div>
<div>Max</div>
<div><br>
</div>
<div>
<div><img alt="image.png" src="cid:16907e1e525cb971f161" width="474" height="317"></div>
<div><br>
</div>
<div>
<div><img alt="image.png" src="cid:16907e1e525cb971f162" width="474" height="368"></div>
<div>
<div><img alt="image.png" src="cid:16907e1e525cb971f163" width="474" height="399"><br>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sat, Feb 16, 2019 at 6:30 AM Gaetano Camarda <<a href="mailto:gaetano_18@live.it" target="_blank">gaetano_18@live.it</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div class="gmail-m_5527032270517882208gmail-m_3223345452502574116WordSection1">
<p class="MsoNormal">Hi guys,</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Yes, the goal for me was to plot the data on Gmsh, but i would like to ask if there’s a way to show sigx, sigy and so on as view in Gmsh Post-processing, it’s look like all the tensor is plotted all together, am I wrong?</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I will do a little brief on my work, maybe it could be better understanding. I’m sorry in advance if my english is not so perfect.</p>
<p class="MsoNormal">I’m working on a Multiscale RVE, not the one I’,ve send but another one more complex, but I’m working in parallel with this smaller one for having more control on the code and on the operations.</p>
<p class="MsoNormal">First thing first, I’ve generated the mesh on Gmsh and I processed the operation on Matlab, part of the code was made by me, and for the solver I’ve used CalFem function (freeware script for matlab and python for studying FEM). The .txt
 file with stress matrix in gauss point and nodal forces for each element are the result of the Calfem function’s I’ve used. The Gauss stress matrix I’ve sent as .txt has 8 point for each element, this is why it has 512 row, cause the element are 64 with 8
 gauss point = 512.</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Now the post-processing part, my goal is to plot the data on Gmsh postprocessing, thanks to Max now I’ve understood a little bit more about the input data Gmsh needed, but I would like to plot stress field separately (as sigx view, sigy
 view and so on); I also have another matrix for strain, but I didn’t sent because I think the process to do is the same for the gauss-points stresses.
</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">About this question from Max:</p>
<p class="MsoNormal">“<span style="font-size:11.5pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33);background:white none repeat scroll 0% 0%">For the displacement, does it seem like it aligns? I don't want to make a mistake in interpreting the data. “<u></u><u></u></span></p>
<p class="MsoNormal">I haven’t understand well your question, my english is a little bit limited, what to you mean with “it seem like it alligns?” If I’ve understood the question you mean the fact that displacement are not the same in faces dx and sx? if that’s
 your meaning, i’ve noticed it too, and I think it depends on the boundary conditions I gave to the element.</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Sincerily,</p>
<p class="MsoNormal">Gaetano</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Inviato da <a href="https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgo.microsoft.com%2Ffwlink%2F%3FLinkId%3D550986&data=02%7C01%7C%7C73310a179a5841aab52408d6942cd9da%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636859317703189999&sdata=G7qtTFXlvkUpFI8GnWPOiN4%2B6AtBN1Baj1qFK%2F%2F3sfE%3D&reserved=0" target="_blank">
Posta</a> per Windows 10</p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5527032270517882208gmail-m_3223345452502574116divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>Da:</b> Jeremy Theler <<a href="mailto:jeremy@seamplex.com" target="_blank">jeremy@seamplex.com</a>><br>
<b>Inviato:</b> Saturday, February 16, 2019 11:29:03 AM<br>
<b>A:</b> Max Orok<br>
<b>Cc:</b> Gaetano Camarda; <a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
<b>Oggetto:</b> Re: [Gmsh] Stress on Gauss Points</font>
<div> </div>
</div>
<div>
<div>Hey Max</div>
<div><br>
</div>
<div>On Fri, 2019-02-15 at 11:01 -0500, Max Orok wrote:</div>
<div><br>
</div>
<div><br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>My impression going in was that the Gauss point steps were part of the simulation and the goal was just to plot the data.
<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Well, there is a fine line here but it can be considered that the modeling [1] is finished when the stresses (or the heat fluxes, etc) are computed at the gauss points.</div>
<div><br>
</div>
<div>Actually, I think the original question was: "given that I have data at the Gauss points (and not at the nodes nor at the cell centers), how can I plot them on Gmsh?"</div>
<div><br>
</div>
<div>Indeed, I think Gaetano wanted to have Gmsh to do something like to re-mesh the base grid and have new cells around the gauss points and plot the data as cell-centered values. I do not know if I am clear, but for a 8-node hexa with 2 Gauss points per dimensions
 that would be to replace the original element with 8 little hexas, one in each octant and then assign a cell value from the dataset to each of them. For other elements it would be more complex.</div>
<div><br>
</div>
<div>This is definetely in the post-processing step of the worflow.</div>
<div><br>
</div>
<div><br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>I think for this data set, the shared node data does repeat, but a single element will vary as a function of its nodes.
<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>As far as I can see, the data corresponding to the shared nodes does not repeat. Again, this is what is expected from secondary fields computed as derivatives of a primary field at the Gauss points.</div>
<div><br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div><br>
</div>
<div>For any averaging etc., perhaps this could be done during the simulation? </div>
</div>
</blockquote>
<div><br>
</div>
<div>Yes, it could. This is why I am interested in this issue, and raised the question on how to average if the solver has to output nodal values of stresses. Weighted average? With the element's volume? With the element's quality?</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>[1] <a href="https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.seamplex.com%2Fblog%2Fsay-modeling-not-simulation.html&data=02%7C01%7C%7C73310a179a5841aab52408d6942cd9da%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636859317703200010&sdata=HyJU40M0h8b0L3zxRw0c6j%2FCJa%2FpwWvLqSSaoxKYNHA%3D&reserved=0" target="_blank">
https://www.seamplex.com/blog/say-modeling-not-simulation.html</a></div>
<div><br>
</div>
<div><br>
</div>
<div>--</div>
<div>Jeremy Theler</div>
<div><a href="https://eur02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.seamplex.com&data=02%7C01%7C%7C73310a179a5841aab52408d6942cd9da%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636859317703210021&sdata=xGwsNp%2BFJ0Z6vASU1zUR5eXqaPNz4pRMgTya%2F%2FxNFks%3D&reserved=0" target="_blank">www.seamplex.com</a></div>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<br>
-- <br>
<div dir="ltr" class="gmail-m_5527032270517882208gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>Max Orok</div>
<div>Contractor<br>
</div>
</div>
<div dir="ltr"><a href="https://eur02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.mevex.com&data=02%7C01%7C%7C73310a179a5841aab52408d6942cd9da%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636859317703220026&sdata=v%2FR9k2G3Zec8F1qQcL%2BvNdMB89wmC%2B%2FgNEmzSV%2F3UCQ%3D&reserved=0" target="_blank">www.mevex.com</a><br>
<br>
<img width="164" height="42"><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>

</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Max Orok</div><div>Contractor<br></div></div><div dir="ltr"><a href="http://www.mevex.com" target="_blank">www.mevex.com</a><br><br><img src="https://docs.google.com/uc?export=download&id=1fHTIiW4OMUjQr1iOkspQ7wiEsxunoOs0&revid=0B6x5w-5zVaEjSkpwbm5oY29jbG1XMzJoYldXTmJpNGFtb3dVPQ" width="164" height="42"><br></div></div></div></div></div></div>