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),
}
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.
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
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.
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
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.
>
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.
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
