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
0
votes
1 answer

Calculating thickness of 3D geometry of a model

I am doing finite element modelling. I have a geometry which does not have a uniform thickness and I want to calculate the thickness of that in 4 distinct nodes and compare it with another model. Is there any method to calculate the thickness of the…
0
votes
1 answer

Generate unstructured surface mesh from coordinates

Problem summary I have a rectangular prism like in the image below where I know the coordinates of the vertices. I want to discretize the lateral faces of the prism into an unstructured surface mesh using quadrilateral elements with a certain…
fma
  • 219
  • 2
  • 14
0
votes
0 answers

Earthquake simulation - finite elements method (Matlab)

do you know is there's a way to simulate a earthquake with the finite elements method using the Partial Differential Equation Toolbox on Matlab ? My goal is to visualize the amplitude of a (Rayleigh) wave in all points of a forest for different…
Steph
  • 1
0
votes
0 answers

How to avoid accumulated error over time and "drift" in numerical time integration?

I'm coding a few programs to study structural dynamics problems, and I'm trying to understand how to estimate and minimize accumulated error over time for large-time scale calculations. An explicit time integrator uses the knowledge of the current…
0
votes
2 answers

Is it possible to write a 1D subroutine in ABAQUS?

I'm trying to write a subroutine for ABAQUS in fortran and during my research on how to write subroutines I only encountered subroutines for 3D models. So I was wandering if it's possible to make 1D subroutines.
0
votes
0 answers

Setting Boundary Conditions in Finite Element

EDIT: See solution in the full code, following @bgb2's comment. I'm currently trying to code a Finite Element Analysis to solve a 2D heat conduction problem. For now I'm looking at the steady state system: where k is the thermal conductivity of the…
Sorade
  • 915
  • 1
  • 14
  • 29
0
votes
0 answers

Is it possible to perform a geometric non-linear analysis in Abaqus with multiple independent steps?

I am trying to analyse a large structure with multiple loading scenarios (independent - different loading directions), while considering geometric non-linearity. Using different steps is not sufficient, since each step starts off with the geometry…
Aravind
  • 3
  • 3
0
votes
0 answers

Application of the formula SymPy

I am trying to use this formula to calculate integrals with SymPy, but i can't find a way how to do it. This formula is used in FEM to calculate matrices and vectors on triangles. Maybe you know another way how to use symbolic calculus for this…
0
votes
1 answer

How to use two UEL in same problem in ABAQUS

I am solving a contact problem. A footing under compressive load placed over soil. I would like to use two different types of user-defined elements one for the soil and other for the footing. Please note: I would like to use two different UEL. Not…
0
votes
1 answer

Is three.js able to show FEA results?

I’m trying to develop an real-time wind turbine simulation Finite Element Analysis software using three.js to display the caculated 3D FEA results. Like the image one frame of FEA result.the showed 3D wind turbine rotates like the real physical one…
alex zhang
  • 23
  • 4
0
votes
1 answer

How to sum the duplicates in coo matrix efficiently in fortran?

I have a sparse matrix in the coo format obtained from finite element assembly. It has many duplicate entries in the triplet format. I want to sum them up efficiently using Fortran. I went through the sparse kit package but I didn't find a…
bharat
  • 21
  • 1
0
votes
1 answer

How to parallelize the nested do loops with a counter variable

jx=1 (Counter variable) ! squeeze out zeros from k !$omp parallel do do ix=1,tdof do iy=1,tdof if (iy.ge.ix) then ! Upper triangular elements if ((abs(K(ix,iy))-0.d0)>1e-10) then ! Non zero entries …
bharat
  • 21
  • 1
0
votes
1 answer

Applying displacement control loading using Lagrange multipliers in the material non-linear finite element method

Hi I am trying to implement a simple plasticity based finite element code. I am not clear how to set up displacement control applied through Lagrange multipliers. In case of a linear problem, I did the following to apply displcament control and it…
0
votes
1 answer

Variable dependent flux boundary condition in fipy

I have a case where reaction flux is present at a boundary and flux is proportional to the concentration variable c. Will the following correctly apply the boundary condition: c.faceGrad.constrain(k*c.faceValue, mesh.physicalFaces["ReactingFace"]);
DKS
  • 188
  • 10
0
votes
0 answers

How to find number of nonzero terms by using element connectivity in FEM (language : C)?

In Finite element method, for making sparse matrix of tangent stiffness matrix we should know the number of nonzero terms because, I am using language : C. I could not find an optimal algorithm to find the number of nonzero terms using "element…
Song
  • 13
  • 3