3.6.20 · D1Spacecraft Structures & Systems Engineering

Foundations — FEM software — NASTRAN, ABAQUS (concepts and use)

3,403 words15 min readBack to topic

Everything in the parent note (FEM software — NASTRAN, ABAQUS (concepts and use)) is built on a handful of ideas. This page defines every symbol and word from absolute zero, in an order where each one leans only on the ones before it.


1. Displacement — "how far did a point move?"

Imagine drawing a dot on a metal bracket, then pushing on the bracket. The dot slides to a new spot. The little arrow from where the dot was to where it is now is called the displacement.

Figure — FEM software — NASTRAN, ABAQUS (concepts and use)
  • Picture: the red arrow in the figure — tail at the old position, head at the new one.
  • Why the topic needs it: displacement is the thing we solve for. Everything else (stress, whether the bracket cracks) is computed from how points moved. When the parent writes in , this is it.

The bold letter means it is not a single number but a list of numbers — one for each direction a point can move. That leads straight to the next idea.


2. Degrees of Freedom (DOF) — "how many ways can it move?"

A single point in 3D space can slide three ways: left–right (), forward–back (), up–down (). Each independent way of moving is one degree of freedom.

Symbol check
The bold is a column of all the DOF displacements of all the nodes, stacked one on top of another.

3. Nodes and Elements — "chop the smooth thing into pieces"

You cannot store a smooth continuous shape in a computer. So FEM covers the structure with a grid.

Figure — FEM software — NASTRAN, ABAQUS (concepts and use)
  • Picture: the smooth L-bracket on the left, the same bracket covered in red quad elements on the right. The mesh approximates the smooth shape.
  • Why the topic needs it: we only know how to write physics equations for a simple shape. So we solve simple shapes and glue the answers together. This gluing is called assembly: the element stiffness tables are stacked into one global by adding contributions at shared nodes (Section 7 shows where comes from). This is literally the "finite" in Finite Element — finitely many pieces instead of infinitely smooth material.
Why quads for thin walls?
A thin bracket is basically a surface, so flat 4-node shell elements (CQUAD4) capture bending with far fewer DOF than solid bricks.

4. Force — "the push"

  • Picture: an arrow pushing on a node; its length is the force size.
  • Why bold? Same reason as — one force number per DOF, stacked into a column.

5. Stiffness — "how hard is it to push?"

Push a soft spring and it moves a lot; push a stiff spring and it barely budges. The ratio "force needed per unit of movement" is stiffness.

Figure — FEM software — NASTRAN, ABAQUS (concepts and use)
  • Picture: the straight red line . Its slope is . Steeper line = stiffer.
  • Why a matrix, not a number? A real structure has many DOF, and pushing node 1 also moves node 2 (they are connected). So we need a table where entry says "force at DOF caused by a unit move of DOF ." That table is the stiffness matrix, and the many-spring version of is exactly .

6. Boundary conditions — "pin it down or it floats away"

A structure floating free in space has nothing to push against: give it any force and the whole thing drifts, so no unique displacement exists. Mathematically, the raw is singular (its determinant is zero — you cannot invert it, and has infinitely many solutions).

Why does a free-floating model give no unique ?
With nothing pinned, rigid-body drift adds any constant motion for free, so is singular and solutions are non-unique.

7. Stress , Strain , and the matrices

Displacement tells you how points move. But to know if the metal cracks, you need what is happening inside it.

Now we need to connect the inside quantities (strain, stress) back to the nodal displacements we actually solve for. Two matrices do this, and both must be defined before use.

  • Why the colon in (parent's virtual-work integral)? Stress and strain are 2D tables (tensors), and "" just means "multiply matching entries and add them up" — it produces energy per volume.

The chain the parent uses is: move the nodes → get strain → get stress → check margin.


8. Young's Modulus and Density — "the material facts"

  • Why both and ? (via ) builds the stiffness ; builds the mass (next). Stiffness resists bending; mass resists acceleration. Vibration is a fight between the two.

9. Mass , Frequency , and Mode Shape — "how it vibrates"

Pluck a structure and it wobbles at special rates. FEM finds those.

  • Why it matters for the parent's Hz: if launch shaking contains energy at 850 Hz, the bracket resonates and can break. This is the bridge to Vibration and Modal Analysis and Launch Vehicle Loads.
Convert the parent's
rad/s.
Why do mode shapes exist only at certain frequencies?
Because needs singular, i.e. , which only specific satisfy.

10. Linear vs Nonlinear — "is fixed or moving?"

There are three flavours of nonlinearity the parent relies on, and each one is just " is no longer constant" for a different physical reason:

  • Material nonlinearity: past a certain stress the material stops springing back proportionally (yielding/plasticity, composite damage), so its stiffness drops as load grows.
  • Geometric nonlinearity: when parts bend or rotate a lot (deployable booms), the geometry itself changes enough that the stiffness depends on the current shape.
  • Contact nonlinearity: two surfaces touch or separate (bolted joints), and stiffness jumps the instant they come into contact.
Why iterate at all in nonlinear?
Because (hence ) depends on the answer , you cannot solve in one shot — you guess, measure the imbalance , correct, and repeat until .

11. Factor of Safety & Margin — "how much did we over-build?"

Connects to Factors of Safety and Stress Analysis and Margins — the decision stage after , are known.


How these foundations feed the topic

Nodes and Elements

Degrees of Freedom

Displacement vector u

Youngs modulus E

Material matrix D

Poisson ratio nu

Stiffness matrix K

Strain-displacement B

Boundary conditions

K u = F

External force F ext

Strain then Stress

Margins and Failure index

det K minus omega2 M = 0

Mass matrix M

Density rho

Nonlinear R and tangent Kt

NASTRAN linear and modal

ABAQUS nonlinear


Equipment checklist

Cover the right side and answer before revealing.

What does the arrow physically represent?
How far and in which direction a point of the structure moved from rest.
How many DOF does a shell-corner node with translation and rotation have?
6 (three translations, three rotations).
What is the difference between a node and an element?
A node is a labelled corner point; an element is a simple shape whose corners are nodes.
In , what is the slope of the line?
The stiffness .
Why is a matrix rather than a single number?
Because moving one DOF affects connected DOF, so we need a table of cross-couplings.
Name two mathematical properties of and why they matter.
Symmetric () and, once constrained, positive definite — giving a unique solution and fast Cholesky-type solves.
Why must you apply boundary conditions before solving ?
Without pinned DOF, is singular (rigid-body drift) and has no unique solution.
What is strain, and what are its units?
Fractional stretch (new length change ÷ original length); dimensionless.
What is Poisson's ratio , and a typical metal value?
Sideways contraction strain ÷ lengthwise stretch strain; about 0.33 for aluminium.
What do the shape functions do?
Blend the corner-node displacements to give the displacement at every interior point ().
What is ?
The short vector of just one element's nodal DOF, a slice of the global .
What does the strain–displacement matrix do?
Turns element nodal displacements into strain by differentiating the shape functions ().
What two material facts build , and what builds ?
and build (hence ); density builds .
Convert to ordinary frequency .
.
What condition makes nonzero mode shapes possible?
— the matrix must be singular.
Name the three flavours of nonlinearity.
Material (yielding/damage), geometric (large deformation), and contact (surfaces touching/separating).
Distinguish , , and .
= applied load (); = elastic push-back; = imbalance, driven to zero.
Why does Newton–Raphson use a derivative ?
The slope of the curved force line changes everywhere; the derivative gives the local slope to take the next step.
Write the Margin of Safety in terms of FoS.
.
A failure index of 1.0 means what?
The material is exactly at its breaking point.