Biomechanics and Modeling in Mechanobiology最新文献

{"title":"Computational model of coarctation of the aorta in rabbits suggests persistent ascending aortic remodeling post-correction.","authors":"Ashley A Hiebing, Matthew A Culver, John F LaDisa, Colleen M Witzenburg","doi":"10.1007/s10237-025-01933-y","DOIUrl":"https://doi.org/10.1007/s10237-025-01933-y","url":null,"abstract":"<p><p>Coarctation of the aorta (CoA) is a common congenital cardiovascular lesion that presents as a localized narrowing of the proximal descending aorta. While improvements in surgical and catheter-based techniques have increased short-term survival, there is a high long-term risk of hypertension and a reduced average lifespan despite correction. Computational models can be used to estimate aortic remodeling and peripheral vascular compensation, potentially serving as key tools in developing a mechanistic understanding of the interplay between pre-treatment dynamics, post-treatment recovery, and long-term hypertension risk. In this study, we developed a lumped-parameter model of the heart and circulation to simulate CoA. After fitting model parameters using imaging and catheterization data from healthy rabbits, we then used the model to estimate differences in ascending aortic compliance and peripheral resistance between the healthy group and rabbits with both untreated and corrected CoA using their imaging and catheterization data. CoA was defined by the current putative clinical treatment threshold (a pressure gradient > 20 mm Hg). Model inputs were fitted such that outputs matched reported stroke volume, ejection fraction, systolic and diastolic aortic pressure, peak aortic flow, mean and peak blood pressure gradients, and upper-to-lower body flow split, with all results falling within one standard deviation of the data for all groups. In the untreated CoA and corrected simulations, a decrease in ascending aortic compliance was necessary to match reported hemodynamics. This suggests exposure to a pressure gradient > 20 mm Hg results in vascular remodeling that persists after repair, a process strongly correlated with hypertension.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143668568","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical simulation of voluntary respiration in a model of the whole human lower airway.","authors":"Xinying Ou, Jiahuan Meng, Chen Ma, Huajing Wan, Yu Chen, Fengming Luo","doi":"10.1007/s10237-025-01932-z","DOIUrl":"https://doi.org/10.1007/s10237-025-01932-z","url":null,"abstract":"<p><p>The lung model construction is limited to the local scale, and the numerical simulation of autonomous breathing is mostly computed from top to bottom in recent research. In this study, models of the entire lower airway from G0-G23 were constructed, and computational simulations were performed for the alveolar model using coupled fluid-solid analysis with pressure changes on the wall and for the rigid bronchial model using computational fluid dynamics by transmitting the boundary conditions step from bottom to top. This paper provides the results under spontaneous respiration, including the ventilation volume of the tracheobronchial tree, the situation of the internal flow field, and the mechanical characteristics of the lung tissues. The mechanical characteristics and the lung functions computed by the models were consistent with clinical or experimental data. This model could provide quantitative analysis results of respiratory mechanics in the lower respiratory tract of the human, which offers a reference for mechanical studies, such as the morphological changes and differentiation of cell types induced by force stimulation and tumor induction. Furthermore, various pathological models can be developed based on this model.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143646768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Advances in computational modeling of cytokine and growth factor dynamics in bone healing: a scoping review.","authors":"Ahmad Hedayatzadeh Razavi, Nazanin Nafisi, Maria Velasquez-Hammerle, Mohammad Javad Shariyate, Mohammad Khak, Alireza Mirahmadi, Megan McNichol, Edward K Rodrogiuez, Ara Nazarian","doi":"10.1007/s10237-025-01938-7","DOIUrl":"https://doi.org/10.1007/s10237-025-01938-7","url":null,"abstract":"<p><p>Bone healing is a complex process regulated by intricate biological and mechanical factors and spatially varied regions over time. This scoping review synthesizes current computational models that incorporate cytokines and growth factors, examining their role in bone healing. Through a systematic analysis of 71 studies, this review identifies and categorizes the modeling methodologies used, including mathematical, finite element, agent-based, mechanobiological, pharmacobiological, and hybrid approaches. The findings highlight the predominant use of mathematical models while noting a recent shift toward more sophisticated techniques like finite element and agent-based models. Key cytokines and growth factors, such as TGF-β, RANK-RANKL-OPG, and PTH, are repeatedly used, underscoring their essential roles in regulating cellular processes. This review also analyzes parameter selection and validation strategies, identifying gaps in current practices and emphasizing the need for high-quality experimental validation to improve model reliability. Some bibliometric analyses provide insights into citation networks and keyword co-occurrence, illustrating influential studies in the field and central themes. The findings offer a foundation for future research to enhance model accuracy, aiming toward more predictive and clinically relevant models accounting for biology and mechanics in bone healing.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143630137","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The order of precedence in treatment of multiple intracranial aneurysms: insights from a fluid-structure interaction study.","authors":"Kornelia M Kliś, Jerzy Gąsowski, Antoni Cierniak, Borys M Kwinta, Krzysztof Stachura, Tadeusz J Popiela, Igor Szydłowski, Bartłomiej Łasocha, Karolina Piotrowicz, Tomasz Grodzicki, Roger M Krzyżewski","doi":"10.1007/s10237-025-01928-9","DOIUrl":"https://doi.org/10.1007/s10237-025-01928-9","url":null,"abstract":"<p><p>Treatment strategy for multiple intracranial aneurysms is challenging, as in many cases the choice of the order in which to treat aneurysms is not based on high-quality evidence. We aimed to digitally simulate clinical scenarios of two different orders in which multiple intracranial aneurysms were treated and analyze changes in hemodynamics after first stage of treatment. We prospectively included patients with two intracranial aneurysms, with order of treatment difficult to determine based on clinical data alone. For each patient we prepared three models of arteries harboring aneurysms: with both aneurysms present and with one of them removed. Computational modeling of blood flow using fluid-structure interaction methodology was performed for each model. Hemodynamic parameters of aneurysm domes were compared between models with both aneurysm present, and models in with aneurysms were removed in changing order. In 25 included patients, the calculated hemodynamic parameters such as Time-Averaged Wall Shear Stress (0.46 ± 0.40 vs. 0.54 ± 0.44 Pa; p < 0.01) and surface vortex fraction (12.73% ± 7.92% vs. 14.26% ± 7.46%; p = 0.02) decreased after first stage of treatment, while Time-Averaged Wall Shear Stress Gradient (1.44 ± 0.41 vs. 1.34 ± 0.46 Pa; p = 0.04) and percentage of wall shear stress < 0.5 Pa (50.13% ± 33.01% vs. 44.08% ± 34.16%; p < 0.01) increased. Changes of wall shear stress in remaining aneurysm dome were independently correlated with dome-to-neck ratio of both removed and remaining aneurysms. Hemodynamics of untreated aneurysm worsens after first stage of treatment. Dome-to-neck ratio of both treated and untreated aneurysm was the strongest and independent predictor of that change.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143630140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A combined 4D flow MR imaging and fluid-structure interaction analysis of ascending thoracic aortic aneurysms.","authors":"Yu Zhu, Chlöe Armour, Binghuan Li, Selene Pirola, Yousuf Salmasi, Thanos Athanasiou, Declan P O'Regan, Xiao Yun Xu","doi":"10.1007/s10237-025-01939-6","DOIUrl":"https://doi.org/10.1007/s10237-025-01939-6","url":null,"abstract":"<p><p>This study aimed to characterize the altered hemodynamics and wall mechanics in ascending thoracic aortic aneurysms (ATAA) by employing fully coupled two-way fluid-structure interaction (FSI) analyses. Our FSI models incorporated hyperelastic wall mechanical properties, prestress, and patient-specific inlet velocity profiles (IVP) extracted from 4D flow magnetic resonance imaging (MRI). By performing FSI analyses on 7 patient-specific ATAA models and 6 healthy aortas, the primary objective of the study was to compare hemodynamic and biomechanical features in ATAA versus healthy controls. A secondary objective was to examine the need for 4D flow MRI-derived IVP in FSI simulations by comparing results with those using two commonly adopted idealized IVPs: Flat-IVP and Para-IVP for selected cases. Our results show that, compared to the healthy aortas, the ATAA models exhibited highly disturbed blood flow in the ascending aorta. Consequently, maximum turbulent kinetic energy (TKE) at peak systole (155.0 ± 188.4 Pa) and maximum time-averaged wall shear stress (TAWSS) (8.6 ± 6.5 Pa) were significantly higher in the ATAA cohort, compared to 0.6 ± 0.5 Pa and 2.8 ± 0.7 Pa in the healthy aortas. Peak wall stress was also nearly doubled in the ATAA group (414 ± 108 kPa vs. 215 ± 31 kPa). Additionally, comparisons of simulation results across models with different IVPs underscore the importance of prescribing 3D-IVP at the inlet, especially for ATAA cases. Using idealized IVPs in two selected ATAA models (P1 and P7) substantially reduced the maximum TKE from 571 Pa to 0.01 Pa (Flat-IVP) and 0.02 Pa (Para-IVP) in P1 and from 73 Pa to 0.01 Pa (Flat-IVP) and 0.08 Pa (Para-IVP) in P7, while the maximum TAWSS in the ascending aorta decreased from 9.6 Pa to 0.7 Pa (Flat-IVP) and 0.9 Pa (Para-IVP) in P1, and from 3.6 Pa to 1.2 Pa and 0.9 Pa, respectively, in P7. Moreover, idealized IVPs also caused the peak wall stress to reduce by up to 11.5% in P1 with severe aortic valve stenosis, and by up to 2% in P7 with mild aortic regurgitation. These results highlight the importance of FSI simulations combined with 4D flow MRI in capturing realistic hemodynamic and biomechanical changes in aneurysmal aortas.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143603215","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Surface-based versus voxel-based finite element head models: comparative analyses of strain responses.","authors":"Zhou Zhou, Xiaogai Li, Svein Kleiven","doi":"10.1007/s10237-025-01940-z","DOIUrl":"https://doi.org/10.1007/s10237-025-01940-z","url":null,"abstract":"<p><p>Finite element (FE) models of the human head are important injury assessment tools but developing a high-quality, hexahedral-meshed FE head model without compromising geometric accuracy is a challenging task. Important brain features, such as the cortical folds and ventricles, were captured only in a handful of FE head models that were primarily developed from two meshing techniques, i.e., surface-based meshing with conforming elements to capture the interfacial boundaries and voxel-based meshing by converting the segmented voxels into elements with and without mesh smoothing. Despite these advancements, little knowledge existed of how similar the strain responses were between surface- and voxel-based FE head models. This study uniquely addressed this gap by presenting three anatomically detailed models - a surface-based model with conforming meshes to capture the cortical folds-subarachnoid cerebrospinal fluid and brain-ventricle interfaces, and two voxel-based models (with and without mesh smoothing) - derived from the same imaging dataset. All numerical settings in the three models were exactly the same, except for the meshes. These three models were employed to simulate head impacts. The results showed that, when calculating commonly used injury metrics, including the percentile strains below the maximum (e.g., 99 percentile strain) and the volume of brain element with the strain over certain thresholds, the responses of the three models were virtually identical. Different strain patterns existed between the surface- and the voxel-based models at the interfacial boundary (e.g., sulci and gyri in the cortex, regions adjacent to the falx and tentorium) with strain differences exceeding 0.1, but remarkable similarities were noted at the non-interfacial region. The mesh smoothing procedure marginally reduced the strain discrepancies between the voxel- and surface-based model. This study yielded new quantitative insights into the general similarity in the strain responses between the surface- and voxel-based FE head models and underscored that caution should be exercised when using the strain at the interface to predict injury.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143603230","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Static and dynamic optimisation of fluid-filled responsive orthotic insoles.","authors":"Dayna Cracknell, Mark Battley, Justin Fernandez, Maedeh Amirpour","doi":"10.1007/s10237-025-01935-w","DOIUrl":"https://doi.org/10.1007/s10237-025-01935-w","url":null,"abstract":"<p><p>This study was focused on developing an optimisation-based methodology to create customised solid-liquid composite (SLC) orthotic insoles. The goal was to reduce peak plantar pressures through gait through a dynamic numerical optimisation. A gait simulation was developed through a series of numerical models with increasing complexity. These models were validated against experimental analyses. The insole was designed based on numerical optimisation techniques that regionally tailored the insole with the aim to reduce temporal peak pressures. A prototype of the optimised insole was created using additive manufacturing and tested experimentally. The numerical gait simulation showed good correlation with experimental results. The largest differences are attributed to the bone geometry adopted from a previous study from a subject of different age, gender and size demographics. The optimisation process showed significant reductions in peak plantar pressures in the static peak pressures by approximately 8% and in the summation of dynamic peak pressures by 50%. Experimental validation confirmed the numerical predictions, highlighting the effectiveness of the optimised insole. The findings suggest that the optimised insoles can improve plantar pressure distributions and reduce peak pressures, making them a viable alternative to traditional orthotic insoles. Future research should focus on more accurate geometry for the numerical models and clinical trials.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143539864","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Phase-field simulation of crack growth in cortical bone microstructure: parameter identification and comparison against experiments.","authors":"Jenny Carlsson, Olivia Karlsson, Hanna Isaksson, Anna Gustafsson","doi":"10.1007/s10237-025-01929-8","DOIUrl":"https://doi.org/10.1007/s10237-025-01929-8","url":null,"abstract":"<p><p>Computational models are commonly used to investigate how the cortical bone microstructure affects fracture resistance; recently, phase-field models have been introduced for this purpose. However, experimentally measured material parameters for the microstructural tissues are lacking. Moreover, as no validation studies have been published, it remains unclear to what extent classical phase-field methods, assuming linear-elastic, brittle fracture, accurately represent bone. In this study, we address both these shortcomings by first applying a design-of-experiments methodology to calibrate a set of material parameters for a two-dimensional phase-field finite element model of bovine osteonal microstructure. This was achieved by comparing the outcomes from simulation to data from single-edge notched bending experiments on bovine osteonal bone and subsequent imaging of the crack path. Second, we used these parameters in new bone geometries to evaluate the parameters and the predictive performance of the model. Reasonable agreement was achieved between prediction and experiments in terms of peak load, crack initiation toughness and crack path. However, the model is unable to capture the experimentally observed gradual evolution of damage, leading to a nonlinear force response before the onset of visible crack extension. Nor does it capture the similarly observed increase in toughness with increasing crack length. These limitations are inherent to all classical phase-field methods since they originate from theories of brittle fracture, and alternative formulations are discussed. This is the first study attempting to validate classical phase-field methods in simulation of cortical bone fracture, and it highlights both potential and limitations to be addressed in future work.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143536306","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Homogenized multiscale modelling of an electrically active double poroelastic material representing the myocardium.","authors":"Laura Miller, Raimondo Penta","doi":"10.1007/s10237-025-01931-0","DOIUrl":"https://doi.org/10.1007/s10237-025-01931-0","url":null,"abstract":"<p><p>In this work, we present the derivation of a novel model for the myocardium that incorporates the underlying poroelastic nature of the material constituents as well as the electrical conductivity. The myocardium has a microstructure consisting of a poroelastic extracellular matrix with embedded poroelastic myocytes, i.e. a double poroelastic material. Due to the sharp length scale separation that exists between the microscale, where the individual myocytes are clearly resolved from the surrounding matrix, and the length of the entire heart muscle, we can apply the asymptotic homogenization technique. The novel PDE model accounts for the difference in the electric potentials, elastic properties as well as the differences in the hydraulic conductivities at different points in the microstructure. The differences in these properties are encoded in the coefficients and are to be computed by solving differential cell problems arising when applying the asymptotic homogenization technique. We present a numerical analysis of the obtained Biot's modulus, Young's moduli as well as shears and the effective electrical activity. By investigating the poroelastic and electrical nature of the myocardium in one model, we can understand how the differences in elastic displacements between the extracellular matrix and the myocytes affect mechanotransduction and the influence of disease.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143497565","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Constitutive neural networks for main pulmonary arteries: discovering the undiscovered.","authors":"Thibault Vervenne, Mathias Peirlinck, Nele Famaey, Ellen Kuhl","doi":"10.1007/s10237-025-01930-1","DOIUrl":"https://doi.org/10.1007/s10237-025-01930-1","url":null,"abstract":"<p><p>Accurate modeling of cardiovascular tissues is crucial for understanding and predicting their behavior in various physiological and pathological conditions. In this study, we specifically focus on the pulmonary artery in the context of the Ross procedure, using neural networks to discover the most suitable material model. The Ross procedure is a complex cardiac surgery where the patient's own pulmonary valve is used to replace the diseased aortic valve. Ensuring the successful long-term outcomes of this intervention requires a detailed understanding of the mechanical properties of pulmonary tissue. Constitutive artificial neural networks offer a novel approach to capture such complex stress-strain relationships. Here, we design and train different constitutive neural networks to characterize the hyperelastic, anisotropic behavior of the main pulmonary artery. Informed by experimental biaxial testing data under various axial-circumferential loading ratios, these networks autonomously discover the inherent material behavior, without the limitations of predefined mathematical models. We regularize the model discovery using cross-sample feature selection and explore its sensitivity to the collagen fiber distribution. Strikingly, we uniformly discover an isotropic exponential first-invariant term and an anisotropic quadratic fifth-invariant term. We show that constitutive models with both these terms can reliably predict arterial responses under diverse loading conditions. Our results provide crucial improvements in experimental data agreement, and enhance our understanding into the biomechanical properties of pulmonary tissue. The model outcomes can be used in a variety of computational frameworks of autograft adaptation, ultimately improving the surgical outcomes after the Ross procedure.</p>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":" ","pages":""},"PeriodicalIF":3.0,"publicationDate":"2025-02-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143481893","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}