For much more rotations, a powerful numeric continuation method is your friend.
The "speed" of the rotation is a result of the adaptive step-size of contique (one frame per converged step, constant framerate).
https://github.com/adtzlr/contique
#computationalmechanics #scientificcomputing #FiniteElementMethod #fem #fea #python #opensource
🚀🚀 FElupe is now available on conda-forge! 🚀🚀 Unfortunately, it doesn't handle optional dependencies like pip. Anyway here's a nice animation of the #hyperelastic "twisting column" example from the book of Bonet & Wood. This one works without any fancy continuation techniques, it's just Newton's method.
https://github.com/adtzlr/felupe
#computationalmechanics #FiniteElementMethod #scientificcomputing #fem #fea #python #opensource
Sitting here running #FEA simulations in my hotel room, and my model has 1.24 MILLION degrees of freedom. I’m sure any of you that do FEA probably think that’s pretty pedestrian, but I don’t think many people would realize the size of the matrices being pushed around. This simulation takes about 25 minutes to run because of large displacement and surface contact, and it’s rather amazing that it only takes 25 minutes on modern hardware.
The result of the getting started notebook of #felupe, which is deployed by #jupyterlite on #github, is now ready-to-visualize directly in the browser, thanks to #jupyterview!
https://adtzlr.github.io/felupe-web/lab?path=01_Getting-Started.ipynb
#python #webassembly #scientificcomputing #computationalmechanics #finiteelementmethod
Merry Christmas everyone (...if you're living in a culture where this is celebrated). 🎄🎄 ☃️ ❄️ ⛷️ It's time to take a break from business as usual - get your batteries reloaded! 🔋
(This is a hyperelastic Mooney-Rivlin christmas tree under varying body forces.)
#computationalmechanics #scientificcomputing #python #finiteelementmethod #holiday #simulation #hyperelasticity
Hello, world! Wait, for a nonlinear simulation #python package? If you'd like to use FElupe for your #hyperelastic #simulation but don't know where to start & too lazy to read the docs, you may use the hello_world() function to print a minimal-working example to the console.
Initially, I wasn't sure if that kind of function really makes sense. But it turns out that I'm using this one quite frequently. 📝 (shhh! Author here... 🤓)
https://github.com/adtzlr/felupe
#scientificcomputing #computationalmechanics #finiteelementmethod #fem #fea
Here's an example: let's take the Neo-Hookean [1] isotropic #hyperelastic material model formulation. Just define its strain energy density function and use it in your #simulation [2].
[1] https://en.wikipedia.org/wiki/Neo-Hookean_solid
[2] https://github.com/adtzlr/felupe
#python #scientificcomputing #finiteelementmethod #computationalmechanics #fea #fem #jax
Have you ever thought 💡 of using JAX as 🧮 #automaticdifferentiation engine in 💻 finite element simulations? Boost the performance 🏇 of computationally-expensive hyperelastic material models with #jit in 🔍 FElupe! 🚀 🚀
https://github.com/adtzlr/felupe
#python #jax #finiteelementmethod #scientificcomputing #computationalmechanics #fea #fem #hyperelasticity
With 🔍 FElupe, it's easier than ever to simulate the deformation of nonlinear hyperelastic solid bodies!
https://github.com/adtzlr/felupe
#python #numpy #scipy #fea #finiteelementmethod #numericsimulation #scientificcomputing #opensource #computationalmechanics #hyperelasticity
In my Python FEA package FElupe I've noticed that the time spent on assembling (especially large) finite element matrices takes longer than in recent releases. The slowdown is about 15% to 100%, depending on the number of degrees of freedom. From a computational-cost point of view, there are some cheap arrays stored inside a numeric region. A somewhat massive einsum-call further evaluates the values for the sparse finite element matrix. Due to some recent code changes, these cheap region arrays aren't contiguous anymore. That slows down the finite element matrix assembly up to 100%! By ensuring C-contiguous arrays, performance is back to normal (NumPy array flags are available as attribute ndarray.flags).
https://github.com/adtzlr/felupe
#python #numpy #FiniteElementMethod #finiteelements #fem #fea #scientificcomputing #computationalmechanics
FElupe 🔍 now has visual previews for the schemes of most numeric regions (finite element formulation + quadrature rule) in its HTML API docs - thanks to PyVista. There are also a lot of helpful static and interactive plots in other sections of the docs too. Go check it out!
Simulating the nonlinear deformation of hyperelastic structures has never been easier in a few lines of Python code.
https://felupe.readthedocs.io/
#fea #fem #scientificcomputing #computationalmechanics #FiniteElementMethod #finiteelements #python #plot
Interesting paper by @jedbrown et al.
https://doi.org/10.48550/arXiv.2401.13196
For computational mechanics/physics, if you code by just punching in the equations from the textbooks directly, the physics should work, but computationally the way you evaluate the quantities may be unstable. This paper lists some recipes to avoid these.
Mostly small strain problem, but still feels icky to leave in.
![Figure 2 from the paper with the caption: "Relative error of standard computation of J–1 and its stable [equivalent]"](https://files.mastodon.social/cache/media_attachments/files/111/832/626/532/693/133/original/555b79129213ff4e.png)
Just finished the "beta" of my javascript sparse matrix solver. For testing, I first created a 20000 x 20000 sparse test matrix A and the vector c. Solving the equation Ab=c took 20 seconds. Quite feasible for a javascript code I guess. With the solver I want to carry out some finite element analysis in a web browser.
#StructuralEngineering #FiniteElementMethod #NumericalMethods #javascript
The docs of the Tensor Toolbox for Modern Fortran have finally arrived in the 21st century. New look, search included. Useful for constitutive material formulations e.g. in #FEA software packages like #Abaqus #umat or #Marc #hypela2 subroutines.
#Fortran #FEM #hyperelasticity #scientificcomputing #computationalmechanics #FiniteElementMethod #opensource #openscience
Last week Andrei gave his first #keynote #plenary talk at the #ECCOMAS Young Investigators Conference 2023 in the Faculdade de Engenharia da Universidade do Porto. Many thanks to the organisers, particularly Rodrigo Pinto Carvalho and Igor Lopes, for inviting. Advanced #opensource #finiteelementmethod tools are becoming increasingly popular in the ECCOMAS community!
