Questions tagged [solid-mechanics]

Studies the behavior of solid materials, and deformation under the action of forces

Filter by
Sorted by
Tagged with
0 votes
0 answers
55 views

Derivation of the second Piola-Kirchhoff tensor

I tried many formulas to find the components of the second Piola Kirchhoff. I need help to derive equations 27-29 on reference 1. [![enter image description here][1]][1] we have $$ \mathbf{C}=C_{11} \...
user avatar
1 vote
1 answer
48 views

Applying Stress Boundary Conditions in Commercial Finite Element Analysis Codes

I am trying to replicate a finite element analysis given in a research paper titled On the Detection of Stress Singularities in Finite Element Analysis 1 by G.B.Sinclair et. al. The geometry of the ...
user avatar
1 vote
0 answers
65 views

Stretched Elastic Sheet with Horizontal Cut in Interior

Given a rectangular finite elastic sheet $ABCD$ containing an arbitrary horizontal cut (or discontinuity) $EF$ in its interior. The edge $AB$ of the sheet is tethered but $CD$ is stretched by $\Delta$ ...
user avatar
  • 215
0 votes
0 answers
34 views

3D Euler Bernoulli Beam for Nonlinear FEA

Anybody has any experience in coding 3D beam elements? I am trying to write a C++ code for a 3D euler bernoulli beam. For 2D, I used Reddy for coding 2D for non linear FEA. How should I proceed with ...
user avatar
0 votes
1 answer
70 views

Beam theory: does finer meshing make any difference if the shape functions used are 3rd degree polynomials?

I'm studying the Finite element method for structural mechanics, and I'm reading this source. 3rd degree polynomials are used as shape functions for beam theory, since they are required to be ...
user avatar
  • 123
2 votes
1 answer
142 views

Total stored potential energy of finite element mesh from nodal point displacements and strain energy density function only

I am interested in calculating the total potential energy stored in a finite element mesh given its nodal point displacements alone. The forces that created the displacements are irrelevant because ...
user avatar
  • 215
3 votes
2 answers
295 views

Second Piola-Kirchoff Stress Tensor of Neo-Hookean solid at "zero deformation"

The strain energy of an incompressible Neo-Hookean solid is given as: $$ W = C_{10}(I_1 - 3) $$ Implying that at zero deformation $W = 0$, because $F = I \implies C = F^TF = I \implies I_1 = 3$ ...
user avatar
  • 215
2 votes
2 answers
390 views

Modelling question: example of a physical phenomenon with this jump condition at an interface?

in our finite element class we were talking about interface problems our teacher came up with the following, where $K_i$ are two given functions and $u_i$ is the restriction of the solution $u$ to $\...
user avatar
3 votes
2 answers
201 views

How to enforce fluid and solid dynamic coupling in fluid-structure interactions using the finite element method?

I apologize in advance if the question has been posted before or if it sounds a bit naive. I am writing my own code in MATLAB for a staggered finite element solver for fluid-structure interaction ...
user avatar
2 votes
2 answers
45 views

How does RFEM give non-linear results with a two-node mesh?

I did this simple analysis on RFEM, of a rigid-supported beam loaded with a point moment. Before analysis, I didn't assign any kind of mesh manually. When I turn on the FE Mesh visible on Project ...
user avatar
  • 123
1 vote
0 answers
95 views

Varying Young modulus in FEM simulation

I'm working on a project for which I have inherited some FEM code. This implemented FEM calculates, given some force field, displacements on a discretised square using square elements and assumes a ...
user avatar
2 votes
0 answers
50 views

Extracting a mid-plane for thick shell analysis

I have a complex part that contains features of the form shown in the figure below. Because of the cost of 3D finite element simulation of the part, I want to try an analysis with 2D thick shells. ...
user avatar
1 vote
2 answers
291 views

Dyadic operations, fourth order tensors and Tensor algebra

I am trying to understand the dyadic operation for a while since I am interested in Elasticity problems. I believe an intuitive understanding (rather than assuming) will give me good problem solving ...
user avatar
0 votes
0 answers
152 views

How to interpolate stress at unknown points from the stress values available based on geometrical position for constant load?

I am working on a combined contact, bending, and torsion problem. I have data on geometrical points and their instantaneous stress components. However, based on the available data, I have to ...
user avatar
1 vote
0 answers
32 views

Approximation in the derivation of the Arc Length method

I am studying the proof of the Arc Length method in section 2.2 of this thesis. In equation (2.2) the author introduces the supplementary conditions $$ (\Delta {\bf u} + \delta {\bf u})^T \cdot (\...
user avatar
  • 215
0 votes
1 answer
88 views

relation between different tangent stiffness

I need to find a relation between the tangent stiffness $L_1$ of the first Piola-Kirchhoff stress tensor with the tangents stiffness $L_2$ of the second Piola-kirchoff stress tensor. They are defined ...
user avatar
1 vote
1 answer
117 views

Transition from 2D to 3D finite element code, what are the inevitable modifications to be implemented?

Imagine we have a simple 2D FEM solver (we are dealing with solid mechanics) and we would like to develop it to a 3D FEM solver (let's say for the same solid mechanics problem) in this case what are ...
user avatar
  • 418
2 votes
1 answer
292 views

How to solve a linear problem A x = b in PETSC when matrix A has zero diagonal enteries?

I am solving a structural mechanics problem that involves setting constraints, and I use Lagrange multipliers to set it. Consequently, some diagonal entries of the tangent stiffness matrix vanish, and ...
user avatar
3 votes
3 answers
111 views

How are the classical set of equilibrium equations for linear elasticity derived?

In linear elasticity, the governing PDE is the equilibrium equations (absent of vibration considerations): $$ -\nabla \cdot \sigma = F $$ Is this equation simply derived from the sum of forces and ...
user avatar
0 votes
0 answers
45 views

Hi I am trying to model a 2D Lug angle using Gmsh 4.6. How can I combine transfinite quad and regular full quad meshes in the following geo file?

I need transfinite mesh a small section of the bolt hole to insert a crack. However, The transfinite mesh and regular full quad mesh seem being incompatible and throwing errors. How can I combine ...
user avatar
0 votes
1 answer
642 views

Global stiffness matrix from element stiffness matrices for a thin rectangular plate (Kirchhoff plate)

I have the element stiffness matrix for a thin "kirchhoff" plate. The plate is 3 [m] x 5 [m] and is simply supported on all edges. It's thickness is 0,2 [m]. On the plate there acts a ...
user avatar
1 vote
0 answers
56 views

implementation of shell elements in a topology optimization algorithm

I am working on developing a topology optimization solver based on the finite element method and I want to add a triangular shell element in it. I used the classic finite element method but I didn’t ...
user avatar
  • 11
1 vote
3 answers
373 views

How to determine global stiffness matrix is constrained or not

Background In solid fem, we often solve $$\mathbf{Ku}=\mathbf{p}$$ where $\mathbf{K}$ is global stiffness matrix, $\mathbf{u}$ is displacement, $\mathbf{p}$ is global load vector. If displacement not ...
user avatar
  • 323
5 votes
1 answer
1k views

4th order tensor rotation - sources to refer

I am trying to model a linear elastic material in Abaqus using a UMAT. For my application, I need to rotate the 6x6 compliance matrix for a given set of eigenvectors (or a rotation matrix). I came ...
user avatar
0 votes
2 answers
74 views

Uniaxial stretching solution not uniform in FEM code

I am trapped here for a long time. I wrote a toy Matlab FEM code. I want to run the follow simulation. Mesh Suppose we have a cube, and we divide it into subcube along $x,y,z$ axis, then each subcube ...
user avatar
  • 323
3 votes
1 answer
156 views

FEM with elastic inhomogeneous properties leads to mesh-induced anisotropy

I'm solving an elastic homogenization problem and I'm having problems with mesh artifacts. I would like to first give a brief summary of what I do: I have a system with inhomogeneous (but isotropic) ...
user avatar
  • 131
4 votes
3 answers
225 views

I wrote a 2D Finite Element program for Axial Loaded Plates, but the results are unexpected

TLDR: I used Python to write a 2D Finite Element program using 'Constant Strain Triangles' and my beam keeps pointing slightly upwards instead of straight sideways (like the force). I'm new to FEA and ...
user avatar
  • 43
0 votes
0 answers
32 views

At what l/d ratio will a frame element start to behave as a shell element?

I'm working in ETABS. There are few columns of dimension 300mmx1400mm. The height of building is 36.6 meter above ground level and the dimension of building is 26mx68m. I'm getting the time period of ...
user avatar
  • 1
4 votes
1 answer
142 views

Why does smoothed aggregation multigrid method used as preconditioner in conjugate gradient slows down the solution time?

I'm solving a system of linear equations obtained from the FEM discretization of a simple linear elasticity problem on a cube with zero displacements at one plane and a load on the opposite one. The ...
user avatar
1 vote
0 answers
38 views

Well posedness dirichlet BC in elasticity

I always thought that a Dirichlet BC need to be prescribed on at least one part of the domain in order for linear elasticity to be well-posed. Is this correct? But I was recently told by a colleague ...
user avatar
0 votes
1 answer
153 views

Determining Displacement Field on a Sphere

I need to find the displacement field for this sphere shape in terms of $\delta$. So far, by applying boundary conditions, I know $u_r = u_\theta = u_\phi = 0$ at $r = a$ and $u_r = \delta, u_\theta = ...
user avatar
1 vote
0 answers
76 views

When does reduced integration lead to artificial zero energy modes in stiffness matrix?

This question relates to the topic of locking free finite element development. In the case of application of reduced integration to global stiffness matrix for the Timoshenko beam element with ...
user avatar
4 votes
1 answer
95 views

Why the solid FEM problem can not be solved after constraining 3 degrees of freedom?

I write a simple MATLAB code for solving solid FEM problem. The problem looks like that (1) (2) x-------x | / | | / | | / | x-------x (3) (4) ...
user avatar
  • 323
3 votes
1 answer
97 views

Dynamic linear elasticity problem is stiff (numerically)

I am considering a dynamic linear elasticity problem applied to a simple structure such as a beam. In system form, the PDE can be written as $M \ddot{X} + D \dot{X} + XK = F$ where $X$ represents the ...
user avatar
2 votes
0 answers
322 views

Extracting FEM matrices in matlab pde toolbox

I am trying to follow the dynamic linear elasticity in Matlab, link here. My goal is to extract the FE Matrices using the function assembleFEMatrices in matlab and solve the resulting system of second-...
user avatar
1 vote
2 answers
691 views

Evaluation of slope at iteration ith - Newton-Raphson method

I'd like to know how Ansys computes the slope (=stiffness matrix) at point x1 in figure. I'm studying the way in which Ansys uses the Newton-Raphson method when there are nonlinearities. In the slide ...
user avatar
3 votes
1 answer
221 views

Reference request: Riks method (Nonlinear FEM)

I'm struggling to find a good detailed reference explaining the Arc-length method or, more generally, Riks method and its derivations. I looked for the classical books in nonlinear mechanics (the ones ...
user avatar
3 votes
0 answers
108 views

Strain propagation from surface to bulk in COMSOL

I am trying to simulate strain propagation from the surface into the bulk. I have a rectangular semiconductor block (~2 μm thick) on top of which metal gates (~25 nm thick) are deposited as seen in ...
user avatar
1 vote
1 answer
137 views

Abaqus, ANSYS, and FVM solver for thermal expansion problem converges to different values

Is it reasonable for a FEM and FVM code to converge to slightly different solutions for the same physical problem (identical BCs, geometry, properties, etc...), provided stability constraints are ...
user avatar
1 vote
1 answer
130 views

Split solution of FEM problem depending on number of DOF

Assume we have a 3D finite element structural problem discretized with hexahedral elements with 8 nodes and 3 degrees of freedom per node. Instead of solving the global stiffness matrix system for all ...
user avatar
0 votes
0 answers
30 views

Numerically computing deflection due to thermal expansion

Using linear elasticity formulation, I am attempting to numerically compute the displacement due to thermal expansion. This is done for a 3-D isotropic material. The governing equations are simply: ...
user avatar
1 vote
1 answer
234 views

Isotropic thermal expansion

I frequently see the equation $$ \sigma_t = E\alpha \Delta T $$ as the equation for thermal stress. Where $E$ is Young's modulus, $\alpha$ is the CTE, and $\Delta T$ is the change in temperature. ...
user avatar
0 votes
0 answers
56 views

Strange spectral (finite) element results of a solid plate

I tried to implement the spectral element method1 as proposed in [1]. I simulated an aluminum plate and a vertical concentrated force was applied at the middle of the top left quarter of the plate2. I ...
user avatar
2 votes
1 answer
117 views

Linear elasticity modeling load using traction vs. mixed BC

In classical linear elasticity, when modeling a force/load boundary condition, it appears that we could either: Use a pure Neumann boundary condition, where the 3 traction components are specified. ...
user avatar
6 votes
3 answers
903 views

Why is the FVM traditionally used in CFD, and FEM in computational structures?

Most CFD codes use FVM. Most computational structures codes use FEM. Why is the FEM not frequently used in CFD, and why is FVM not frequently used in FEM?
user avatar
1 vote
1 answer
105 views

Structural boundary conditions - rotational/translational DoFs and displacement/tractions BCs

I am a little bit confused over the concept of translational and rotational degrees of freedom (DoFs) in structures, and their relation to displacement/traction BCs. Do displacement boundary ...
user avatar
3 votes
0 answers
93 views

Structural Analysis Library

Can anyone recommend a structural analysis library that satisfies the following requirements: C++ API Simulate both beam elements and shell (slab) elements Both static and dynamic analysis Free and/...
user avatar
0 votes
1 answer
83 views

Given co-ordinates of 8 vertices, how to calculate the outward normal and surface area for each face of a irregular hexahedron?

I am working on an FEA mesh of hexahedron elements. The elemental level calculations involve finding the surface normals and area for each surface of a hex element. I preferred the vector cross ...
user avatar
2 votes
1 answer
157 views

Is steady linear elasticity inherently ill-conditioned?

Compared to the transient PDE for linear elasticity, the steady equations appear to less well-conditioned. Are they inherently ill-conditioned without the transient term? The condition number for ...
user avatar
1 vote
1 answer
214 views

Strain from FEM simulations to strain gauge measurements

I am looking for some intuition in making comparisons of FEM simulations to experimental measurements. In particular, I am interested in comparisons to strain gauge readings, and perhaps even LVDTs. ...
user avatar