Central Difference

• Command_Manual
  • Tcl Commands
  • Modeling_Commands
    • model
    • uniaxialMaterial
    • ndMaterial
    • frictionModel
    • section
    • geometricTransf
    • element
    • node
    • sp commands
      • fix
      • fixX
      • fixY
      • fixZ
    • mp commands
      • equalDOF
      • rigidDiaphragm
      • rigidLink
    • timeSeries
    • pattern
    • mass
    • block commands
      • block2D
      • block3D
    • region
    • rayleigh
  • Analysis Commands
    • constraints Command
    • numberer Command
    • system Command
    • test Command
    • algorithm Command
    • integrator Command
    • analysis Command
    • eigen Command
    • analyze Command
  • Output Commands
    • recorder
    • record
    • print
    • printA
    • logFile
    • RealTime_Output_Commands
  • Misc Commands
    • eleResponse Command
    • getTime Command
    • getEleTags Command
    • getNodeTags
    • loadConst Command
    • nodeCoord Command
    • nodeDisp Command
    • nodeVel Command
    • nodeAccel Command
    • nodeEigenvector Command
    • nodeBounds
    • print Command
    • remove Command
    • reset Command
    • testIter
    • testNorms
    • wipe Command
    • exit Command
  • DataBase Commands
    • Template:DataBaseCommandsMenu

---

This command is used to construct a Central Difference integrator object.

integrator CentralDifference

---

EXAMPLE:

integrator CentralDifference

---

NOTES:

1. The calculation of <math>U_{t+\Delta t}</math>, as shown below, is based on using the equilibrium equation at time t. For this reason the method is called an explicit integration method.
2. If there is no rayleigh damping and the C matrix is 0, for a diagonal mass matrix a diagonal solver may and should be used.
3. For stability, <math>\frac{\Delta t}{T_n} < \frac{1}{\pi}</math>

---

REFERENCES

---

THEORY:

The Central difference approximations for velocity and acceleration:

<math> v_n = \frac{d_{n+1} - d_{n-1}}{2 \Delta t}</math>

<math> a_n = \frac{d_{n+1} - 2 d_n + d_{n-1}}{\Delta t^2}</math>

In the Central Difference method we determine the displacement solution at time <math>t+\delta t</math> by considering the the eqilibrium equation for the finite element system in motion at time t:

<math>M \ddot U_t + C \dot U_t + K U_t = R_t</math>

which when using the above two expressions of becomes:

<math> \left ( \frac{1}{\Delta t^2} M + \frac{1}{2 \Delta t} C \right ) U_{t+\Delta t} = R_t - \left (K - \frac{2}{\Delta t^2}M \right )U_t - \left (\frac{1}{\Delta t^2}M - \frac{1}{2 \Delta t} C \right) U_{t-\Delta t} </math>

---

Code Developed by: fmk