A mechanical modeling framework to study endothelial permeability

The inner lining of blood vessels, the endothelium, is made up of endothelial cells. Vascular endothelial (VE)-cadherin protein forms a bond with VE-cadherin from neighboring cells to determine the size of gaps between the cells and thereby regulate the size of particles that can cross the endothelium. Chemical cues such as thrombin, along with mechanical properties of the cell and extracellular matrix are known to affect the permeability of endothelial cells. Abnormal permeability is found in patients suffering from diseases including cardiovascular diseases, cancer, and COVID-19. Even though some of the regulatory mechanisms affecting endothelial permeability are well studied, details of how several mechanical and chemical stimuli acting simultaneously affect endothelial permeability are not yet understood. In this article, we present a continuum-level mechanical modeling framework to study the highly dynamic nature of the VE-cadherin bonds. Taking inspiration from the catch-slip behavior that VE-cadherin complexes are known to exhibit, we model the VE-cadherin homophilic bond as cohesive contact with damage following a traction-separation law. We explicitly model the actin cytoskeleton and substrate to study their role in permeability. Our studies show that mechanochemical coupling is necessary to simulate the influence of the mechanical properties of the substrate on permeability. Simulations show that shear between cells is responsible for the variation in permeability between bicellular and tricellular junctions, explaining the phenotypic differences observed in experiments. An increase in the magnitude of traction force due to disturbed flow that endothelial cells experience results in increased permeability, and it is found that the effect is higher on stiffer extracellular matrix. Finally, we show that the cylindrical monolayer exhibits higher permeability than the planar monolayer under unconstrained cases. Thus, we present a contact mechanics-based mechanochemical model to investigate the variation in the permeability of endothelial monolayer due to multiple loads acting simultaneously.</p

Experimental and numerical investigations of stress fibre reorientation in biological cells

Cells are the fundamental units of living organisms controlling the behaviour of tissues, organs and thereby the organ system. The mechanical stimuli has been found to contribute towards changes in mechanical properties of cells, sometimes resulting in diseases. In order to quantify the response of cells to mechanical stimuli, a variety of experiments have been carried out. They have helped us to understand the effects of different types of stimuli, responses in different length and time scales, effect of different environmental factors, and the relation between different components of the cell. Focal adhesions present on the cell membrane contain mechanosensitive proteins called integrins which can form a connection between the cell and the extra cellular matrix, and thus sense the properties of the substrate and thereby the external stimuli. This creates a chain of biochemical reactions within the cytoplasm, leading to a cross bridge between the actin and myosin proteins in the presence of calcium ions, forming stress fibres. In this regard, cyclic loading experiments have been performed to understand the nature of connection between the growth of focal adhesions and stress fibres. In this thesis, a novel DIY design of an uniaxial cell stretcher has been designed and the manufacturing process using 3D printing technology has been explained. The device has been used to apply uniaxial cyclic load with different amplitudes, keeping the frequency constant to study the response of cells to changing strains. The experiments have been performed on two types of cells, fibroblasts and osteoblasts. The results are analysed quantitatively and the stress fibre orientation is studied for varying loading conditions for each cell type. Since the experiments performed are in vitro, the numerical models are developed in order to apply in vivo type loading and study the response of cells. In this regard, a numerical model is developed wherein the already existing stress fibre and focal adhesion growth models are extended and coupled through a feedback loop involving cytoplasmic calcium concentration. Stress fibre is assumed to depend on the calcium concentration and the active stress, while focal adhesion is modelled by assuming that integrins which exist in two states are in thermodynamic equilibrium. The active stress is taken as a product of strain-rate and strain dependent functions. The focal adhesion forming a bond with the substrate, and the cell provides a traction force to the cell. To consider the variation of calcium concentration depending on the focal adhesion growth, a feedback loop is introduced. The effect of substrate stiffness on the response of cells is analysed. The model thus developed is used to obtain solutions to numerical problems simulating biological phenomenon including stress fibre reorientation due to changing amplitude of cyclic loading. The mathematical model developed results in a coupled system of equations, for which the solution scheme needs special consideration. Noting the limitation on the time step that could be used with the staggered coupling scheme, a monolithic coupling scheme is developed where the system of equations are solved simultaneously. The variational formulation of the constitutive equations have been derived along with the algorithmic aspects of the solution schemes. The solution is obtained by following the finite element method with large strain formulation. Numerical problems have been solved to compare the coupling schemes for space and time refinements, and for changing parameter values. The convergence behaviour of the solution schemes have been analysed. Finally, the solution obtained from the numerical scheme have been compared with the experimental observations, and the ways this project could be extended has been discussed. Some of the important results that could be drawn from this thesis are: 1) A novel compact 3D-printed cell stretcher which has the size of a standard 96-well plate. 2) Stress fibres are formed predominantly in directions away from the direction of loading, when cells are subjected to uni-axial cyclic stretch. 3) Reorientation is higher with increasing amplitude of the cyclic loading. 4) Mathematical model for stress fibre growth has been coupled with focal adhesions' and the cytoplasmic calcium is controlled through a feedback loop. 5) Feedback loop plays a crucial role in simulating experiments such as optical tweezers, and ROCK inhibition. 6) Uni-axial cyclic stretch simulation results could be qualitatively compared to the average results obtained from experiments. 7) Coupling the system of equations can be performed following either a staggered or a monolithic approach. 8) Monolithic approach allowed using higher time step values, while staggered approach had limitations in the time step, and the parameter values that could be used

Progressive Damage and Failure Analysis of Bonded Composite Joints at High Energy Dynamic Impacts

Both wing and fuselage structures utilize bonded composite joints for structural efficiency in modern commercial and military aircraft. To ensure compliance with certification requirements mechanical fasteners are typically used as a failsafe mechanism for appropriate strength in the event of complete stiffener disbond. However, the use of fasteners decreases the structural efficiency of the structure by adding weight. This establishes the requirement to better exploit the efficiency of bonded structures and fully understand the failure behavior of adhesively bonded composite structures, particularly when subjected to elevated loading rates due to high energy dynamic impacts (HEDI). For this reason, the NASA Advanced Composite Consortium (ACC) HEDI team developed an experimentation and numerical modeling program for high rate loading of composite joints. In the present work, the response of adhesively bonded composite joints subjected to elevated loading rates is studied numerically and validated against experimental results. Due to dynamic considerations of experiments, the idea of wedge insert was extended to use with Split Hopkinson Pressure Bar (SHPB) testing techniques. Mode-I and Mode-II test configurations were simulated to evaluate the capability of two continuum damage material (CDM) models in LS-DYNA, namely MAT162 and MAT261. Three different levels of fidelity were considered to investigate the level of detail required to numerically predict the failure behavior and the results from high fidelity analysis are presented

Evaluation of LS-DYNA MAT162 for Modeling Composite Fastener Joints for High Rates of Loading

In the present work, the behavior of composite-fastener joints in bearing failure at dynamic stroke rates of 500 in/s, 300 in/s and 100 in/s has been evaluated through progressive damage analysis (PDA) material model in LS-DYNA, namely MAT162. Two joint types: titanium pin and Hi-Lok fastener were analyzed to identify the differences between without and with preload conditions. A meso-level approach where each lamina was modeled separately was employed and a contact definition based on fracture toughness data was defined to represent composite delamination behavior. Test fixture had been modeled in a detailed manner to account for the dynamic effects and the simulation results were validated against experimental data. Preliminary test-analysis correlation indicated that MAT162 predicted results conservatively when compared to tests. Debris accumulation were observed to greatly affect the test results which were not considered in the current modelling strategies

On the monolithic and staggered solution of cell contractility and focal adhesion growth

The mechanical response of cells to stimuli tightly couples biochemical and biomechanical processes, which describe fundamental properties such as cell growth and reorientation. Cells interact continuously with their external surroundings, the extracellular matrix (ECM), by establishing a bond between cell and ECM through the formation of focal adhesions. Focal adhesions are made up of integrins, which are mechanosensitive proteins and responsible for the communication between the cell and the ECM. The governing biochemomechanical processes can be modeled by means of a continuum approach considering mechanical and thermodynamic equilibrium to describe cell contractility and focal adhesion growth. The immanent multiphysical character of cell mechanics involves important aspects such as the coupling of fields of different scales and corresponding interface conditions that are sensitive to the solution of the governing numerical problem. These aspects become even more relevant when considering a feedback loop among the multiphysical solutions fields. In this contribution, we consider solution properties and sensitivity aspects of a nonlinear mechanical continuum model for the prognosis of stress fiber growth and reorientation incorporating a mechanosensitive feedback loop. We provide the governing equations of a Hill model-based stress fiber growth, which is coupled to a thermodynamical approach modeling the focal adhesions. Furthermore, a variational formulation including the algebraic equations is derived for staggered and monolithic solution approaches and the reaction-diffusion equation that models the feedback mechanism. We test both schemes with regard to reliability, accuracy, and numerical efficiency for different model parameters and loading scenarios. We present algorithmic aspects of the considered solution schemes and reveal their robustness with regard to model refinement in space and time and finally provide an assessment of their overall solution performance for multiphysics problems in the context of cell mechanics

Experimental and numerical analysis of dynamic compressive response of Nomex honeycombs

Lightweight phenolic resin-impregnated aramid paper honeycombs, commercially known as Nomex®honeycombs, are promising cores for sandwich structures in aerospace applications due to their high ratios of stiffness and strength to density. The out-of-plane compressive properties of the Nomex honeycombs have been widely investigated under quasi-static and low strain rates (up to 300 s-1). There is a need to understand the behaviour of this structure under higher strain rate compression. This will widen the applicability of these structures to more areas such as debris impact and other impacts which induce high strain rates. This paper reports the out-of-plane compressive responses of Nomex honeycombs subject to quasi-static loading and high strain rate dynamic loading up to 1500 s-1. The work involves experimental measurements and numerical modelling and validation. The compressive responses of the honeycombs were measured using a sensitive magnesium alloy Kolsky bar setup with front and back face impacts. The failure modes of the Nomex honeycombs were identified to be different under quasi-static and dynamic compressions. Under quasi-static compression, the honeycombs failed with local phenolic resin fracture after the elastic buckling of the honeycomb walls. For the dynamic compression, the honeycombs failed with the stubbing of cell walls at the ends of specimens. A finite element (FE) numerical model was devised and validated with the experimental data. The FE model considered the strain rate effect of phenolic resin material. The model predictions were in good agreement with the experimental measurements and facilitated interpreting the out-of-plane compressive response of the Nomex honeycombs. It was shown that there was a linear compressive strength enhancement up to 30% from quasi-static to strain rate of 1500 s-1. The strength enhancement was governed by two mechanisms: the strain rate effect of the phenolic resin and inertial stabilization of the honeycomb unit cell walls, where 61%-74% of the enhancement was contributed by the inertial stabilization of the unit cell walls. In addition, it was shown that the impact method and initial imperfections had negligible effect on the compressive response of the Nomex honeycombs

A feedback-loop extended stress fiber growth model with focal adhesion formation

Contractile cells play a prominent role in the adaptive nature of biological tissues. Contractility is mainly attributed to the growth of the tension dependent actomyosin bundles called stress fibers within the cytoskeleton. Stress fibers extend along the length of the cell and end at focal adhesions on the cell membrane. At the focal adhesion junctions on the cell membrane the integrin proteins are capable of sensing the environment, thereby making the cellular behavior dependent on the cell supporting substrate. It has been observed that the growth of stress fibers influences focal adhesions and vice-versa, resulting in a continuous cross-talk between different processes in the cell. Recent experiments have shown that cells subjected to uni-axial cyclic loading, depending on the substrate properties reorient themselves in a direction away from the loading direction, exhibiting strain avoidance. Mathematical models are important to understand the dependence of the cellular behavior on the substrate properties along with feedback mechanisms and are further used in designing in-vitro experiments. The coupling of the models for stress fibers and focal adhesions results in a non-linear bio-chemo-mechanical problem. In this contribution, we present the positive influence of the growth of focal adhesions along with a mechanosensitive feedback loop on the stress fiber growth and further reveal the characteristics of the re-orientation process due to cyclic loading. We use a non-linear Hill-type model to capture the growth of the active stress involved in the evolution law for the stress fibers and a thermodynamical approach to model the focal adhesions. A highly stable and reliable monolithic solution scheme is used to solve the governing system of coupled equations. Finally, we validate our simulation results with experimental results in regard to different loading conditions