convergence issue

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
lvy
Posts: 26
Joined: Tue Dec 14, 2004 8:22 pm
Location: Shanghai, China

convergence issue

Post by lvy »

I built a bridge model with fiber element and applied uniform excitation. I set g=4.8 and get convergent result. When I increase the value of g, from 4.9 to 20.0, I couldn't get the convergent result. The program stoped earlier and earlier because of the increcement of g. The script and data file are as follows.

Could you like give me some suggestion? Thank for any answer!

Code: Select all

# Hollow High Pier
# ------------------------
#  Reinforced concrete bridge with 30m tall pier 
# -----------------
#  Nonlinear beam-column elements
# 
# Units: kN, m, sec
#
# Written: Liang Zhiyao
#Tongji University
# Date: May 19, 2005


# ------------------------------
# Start of model generation
# ------------------------------

# Create ModelBuilder (with three-dimensions and 6 DOF/node)
model basic -ndm 3 -ndf 6

#record all information
logFile lvy.log

# Create nodes
# ------------

# Set parameters for overall model geometry

# Create nodes
# 1#pier
#    tag        X        Y       Z 
node   1        40.0    0.00    0.0 
node   2        40.0    0.00    0.5 
node   3        40.0    0.00    1.0 
node   4        40.0    0.00    1.5 
node   5        40.0    0.00    2.0 
node   6        40.0    0.00    2.5 
node   7        40.0    0.00    3.0 
node   8        40.0    0.00    3.5 
node   9        40.0    0.00    4.0 
node   10       40.0    0.00    5.0 
node   11       40.0    0.00    6.0 
node   12       40.0    0.00    8.0 
node   13       40.0    0.00    10.0
node   14       40.0    0.00    12.0
node   15       40.0    0.00    14.0 
node   16       40.0    0.00    16.0 
node   17       40.0    0.00    18.0 
node   18       40.0    0.00    20.0 
node   19       40.0    0.00    22.0 
node   20       40.0    0.00    24.0 
node   21       40.0    0.00    26.0 
node   22       40.0    0.00    28.0 
node   23       40.0    0.00    30.0 

# 2#pier                                
node   24       80.0    0.00    0.0 
node   25       80.0    0.00    0.5 
node   26       80.0    0.00    1.0 
node   27       80.0    0.00    1.5 
node   28       80.0    0.00    2.0 
node   29       80.0    0.00    2.5 
node   30       80.0    0.00    3.0 
node   31       80.0    0.00    3.5 
node   32       80.0    0.00    4.0 
node   33       80.0    0.00    5.0 
node   34       80.0    0.00    6.0 
node   35       80.0    0.00    8.0 
node   36       80.0    0.00    10.0 
node   37       80.0    0.00    12.0 
node   38       80.0    0.00    14.0 
node   39       80.0    0.00    16.0 
node   40       80.0    0.00    18.0 
node   41       80.0    0.00    20.0 
node   42       80.0    0.00    22.0 
node   43       80.0    0.00    24.0 
node   44       80.0    0.00    26.0 
node   45       80.0    0.00    28.0 
node   46       80.0    0.00    30.0 

# 3#pier                                
node   47       120.0   0.00    0.0 
node   48       120.0   0.00    0.5 
node   49       120.0   0.00    1.0 
node   50       120.0   0.00    1.5 
node   51       120.0   0.00    2.0 
node   52       120.0   0.00    2.5 
node   53       120.0   0.00    3.0 
node   54       120.0   0.00    3.5 
node   55       120.0   0.00    4.0 
node   56       120.0   0.00    5.0 
node   57       120.0   0.00    6.0 
node   58       120.0   0.00    8.0 
node   59       120.0   0.00    10.0 
node   60       120.0   0.00    12.0 
node   61       120.0   0.00    14.0 
node   62       120.0   0.00    16.0 
node   63       120.0   0.00    18.0 
node   64       120.0   0.00    20.0 
node   65       120.0   0.00    22.0 
node   66       120.0   0.00    24.0 
node   67       120.0   0.00    26.0 
node   68       120.0   0.00    28.0 
node   69       120.0   0.00    30.0 

node   70        0.0    0.00    30.0
node   71       10.0    0.00    30.0
node   72       20.0    0.00    30.0
node   73       30.0    0.00    30.0
node   74       50.0    0.00    30.0
node   75       60.0    0.00    30.0
node   76       70.0    0.00    30.0
node   77       90.0    0.00    30.0
node   78      100.0    0.00    30.0
node   79      110.0    0.00    30.0
node   80      130.0    0.00    30.0
node   81      140.0    0.00    30.0
node   82      150.0    0.00    30.0
node   83      160.0    0.00    30.0

#     tag     MX     MY     MZ     RX     RY     RZ
mass   2    5.6    5.6     0.0     0.0    0.0    0.0
mass   3    5.6    5.6     0.0     0.0    0.0    0.0
mass   4    5.6    5.6     0.0     0.0    0.0    0.0
mass   5    5.6    5.6     0.0     0.0    0.0    0.0
mass   6    5.6    5.6     0.0     0.0    0.0    0.0
mass   7    5.6    5.6     0.0     0.0    0.0    0.0
mass   8    5.6    5.6     0.0     0.0    0.0    0.0
mass   9    8.4    8.4     0.0     0.0    0.0    0.0
mass   10   11.2   11.2    0.0     0.0    0.0    0.0
mass   11   16.8   16.8    0.0     0.0    0.0    0.0
mass   12   22.4   22.4    0.0     0.0    0.0    0.0
mass   13   22.4   22.4    0.0     0.0    0.0    0.0
mass   14   22.4   22.4    0.0     0.0    0.0    0.0
mass   15   22.4   22.4    0.0     0.0    0.0    0.0
mass   16   22.4   22.4    0.0     0.0    0.0    0.0
mass   17   22.4   22.4    0.0     0.0    0.0    0.0
mass   18   22.4   22.4    0.0     0.0    0.0    0.0
mass   19   22.4   22.4    0.0     0.0    0.0    0.0
mass   20   22.4   22.4    0.0     0.0    0.0    0.0
mass   21   22.4   22.4    0.0     0.0    0.0    0.0
mass   22   22.4   22.4    0.0     0.0    0.0    0.0
mass   23   198.7  198.7   0.0     0.0    0.0    0.0

mass   25   5.6    5.6     0.0     0.0    0.0    0.0
mass   26   5.6    5.6     0.0     0.0    0.0    0.0
mass   27   5.6    5.6     0.0     0.0    0.0    0.0
mass   28   5.6    5.6     0.0     0.0    0.0    0.0
mass   29   5.6    5.6     0.0     0.0    0.0    0.0
mass   30   5.6    5.6     0.0     0.0    0.0    0.0
mass   31   5.6    5.6     0.0     0.0    0.0    0.0
mass   32   8.4    8.4     0.0     0.0    0.0    0.0
mass   33   11.2   11.2    0.0     0.0    0.0    0.0
mass   34   16.8   16.8    0.0     0.0    0.0    0.0
mass   35   22.4   22.4    0.0     0.0    0.0    0.0
mass   36   22.4   22.4    0.0     0.0    0.0    0.0
mass   37   22.4   22.4    0.0     0.0    0.0    0.0
mass   38   22.4   22.4    0.0     0.0    0.0    0.0
mass   39   22.4   22.4    0.0     0.0    0.0    0.0
mass   40   22.4   22.4    0.0     0.0    0.0    0.0
mass   41   22.4   22.4    0.0     0.0    0.0    0.0
mass   42   22.4   22.4    0.0     0.0    0.0    0.0
mass   43   22.4   22.4    0.0     0.0    0.0    0.0
mass   44   22.4   22.4    0.0     0.0    0.0    0.0
mass   45   22.4   22.4    0.0     0.0    0.0    0.0
mass   46   198.7  198.7   0.0     0.0    0.0    0.0

mass   48   5.6    5.6     0.0     0.0    0.0    0.0
mass   49   5.6    5.6     0.0     0.0    0.0    0.0
mass   50   5.6    5.6     0.0     0.0    0.0    0.0
mass   51   5.6    5.6     0.0     0.0    0.0    0.0
mass   52   5.6    5.6     0.0     0.0    0.0    0.0
mass   53   5.6    5.6     0.0     0.0    0.0    0.0
mass   54   5.6    5.6     0.0     0.0    0.0    0.0
mass   55   8.4    8.4     0.0     0.0    0.0    0.0
mass   56   11.2   11.2    0.0     0.0    0.0    0.0
mass   57   16.8   16.8    0.0     0.0    0.0    0.0
mass   58   22.4   22.4    0.0     0.0    0.0    0.0
mass   59   22.4   22.4    0.0     0.0    0.0    0.0
mass   60   22.4   22.4    0.0     0.0    0.0    0.0
mass   61   22.4   22.4    0.0     0.0    0.0    0.0
mass   62   22.4   22.4    0.0     0.0    0.0    0.0
mass   63   22.4   22.4    0.0     0.0    0.0    0.0
mass   64   22.4   22.4    0.0     0.0    0.0    0.0
mass   65   22.4   22.4    0.0     0.0    0.0    0.0
mass   66   22.4   22.4    0.0     0.0    0.0    0.0
mass   67   22.4   22.4    0.0     0.0    0.0    0.0
mass   68   22.4   22.4    0.0     0.0    0.0    0.0
mass   69   198.7  198.7   0.0     0.0    0.0    0.0

mass   70   93.75  93.75  0.0     0.0    0.0    0.0
mass   71   187.5  187.5  0.0     0.0    0.0    0.0
mass   72   187.5  187.5  0.0     0.0    0.0    0.0
mass   73   187.5  187.5  0.0     0.0    0.0    0.0
mass   74   187.5  187.5  0.0     0.0    0.0    0.0
mass   75   187.5  187.5  0.0     0.0    0.0    0.0
mass   76   187.5  187.5  0.0     0.0    0.0    0.0
mass   77   187.5  187.5  0.0     0.0    0.0    0.0
mass   78   187.5  187.5  0.0     0.0    0.0    0.0
mass   79   187.5  187.5  0.0     0.0    0.0    0.0
mass   80   187.5  187.5  0.0     0.0    0.0    0.0
mass   81   187.5  187.5  0.0     0.0    0.0    0.0
mass   82   187.5  187.5  0.0     0.0    0.0    0.0
mass   83   93.75  93.75  0.0     0.0    0.0    0.0

# Fix supports at base of columns
#    tag   DX   DY   DZ    RX   RY   RZ
fix   1     1    1    1    1    1    1
fix   24    1    1    1    1    1    1
fix   47    1    1    1    1    1    1
fix   70    0    1    1    1    0    1
fix   85    0    1    1    1    0    1


# Define materials for nonlinear columns
# ------------------------------------------
# CONCRETE                  tag   f'c        ec0   f'cu        ecu
# Core concrete (confined)
uniaxialMaterial Concrete01 1  -44540.0  -0.0051  -28500.0   -0.0254
# Cover concrete (unconfined)
uniaxialMaterial Concrete01 2  -34000.0  -0.002   0.0     -0.006

# STEEL
# Reinforcing steel 
set fy 340000.0;      # Yield stress
set E 200000000.0;    # Young's modulus
#                        tag  fy E0    b
uniaxialMaterial Steel01  3  $fy $E 0.001

# Define cross-section for nonlinear columns
# ------------------------------------------

# set some paramaters
set y  1.55
set y2 1.25
set z  1.55
set z2 1.25
set cover  0.05
set As    0.000490625;     # area of 25mm bars


section Fiber 1 {

    # Create the concrete core fibers
    patch rect 1 2  10 [expr -$y] [expr -$z] [expr -$y2] [expr $z]
    patch rect 1 2  10 [expr $y2] [expr -$z] [expr $y] [expr $z]
    patch rect 1 10  2 [expr -$y2] [expr $z2] [expr $y2] [expr $z]
    patch rect 1 10  2 [expr -$y2] [expr -$z] [expr $y2] [expr -$z2]

    # Create the concrete cover fibers (top, bottom, left, right)
    patch rect 2 10 1  [expr -$y-$cover] [expr $z] [expr $y+$cover] [expr $z+$cover]
    patch rect 2 10 1  [expr -$y-$cover] [expr -$z-$cover] [expr $y+$cover] [expr -$z]
    patch rect 2 1 10  [expr -$y-$cover] [expr -$z] [expr -$y] [expr $z]
    patch rect 2 1 10  [expr $y] [expr -$z] [expr $y+$cover] [expr $z]
    
    patch rect 2 10 1  [expr -$y2] [expr $z2-$cover] [expr $y2] [expr $z2]
    patch rect 2 10 1  [expr -$y2] [expr -$z2] [expr $y2] [expr -$z2+$cover]
    patch rect 2 1 10  [expr -$y2] [expr -$z2+$cover] [expr -$y2+$cover] [expr $z2-$cover]
    patch rect 2 1 10  [expr $y2-$cover] [expr -$z2+$cover] [expr $y2] [expr $z2-$cover]

    # Create the reinforcing fibers (left, middle, right)
    layer straight 3 21 $As [expr -$y] [expr $z] [expr -$y] [expr -$z]
    layer straight 3 21 $As [expr $y] [expr $z] [expr $y] [expr -$z]
    layer straight 3 10 $As [expr -$y2] [expr $z2] [expr -$y2] [expr -$z2]
    layer straight 3 10 $As [expr $y2] [expr $z2] [expr $y2] [expr -$z2]
    layer straight 3 19 $As [expr -$y+0.12] [expr $z] [expr $y-0.12] [expr $z]
    layer straight 3 19 $As [expr -$y+0.12] [expr -$z] [expr $y-0.12] [expr -$z]
    layer straight 3  8 $As [expr -$y2+0.2] [expr -$z2] [expr $y2-0.2] [expr -$z2]
    layer straight 3  8 $As [expr -$y2+0.2] [expr $z2] [expr $y2-0.2] [expr $z2]
}    

#Linear elastic torsion
uniaxialMaterial Elastic  101  [expr 11.6942*13000000.0]

#Attach torsion to the RC section
section Aggregator 2 101 T -section 1

# Define column elements
# ----------------------

# Geometry of column elements
#                tag 
geomTransf Linear 1 0 1 0
geomTransf Linear 2 0 0 1  

# Number of integration points along length of element
set np 5

# Create the coulumns using Beam-column elements
#                           tag ndI ndJ  nsecs secID transfTag  mass
element nonlinearBeamColumn  1   1   2    $np    2       1 
element nonlinearBeamColumn  2   2   3    $np    2       1 
element nonlinearBeamColumn  3   3   4    $np    2       1 
element nonlinearBeamColumn  4   4   5    $np    2       1 
element nonlinearBeamColumn  5   5   6    $np    2       1 
element nonlinearBeamColumn  6   6   7    $np    2       1
element nonlinearBeamColumn  7   7   8    $np    2       1  
element nonlinearBeamColumn  8   8   9    $np    2       1  
element nonlinearBeamColumn  9   9   10   $np    2       1   
element nonlinearBeamColumn  10  10  11   $np    2       1  
element nonlinearBeamColumn  11  11  12   $np    2       1    
element nonlinearBeamColumn  12  12  13   $np    2       1  
element nonlinearBeamColumn  13  13  14   $np    2       1   
element nonlinearBeamColumn  14  14  15   $np    2       1   
element nonlinearBeamColumn  15  15  16   $np    2       1    
element nonlinearBeamColumn  16  16  17   $np    2       1   
element nonlinearBeamColumn  17  17  18   $np    2       1   
element nonlinearBeamColumn  18  18  19   $np    2       1  
element nonlinearBeamColumn  19  19  20   $np    2       1   
element nonlinearBeamColumn  20  20  21   $np    2       1   
element nonlinearBeamColumn  21  21  22   $np    2       1    
element nonlinearBeamColumn  22  22  23   $np    2       1   

element nonlinearBeamColumn  23  24  25   $np    2       1 
element nonlinearBeamColumn  24  25  26   $np    2       1    
element nonlinearBeamColumn  25  26  27   $np    2       1   
element nonlinearBeamColumn  26  27  28   $np    2       1   
element nonlinearBeamColumn  27  28  29   $np    2       1  
element nonlinearBeamColumn  28  29  30   $np    2       1   
element nonlinearBeamColumn  29  30  31   $np    2       1   
element nonlinearBeamColumn  30  31  32   $np    2       1    
element nonlinearBeamColumn  31  32  33   $np    2       1 
element nonlinearBeamColumn  32  33  34   $np    2       1   
element nonlinearBeamColumn  33  34  35   $np    2       1 
element nonlinearBeamColumn  34  35  36   $np    2       1    
element nonlinearBeamColumn  35  36  37   $np    2       1   
element nonlinearBeamColumn  36  37  38   $np    2       1   
element nonlinearBeamColumn  37  38  39   $np    2       1  
element nonlinearBeamColumn  38  39  40   $np    2       1   
element nonlinearBeamColumn  39  40  41   $np    2       1   
element nonlinearBeamColumn  40  41  42   $np    2       1    
element nonlinearBeamColumn  41  42  43   $np    2       1
element nonlinearBeamColumn  42  43  44   $np    2       1     
element nonlinearBeamColumn  43  44  45   $np    2       1 
element nonlinearBeamColumn  44  45  46   $np    2       1    

element nonlinearBeamColumn  45  47  48   $np    2       1   
element nonlinearBeamColumn  46  48  49   $np    2       1  
element nonlinearBeamColumn  47  49  50   $np    2       1   
element nonlinearBeamColumn  48  50  51   $np    2       1   
element nonlinearBeamColumn  49  51  52   $np    2       1    
element nonlinearBeamColumn  50  52  53   $np    2       1 
element nonlinearBeamColumn  51  53  54   $np    2       1   
element nonlinearBeamColumn  52  54  55   $np    2       1 
element nonlinearBeamColumn  53  55  56   $np    2       1    
element nonlinearBeamColumn  54  56  57   $np    2       1   
element nonlinearBeamColumn  55  57  58   $np    2       1   
element nonlinearBeamColumn  56  58  59   $np    2       1  
element nonlinearBeamColumn  57  59  60   $np    2       1   
element nonlinearBeamColumn  58  60  61   $np    2       1   
element nonlinearBeamColumn  59  61  62   $np    2       1    
element nonlinearBeamColumn  60  62  63   $np    2       1
element nonlinearBeamColumn  61  63  64   $np    2       1     
element nonlinearBeamColumn  62  64  65   $np    2       1 
element nonlinearBeamColumn  63  65  66   $np    2       1 
element nonlinearBeamColumn  64  65  67   $np    2       1    
element nonlinearBeamColumn  65  67  68   $np    2       1
element nonlinearBeamColumn  66  68  69   $np    2       1 

#                         $eleTag $iNode $jNode $A  $E $G  $J  $Iy $Iz $transfTag
element elasticBeamColumn 67  70 71 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 68  71 72 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 69  72 73 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 70  73 23 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 71  23 74 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 72  74 75 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 73  75 76 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 74  76 46 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 75  46 77 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 76  77 78 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 77  78 79 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 78  79 69 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 79  69 80 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 80  80 81 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 81  81 82 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2
element elasticBeamColumn 82  82 83 7.5 34500000.0 14375000.0 6.034 2.0 107.5  2

# initialize in case we need to do an initial stiffness iteration
initialize

# ------------------------------
# End of model generation
# ------------------------------

set f [open mode.out a+]
#statement line
puts $f "    Period           Frequency"
puts $f "    T/sec              Hz" 
set lambda [eigen 10]
for {set i 0} {$i<10} {incr i 1} {
set lambdaa [lindex $lambda $i]
set omega [expr pow($lambdaa,0.5)]
set Tperiod [expr 2*3.14/$omega]
set Freq [expr 1/$Tperiod]
# period (sec.) and frequency 
puts $f "$Tperiod\t   $Freq"
}
close $f

for {set j 1} {$j<11} {incr j 1} {
recorder Node -file modevector$j.out -time -nodeRange 1 83 -dof 1 "eigen $j"
}

#***********************************************************************************
#***********************************************************************************
# Start of Model Generation & Initial Gravity Analysis
# ----------------------------------------------------

# ----------------------------------------------------
# Start of additional modeling for dynamic loads
# ----------------------------------------------------

# Define dynamic loads
# --------------------
# Set some parameters
set inputFile elcentro.txt
set dt 0.01
set g 4.8 

# Set time series to be passed to uniform excitation
set accelSeries "Path -filePath $inputFile -dt $dt -factor $g"

# Create UniformExcitation load pattern
#                        tag dir 
pattern UniformExcitation 1 1 -accel $accelSeries

# set the rayleigh damping factors for nodes & elements
rayleigh 0.0 0.0 0.0 0.0
# ----------------------------------------------------
# End of additional modeling for dynamic loads
# ----------------------------------------------------

# Create the convergence test, the norm of the residual with a tolerance of 
# 1e-12 and a max number of iterations of 10
test EnergyIncr 1.0e-12 15 0 

# Create the solution algorithm, a Newton-Raphson algorithm
algorithm Newton

# Create the integration scheme, Newmark with gamma = 0.5 and beta =  0.25
integrator Newmark  0.5  0.25 

# Create the system of equation, a banded general storage scheme
system UmfPack

# Create the constraint handler, a plain handler as homogeneous boundary conditions
constraints Transformation

# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer RCM

# Create the analysis object
analysis Transient
# ---------------------------------------------------------
# End of modifications to analysis for transient analysis
# ---------------------------------------------------------


# ------------------------------
# Start of recorder generation
# ------------------------------
# Create a recorder which writes to Node.out and prints
# the current load factor (pseudo-time) and dof 1 displacements at node 38
recorder Node -file PierNode-1.out -time -nodeRange 1  23 -dof 1 disp
recorder Node -file PierNode-2.out -time -nodeRange 24 46 -dof 1 disp
recorder Node -file PierNode-3.out -time -nodeRange 47 69 -dof 1 disp

recorder Element -file PierElement-1.out -time -eleRange 1  22 section 1 deformation
recorder Element -file PierElement-2.out -time -eleRange 23 44 section 1 deformation
recorder Element -file PierElement-3.out -time -eleRange 45 66 section 1 deformation

# a window to plot the nodal displacements versus load for node 38
recorder plot PierNode-2.out Node_46_Xdisp 10 100 800 300 -columns 1 24 -dT 0.1
recorder display DynamicColumn 10 410 500 300 -wipe 
prp 100 100 120.5 
vup 0 0 1 
display 1 1 1
# ------------------------------
# End of recorder generation
# ------------------------------

# set some variables
set tFinal [expr 2000*0.01]
set tCurrent [getTime]
set ok 0

# Perform the transient analysis
while {$ok == 0 && $tCurrent < $tFinal} {
    
    set ok [analyze 1 0.01]
    
    # if the analysis fails try initial tangent iteration
    if {$ok != 0} {
	puts "regular newton failed .. lets try an initail stiffness for this step"
	test NormDispIncr 1.0e-12  100 0
	algorithm ModifiedNewton -initial
	set ok [analyze 1 0.01]
	if {$ok == 0} {puts "that worked .. back to regular newton"}
	test NormDispIncr 1.0e-12  10 
	algorithm Newton
    }
    
    set tCurrent [getTime]
}

# Print a message to indicate if analysis succesfull or not
if {$ok == 0} {
   puts "Transient analysis completed SUCCESSFULLY";
} else {
   puts "Transient analysis completed FAILED";    
}

# output the eigenvalue analysis
set f2 [open modevector.out a+]
for {set j 1} {$j<11} {incr j 1} {
set f3 [open modevector$j.out r]
gets $f3 line
puts $f2 $line
close $f3
}

remove recorders
for {set j 1} {$j<11} {incr j 1} {
file delete -force modevector$j.out
}

Code: Select all

0.0
-0.006403
-0.006029
0.00053
0.007743
0.00689
0.006877
0.006705
0.006473
0.006274
0.006153
0.006134
0.0062
0.006279
0.006305
0.006276
0.00623
0.006201
0.006189
0.006183
0.006171
0.006149
0.006117
0.006082
0.006058
0.006054
0.006071
0.006099
0.006132
0.006176
0.006223
0.006247
0.00624
0.006221
0.006212
0.006215
0.006221
0.006217
0.006186
0.006116
0.006012
0.005923
0.005934
0.006085
0.006319
0.006557
0.006773
0.00695
0.006975
0.006701
0.006155
0.005509
0.004896
0.004341
0.003941
0.004019
0.004751
0.005856
0.006661
0.006626
0.005996
0.005489
0.005539
0.006114
0.006905
0.007611
0.008139
0.008484
0.008457
0.007841
0.006888
0.006233
0.006218
0.006617
0.006992
0.007113
0.006989
0.006582
0.005808
0.004758
0.003621
0.002485
0.001265
-0.000088
-0.001298
-0.00201
-0.002243
-0.002483
-0.003321
-0.004912
-0.006945
-0.008873
-0.009937
-0.009639
-0.008336
-0.00689
-0.005847
-0.005167
-0.00448
-0.003587
-0.002955
-0.003375
-0.004967
-0.00725
-0.009874
-0.0125
-0.0149
-0.0169
-0.0189
-0.021
-0.023
-0.0238
-0.0232
-0.022
-0.0206
-0.0186
-0.0159
-0.0135
-0.0124
-0.0112
-0.006375
0.003818
0.0175
0.031
0.0398
0.0404
0.0343
0.0268
0.0217
0.0191
0.017
0.0133
0.007953
0.002115
-0.00326
-0.007974
-0.0124
-0.0171
-0.0222
-0.0279
-0.0344
-0.042
-0.0501
-0.0583
-0.0663
-0.0743
-0.0826
-0.0891
-0.0897
-0.0829
-0.0724
-0.0624
-0.055
-0.0516
-0.0528
-0.0565
-0.0593
-0.0591
-0.0565
-0.054
-0.054
-0.058
-0.0651
-0.0729
-0.0791
-0.0764
-0.0499
0.004302
0.0726
0.1337
0.1667
0.1657
0.1468
0.1319
0.1298
0.1356
0.14
0.1381
0.1324
0.1284
0.1287
0.1318
0.1354
0.138
0.1396
0.1391
0.134
0.1238
0.1114
0.0996
0.0893
0.0796
0.0692
0.057
0.0422
0.0251
0.007266
-0.0101
-0.0269
-0.0436
-0.0607
-0.0783
-0.0965
-0.115
-0.1335
-0.1519
-0.1703
-0.1887
-0.2073
-0.2262
-0.2443
-0.2596
-0.2716
-0.2822
-0.2934
-0.3054
-0.3129
-0.3082
-0.2902
-0.2624
-0.2215
-0.1455
-0.00573
0.165
0.2712
0.266
0.1963
0.1343
0.1131
0.121
0.129
0.1185
0.0893
0.0528
0.0211
-0.000354
-0.0142
-0.0269
-0.0429
-0.0623
-0.0832
-0.1034
-0.122
-0.1397
-0.159
-0.1726
-0.1364
-0.002874
0.1817
0.2956
0.2819
0.1939
0.1148
0.0837
0.0877
0.0936
0.0793
0.0432
-0.000885
-0.025
0.002765
0.0847
0.1789
0.2299
0.219
0.1763
0.1415
0.1313
0.1374
0.1431
0.1381
0.1228
0.1037
0.0867
0.0738
0.0635
0.0535
0.042
0.0293
0.0161
0.003169
-0.009445
-0.0203
-0.023
-0.0112
0.012
0.0374
0.0587
0.0761
0.0901
0.0958
0.0896
0.0761
0.0628
0.0532
0.0459
0.038
0.0272
0.0134
-0.001676
-0.0164
-0.03
-0.0429
-0.0561
-0.0664
-0.0609
-0.0288
0.0185
0.0509
0.0489
0.0222
-0.006025
-0.0221
-0.0269
-0.0296
-0.0373
-0.0516
-0.0693
-0.0866
-0.1017
-0.1149
-0.1284
-0.1395
-0.1274
-0.0749
-0.0117
0.0122
-0.0124
-0.0559
-0.0885
-0.1017
-0.1045
-0.11
-0.1243
-0.1458
-0.1704
-0.1892
-0.161
-0.0388
0.1277
0.2175
0.1879
0.0998
0.0286
0.003559
0.007744
0.0102
-0.00835
-0.042
-0.0487
0.0159
0.1243
0.1943
0.1881
0.1379
0.0934
0.0776
0.0831
0.0902
0.0859
0.0692
0.0472
0.027
0.0113
0.006137
0.0225
0.053
0.0708
0.0637
0.0425
0.0236
0.0143
0.0121
0.0104
0.004291
-0.007184
-0.0216
-0.0357
-0.0482
-0.0585
-0.0627
-0.0548
-0.0358
-0.0134
0.00617
0.0211
0.0277
0.0216
0.00553
-0.0121
-0.0258
-0.0359
-0.0461
-0.0596
-0.0768
-0.0956
-0.1139
-0.1308
-0.147
-0.1635
-0.1723
-0.153
-0.0965
-0.0245
0.0236
0.025
-0.006335
-0.0394
-0.0536
-0.0427
-0.0126
0.0196
0.0335
0.0226
-0.001254
-0.0221
-0.033
-0.037
-0.0408
-0.0487
-0.0609
-0.0752
-0.0895
-0.1026
-0.1146
-0.1263
-0.1384
-0.1511
-0.1646
-0.1786
-0.1909
-0.1979
-0.1983
-0.1957
-0.1941
-0.1934
-0.1888
-0.1768
-0.1484
-0.0745
0.0445
0.1374
0.1443
0.086
0.0204
-0.0137
-0.004913
0.04
0.0988
0.1485
0.1818
0.2072
0.2304
0.2419
0.2319
0.2067
0.1802
0.1597
0.1438
0.1273
0.1067
0.0811
0.0516
0.0202
-0.0109
-0.0403
-0.0688
-0.0975
-0.1271
-0.1579
-0.1781
-0.1557
-0.0755
0.022
0.0706
0.0481
-0.0101
-0.0585
-0.079
-0.0815
-0.0852
-0.1011
-0.1287
-0.1604
-0.1894
-0.2138
-0.2366
-0.2539
-0.2299
-0.1202
0.0498
0.1818
0.2062
0.148
0.0803
0.0504
0.0574
0.0735
0.076
0.0716
0.0908
0.1461
0.2058
0.2279
0.2069
0.1703
0.1447
0.1365
0.1366
0.1335
0.122
0.1041
0.0847
0.0677
0.0533
0.043
0.0448
0.0663
0.1014
0.132
0.142
0.132
0.1159
0.1062
0.1041
0.1011
0.0902
0.0716
0.0504
0.0307
0.0121
-0.000042
0.0241
0.0979
0.1685
0.1802
0.1401
0.09
0.0599
0.052
0.0517
0.0447
0.0264
0.001223
-0.0241
-0.0457
-0.0641
-0.0803
-0.0898
-0.0845
-0.0647
-0.0401
-0.0186
-0.002672
0.003997
-0.002786
-0.0197
-0.0377
-0.0514
-0.0615
-0.0715
-0.084
-0.0993
-0.1157
-0.1316
-0.1463
-0.1603
-0.1745
-0.1861
-0.1838
-0.1583
-0.1159
-0.0718
-0.0367
-0.0164
-0.0141
-0.0254
-0.0387
-0.0462
-0.0486
-0.0485
-0.044
-0.0318
-0.0138
0.004902
0.0194
0.0246
0.0202
0.0128
0.0103
0.0149
0.023
0.0289
0.0287
0.0224
0.0139
0.006441
-0.000287
-0.008797
-0.0201
-0.0331
-0.0464
-0.0563
-0.0557
-0.0409
-0.0179
0.004536
0.024
0.0393
0.043
0.0306
0.00912
-0.0114
-0.0254
-0.0288
-0.0181
0.002097
0.0195
0.0226
0.0109
-0.006024
-0.0189
-0.0252
-0.0283
-0.0324
-0.04
-0.0479
-0.0467
-0.0311
-0.009104
0.008346
0.0171
0.017
0.0104
0.001666
-0.005125
-0.009063
-0.009655
-0.005717
0.001803
0.009268
0.0126
0.0104
0.005043
0.000218
-0.00231
-0.003245
-0.00424
-0.005961
-0.005873
0.000031
0.012
0.0264
0.0395
0.0469
0.0444
0.0332
0.0219
0.0171
0.0184
0.0209
0.0207
0.0174
0.0126
0.005942
-0.004817
-0.016
-0.0177
-0.005781
0.0113
0.0206
0.0169
0.007134
0.001021
0.002373
0.008279
0.0139
0.0161
0.0142
0.00903
0.002882
-0.002306
-0.006678
-0.0117
-0.0181
-0.0255
-0.0338
-0.0447
-0.06
-0.0774
-0.0866
-0.0749
-0.0421
-0.00142
0.0308
0.0417
0.0305
0.0102
-0.005413
-0.012
-0.0139
-0.0138
-0.008906
0.002902
0.0184
0.0335
0.0463
0.0539
0.0517
0.0403
0.0257
0.0137
0.004893
-0.002469
-0.009946
-0.0181
-0.027
-0.0362
-0.0422
-0.0387
-0.024
-0.003848
0.015
0.0272
0.0284
0.0196
0.009157
0.005683
0.0107
0.0187
0.0226
0.019
0.0103
0.001096
-0.005432
-0.009396
-0.0126
-0.0166
-0.022
-0.027
-0.0276
-0.0218
-0.0122
-0.002367
0.003657
0.001327
-0.0106
-0.0272
-0.0429
-0.0539
-0.0561
-0.0462
-0.0283
-0.011
-0.000932
0.001049
-0.000135
0.000371
0.001529
-0.000118
-0.003501
-0.006005
-0.008432
-0.013
-0.0202
-0.0283
-0.0361
-0.043
-0.0497
-0.0554
-0.0572
-0.0538
-0.0474
-0.0411
-0.0378
-0.0407
-0.0509
-0.0657
-0.0783
-0.0776
-0.0591
-0.0367
-0.0263
-0.0262
-0.0228
-0.00824
0.0117
0.0261
0.0315
0.0322
0.0335
0.0371
0.0414
0.0445
0.0462
0.0466
0.0459
0.0452
0.0443
0.0383
0.022
-0.003392
-0.024
-0.0231
-0.000445
0.0294
0.0534
0.0651
0.062
0.0479
0.0339
0.0285
0.0312
0.0358
0.0374
0.0354
0.0315
0.0281
0.0262
0.0256
0.0253
0.0252
0.0261
0.0283
0.0311
0.0333
0.0331
0.0291
0.0223
0.015
0.008158
0.000869
-0.007488
-0.0168
-0.0261
-0.0334
-0.0365
-0.0363
-0.0354
-0.0357
-0.037
-0.0381
-0.0378
-0.0363
-0.0344
-0.0328
-0.0307
-0.0269
-0.0216
-0.0166
-0.0152
-0.0179
-0.0227
-0.0272
-0.0301
-0.0153
0.0374
0.107
0.1523
0.1612
0.1477
0.1302
0.1205
0.1199
0.1221
0.1207
0.113
0.1008
0.0876
0.0758
0.0661
0.0575
0.0485
0.0387
0.0283
0.0177
0.007809
0.001512
0.002023
0.008365
0.0161
0.021
0.0195
0.0101
-0.003821
-0.0168
-0.0238
-0.0211
-0.0114
-0.005194
-0.0122
-0.0319
-0.0554
-0.075
-0.0893
-0.1018
-0.1165
-0.1349
-0.1567
-0.1753
-0.1598
-0.0743
0.0525
0.1347
0.1284
0.068
0.0116
-0.0112
0.000514
0.0384
0.0844
0.1106
0.1034
0.076
0.0508
0.0419
0.0513
0.0675
0.0715
0.0534
0.0208
-0.0123
-0.0402
-0.066
-0.0945
-0.1276
-0.1606
-0.1736
-0.1445
-0.0838
-0.0329
-0.0218
-0.0418
-0.0611
-0.056
-0.0276
0.00752
0.0352
0.0524
0.0644
0.0776
0.0936
0.1091
0.1207
0.1246
0.1154
0.0935
0.0665
0.0419
0.0215
0.002584
-0.0177
-0.0413
-0.0701
-0.1042
-0.1416
-0.1652
-0.1372
-0.0501
0.0374
0.0659
0.038
-0.004124
-0.0267
-0.0247
-0.0127
-0.006331
-0.014
-0.0362
-0.0664
-0.087
-0.0748
-0.0248
0.0401
0.0832
0.0813
0.0451
0.006523
-0.0109
-0.006159
0.005675
0.0105
0.003911
-0.0116
-0.0312
-0.0495
-0.0621
-0.0688
-0.0732
-0.0787
-0.0829
-0.0799
-0.0686
-0.0532
-0.0363
-0.0187
-0.004139
0.000585
-0.001729
0.002033
0.0194
0.0437
0.0582
0.0509
0.0262
-0.001661
-0.0228
-0.0372
-0.0487
-0.0502
-0.0315
-0.001749
0.018
0.0186
0.006896
-0.003869
-0.003272
0.009108
0.0257
0.0389
0.0457
0.0457
0.0395
0.0306
0.0228
0.0173
0.0121
0.002921
-0.0122
-0.0308
-0.0495
-0.0666
-0.0826
-0.0987
-0.1126
-0.1142
-0.0965
-0.0648
-0.0297
0.002926
0.0305
0.0462
0.0427
0.0252
0.0118
0.021
0.0549
0.096
0.1188
0.1094
0.0773
0.0438
0.0214
0.009212
-0.000098
-0.0106
-0.0225
-0.035
-0.0472
-0.0546
-0.0506
-0.0351
-0.0155
0.000026
0.007935
0.0106
0.0133
0.0191
0.0261
0.0295
0.0274
0.0246
0.0276
0.0378
0.0508
0.0623
0.0712
0.0783
0.0813
0.0757
0.0618
0.0452
0.0306
0.0187
0.007168
-0.004424
-0.0113
-0.009452
-0.002189
0.003773
0.005519
0.004416
0.005598
0.0163
0.0327
0.0387
0.0269
0.00558
-0.014
-0.0273
-0.0367
-0.0461
-0.0574
-0.069
-0.074
-0.0651
-0.0432
-0.0173
0.00543
0.0238
0.0392
0.0479
0.0451
0.0338
0.0214
0.0123
0.005856
-0.000495
-0.008919
-0.0196
-0.0317
-0.0409
-0.0383
-0.0192
0.009329
0.035
0.0501
0.0566
0.0613
0.066
0.0679
0.0656
0.0611
0.0564
0.0521
0.0474
0.042
0.0367
0.0316
0.0269
0.0251
0.0305
0.0436
0.0585
0.069
0.0748
0.082
0.0956
0.1092
0.1116
0.1012
0.0858
0.0731
0.065
0.0601
0.056
0.051
0.0447
0.0374
0.0301
0.0227
0.014
0.003443
-0.007961
-0.0191
-0.0297
-0.0399
-0.0501
-0.0596
-0.0681
-0.0758
-0.0834
-0.0914
-0.0995
-0.107
-0.1124
-0.1153
-0.1171
-0.1193
-0.1224
-0.1259
-0.1294
-0.1344
-0.1431
-0.1555
-0.1694
-0.1826
-0.1953
-0.2047
-0.1985
-0.1634
-0.1059
-0.0487
-0.008883
0.0151
0.035
0.0545
0.0663
0.0667
0.061
0.0595
0.0698
0.0901
0.1125
0.1298
0.1365
0.1308
0.1179
0.1057
0.0982
0.0941
0.0896
0.0833
0.078
0.0782
0.0835
0.0893
0.092
0.0917
0.0899
0.0885
0.0885
0.0895
0.0906
0.0916
0.092
0.0914
0.0863
0.0707
0.0431
0.009752
-0.018
-0.0301
-0.0269
-0.0206
-0.024
-0.0388
-0.0536
-0.0533
-0.0364
-0.0193
-0.0179
-0.0323
-0.0514
-0.0656
-0.0693
-0.0617
-0.0478
-0.0351
-0.0266
-0.021
-0.0172
-0.0185
-0.0276
-0.0411
-0.0512
-0.0494
-0.0366
-0.0247
-0.0236
-0.0316
-0.0391
-0.0365
-0.0231
-0.007839
0.000784
0.002265
0.002623
0.007474
0.0143
0.0165
0.0122
0.005243
-0.000824
-0.004878
-0.007035
-0.005462
0.001752
0.013
0.0248
0.0347
0.0394
0.0353
0.0234
0.009282
-0.003055
-0.0142
-0.0262
-0.0373
-0.044
-0.046
-0.0458
-0.0457
-0.0464
-0.047
-0.047
-0.0458
-0.0408
-0.0293
-0.0169
-0.0122
-0.0166
-0.025
-0.0329
-0.0387
-0.0433
-0.0479
-0.0526
-0.0558
-0.0571
-0.0582
-0.063
-0.0724
-0.0815
-0.086
-0.0878
-0.0891
-0.0824
-0.0589
-0.0249
0.002546
0.0145
0.0166
0.0188
0.0258
0.0361
0.0457
0.0507
0.0489
0.0417
0.0328
0.0262
0.0234
0.021
0.0102
-0.0143
-0.0388
-0.0396
-0.0102
0.0319
0.0602
0.0576
0.0302
-0.001537
-0.023
-0.0339
-0.0401
-0.0385
-0.0191
0.0174
0.0592
0.0929
0.1031
0.0823
0.0431
0.009392
-0.00432
-0.002533
0.002529
0.006473
0.0119
0.0195
0.0275
0.0341
0.0376
0.0365
0.0325
0.0292
0.0284
0.029
0.0274
0.0213
0.0124
0.00532
0.002678
0.001742
-0.002907
-0.0127
-0.0235
-0.0314
-0.0372
-0.0481
-0.0712
-0.0995
-0.1103
-0.0878
-0.0419
0.006185
0.0438
0.0719
0.0928
0.1016
0.0957
0.0826
0.0714
0.0649
0.0602
0.0538
0.044
0.0317
0.0192
0.00862
0.000976
-0.004722
-0.0111
-0.0219
-0.0373
-0.0453
-0.0281
0.0131
0.0523
0.0668
0.0589
0.0476
0.0476
0.0577
0.0647
0.0601
0.0461
0.0282
0.005421
-0.024
-0.0533
-0.074
-0.0865
-0.097
-0.1113
-0.1296
-0.1439
-0.1441
-0.1299
-0.1087
-0.0838
-0.0549
-0.027
-0.006351
0.007828
0.0189
0.0236
0.0179
0.004509
-0.009868
-0.0217
-0.0315
-0.0417
-0.0511
-0.051
-0.035
-0.0103
0.008371
0.0167
0.0195
0.0184
0.0134
0.006704
0.000767
-0.002207
0.000423
0.008594
0.0183
0.0253
0.0287
0.031
0.0363
0.0463
0.0592
0.0718
0.0787
0.073
0.0548
0.0334
0.0189
0.0133
0.0112
0.006921
-0.000283
-0.007248
-0.0113
-0.0127
-0.0134
-0.0164
-0.0248
-0.037
-0.0422
-0.028
0.005105
0.0444
0.0738
0.0794
0.0605
0.032
0.01
-0.000731
-0.005576
-0.0124
-0.0247
-0.0372
-0.0401
-0.0296
-0.0121
0.004591
0.0173
0.0253
0.0249
0.0133
-0.005169
-0.0232
-0.0367
-0.0469
-0.0574
-0.0696
-0.0781
-0.0745
-0.0582
-0.037
-0.0172
0.000286
0.0141
0.0217
0.0239
0.0245
0.0247
0.0234
0.0201
0.0153
0.009236
0.00105
-0.0102
-0.0206
-0.0212
-0.009086
0.009208
0.0259
0.0383
0.0472
0.0508
0.046
0.0346
0.0221
0.012
0.004134
-0.003819
-0.0112
-0.0121
-0.002734
0.009845
0.0131
0.002614
-0.0149
-0.0306
-0.0411
-0.0482
-0.0544
-0.0584
-0.0569
-0.0506
-0.043
-0.0367
-0.0333
-0.0356
-0.0441
-0.0553
-0.0655
-0.073
-0.0735
-0.0603
-0.0339
-0.003302
0.0221
0.0385
0.0495
0.0602
0.0682
0.0658
0.0522
0.0347
0.0207
0.0126
0.007404
0.001035
-0.005899
-0.006775
0.002582
0.0186
0.0351
0.0484
0.0554
0.0529
0.0423
0.0297
0.0197
0.0126
0.006227
-0.00172
-0.0105
-0.0155
-0.0129
-0.004395
0.00509
0.0121
0.0158
0.0186
0.0247
0.0339
0.0395
0.0344
0.02
0.003572
-0.00975
-0.019
-0.0265
-0.035
-0.0444
-0.0528
-0.0592
-0.0637
-0.0635
-0.0554
-0.041
-0.0237
-0.001636
0.0288
0.0635
0.0871
0.0863
0.0657
0.0418
0.0263
0.0195
0.0148
0.006378
-0.007191
-0.0238
-0.038
-0.0408
-0.0289
-0.0131
-0.008626
-0.0182
-0.033
-0.0422
-0.0407
-0.0316
-0.022
-0.0169
-0.0169
-0.0195
-0.0215
-0.0209
-0.0163
-0.007654
0.003461
0.0141
0.0229
0.0306
0.0367
0.0374
0.0304
0.0192
0.008523
0.00023
-0.006749
-0.013
-0.014
-0.004192
0.0136
0.0288
0.0343
0.0328
0.0309
0.032
0.0336
0.0315
0.0259
0.0207
0.0185
0.0187
0.0193
0.0188
0.0162
0.009744
0.001413
-0.001915
0.003085
0.0128
0.021
0.0218
0.0126
-0.002931
-0.0191
-0.0337
-0.0478
-0.0628
-0.0784
-0.0911
-0.0965
-0.0932
-0.0849
-0.0771
-0.0703
-0.0606
-0.0458
-0.0282
-0.0124
-0.001259
0.006461
0.0133
0.0183
0.0186
0.014
0.007503
0.001788
-0.002583
-0.006563
-0.0112
-0.0168
-0.0212
-0.0201
-0.0116
0.000918
0.013
0.0227
0.0308
0.0392
0.0487
0.0591
0.0685
0.0722
0.0673
0.0568
0.046
0.0378
0.0318
0.0258
0.0184
0.009443
0.001669
0.0003
0.007495
0.0194
0.0309
0.0395
0.0448
0.0454
0.0417
0.0364
0.0323
0.03
0.0303
0.0343
0.0409
0.0478
0.0527
0.0538
0.0509
0.0461
0.042
0.0393
0.0367
0.0336
0.0307
0.028
0.0254
0.0225
0.0187
0.0142
0.00815
-0.004998
-0.0303
-0.0607
-0.082
-0.0892
-0.0887
-0.086
-0.0818
-0.077
-0.0741
-0.0711
-0.0601
-0.0381
-0.0147
-0.000556
0.003741
0.004675
0.006049
0.006647
0.00492
0.001367
-0.00285
-0.006689
-0.009229
-0.0104
-0.011
-0.0123
-0.0154
-0.0203
-0.0237
-0.0215
-0.0138
-0.004242
0.004042
0.009917
0.0111
0.005291
-0.004184
-0.0114
-0.0143
-0.015
-0.0143
-0.0106
-0.003229
0.005822
0.0145
0.0214
0.0233
0.0173
0.005679
-0.006545
-0.0161
-0.0207
-0.0173
-0.006095
0.008127
0.0206
0.0301
0.0358
0.0347
0.0259
0.0137
0.004011
-0.000225
0.000011
0.001482
0.001481
-0.000602
-0.0036
-0.006136
-0.007557
-0.008065
-0.008417
-0.009465
-0.0115
-0.0141
-0.0167
-0.0191
-0.0211
-0.023
-0.0248
-0.0267
-0.0288
-0.0309
-0.033
-0.0351
-0.0371
-0.0392
-0.0422
-0.0463
-0.049
-0.0466
-0.0384
-0.0282
-0.0192
-0.0122
-0.005999
0.001251
0.009186
0.0133
0.009956
0.000711
-0.009111
-0.0126
-0.004064
0.0149
0.0366
0.0548
0.0664
0.069
0.0641
0.0573
0.0519
0.0451
0.0329
0.0155
-0.00387
-0.0226
-0.036
-0.0371
-0.0245
-0.005772
0.009149
0.0144
0.0115
0.007038
0.006865
0.0134
0.0252
0.038
0.0488
0.0544
0.0499
0.0342
0.0132
-0.00458
-0.0136
-0.0139
-0.0108
-0.00945
-0.0117
-0.0162
-0.0215
-0.0264
-0.0303
-0.0337
-0.0371
-0.0401
-0.0386
-0.0283
-0.0109
0.008114
0.0243
0.0342
0.0357
0.0309
0.0255
0.0231
0.0243
0.0278
0.0318
0.0344
0.0337
0.0284
0.0196
0.009486
-0.000974
-0.0116
-0.0225
-0.0337
-0.0433
-0.0466
-0.0413
-0.0313
-0.0219
-0.0151
-0.009756
-0.003749
0.003706
0.0116
0.0187
0.0246
0.0299
0.0323
0.0282
0.0172
0.003234
-0.009863
-0.0211
-0.0316
-0.043
-0.0548
-0.0615
-0.0561
-0.0396
-0.0201
-0.004715
0.003565
0.00642
0.008513
0.0133
0.021
0.0292
0.035
0.0371
0.0369
0.0365
0.0372
0.0389
0.0409
0.0432
0.0465
0.0501
0.0514
0.0486
0.0436
0.0399
0.0389
0.0395
0.0389
0.0341
0.0251
0.0144
0.005315
-0.000936
-0.004609
-0.005125
-0.002143
0.002717
0.007314
0.0108
0.0138
0.0159
0.0141
0.005673
-0.007782
-0.0226
-0.0365
-0.0496
-0.0624
-0.0726
-0.0772
-0.0775
-0.0769
-0.077
-0.0769
-0.0748
-0.0697
-0.0623
-0.0537
-0.0443
-0.0331
-0.0198
-0.004259
0.0124
0.0265
0.0348
0.0386
0.0415
0.045
0.0472
0.0467
0.044
0.0408
0.0375
0.0322
0.0233
0.011
-0.003788
-0.0194
-0.0331
-0.0432
-0.051
-0.0576
-0.0602
-0.0547
-0.0429
-0.0303
-0.0208
-0.0142
-0.008016
-0.000374
0.009258
0.021
0.034
0.0435
0.042
0.029
0.0117
-0.003145
-0.0138
-0.0225
-0.0322
-0.0435
-0.0554
-0.065
-0.068
-0.0635
-0.0547
-0.0446
-0.0342
-0.0238
-0.013
-0.001342
0.0112
0.0243
0.0373
0.0492
0.0594
0.0685
0.0776
0.0873
0.0947
0.0943
0.0841
0.069
0.0544
0.0403
0.0234
0.002989
-0.0173
-0.0329
-0.0437
-0.053
-0.0621
-0.0662
-0.0606
-0.0479
-0.0355
-0.03
-0.0312
-0.0342
-0.035
-0.0315
-0.0227
-0.009437
0.005245
0.018
0.0249
0.0225
0.0127
0.001527
-0.006787
-0.0123
-0.0173
-0.0238
-0.032
-0.0407
-0.0486
-0.0555
-0.0606
-0.0605
-0.0528
-0.0401
-0.0274
-0.0177
-0.0107
-0.004161
0.003807
0.0128
0.0207
0.0262
0.03
0.0336
0.0378
0.0426
0.0476
0.0519
0.054
0.0517
0.0455
0.0383
0.0326
0.0285
0.025
0.0213
0.0184
0.0169
0.0162
0.0155
0.0141
0.0115
0.007993
0.004336
0.00116
-0.001605
-0.004413
-0.006889
-0.007098
-0.003625
0.00239
0.008626
0.0138
0.0181
0.0224
0.0274
0.0329
0.0374
0.0385
0.0357
0.0311
0.0272
0.0258
0.028
0.0333
0.0393
0.0442
0.0479
0.0511
0.0544
0.0584
0.0625
0.0643
0.0613
0.0545
0.0461
0.0338
0.0133
-0.0112
-0.0295
-0.0384
-0.0427
-0.0487
-0.058
-0.0666
-0.0683
-0.062
-0.0522
-0.0436
-0.0389
-0.0382
-0.0398
-0.0408
-0.0401
-0.0383
-0.0361
-0.0304
-0.0159
0.005278
0.0247
0.0369
0.0439
0.0484
0.048
0.0402
0.0277
0.015
0.004416
-0.004338
-0.0127
-0.0218
-0.032
-0.0428
-0.051
-0.0522
-0.0457
-0.0359
-0.0273
-0.0215
-0.0177
-0.0139
-0.008895
-0.00288
0.001773
0.002547
-0.000294
-0.004214
-0.00718
-0.008851
-0.01
-0.0116
-0.0135
-0.0142
-0.0124
-0.009265
-0.006338
-0.004334
-0.002959
-0.001519
0.00043
0.002432
0.003218
0.002131
-0.0001
-0.001768
-0.00087
0.003174
0.008995
0.0153
0.022
0.029
0.0363
0.0438
0.0514
0.0594
0.0666
0.0663
0.051
0.0249
0.001902
-0.009653
-0.0131
-0.0163
-0.0238
-0.035
-0.0466
-0.0551
-0.0569
-0.0513
-0.0418
-0.0339
-0.0307
-0.0317
-0.0336
-0.0341
-0.0325
-0.0302
-0.0282
-0.0258
-0.0214
-0.0144
-0.004849
0.006091
0.0172
0.0278
0.0359
0.0385
0.0354
0.0301
0.0268
0.0264
0.0274
0.0277
0.0265
0.0249
0.0246
0.0267
0.0302
0.0337
0.0365
0.0384
0.0383
0.0358
0.0321
0.028
0.0224
0.0138
0.00285
-0.008793
-0.0189
-0.0241
-0.0224
-0.0163
-0.0113
-0.0105
-0.0126
-0.0147
-0.0151
-0.0137
-0.0118
-0.0104
-0.009933
-0.009902
-0.009894
-0.009038
-0.005747
0.000677
0.008861
0.017
0.0243
0.0314
0.0389
0.0456
0.0477
0.0429
0.034
0.0251
0.0184
0.0133
0.008101
0.001806
-0.005611
-0.0128
-0.0174
-0.0181
-0.0164
-0.0148
-0.0145
-0.0144
-0.0119
-0.005614
0.003241
0.0123
0.0203
0.0268
0.0285
0.0205
0.003636
-0.0148
-0.0284
-0.0347
-0.0324
-0.0228
-0.0111
-0.00137
0.00564
0.012
0.0182
0.0204
0.0158
0.006338
-0.003535
-0.0113
-0.0169
-0.0222
-0.0284
-0.0359
-0.0441
-0.0512
-0.0524
-0.0433
-0.0262
-0.007768
0.007732
0.0204
0.0327
0.0443
0.0495
0.0442
0.0319
0.019
0.008334
-0.001418
-0.0124
-0.0238
-0.03
-0.0267
-0.0167
-0.005803
0.002268
0.008031
0.0142
0.0228
0.0339
0.0459
0.0573
0.0677
0.0768
0.0814
0.0752
0.0574
0.0345
0.0134
-0.003704
-0.0188
-0.035
-0.0534
-0.073
-0.091
-0.1059
-0.1141
-0.1084
-0.0872
-0.0587
-0.0323
-0.0111
0.007715
0.0279
0.0512
0.077
0.1033
0.1233
0.1262
0.1095
0.0837
0.0604
0.0438
0.0307
0.0162
-0.002144
-0.0235
-0.0421
-0.0496
-0.0435
-0.0307
-0.0193
-0.0124
-0.009174
-0.0101
-0.0161
-0.0251
-0.0337
-0.0407
-0.0468
-0.0533
-0.0585
-0.0592
-0.0552
-0.0506
-0.0498
-0.0534
-0.0583
-0.062
-0.0635
-0.0627
-0.06
-0.0562
-0.05
-0.0388
-0.023
-0.005751
0.008917
0.0186
0.0241
0.0278
0.0307
0.0328
0.0346
0.0362
0.0384
0.0423
0.0485
0.0573
0.0672
0.0749
0.0785
0.0794
0.0806
0.0836
0.0872
0.0859
0.0741
0.0529
0.0293
0.00854
-0.008295
-0.0214
-0.0319
-0.041
-0.0472
-0.0467
-0.0399
-0.0308
-0.0228
-0.0166
-0.0136
-0.017
-0.0282
-0.0429
-0.0568
-0.0666
-0.0682
-0.0589
-0.0427
-0.0272
-0.0166
-0.0101
-0.005631
-0.002508
-0.000555
0.001108
0.004629
0.0127
0.0248
0.0369
0.0448
0.0478
0.048
0.0486
0.052
0.0598
0.0719
0.0827
0.0813
0.0626
0.0337
0.005507
-0.0167
-0.033
-0.045
-0.0551
-0.0657
-0.0753
-0.078
-0.0716
-0.0602
-0.049
-0.0415
-0.0404
-0.0462
-0.0553
-0.0634
-0.0686
-0.072
-0.0755
-0.0804
-0.0861
-0.091
-0.0945
-0.0966
-0.0965
-0.0938
-0.0895
-0.0828
-0.0682
-0.042
-0.008257
0.0257
0.0562
0.0843
0.1125
0.1354
0.1383
0.1159
0.0804
0.0473
0.0224
0.001945
-0.0198
-0.0442
-0.0672
-0.0821
-0.0819
-0.0669
-0.0456
-0.0268
-0.0128
-0.000784
0.0129
0.0294
0.0439
0.0465
0.0338
0.0131
-0.006368
-0.0208
-0.0326
-0.048
-0.0708
-0.0956
-0.1123
-0.1182
-0.1169
-0.1116
-0.1043
-0.0969
-0.09
-0.0829
-0.0751
-0.0666
-0.0584
-0.0513
-0.0434
-0.0306
-0.0112
0.012
0.0347
0.0552
0.0743
0.0943
0.1159
0.1341
0.1407
0.1353
0.1254
0.118
0.1139
0.1089
0.0995
0.0862
0.0716
0.059
0.0515
0.049
0.0485
0.0466
0.0418
0.0343
0.0261
0.0191
0.0149
0.0136
0.0135
0.0126
0.0105
0.007293
0.003423
-0.000607
-0.004349
-0.007452
-0.008654
-0.006652
-0.002208
0.002533
0.005826
0.006994
0.006252
0.004282
0.00205
0.000312
-0.00145
-0.004873
-0.0108
-0.0182
-0.0257
-0.0325
-0.0385
-0.0443
-0.0503
-0.0565
-0.0611
-0.0613
-0.0569
-0.05
-0.042
-0.0332
-0.0241
-0.0149
-0.006474
-0.000441
0.002368
0.002955
0.003261
0.004393
0.006307
0.008922
0.0125
0.0171
0.0215
0.0223
0.0161
0.004091
-0.008566
-0.0173
-0.0215
-0.0238
-0.0258
-0.0274
-0.0279
-0.0276
-0.026
-0.0212
-0.0124
-0.00236
0.005407
0.0101
0.0136
0.0171
0.0196
0.019
0.0155
0.0112
0.007463
0.005161
0.005438
0.008759
0.0133
0.0167
0.0179
0.0183
0.019
0.0209
0.0232
0.0243
0.0227
0.0189
0.0145
0.0111
0.0105
0.0138
0.0195
0.0254
0.0297
0.0318
0.0323
0.0325
0.0343
0.0389
0.0444
0.0463
0.0424
0.035
0.0276
0.0211
0.0145
0.007942
0.002273
-0.002553
-0.005639
-0.004662
0.000643
0.006865
0.0103
0.0105
0.009828
0.0102
0.0121
0.0147
0.0174
0.0205
0.0246
0.0295
0.0345
0.0393
0.0436
0.0459
0.0448
0.041
0.0367
0.0332
0.0302
0.0269
0.0228
0.0181
0.0123
0.004527
-0.004687
-0.0142
-0.023
-0.0298
-0.0343
-0.0373
-0.039
-0.0386
-0.0358
-0.0312
-0.0246
-0.0159
-0.006381
0.001773
0.007325
0.011
0.0147
0.0194
0.024
0.0253
0.0217
0.0153
0.009034
0.004752
0.002169
-0.000172
-0.003374
-0.007769
-0.0131
-0.0187
-0.0241
-0.0291
-0.0335
-0.0367
-0.0381
-0.0363
-0.0307
-0.0226
-0.0141
-0.005872
0.002284
0.0108
0.02
0.0284
0.032
0.0283
0.0196
0.0105
0.003737
-0.001189
-0.006135
-0.0123
-0.0198
-0.0276
-0.0325
-0.0311
-0.0239
-0.0155
-0.01
-0.008303
-0.008058
-0.006884
-0.004045
-0.000769
0.001334
0.001942
0.00181
0.001882
0.00248
0.003385
0.004858
0.007271
0.0104
0.0138
0.0165
0.0171
0.0143
0.009359
0.004221
0.000065
-0.002272
-0.00101
0.004313
0.0116
0.0182
0.0233
0.0278
0.0327
0.0375
0.0391
0.0352
0.028
0.0205
0.0146
0.009756
0.004239
-0.003209
-0.0122
-0.0215
-0.0298
-0.0361
-0.0401
-0.0422
-0.0425
-0.0412
-0.0395
-0.0381
-0.0357
-0.0308
-0.0237
-0.0165
-0.0119
-0.009959
-0.008988
-0.007086
-0.003507
0.001018
0.005358
0.009065
0.012
0.0136
0.0136
0.0123
0.011
0.0102
0.009829
0.009438
0.008955
0.0094
0.0118
0.0158
0.02
0.0237
0.0269
0.0299
0.0333
0.0365
0.0379
0.0361
0.0321
0.0279
0.0234
0.0162
0.005171
-0.005952
-0.0132
-0.0164
-0.0185
-0.0217
-0.0251
-0.0262
-0.024
-0.0201
-0.0164
-0.0136
-0.0118
-0.0118
-0.014
-0.0175
-0.0207
-0.0233
-0.0248
-0.0235
-0.0188
-0.0121
-0.005605
-0.00022
0.004514
0.008743
0.0108
0.009293
0.005086
0.000585
-0.002859
-0.005462
-0.008049
-0.0104
-0.0111
-0.009148
-0.005709
-0.002301
0.000319
0.002461
0.004724
0.007299
0.009613
0.0108
0.011
0.0109
0.0109
0.0114
0.0123
0.0139
0.016
0.0184
0.0201
0.0194
0.0151
0.008299
0.001687
-0.00287
-0.005424
-0.007367
-0.01
-0.0138
-0.018
-0.0212
-0.0218
-0.02
-0.0172
-0.015
-0.0137
-0.0127
-0.0116
-0.011
-0.0116
-0.0132
-0.015
-0.0161
-0.0156
-0.0139
-0.0121
-0.0108
-0.01
-0.009254
-0.007739
-0.004494
0.000191
0.004358
0.006572
0.007252
0.007812
0.009891
0.0148
0.0223
0.0305
0.0378
0.0434
0.0476
0.0513
0.0534
0.052
0.0473
0.0415
0.0364
0.0316
0.0252
0.0168
0.007036
-0.002273
-0.009403
-0.0134
-0.0148
-0.0149
-0.0152
-0.0163
-0.0177
-0.0192
-0.0213
-0.0243
-0.028
-0.0319
-0.035
-0.0355
-0.0325
-0.0273
-0.0222
-0.0177
-0.0129
-0.007374
-0.000974
0.006052
0.0131
0.0187
0.0219
0.0236
0.0256
0.0283
0.0309
0.0314
0.029
0.0252
0.0215
0.0187
0.0165
0.0145
0.0134
0.0134
0.0142
0.0149
0.0142
0.0117
0.007899
0.003993
0.000742
-0.001956
-0.00427
-0.005674
-0.005728
-0.004933
-0.00416
-0.003786
-0.003852
-0.004794
-0.007165
-0.0106
-0.0142
-0.0174
-0.0201
-0.0229
-0.0259
-0.0291
-0.0315
-0.0316
-0.0291
-0.0258
-0.0229
-0.0209
-0.0192
-0.0173
-0.0161
-0.0169
-0.0195
-0.0227
-0.0254
-0.0276
-0.0295
-0.0302
-0.028
-0.0227
-0.0162
-0.009926
-0.004105
0.001893
0.008455
0.0155
0.0229
0.0302
0.0364
0.0384
0.0345
0.0271
0.02
0.0147
0.0106
0.006314
0.001674
-0.001098
0.000017
0.00418
0.008903
0.0124
0.014
0.0142
0.014
0.0143
0.015
0.015
0.014
0.0125
0.0108
0.00914
0.007121
0.00459
0.001761
-0.001075
-0.00379
-0.006442
-0.009302
-0.0128
-0.017
-0.0215
-0.0259
-0.03
-0.0332
-0.0344
-0.0338
-0.032
-0.0287
-0.0232
-0.0161
-0.008658
-0.00233
0.001403
0.002383
0.002242
0.002675
0.004149
0.006212
0.00847
0.0108
0.013
0.0151
0.017
0.0176
0.0158
0.0121
0.008058
0.004689
0.002341
0.00132
0.001514
0.00212
0.002374
0.001964
0.000767
-0.001099
-0.003156
-0.004932
-0.006237
-0.007006
-0.007375
-0.00765
-0.008109
-0.008985
-0.0103
-0.0118
-0.0132
-0.0146
-0.0158
-0.017
-0.0184
-0.0204
-0.0233
-0.027
-0.0306
-0.0325
-0.0314
-0.0278
-0.0234
-0.0185
-0.0126
-0.005961
-0.000012
0.00336
0.003902
0.003315
0.003301
0.004426
0.006396
0.008599
0.0104
0.0108
0.009224
0.006015
0.002542
-0.000275
-0.002492
-0.004699
-0.007386
-0.0106
-0.0142
-0.0176
-0.0202
-0.0205
-0.0174
-0.0123
-0.007182
-0.003363
-0.001059
0.000368
0.001932
0.004129
0.006769
0.009411
0.0113
0.0115
0.0101
0.007914
0.006059
0.004872
0.004121
0.00356
0.003013
0.002346
0.00159
0.001004
0.000857
0.001143
0.001598
0.001969
0.00221
0.002327
0.00225
0.001949
0.001541
0.001174
0.000883
0.000572
-0.000086
-0.001474
-0.003564
-0.005888
-0.008007
-0.009861
-0.0116
-0.0135
-0.0151
-0.0155
-0.0144
-0.0125
-0.0105
-0.008578
-0.00651
-0.004294
-0.001921
0.000559
0.003106
0.005666
0.008221
0.0103
0.0108
0.009047
0.005933
0.002871
0.000532
-0.001157
-0.00231
-0.003088
-0.003821
-0.004793
-0.005842
-0.006388
-0.006119
-0.005281
-0.004328
-0.003387
-0.002333
-0.001144
0.000132
0.001474
0.002796
0.003759
0.003895
0.003174
0.002104
0.001202
0.000614
0.000307
0.000394
0.000861
0.001447
0.001941
0.002274
0.002543
0.002877
0.003255
0.00329
0.00251
0.001028
-0.000618
-0.002077
-0.003341
-0.004527
-0.005748
-0.006877
-0.0075
-0.007344
-0.006633
-0.005669
-0.004185
-0.001885
0.000929
0.003789
0.006334
0.007936
0.008078
0.007092
0.005907
0.005461
0.006503
0.009162
0.0126
0.0158
0.0183
0.0191
0.0174
0.0138
0.009971
0.006896
0.004901
0.003923
0.003551
0.003118
0.002249
0.001011
-0.000301
-0.001402
-0.002379
-0.00374
-0.0059
-0.008698
-0.0115
-0.0134
-0.0131
-0.0106
-0.007265
-0.004343
-0.002204
-0.00045
0.001094
0.002157
0.002647
0.002817
0.003002
0.00332
0.00366
0.00392
0.004173
0.004933
0.006827
0.009606
0.012
0.0127
0.0117
0.0102
0.009108
0.008743
0.008688
0.008475
0.007929
0.007518
0.008024
0.009643
0.0118
0.0137
0.0153
0.0163
0.0163
0.0152
0.0137
0.0123
0.0114
0.0108
0.0104
0.0101
0.009762
0.009261
0.008721
0.008223
0.007569
0.006051
0.003242
-0.000311
-0.003767
-0.006825
-0.009509
-0.0119
-0.0141
-0.0163
-0.0185
-0.0207
-0.0217
-0.0207
-0.0183
-0.0157
-0.0137
-0.0122
-0.0107
-0.008951
-0.006699
-0.00377
-0.000148
0.003792
0.007426
0.009908
0.0109
0.011
0.0113
0.012
0.0132
0.0149
0.0172
0.0199
0.0227
0.0253
0.0274
0.0279
0.0261
0.0227
0.0193
0.0167
0.0146
0.0125
0.009846
0.006785
0.003537
0.000371
-0.002615
-0.005358
-0.007794
-0.009977
-0.0121
-0.0143
-0.0167
-0.019
-0.0214
-0.0237
-0.0259
-0.0276
-0.0283
-0.0278
-0.0267
-0.0249
-0.0225
-0.0199
-0.0172
-0.0146
-0.012
-0.00884
-0.004403
0.001174
0.007129
0.0127
0.0173
0.0204
0.0223
0.0239
0.025
0.0252
0.0248
0.0241
0.0235
0.0226
0.0211
0.0187
0.0159
0.0133
0.0109
0.008387
0.005458
0.002035
-0.001538
-0.004779
-0.007441
-0.009538
-0.0109
-0.0113
-0.0111
-0.0109
-0.011
-0.0112
-0.0114
-0.0117
-0.0126
-0.014
-0.0152
-0.015
-0.0135
-0.0115
-0.009974
-0.008678
-0.007037
-0.00471
-0.001938
0.000903
0.0036
0.006128
0.008627
0.0112
0.0139
0.0167
0.0191
0.0199
0.0182
0.0149
0.0114
0.008771
0.006682
0.004795
0.003378
0.002917
0.003199
0.003604
0.003772
0.003722
0.003691
0.003809
0.003836
0.003427
0.002522
0.00148
0.00082
0.000682
0.000813
0.000897
0.000786
0.000557
0.000351
0.000229
0.000207
0.000131
-0.000538
-0.00225
-0.004777
-0.007466
-0.009878
-0.0119
-0.0129
-0.0121
-0.00959
-0.00665
-0.004178
-0.002294
-0.000602
0.001373
0.003741
0.006251
0.008341
0.009326
0.009174
0.008606
0.008307
0.008398
0.008624
0.008764
0.008713
0.008312
0.007604
0.005913
0.003914
0.002045
0.000393
-0.001165
-0.002719
-0.004021
-0.004551
-0.004268
-0.003701
-0.003305
-0.003166
-0.003189
-0.003401
-0.00395
-0.004801
-0.005719
-0.006533
-0.007258
-0.007982
-0.008787
-0.009525
-0.009786
-0.009391
-0.008643
-0.007924
-0.007392
-0.006999
-0.006699
-0.0067
-0.007182
-0.007971
-0.008781
-0.009491
-0.0101
-0.0108
-0.0112
-0.011
-0.01
-0.008718
-0.007509
-0.006526
-0.005717
-0.00524
-0.005476
-0.006441
-0.00763
-0.008631
-0.009398
-0.0101
-0.0109
-0.0119
-0.0127
-0.0127
-0.0113
-0.009236
-0.007162
-0.005497
-0.004153
-0.002854
-0.001401
0.000232
0.001883
0.003289
0.004372
0.005273
0.006181
0.007198
0.008304
0.009506
0.0109
0.0127
0.0146
0.0165
0.0183
0.02
0.0215
0.0227
0.0237
0.0244
0.0247
0.0249
0.025
0.0251
0.0253
0.0252
0.0248
0.0239
0.023
0.0223
0.0217
0.0211
0.0204
0.0196
0.0181
0.0159
0.0134
0.011
0.008906
0.00729
0.006414
0.006174
0.006109
0.005855
0.005374
0.004829
0.004288
0.003653
0.002895
0.002078
0.001273
0.000512
-0.000226
-0.000969
-0.001813
-0.002969
-0.004538
-0.006151
-0.00719
-0.007469
-0.007333
-0.007221
-0.007282
-0.007375
-0.007362
-0.007225
-0.007014
-0.006798
-0.006622
-0.006483
-0.006361
-0.006223
-0.006106
-0.006192
-0.006636
-0.007357
-0.008118
-0.008791
-0.009425
-0.0101
-0.0109
-0.0118
-0.0126
-0.0134
-0.0142
-0.0149
-0.0155
-0.0156
-0.0152
-0.0145
-0.0138
-0.0133
-0.013
-0.0126
-0.0121
-0.0111
-0.009655
-0.007929
-0.00626
-0.004846
-0.003735
-0.002835
-0.002147
-0.001905
-0.002181
-0.002693
-0.003132
-0.003426
-0.003654
-0.003913
-0.00426
-0.004674
-0.005139
-0.005691
-0.006329
-0.006988
-0.007597
-0.008156
-0.008681
-0.009031
-0.00893
-0.008332
-0.007524
-0.006826
-0.00645
-0.006439
-0.006595
-0.006681
-0.006628
-0.006484
-0.006346
-0.006289
-0.006216
-0.005874
-0.005101
-0.004045
-0.002984
-0.002072
-0.001304
-0.000568
0.000114
0.000639
0.001014
0.001324
0.00164
0.001985
0.002318
0.002452
0.002147
0.001423
0.000574
-0.000165
-0.00075
-0.001252
-0.00177
-0.002354
-0.002968
-0.003569
-0.004047
-0.004075
-0.003369
-0.002123
-0.000826
0.00022
0.001039
0.001806
0.002648
0.003581
0.004546
0.005338
0.005575
0.005078
0.004133
0.003204
0.002531
0.00205
0.001638
0.001312
0.001122
0.001047
0.000995
0.000892
0.000764
0.000659
0.000568
0.000445
0.000222
-0.000103
-0.000445
-0.000741
-0.000994
-0.001212
-0.001425
-0.00165
-0.001855
-0.001995
-0.001987
-0.001821
-0.001603
-0.001416
-0.001282
-0.00117
-0.001045
-0.000897
-0.00074
-0.00052
-0.000055
0.000769
0.001813
0.002821
0.003623
0.004125
0.004317
0.004353
0.004412
0.004569
0.004807
0.005112
0.005471
0.00582
0.006054
0.006041
0.005749
0.005324
0.004945
0.004672
0.004461
0.004243
0.003982
0.003689
0.003393
0.003116
0.002873
0.002675
0.002525
0.002399
0.002269
0.002128
0.001986
0.001853
0.001732
0.001624
0.001488
0.001257
0.00091
0.000516
0.000157
-0.00012
-0.000265
-0.000262
-0.00018
-0.000108
-0.000086
-0.000096
-0.000103
-0.000083
-0.000047
-0.000034
-0.000079
-0.000175
-0.000278
-0.000358
-0.000412
-0.000452
-0.000492
-0.000526
-0.000522
-0.000464
-0.000373
-0.000282
-0.000204
-0.000135
-0.00007
-0.000008205
0.000053
0.000111
0.000161
0.000188
0.000176
0.00013
0.000078
0.000039
0.000017
0.000002885
-0.000012
-0.00003
-0.000047
-0.00006
-0.000066
-0.000064
-0.000054
-0.00004
-0.000025
-0.000012
0.00000207
0.000017
0.000033
0.000051
0.000068
0.000076
0.000068
0.000054
0.000049
0.000073
0.000138
0.000223
0.000299
0.000351
0.000382
0.000401
0.000415
0.000428
0.000438
0.000442
0.00044
0.000432
0.000419
0.000402
0.000383
0.000358
0.000329
0.000298
0.000271
0.000247
0.000227
0.000208
0.000191
0.000176
0.000163
0.000153
0.000144
0.000137
0.000131
0.000125
0.000121
0.000116
0.00011
0.000103
0.000095

[/code]
berktaftali
Posts: 68
Joined: Fri Jul 02, 2004 6:10 am
Location: Computers and Structures, Inc.

Post by berktaftali »

Try reducing the time step by a factor (say 10) when a step does not converge. You can go back to your original time step after that. It is similar to what you have for switching between regular and modified Newton algorithms in the input file you posted...

PS. If g stands for the PGA of your ground acceleration record, do you really expect your bridge to survive a 20g earthquake?
Berk Taftali
Georgia Institute of Technology
Ph.D. Candidate, Structural Engineering, Mechanics, and Materials
School of Civil and Environmental Engineering
Atlanta, GA 30332 USA
Email: gte994y@mail.gatech.edu
lvy
Posts: 26
Joined: Tue Dec 14, 2004 8:22 pm
Location: Shanghai, China

Post by lvy »

Thank you for your answer, Berk Taftali. Perhaps I didn't make you understand.

I want to apply IDA(increment Dynamic Analysis) that scale the sesmic acceleration, so it may reach about 5g. I had try reducing the time step to 1/10 present time step for convegence, but it is unsuccessful. The 'g' in my script just is scale factor, not gravity acceleration. When 'g' equal 4.9, PGA equals about 0.25g. It is reasonable.
ozancemcelik
Posts: 3
Joined: Fri May 13, 2005 3:26 pm
Location: Georgia Tech

Post by ozancemcelik »

Try using displacement-based beam column element (dispBeamColumn) with a more refined mesh. In doing this, you may want to reduce the number of integration points as well.

Ozan.
lvy
Posts: 26
Joined: Tue Dec 14, 2004 8:22 pm
Location: Shanghai, China

Post by lvy »

Displacement based element suppose distributed plasticity has linear curvature distribution. The purpose of my research find the actual curvature distribution, so it couldn't meet. I will refine the element according your suggestion, Thank you!
Post Reply