Dynamic relaxation

<h2 id="main-equations-used">Main equations used</h2>
<p>Considering <a href="/facts/Newton%2527s_second_law_of_motion/XKzGTXK1">Newton's second law of motion</a> (force is mass multiplied by acceleration) in the 
  
    
      
        x
      
    
    {\displaystyle x}
  
 direction at the 
  
    
      
        i
      
    
    {\displaystyle i}
  
th node at time 
  
    
      
        t
      
    
    {\displaystyle t}
  
:
</p>

R
          
            i
            x
          
        
        (
        t
        )
        =
        
          M
          
            i
          
        
        
          A
          
            i
            x
          
        
        (
        t
        )
        
          
            
            
          
        
      
    
    {\displaystyle R_{ix}(t)=M_{i}A_{ix}(t){\frac {}{}}}

<p>Where:
</p>

R
      
    
    {\displaystyle R}
  
 is the residual force

M
      
    
    {\displaystyle M}
  
 is the nodal mass

A
      
    
    {\displaystyle A}
  
 is the nodal acceleration
<p>Note that fictitious nodal masses may be chosen to speed up the process of form-finding.
</p><p>The relationship between the speed 
  
    
      
        V
      
    
    {\displaystyle V}
  
, the geometry 
  
    
      
        X
      
    
    {\displaystyle X}
  
 and the residuals can be obtained by performing a double numerical integration of the acceleration (here in <a href="/facts/Central_difference/aY1txFhj">central finite difference</a> form<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a>), :
</p>

V
          
            i
            x
          
        
        
          (
          
            t
            +
            
              
                
                  Δ
                  t
                
                2
              
            
          
          )
        
        =
        
          V
          
            i
            x
          
        
        
          (
          
            t
            −
            
              
                
                  Δ
                  t
                
                2
              
            
          
          )
        
        +
        
          
            
              Δ
              t
            
            
              M
              
                i
              
            
          
        
        
          R
          
            i
            x
          
        
        (
        t
        )
      
    
    {\displaystyle V_{ix}\left(t+{\frac {\Delta t}{2}}\right)=V_{ix}\left(t-{\frac {\Delta t}{2}}\right)+{\frac {\Delta t}{M_{i}}}R_{ix}(t)}

X
          
            i
          
        
        (
        t
        +
        Δ
        t
        )
        =
        
          X
          
            i
          
        
        (
        t
        )
        +
        Δ
        t
        ×
        
          V
          
            i
            x
          
        
        
          (
          
            t
            +
            
              
                
                  Δ
                  t
                
                2
              
            
          
          )
        
      
    
    {\displaystyle X_{i}(t+\Delta t)=X_{i}(t)+\Delta t\times V_{ix}\left(t+{\frac {\Delta t}{2}}\right)}

<p>Where:
</p>

Δ
        t
      
    
    {\displaystyle \Delta t}
  
 is the time interval between two updates.
<p>By the principle of equilibrium of forces, the relationship between the residuals and the geometry can be obtained:
</p>

R
          
            i
            x
          
        
        (
        t
        +
        Δ
        t
        )
        =
        
          P
          
            i
            x
          
        
        (
        t
        +
        Δ
        t
        )
        +
        ∑
        
          
            
              
                T
                
                  m
                
              
              (
              t
              +
              Δ
              t
              )
            
            
              
                l
                
                  m
                
              
              (
              t
              +
              Δ
              t
              )
            
          
        
        ×
        (
        
          X
          
            j
          
        
        (
        t
        +
        Δ
        t
        )
        −
        
          X
          
            i
          
        
        (
        t
        +
        Δ
        t
        )
        )
      
    
    {\displaystyle R_{ix}(t+\Delta t)=P_{ix}(t+\Delta t)+\sum {\frac {T_{m}(t+\Delta t)}{l_{m}(t+\Delta t)}}\times (X_{j}(t+\Delta t)-X_{i}(t+\Delta t))}

<p>where:
</p>

P
      
    
    {\displaystyle P}
  
 is the applied load component

T
      
    
    {\displaystyle T}
  
 is the tension in link 
  
    
      
        m
      
    
    {\displaystyle m}
  
 between nodes 
  
    
      
        i
      
    
    {\displaystyle i}
  
 and 
  
    
      
        j
      
    
    {\displaystyle j}

l
      
    
    {\displaystyle l}
  
 is the length of the link.
<p>The sum must cover the forces in all the connections between the node and other nodes.
By repeating the use of the relationship between the residuals and the geometry, and the relationship between the geometry and the residual, the pseudo-dynamic process is simulated.
</p>
<h2 id="iteration-steps">Iteration steps</h2>
<p>1. Set the initial kinetic energy and all nodal velocity components to zero:
</p>

E
          
            k
          
        
        (
        t
        =
        0
        )
        =
        0
        
          
            
            
          
        
      
    
    {\displaystyle E_{k}(t=0)=0{\frac {}{}}}

V
          
            i
          
        
        (
        t
        =
        0
        )
        =
        0
        
          
            
            
          
        
      
    
    {\displaystyle V_{i}(t=0)=0{\frac {}{}}}

<p>2. Compute the geometry set and the applied load component:
</p>

X
          
            i
          
        
        (
        t
        =
        0
        )
        
          
            
            
          
        
      
    
    {\displaystyle X_{i}(t=0){\frac {}{}}}

P
          
            i
          
        
        (
        t
        =
        0
        )
        
          
            
            
          
        
      
    
    {\displaystyle P_{i}(t=0){\frac {}{}}}

<p>3. Compute the residual:
</p>

T
          
            m
          
        
        (
        t
        )
        
          
            
            
          
        
      
    
    {\displaystyle T_{m}(t){\frac {}{}}}

R
          
            i
          
        
        (
        t
        )
        
          
            
            
          
        
      
    
    {\displaystyle R_{i}(t){\frac {}{}}}

<p>4. Reset the residuals of constrained nodes to zero
</p><p>5. Update velocity and coordinates:
</p>

V
          
            i
          
        
        (
        t
        +
        
          
            
              Δ
              t
            
            2
          
        
        )
        
          
            
            
          
        
      
    
    {\displaystyle V_{i}(t+{\frac {\Delta t}{2}}){\frac {}{}}}

X
          
            i
          
        
        (
        t
        +
        Δ
        t
        )
        
          
            
            
          
        
      
    
    {\displaystyle X_{i}(t+\Delta t){\frac {}{}}}

<p>6. Return to step 3 until the structure is in static <a href="/facts/Mechanical_equilibrium/N2XKaK6s">equilibrium</a>
</p>
<h2 id="damping">Damping</h2>
<p>It is possible to make dynamic relaxation more computationally efficient (reducing the number of iterations) by using damping.<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a>
There are two methods of damping:
</p>
<ul><li>Viscous damping, which assumes that connection between the nodes has a viscous force component.</li>
<li>Kinetic energy damping, where the coordinates at peak kinetic energy are calculated (the equilibrium position), then updates the geometry to this position and resets the velocity to zero.</li></ul>
<p>The advantage of viscous damping is that it represents the reality of a cable with viscous properties. Moreover, it is easy to realize because the speed is already computed.
The kinetic energy damping is an artificial damping which is not a real effect, but offers a drastic reduction in the number of iterations required to find a solution. However, there is a computational penalty in that the kinetic energy and peak location must be calculated, after which the geometry has to be updated to this position.
</p>
<h2 id="see-also">See also</h2>
<ul><li><a href="/facts/Tensile_structure/MUlT2xgs">Tensile structures</a></li>
<li><a href="/facts/Optimization_(mathematics)/oRn8Iv5I">Optimization (mathematics)</a></li></ul>
<h2 id="further-reading">Further reading</h2>
<ul><li>A S Day, <i>An introduction to dynamic relaxation.</i> The Engineer 1965, 219:218–221</li>
<li>H.A. BUCHHOLDT, <i>An introduction to cable roof structures</i>, 2nd ed, London, Telford, 1999</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1"><p>W. J. Lewis, Tension Structures: Form and behaviour, London, Telford, 2003 <a href="/wiki/Wanda_Lewis" target="_blank">/wiki/Wanda_Lewis</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li>
<li id="fn:2"><p>D S WAKEFIELD, Engineering analysis of tension structures: theory and practice, Bath, Tensys Limited, 1999 <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li>
<li id="fn:3"><p>W. J. Lewis, Tension Structures: Form and behaviour, London, Telford, 2003 <a href="/wiki/Wanda_Lewis" target="_blank">/wiki/Wanda_Lewis</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li>
</ol>

Dynamic relaxation open-in-new

Dynamic relaxation