Mailman 3 August 2013 - SfePy

Re: Torque
by Robert Cimrman 09 Mar '18

09 Mar '18

In examples/large_deformation/hyperelastic.py a rotation by displacements is applied. By using a similar function the vectors defining the force couples could be defined for dw_surface_ltr (IMHO). Does it make sense? r. ----- Reply message ----- From: "Andre Smit" <freev...(a)gmail.com> To: <sfepy...(a)googlegroups.com> Subject: Torque Date: Sat, Dec 18, 2010 05:10 What is the best way to apply a torque load to a model? -- Andre -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

4 7

1D line elements in SfePy
by Nimish 22 Jan '16

22 Jan '16

I am currrently looking for FEM packages to help me solve a system of beams and columns, basically a collection of 1D bernoulli/timoshenko line elements. I started reading SfePy docs and i am getting the idea that doing the above is not really possible here, am i right? Are only 2D area elements permitted in SfePy? Or is there any direct support for solving 1D line elements too.. Cheers Nimish

4 9

Python 3 compatibility
by Robert Cimrman 11 Aug '14

11 Aug '14

FYI: As SciPy 0.12.0 is out and one of the release highlights is "Support for Python 2 and Python 3 from the same code base (no more 2to3)", we can think seriously about updating SfePy in this respect as well, cf. [1]. r. [1] https://github.com/sfepy/sfepy/issues/164

4 5

Evaluate a solution in the arbitrary point in the domain
by Alec Kalinin 03 Apr '14

03 Apr '14

Dear SfePy users, Is it possible to evaluate a solution not only in the FEM mesh node, but in any arbitrary point in the domain with the given (x, y, z) coordinates? For example, consider Dirichlet problem for Poisson equation. We apply essential boundary conditions on the surface nodes and after the problem has been solved we have the solution vector, i.e. vector of values in the FEM mesh nodes. But I want to know the solution in point v(x, y, z) that is not FEM mesh node. What is the best way to obtain solution in this point v? Sincerely, Alec Kalinin

3 25

elliptical hole under isotropic stress
by David N. Mashburn 13 Jan '14

13 Jan '14

Hello sfepy developers and users! I am modelling a simple linear elastic sheet under isotropic stress with an elliptical hole in the center (and I have it working under sfepy, great little platform!). It is obvious the model should initially yield more easily in the direction of the short axis of the ellipse. What is not so obvious to me is what should happen in the limit as stress goes to infinity. Part of me wants to believe that the hole should eventually become a circular, but the results of the simulation show that the ellipse eventually switches its aspect ratio with what was the the short axis becoming the long axis and vice-versa. My question is whether: A: The finite element result is the product of a small-displacement/non-moving mesh artifact (and if so, if there is a way to get the correct behavior using sfepy...) OR B: My intuition about the physical behavior of this ideal system is incorrect and the ellipse really wouldn't round out into a circle under increasingly large stress (aka, the FE model is still physical/correct with large displacements). This might be obvious to people who have done more finite element modeling than I have, but thanks anyway! I'm attaching a picture to make it easier to see at a glance (quarter-ellipse with x and y symmetry boundary conditions and equal tractions applied at the top and right boundaries). Thanks! -David Mashburn

4 12

post process energy calculation
by Anastasia Konovalova 15 Nov '13

15 Nov '13

Hello, sfepy group, I have a question about post-process calculations. For my problem I used the modified examples linear_elastic.py and hyperelastic.py. After the calculation of strain and stress tensors I need to get energy of deformation \int_{V} S : C, where C = 2*E+I. I think I need to calculate the multiplication of tensors first, as some material variable, and then integrate it. So I have 2 questions: 1. How can I calculate the double scalar multiplication of two arrays strain and stress that have shape (number of elements, 1, 6, 1)? (Double scalar multiplication of A and B is a scalar that can be got as A_{ij}B_{ji}). Because of the shape of arrays I can't use numpy.inner(stress,strain). I made some function for these calculation, but perhaps, the simplier way exsists. 2. Then I get the multiplication - the scalar function of energy I should integrate it. For this I get the dw_volume_integrate term. I'm trying to integrate material parameter, where energy has been saved. solid.update({ 'E' : get_W(stress,strain),#energy in elements }) U = problem.evaluate('dw_volume_integrate.i1.Omega(solid.E, v)',) Or may be I should use some new variable and initialize it in post-proc. I'm asking for help in correct formulation of this integral. Should I define energy as material parameter or some new variable? And how I can do it. Thanks in advance,

2 15

GSoC Phase II
by Ankit Mahato 19 Sep '13

19 Sep '13

Hi R, I tried mem_checkIntegrity inside the functions but it is not printing any debugging result. mem_checkIntegrity(421,"convect_build_vtg","termsNavierStokes.c","~/sfepy/sfepy/terms/extmods"); mem_checkIntegrity(546,"term_ns_asm_convect","termsNavierStokes.c","~/sfepy/sfepy/terms/extmods"); Am I adding it correctly? regards.

2 25

Re: [sfepy-devel] GSoC Phase II
by Robert Cimrman 12 Aug '13

12 Aug '13

Probably yes but not today. Try putting the check statement to the beginning of the function and various other places. Also try printf/fmf_print, and print sizes of variable, to help determining where the problem actually is. r. ----- Reply message ----- From: "Ankit Mahato" <aryan...(a)gmail.com> To: <sfepy...(a)googlegroups.com> Subject: [sfepy-devel] GSoC Phase II Date: Mon, Aug 12, 2013 15:28 Hi R, Is it possible for you to come on Teamviewer - http://www.teamviewer.com/hi/index.aspx Regards. On Monday, 12 August 2013 18:52:29 UTC+5:30, Ankit Mahato wrote: On Sunday, 11 August 2013 23:21:13 UTC+5:30, Robert Cimrman wrote:Hi Ankit, On Sat, 10 Aug 2013, Ankit Mahato wrote: > Hi R, > > I tried mem_checkIntegrity inside the functions but it is not printing > any debugging result.ï¿½ > mem_checkIntegrity(421,"convect_build_vtg","termsNavierStokes.c","~/sfepy/sfepy/ > terms/extmods"); > mem_checkIntegrity(546,"term_ns_asm_convect","termsNavierStokes.c","~/sfepy/sfep > y/terms/extmods"); > > Am I adding it correctly? Use the macro without arguments (see sfepy/fem/extmods/common.h for definition): check_memory_integrity(); and do not forget to rebuild the sources. After adding it to convect_build_vtg() ... default: errput( ErrHead "ERR_Switch\n" ); return( RET_Fail ); } check_memory_integrity(); return( RET_OK ); } ... and typing "make", I get tons of checking memory integrity in sfepy/terms/extmods, sfepy/terms/extmods/termsNavierStokes.c, convect_build_vtg(), 476: allocated memory: 8 records, usage: 52032, max: 52032 memory OK. - trying with the examples/navier_stokes/navier_stokes2d.py Also try changing the debug flags in site_cfg.py: debug_flags = '-DDEBUG_FMF' or even debug_flags = '-DDEBUG_FMF -DDEBUG_MESH' It printed the debug results for examples/navier_stokes/navier_stokes2d.py Then I tried running the convective_difffusive.py and it did not print any line containing allocated memory report: .. .. .. convdiff: matrix structural nonzeros: 1747752 (8.65e-04% fill) convdiff: updating materials... convdiff: m convdiff: ...done in 0.01 s convdiff: nls: iter: 0, residual: 2.002082e-02 (rel: 1.000000e+00) convect_build_vtg(): ERR_Switch Traceback (most recent call last): File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 131, in <module> main() File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 113, in main flow = problem.solve() File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/problemDef.py", line 933, in solve vec = solvers.nls(vec0) File "/usr/local/lib/python2.7/dist-packages/sfepy/solvers/nls.py", line 345, in __call__ mtx_a = fun_grad(vec_x) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/evaluate.py", line 66, in eval_tangent_matrix mtx = pb.equations.eval_tangent_matrices(vec, mtx) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 640, in eval_tangent_matrices self.evaluate(mode='weak', dw_mode='matrix', asm_obj=tangent_matrix) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 526, in evaluate asm_obj=asm_obj) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 766, in evaluate ret_status=True) File "/usr/local/lib/python2.7/dist-packages/sfepy/terms/terms.py", line 1474, in evaluate vals.append(self.sign * val) RuntimeError: ccore error (see above) r. -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sfepy-devel...(a)googlegroups.com. To post to this group, send email to sfepy...(a)googlegroups.com. Visit this group at http://groups.google.com/group/sfepy-devel.

2 1

# triangular faces have -1 as 4. point.
by Hendrik Mohrmann 02 Aug '13

02 Aug '13

Hello. I got the following error message, when I run sfepy with my mesh File "/home/rasputin/Uni/sfepy/sfepy/fem/domain.py", line 73, in _region_leaf if aux[0] == -1: # Triangular faces have -1 as 4. point. IndexError: index 0 is out of bounds for axis 0 with size 0 How would I trace the error? Is it an error from the mesh, so do I have to change the mesh in some way? I checked the .mesh file and there are at least no negative values as suggested by the comment. I have very basic knowledge about the stuff I am trying to do and I constantly fear that the answers to my questions are very trivial. I just modified one of the given navier-stokes examples and try to run it with my own geometry of a micromixing device. I created the mesh with gmsh and tried the optimize 3D and optimize 3D (Netgen) (without really knowing what they do), but that did not change the error. Your help would be very much appreciated. Best, Hendrik

2 11