Mechanism of modulus improvement for epoxy resin matrices: A molecular dynamics simulation