  Metallic glasses are frequently used as structural materials. Therefore, it
is important to develop methods to predict their mechanical response as a
function of the microstructure prior to loading. We develop a coarse-grained
spring network model, which describes the mechanical response of metallic
glasses using an equivalent series network of springs, which can break and
re-form to mimic atomic rearrangements during deformation. To validate the
model, we perform simulations of quasistatic, uniaxial tension of Lennard-Jones
and embedded atom method (EAM) potentials for Cu$_{50}$Zr$_{50}$ metallic
glasses. We consider samples prepared using a wide range of cooling rates and
with different amounts of crystalline order. We show that both the
Lennard-Jones and EAM models possess qualitatively similar stress $\sigma$
versus strain $\gamma$ curves. By specifying five parameters in the spring
network model (ultimate strength, strain at ultimate strength, slopes of
$\sigma(\gamma)$ at $\gamma=0$ and at large strain, and strain at fracture
where $\sigma=0$), we can accurately describe the form of the stress-strain
curves during uniaxial tension for the computational studies of
Cu$_{50}$Zr$_{50}$, as well as recent experimental studies of several Zr-based
metallic glasses. For the computational studies of Cu$_{50}$Zr$_{50}$, we find
that the yield strain distribution is shifted to larger strains for slowly
cooled glasses compared to rapidly cooled glasses. In addition, the average
number of new springs and their rate of formation decreases with decreasing
cooling rate. These effects offset each other at large strains, causing the
stress-strain curve to become independent of the sample preparation protocol in
this regime. In future studies, we will extract the parameters that define the
spring network model directly from atomic rearrangements that occur during
uniaxial deformation.
