Sharif Digital Repository

Controlling topographic features at all length scales is of great importance for the interaction of cells with tissue regenerative materials. We utilized an indirect three-dimensional printing method to fabricate polymeric scaffolds with pre-defined and controlled external and internal architecture that had an interconnected structure with macro- (400-500 μm) and micro- (∼25 μm) porosity. Polycaprolactone (PCL) was used as model system to study the kinetics of tissue growth within porous scaffolds. The surface of the scaffolds was decorated with TiO2 and bioactive glass (BG) nanoparticles to the better match to nanoarchitecture of extracellular matrix (ECM). Micrometric BG particles were... 

In this study, a nonlinear poroelastic model of intervertebral disc as an infrastructure was developed. Moreover, a new element was defined consisting a disc (Viscoelastic Euler Beam Element) and a vertebra (Rigid Link) as a unit element. Using the new element, three different viscoelastic finite element models were prepared for lumbar motion segment (L4/L5). Prolonged loading (short-term and long-term creep) and cyclic loading were applied to the models and the results were compared with results of in vivo tests. Simplification of the models by using the new element leads to reduction of the runtime of the models in dynamic analyses to few minutes without losing the accuracy in the results  

Finite element analysis is an effective tool to evaluate the material properties of living tissue. For an interactive optimization procedure, the finite element analysis usually needs many simulations to reach a reasonable solution. The metamodel analysis of finite element simulation can be used to reduce the computation of a structure with complex geometry or a material with composite constitutive equations. The intervertebral disc is a complex, heterogeneous, and hydrated porous structure. A poroelastic finite element model can be used to observe the fluid transferring, pressure deviation, and other properties within the disc. Defining reasonable poroelastic material properties of the... 

Two artificial neural networks (ANNs) are constructed, trained, and tested to map inputs of a complex trunk finite element (FE) model to its outputs for spinal loads and muscle forces. Five input variables (thorax flexion angle, load magnitude, its anterior and lateral positions, load handling technique, i.e., one- or two-handed static lifting) and four model outputs (L4-L5 and L5-S1 disc compression and anterior-posterior shear forces) for spinal loads and 76 model outputs (forces in individual trunk muscles) are considered. Moreover, full quadratic regression equations mapping input-outputs of the model developed here for muscle forces and previously for spine loads are used to compare the... 

An analytical approach which is popular in micromechanical studies has been extended to the solution for the interference fit problem of the femoral stem in cementless total hip arthroplasty (THA). The multiple inhomogeneity problem of THA in transverse plane, including an elliptical stem, a cortical wall, and a cancellous layer interface, was formulated using the equivalent inclusion method (EIM) to obtain the induced interference elastic fields. Results indicated a maximum interference fit of about 210 μm before bone fracture, predicted based on the Drucker-Prager criterion for a partially reamed section. The cancellous layer had a significant effect on reducing the hoop stresses in the... 

This paper introduces three geometric features, from deformed shape of a soft tissue, which demonstrate good correlation with probing force and maximum local stress. Using FEM simulation, 2D and 3D model of an in vivo porcine liver was built for different probing tasks. Maximum deformation angle, maximum deformation depth and width of displacement constraint of the reconstructed shape of the deformed body were calculated. Two neural networks were trained from these features and the calculated interaction forces. The features are shown to have high potential to provide force estimation either for haptic devices or to assess the damage to the tissue in large deformations of up to 40%  

The aim of this study was to use finite element analysis (FEA) to assess the influence of microthread design at the implant neck on stress distribution in the surrounding bone. A commercially available implant with 3.5 mm diameter and 10.5 mm length was selected and used as a model. For the purpose of designing the microthread implant model, microthreads were added to the implant neck in a computerized model. A force measuring 100 N was then applied to the entire surface of the abutment in the vertical direction. The results showed that in both models, stress was mainly concentrated at the cortical bone adjacent to the neck of the implant. Maximum stress values in the cortical bone... 

The performance of cellular solids in biomedical applications relies strongly on a detailed understanding of the effects of pore topology on mechanical properties. This study aims at characterizing the failure mechanism of scaffolds based on nodal connectivity (number of struts that meet in joints) and geometry of the pores. Plastic models of scaffolds having the same relative density but different cubic and trigonal unit cells were designed and then fabricated via three dimensional (3-D) printing. Unit cells were repeated in different arrangements in 3-D space. An in-situ imaging technique was utilized to study the progressive deformation of the scaffold models. Different nodal... 

Background: Asymmetric lifting activities are associated with low back pain. Methods: A finite element biomechanical model is used to estimate spinal loads during one- and two-handed asymmetric static lifting activities. Model input variables are thorax flexion angle, load magnitude as well as load sagittal and lateral positions while response variables are L4-L5 and L5-S1 disc compression and shear forces. A number of levels are considered for each input variable and all their possible combinations are introduced into the model. Robust yet user-friendly predictive equations that relate model responses to its inputs are established. Findings: Predictive equations with adequate... 

The correlations between the hardness, yield stress, and modulus of elasticity of isotactic polypropylene (iPP) were evaluated on the local and global scales. Nanoindentation and traditional macromechanical tests were incorporated for this purpose. Thus, local and global mechanical properties were measured at various temperatures and strain rates. A certain relation was found between the local and global mechanical properties. Moreover, Johnson's model (developed according to the expanding cavity model) was also evaluated at various temperatures and strain rates. The Johnson model was valid only for the yield stresses obtained by nanoindentation and compressive tests and also the elastic... 

Muscle force partitioning methods and musculoskeletal system simplifications are key modeling issues that can alter outcomes, and thus change conclusions and recommendations addressed to health and safety professionals. A critical modeling concern is the use of single-joint equilibrium to estimate muscle forces and joint loads in a multi-joint system, an unjustified simplification made by most lumbar spine biomechanical models. In the context of common occupational tasks, an EMG-assisted optimization method (EMGAO) is modified in this study to simultaneously account for the equilibrium at all lumbar joints (M-EMGAO). The results of this improved approach were compared to those of its... 

Porous shape memory alloys (SMAs) exhibit the interesting characteristics of porous metals together with shape memory effect and pseudo-elasticity of SMAs that make them appropriate for biomedical applications. In this paper, a 3-D phenomenological constitutive model for the pseudo-elastic behavior and shape memory effect of porous SMAs is developed within the framework of irreversible thermodynamics. Comparing to micromechanical and computational models, the proposed model is computationally cost effective and predicts the behavior of porous SMAs under proportional and non-proportional multiaxial loadings. Considering the pressure dependency of phase transformation in porous SMAs, proper... 

Nafion membrane encounters many different thermal conditions and mechanical loadings because of its wide range of applications as a proton exchange membrane (PEM). Molecular dynamics simulation of hydrated Nafion at different temperatures is carried out to investigate the alteration of the physical properties of Nafion under uniaxial loading over a wide range of temperatures. According to the simulation results, increase of the temperature reduces the yield stress. The results also show that the polymer chains ordering increases the glass transition temperature and enhances the self-diffusion coefficient of water in hydrated Nafion. Comparisons show that the elastic modulus and viscosity... 

Axon is an important part of the neuronal cells and axonal microtubules are bundles in axons. In axons, microtubules are coated with microtubule-associated protein tau, a natively unfolded filamentous protein in the central nervous system. These proteins are responsible for cross-linking axonal microtubule bundles. Through complimentary dimerization with other tau proteins, bridges are formed between nearby microtubules creating bundles. Formation of bundles of microtubules causes their transverse reinforcement and has been shown to enhance their ability to bear compressive loads. Though microtubules are conventionally regarded as bearing compressive loads, in certain circumstances during... 

This study investigates the effect of body acceleration on human cardiac function. Finite element analysis is conducted to simulate geometrical and mechanical properties of human heart. Heart geometrical modeling in three-dimension is performed by segmentation of cardiac MRI images. The nonlinear mechanical behavior of myocardium is modeled by Mooney-Rivlin, Polynomial, Ogden and Yeoh hyperelastic material models. Stress-strain curves of myocardial tissue are obtained from experimental compression tests on bovine heart samples. The experimental results are employed for the evaluation of material coefficients by the nonlinear least squares method. Among hyperelastic models, the Yeoh model... 

In this research, the characteristics of the prepared samples in epoxy matrix by means of X-ray diffraction (XRD), energy dispersive X-ray spectroscopy (EDS), as well as scanning electron microscope (SEM) are evaluated. Meanwhile, the obtained mechanical properties of the specimen are investigated. Thermogravimetric analysis (TGA) is also employed to evaluate the thermal degradation of manufactured nanocomposites. The thermal neutron absorption properties of nanocomposites containing 3 wt% of montmorillonite nanoclay (closite30B) have been studied experimentally, using an Am-Be point source. Mechanical tests reveal that the higher B4C concentrations, the more tensile strengths, but lower... 

In the present study, we investigate the effect of the hemodynamic factors of the blood flow on the cerebral aneurysms. To this end, a hypothetical geometry of the aneurysm in the circle of Willis, located in the bifurcation point of the anterior cerebral artery (ACA) and anterior communicating artery (ACoA) is modeled in a three-dimensional manner. Three cases are chosen in the current study: an untreated thin wall (first case), untreated thick wall (second case), and a treated aneurysm (third case). The effect of increasing the aneurysm wall thickness on the deformation and stress distribution of the walls are studied. The obtained results showed that in the second case, a reduction in the... 

In this paper, a constitutive and micromechanical model for prediction of rate-dependent behavior of connective tissues (CTs) is presented. Connective tissues are considered as nonlinear viscoelastic material. The rate-dependent behavior of CTs is incorporated into model using the well-known quasi-linear viscoelasticity (QLV) theory. A planar wavy representative volume element (RVE) is considered based on the tissue microstructure histological evidences. The presented model parameters are identified based on the available experiments in the literature. The presented constitutive model introduced to ABAQUS by means of UMAT subroutine. Results show that, monotonic uniaxial test predictions of... 

Obstruction of left anterior descending artery (LAD) due to the thrombosis or atherosclerotic plaques is the leading cause of death worldwide. Targeted delivery of drugs through micro- and nanoparticles is a very promising approach for developing new strategies in clot-busting or treating restenosis. In this work, we modelled the blood flow characteristics in a patient-specific reconstructed LAD artery by the fluid–solid interaction method and based on physiological boundary conditions. Next, we provided a Lagrangian description of micro- and nanoparticles dynamics in the blood flow considering their Brownian motion and the particle–particle interactions. Our results state that the number of... 

The effect of the interface stresses is studied upon the size-dependent elastic deformation of an elastic half-plane having a cylindrical inclusion with distinct elastic properties. The elastic half-plane is subjected to either a uniaxial loading at infinity or a uniform non-shear eigenstrain in the inclusion. The straight edge of the half-plane is either traction-free, or rigid-slip, or motionless, which represents three practical situations of mechanical structures. Using two-dimensional Papkovich-Neuber potentials and the theory of surface/interface elasticity, the elastic field in the elastic half-plane is obtained. Comparable with classical result, the new formulation renders the...