Hi,
I'm trying to solve a large deformation problem interactively using updated
lagrangian formulation. I have some trouble in updating the mesh. In the
code I attached here, I attempted to define new_ulf_iteration in
problem.nls_iter_hook (line 149) and nls.iter_hook (line 152). But both
methods failed to update the msh. How do you update mesh in
interactive mode? Thanks!
Regarls,
Peixun
Hi,
My question is about order of approximation in field as well as integral. My understanding is that the order of field is the order of the polynomial shape functions in element and order of integral is for the integration in qp. For a 3 point triangular element as well as 4 point quadrilateral element, I think order of field should be 1 and using any other value is "wrong". Is this correct ? Also, I noticed from user manual that the order of integral = 2* order of field. I was curious what is the basis for this assumption ?
Best,
Ali
Hi
With the same element size, Poisson's ratio and elastic modulus, why is the arrangement order of the element stiffness matrix obtained in ANSYS different from that in sfepy? Can you tell me which finite element book you are referring to? Secondly, I want to set different materials for different finite element elements. What functions do I need to call?
Best，
Lee
Dear all,
I notice something in my simulations and I wanted to share it with you. I appreciate if I could have your thoughts on it: In one simulation, when I uniaxially strectch a plate and change the mesh size, I noticed that the nodal stresses and strains remain in the same range. However, once I added an elliptical hole in its center, it seems that as I decrease the mesh size, stress and strain especially at the sharp edges of the ellipse goes up and do not reach to a plateau. My understanding is that it shouldn't be real and in fact should reach to a steady state. Please let me know what you think.
Best,
Ali
I am pleased to announce the release of SfePy 2021.3.
Description
-----------
SfePy (simple finite elements in Python) is a software for solving systems of
coupled partial differential equations by finite element methods. It is
distributed under the new BSD license.
Home page: https://sfepy.org
Mailing list: https://mail.python.org/mm3/mailman3/lists/sfepy.python.org/
Git (source) repository, issue tracker: https://github.com/sfepy/sfepy
Highlights of this release
--------------------------
- unified volume and surface integration terms
- improved pyvista-based visualization script
For full release notes see [1].
Cheers,
Robert Cimrman
[1] http://docs.sfepy.org/doc/release_notes.html#id1
---
Contributors to this release in alphabetical order:
Robert Cimrman
Hugo Levy
Vladimir Lukes
Hello,
I am using SfePy for a project in aerothermal coupling. Congratulations for
the amazing job done on SfePy and thank you for making it open source.
I have a question regarding high order.
You have a parameter integral,
integrals = {
'i' : 2,
}
That we can set to the order we want.
If I am correct, it represents the order of the polynomial shape functions
and is used for the interpolation of the unknown variable in an element.
I can't find anything about the shape functions, nor in the code, nor in
the documentation.
Could you help me and point out where it is coded ?
Also, how high could SfePy support to go in order ?
I tried to put different values to this parameter but it didn't seem to
affect my solution in temperature.
I thank you in advance if you can help me.
I would like to talk to my supervisor about this during the week so I would
be really grateful if you could help me as soon as possible.
Sincerely,
Thomas Sharifi,
Sorbonne University
Dear all,
I'm observing something suspicious in a simulation that I wanted to share with you. I appreciate if you could help me out with it please: I'm simulating a 2D plate which is fixed in one end and I'm displacing the other end mimicking a uniaxial test. Now, when I average the nodal stress and divide by strain, I got a Young's modulus that shows 5% difference comparing with initial input E. So, is there any plausible explanation for it or there is something I'm missing in my simulation ? I tried different types of mesh, mesh size, mesh refinement, different integration order but all ended up into same observation.
Best,
Ali
Hi everyone,
Just a quick question on .mesh file generated with GMSH. I notice that when I generate .mesh file ( for instance with [1]) it contains EDGES section and thus if I use the generated .mesh file in sfepy if gives me error while if I manually remove that section, everything goes well. So, I was curious is there any trick so that I can avoid writing EDGES section in .mesh file or is there a shortcut in sfepy for it ? I need to avoid manually removing it.
Best,
Ali
[1] https://gitlab.onelab.info/gmsh/gmsh/blob/gmsh_4_8_4/tutorial/python/t1.py
