In the current study, we developed a new computational methodology to simulate wound healing in soft tissues. We assumed that the injured tissue recovers partially its mechanical strength and stiffness by gradually increasing the volume fraction of collagen fibers. Following the principles of the constrained mixture theory, we assumed that new collagen fibers are deposited at homeostatic tension while the already existing tissue undergoes a permanent deformation due to the effects of remodeling. The model was implemented in the finite-element software Abaqus® through a VUMAT subroutine and applied to a complex and realistic case: simulating wound healing following midline laparotomy closure. The incidence of incisional hernia is still quite significant clinically, and our goal was to investigate different conditions hampering the success of these procedures. We simulated wound healing over periods of 6 months on a patient-specific geometry. One of the outcomes of the finite-element simulations was the width of the wound tissue, which was found to be clinically correlated with the development of incisional hernia after midline laparotomy closure. We studied the impact of different suturing modalities and the effects of situations inducing increased intra-abdominal pressure or its intermittent variations such as coughing. Eventually, the results showed that the main risks of developing an incisional hernia mostly depend on the elastic strains reached in the wound tissue after degradation of the suturing wires. Despite the need for clinical validation, these results are promising for establishing a digital twin of wound healing in midline laparotomy incision.