Re: Torque
by Robert Cimrman
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.
1 year, 1 month
1D line elements in SfePy
by Nimish
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
3 years, 2 months
ANN: SfePy 2014.4
by Robert Cimrman
I am pleased to announce release 2014.4 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 (preliminary support). It is distributed under the new
BSD license.
Home page: http://sfepy.org
Mailing list: http://groups.google.com/group/sfepy-devel
Git (source) repository, issue tracker, wiki: http://github.com/sfepy
Highlights of this release
--------------------------
- preliminary support for 1D problems
- data probes using pyVTK library
For full release notes see http://docs.sfepy.org/doc/release_notes.html#id1
(rather long and technical).
Best regards,
Robert Cimrman and Contributors (*)
(*) Contributors to this release (alphabetical order):
Lubos Kejzlar, Vladimir Lukes
4 years, 4 months
"MemoryError" when using large mesh and how does integral order change prestress shape?
by Ronghai Wu
Dear sfepy users and developers,
I am using sfepy for solving an 2D inclusion elastic problem. So, I
implemented based on "linear_elasticity/prestress_fibres" example in
interactive way. However, the prestress in my problem is not constant, so I
manually assign values.
I have two questions:
1. why when the integral order is 3, the array shape of prestress.datas[(
'Omega', 'i')][0]['val'] is (nx*ny,4,3,1), so when I assign by
prestress_value, I have to make it the same shape. If change integral order
to 1, the shape changes to (nx*ny,8,3,1), why?
2. When I use mesh smaller than 70*70, no problem. When 90*90 very slow.
When 120*120 "MemoryError" comes out.
Thanks a lot.
My codes and errors are as follow:
field = Field.from_args('displacement', np.float64, 'vector', omega, 1)
u = FieldVariable('u', 'unknown', field)
v = FieldVariable('v', 'test', field, primary_var_name='u')
prestress_value = np.zeros((nx*ny,3,1))
m = Material('m', lam=lambda_, mu=mu)
prestress = Material('prestress', val=prestress_value)
integral = Integral('i', order=3)
t1 = Term.new('dw_lin_elastic_iso(m.lam, m.mu, v, u )', integral, omega, m=m
, v=v, u=u)
t2 = Term.new('dw_lin_prestress(prestress.val, v )',integral, omega,
prestress=prestress, v=v)
eq = Equation('balance', t1-t2 )
eqs = Equations([eq])
ls = ScipyDirect({})
nls_status = IndexedStruct()
nls = Newton({'i_max' : 1,'eps_a' : 1e-6,'problem' : 'nonlinear'},
lin_solver=ls, status=nls_status)
pb = Problem('elasticity', equations=eqs, nls=nls, ls=ls)
pb.update_materials()
prestress_value = np.zeros((nx*ny,4,3,1))
#######################################################################
this is how I assign values to prestress, eta is a known scalar field.
#######################################################################
for i in range(nx*ny):
for j in range(4):
prestress_value[i][j] = np.array([[stress00_11], [stress00_22], [
stress00_12]]) * eta.value[i]
prestress.datas[('Omega', 'i')][0]['val'] = prestress_value
pb.update_materials()
########################################################################
following is error:
sfepy: updating materials...
sfepy: m
sfepy: prestress
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File
"/usr/local/lib/python2.7/dist-packages/spyderlib/widgets/externalshell/sitecustomize.py"
, line 580, in runfile
execfile(filename, namespace)
File "/home/ronghaiwu/temporary/test.py", line 304, in <module>
main()
File "/home/ronghaiwu/temporary/test.py", line 179, in main
pb.update_materials()
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/problem.py"
, line 522, in update_materials
problem=self, verbose=verbose)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/equations.py"
, line 313, in time_update_materials
verbose=verbose)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py"
, line 70, in time_update
mat.time_update(ts, equations, mode=mode, problem=problem)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py"
, line 339, in time_update
self.update_data(key, ts, equations, term, problem=problem)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py"
, line 251, in update_data
**self.extra_args)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/functions.py"
, line 34, in __call__
return self.function(*args, **_kwargs)
File
"/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/functions.py"
, line 66, in get_constants
out[key] = nm.tile(val, (coors.shape[0], 1, 1))
File "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py",
line 860, in tile
c = c.reshape(-1, n).repeat(nrep, 0)
MemoryError
Regards
Ronghai
4 years, 4 months
Problem with a hexahedral mesh
by patricia...＠upf.edu
Hi!
I am trying to simulate a hyperelastic 3D cylinder with hexahedral
elements. I did the mesh with gmsh and save it as a .mesh file. When I do
it with tetrahedral elments everything works perfect, but when I try to
simulate the mesh with hexahredal elemnts, SfePy can not read well the mesh
and give me this error:
sfepy: warning: bad element orientation, trying to correct...
sfepy: warning: bad element orientation, trying to correct...
RuntimeError: elements cannot be oriented! (0, 3_4)
I've attached the .mesh file with hexahedral elements
Thanks a lot for the help!
4 years, 4 months
Installation test failed
by Serge-Étienne Parent
Hi,
I tried to install SfePy on my openSUSE 13.2 64bit machine. I had a hard
time installing PyTables, but I could handle this by installing it using
/usr/bin/pip2.7 install tables
...after installing python-pip, hdf5-devel, python-devel and
python-numpy-devel.
Everything else went smooth and SfePy compiled successfully. But most tests
failed and importing in python did not work.
*Importing in Python:*
Python 2.7.8 (default, Sep 30 2014, 15:34:38) [GCC] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> from sfepy.base.base import *
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ImportError: No module named sfepy.base.base
*Test:*
useraccount@linux:~/Programmes/sfepy> ./run_tests.py
<<< directory: tests, test files: 82
tests/test_assembling.py
--- test_assemble_matrix: failed!
--- test_assemble_matrix_complex: failed!
--- test_assemble_vector: failed!
--- test_assemble_vector_complex: failed!
tests/test_base.py
+++ test_resolve_deps: ok
+++ test_struct_add: ok
+++ test_struct_i_add: ok
+++ test_verbose_output: ok
tests/test_cmesh.py
--- test_cmesh_counts: failed!
tests/test_domain.py
--- test instance creation failed
tests/test_ebcs.py
--- test instance creation failed
tests/test_elasticity_small_strain.py
--- test_get_solution: failed!
--- test_linear_terms: failed!
tests/test_fem.py
--- test instance creation failed
tests/test_functions.py
--- test instance creation failed
tests/test_high_level.py
--- test instance creation failed
tests/test_homogenization_perfusion.py
--- test_solution: failed!
tests/test_hyperelastic_tlul.py
--- test_solution: failed!
tests/test_input_acoustics.py
--- test instance creation failed
tests/test_input_acoustics3d.py
--- test instance creation failed
tests/test_input_active_fibres.py
--- test instance creation failed
tests/test_input_biot.py
--- test instance creation failed
tests/test_input_biot_npbc.py
--- test instance creation failed
tests/test_input_biot_npbc_lagrange.py
--- test instance creation failed
tests/test_input_cube.py
--- test instance creation failed
tests/test_input_elastic_contact_planes.py
--- test instance creation failed
tests/test_input_elastic_contact_sphere.py
--- test instance creation failed
tests/test_input_elastic_shifted_periodic.py
--- test instance creation failed
tests/test_input_its2D_2.py
--- test instance creation failed
tests/test_input_laplace_coupling_lcbcs.py
--- test instance creation failed
tests/test_input_laplace_time_ebcs.py
--- test instance creation failed
tests/test_input_linear_elastic.py
--- test instance creation failed
tests/test_input_linear_elastic_damping.py
--- test instance creation failed
tests/test_input_linear_elastic_mM.py
--- test instance creation failed
tests/test_input_linear_elastic_probes.py
--- test instance creation failed
tests/test_input_linear_elastic_tractions.py
--- test instance creation failed
tests/test_input_linear_elastic_up.py
--- test instance creation failed
tests/test_input_linear_viscoelastic.py
--- test instance creation failed
tests/test_input_material_nonlinearity.py
--- test instance creation failed
tests/test_input_navier_stokes.py
--- test instance creation failed
tests/test_input_navier_stokes2d.py
--- test instance creation failed
tests/test_input_perfusion_tl.py
--- test instance creation failed
tests/test_input_piezo.py
--- test instance creation failed
tests/test_input_poisson.py
--- test instance creation failed
tests/test_input_poisson_field_dependent_material.py
--- test instance creation failed
tests/test_input_poisson_functions.py
--- test instance creation failed
tests/test_input_poisson_neumann.py
--- test instance creation failed
tests/test_input_poisson_periodic_boundary_condition.py
--- test instance creation failed
tests/test_input_prestress_fibres.py
--- test instance creation failed
tests/test_input_sinbc.py
--- test instance creation failed
tests/test_input_stabilized_navier_stokes.py
--- test instance creation failed
tests/test_input_stokes.py
--- test instance creation failed
tests/test_input_stokes_slip_bc.py
--- test instance creation failed
tests/test_input_thermo_elasticity.py
--- test instance creation failed
tests/test_input_thermo_elasticity_ess.py
--- test instance creation failed
tests/test_input_time_poisson.py
--- test instance creation failed
tests/test_input_vibro_acoustic3d.py
--- test instance creation failed
tests/test_io.py
+++ test_recursive_dict_hdf5: ok
+++ test_sparse_matrix_hdf5: ok
tests/test_laplace_unit_disk.py
--- test instance creation failed
tests/test_laplace_unit_square.py
--- test instance creation failed
tests/test_lcbc_2d.py
--- test_linear_rigid_body_bc: failed!
tests/test_lcbc_3d.py
--- test_linear_rigid_body_bc: failed!
tests/test_lcbcs.py
--- test_laplace_shifted_periodic: failed!
tests/test_linalg.py
--- test_assemble1d: failed!
--- test_geometry: failed!
--- test_tensors: failed!
--- test_unique_rows: failed!
tests/test_linear_solvers.py
--- test instance creation failed
tests/test_linearization.py
--- test_linearization: failed!
tests/test_matcoefs.py
+++ test_conversion_functions: ok
+++ test_elastic_constants: ok
+++ test_stiffness_tensors: ok
tests/test_mesh_generators.py
--- test_gen_block_mesh: failed!
--- test_gen_cylinder_mesh: failed!
--- test_gen_extended_block_mesh: failed!
tests/test_mesh_interp.py
--- test_interpolation: failed!
--- test_interpolation_two_meshes: failed!
--- test_invariance: failed!
--- test_invariance_qp: failed!
tests/test_mesh_smoothing.py
--- test_mesh_smoothing: failed!
tests/test_meshio.py
--- test instance creation failed
tests/test_msm_laplace.py
--- test instance creation failed
tests/test_msm_symbolic.py
--- test instance creation failed
tests/test_normals.py
--- test_normals: failed!
tests/test_parsing.py
--- test_parse_equations: failed!
--- test_parse_regions: failed!
tests/test_periodic_bc_2d.py
--- test instance creation failed
tests/test_periodic_bc_3d.py
--- test instance creation failed
tests/test_permutations.py
--- test_rcm: failed!
tests/test_poly_spaces.py
--- test instance creation failed
tests/test_projections.py
--- test instance creation failed
tests/test_quadratures.py
--- test_quadratures: failed!
--- test_weight_consistency: failed!
tests/test_regions.py
--- test instance creation failed
tests/test_semismooth_newton.py
--- test_semismooth_newton: failed!
tests/test_sparse.py
--- test_compose_sparse: failed!
tests/test_tensors.py
--- test_stress_transform: failed!
--- test_tensors: failed!
--- test_transform_data: failed!
--- test_transform_data4: failed!
tests/test_term_call_modes.py
--- test instance creation failed
tests/test_term_consistency.py
--- test instance creation failed
tests/test_units.py
+++ test_consistent_sets: ok
+++ test_units: ok
tests/test_volume.py
--- test instance creation failed
82 test file(s) executed in 0.46 s, 95 failure(s) of 106 test(s)
*Test debug for test_term_consistency.py:*
useraccount@linux:~/Programmes/sfepy> python run_tests.py --debug
tests/test_term_consistency.py
<<< directory: tests, test files: 1
<<< tests/test_term_consistency.py
sfepy: left over: ['verbose', '__builtins__', '_integrate', 'test_terms',
'__doc__', '__name__', 'get_pars', 'TestCommon', 'dim', 'data_dir', 'nm',
'__package__', '_filename', '__file__']
--- test instance creation failed
Traceback (most recent call last):
File "run_tests.py", line 206, in <module>
main()
File "run_tests.py", line 197, in main
run_tests(stats, dirname, [filename])
File "run_tests.py", line 131, in run_tests
n_fail, n_total, test_time = run_test(conf_name, options)
File "run_tests.py", line 83, in run_test
test = conf.funmod.Test.from_conf(conf, options)
File
"/home/sergeetienne/Programmes/sfepy/tests/test_term_consistency.py", line
94, in from_conf
from sfepy.discrete import Problem
File "/home/sergeetienne/Programmes/sfepy/sfepy/discrete/__init__.py",
line 8, in <module>
from variables import Variables, Variable, FieldVariable,
create_adof_conns
File "/home/sergeetienne/Programmes/sfepy/sfepy/discrete/variables.py",
line 12, in <module>
import sfepy.linalg as la
File "/home/sergeetienne/Programmes/sfepy/sfepy/linalg/__init__.py", line
1, in <module>
from extmods.crcm import rcm, permute_in_place
ImportError: No module named crcm
4 years, 5 months
Re: [sfepy-devel] How to implement non-constant prestress problem?
by Robert Cimrman
Hi Ronghai,
On 11/05/2014 09:42 PM, Ronghai Wu wrote:
> Dear Sfepy developers and users,
>
> I am wondering if it is possible to solve a non-constant prestress problem
> with sfepy, as show
>
>
>
> where all fields exist in the whole domain(no subdomain is considered).
> is a constant stress field, however is a non-constant scale field(having
> different values at different nodes, but values are known). So, therefore
> prestress is also non-constant.
> I know the "linear_elasticity/prestress_fibres" example in Sfepy
> documentation. But the prestress in that example is constant(having same
> values at different nodes). So, I am wondering if it is possible to
> implement non-constant prestress problem. Thanks.
It should work out of the box. Look at [1] to see how to define a non-constant
material parameter (the prestress here) by using functions.
r.
[1] http://sfepy.org/doc-devel/users_guide.html#defining-material-parameters
4 years, 5 months
