Sharif Digital Repository

In this paper, an adaptive finite element analysis is presented for 3D modeling of non-planar curved crack growth. The fracture mechanical evaluation is performed based on a general technique for non-planar curved cracks. The Schollmann's crack kinking criterion is used for the process of crack propagation in 3D problems. The Zienkiewicz-Zhu error estimator is employed in conjunction with a weighted SPR technique at each patch to improve the accuracy of error estimation. Applying the proposed technique to 3D non-planar curved crack growth problems shows significant improvements particularly at the boundaries and near crack tip regions. Several numerical examples are presented to illustrate... 

Cell survival during the early stages of transplantation and before new blood vessels formation is a major challenge in translational applications of 3D bioprinted tissues. Supplementing oxygen (O2) to transplanted cells via an O2 generating source such as calcium peroxide (CPO) is an attractive approach to ensure cell viability. Calcium peroxide also produces calcium hydroxide that reduces the viscosity of bioinks, which is a limiting factor for bioprinting. Therefore, adapting this solution into 3D bioprinting is of significant importance. In this study, a gelatin methacryloyl (GelMA) bioink that is optimized in terms of pH and viscosity is developed. The improved rheological properties... 

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.... 

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... 

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... 

Alumina is widely used as a vacuum ceramic separator in the pillbox-type output window of high-power electron tubes. The fracture of alumina is one of the most important factors in the failure of these tubes. In this article, the fracture process of alumina in the pillbox-type window of high-power electron tubes is assessed microscopically. The investigation shows the twinning phenomenon in the fractured alumina. The reason for initiating and propagating twins is also investigated. An increase in the alumina temperature decreases the required stress level for generating twins. The overheating condition occurs due to an impedance mismatch in the pillbox-type window. Experimental and numerical... 

A mathematical model is proposed for the determination of the thermo-mechanical stresses in work-rolls during hot-strip rolling. The model describes the evolution of the temperature fields in the work-roll and in the work-piece, and in the latter the plastic heat generation is taken into account. The frictional heat generated on the contact surface is also included. The problem is treated in two steps. First, a numerical method is developed for the analysis of the coupled thermal problem of the temperature distributions within the work-roll and metal being rolled, while an admissible velocity field is employed to estimate the heat of deformation. In the second step, the finite-element method... 

A modified two-surface critical state plasticity model for saturated and unsaturated soil is presented in this study. The key modification in new model is inclusion of an alternative yield surface used to simulate the behavior of unsaturated soils in addition to corresponding saturated conditions. Moreover, a numerical technique is used to obtain an incremental stress–strain response from loading curves. Modification is applied continuously in each incremental step to return the final stress states and hardening parameters to the yield surface. Results revealed that the adopted modeling approach can predict two independent sets of laboratory unsaturated experiments under various conditions... 

This paper aims to identify an optimum bone fracture stabilizer. For this purpose, three design variables including the ratio of the screw diameter to the plate width at three levels, the ratio of the plate thickness to the plate width at three levels, and the diameter of the bone at two levels were selected for analysis. Eighteen 3D verified finite element models were developed to examine the effects of these parameters on the weight, maximum displacement and maximum von Mises stress of the fixation structure. Considering the relations between the inputs and outputs using multivariate regression, a genetic algorithm was used to find the optimal choices. Results showed that the diameter of... 

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  

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... 

Based on the force-deflection equation for a beam subjected to lateral point loads, a C2 continuous piecewise bicubic mathematical representation was proposed to model complicated geometrical surfaces, e.g. the articular surfaces of human joints. The method was then extended so that it could be used for mathematical modelling of incomplete nets of data points, as well as smoothing of noisy and/or filtering of erroneous data points. Mathematical techniques were also developed to calculate the required unknown parameters explicitly, with no need to solve the system of equations simultaneously. The performance of the proposed method was evaluated on a number of surface modelling problems,... 

In this paper an integrated model has been developed that is capable of considering the effects of process parameters on the thermo-mechanical behaviour of the work-roll and the strip during hot rolling operations. The commercially finite element package ABAQUS is utilized to solve the governing equations, while an elasto-viscoplastic material behaviour is utilized in the model. The model may be employed in single-pass as well as multi-pass rolling layouts. Various aspects of the hot rolling process are predicted including temperature field and thermo-mechanical stresses within the work-roll as well as temperature, strain, and residual stresses distributions within the strip. To assess the... 

Patellar lateral stability was studied using a 2D transverse plane model with deformable articular surfaces. Quadriceps muscles and patellar tendon were considered as strings with predefined forces and lateral and medial retinaculum as tensile springs. Deformation behavior of articular cartilage was modeled by a set of compression springs perpendicular to articular surfaces, based on rigid body spring model method (RBSM). Patellar lateral stability was investigated using restraining force method (the external force required to cause up to 10 mm lateral displacement on patella). The results were in good agreement with experimental reports for normal joint, vastus lateralis and vastus medialis... 

Atherosclerosis as a common cardiovascular disease is a result of both adverse hemodynamics conditions and monocyte deposition within coronary arteries. It is known that the adhesion of monocytes on the arterial wall and their interaction with the vascular surface are one of the main parameters in the initiation and progression of atherosclerosis. In this work, hemodynamic parameters and monocyte deposition have been investigated in a 3D computational model of the Left Anterior Descending coronary artery (LAD) and its first diagonal branch (D1) under the heart motion. A one-way Lagrangian approach is performed to trace the monocyte particles under different blood flow regimes and heart motion... 

This study presents a novel statistical volume element (SVE) for micromechanical modeling of the white matter structures, with histology-informed randomized distribution of axonal tracts within the extracellular matrix. The model was constructed based on the probability distribution functions obtained from the results of diffusion tensor imaging as well as the histological observations of scanning electron micrograph, at two structures of white matter susceptible to traumatic brain injury, i.e. corpus callosum and corona radiata. A simplistic representative volume element (RVE) with symmetrical arrangement of fully alligned axonal fibers was also created as a reference for comparison. A... 

Background: The development of endoscopic sinus surgery (ESS) training simulators for clinical environment applications has reduced the existing shortcomings in conventional teaching methods, creating a standard environment for trainers and trainees in a more accurate and repeatable fashion. Materials and methods: In this research, the validation study of an ESS training simulator has been addressed. It is important to consider components that guide trainees to improve their hand movements control in the orbital floor removal in an ESS operation. Therefore, we defined three tasks to perform: pre-experiment learning, training, and evaluation. In these tasks, the critical regions introduced in... 

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... 

In this paper, some basis-free expressions for the material time derivative of Lagrangian stress tensors are presented which are generally valid in all cases of coalescent principal stretches. The material is assumed to be elastic and isotropic  

If the double contraction of a stress tensor such as T and rate of a Lagrangean strain tensor such as E, i.e. T : Ė produces the stress power then these stress and strain tensors are called a conjugate pair. The applications of the conjugate stress and strain measures are in the development of the basic relations in nonlinear continuum mechanics analysis such as modeling of constitutive equations of elastic-plastic materials. In this paper relations for stress tensors conjugate to an arbitrary Lagrangean strain measure of Hill's class are obtained. The results of this paper are more compact and simpler in compare with those available in the literature. The results are valid for the three...