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

Polynomial functions as matrix elements

Matrix N is a 2x2 matrix that contains polynomial functions: N=[(1-x).*(1-y), 0 ; 0 , (1-x).*(1-y) ]; where x and y are defined over the interval [0:1] to make a square element I also have an operator matrix, D D= [d/dx, 0 0, d/dy]; I want…
0
votes
0 answers

Python not applying regex on large files

import re import sys, getopt import mmap shakes = open(sys.argv[1:][0],'r') love = open(sys.argv[1:][1], "w") #moreLove = open (sys.argv[1:][2], "w") #HardLove = open (sys.argv[1:][3], "w") node = re.compile('\*NODE[a-zA-Z,…
0
votes
0 answers

Meshing options to generate number of the sides of and element (tetgen-triangle)

I wrote a finite element code in fortran 90. This code is really fast, except the meshing process. I used triangle and tetgen for meshing in 2D and 3D, respectively, so this process is fast, of course. For example, for the unit square [0,1]x[0,1] in…
user106306
  • 315
  • 2
  • 10
0
votes
0 answers

ABAQUS: Balancing a thin plate by two fixtures on edge centres under self weight

I've created a plate of dimensions(200x200x5mm). Density=5000. The plate is inclined at 30 degrees with x axis. Now, just under the influence of self weight, the plate is put in equilibrium using two supports at the midpoints of the faces shown in…
0
votes
0 answers

Ansys LS-DYNA explicit dynamics element COMBI165 length is affecting the solution

I'm analysing some dynamic models using Ansys LS-DYNA explicit dynamic models. My model only has beam elements (BEAM161), concentrated masses (MASS166) and springs/dampers (COMBI165). The problem is that the length of the springs/dampers is…
0
votes
0 answers

Abaqus composite ply strain extraction

I'm trying to extract ply strains from a composite shell model in Abaqus. I would like to be able to export (to a .dat (preferably) or any sort of text file) every ply strain in every element. So far I have managed to find a method to extract the…
user2739143
  • 591
  • 2
  • 8
  • 13
0
votes
1 answer

Using implicit difference method to solve the heat equation

I trying to make a Matlab code to plot a discrete solution of the heat equation using the implicit method. The information I am given about the heat equation is the following: d^2u/d^2x=du/dt Initial conditions (t=0): u=0 if x>0 Otherwise u=1 (when…
David
  • 159
  • 8
0
votes
0 answers

Meshing a construcitve solid geometry (CSG)

I have many geometric primitives (100) which cut a cube. This could be modelled using CSG. Now I need a relatively coarse tetra mesh for finite element simulations. For generating the mesh I came across netgen (see the picture). But netgen will give…
user3528789
  • 103
  • 7
0
votes
1 answer

Matlab - FEM: Meshing polygon

I have problem with meshing this polygon: I need fill this polygon with one of this finite element: I don't how to type a code for this problem and then plot this polygon with mesh, all without using PDEtool. Thank you.
MrPitivier
  • 71
  • 9
0
votes
1 answer

Solving nonlinear FEM in MATLAB

I try to solve a heat diffusion problem on tetrahedral finite elements with nodal heat sources, which depend on the solution vector, in MATLAB. The nonlinear equation system looks like this: BU' + AU = q(T) with B being the heat capactiy matrix, A…
0
votes
1 answer

2d linear convection in c++ (segmentation fault)

im following Dr. Lorena Barbara's 12 steps to Navier stokes equation (http://lorenabarba.com/blog/cfd-python-12-steps-to-navier-stokes/) structure but I'm not sure how to do the 2d linear convection. I was wondering if anyone here would be familiar…
user2804865
  • 976
  • 2
  • 9
  • 15
0
votes
1 answer

Indentation error

I receive an indentation error that I cannot figure out the reason. The error is ('unexpected indent', ('C:/Hamid/Failure_index.py',15,1,'\tSDV2=xyList[0]\n')). My code is from abaqusConstants import * from odbAccess import * from…
0
votes
0 answers

MPI double a map from a given GID and LID list

I am working with a finite element code, where the nodes of the mesh are distributed according to a map (let us call it single_mesh_nodes) with local and global identifier. I need to make a new map (let us call it double_mesh_nodes) which should…
0
votes
1 answer

Check for edge connectivity in a 2D mesh

I'm trying to optimize structures by using FEM and genetic algorithms (GA). The FEM solver is a commercial one, and I'm programming the GA. Eventually I hope to get something like this. My first approach is simple, just erasing elements and…
0
votes
0 answers

I am working with a code written by another person and there are things that I can not understand

I am working with a parallel finite element code that was written by another person, and there are some pieces of code that I can not understand. In particular there is the map creation for the degrees of freedom belonging to each node of the…