Mailman 3 May 2017 - SfePy

Material Parameters of Individual Nodes
by derj...＠web.de April 12, 2022

April 12, 2022

Hello! It's a great tool you built there! I'd like to use it for some kind of topology optimisation. Is there an "idiomatic" way to set the material parameters of individual nodes? Thank you very much! JonnyB

3 5

Calculate total force from cauchy stress
by Patricia Garcia Cañadilla Jan. 17, 2019

Jan. 17, 2019

Dear Robert, I have a 2D body of hyperelastic material which contracts and I would like to compute the total force developed by the body from the cauchy stress. I am trying to follow some of your indications I found in this group, but I still couldn't make it works. Could you please help me to fix the problem? I am getting the following error: key = (region.name, integral.order, integration) AttributeError: 'dict' object has no attribute 'name' I am trying to do the following, inside stress_strain post-processing function: def stress_strain(out, problem, state, extend = False ): from sfepy.base.base import Struct from sfepy.mechanics.tensors import StressTransform from sfepy.mechanics.tensors import transform_data from sfepy.discrete.common.mappings import get_normals ev = problem.evaluate field = problem.fields['displacement'] region = problem.domain.regions['Gamma'] integral = problem.integrals['i2'] n = get_normals(field,integral,regions) stress = ev('dw_tl_fib_a.1.Omega(f1.fmax, f1.eps_opt, f1.s, f1.fdir, f1.act, v, u )',mode='qp', term_mode= 'stress'); F = ev('ev_def_grad.1.Omega(u)',mode='el_avg'); transform = StressTransform(F) Cstress = transform.get_cauchy_from_2pk(stress) T = Cstress*n; Force = ev('ev_surface_integrate.2.Gamma(T)') And here it is part of the problem configuration file. fields = { 'displacement': ('real', 'vector', 'Omega', 1), } materials = { 'solid' : (None, 'get_elastic_pars'), 'load' : (None, 'linear_tension'), 'f1' : 'get_pars_fibres1', } variables = { 'u': ('unknown field', 'displacement', 0), 'v': ('test field', 'displacement', 'u'), } regions = { 'Omega' : 'all', 'Fix1' : ('vertices in x < %.10f' % (fix_point + eps2), 'facet'), 'Fix2' : ('vertices in x > %.10f' % (fix_point - eps2), 'facet'), 'Fix' : ('r.Fix1 *v r.Fix2', 'facet'), 'Gamma' : ('vertices of surface','edge'), } ebcs = { 'fixb' : ('Fix', {'u.all' : 0.0}), } integrals = { 'i1' : ('v', 1), 'i2' : ('s', 2), }

4 5

Re: Torque
by Robert Cimrman March 9, 2018

March 9, 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

multiphysics: piezoelectric + solid mechanics + acoustic
by firdaus prabowo May 22, 2017

May 22, 2017

Hi, I'd like to try solving multiphysics problem with Sfepy, in particular when a piezoelectric disc attached on a solid vibrates and propagating waves from the solid part to a water domain (acoustic). Is there any example available for a similar problem? I would like to compare the results with COMSOL results (readily available). Thanks. Regards, Firdaus

2 1

question about non-linear solver
by Patricia Garcia Cañadilla May 20, 2017

May 20, 2017

Dear all, I have some questions about the internal behaviour of Sfepy that I hope you can answer me. - What is it the convergence criterion of the non-linear Newton solver? - When modelling a Nearly incompressible Mooney-Rivlin hyperelastic material, like the example shown here: http://sfepy.org/doc-devel/examples/large_deformation/hyperelastic.html, how the nearly-incompressibility constraint is stabilised? Thank you very much in advance.

2 1

ANN: SfePy 2017.2
by Robert Cimrman May 19, 2017

May 19, 2017

I am pleased to announce release 2017.2 of SfePy. Description ----------- SfePy (simple finite elements in Python) is a software for solving systems of coupled partial differential equations by the finite element method or by the isogeometric analysis (limited support). It is distributed under the new BSD license. Home page: http://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 -------------------------- - simplified and unified implementation of some homogenized coefficients - support for saving custom structured data to HDF5 files - new tutorial on preparing meshes using FreeCAD/OpenSCAD and Gmsh For full release notes see http://docs.sfepy.org/doc/release_notes.html#id1 (rather long and technical). Cheers, Robert Cimrman --- Contributors to this release in alphabetical order: Robert Cimrman Jan Heczko Lubos Kejzlar Vladimir Lukes Matyas Novak

1 0

Re: Problems with Linear Elasticity Example
by Robert Cimrman May 9, 2017

May 9, 2017

On 05/09/2017 01:36 AM, Jeshwanth.Kundem(a)wagner-group.com wrote: > You are right about scikit-umpfack being installed somewhere else. I removed > it and it worked as you mentioned in the previous email. OK. Now if you do a fresh scikit-umfpack install, does it help? > I will check the example. > > I have to disagree with the Anaconda statement. I tried, tried and tried > with Anaconda. It might be an easy install, but what is happening is the > Anaconda is at 4.3.1. When I install sfepy through conda forge channel, the > sfepy or conda forge wants the anaconda to be 4.2.13 or something. For that > reason, it gave me an option mentioning "The packages will be superseded by > a higher priority channel" and it is actually trying to get back to 4.2.13. > Because of this lot of anaconda packages are being unlinked with conda. They > are causing other problems down the line, with mayavi, wxpython, serial etc. > I would like to know what you think about my statement. I am not aware of any such anaconda version dependence in the sfepy recipe - maybe some of the dependencies need that. Anyway, I am not sure how to solve that. If you happen to find some more information, let us know. > Are you mentioning that if I use conda, I would not have the scikit-umfpack > problem and UMFPACK would work without any errors? If so let me know, I will > give it a try once again. It might. I did not try it recently, but that was the idea. If it does not work, submit an issue to the package that fails, please. r. > Because of this I intentionally deleted anaconda and installed python and > packages one by one with pip. > > > > > > -----Original Message----- From: Robert Cimrman [mailto:cimrman3@ntc.zcu.cz] > Sent: Monday, May 08, 2017 5:03 PM To: Kundem, Jeshwanth Cc: > sfepy(a)python.org Subject: Re: [sfepy] Problems with Linear Elasticity > Example > > On 8.5.2017 23:34, Jeshwanth.Kundem(a)wagner-group.com wrote: >> Hi Robert >> >> The example works when I change the ls solver from scipy_direct to >> scipy_iterative. Not when I uninstall scikit-umfpack. It is apparently not >> being defaulted to scipy LU solver. Any thoughts on this? However, I am > > If you remove scikit-umfpack and still see the same error, I would guess you > have another version of scikit-umfpack installed somewhere. > >> unable to specifiy the number of iterations if I change to >> scipy_iterative. > > See examples/diffusion/cube.py for an example of how to set the number of > iterations. > >> The UMFPACK has to be installed before using sckit-umfpack. There is no >> good documentation on that . I got onto this website: >> http://faculty.cse.tamu.edu/davis/suitesparse.html . However it is only >> available for MATLAB and C. > > How did you install python and other sfepy dependencies? (Which versions?) > The easiest way nowadays is to use anaconda - there is even a recipe for > sfepy on conda-forge ("conda install -c conda-forge sfepy" should work). > > r. >

1 0

Re: Problems with Linear Elasticity Example
by Robert Cimrman May 8, 2017

May 8, 2017

On 8.5.2017 23:34, Jeshwanth.Kundem(a)wagner-group.com wrote: > Hi Robert > > The example works when I change the ls solver from scipy_direct to > scipy_iterative. Not when I uninstall scikit-umfpack. It is apparently not > being defaulted to scipy LU solver. Any thoughts on this? However, I am If you remove scikit-umfpack and still see the same error, I would guess you have another version of scikit-umfpack installed somewhere. > unable to specifiy the number of iterations if I change to scipy_iterative. See examples/diffusion/cube.py for an example of how to set the number of iterations. > The UMFPACK has to be installed before using sckit-umfpack. There is no good > documentation on that . I got onto this website: > http://faculty.cse.tamu.edu/davis/suitesparse.html . However it is only > available for MATLAB and C. How did you install python and other sfepy dependencies? (Which versions?) The easiest way nowadays is to use anaconda - there is even a recipe for sfepy on conda-forge ("conda install -c conda-forge sfepy" should work). r.

1 0

Problems with Linear Elasticity Example
by Jeshwanth.Kundem＠wagner-group.com May 8, 2017

May 8, 2017

Hi I am trying to work out the example located on the following webpage http://sfepy.org/doc-devel/primer.html#sec-primer I am only trying to do the first piece of code, I am getting the following error and I am completely unsure what the problem is. I tried uninstalling and installing scikit-umfpack like bunch of times. Still do not know what the problem is. Any help would be appreciated on how to install umfpack. Other examples are running fine, except this one and elasticity is the only important problem for me. [cid:image001.png@01D2C813.8ECF9BB0] Regards Jeshwanth

2 1