Sharif Digital Repository

This paper presents a bi-directional closed-form analytical solution, in the framework of nonlinear soft composites mechanics, for top-down hyperelastic characterization of brain white matter tissue components, based on the directional homogenized responses of the tissue in the axial and transverse directions. The white matter is considered as a transversely isotropic neo-Hookean composite made of unidirectional distribution of axonal fibers within the extracellular matrix. First, two homogenization formulations are derived for the homogenized axial and transverse shear moduli of the tissue, based on definition of the strain energy density function. Next, the rule of mixtures and... 

The mechanical stability of conductive, single-walled carbon nanotubes (SWCNTs) under applied electric field and compressive loading is investigated. The distribution of electric charges on the nanotube surface is determined by employing a method based on the classical electrostatic theory. For mechanical stability analysis, a hybrid atomistic-structural element is proposed, which takes into account the nonlinear features of the stability. Nonlinear stability analysis based on an iterative solution procedure is used to determine the buckling force. The coupling between electrical and mechanical models is accomplished by adding Coulomb interactions to the mechanical model. The results show... 

Pregnant vehicle occupants experience relatively large acceleration when the vehicle passes a speed-bump. In this paper, the effect of such sudden acceleration on a pregnant uterus is investigated. A biomechanical model representing the fundamental dynamic behaviors of a pregnant uterus has been developed. The model relates to the 32nd week of gestation when the fetus is in head-down, occipito-anterior position. Considering the drag and squeeze effects of the amniotic fluid, we derive a comprehensive differential equation that represents the interaction of the uterus and fetus. Solving the governing equation, we obtain the system response to different speed-bump excitations. Using the fetal... 

Objective: Spine kinematics, kinetics, and trunk muscle activities were evaluated during different stages of a fatigue-induced symmetric lifting task over time. Background: Due to neuromuscular adaptations, postural behaviors of workers during lifting tasks are affected by fatigue. Comprehensive aspects of these adaptations remain to be investigated. Method: Eighteen volunteers repeatedly lifted a box until perceived exhaustion. Body center of mass (CoM), trunk and box kinematics, and feet center of pressure (CoP) were estimated by a motion capture system and force-plate. Electromyographic (EMG) signals of trunk/abdominal muscles were assessed using linear and nonlinear approaches. The L5-S1... 

Carbon dioxide storage in geological formations is one of the mature strategies developed for controlling global warming. This paper represents a comprehensive experimental and geochemical modeling study to analyze CO2-brine-rock interactions in a carbonate rock containing calcite and dolomite minerals. PHREEQC geochemical package has been applied for modeling the geochemical reactions in the studied porous media. Firstly, dynamic experiments are performed to calibrate the geochemical model. Then, static experiments are conducted to study the geochemical reactions in the CO2-brine-rock interaction system. This study contributes to analyzing the precipitation-dissolution and ion exchange... 

A diffusion tensor imaging (DTI) -based statistical micromechanical model was developed to study the effect of axonal fiber architecture on the inter- and intra-regional mechanical heterogeneity of the white matter. Three characteristic regions within the white matter, i.e., corpus callosum, brain stem, and corona radiata, were studied considering the previous observations of locations of diffuse axonal injury. The embedded element technique was used to create a fiber-reinforced model, where the fiber was characterized by a Holzapfel hyperelastic material model with variable dispersion of axonal orientations. A relationship between the fractional anisotropy and the dispersion parameter of... 

Muscle fatigue affects the underlying EMG-force relationship on which EMG-assisted biomechanical models rely. The aim of this study was to evaluate the impact of short duration muscle fatigue on the muscle gain value. Participants performed controlled, isometric trunk extension exertions at 10, 20, and 30 degrees of trunk flexion and controlled isokinetic trunk extension exertions at 5 and 15°/sec on five separate days. Fatigue of the lumbar extensors was generated by moderate-intensity, trunk extension exertions. Participants performed controlled test contractions at defined intervals throughout the fatiguing bout and the EMG activities of trunk muscles were collected. These EMG data were... 

In this work, thermo-mechanical behavior and microstructural evolution in similar and dissimilar friction stir welding of AA6061-T6 and AA5086-O have been investigated. Firstly, the thermo-mechanical behaviors of materials during similar and dissimilar FSW operations have been predicted using three-dimensional finite element software, ABAQUS, then, the mechanical properties and the developed microstructures within the welded samples have been studied with the aid of experimental observations and model predictions. It is found that different strengthening mechanisms in AA5086 and AA6061 result in complex behaviors in hardness of the welded cross section where the hardness variation in similar... 

In this paper, the continuum damage mechanics model originally proposed by Lemaitre (Journal of Engineering Materials and Technology. 1985; 107: 83-89) is presented through an adaptive finite element method for three-dimensional ductile materials. The macro-crack initiation-propagation criterion is used based on the distribution of damage variable in the continuum damage model. The microcrack closure effect is incorporated to simulate the damage evolution more realistic. The Zienkiewicz-Zhu posteriori error estimator is employed in conjunction with a weighted Superconvergence Patch Recovery (SPR) technique at each patch to improve the accuracy of error estimation and data transfer process.... 

A new, robust homogenization scheme for determination of the effective properties of a periodic piezoelectric composite with general multi-coated inhomogeneities is developed. In this scheme the coating does not have to be thin, the shape and orientation of the inclusion and coatings do not have to be identical, their centers do not have to coincide, their properties do not have to remain uniform, and the microstructure can be with the 2D elliptic or the 3D ellipsoidal inclusions. The development starts from the local electromechanical equivalent inclusion principle through the introduction of the position-dependent equivalent eigenstrain and electric field. Then with a Fourier series... 

A hallmark of more advanced models is their higher details of trunk muscles represented by a larger number of muscles. The question is if in reality we control these muscles individually as independent agents or we control groups of them called "synergy". To address this, we employed a 3-D biomechanical model of the spine with 18 trunk muscles that satisfied equilibrium conditions at L4/5, with different cost functions. The solutions of several 2-D and 3-D tasks were arranged in a data matrix and the synergies were computed by using non-negative matrix factorization (NMF) algorithms. Variance accounted for (VAF) was used to evaluate the number of synergies that emerged by the analysis, which... 

This paper presents a three-dimensional micromechanical model of brain white matter tissue as a transversely isotropic soft composite described by the generalized Ogden hyperelastic model. The embedded element technique, with corrected stiffness redundancy in large deformations, was used for the embedment of a histology-informed probabilistic distribution of the axonal fibers in the extracellular matrix. The model was linked to a multi-objective, multi-parametric optimization algorithm, using the response surface methodology, for characterization of material properties of the axonal fibers and extracellular matrix in an inverse finite element analysis. The optimum hyperelastic... 

In this research, the feasibility of friction-stir welding (FSW) for dissimilar lap-joining of an aluminium-magnesium alloy (AA5058) and poly-methyl-methacrylate sheets to attain sound and defect-free joints was examined. The inter-mixing flow patterns between the metal and polymer counterparts during FSW were predicted by employing three-dimensional finite element models. It is shown that the bonding mechanism between the dissimilar materials is mechanical interlocking at the interface which controls the joint strength depending on the processing parameters. The most suitable dissimilar lap-joining regarding microstructural soundness is attained at w= 1600 rev min−1 and v = 25 mm min−1.... 

The Calcium ion Ca2+ plays a critical role as an initiator and preserving agent of the cross-bridge cycle in the force generation of skeletal muscle. A new multi-scale chemo-mechanical model is presented in order to analyze the role of Ca2+ in muscle fatigue and to predict fatigue behavior. To this end, a cross-bridge kinematic model was incorporated in a continuum based mechanical model, considering a thermodynamic compatible framework. The contractile velocity and the generated active force were directly related to the force-bearing states that were considered for the cross-bridge cycle. In order to determine the values of the model parameters, the output results of an isometric simulation... 

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase... 

The 1991 NIOSH Lifting Equation (NLE) is widely used to assess the risk of injury to spine by providing estimates of the recommended weight limit (RWL) in hands. The present study uses the predictive equations developed based on a detailed trunk musculoskeletal biomechanical model to verify whether the RWL generates L5-S1 loads within the limits (e.g., 3400N for compression recommended by NIOSH and 1000N for shear recommended in some studies). Fifty lifting activities are simulated here to evaluate the RWL by the NLE and the L5-S1 loads by the predictive equations. In lifting activities involving moderate to large forward trunk flexion, the estimated RWL generates L5-S1 spine loads exceeding... 

Epidemiological studies are divided over the causative role of body weight (BW) in low back pain. Biomechanical modeling is a valuable approach to examine the effect of changes in BW on spinal loads and risk of back pain. Changes in BW have not been properly simulated by previous models as associated alterations in model inputs on the musculature and moment arm of gravity loads have been neglected. A detailed, multi-joint, scalable model of the thoracolumbar spine is used to study the effect of BW (varying at five levels, i.e., 51, 68, 85, 102, and 119kg) on the L5-S1 spinal loads during various static symmetric activities while scaling moment arms and physiological cross-sectional areas of... 

Different lifting analysis tools are commonly used to assess spinal loads and risk of injury. Distinct musculoskeletal models with various degrees of accuracy are employed in these tools affecting thus their relative accuracy in practical applications. The present study aims to compare predictions of six tools (HCBCF, LSBM, 3DSSPP, AnyBody, simple polynomial, and regression models) for the L4-L5 and L5-S1 compression and shear loads in twenty-six static activities with and without hand load. Significantly different spinal loads but relatively similar patterns for the compression (R2>0.87) were computed. Regression models and AnyBody predicted intradiscal pressures in closer agreement with... 

The human motor system is organized for execution of various motor tasks in a different and flexible manner. The kinetic redundancy in the human musculoskeletal system is a significant property by which the central nervous system achieves many complementary goals. An equilibrium-based biomechanical model of isometric three-dimensional exertions of trunk muscles has been developed. Following the definition and role of the uncontrolled manifold, the kinetic redundancy concept is explored in mathematical terms. The null space of the kinetically redundant system when a certain joint moment and/or stiffness are needed is derived and discussed. The aforementioned concepts have been illustrated,... 

Range of motion (ROM) of the thoracic spine has implications in patient discrimination for diagnostic purposes and in biomechanical models for predictions of spinal loads. Few previous studies have reported quite different thoracic ROMs. Total (T1-T12), lower (T5-T12) and upper (T1-T5) thoracic, lumbar (T12-S1), pelvis, and entire trunk (T1) ROMs were measured using an inertial tracking device as asymptomatic subjects flexed forward from their neutral upright position to full forward flexion. Correlations between body height and the ROMs were conducted. An effect of measurement errors of the trunk flexion (T1) on the model-predicted spinal loads was investigated. Mean of peak voluntary total...