1D Nonlinear Numerical Methods Page 1 1D Nonlinear Numerical Methods Reading Assignment Lecture Notes Pp. 275-280 Kramer DEEPSOIL.pdf 2001 Darendeli, Ch. 10 Other Materials DeepSoil User's Manual 2001 Darendeli Homework Assignment #5 1. 2. 3. 4. Obtain the scaled Matahina Dam, New Zealand record from the course website and plot the following: (10 points) a. Plot the scaled acceleration time history b. Plot the scaled response spectrum Develop a soil profile for ground response analysis using soil properties for the I-15 project at 600 South Street (see attached) and the shear wave velocities found in SLC Vs profile.xls. (20 points) a. For sands, Darendeli, 2001 curves b. For silts, use Darendeli, 2001 with PI = 0 c. For clays, use Darendeli, 2001 curves with PI = 20 d. Treat layer 18 as a clay with PI = 20 and use Darendeli, 2001 curves e. Treat layer 19 as a sand and use Darendeli, 2001 curves f. For the bedrock velocity, use the velocity corresponding to the deepest Vs measurement in the soil profile with 2 percent damping Perform a site-specific, non-linear time domain ground response analysis for this soil profile using the pressure dependent hyperbolic model and Masing critera. Provide the following plots of the results: (15 points) a. Response spectrum summary b. Acceleration time histories for layer 1 c. pga profile Repeat problem 3 but perform a EQL analysis using the directions given in HW#3 problem 3. Plot a comparative plot of the response spectra using the spectrum from the nonlinear pressure dependent model (previous problem) versus the EQL pressure independent model (HW3 problem 4). (10 points). (SEE NEXT PG.)
1D Nonlinear Numerical Methods Page 2 Nonlinear Methods Homework Assignment #5 (cont.) 5. Modify the finite difference spreadsheet provided on the course website to include (20 points): a. Heterogeneous layers i. Varying thickness ii. Varying unit weight iii. Varying shear modulus b. Damping c. Given the information below, use the modified spreadsheet to perform a dynamic analysis for a duration of 2.0 s. Plot the response of the surface node versus time for verification: Layer # layer thickness unit weight Vs Damping (m) kn/m^3 (m/s) 1 1 19 150 5 2 1 19 170 5 3 1 19 190 5 4 0.5 20 150 5 5 1 20 150 5 6 0.5 20 150 5 7 2 20 150 5 8 1 21 170 5 9 1 21 170 5 10 1 21 170 5 Poisson ratio = 0.35 v(t) = A cos( t + ) A = 0.3 6.283 0.000 6. Verify your solution in 5 by performing an linear elastic analysis in DEEPSoil or FLAC for the same soil properties and velocity input (10 points).
1D Nonlinear Numerical Methods Page 3 Nonlinear Methods Homework Assignment #5 (cont.) 5. Solution (Excel) for uniform Vs = 80 m/s and 10 damping 6. Solution (FLAC)
1D Nonlinear Numerical Methods Page 4 Nonlinear Methods Homework Assignment #5 (cont.) 5. Solution (Excel) (first 5 time steps)
1D Nonlinear Numerical Methods Page 5 Comparison of 1D Equivalent Liner vs. 1D Nonlinear Methods EQL Method Nonlinear Methods
1D Nonlinear Numerical Methods Page 6 EQL vs NL Comparisons Target Spectrum for Comparisons
1D Nonlinear Numerical Methods Page 7 EQL vs NL Comparisons (cont.) EQL (Shake) Results at Surface from 5 km Convolution Nonlinear Results (DEEPSoil at Surface from 5 km Convolution
1D Nonlinear Numerical Methods Page 8 Lumped Mass System used in DeepSoil Fundamental Equation of Motion
1D Nonlinear Numerical Methods Page 9 DEEPSoil - Hyperbolic Model Modified Soil Hyperbolic Model used in DeepSoil
1D Nonlinear Numerical Methods Page 10 DEEPSoil (cont.) Introducing Pressure Dependency (Important for Deep Sediments)
1D Nonlinear Numerical Methods Page 11 DeepSoil (cont.) Incorporating Pressure Dependency in Damping [K] = stiffness matrix small strain viscous damping hysteretic damping incorporated by the hysteretic behavior of the soil
1D Nonlinear Numerical Methods Page 12 DEEPSoil (cont.) Pressure-dependent parameters b and d used to adjust curves in DEEPSoil. However, DARENDELI, 2001 has published newer curves based on confining pressure and PI. These are also incorporated in DEEPSoil.
1D Nonlinear Numerical Methods Page 13 Shear Modulus and Damping Curves from DARENDELI, 2001 As part of various research projects [including the SRS (Savannah River Site) Project AA891070. EPRI (Electric Power Research Institute) Project 3302. and ROSRINE (Resolution of Site Response Issues from the Northridge Earthquake) Project], numerous geotechnical sites were drilled and sampled. Intact soil samples over a depth range of several hundred meters were recovered from 20 of these sites. These soil samples were tested in the laboratory at The University of Texas at Austin (UTA) to characterize the materials dynamically. The presence of a database accumulated from testing these intact specimens motivated a re-evaluation of empirical curves employed in the state of practice. The weaknesses of empirical curves reported in the literature were identified and the necessity of developing an improved set of empirical curves was recognized. This study focused on developing the empirical framework that can be used to generate normalized modulus reduction and material damping curves. This framework is composed of simple equations. which incorporate the key parameters that control nonlinear soil behavior. The data collected over the past decade at The University of Texas at Austin are statistically analyzed using First-order. Second-moment Bayesian Method (FSBM). The effects of various parameters (such as confining pressure and soil plasticity on dynamic soil properties are evaluated and quantified within this framework. One of the most important aspects of this study is estimating not only the mean values of the empirical curves but also estimating the uncertainty associated with these values. This study provides the opportunity to handle uncertainty in the empirical estimates of dynamic soil properties within the probabilistic seismic hazard analysis framework. A refinement in site-specific probabilistic seismic hazard assessment is expected to materialize in the near future by incorporating the results of this study into the state of practice.
1D Nonlinear Numerical Methods Page 14 Effects of Mean Effective Stress on Shear Modulus and Damping Curves DARENDELI, 2001
1D Nonlinear Numerical Methods Page 15 Effects of Mean Effective Stress on Shear Modulus and Damping Curves (cont.) Shearing Strain (%) σo' = 0.25 atm σo' = 1.0 atm σo' = 4.0 atm σo' = 16 atm 1.00E-05 0.999 0.999 1.000 1.000 2.20E-05 0.998 0.999 0.999 1.000 4.84E-05 0.996 0.998 0.998 0.999 1.00E-04 0.993 0.995 0.997 0.998 2.20E-04 0.986 0.991 0.994 0.996 4.84E-04 0.971 0.981 0.988 0.992 1.00E-03 0.944 0.964 0.976 0.985 2.20E-03 0.891 0.928 0.952 0.969 4.84E-03 0.799 0.861 0.906 0.938 1.00E-02 0.671 0.761 0.832 0.885 2.20E-02 0.497 0.607 0.706 0.789 4.84E-02 0.324 0.428 0.538 0.645 1.00E-01 0.197 0.277 0.374 0.482 2.20E-01 0.107 0.157 0.225 0.311 4.84E-01 0.055 0.083 0.123 0.179 1.00E+00 0.029 0.044 0.067 0.101 Shearing Strain (%) σo' = 0.25 atm σo' = 1.0 atm σo' = 4.0 atm σo' = 16 atm 1.00E-05 1.201 0.804 0.539 0.361 2.20E-05 1.207 0.808 0.541 0.362 4.84E-05 1.226 0.820 0.548 0.367 1.00E-04 1.257 0.839 0.560 0.374 2.20E-04 1.330 0.884 0.588 0.391 4.84E-04 1.487 0.982 0.649 0.429 1.00E-03 1.792 1.174 0.769 0.503 2.20E-03 2.458 1.602 1.039 0.673 4.84E-03 3.762 2.474 1.607 1.035 1.00E-02 5.821 3.953 2.618 1.702 2.20E-02 9.097 6.579 4.572 3.075 4.84E-02 12.993 10.184 7.621 5.449 1.00E-01 16.376 13.788 11.134 8.573 2.20E-01 19.181 17.199 14.946 12.483 4.84E-01 20.829 19.565 17.990 16.070 1.00E+00 21.393 20.716 19.792 18.528 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 16 Effects of Mean Effective Stress on Shear Modulus and Damping Curves (cont.) Curve 1 Curve 2 Curve 1 - Sand Darendeli, 2001 v' (psf) = 11357 OCR = 1 Ko = 0.4 N = 10 F = 1 Hz Curve 2 - Sand Darendeli, 2001 v' (psf) = 576 OCR = 1 Ko = 0.4 N = 10 F = 1 Hz Curve 2 Curve 1 DEEPSoil V4.0
1D Nonlinear Numerical Methods Page 17 Effects of Plasticity on Shear Modulus and Damping Curves DARENDELI, 2001
1D Nonlinear Numerical Methods Page 18 Effects of Plasticity on Shear Modulus and Damping Curves (cont.) Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.999 1.000 1.000 1.000 1.000 2.20E-05 0.999 0.999 0.999 1.000 1.000 4.84E-05 0.998 0.998 0.999 0.999 0.999 1.00E-04 0.995 0.997 0.997 0.998 0.999 2.20E-04 0.991 0.993 0.995 0.996 0.997 4.84E-04 0.981 0.986 0.989 0.992 0.994 1.00E-03 0.964 0.973 0.979 0.984 0.989 2.20E-03 0.928 0.947 0.958 0.967 0.978 4.84E-03 0.861 0.896 0.917 0.934 0.956 1.00E-02 0.761 0.816 0.849 0.878 0.917 2.20E-02 0.607 0.682 0.732 0.778 0.843 4.84E-02 0.428 0.509 0.569 0.629 0.722 1.00E-01 0.277 0.348 0.404 0.465 0.571 2.20E-01 0.157 0.205 0.248 0.296 0.392 4.84E-01 0.083 0.111 0.137 0.169 0.238 1.00E+00 0.044 0.060 0.076 0.095 0.138 Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.804 0.997 1.191 1.450 2.096 2.20E-05 0.808 1.000 1.193 1.451 2.097 4.84E-05 0.820 1.008 1.199 1.456 2.100 1.00E-04 0.839 1.021 1.209 1.464 2.105 2.20E-04 0.884 1.053 1.234 1.482 2.117 4.84E-04 0.982 1.122 1.287 1.523 2.143 1.00E-03 1.174 1.257 1.392 1.603 2.193 2.20E-03 1.602 1.562 1.628 1.786 2.309 4.84E-03 2.474 2.198 2.128 2.175 2.560 1.00E-02 3.953 3.317 3.028 2.888 3.029 2.20E-02 6.579 5.440 4.803 4.343 4.029 4.84E-02 10.184 8.650 7.664 6.824 5.876 1.00E-01 13.788 12.217 11.092 10.024 8.541 2.20E-01 17.199 15.951 14.966 13.941 12.279 4.84E-01 19.565 18.829 18.185 17.458 16.132 1.00E+00 20.716 20.460 20.178 19.815 19.069 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 19 Shear Modulus and Damping Curves ( ' = 0.25 atm) DARENDELI, 2001
1D Nonlinear Numerical Methods Page 20 Shear Modulus and Damping Curves ( ' = 0.25 atm) Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.999 0.999 1.000 1.000 1.000 2.20E-05 0.998 0.999 0.999 0.999 1.000 4.84E-05 0.996 0.997 0.998 0.998 0.999 1.00E-04 0.993 0.995 0.996 0.997 0.998 2.20E-04 0.986 0.990 0.992 0.994 0.996 4.84E-04 0.971 0.979 0.983 0.987 0.991 1.00E-03 0.944 0.959 0.968 0.975 0.983 2.20E-03 0.891 0.919 0.936 0.949 0.966 4.84E-03 0.799 0.847 0.876 0.900 0.932 1.00E-02 0.671 0.739 0.783 0.822 0.876 2.20E-02 0.497 0.579 0.637 0.692 0.774 4.84E-02 0.324 0.400 0.459 0.521 0.625 1.00E-01 0.197 0.255 0.303 0.358 0.461 2.20E-01 0.107 0.142 0.174 0.213 0.293 4.84E-01 0.055 0.074 0.093 0.116 0.167 1.00E+00 0.029 0.040 0.050 0.063 0.093 Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 1.201 1.489 1.778 2.164 3.129 2.20E-05 1.207 1.493 1.781 2.166 3.131 4.84E-05 1.226 1.506 1.791 2.174 3.136 1.00E-04 1.257 1.528 1.808 2.187 3.144 2.20E-04 1.330 1.579 1.848 2.217 3.163 4.84E-04 1.487 1.690 1.933 2.282 3.204 1.00E-03 1.792 1.906 2.101 2.411 3.286 2.20E-03 2.458 2.387 2.476 2.702 3.472 4.84E-03 3.762 3.358 3.249 3.310 3.868 1.00E-02 5.821 4.977 4.581 4.386 4.593 2.20E-02 9.097 7.778 7.010 6.441 6.070 4.84E-02 12.993 11.489 10.477 9.589 8.579 1.00E-01 16.376 15.064 14.088 13.137 11.798 2.20E-01 19.181 18.334 17.640 16.904 15.716 4.84E-01 20.829 20.515 20.208 19.849 19.213 1.00E+00 21.393 21.507 21.542 21.547 21.544 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 21 Shear Modulus and Damping Curves ( ' = 1 atm) DARENDELI, 2001
1D Nonlinear Numerical Methods Page 22 Shear Modulus and Damping Curves ( ' = 1 atm) Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.999 1.000 1.000 1.000 1.000 2.20E-05 0.999 0.999 0.999 1.000 1.000 4.84E-05 0.998 0.998 0.999 0.999 0.999 1.00E-04 0.995 0.997 0.997 0.998 0.999 2.20E-04 0.991 0.993 0.995 0.996 0.997 4.84E-04 0.981 0.986 0.989 0.992 0.994 1.00E-03 0.964 0.973 0.979 0.984 0.989 2.20E-03 0.928 0.947 0.958 0.967 0.978 4.84E-03 0.861 0.896 0.917 0.934 0.956 1.00E-02 0.761 0.816 0.849 0.878 0.917 2.20E-02 0.607 0.682 0.732 0.778 0.843 4.84E-02 0.428 0.509 0.569 0.629 0.722 1.00E-01 0.277 0.348 0.404 0.465 0.571 2.20E-01 0.157 0.205 0.248 0.296 0.392 4.84E-01 0.083 0.111 0.137 0.169 0.238 1.00E+00 0.044 0.060 0.076 0.095 0.138 Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.804 0.997 1.191 1.450 2.096 2.20E-05 0.808 1.000 1.193 1.451 2.097 4.84E-05 0.820 1.008 1.199 1.456 2.100 1.00E-04 0.839 1.021 1.209 1.464 2.105 2.20E-04 0.884 1.053 1.234 1.482 2.117 4.84E-04 0.982 1.122 1.287 1.523 2.143 1.00E-03 1.174 1.257 1.392 1.603 2.193 2.20E-03 1.602 1.562 1.628 1.786 2.309 4.84E-03 2.474 2.198 2.128 2.175 2.560 1.00E-02 3.953 3.317 3.028 2.888 3.029 2.20E-02 6.579 5.440 4.803 4.343 4.029 4.84E-02 10.184 8.650 7.664 6.824 5.876 1.00E-01 13.788 12.217 11.092 10.024 8.541 2.20E-01 17.199 15.951 14.966 13.941 12.279 4.84E-01 19.565 18.829 18.185 17.458 16.132 1.00E+00 20.716 20.460 20.178 19.815 19.069 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 23 Shear Modulus and Damping Curves ( ' = 4 atm) DARENDELI, 2001
1D Nonlinear Numerical Methods Page 24 Shear Modulus and Damping Curves ( ' = 4 atm) Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 1.000 1.000 1.000 1.000 1.000 2.20E-05 0.999 1.000 1.000 1.000 1.000 4.84E-05 0.998 0.999 0.999 0.999 1.000 1.00E-04 0.997 0.998 0.998 0.999 0.999 2.20E-04 0.994 0.996 0.997 0.997 0.998 4.84E-04 0.988 0.991 0.993 0.995 0.996 1.00E-03 0.976 0.983 0.986 0.989 0.993 2.20E-03 0.952 0.965 0.972 0.978 0.986 4.84E-03 0.906 0.931 0.945 0.956 0.971 1.00E-02 0.832 0.873 0.898 0.918 0.945 2.20E-02 0.706 0.770 0.810 0.845 0.893 4.84E-02 0.538 0.618 0.673 0.725 0.802 1.00E-01 0.374 0.454 0.514 0.575 0.675 2.20E-01 0.225 0.287 0.339 0.396 0.501 4.84E-01 0.123 0.163 0.199 0.241 0.327 1.00E+00 0.067 0.091 0.113 0.140 0.200 Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.539 0.668 0.798 0.971 1.404 2.20E-05 0.541 0.670 0.799 0.972 1.405 4.84E-05 0.548 0.675 0.803 0.975 1.407 1.00E-04 0.560 0.683 0.809 0.980 1.410 2.20E-04 0.588 0.703 0.824 0.991 1.417 4.84E-04 0.649 0.745 0.857 1.016 1.433 1.00E-03 0.769 0.829 0.922 1.066 1.464 2.20E-03 1.039 1.021 1.070 1.180 1.537 4.84E-03 1.607 1.428 1.388 1.426 1.693 1.00E-02 2.618 2.173 1.977 1.886 1.991 2.20E-02 4.572 3.684 3.206 2.871 2.648 4.84E-02 7.621 6.235 5.387 4.693 3.934 1.00E-01 11.134 9.482 8.357 7.333 5.972 2.20E-01 14.946 13.400 12.231 11.056 9.226 4.84E-01 17.990 16.866 15.935 14.917 13.118 1.00E+00 19.792 19.158 18.571 17.876 16.513 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 25 Shear Modulus and Damping Curves ( ' = 16 atm) DARENDELI, 2001
1D Nonlinear Numerical Methods Page 26 Shear Modulus and Damping Curves ( ' = 16 atm) Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 1.000 1.000 1.000 1.000 1.000 2.20E-05 1.000 1.000 1.000 1.000 1.000 4.84E-05 0.999 0.999 0.999 1.000 1.000 1.00E-04 0.998 0.999 0.999 0.999 0.999 2.20E-04 0.996 0.997 0.998 0.998 0.999 4.84E-04 0.992 0.994 0.996 0.997 0.998 1.00E-03 0.985 0.989 0.991 0.993 0.996 2.20E-03 0.969 0.977 0.982 0.986 0.991 4.84E-03 0.938 0.954 0.964 0.972 0.981 1.00E-02 0.885 0.915 0.932 0.946 0.964 2.20E-02 0.789 0.839 0.869 0.895 0.929 4.84E-02 0.645 0.716 0.763 0.804 0.863 1.00E-01 0.482 0.564 0.623 0.679 0.764 2.20E-01 0.311 0.386 0.444 0.506 0.610 4.84E-01 0.179 0.233 0.279 0.331 0.431 1.00E+00 0.101 0.135 0.166 0.203 0.280 Shearing Strain (%) PI = 0 % PI = 15 % PI = 30 % PI = 50 % PI = 100 % 1.00E-05 0.361 0.448 0.534 0.650 0.941 2.20E-05 0.362 0.449 0.535 0.651 0.941 4.84E-05 0.367 0.452 0.538 0.653 0.942 1.00E-04 0.374 0.457 0.541 0.656 0.944 2.20E-04 0.391 0.469 0.551 0.663 0.949 4.84E-04 0.429 0.495 0.571 0.678 0.958 1.00E-03 0.503 0.547 0.611 0.709 0.978 2.20E-03 0.673 0.667 0.704 0.780 1.023 4.84E-03 1.035 0.924 0.903 0.934 1.120 1.00E-02 1.702 1.407 1.281 1.227 1.308 2.20E-02 3.075 2.433 2.100 1.871 1.729 4.84E-02 5.449 4.318 3.659 3.138 2.589 1.00E-01 8.573 7.021 6.022 5.151 4.049 2.20E-01 12.483 10.780 9.557 8.381 6.651 4.84E-01 16.070 14.619 13.472 12.268 10.241 1.00E+00 18.528 17.522 16.655 15.677 13.847 DARENDELI, 2001
1D Nonlinear Numerical Methods Page 27 Finite Difference Approach Wednesday, August 17, 2011 12:45 PM Note that with this approach we can approximate the change of things that vary either in space or time, or both. In regards to time, we will use the forward differencing approach in formulating the finite difference approach.
1D Nonlinear Numerical Methods Page 28 Finite Difference Approach (cont.) Wednesday, August 17, 2011 12:45 PM Finite difference calculation loop written with differential calculus
1D Nonlinear Numerical Methods Page 29 Finite Difference Approach (cont.) Finite difference calculation loop written with incremental approach
1D Nonlinear Numerical Methods Page 30 1D Finite Difference Solution for Wave Propagation Wednesday, August 17, 2011 12:45 PM
1D Nonlinear Numerical Methods Page 31 1D Finite Difference Solution for Wave Propagation (cont.) Tuesday, March 04, 2014 11:45 AM Steven F. Bartlett, 2014
1D Nonlinear Numerical Methods Page 32 1D Finite Difference Solution for Wave Propagation (cont.) Wednesday, August 17, 2011 12:45 PM
1D Nonlinear Numerical Methods Page 33 1D Finite Difference Solution for Wave Propagation (cont.) Wednesday, March 05, 2014 11:45 AM Steven F. Bartlett, 2014
1D Nonlinear Numerical Methods Page 34 1D Finite Difference Solution for Wave Propagation (cont.) Wednesday, August 17, 2011 12:45 PM 1 Displacement of Top Node vs Time 0.8 0.6 0.4 0.2 0 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5-0.2-0.4-0.6-0.8-1
1D Nonlinear Numerical Methods Page 35 1D Finite Difference Solution for Wave Propagation (cont.) Wednesday, August 17, 2011 12:45 PM ;FLAC verification of solution without damping config dynamic extra 5 grid 1 10 model elastic ini y mul 1 ;set dy_damp rayl 0.05 5; 5 percent damping at 5 hz fix y prop dens 2000 bulk 9.6E6 shear 3.2E6 def wave wave=amp*cos(omega*dytime) if dytime>=100 wave=0 endif end set amp=0.3 set omega = 6.283 apply xvel 1 hist wave yvel=0 j=1 his 1 xdisp i 1 j 1 his 2 xdisp i 1 j 11 his 3 xvel i 1 j 1 his 4 dytime set dytime = 0 ;set dydt = 0.0002; Flac can calc automatically solve dytime 5.01 save model2.sav 'last project state'
1D Nonlinear Numerical Methods Page 36 Incorporating Damping Note that shear resistance has two components: elastic and damping.
1D Nonlinear Numerical Methods Page 37 Incorporating Damping (cont.) Stiffness due to viscous damping
1D Nonlinear Numerical Methods Page 38 Hysteretic Damping as Implemented in FLAC Background The equivalent-linear method (see Section 3.2) has been in use for many years to calculate the wave propagation (and response spectra) in soil and rock, at sites subjected to seismic excitation. The method does not capture directly any nonlinear effects because it assumes linearity during the solution process; strain-dependent modulus and damping functions are only taken into account in an average sense, in order to approximate some effects of nonlinearity (damping and material softening). Although fully nonlinear codes such as FLAC are capable in principle of modeling the correct physics, it has been difficult to convince designers and licensing authorities to accept fully nonlinear simulations. One reason is that the constitutive models available to FLAC are either too simple (e.g., an elastic/plastic model, which does not reproduce the continuous yielding seen in soils), or too complicated (e.g., the Wang model [Wang et al. 2001], which needs many parameters and a lengthy calibration process). Further, there is a need to accept directly the same degradation curves used by equivalent-linear methods (see Figure 3.23 for an example), to allow engineers to move easily from using these methods to using fully nonlinear methods.
1D Nonlinear Numerical Methods Page 39 Hysteretic Damping (cont.) Formulation Modulus degradation curves, as illustrated in Figure 3.23, imply a nonlinear stress/strain curve. If we assume an ideal soil, in which the stress depends only on the strain (not on the number of cycles, or time), we can derive an incremental constitutive relation from the degradation curve, described by τe/γ = Ms, where τe is the normalized shear stress, γ the shear strain and Ms the normalized secant modulus. τe = Msγ (elastic component) Mt = dτe / dγ = Ms + γ dms / dγ (elastic and viscous component) where Mt is the normalized tangent modulus. The incremental shear modulus in a nonlinear simulation is then given by G Mt, where G is the small-strain shear modulus of the material.
1D Nonlinear Numerical Methods Page 40 Hysteretic Damping (cont.) FLAC code for single zone model with hysteretic damping conf dyn ext 5 grid 1 1 model elas prop dens 1000 shear 5e8 bulk 10e8 fix x y set dydt 1e-4 ini dy_damp hyst default -3.5 1.3 his sxy i 1 j 1 his xdis i 1 j 2 his nstep 1 ini xvel 1e-2 j=2 cyc 1000 ini xvel mul -1 cyc 250 ini xvel mul -1 cyc 500
1D Nonlinear Numerical Methods Page 41 Hysteretic Damping - Types of Tangent-Modulus Functions
1D Nonlinear Numerical Methods Page 42 Hysteretic Damping - Tangent-Modulus Functions (cont.) Default Model (cont)
1D Nonlinear Numerical Methods Page 43 Hysteretic Damping - Tangent-Modulus Functions (cont.)
1D Nonlinear Numerical Methods Page 44 Hysteretic Damping - Tangent-Modulus Functions (cont.)
1D Nonlinear Numerical Methods Page 45 Hysteretic Damping - Tangent-Modulus Functions (cont.)
1D Nonlinear Numerical Methods Page 46 Hysteretic Damping - Tangent-Modulus Functions (cont.)
1D Nonlinear Numerical Methods Page 47 Hysteretic Damping - Tangent-Modulus Functions (cont.) The parameters for the various tangent-modulus functions can be changed to fit or types of modulus reduction and damping data. To judge the fit of the function parameters to the experimental data, the following FLAC subroutine can be used. conf dy def setup givenshear = 1e8 ; shear modulus CycStrain = 0.01 ; cyclic strain (%) / 10 ;---- derived.. setvel = 0.1 * min(1.0,cycstrain/0.1) givenbulk = 2.0 * givenshear timestep = min(1e-4,1e-5 / CycStrain) nstep1 = int(0.5 + 1.0 / (timestep * 10.0)) nstep2 = nstep1 * 2 nstep3 = nstep1 + nstep2 nstep5 = nstep1 + 2 * nstep2 end setup ; gri 1 1 ;m mohr m elastic prop den 1000 sh givenshear bu givenbulk cohesion = 50e3 fix x y ini xvel setvel j=2 set dydt 1e-4 ini dy_damp hyst default -3.325 0.823; hysteretic damping his sxy i 1 j 1 his xdis i 1 j 2 his nstep 1 cyc nstep1 ini xv mul -1 cyc nstep2 ini xv mul -1 cyc nstep2 his write 1 vs 2 tab 1
1D Nonlinear Numerical Methods Page 48 Hysteretic Damping - Tangent-Modulus Functions (cont.) def HLoop emax = 0.0 emin = 0.0 tmax = 0.0 tmin = 0.0 loop n (1,nstep5) emax = max(xtable(1,n),emax) emin = min(xtable(1,n),emin) tmax = max(ytable(1,n),tmax) tmin = min(ytable(1,n),tmin) endloop slope = ((tmax - tmin) / (emax - emin)) / givenshear oo = out(' strain = '+string(emax*100.0)+'% G/Gmax = '+string(slope)) Tbase = ytable(1,nstep3) Lsum = 0.0 loop n (nstep1,nstep3-1) meant = (ytable(1,n) + ytable(1,n+1)) / 2.0 Lsum = Lsum + (xtable(1,n)-xtable(1,n+1)) * (meant - Tbase) endloop Usum = 0.0 loop n (nstep3,nstep5-1) meant = (ytable(1,n) + ytable(1,n+1)) / 2.0 Usum = Usum + (xtable(1,n+1)-xtable(1,n)) * (meant - Tbase) endloop Wdiff = Usum - Lsum Senergy = 0.5 * xtable(1,nstep1) * ytable(1,nstep1) Drat = Wdiff / (Senergy * 4.0 * pi) oo = out(' damping ratio = '+string(drat*100.0)+'%') end HLoop save singleelement.sav 'last project state'
1D Nonlinear Numerical Methods Page 49 Rayleigh Damping
1D Nonlinear Numerical Methods Page 50 Rayleigh Damping (cont.)
1D Nonlinear Numerical Methods Page 51 Hysteretic vs Rayleigh Damping Comparison
1D Nonlinear Numerical Methods Page 52
1D Nonlinear Numerical Methods Page 53 Blank
1D Nonlinear Numerical Methods Page 54
1D Nonlinear Numerical Methods Page 55 Blank
1D Nonlinear Numerical Methods Page 56
1D Nonlinear Numerical Methods Page 57 Blank