Skip to content

1_2_0
Compare
Choose a tag to compare
@dunyuliu dunyuliu released this 24 Sep 19:40
1_2_0
66958e3

Major Achievements:

  • Combining the first cycle and later cycles calculated by EQquasi in one code;
  • Simplification of code structure. No variable transference.
  • Simulation of BP5-QD (quasi-dynamic with a radiated damping term).

Version: M-1.2.0
Author: DL
Branch: Master Parent version: M-1.1.0

Problems:

  • The first earthquake starts quickly and yields reasonable ‘rupture’ time contour. However, later cycles show abrupt change of maximum slip rate in the transition zone of frictional parameter a. A 1e-5 m/s will suddenly jump out of a sequence of 1e-9 m/s slip rates. Will model sizes affect it? The old models have 120 by 60 by 120 in x, y, z directions. Now it is changed to 160 by 70 by 50 instead.
    [Solved] Because the critical slip distance L is increased to 0.13/0.14m in BP5 compared to the 0.04 m in BP4, so the previous formula to calculate time step dtev = ksi * L/maxsliprate yields too rough time step and solution blows up. Ksi is reduced to 0.05 from 0.2 in the current model and it works. [Continued] Ksi is reduced to 0.025 based on equations in Lapusta et al. (2000). A matlab script is developed to process it.

Changes:

  • New: A license is created.
  • New: Use one code to simulate all earthquake cycles. For the first cycle, initial conditions are assigned. For later ones, initial conditions are imported from cplot_EQquasi.txt.
  • New: Significant simplification of the code structure. No variable is transferred between major subroutines.
  • New: files are renamed.
  • New: library_bound.f90 is created. It contains two subroutines bound_load and bound_ft_ku to handle far_field loading and calculate KU on faults for use in subroutine faulting.
  • New: library_output.f90 is created. It contains several subroutines output_onfault_st, output_offfault_st, output_onfault_transfer, output_timedy, output_globaldat, and output_prof. Future outputs could be implemented here.
  • New: Read_Input_Files.f90 is copied from EQdyna 3D version 5.2.0 to import station locations and to read in icstart from currentcycle.txt, which is created by the batch script. Other systematic parameters should be implemented in a later time.
  • New: crs.f90 is created to handle the CRS format of stiffness structure construction.
  • New: elemcal.f90 is created to handle the construction of the stiffness matrix and mass array in the CRS format.
  • New: ./script is created for MATLAB scripts for post processing.
    Improvement: on-screen information printing is improved.
  • [Verification] SCEC SEAS benchmark BP-5-QD.
Assets 2
Loading