Mailman 3 March 2014 - SfePy

Re: Torque
by Robert Cimrman March 10, 2018

March 10, 2018

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 Jan. 22, 2016

Jan. 22, 2016

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 Aug. 11, 2014

Aug. 11, 2014

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 April 3, 2014

April 3, 2014

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

thermal diffusion Dirichlet vs. Neumann BCs
by seismo_phil March 26, 2014

March 26, 2014

Hi, Thanks for inviting me to contribute usage issues, Robert. I'm using sfepy to solve a relatively simple, 2D, thermal diffusion problem. i can solve it useing the following ebcs and equations (i.e., implied zero heat flux at left and right boundaries): ebcs = { 't1' : ('Gamma_Top', {'t.0' : 0.0}), 't2' : ('Gamma_Bottom', {'t.0' : 1000.0}), } equations = { 'Temperature' : """dw_laplace.2.Omega( coef.val, s, t ) = 0 """ } This works perfectly. I actually need to solve the same problem with an ebc at the top but a Neumann boundary condition at the bottom. So I use: ebcs = { 't1' : ('Gamma_Top', {'t.0' : 0.0}), } equations = { 'Temperature' : """dw_laplace.2.Omega( coef.val, s, t ) = dw_surface_integrate.2.Gamma_Bottom(flux.val, s)""" } where flux is defined in materials as: materials = { 'flux' : ({'val' : 1.0},), 'coef' : ({'val' : 1.0},), } This seems to be exactly like what is done in the example poisson_neumann.py But now, I find that the solver fails: sfepy: updating materials... sfepy: coef sfepy: flux sfepy: ...done in 0.01 s sfepy: nls: iter: 0, residual: 9.332878e+04 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.00 [s] sfepy: matrix: 0.00 [s] sfepy: linear system not solved! (err = 5.659775e-10 < 1.000000e-10) sfepy: nls: iter: 1, residual: 5.748516e-10 (rel: 6.159425e-15) I'm confused by this, because I thought that this is still a well-posed problem. Can anyone please suggest to me what I may be doing wrong? Regards, - Phil

2 5

usage issues
by seismo_phil March 24, 2014

March 24, 2014

Hi, I have just started with sfepy and have a question about using it. I have browsed the tutorial and primer but haven't found what I need. The question is about the use of dw_volume_integrate. I will be happy to go into specifics but first wanted to check that it's OK to send usage questions to this list? Thanks, - Phil

2 1