Questions tagged [finite-element-analysis]

Anything related to finite element analysis or finite element methods (FEMs), i.e. a number of advanced numerical analysis techniques for finding approximate solutions of differential equations. FEMs are heavily used in many branches of engineering, for example to simplify the design phase of a project.

The methods transform the differential equation into a weak form over the problem domain, which is then discretized by a set of elements (e.g. a triangle-shaped plate or quadrilateral shells). This allows for a per-element approximation of the solution and an assembly of global quantities by local quadrature rules, resulting in a possibly nonlinear system of equations. For large problem statements, finding a solution might require parallel computing strategies.

Many areas of engineering profit from these techniques, e.g. solid and structural mechanics on various time- and size-scales, electrical field theory or diffusion problems, and there are plenty of (free) software solutions available. In scientific computing, users regularly need to extend the method by custom element formulations or new solution procedures, which is often done in fortran, c++ or python.

See also the Wikipedia page on FEMs.

235 questions
2
votes
2 answers

re-grid/re-mesh unstructured 3D array data into a pre-existing grid

I have a map of values for a specific grid: X Y Z Value 0.555 0.334 0.472 1.361 0.674 0.729 0.711 2.114 0.538 0.982 0.886 2.406 0.775 0.279 0.789 1.843 0.292 0.206 0.984 1.482 I would like to convert…
umn
  • 431
  • 6
  • 17
2
votes
1 answer

parallelizing python finite element analysis

I am trying to parallelize the python function below, which is involved in finite element analysis. In particular, I am trying to make the for loop in the function runs in parallel. I have done it in Matlab using parfor, and I am trying to do the…
2
votes
0 answers

How to import non-Newtonian model in FEATOOL Multiphysics?

I am using FEATool Multiphysics. I have a MATLAB function named ex_nonnewtonian1.m, which is supported by FEATool itself and contains the definition of a power-law non-Newtonial model. This file wasn't added to FEATOOL, so I need to add it…
2
votes
2 answers

Create triangular meshing of a non rectangular shape in 2D in Python

Suppose I have a set of points defining a perimeter of a non rectangular shape in the 2d plane. I need a function to create a triangular meshing where I can modify the number of triangle cells and return (x,y) coordinates of each cell. Thank you.
2
votes
2 answers

Solving a large (n>1000000) linear system of equations

I have reached my limit on the following problem: As part of my FEA code (in MATLAB) I need to find x, x=A\b. Both A and b are sparse, complex, double precision matrix and vector respectively. The size of A is (n,n) and b is (n,1) where n is 850000…
user5489
  • 123
  • 2
  • 8
2
votes
2 answers

Add internal boundary or crack in PyGmsh / Gmsh

I am trying to generate a finite element mesh using PyGmsh, using the following code: import pygmsh geom = pygmsh.opencascade.Geometry( characteristic_length_min=0.1, characteristic_length_max=0.1, ) rectangle =…
Carlos
  • 586
  • 1
  • 9
  • 19
2
votes
2 answers

How to merge values from dictionary on different keys while iterating through it ? Finite element mesh merge algorithm

I am working on my phd and I am stuck on this step. The problem consists of implementing a finite element mesh merging algorithm and maybe my solution is not the best, so if you think of a better one I am open to suggestions. Regarding the problem:…
2
votes
0 answers

How to work around GMSH and Calculix GUI in FreeCAD to automate Finite Element Analysis using Macro?

I created a macro (python script) in FreeCAD that generates a cantilever beam, creates an Analysis, selects a material, and applies constraints to the body. How do I now extend the script, so I can generate the mesh (using GMSH) and run the analysis…
avgJoe
  • 832
  • 7
  • 24
2
votes
2 answers

Matlab: associating an ID with a dataset (e.g. struct)?

I am developing a certain feature for a high-order finite element simulation algorithm in Matlab and I am wondering what is a good way of implementing a certain task. I believe I am facing a somewhat common problem, but after doing some digging, I'm…
2
votes
2 answers

efficiently insert diagonal in sparse matrix in julia

I'm studying EM computational methods on this resource. These methods use a lot of big sparse matrices with only few diagonals set to non-zero. So my question is this: how do I efficiently set a diagonal of the existing matrix in place in julia?
Sergey
  • 1,166
  • 14
  • 27
2
votes
1 answer

Abaqus C++ API syntax

I am new to Abaqus C++ API and have question about the syntax. db_Odb& odb = openOdb("file.odb"); odb_Step& step = odb.steps()["Step-1"]; odb_Instance& instance = odb.rootAssembly().instances()["PART-1-1"]; There are several instances that…
dfrw4r3
  • 25
  • 4
2
votes
2 answers

Solving a linear system of equations

I am working on using the finite element method to calculate heat flow through elements. I'm currently stuck on solving a system of equations where there is variables on both side of the equality. A simple example could be something like this | 1 …
2
votes
1 answer

Riding the wave Numerical schemes for hyperbolic PDEs, lorena barba lessons, assistance needed

I am a beginner python user who is trying to get a feel for computer science, I've been learning how to use it by studying concepts/subjects I'm already familiar with, such as Computation Fluid Mechanics & Finite Element Analysis. I got my degree in…
2
votes
1 answer

Boundary markers of nodes and faces in Tetgen

Let's say I create a .poly detailing the information of the my domain geometry. I also add node boundary markers and face boundary markers. Question 1: Do I have to be consistent between the boundary markers of nodes and faces? I mean, if a node…
2
votes
0 answers

Numerical integration over a mesh in Python with existing mesh

I found this discussion in my previous search which shows how to use Python for numerical integration of a function on points in space. Looking for Python package for numerical integration over a tessellated domain My first question is will the…
1 2
3
15 16