An efficient algorithm for sampling the shear-modulus reduction curve in the context of wave propagation using the elastoplastic Iwan model

The elastoplastic Iwan model has been used since the end of the 1970s to simulate nonlinear soil behaviour in seismic wave propagation. In this work, we present an automatic algorithm to efficiently sample the shear-modulus reduction curve in function of shear deformation, which constitutes the exclusive ingredient of the elastoplastic model. This model requires the data from the shear- modulus reduction as a function of shear deformation, which are readily available in the literature and from specific laboratory tests. The method involves a discretization and interpolation of these data to be used. The quality of the solution depends on the number of interpolated points. However, a larger number of them produce an increase of the computational time. To overcome this, we present an automatic algorithm to efficiently sample the shear-modulus reduction curve. We numerically prove that the chosen discretization of the curve has a strong impact on the calculation load, in addition to the well-known dependence on the input motion amplitude level. Two tests of nonlinear wave propagation in 1-D and 3-D media show the clear gain in computation time when using the proposed automatic sampling algorithm.