Questions tagged [solid-mechanics]

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

Filter by
Sorted by
Tagged with
2 votes
0 answers
87 views

Iterative solvers for problems in solid and structural mechanics

I am looking for comprehensive literature (papers, books, reports etc..) on iterative solvers for solid and structural mechanics problems to understand the best iterative solvers and preconditioners ...
Chenna K's user avatar
  • 934
0 votes
1 answer
51 views

How to implement the interface extension of fluid "displacement" in ALE?

In ALE, we first set a referenced space for fluid, then we extend the boundary fluid displacement to the whole fluid region, take harmonic extension as an example, we need $$\Delta \left ( \hat{u} \...
吴yuer's user avatar
  • 193
1 vote
0 answers
33 views

Recommendations for some new books about computational contact mechanics in solid mechanics

I want to simulate some frictional contact problems, but I'm not familiar with this field, could you please recommend some new books as introductions?thank you
吴yuer's user avatar
  • 193
0 votes
0 answers
19 views

Is there any papers about simulating material damage under supersonic combustion?

I'm quite interested in this field, Is there any papers about simulating the damage of material(such as alloy) under supersonic combustion, after searching for a while, I didn't find anything
吴yuer's user avatar
  • 193
2 votes
1 answer
111 views

FENICSx or deal.ii for modeling fluid structure interaction of cylindrical energy harvesters in various arrangements

I'm trying to model fluid-structure interaction of cylindrical energy harvesters in various arrangements. I have chosen OpenFOAM for the fluid part but I have no clue which framework to use for the ...
Mechanic7's user avatar
  • 123
2 votes
1 answer
190 views

Neumann BC in the current configuration in a finite-strain problem

For a hyperelastic problem, I understand the variational formulation can be written as the minimisation of $\Pi$ with $\Pi = \int_{\Omega} \psi( \pmb{u} )dx - \int_{\partial\Omega} \pmb{T}\cdot \pmb{u}...
Zed's user avatar
  • 23
0 votes
2 answers
166 views

Determination of the domain of nonlinearity in a Neo-Hook solid model (Finite elements)

For a FEM simulation of a Neo-Hook solid model, how do we know we are in the "regime" of nonlinearity of the solid? In other words, how do I know the hyperelastic material law is really used,...
Zed's user avatar
  • 23
0 votes
0 answers
33 views

RFEM: Equilibrium over plate section

I'm calculating a simple plate arrangement in RFEM: The plate has dimensions of 1m X 1m, thickness of 50mm. Elastic modulus of the material is 205000.0 N/mm2, Poisson's ratio is 0.3. The plate is ...
S. Rotos's user avatar
  • 133
4 votes
1 answer
184 views

Estimating forces on a model from the displacements of nodes

In any FEM problem involving mechanics, we try to solve the differential equation for the displacement field, $u$ given the force vector in the nodes, $F$. In industry, we often see our automobiles ...
s6292_1997's user avatar
0 votes
0 answers
78 views

How to solve evolution equation numerically?

How do I solve the following evolution equation numerically: $$ \dot{\boldsymbol{\mathcal A}} = -\lambda A_1|\boldsymbol{s}-\boldsymbol{\alpha}|^2\left[(\boldsymbol{n}_r:\boldsymbol{\alpha})\...
Mahmudur Rahman's user avatar
3 votes
1 answer
274 views

Is it really necessary to solve a system of linear equations in the Finite Element Method?

When we solve some boundary value problem by Finite Element Method, the appropriate system of linear equations is built, $$Ax=b.$$ Usually we use the solution x just for plugging it into some ...
Fidel Pestrukhine's user avatar
2 votes
1 answer
368 views

Finite Element Modelling of Hyperelastic Material under 2D Plane Strain Conditions

I am currently working on writing a MATLAB code for running a finite element simulation of a hyperelastic material in 2D. Since I am building this simulation as a part of a fluid-structure interaction ...
Mohamed Abdelhamid's user avatar
0 votes
1 answer
97 views

Calculate strains based on X,Y Deformation gradient with time

Recently I have obtained a csv value form and experiment I have computed. However, I was trying to understand how each individual component being calculated. The image attached shows sample point ...
Abdullah's user avatar
1 vote
3 answers
501 views

If FEM is exact at the nodes, why do first and second-order elements give very different results?

I'm looking at the solution to a structural mechanics problem that is modeled with first-order elements and then as a comparison with second-order elements. It is clear that the first-order elements ...
Tepa's user avatar
  • 31
6 votes
1 answer
269 views

Why FEM for incompressible materials is ill-posed?

I am an engineer who is trying to get a deeper understanding of FEM. I have been using the Zienkiewicz texts as my bible. It touches on the issue of incompressibility but I need a more intuitive way ...
TheCodeNovice's user avatar
1 vote
2 answers
342 views

What is the difference between non-linear elastic simulation and linear elastic simulation with plasticity?

I'm learning how to do Finite Element calculations using Comsol Multiphysics. In Comsol, Linear Elastic Material and Nonlinear Elastic Material are available as material models: Using Linear Elastic ...
S. Rotos's user avatar
  • 133
1 vote
1 answer
136 views

What are the prerequisites and resources to self-learn the Boundary Element Method for Contact Mechanics problems?

What are the prerequisites to learning BEM? In which order is it advisable to learn BEM and FEM - either one before the other, or does it not matter? What are some good resources to self-learn BEM? P....
SNIreaPER's user avatar
1 vote
0 answers
127 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} \...
farid1889's user avatar
1 vote
1 answer
106 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 ...
Ali Baig's user avatar
1 vote
0 answers
77 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$ ...
Olumide's user avatar
  • 265
0 votes
1 answer
75 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 ...
S. Rotos's user avatar
  • 133
2 votes
1 answer
310 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 ...
Olumide's user avatar
  • 265
3 votes
2 answers
815 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$ ...
Olumide's user avatar
  • 265
3 votes
2 answers
413 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 $\...
FEGirl's user avatar
  • 333
3 votes
2 answers
277 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 ...
Mohamed Abdelhamid's user avatar
2 votes
2 answers
55 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 ...
S. Rotos's user avatar
  • 133
1 vote
0 answers
204 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 ...
QuestionsQuestion5839's user avatar
2 votes
0 answers
73 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. ...
Biswajit Banerjee's user avatar
1 vote
2 answers
795 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 ...
Bruce Lee Jun Fan's user avatar
0 votes
0 answers
237 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 ...
Srikumar Gopalakrishnan's user avatar
1 vote
0 answers
37 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 (\...
Olumide's user avatar
  • 265
0 votes
1 answer
138 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 ...
Donald Zeka's user avatar
1 vote
1 answer
143 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 ...
Dude's user avatar
  • 580
2 votes
1 answer
608 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 ...
akr's user avatar
  • 55
3 votes
3 answers
135 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 ...
student010101's user avatar
0 votes
0 answers
52 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 ...
Bruce Lee Jun Fan's user avatar
0 votes
1 answer
1k 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 ...
Ahmed Hossam's user avatar
1 vote
0 answers
65 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 ...
Moum's user avatar
  • 11
1 vote
3 answers
557 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 ...
Xu Hui's user avatar
  • 323
7 votes
1 answer
3k 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 ...
Sagar Shah's user avatar
0 votes
2 answers
89 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 ...
Xu Hui's user avatar
  • 323
3 votes
1 answer
174 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) ...
rasmodius's user avatar
  • 131
4 votes
3 answers
557 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 ...
Ben's user avatar
  • 43
0 votes
0 answers
37 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 ...
John's user avatar
  • 1
4 votes
1 answer
184 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 ...
Alexander's user avatar
1 vote
0 answers
43 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 ...
user5965026's user avatar
0 votes
1 answer
228 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 = ...
HokieFan7's user avatar
1 vote
0 answers
99 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 ...
Yirmidokuz's user avatar
4 votes
1 answer
110 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) ...
Xu Hui's user avatar
  • 323
3 votes
1 answer
114 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 ...
secondrate's user avatar