Home About us Contact
|
Home About us Contact | ||
|
|
||
Element Analysis (element + analysis)
Kinds of Element Analysis Selected AbstractsSTRESS REDUCTION IN A GEAR TOOTH USING PHOTOELASTICITY AND FINITE ELEMENT ANALYSISEXPERIMENTAL TECHNIQUES, Issue 5 2001A. Ahlqvist No abstract is available for this article. [source] Recent Developments in Trace Element Analysis by ICP-AES and ICP-MS with Particular Reference to Geological and Environmental SamplesGEOSTANDARDS & GEOANALYTICAL RESEARCH, Issue 1 2005Kathryn L. Linge This review describes recent developments in trace element analysis using inductively coupled plasma-atomic emission spectrometry (ICP-AES) and inductively coupled plasma-mass spectrometry (ICP-MS). It aims to focus on the application of ICP techniques to geological and environmental samples. Therefore, fundamental studies in ICP-MS and ICP-AES instrumentation have largely been ignored. Whereas the majority of literature reviewed related to ICP-MS, indicating that ICP-MS is now the preferred technique for all geological analysis, there is still a steady development of ICP-AES to environmental applications. It is clear that true flexibility in elemental analysis can only be achieved by combining the advantages of both ICP-AES and ICP-MS. Two particular groups of elements (long-lived radionuclide and the platinum-group elements) stood out as warranting dedicated sections describing analytical developments these areas. [source] Surgical correction of scoliosis: Numerical analysis and optimization of the procedureINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 9 2010J. F. Aguilar Madeira Abstract A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ,, n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed. Copyright © 2010 John Wiley & Sons, Ltd. [source] Finite Element Analysis of the Proximal Femur and Hip Fracture Risk in Older Men,JOURNAL OF BONE AND MINERAL RESEARCH, Issue 3 2009Eric S Orwoll Abstract Low areal BMD (aBMD) is associated with increased risk of hip fracture, but many hip fractures occur in persons without low aBMD. Finite element (FE) analysis of QCT scans provides a measure of hip strength. We studied the association of FE measures with risk of hip fracture in older men. A prospective case-cohort study of all first hip fractures (n = 40) and a random sample (n = 210) of nonfracture cases from 3549 community-dwelling men ,65 yr of age used baseline QCT scans of the hip (mean follow-up, 5.6 yr). Analyses included FE measures of strength and load-to-strength ratio and BMD by DXA. Hazard ratios (HRs) for hip fracture were estimated with proportional hazards regression. Both femoral strength (HR per SD change = 13.1; 95% CI: 3.9,43.5) and the load-to-strength ratio (HR = 4.0; 95% CI: 2.7,6.0) were strongly associated with hip fracture risk, as was aBMD as measured by DXA (HR = 5.1; 95% CI: 2.8,9.2). After adjusting for age, BMI, and study site, the associations remained significant (femoral strength HR = 6.5, 95% CI: 2.3,18.3; load-to-strength ratio HR = 4.3, 95% CI: 2.5,7.4; aBMD HR = 4.4, 95% CI: 2.1,9.1). When adjusted additionally for aBMD, the load-to-strength ratio remained significantly associated with fracture (HR = 3.1, 95% CI: 1.6,6.1). These results provide insight into hip fracture etiology and demonstrate the ability of FE-based biomechanical analysis of QCT scans to prospectively predict hip fractures in men. [source] Structure, Orientation and Finite Element Analysis of the Tail Club of Mamenchisaurus hochuanensisACTA GEOLOGICA SINICA (ENGLISH EDITION), Issue 6 2009XING Lida Abstract: The structure and orientation of the posterior extremity (tail club) of the caudal vertebrae of Mamenchisaurus hochuanensis Young and Chao, 1972 from the Upper Jurassic Shangshaximiao Formation has been analyzed to determine the tail club function using Finite Element Analysis. Of the four caudal vertebrae composing the tail club, the second largest (C"1") was probably the most proximal, and is fixed with the preceding sequence of the caudal vertebrae, whereas the smallest (C"4") is free and forms the termination of the tail club. Our analysis also suggests that the tail club is more efficient in lateral swinging rather than up-and-down motion, and that the best region for the tail club to impact is at the spine of the largest of the four caudals (C"2"), with a maximum load for impact at about 450 N. The tail club of Mamenchisaurus hochuanensis probably also had limitations as a defense weapon and was more possibly a sensory organ to improve nerve conduction velocity to enhance the capacity for sensory perception of its surroundings. [source] Micromotion and Stress Distribution of Immediate Loaded Implants: A Finite Element AnalysisCLINICAL IMPLANT DENTISTRY AND RELATED RESEARCH, Issue 4 2009A. Fazel DDS ABSTRACT Background: Primary stability and micromotion of the implant fixture is mostly influenced by its macrodesign. Purpose: To assess and compare the peri-implant stress distribution and micromotion of two types of immediate loading implants, immediate loaded screw (ILS) Nisastan and Xive (DENTSPLY/Friadent, Monnheim, Germany), and to determine the best macrodesign of these two implants by finite element analysis. Methods: In this experimental study, the accurate pictures of two fixtures (ILS: height = 13, diameter = 4 mm and Xive: height = 13, diameter = 3.8 mm) were taken by a new digital camera (Nikon Coolpix 5700 [Nikon, Japan], resolution = 5.24 megapixel, lens = 8× optical, 4× digital zoom). Following accurate measurements, the three-dimensional finite element computer model was simulated and inserted in simulated mandibular bone (D2) in SolidWorks 2003 (SolidWork Corp., MA, USA) and Ansys 7.1 (Ansys, Inc., Canonsburg, PA, USA). After loading (500 N, 75° above horizon), the displacement was displayed and von Mises stress was recorded. Results: It was found that the primary stability of ILS was greater (152 µm) than Xive (284 µm). ILS exhibited more favorable stress distribution. Maximum stress concentration found in periapical bone around Xive (,30 MPa) was lesser than Nisastan (,37 MPa). Conclusions: Macrodesign of ILS leads to better primary stability and stress distribution. Maximum stress around Xive was less. [source] Global and local linear buckling behavior of a chiral cellular structurePHYSICA STATUS SOLIDI (B) BASIC SOLID STATE PHYSICS, Issue 3 2005A. Spadoni Abstract This paper investigates the flat-wise compression behavior of an innovative cellular structure configuration. The considered layout has a hexagonal chiral geometry featuring cylinders, or nodes, joined by ligaments, or ribs. The resulting assembly is characterized by a number of interesting properties that can be exploited for the design of alternative honeycombs or cellular topologies to be used in sandwich construction. The flat-wise strength of the chiral geometry is investigated through classical analytical formulas for the linear buckling of thin plates and shells and a bifurcation analysis performed on a Finite Element model. The analytical expressions predict the global buckling behavior and the resulting critical loads, and can be directly compared with the results obtained from the Finite Element analysis. In addition, the Finite Element model predicts local buckling modes, which should be considered to evaluate the possible development of localized plasticity. A sensitivity study is performed to evaluate the influence of the geometry of the chiral structure on its buckling strength. The study shows that the considered topology can offer great design flexibility, whereby several parameters can be selected and modified to improve the flat-wise performance. The comparison with traditional, hexagonal centro-symmetric structural configurations concludes the paper and demonstrates the enhanced performance and the potentials of chiral noncentro-symmetric designs. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source] Weathering the escarpment: chemical and physical rates and processes, south-eastern AustraliaEARTH SURFACE PROCESSES AND LANDFORMS, Issue 6 2009Benjamin C. Burke Abstract Differences in chemical weathering extent and character are expected to exist across topographic escarpments due to spatial gradients of climatic and/or tectonic forcing. The passive margin escarpment of south-eastern Australia has a debated but generally accepted model of propagation in which it retreated (within 40 Ma) to near its current position following rifting between Australia and New Zealand 85,100 Ma before present. We focus on this escarpment to quantify chemical weathering rates and processes and how they may provide insight into scarp evolution and retreat. We compare chemical weathering extents and rates above and below the escarpment using a mass balance approach coupling major and trace element analyses with previous measurements of denudation rates using cosmogenic nuclides (10Be and 26Al). We find a slight gradient in saprolite chemical weathering rate as a percentage of total weathering rate across the escarpment. The lowlands area, encompassing the region extending from the base of the escarpment to the coast, experiences a greater extent of chemical weathering than the highland region above the escarpment. Percents of denudation attributable to saprolite weathering average 57 ± 6% and 47 ± 7% at low and high sites respectively. Furthermore, the chemical index of alteration (CIA), a ratio of immobile to mobile oxides in granitic material that increases with weathering extent, have corresponding average values of 73·7 ± 3·9 and 65·5 ± 3·4, indicating lower extents of weathering above the escarpment. Finally, we quantify variations in the rates and extent of chemical weathering at the hillslope scale across the escarpment to suggest new insight into how climate differences and hillslope topography help drive landscape evolution, potentially overprinting longer term tectonic forcing. Copyright © 2009 John Wiley & Sons, Ltd. [source] An efficient three-dimensional solid finite element dynamic analysis of reinforced concrete structuresEARTHQUAKE ENGINEERING AND STRUCTURAL DYNAMICS, Issue 2 2006K. V. Spiliopoulos Abstract Most of the finite element analyses of reinforced concrete structures are restricted to two-dimensional elements. Three-dimensional solid elements have rarely been used although nearly all reinforced concrete structures are under a triaxial stress state. In this work, a three-dimensional solid element based on a smeared fixed crack model that has been used in the past mainly for monotonic static loading analysis is extended to cater for dynamic analysis. The only material parameter that needs to be input for this model is the uniaxial compressive strength of concrete. Steel bars are modelled as uniaxial elements and an embedded formulation allows them to have any orientation inside the concrete elements. The proposed strategy for loading or unloading renders a numerical procedure which is stable and efficient. The whole process is applied to two RC frames and compared against existing experiments in the literature. Results show that the proposed approach may adequately be used to predict the dynamic response of a structure. Copyright © 2005 John Wiley & Sons, Ltd. [source] A parallel multigrid solver for high-frequency electromagnetic field analyses with small-scale PC clusterELECTRONICS & COMMUNICATIONS IN JAPAN, Issue 9 2008Kuniaki Yosui Abstract Finite element analyses of electromagnetic fields are commonly used for designing various electronic devices. The scale of the analyses becomes larger and larger, therefore, a fast linear solver is needed to solve linear equations arising from the finite element method. Since a multigrid solver is the fastest linear solver for these problems, parallelization of a multigrid solver is quite a useful approach. From the viewpoint of industrial applications, an effective usage of a small-scale PC cluster is important due to initial cost for introducing parallel computers. In this paper, a distributed parallel multigrid solver for a small-scale PC cluster is developed. In high-frequency electromagnetic analyses, a special block Gauss, Seidel smoother is used for the multigrid solver instead of general smoothers such as a Gauss, Seidel or Jacobi smoother in order to improve the convergence rate. The block multicolor ordering technique is applied to parallelize the smoother. A numerical example shows that a 3.7-fold speed-up in computational time and a 3.0-fold increase in the scale of the analysis were attained when the number of CPUs was increased from one to five. © 2009 Wiley Periodicals, Inc. Electron Comm Jpn, 91(9): 28, 36, 2008; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/ecj.10160 [source] Stress-life fatigue assessment of pipelines with plain dentsFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 12 2009S. B. CUNHA ABSTRACT This paper presents a new algorithm for assessing the fatigue life of dented pipelines. The proposed methodology was conceived according to the current stress-life fatigue theory and design practice: it employs S,N curves inferred from tensile test material properties and uses well established methodologies to deal with the stress concentration, the mean stress and the multi-axial stress state that characterizes a dented pipe. Finite element analyses are carried out to model the denting process and to determine the stress concentration factors of several pipe-dent geometries. Using dimensional analysis over the numerical results, a non-dimensional number to characterize the pipe-dent geometry is determined and linear interpolation expressions for the stress concentration factors of dented pipelines are developed. Fatigue tests are conducted with the application of cyclic internal pressure on small-scale dented steel pipe models. In view of the fatigue test results, the more appropriate S,N curve and mean stress criteria are selected. [source] Numerical simulation of rotating bending process for U-tubes in heat exchangersFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 10 2009H.-S. KIM ABSTRACT Heat exchangers comprise thousands of tubes having U-shaped portions. Rotating bending method has been widely utilized to make U-bends. Although this method shows an excellent performance, cracks have been frequently detected in the U-bends due to residual stresses induced by bending. In this paper, the bending process is simulated based on elastic,plastic finite element analyses in order to investigate the magnitude and distribution of the residual stresses including the effects of operating pressure. Analyses results show that the residual stress increases as the radius of U-bend decreases and that operating pressure has a detrimental effect in terms of stress corrosion cracking at the intrados of U-bend. It is thought that these results can be utilized for the estimations of fracture mechanics parameters such as limit load, stress intensity factor and J-integral, prevention of the cracking, and establishment of the optimum inspection strategy for the heat exchanger tubes. [source] Significance of the elastic peak stress evaluated by FE analyses at the point of singularity of sharp V-notched componentsFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 2 2007G. MENEGHETTI ABSTRACT The paper presents an expression useful to estimate the notch stress intensity factor (NSIF) from finite element analyses carried out by using a mesh pattern with a constant element size. The evaluation of the NSIF from a numerical analysis of the local stress field usually requires very refined meshes and then large computational effort. The usefulness of the presented expression is that (i) only the elastic peak stress numerically evaluated at the V-notch tip is needed and no longer the whole stress,distance set of data; (ii) the adopted meshes are rather coarse if compared to those necessary for the evaluation of the whole local stress field. The proposed expression needs the evaluation of a virtual V-notch tip radius, i.e. the radius which would produce the same elastic peak stress than that calculated by FEM at the sharp V-notch tip by means of a given mesh pattern. Once such a radius has been theoretically determined for a given geometry, the expression can be applied in a wide range of notch depths and opening angles. [source] ,Arching' effect in elastic polycrystals: implications for the variability of fatigue livesFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 4 2002S. Pommier ABSTRACT The paper deals with a study of heterogeneous stress and strain distribution in polycrystals in relation with elastic anisotropy of grains. A similitude with the arching effect widely studied in granular materials is proposed and this concept is extended to heterogeneous polycrystals in which the load transfer is not binary in the way it is in granular media but may vary significantly and suddenly from one grain to another according to the crystal orientation to the load direction. Experiments and 3D finite element analyses show that though the individual orientation of grains is random, the strain and stress distribution is not. A network is formed inside the polycrystal whose scale is larger than the grain size. The load percolation network consists in heavily loaded links whose direction is coincident with the direction of the principal stresses. So, the typical scale for the variability of the local stresses is not the grain size but the size of the load percolation network. Since this scale is found to be rather large in particular for iron, zinc and copper, this effect should contribute significantly to the variability of the fatigue lives of notched vs. smooth components. [source] A shear stress-based parameter for fretting fatigue crack initiationFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 7 2001C. D. Lykins The purpose of this study was to investigate the fretting fatigue crack initiation behaviour of titanium alloy, Ti,6Al,4V. Fretting contact conditions were varied by using different geometries of the fretting pad. Applied forces were also varied to obtain fretting fatigue crack initiation lives in both the low- and high-cycle fatigue regimes. Fretting fatigue specimens were examined to determine the crack location and the crack angle orientation along the contact surface. Salient features of fretting fatigue experiments were modelled and analysed with finite element analysis. Computed results of the finite element analyses were used to formulate a shear stress-based parameter to predict the fretting fatigue crack initiation life, location and orientation. Comparison of the analytical and experimental results showed that fretting fatigue crack initiation was governed by the maximum shear stress, and therefore a parameter involving the maximum shear stress range on the critical plane with the correction factor for the local mean stress or stress ratio effect was found to be effective in characterizing the fretting fatigue crack initiation behaviour in titanium alloy, Ti,6Al,4V. [source] Multiaxial fatigue of welded joints under constant and variable amplitude loadingsFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 5 2001C. M. Sonsino Flange-tube joints from fine grained steel StE 460 with unmachined welds were investigated under biaxial constant and variable amplitude loading (bending and torsion) in the range of 103 to 5,×,106 cycles to crack initiation and break-through, respectively. In order not to interfere with residual stresses they were relieved by a heat treatment. In-phase loading can be treated fairly well using the conventional hypotheses (von Mises or Tresca) on the basis of nominal, structural or local strains or stresses. But the influence of out-of-phase loading on fatigue life is severely overestimated if conventional hypotheses are used. However, the hypothesis of the effective equivalent stress that is introduced leads to fairly good predictions for constant as well as for random variable amplitude loads. Therefore, the knowledge of local strains or stresses is necessary. They are determined by boundary element analyses that are dependent on weld geometry. This hypothesis considers the fatigue-life-reducing influence of out-of-phase loading by taking into account the interaction of local shear stresses acting in different surface planes of the material. Further, size effects resulting from weld geometry and loading mode were included. Damage accumulation under a Gaussian spectrum can be assessed for in- and out-of-phase combined bending and torsion using an allowable damage sum of 0.35. [source] Ductile fracture of commercial purity titanium at room temperatureFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 10 2000G. Shatil An experimental and numerical program was carried out to examine and assess the deformation and fracture behaviour of alloys of commercial purity (CP) titanium. The material rate-dependent deformation under constant displacement rates and under sustained loads was directly simulated in finite element analyses using an implemented unified material model. The simulations predicted the fracture of compact tension specimens subjected to J,R tests and sustained load tests employing a dimensional analysis and strain-hardening approach. Differences between two batches with different oxygen contents were examined and the limitation of the material model was investigated. [source] Bonded aircraft repairs under variable amplitude fatigue loading and at low temperaturesFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 1 2000Vlot Bonded repairs can replace mechanically fastened repairs for aircraft structures. Compared to mechanical fastening, adhesive bonding provides a more uniform and efficient load transfer into the patch, and can reduce the risk of high stress concentrations caused by additional fastener holes necessary for riveted repairs. Previous fatigue tests on bonded Glare (glass-reinforced aluminium laminate) repairs were performed at room temperature and under constant amplitude fatigue loading. However, the realistic operating temperature of ,40 °C may degrade the material and will cause unfavourable thermal stresses. Bonded repair specimens were tested at ,40 °C and other specimens were tested at room temperature after subjecting them to temperature cycles. Also, tests were performed with a realistic C-5A Galaxy fuselage fatigue spectrum at room temperature. The behaviour of Glare repair patches was compared with boron/epoxy ones with equal extensional stiffness. The thermal cycles before fatigue cycling did not degrade the repair. A constant temperature of ,40 °C during the mechanical fatigue load had a favourable effect on the fatigue crack growth rate. Glare repair patches showed lower crack growth rates than boron/epoxy repairs. Finite element analyses revealed that the higher crack growth rates for boron/epoxy repairs are caused by the higher thermal stresses induced by the curing of the adhesive. The fatigue crack growth rate under spectrum loading could be accurately predicted with stress intensity factors calculated by finite element modelling and cycle-by-cycle integration that neglected interaction effects of the different stress amplitudes, which is possible because stress intensities at the crack tip under the repair patch remain small. For an accurate prediction it was necessary to use an effective stress intensity factor that is a function of the stress ratio at the crack tip Rcrack tip including the thermal stress under the bonded patch. [source] Finite element analyses of layered visco-elastic system under vertical circular loadingINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 8 2008Q. Xu Abstract Analyses for the response of a linear visco-elastic system subjected to axi-symmetric vertical circular loading are presented. Hankel transforms with respect to the radial spatial coordinate are used to reduce the three-dimensional problem to that involving only a single spatial dimension, which is then discretized using the finite element method. Three techniques are employed to handle the time factor in the visco-elastic material: (i) direct time integration; (ii) Fourier transforms; and (iii) Laplace transforms. These methods are compared and evaluated through their numerical results. Copyright © 2007 John Wiley & Sons, Ltd. [source] Three-dimensional finite element analyses of passive pile behaviourINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 7 2006L. F. Miao Abstract Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction-related activities such as embankment construction or excavation operations. Three-dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20-node quadrilateral brick elements with reduced integration. For compatibility between the soil,pile interface elements, 27-node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil,pile interface. A Mohr,Coulomb elastic,plastic constitutive model with large-strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 John Wiley & Sons, Ltd. [source] Application of micropolar plasticity to post failure analysis in geomechanicsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 10 2004Majid T. Manzari Abstract A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy,Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd. [source] Curvature- and displacement-based finite element analyses of flexible slider crank mechanismsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 10 2010Y. L. Kuo Abstract The paper presents the applications of the curvature- and displacement-based finite element methods to flexible slider crank mechanisms. The displacement-based method usually needs more elements or high-degree polynomials to obtain highly accurate solutions. The curvature-based method assumes a polynomial to approximate a curvature distribution, and the expressions are investigated to obtain the displacement and rotation distributions. During the process, the boundary conditions associated with displacement, rotation, and curvature are imposed, which leads the great reduction of the number of degrees of freedom that are required. The numerical results demonstrate that the errors obtained by applying the curvature-based method are much smaller than those by applying the displacement-based method, based on the comparison of the same number of degrees of freedom. Copyright © 2008 John Wiley & Sons, Ltd. [source] Surgical correction of scoliosis: Numerical analysis and optimization of the procedureINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 9 2010J. F. Aguilar Madeira Abstract A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ,, n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed. Copyright © 2010 John Wiley & Sons, Ltd. [source] Nonlinear transient dynamic analysis by explicit finite element with iterative consistent mass matrixINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 3 2009Shen Rong Wu Abstract Various mass matrices in the explicit finite element analyses of nonlinear transient dynamic problems are investigated. The matrices are obtained as a linear combination of lumped and consistent mass matrices. An iterative procedure to calculate the inverse of the consistent and the mixed mass matrices in the framework of explicit finite element method is presented. The convergence of the iterative procedure is proved. The inverse of the consistent and mixed mass matrices is approximated by the iteration and is used to compare the results from the lumped mass matrix. For the impact of a structural component and a vehicle, some difference in the results by using coarse mesh is observed. For the component using fine mesh, no significant difference is found. Copyright © 2008 John Wiley & Sons, Ltd. [source] Comparison between cohesive zone modelsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 11 2004K. Y. Volokh Cohesive zone models (CZMs) are widely used for numerical simulation of the fracture process. Cohesive zones are surfaces of discontinuities where displacements jump. A specific constitutive law relating the displacement jumps and proper tractions defines the cohesive zone model. Within the cohesive zone approach crack nucleation, propagation, and arrest are a natural outcome of the theory. The latter is in contrast to the traditional approach of fracture mechanics where stress analysis is separated from a description of the actual process of material failure. The common wisdom says that only cohesive strength,the maximum stress on the traction,separation curve,and the separation work,the area under the traction,separation curve,are important in setting a CZM while the shape of the traction,separation curve is subsidiary. It is shown in our note that this rule may not be correct and a specific shape of the cohesive zone model can significantly affect results of the fracture analysis. For this purpose four different cohesive zone models,bilinear, parabolic, sinusoidal, and exponential,are compared by using a block-peel test, which allows for simple analytical solutions. Numerical performance of the cohesive zone models is considered. It appears that the convergence properties of nonlinear finite element analyses are similar for all four CZMs in the case of the block-peel test. Copyright © 2004 John Wiley & Sons, Ltd. [source] Guaranteed computable error bounds for conforming and nonconforming finite element analyses in planar elasticityINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2010Mark Ainsworth Abstract We obtain fully computable a posteriori error estimators for the energy norm of the error in second-order conforming and nonconforming finite element approximations in planar elasticity. These estimators are completely free of unknown constants and give a guaranteed numerical upper bound on the norm of the error. The estimators are shown to also provide local lower bounds, up to a constant and higher-order data oscillation terms. Numerical examples are presented illustrating the theory and confirming the effectiveness of the estimator. Copyright © 2009 John Wiley & Sons, Ltd. [source] The stability of stars of triangular equilibrium plate elementsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 7 2009E. A. W. Maunder Abstract Equilibrium models for finite element analyses are becoming increasingly important in complementary roles to those from conventional conforming models, but when formulating equilibrium models questions of stability, or admissibility of loads, are of major concern. This paper addresses these questions in the context of flat plates modelled with triangular hybrid elements involving membrane and/or flexural actions. Patches of elements that share a common vertex are considered, and such patches are termed stars. Stars may be used in global analyses as assemblies of elements forming macro-elements, or in local analyses. The conditions for stability, or the existence and number of spurious kinematic modes, are determined in a general algebraic procedure for any degree of the interpolation polynomials and for any geometric configuration. The procedure involves the determination of the rank of a compatibility matrix by its transformation to row echelon form. Examples are presented to illustrate some of the characteristics of spurious kinematic modes when they exist in stars with open or closed links. Copyright © 2008 John Wiley & Sons, Ltd. [source] A random field model for generating synthetic microstructures of functionally graded materialsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 7 2008Sharif Rahman Abstract This article presents a new level-cut, inhomogeneous, filtered Poisson random field model for representing two-phase microstructures of statistically inhomogeneous, functionally graded materials with fully penetrable embedded particles. The model involves an inhomogeneous, filtered Poisson random field comprising a sum of deterministic kernel functions that are scaled by random variables and a cut of the filtered Poisson field above a specified level. The resulting level-cut field depends on the Poisson intensity, level, kernel functions, random scaling variables, and random rotation matrices. A reconstruction algorithm including model calibration and Monte Carlo simulation is presented for generating samples of two-phase microstructures of statistically inhomogeneous media. Numerical examples demonstrate that the model developed is capable of producing a wide variety of two- and three-dimensional microstructures of functionally graded composites containing particles of various sizes, shapes, densities, gradations, and orientations. An example involving finite element analyses of random microstructures, leading to statistics of effective properties of functionally graded composites, illustrates the usefulness of the proposed model. Copyright © 2008 John Wiley & Sons, Ltd. [source] Locking-free finite elements for shear deformable orthotropic thin-walled beamsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 7 2007F. Minghini Abstract Numerical models for finite element analyses of assemblages of thin-walled open-section profiles are presented. The assumed kinematical model is based on Timoshenko,Reissner theory so as to take shear strain effects of non-uniform bending and torsion into account. Hence, strain elastic-energy coupling terms arise between bending in the two principal planes and between bending and torsion. The adopted model holds for both isotropic and orthotropic beams. Several displacement interpolation fields are compared with the available numerical examples. In particular, some shape functions are obtained from ,modified' Hermitian polynomials that produce a locking-free Timoshenko beam element. Analogously, numerical interpolation for torsional rotation and cross-section warping are proposed resorting to one Hermitian and six Lagrangian formulation. Analyses of beams with mono-symmetric and non-symmetric cross-sections are performed to verify convergence rate and accuracy of the proposed formulations, especially in the presence of coupling terms due to shear deformations, pointing out the decay length of end effects. Profiles made of both isotropic and fibre-reinforced plastic materials are considered. The presented beam models are compared with results given by plate-shell models. Copyright © 2007 John Wiley & Sons, Ltd. [source] Localized remeshing techniques for three-dimensional metal forming simulations with linear tetrahedral elementsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2006Il-Heon Son Abstract The localized remeshing technique for three-dimensional metal forming simulations is proposed based on a mixed finite element formulation with linear tetrahedral elements in the present study. The numerical algorithm to generate linear tetrahedral elements is developed for finite element analyses using the advancing front technique with local optimization method which keeps the advancing fronts smooth. The surface mesh generation using mesh manipulations of the boundary elements of the old mesh system was made to improve mesh quality of the boundary surface elements, resulting in reduction of volume change in forming simulations. The mesh quality generated was compared with that obtained from the commercial CAD package for the complex geometry like lumbar. The simulation results of backward extrusion and bevel gear and spider forgings indicate that the currently developed simulation technique with the localized remeshing can be used effectively to simulate the three-dimensional forming processes with a reduced computation time. Copyright © 2006 John Wiley & Sons, Ltd. [source] | |||||||||||||||||||||||||||||||