Stress fibers are actomyosin bundles that regulate cellular mechanosensation and force transduction. Interacting with the extracellular matrix through focal adhesion complexes, stress fibers are highly dynamic structures regulated by myosin motors and crosslinking proteins. Under external mechanical stimuli such as tensile forces, the stress fiber remodels its architecture to adapt to external cues, displaying properties of viscoelastic materials. How the structural remodeling of stress fibers is related to the generation of contractile force is not well understood. In this work, we simulate mechanochemical dynamics and force generation of stress fibers using the molecular simulation platform MEDYAN. We model stress fiber as two connecting bipolar bundles attached at the ends to focal adhesion complexes. The simulated stress fibers generate contractile force that is regulated by myosin motors and alpha-actinin crosslinkers. We find that stress fibers enhance contractility by reducing the distance between actin filaments to increase crosslinker binding, and this structural remodeling ability depends on the crosslinker turnover rate. Under tensile pulling force, the stress fiber shows an instantaneous increase of the contractile forces followed by a slow relaxation into a new steady state. While the new steady state contractility after pulling depends only on the overlap between actin bundles, the short-term contractility enhancement is sensitive to the tensile pulling distance. We further show that this mechanical response is also sensitive to the crosslinker turnover rate. Our results provide new insights into the stress fiber mechanics that have significant implications for understanding cellular adaptation to mechanical signaling.