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

How to solve "Cannot open include file 'odb_API.h'" ABAQUS 2018

The problem is solved by watching this video: https://www.youtube.com/watch?v=qGNAdlDoeTs I've tried linking Abaqus 2018 with Fortran and the subroutine check passed. But I got this error when running my UMAT: "Error in job UMAT: Problem during…
0
votes
2 answers

What are some software / techniques I can use to take a chunk out of a volumetric 3D mesh?

For some work using the Finite Element Method, I manually generated a volumetric 3D mesh with tetrahedral elements, from medical image data. The object has an arbitrary shape that is almost cylindrical. I also have a separate mesh file for only the…
Cogicero
  • 1,514
  • 2
  • 17
  • 36
0
votes
1 answer

Polygon Mesh of a 3D Model for Finite Element Analysis

I have a complicated 3d shape which I have created in Autocad that needs to be analyzed in SAP2000. SAP2000 is only capable of reading a 3D solid made out of Polygon Mesh. The present Autocad Versions can only mesh as Polyface meshes. Is there any…
0
votes
2 answers

For a buckling analysis, must all forces be multiplied by the resulting eigenvalue, or only the compressive load?

I am trying to do a linear buckling analysis (sol 105) with Nastran on a cylindrical shell structure. My understanding is that the compressive load that I apply to the structure must be multiplied by the resulting eigenvalue to get the buckling…
Jeroen
  • 119
  • 2
  • 2
  • 8
0
votes
0 answers

Overlaying FEA on video stream

I got a video stream from a camera. I would like to overlay an FEA like 3D of one object that showing the stresses on the object dynamically changes. Do you guys have any idea how to implement this?
0
votes
1 answer

ParaView shows displacement as color

In my current "Partial Differential Equations" lecture we write software to solve PDEs numerically. In the process we had to approximate a function f by a continuous, piecewise linear function. Our software has the capabilities to export the…
s1624210
  • 627
  • 4
  • 11
0
votes
0 answers

Slow computation on google colab while solving partial differential equation

I 'm using google colab to solve the homogeneous heat equation. I had made a program earlier with scipy using sparse matrices which worked upto N = 10(hyperparameter) but I need to run it for like N = 4... 1000 and thus it won't work on my pc. I…
0
votes
1 answer

Is that way of writing/reading a solution safe?

I need to store a solution of expensive FEM calculation to use it in further analysis. Browsing through tutorials I have so far discovered, that I may store my results like this: from fenics import * mesh = Mesh('mesh/UnitSquare8x8.xml') V =…
abukaj
  • 2,582
  • 1
  • 22
  • 45
0
votes
1 answer

Finite element analysis: How can I find the boundary conditions and settings for the knee joints?

I am not an expert in finite element models (FEMs), however, I am reading related articles and trying to find the boundary conditions for knee joint including the bone and cartilage. Does anyone have any suggestion? How can I find the boundary…
S.EB
  • 1,966
  • 4
  • 29
  • 54
0
votes
0 answers

An improvement of a conjugated gradient-like approach?

hello everyone the problem is slightly complicated to describe in details but i will provide as much details as i can, just ask if you need more information. In the project i am involved in, i had to work on an algorithm that couples 2 phenomena: a…
0
votes
1 answer

Cantilever Beam Height versus Deflection

I have a problem whereby in a LINEAR analysis, a cantilever is vertically bent through a tip load. A parametric study is required between... Change in Tip Deflection (dZ) versus Change in Beam Height (dH) I want to see the SENESITIVITY is more…
0
votes
0 answers

generating a 2D heptagonal finite element mesh in a rectangular domain

I want to generate a 2D finite element mesh which consists of heptagons and other polygons except triangles, bounded by a rectangle. I am very new to coding I don't know how to start. Can anyone please help to devise an algorithm to produce the…
0
votes
0 answers

No speedup with parfor

I wrote a finite element toolbox in Matlab which turned to be rather slow for large meshes, so I decided to parallelize the element matrices assembly. Hence, after initiating the workers pools I use cell arrays in order to build the global matrices,…
0
votes
1 answer

GMSH 3D surface mesh

I've been trying without success to create a 3D SURFACE mesh using GMSH (v.3.0.6). The problem that I am having is, after creating the surface and generating the mesh, when I read the .MSH file, I get a bizarre node numbering, i.e., not all the…
Gabs
  • 292
  • 5
  • 23
0
votes
1 answer

FileNotFoundError: [WinError 2]

I am trying to execute the following code from https://pypi.org/project/pygmsh/ import pygmsh import numpy as np geom = pygmsh.built_in.Geometry() # Draw a cross. poly = geom.add_polygon([ [ 0.0, 0.5, 0.0], [-0.1, 0.1, 0.0], [-0.5, 0.0, 0.0], …