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
0 answers

ABAQUS python issue extracting deformed coordinates using COORD

I have produced a script for ABAQUS CAE using Python to extract the deformed coordinates of a node set called 'SET-1'. A load is applied to the top surface of the cube and the cube deforms. The field output 'COORD' is supposed to give the deformed…
0
votes
0 answers

How to specify current density in Matlab's PDE's Toolbox/FEM solver?

I have a few questions about a specific part of the Matlab FEM solver because the documentation isn't clear from what I can see: For this specific problem statement I cannot figure out what are the units of the current density specified in the…
0
votes
1 answer

Quad mesh for an annular disc using gmsh

I am trying to generate a quad mesh grid for an annular disc geometry. I was able to generate this mesh for a rectangular shape using gmsh: import gmsh import sys gmsh.initialize() gmsh.model.add("t1") lc = 1e-2 p1 = gmsh.model.geo.addPoint(0, 0,…
H.H
  • 188
  • 1
  • 13
0
votes
0 answers

finite element analysis in mathematica

I have been having trouble solving a coupled set of PDEs in mathematica (electrical current flow), and the issue appears to be in the boundary conditions. Specifically, I have PDES as follows: sum_{i} \partial_i J_i = 0 J_i = sum_{j} c_{ij}…
Rain
  • 265
  • 2
  • 16
0
votes
0 answers

Calculation of the nodal displacements using python and direct stiffness matrix method

I analyse the folowing Frame using python and direct stiffness matrix method. I use the following definition of coordinate system: The values of displacement, which I obtain is correct, but the direction of nodal displacements in to x-direction are…
0
votes
0 answers

How can I extract facets and their associated nodes from a 3D mesh using NX Open?

I'm quite new to NX Open so pardon me if this is a trivial question. I have a 3D part file that I meshed using tetrahedrons, and applied some loads and constraints. The solver I will be using is FENICSx. The way FENICSx works in my project is that…
Haridev V
  • 3
  • 3
0
votes
0 answers

The strain displacement matrix for beam element

I want calculate the strain and stresses in a beam element with following DOF {u1,v1,O1,u2,v2,O2} using FEM analysis. Stresses: Strain: [ε] = [B]x[u] with: [B] = strain displacement matrix [u] = calculated displacement of a beam element I can't…
0
votes
0 answers

How to Solve Coupled Equaitons Using Weak Form of Comsol?

I want to solve a set of coupled equations using weak form of Comsol. Some tuitials are available for solve simple uncoupled equations using weak form in comsol, which doesn't seem to work for my case. The equations may look like this (just a…
HC6
  • 21
  • 1
0
votes
0 answers

Efficient method to get node connectivity in FEM

I have information of element connectivity in problem of finite element method. For example, 1, 5, 6, 8, 7, 1, 2, 4, 3 The first column is the index of element, and the last 8 columns are the node index which composed entries for element 1. The goal…
Song
  • 13
  • 3
0
votes
1 answer

How to Create Piecewise Function in SymPy with Intervals

i need to create a piece-wise function inside an interval but sympy piecewise can't use and (&). I read that the function can't recieve Boolean values so I tried to add them together and it doesn't seem to be right. The code is as follows: import…
0
votes
1 answer

Library for numerical Integration of a function over a tetrahedron (3D)

I am looking for å python library that has a function to solve 3D integrals over a tetrahedron. I would like to be able to input four points on the form (x, y, z) and a function f(x, y, z) where f is a polynomial function. Have only found functions…
0
votes
1 answer

"Error using vertcat Dimensions of arrays being concatenated are not consistent." in matlab while using integral2

I'm coding a programme about FEM and here is my code FEM_2D_TRI_QUA_1.m N1=2; N2=2; N=2*N1*N2; top=1; bottom=0; left=0; right=1; h1=(right-left)/N1; h2=(top-bottom)/N2; T=zeros(3,2*N1*N2); N10=N1+1; N20=N2+1; for i=1:N2 for j=1:2*N1 …
YuerWu
  • 11
  • 1
0
votes
1 answer

An algorithm suggestion for interpolation in a 2d plane, when the plane is in 3D space

I'm trying to find value for a point within a triangle or square. By "value" I do not mean coordinates. Suppose there is a value (number) assigned to each node of the square/triangle. The square/triangle is on a plane. How can I interpolate values…
Neg K
  • 17
  • 7
0
votes
1 answer

How to setup MATLAB shape function solver

I am trying to find the shape function for a finite element analysis problem in MATLAB. The basic shape function looks like, with the boundary condition on the left. At Ni=N1, x=x1, y=y1 (the coordinate of node 1), the shape function would be equal…
Faito Dayo
  • 115
  • 4
0
votes
0 answers

FEA Stress plot in Python from 3 1D-arrays

I have 3 1D arrays (node x-coordinates, node y-coordinates and Von-Mises stress scalar) exported from an FEA solver. I want to create 2D contour plots as shown below in Python: Stress plot example I have managed to create such plot as shown…