Arc-Length Control: Difference between revisions
(Created page with 'THEORY: If we write the governing finite element equation at <math>t + \Delta t\!</math>as: :<math> R(U_{t+\Delta t}, \lambda_{t+\Delta t}) = \lambda_{t+\Delta t} F^{ext} - F(U...') |
No edit summary |
||
| Line 1: | Line 1: | ||
{{CommandManualMenu}} | |||
This command is used to construct an ArcLength integrator object. In an analysis step with ArcLength we seek to determine the time step that will result in our constraint equation being satisfied. | |||
{| | |||
| style="background:yellow; color:black; width:800px" | '''integrator ArcLength $s $alpha''' | |||
|} | |||
---- | |||
{| | |||
| style="width:150px" | '''$s''' || <math>s</math> the arcLength. | |||
|- | |||
| '''$alpha''' || <math>\alpha</math> a scaling factor on the reference loads. | |||
|} | |||
---- | |||
EXAMPLE: | |||
integrator ArcLength 1.0 0.1; | |||
---- | |||
THEORY: | THEORY: | ||
| Line 12: | Line 38: | ||
This equation represents n equations in <math> n+1</math> unknowns, and so an additional equation is needed to solve the equation. For displacement control, we introduce a new constraint equation in which in each analysis step we set to ensure that the displacement increment for the degree-of-freedom <math>\text{dof}</math> at the specified node is: | This equation represents n equations in <math> n+1</math> unknowns, and so an additional equation is needed to solve the equation. For displacement control, we introduce a new constraint equation in which in each analysis step we set to ensure that the displacement increment for the degree-of-freedom <math>\text{dof}</math> at the specified node is: | ||
:<math> | :<math> \!</math> | ||
Latest revision as of 20:54, 12 March 2010
- Command_Manual
- Tcl Commands
- Modeling_Commands
- model
- uniaxialMaterial
- ndMaterial
- frictionModel
- section
- geometricTransf
- element
- node
- sp commands
- mp commands
- timeSeries
- pattern
- mass
- block commands
- region
- rayleigh
- Analysis Commands
- Output Commands
- Misc Commands
- DataBase Commands
This command is used to construct an ArcLength integrator object. In an analysis step with ArcLength we seek to determine the time step that will result in our constraint equation being satisfied.
| integrator ArcLength $s $alpha |
| $s | <math>s</math> the arcLength. |
| $alpha | <math>\alpha</math> a scaling factor on the reference loads. |
EXAMPLE:
integrator ArcLength 1.0 0.1;
THEORY:
If we write the governing finite element equation at <math>t + \Delta t\!</math>as:
- <math> R(U_{t+\Delta t}, \lambda_{t+\Delta t}) = \lambda_{t+\Delta t} F^{ext} - F(U_{t+\Delta t}) \!</math>
where <math>F(U_{t+\Delta t})\!</math> are the internal forces which are a function of the displacements <math>U_{t+\Delta t}\!</math>, <math>F^{ext}\!</math> is the set of reference loads and <math>\lambda\!</math> is the load multiplier. Linearizing the equation results in:
- <math> K_{t+\Delta t}^{*i} \Delta U_{t+\Delta t}^{i+1} = \left ( \lambda^i_{t+\Delta t} + \Delta \lambda^i \right ) F^{ext} - F(U_{t+\Delta t})</math>
This equation represents n equations in <math> n+1</math> unknowns, and so an additional equation is needed to solve the equation. For displacement control, we introduce a new constraint equation in which in each analysis step we set to ensure that the displacement increment for the degree-of-freedom <math>\text{dof}</math> at the specified node is:
- <math> \!</math>
MORE TO COME:
Code Developed by: fmk