total shear modulus $G = G_\mathrm{A} -
G_\mathrm{NA}$ shows a low variation (by comparing $G_\mathrm{A} \simeq
G_\mathrm{NA} \simeq 900$ with $G \simeq 24$ at $\epsilon_\text{bend}=10^4$). The bending rigidity increases the affine shear modulus but, at the same time, the non-affine component also increases to cancel the increase in $G_\mathrm{A}$, and as a result, the total shear modulus presents a low increase. The elasticity of the shear deformation is therefore different from that of the bulk deformation, which is obvious when the elastic moduli are decomposed into affine and non-affine components.

Figure \[fig3\] also shows $K$ in (b) and $G$ in (d) for $L=3$. The values of $K$ and $G$ of $L=3$ are smaller than those of $L=50$, due to the bonding energy, $\varepsilon_\mathrm{FENE}$, connecting the monomers along the polymeric chains. The responses of $K$ and $G$ to the variation of $\varepsilon_\mathrm{bend}$ are also weaker for $L = 3$. However, $K$ and $G$, as well as affine $K_\mathrm{A}$ and $G_\mathrm{A}$ and non-affine $K_\mathrm{NA}$ and $G_\mathrm{NA}$, exhibit overall common dependences on $\varepsilon_\mathrm{bend}$ between $L=3$ and $50$. Therefore, the decrease in $\rho$ and increase in $\varepsilon_\mathrm{bend}$ engenders similar effects on the elasticity for $L=3$ and $50$.

Finally, it is remarked that the polymer 