The umbrella sampling (US) method was used to calculate PMF for sliding of the molecular chains from the initial state to the deform state, where the sliding chain is kept straight. (a reasonable schematic is provided in Figure ). All US simulations were 20 ns long, of which the initial 10 ns were discarded as equilibration. Umbrella potentials with a force constant of 5000 kJ mol -1 •nm were applied to the chain center of mass and distributed with 0.2 nm spacing along the sliding direction. The pulling process of three architectures were presented in Figure . In all US calculations, positional restraints were unbiased using the weighted histogram analysis method (WHAM). To identity the simulated stress state in different weaving structures and impact condition, the virial stress per atom 𝝈 has been characterized in this paper, given as: