A finite difference method for off-fault plasticity throughout the earthquake cycle

NASA Astrophysics Data System (ADS)

Erickson, Brittany A.; Dunham, Eric M.; Khosravifar, Arash

2017-12-01

We have developed an efficient computational framework for simulating multiple earthquake cycles with off-fault plasticity. The method is developed for the classical antiplane problem of a vertical strike-slip fault governed by rate-and-state friction, with inertial effects captured through the radiation-damping approximation. Both rate-independent plasticity and viscoplasticity are considered, where stresses are constrained by a Drucker-Prager yield condition. The off-fault volume is discretized using finite differences and tectonic loading is imposed by displacing the remote side boundaries at a constant rate. Time-stepping combines an adaptive Runge-Kutta method with an incremental solution process which makes use of an elastoplastic tangent stiffness tensor and the return-mapping algorithm. Solutions are verified by convergence tests and comparison to a finite element solution. We quantify how viscosity, isotropic hardening, and cohesion affect the magnitude and off-fault extent of plastic strain that develops over many ruptures. If hardening is included, plastic strain saturates after the first event and the response during subsequent ruptures is effectively elastic. For viscoplasticity without hardening, however, successive ruptures continue to generate additional plastic strain. In all cases, coseismic slip in the shallow sub-surface is diminished compared to slip accumulated at depth during interseismic loading. The evolution of this slip deficit with each subsequent event, however, is dictated by the plasticity model. Integration of the off-fault plastic strain from the viscoplastic model reveals that a significant amount of tectonic offset is accommodated by inelastic deformation ( ∼ 0.1 m per rupture, or ∼ 10% of the tectonic deformation budget).

Analysis of tablet compaction. I. Characterization of mechanical behavior of powder and powder/tooling friction.

PubMed

Cunningham, J C; Sinka, I C; Zavaliangos, A

2004-08-01

In this first of two articles on the modeling of tablet compaction, the experimental inputs related to the constitutive model of the powder and the powder/tooling friction are determined. The continuum-based analysis of tableting makes use of an elasto-plastic model, which incorporates the elements of yield, plastic flow potential, and hardening, to describe the mechanical behavior of microcrystalline cellulose over the range of densities experienced during tableting. Specifically, a modified Drucker-Prager/cap plasticity model, which includes material parameters such as cohesion, internal friction, and hydrostatic yield pressure that evolve with the internal state variable relative density, was applied. Linear elasticity is assumed with the elastic parameters, Young's modulus, and Poisson's ratio dependent on the relative density. The calibration techniques were developed based on a series of simple mechanical tests including diametrical compression, simple compression, and die compaction using an instrumented die. The friction behavior is measured using an instrumented die and the experimental data are analyzed using the method of differential slices. The constitutive model and frictional properties are essential experimental inputs to the finite element-based model described in the companion article. Copyright 2004 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 93:2022-2039, 2004

Gravity-driven postseismic deformation following the Mw 6.3 2009 L’Aquila (Italy) earthquake

PubMed Central

Albano, Matteo; Barba, Salvatore; Saroli, Michele; Moro, Marco; Malvarosa, Fabio; Costantini, Mario; Bignami, Christian; Stramondo, Salvatore

2015-01-01

The present work focuses on the postseismic deformation observed in the region of L’Aquila (central Italy) following the Mw 6.3 earthquake that occurred on April 6, 2009. A new, 16-month-long dataset of COSMO-SkyMed SAR images was analysed using the Persistent Scatterer Pairs interferometric technique. The analysis revealed the existence of postseismic ground subsidence in the mountainous rocky area of Mt Ocre ridge, contiguous to the sedimentary plain that experienced coseismic subsidence. The postseismic subsidence was characterized by displacements of 10 to 35 mm along the SAR line of sight. In the Mt Ocre ridge, widespread morphological elements associated with gravitational spreading have been previously mapped. We tested the hypothesis that the postseismic subsidence of the Mt Ocre ridge compensates the loss of equilibrium induced by the nearby coseismic subsidence. Therefore, we simulated the coseismic and postseismic displacement fields via the finite element method. We included the gravitational load and fault slip and accounted for the geometrical and rheological characteristics of the area. We found that the elastoplastic behaviour of the material under gravitational loading best explains the observed postseismic displacement. These findings emphasize the role of gravity in the postseismic processes at the fault scale. PMID:26553120

Simulation of interaction of damage agents of different shape with shaped-charge munition

NASA Astrophysics Data System (ADS)

Radchenko, P. A.; Batuev, S. P.; Radchenko, A. V.; Tukaev, A. M.

2017-01-01

The present paper studies the influence of the shape of projectile (damage agent) on its penetration capability. Steel projectiles of different shape have been considered as damage agents: sphere, regular tetrahedron, cube, cylinder and plate. The weight of projectiles has been kept the same. Antitank grenade has been used as a target. The study has been conducted by means of numerical simulation using finite element analysis. The simulation is three-dimensional. Behavior of materials has been described by elasto-plastic model taking into consideration the fracture and fragmentation of interacting bodies. The speed of interaction has been considered within the range of 800 to 2000 m/s. Research results demonstrated significant influence of the projectile shape on its penetration capability. Projectile in the shape of elongated cylinder has shown better penetration capability. Considering the weight of damage agents (except for sphere and plate) their maximum penetration capability has been reached at the speed of 1400 m/s. Increase of the speed of interaction has been followed by intensive fracture of damage agents and their penetration capability thus has worsened.

Progressive failure of lower San Fernando dam

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gu, W.H.; Morgenstern, N.R.; Robertson, P.K.

1993-02-01

Postearthquake deformation analyses of the lower San Fernando dam were conducted using an incremental finite-element method. In the analyses, an undrained elastoplastic model was used to simulate the collapse of liquefied materials. The model is developed based on the critical-state boundary-surface theory, the concept of steady-state strength, and the undrained behavior of liquefiable soils. A triggering condition in terms of a collapse surface was considered in this model. The hyperbolic strain-softening relationship has been introduced to simulate the postpeak behavior of liquefied materials. The analyses have shown that a progressive failure under undrained conditions may explain the observed response ofmore » the lower San Fernando dam following the 1971 earthquake. Stress redistribution initiated by the strain softening of liquefied materials is the main reason for undrained flow failures of dams, slopes, and foundations and can occur in a short period ranging from a few seconds to a few minutes. The liquefied zone after stress redistribution may be much larger than the initial liquefied zone caused directly by an earthquake. Therefore, a postearthquake deformation analysis may be essential in liquefaction stability evaluations.« less

Evaluation of TBM tunnels with respect to stability against spalling

NASA Astrophysics Data System (ADS)

Shaalan, Heyam; Ismail, Mohd Ashraf Mohd; Azit, Romziah

2017-10-01

As the depth of tunnels and underground construction increases, instability occurs in the form of rock bursting or spalling because of the induced stresses. Spalling may appear as a strong compressive stress causing crack growth behind the excavated surface and buckling of the thin rock slabs. In this paper, we describe how to reduce the rock spalling failure to increase the underground safety and the tunnel stability. Thus, a parametric study is implemented using 2-D Elasto-plastic finite elements stress analysis software to investigate the parameters that can minimize the extent and depth of the failure zone. The critical section of Pahang Selangor Raw Water Transfer Tunnel under high overburden is analyzed. The effect of the shotcrete lining thickness, tunnel size and the removal of fallouts or scaled v-notch on the failure zone depth is investigated. The results demonstrate that the shotcrete lining thickness has less influence on the failure depth, while a small tunnel diameter minimizes the failure depth. In addition, the stability of the tunnel improves by removing the loose rock mass.

Ultrasonic impact peening for the surface properties’ management

NASA Astrophysics Data System (ADS)

Mordyuk, Bohdan N.; Prokopenko, Georgiy I.

2007-12-01

It is demonstrated that the ultrasonic impact peening (UIP) technique is a beneficial method for essential increase in the fatigue durability of metallic materials due to the surface nanocrystallization and hardening process provided for severe plastic deformation of surface via multiple impacts of high velocity impact pins. Nano-scale grain structures were obtained in the surface layers of stainless steel, low carbon steel weld and different titanium alloys using developed equipment for the UIP. Both the surface nanostructure and compressive residual stresses are shown to attribute to the essential hardness increase. It is revealed experimentally using profilometry that new modification of the UIP apparatus providing high velocity "sliding" impacts leads to marked diminution of the surface roughness, which is another important factor affecting to the fatigue cracks initiation process. The two-dimensional finite element model is used to simulate the indent formation process during single impaction. The solid steel pin and the Al alloy plate are modeled as a rigid material and an elasto-plastic material, respectively. It is shown that the surface roughness magnitude depends on the correlation of the vertical and lateral load components.

Simulation for analysis and control of superplastic forming. Final report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zacharia, T.; Aramayo, G.A.; Simunovic, S.

1996-08-01

A joint study was conducted by Oak Ridge National Laboratory (ORNL) and the Pacific Northwest Laboratory (PNL) for the U.S. Department of Energy-Lightweight Materials (DOE-LWM) Program. the purpose of the study was to assess and benchmark the current modeling capabilities with respect to accuracy of predictions and simulation time. Two modeling capabilities with respect to accuracy of predictions and simulation time. Two simulation platforms were considered in this study, which included the LS-DYNA3D code installed on ORNL`s high- performance computers and the finite element code MARC used at PNL. both ORNL and PNL performed superplastic forming (SPF) analysis on amore » standard butter-tray geometry, which was defined by PNL, to better understand the capabilities of the respective models. The specific geometry was selected and formed at PNL, and the experimental results, such as forming time and thickness at specific locations, were provided for comparisons with numerical predictions. Furthermore, comparisons between the ORNL simulation results, using elasto-plastic analysis, and PNL`s results, using rigid-plastic flow analysis, were performed.« less

Frictional interactions in forming processes: New studies with transparent sapphire strip-drawing dies

NASA Astrophysics Data System (ADS)

Rao, R. S.; Lu, C. Y.; Wright, P. K.; Devenpeck, M. L.; Richmond, O.; Appleby, E. J.

1982-05-01

This research is concerned with the frictional interactions at the toolwork interfaces in the machining and strip-drawing processes. A novel feature is that transparent sapphire (single crystal Al2O3) is being used as the tool and die material. This allows the tribological features of the interface to be directly observed and recorded on movie-film. These qualitative studies provide information on the role of lubricants. In addition, techniques are being developed to quantify the velocity gradient along the interface. For example, in the drawing work it has been found that tracer markings (e.g. dye-spots), applied to the undrawn strip, remain intact during drawing and can be tracked along the sapphire/strip interface. Such data will be used as input to a finite-element, elasto-plastic-workhardening model of the deformation process. The latter can compute strip deformation characteristics, drawing forces and local coefficients of friction at the interface. Introductory results will be presented in this paper, obtained from drawing tin-plated mild steel with sapphire and cemented carbide dies. Drawing loads and die-separating forces will be presented and movie-films of the action of tracer markings at the interface shown. In order to demonstrate how this data can be used in an analysis of a large strain deformation process with friction, initial results from running the FIPDEF elasto-plastic code will be discussed. From a commercial viewpoint research on strip-drawing is of special interest to the can-making industry. From a physical viewpoint stripdrawing is of particular interest because it is a symmetrical, plane strain deformation and, in comparison with other metal processing operations, it is more readily modeled. However, until now the elasto-plastic codes that have been developed to predictively model drawing have had limitations: the most notable being that of quantifying the friction conditions at the die-work interface. Hence the specification of the latter is the ultimate goal of this research program.

The sound of friction: Real-time models, playability and musical applications

NASA Astrophysics Data System (ADS)

Serafin, Stefania

Friction, the tangential force between objects in contact, in most engineering applications needs to be removed as a source of noise and instabilities. In musical applications, friction is a desirable component, being the sound production mechanism of different musical instruments such as bowed strings, musical saws, rubbed bowls and any other sonority produced by interactions between rubbed dry surfaces. The goal of the dissertation is to simulate different instrument whose main excitation mechanism is friction. An efficient yet accurate model of a bowed string instrument, which combines the latest results in violin acoustics with the efficient digital waveguide approach, is provided. In particular, the bowed string physical model proposed uses a thermodynamic friction model in which the finite width of the bow is taken into account; this solution is compared to the recently developed elasto-plastic friction models used in haptics and robotics. Different solutions are also proposed to model the body of the instrument. Other less common instruments driven by friction are also proposed, and the elasto-plastic model is used to provide audio-visual simulations of everyday friction sounds such as squeaking doors and rubbed wine glasses. Finally, playability evaluations and musical applications in which the models have been used are discussed.

Vascular stents with submicrometer-scale surface patterning realized via titanium deep reactive ion etching

NASA Astrophysics Data System (ADS)

Gott, Shannon C.; Jabola, Benjamin A.; Rao, Masaru P.

2015-08-01

Herein, we report progress towards realization of vascular stents that will eventually provide opportunity for evaluating cellular response to rationally-designed, submicrometer-scale surface patterning in physiologically-relevant contexts, i.e. those that provide exposure to the complex multicellular milieu, flow-induced shear, and tissue-device interactions present in vivo. Specifically, using our novel titanium deep reactive ion etching technique (Ti DRIE), we discuss recent advances that have enabled: (a) fabrication of precisely-defined, grating-based surface patterns on planar Ti foils with minimum feature sizes as small as 0.15 μm (b) creation of cylindrical stents from micromachined planar Ti foils; and (c) integration of these processes to produce the first submicrometer-scale surface-patterned Ti stents that are compatible with conventional balloon catheter deployment techniques. We also discuss results from elastoplastic finite element simulations and preliminary mechanical testing of these devices to assess their mechanical performance. These efforts represent key steps towards our long-term goal of developing a new paradigm in stenting, where rationally-designed surface patterning provides a physical means for facilitating healing, and thus, improving outcomes in vascular intervention applications.

Optimal slot dimension for skirt support structure of coke drums

NASA Astrophysics Data System (ADS)

Wang, Edward; Xia, Zihui

2018-03-01

The skirt-to-shell junction weld on coke drums is susceptible to fatigue failure due to severe thermal cyclic stresses. One method to decrease junction stress is to add slots near the top of the skirt, thereby reducing the local stiffness close to the weld. The most common skirt slot design is thin relative to its circumferential spacing. A new slot design, which is significantly wider, is proposed. In this study, thermal-mechanical elastoplastic 3-D finite element models of coke drums are created to analyze the effect of different skirt designs on the stress/strain field near the shell-to-skirt junction weld, as well as any other critical stress locations in the overall skirt design. The results confirm that the inclusion of the conventional slot design effectively reduces stress in the junction weld. However, it has also been found that the critical stress location migrates from the shell-to-skirt junction weld to the slot ends. A method is used to estimate the fatigue life near the critical areas of each skirt slot design. It is found that wider skirt slots provide a significant improvement on fatigue life in the weld and slot area.

Numerical modelling of hydrologically-driven slope instability by means of porous media mechanics

NASA Astrophysics Data System (ADS)

Kakogiannou, Evanthia; Sanavia, Lorenzo; Lora, Marco; Schrefler, Bernhard

2015-04-01

Heavy rainfall can trigger slope failure which generally involves shallow soil deposit of different grading and origin usually in a state of partial saturation. In this case of slope instability, the behaviour of the soil slope is closely related not only to the distribution of pore-water pressure but also to the stress state during rainfall infiltration involving both mechanical and hydrological processes. In order to understand better these physical key processes, in this research work, the modelling of rainfall induced slope failure is considered as a coupled variably saturated hydro-mechanical problem. Therefore, the geometrically linear finite element code Comes-Geo for non-isothermal elasto-plastic multiphase solid porous materials is used, as developed by B.A. Schrefler and his co-workers. In this context, a detailed numerical analysis of an experimental slope stability test due to rainfall infiltration is presented. The main goals of this work are to understand the triggering mechanisms during the progressive failure, the effect of using different constitutive models of the mechanical soil behavior on the numerical results and the use of the second order work criterion on the detection of slope instability.

Subsidence and well failure in the South Belridge Diatomite field

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rouffignac, E.P. de; Bondor, P.L.; Karanikas, J.M. Hara, S.K.

1995-12-31

Withdrawal of fluids from shallow, thick and low strength rock can cause substantial reservoir compaction leading to surface subsidence and well failure. This is the case for the Diatomite reservoir, where over 10 ft of subsidence have occurred in some areas. Well failure rates have averaged over 3% per year, resulting in several million dollars per year in well replacement and repair costs in the South Belridge Diatomite alone. A program has been underway to address this issue, including experimental, modeling and field monitoring work. An updated elastoplastic rock law based on laboratory data has been generated which includes notmore » only standard shear failure mechanisms but also irreversible pore collapse occurring at low effective stresses (<150 psi). This law was incorporated into a commercial finite element geomechanics simulator. Since the late 1980s, a network of level survey monuments has been used to monitor subsidence at Belridge. Model predictions of subsidence in Section 33 compare very well with field measured data, which show that water injection reduces subsidence from 7--8 inches per year to 1--2 inches per year, but does not abate well failure.« less

Nonassociative plasticity model for cohesionless materials and its implementation in soil-structure interaction

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hashmi, Q.S.E.

A constitutive model based on rate-independent elastoplasticity concepts is developed and used to simulate the behavior of geologic materials under arbitrary three-dimensional stress paths. The model accounts for various factors such as friction, stress path, and stress history that influence the behavior of geologic materials. A hierarchical approach is adopted whereby models of progressively increasing sophistication are developed from a basic isotropic-hardening associate model. Nonassociativeness is introduced as correction or perturbation to the basic model. Deviation of normality of the plastic-strain increments to the yield surface F is captured through nonassociativeness. The plastic potential Q is obtained by applying amore » correction to F. This simplified approach restricts the number of extra parameters required to define the plastic potential Q. The material constants associated with the model are identified, and they are evaluated for three different sands (Leighton Buzzard, Munich and McCormick Ranch). The model is then verified by comparing predictions with laboratory tests from which the constants were found, and typical tests not used for finding the constants. Based on the above findings, a soil-footing system is analyzed using finite-element techniques.« less

Nucleation and growth of rolling contact failure of 440C bearing steel

NASA Technical Reports Server (NTRS)

Gupta, V.; Bastias, P. C.; Hahn, G. T.; Rubin, C. A.

1992-01-01

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

The method of lines in analyzing solids containing cracks

NASA Technical Reports Server (NTRS)

Gyekenyesi, John P.

1990-01-01

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Mechanical behavior, damage tolerance and durability of fiber metal laminates for aircraft structures

NASA Astrophysics Data System (ADS)

Wu, Guocai

This study systematically explores the mechanical behavior, damage tolerance and durability of fiber metal laminates, a promising candidate materials system for next generation aerospace structures. The experimental results indicated that GLARE laminates exhibited a bilinear deformation behavior under static in-plane loading. Both an analytical constitutive model based on a modified classical lamination theory which incorporates the elasto-plastic behavior of aluminum alloy and a numerical simulation based on finite element modeling are used to predict the nonlinear stress-strain response and deformation behavior of GLARE laminates. The blunt notched strength of GLARE laminates increased with decreasing specimen width and decreasing hole diameter. The notched strength of GLARE laminates was evaluated based on a modified point stress criterion. A computer simulation based on finite element method was performed to study stress concentration and distribution around the notch and verify the analytical and experimental results of notched strength. Good agreement is obtained between the model predictions and experimental results. Experimental results also indicate that GLARE laminates exhibited superior impact properties to those of monolithic 2024-T3 aluminum alloy at low velocity impact loading. The GLARE 5-2/1 laminate with 0°/90°/90°/0° fiber configuration exhibits a better impact resistance than the GLARE 4-3/2 laminate with 0°/90°/0° fiber orientation. The characteristic impact energies, the damage area, and the permanent deflection of laminates are used to evaluate the impact damage resistance. The post-impact residual tensile strength under various damage states ranging from the plastic dent, barely visible impact damage (BVID), clearly visible impact damage (CVID) up to the complete perforation was also measured and compared. The post-impact fatigue behavior under various stress levels and impact damage states was extensively explored. The damage initiation and progression, failure modes and crack propagation under different loading conditions were investigated and identified with microscopy, SEM, X-ray radiography, and by chemically removing outer aluminum layers.

Theoretical investigations into the influence of the position of a breaking line on the tensile failure of flat, round, bevel-edged tablets using finite element methodology (FEM) and its practical relevance for industrial tablet strength testing.

PubMed

Podczeck, Fridrun; Newton, J Michael; Fromme, Paul

2014-12-30

Flat, round tablets may have a breaking ("score") line. Pharmacopoeial tablet breaking load tests are diametral in their design, and industrially used breaking load testers often have automatic tablet feeding systems, which position the tablets between the loading platens of the machine with the breaking lines in random orientation to the applied load. The aim of this work was to ascertain the influence of the position of the breaking line in a diametral compression test using finite element methodology (FEM) and to compare the theoretical results with practical findings using commercially produced bevel-edged, scored tablets. Breaking line test positions at an angle of 0°, 22.5°, 45°, 67.5° and 90° relative to the loading plane were studied. FEM results obtained for fully elastic and elasto-plastic tablets were fairly similar, but they highlighted large differences in stress distributions depending on the position of the breaking line. The stress values at failure were predicted to be similar for tablets tested at an angle of 45° or above, whereas at lower test angles the predicted breaking loads were up to three times larger. The stress distributions suggested that not all breaking line angles would result in clean tensile failure. Practical results, however, did not confirm the differences in the predicted breaking loads, but they confirmed differences in the way tablets broke. The results suggest that it is not advisable to convert breaking loads obtained on scored tablets into tablet tensile strength values, and comparisons between different tablets or batches should carefully consider the orientation of the breaking line with respect to the loading plane, as the failure mechanisms appear to vary. Copyright © 2014 Elsevier B.V. All rights reserved.

Explicit and implicit springback simulation in sheet metal forming using fully coupled ductile damage and distortional hardening model

NASA Astrophysics Data System (ADS)

Yetna n'jock, M.; Houssem, B.; Labergere, C.; Saanouni, K.; Zhenming, Y.

2018-05-01

The springback is an important phenomenon which accompanies the forming of metallic sheets especially for high strength materials. A quantitative prediction of springback becomes very important for newly developed material with high mechanical characteristics. In this work, a numerical methodology is developed to quantify this undesirable phenomenon. This methodoly is based on the use of both explicit and implicit finite element solvers of Abaqus®. The most important ingredient of this methodology consists on the use of highly predictive mechanical model. A thermodynamically-consistent, non-associative and fully anisotropic elastoplastic constitutive model strongly coupled with isotropic ductile damage and accounting for distortional hardening is then used. An algorithm for local integration of the complete set of the constitutive equations is developed. This algorithm considers the rotated frame formulation (RFF) to ensure the incremental objectivity of the model in the framework of finite strains. This algorithm is implemented in both explicit (Abaqus/Explicit®) and implicit (Abaqus/Standard®) solvers of Abaqus® through the users routine VUMAT and UMAT respectively. The implicit solver of Abaqus® has been used to study spingback as it is generally a quasi-static unloading. In order to compare the methods `efficiency, the explicit method (Dynamic Relaxation Method) proposed by Rayleigh has been also used for springback prediction. The results obtained within U draw/bending benchmark are studied, discussed and compared with experimental results as reference. Finally, the purpose of this work is to evaluate the reliability of different methods predict efficiently springback in sheet metal forming.

Thermodynamically consistent constitutive equations for nonisothermal large strain, elasto-plastic, creep behavior

NASA Technical Reports Server (NTRS)

Riff, R.; Carlson, R. L.; Simitses, G. J.

1985-01-01

The paper is concerned with the development of constitutive relations for large nonisothermal elastic-viscoplastic deformations for metals. The kinematics of elastic-plastic deformation, valid for finite strains and rotations, is presented. The resulting elastic-plastic uncoupled equations for the deformation rate combined with use of the incremental elasticity law permits a precise and purely deductive development of elastic-viscoplastic theory. It is shown that a phenomenological thermodynamic theory in which the elastic deformation and the temperature are state variables, including few internal variables, can be utilized to construct elastic-viscoplastic constitutive equations, which are appropriate for metals. The limiting case of inviscid plasticity is examined.

Elastoplastic properties of a low-modulus titanium-based β alloy

NASA Astrophysics Data System (ADS)

Betekhtin, V. I.; Kolobov, Yu. R.; Golosova, O. A.; Kardashev, B. K.; Kadomtsev, A. G.; Narykova, M. V.; Ivanov, M. B.; Vershinina, T. N.

2013-10-01

The elastoplastic properties (elastic modulus, amplitude-independent damping ratio, microplastic flow stress) of a Ti-26Nb-7Mo-12Zr titanium β alloy are determined using an acoustic resonance method. The effect of the strain during thermomechanical treatment on the structural features of the micro-crystalline alloy and, hence, its elastoplastic properties is analyzed.

Application of boundary integral equations to elastoplastic problems

NASA Technical Reports Server (NTRS)

Mendelson, A.; Albers, L. U.

1975-01-01

The application of boundary integral equations to elastoplastic problems is reviewed. Details of the analysis as applied to torsion problems and to plane problems is discussed. Results are presented for the elastoplastic torsion of a square cross section bar and for the plane problem of notched beams. A comparison of different formulations as well as comparisons with experimental results are presented.

Orientation Effects in Fault Reactivation in Geological CO2 Sequestration

NASA Astrophysics Data System (ADS)

Castelletto, N.; Ferronato, M.; Gambolati, G.; Janna, C.; Teatini, P.

2012-12-01

Geological CO2 sequestration remains one of the most promising option for reducing the greenhouse gases emission. The accurate simulation of the complex coupled physical processes occurring during the injection and the post-injection stage represents a key issue for investigating the feasibility and the safety of the sequestration. The fluid-dynamical and geochemical aspects related to sequestering CO2 underground have been widely debated in the scientific literature over more than one decade. Recently, the importance of geomechanical processes has been widely recognized. In the present modeling study, we focus on fault reactivation induced by injection, an essential aspect for the evaluation of CO2 sequestration projects that needs to be adequately investigated to avoid the generation of preferential leaking path for CO2 and the related risk of induced seismicity. We use a geomechanical model based on the structural equations of poroelasticity solved by the Finite Element (FE) - Interface Element (IE) approach. Standard FEs are used to represent a continuum, while IEs prove especially suited to assess the relative displacements of adjacent elements such as the opening and slippage of existing faults or the generation of new fractures [1]. The IEs allow for the modeling of fault mechanics using an elasto-plastic constitutive law based on the Mohr-Coulomb failure criterion. We analyze the reactivation of a single fault in a synthetic reservoir by varying the fault orientation and size, hydraulic conductivity of the faulted zone, initial vertical and horizontal stress state and Mohr-Coulomb parameters (i.e., friction angle and cohesion). References: [1] Ferronato, M., G. Gambolati, C. Janna, and P. Teatini (2008), Numerical modeling of regional faults in land subsidence prediction above gas/oil reservoirs, Int. J. Numer. Anal. Methods Geomech., 32, 633-657.

Mathematical modeling of thermal stresses in basic oxygen furnace hood tubes

NASA Astrophysics Data System (ADS)

Samarasekera, I. V.

1985-06-01

The stress-strain history of Basic Oxygen Furnace hood tubes during thermal cycling has been computed using heat flow and stress analyses. The steady-state temperature distribution in a transverse section of the tube was computed at a location where gas temperature in the hood could be expected to be a maximum. Calculations were performed for peak gas temperatures in the range 1950 to 2480 °C (3500 to 4500 °F). The stress-strain history of an element of material located at the center of the tube hot face was traced for three consecutive cycles using elasto-plastic finite-element analysis. It has been shown that the state of stress in the element alternates between compression and tension as the tube successively heats and cools. Yielding and plastic flow occurs at the end of each half of a given cycle. It was postulated that owing to repctitive yielding, plastic strain energy accumulates causing failure of the tubes by fatigue in the low cycle region. Using fatigue theory a conservative estimate for tube life was arrived at. In-plant observations support this mechanism of failure, and the number of cycles within which tube cracking was observed compares reasonably with model predictions. Utilizing the heat flow and stress models it was recommended that tube life could be enhanced by changing the tube material to ARMCO 17-4 pH or AISI 405 steel or alternatively reconstructing hoods with AISI 316L tubes of reduced thickness. These recommendations were based on the criterion that low-cycle fatigue failure could be averted if the magnitude of the cyclic strain could be reduced or if macroscopic plastic flow could be prevented.

Low velocity impact analysis of composite laminated plates

NASA Astrophysics Data System (ADS)

Zheng, Daihua

2007-12-01

In the past few decades polymer composites have been utilized more in structures where high strength and light weight are major concerns, e.g., aircraft, high-speed boats and sports supplies. It is well known that they are susceptible to damage resulting from lateral impact by foreign objects, such as dropped tools, hail and debris thrown up from the runway. The impact response of the structures depends not only on the material properties but also on the dynamic behavior of the impacted structure. Although commercial software is capable of analyzing such impact processes, it often requires extensive expertise and rigorous training for design and analysis. Analytical models are useful as they allow parametric studies and provide a foundation for validating the numerical results from large-scale commercial software. Therefore, it is necessary to develop analytical or semi-analytical models to better understand the behaviors of composite structures under impact and their associated failure process. In this study, several analytical models are proposed in order to analyze the impact response of composite laminated plates. Based on Meyer's Power Law, a semi-analytical model is obtained for small mass impact response of infinite composite laminates by the method of asymptotic expansion. The original nonlinear second-order ordinary differential equation is transformed into two linear ordinary differential equations. This is achieved by neglecting high-order terms in the asymptotic expansion. As a result, the semi-analytical solution of the overall impact response can be applied to contact laws with varying coefficients. Then an analytical model accounting for permanent deformation based on an elasto-plastic contact law is proposed to obtain the closed-form solutions of the wave-controlled impact responses of composite laminates. The analytical model is also used to predict the threshold velocity for delamination onset by combining with an existing quasi-static delamination criterion. The predictions are compared with experimental data and explicit finite element LS-DYNA simulation. The comparisons show reasonable agreement. Furthermore, an analytical model is developed to evaluate the combined effects of prestresses and permanent deformation based on the linearized elasto-plastic contact law and the Laplace Transform technique. It is demonstrated that prestresses do not have noticeable effects on the time history of contact force and strains, but they have significant consequences on the plate central displacement. For a impacted composite laminate with the presence of prestresses, the contact force increases with the increasing of the mass of impactor, thickness and interlaminar shear strength of the laminate. The combined analytical and numerical investigations provide validated models for elastic and elasto-plastic impact analysis of composite structures and shed light on the design of impact-resistant composite systems.

Calculation analysis of magnetic-pulse compaction of explosively formed high-velocity metal elements used for meteoroid protection testing

NASA Astrophysics Data System (ADS)

Fedorov, Sergey V.; Selivanov, Victor V.; Veldanov, Vladislav A.

2017-06-01

Accumulation of microdamages as a result of intensive plastic deformation leads to a decrease in the average density of the high-velocity elements that are formed at the explosive collapse of the special shape metal liners. For compaction of such elements in tests of their spacecraft meteoroid protection reliability, the use of magnetic-field action on the produced elements during their movement trajectory before interaction with a target is proposed. On the basis of numerical modeling within the one-dimensional axisymmetric problem of continuum mechanics and electrodynamics, the physical processes occurring in the porous conducting elastoplastic cylinder placed in a magnetic field are investigated. Using this model, the parameters of the magnetic-pulse action necessary for the compaction of the steel and aluminum elements are determined.

Effective viscoelastic properties of shales.

NASA Astrophysics Data System (ADS)

Cornet, Jan; Dabrowski, Marcin; Schmid, Daniel

2017-04-01

Shales are often characterized as being elasto-plastic: they deform elastically for stresses below a certain yield and plastically at the limit. This approach dismisses any time dependent behavior that occurs in nature. Our goal is to better understand this time dependency by considering the visco-elastic behavior of shales before plasticity is reached. Shales are also typically heterogeneous and the question arises as to how to derive their effective properties in order to model them as a homogeneous medium. We model shales using inclusion based models due to their versatility and their ability to represent the microstructure. The inclusions represent competent quartz or calcite grains which are set in a viscous matrix made of clay minerals. Our approach relies on both numerical and analytical results in two dimension and we use them to cross check each other. The numerical results are obtained using MILAMIN, a fast-finite element solver for large problems, while the analytical solutions are based on the correspondence principle of linear viscoelasticity. This principle allows us to use the results on effective properties already derived for elastic bodies and to adapt them to viscoelastic bodies. We start by revisiting the problem of a single inclusion in an infinite medium and then move on to consider many inclusions.

Mechanical properties of metal-ceramic nanolaminates: Effect of constraint and temperature

DOE PAGES

Yang, Ling Wei; Mayer, Carl; Li, Nan; ...

2017-09-21

Al/SiC nanolaminates with equal nominal thicknesses of the Al and SiC layers (10, 25, 50 and 100 nm) were manufactured by magnetron sputtering. The mechanical properties were measured at 25 °C and 100 °C by means of nanoindentation and micropillar compression tests and the deformation mechanisms were analyzed by in situ micropillar compression tests in the transmission electron microscope. In addition, finite element simulations of both tests were carried out to ascertain the role played by the strength of the Al layers and by the elastic constraint of the ceramic layers on the plastic flow of Al in the mechanicalmore » response. It was found that the mechanical response was mainly controlled by the constraint during nanoindentation or micropillar compression tests of very thin layered (≈10 nm) laminates, while the influence of the strength of Al layers was not as critical. This behavior was reversed, however, for thick layered laminates (100 nm). Here, these mechanisms point to the different effects of layer thickness during nanoindentation and micropillar compression, at both temperatures, and showed the critical role played by constraint on the mechanical response of nanolaminates made of materials with a very large difference in the elasto-plastic properties.« less

Computational dynamics of soft machines

NASA Astrophysics Data System (ADS)

Hu, Haiyan; Tian, Qiang; Liu, Cheng

2017-06-01

Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and actuation. Hence, soft machines have raised great challenges to computational dynamics. In this review article, recent studies of the authors on the dynamic modeling, numerical simulation, and experimental validation of soft machines are summarized in the framework of multibody system dynamics. The dynamic modeling approaches are presented first for the geometric nonlinearities of coupled overall motions and large deformations of a soft component, the physical nonlinearities of a soft component made of hyperelastic or elastoplastic materials, and the frictional contacts/impacts of soft components, respectively. Then the computation approach is outlined for the dynamic simulation of soft machines governed by a set of differential-algebraic equations of very high dimensions, with an emphasis on the efficient computations of the nonlinear elastic force vector of finite elements. The validations of the proposed approaches are given via three case studies, including the locomotion of a soft quadrupedal robot, the spinning deployment of a solar sail of a spacecraft, and the deployment of a mesh reflector of a satellite antenna, as well as the corresponding experimental studies. Finally, some remarks are made for future studies.

Evaluation of transtension and transpression within contractional fault steps: Comparing kinematic and mechanical models to field data

NASA Astrophysics Data System (ADS)

Nevitt, Johanna M.; Pollard, David D.; Warren, Jessica M.

2014-03-01

Rock deformation often is investigated using kinematic and/or mechanical models. Here we provide a direct comparison of these modeling techniques in the context of a deformed dike within a meter-scale contractional fault step. The kinematic models consider two possible shear plane orientations and various modes of deformation (simple shear, transtension, transpression), while the mechanical model uses the finite element method and assumes elastoplastic constitutive behavior. The results for the kinematic and mechanical models are directly compared using the modeled maximum and minimum principal stretches. The kinematic analysis indicates that the contractional step may be classified as either transtensional or transpressional depending on the modeled shear plane orientation, suggesting that these terms may be inappropriate descriptors of step-related deformation. While the kinematic models do an acceptable job of depicting the change in dike shape and orientation, they are restricted to a prescribed homogeneous deformation. In contrast, the mechanical model allows for heterogeneous deformation within the step to accurately represent the deformation. The ability to characterize heterogeneous deformation and include fault slip - not as a prescription, but as a solution to the governing equations of motion - represents a significant advantage of the mechanical model over the kinematic models.

Phase-field modeling of the beta to omega phase transformation in Zr–Nb alloys

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yeddu, Hemantha Kumar; Lookman, Turab

A three-dimensional elastoplastic phase-field model is developed, using the Finite Element Method (FEM), for modeling the athermal beta to omega phase transformation in Zr–Nb alloys by including plastic deformation and strain hardening of the material. The microstructure evolution during athermal transformation as well as under different stress states, e.g. uni-axial tensile and compressive, bi-axial tensile and compressive, shear and tri-axial loadings, is studied. The effects of plasticity, stress states and the stress loading direction on the microstructure evolution as well as on the mechanical properties are studied. The input data corresponding to a Zr – 8 at.% Nb alloy aremore » acquired from experimental studies as well as by using the CALPHAD method. Our simulations show that the four different omega variants grow as ellipsoidal shaped particles. Our results show that due to stress relaxation, the athermal phase transformation occurs slightly more readily in the presence of plasticity compared to that in its absence. The evolution of omega phase is different under different stress states, which leads to the differences in the mechanical properties of the material. The variant selection mechanism, i.e. formation of different variants under different stress loading directions, is also nicely captured by our model.« less

Modeling the rubbing contact in honeycomb seals

NASA Astrophysics Data System (ADS)

Fischer, Tim; Welzenbach, Sarah; Meier, Felix; Werner, Ewald; kyzy, Sonun Ulan; Munz, Oliver

2018-03-01

Metallic honeycomb labyrinth seals are commonly used as sealing systems in gas turbine engines. Because of their capability to withstand high thermo-mechanical loads and oxidation, polycrystalline nickel-based superalloys, such as Hastelloy X and Haynes 214, are used as sealing material. In addition, these materials must exhibit a tolerance against rubbing between the rotating part and the stationary seal component. The tolerance of the sealing material against rubbing preserves the integrity of the rotating part. In this article, the rubbing behavior at the rotor-stator interface is considered numerically. A simulation model is incorporated into the commercial finite element code ABAQUS/explicit and is utilized to simulate a simplified rubbing process. A user-defined interaction routine between the contact surfaces accounts for the thermal and mechanical interfacial behavior. Furthermore, an elasto-plastic constitutive material law captures the extreme temperature conditions and the damage behavior of the alloys. To validate the model, representative quantities of the rubbing process are determined and compared with experimental data from the literature. The simulation results correctly reproduce the observations made on a test rig with a reference stainless steel material (AISI 304). A parametric study using the nickel-based superalloys reveals a clear dependency of the rubbing behavior on the sliding and incursion velocity. Compared to each other, the two superalloys studied exhibit a different rubbing behavior.

Strain-dependent Damage Evolution and Velocity Reduction in Fault Zones Induced by Earthquake Rupture

NASA Astrophysics Data System (ADS)

Zhong, J.; Duan, B.

2009-12-01

Low-velocity fault zones (LVFZs) with reduced seismic velocities relative to the surrounding wall rocks are widely observed around active faults. The presence of such a zone will affect rupture propagation, near-field ground motion, and off-fault damage in subsequent earth-quakes. In this study, we quantify the reduction of seismic velocities caused by dynamic rup-ture on a 2D planar fault surrounded by a low-velocity fault zone. First, we implement the damage rheology (Lyakhovsky et al. 1997) in EQdyna (Duan and Oglesby 2006), an explicit dynamic finite element code. We further extend this damage rheology model to include the dependence of strains on crack density. Then, we quantify off-fault continuum damage distribution and velocity reduction induced by earthquake rupture with the presence of a preexisting LVFZ. We find that the presence of a LVFZ affects the tempo-spatial distribu-tions of off-fault damage. Because lack of constraint in some damage parameters, we further investigate the relationship between velocity reduction and these damage prameters by a large suite of numerical simulations. Slip velocity, slip, and near-field ground motions computed from damage rheology are also compared with those from off-fault elastic or elastoplastic responses. We find that the reduction in elastic moduli during dynamic rupture has profound impact on these quantities.

Numerical analysis of nuclear power plant structure subjected to aircraft crash

NASA Astrophysics Data System (ADS)

Saberi, Reza; Alinejad, Majid; Mahdavi, Mir Omid; Sepanloo, Kamran

2017-12-01

An aircraft crashing into a nuclear containment may induce a series of disasters related to containment capacity, including local penetration and perforation of the containment, intensive vibrations, and fire ignited after jet fuel leakage. In this study, structural safety of a reinforced concrete containment vessel (RCCV) has been studied against the direct hit of Airbus A320, Boeing 707-320 and Phantom F4 aircrafts. ABAQUS/explicit finite element code has been used to carry out the three-dimensional numerical simulations. The impact locations identified on the nuclear containment structure are mid height of containment, center of the cylindrical portion, junction of dome and cylinder, and over the cylindrical portion close to the foundation level. The loading of the aircraft has been assigned through the corresponding reaction-time response curve. The concrete damaged plasticity model was predicted to simulate the behavior of concrete while the behavior of steel reinforcement was incorporated using elastoplastic material model. Dynamic loading conditions were considered using dynamic increase factor. The mid height of containment and center of cylindrical portion have been found to experience most severe deformation against each aircraft crash. It has also been found that compression damage in concrete is not critical at none of the impact locations.

Capturing strain localization behind a geosynthetic-reinforced soil wall

NASA Astrophysics Data System (ADS)

Lai, Timothy Y.; Borja, Ronaldo I.; Duvernay, Blaise G.; Meehan, Richard L.

2003-04-01

This paper presents the results of finite element (FE) analyses of shear strain localization that occurred in cohesionless soils supported by a geosynthetic-reinforced retaining wall. The innovative aspects of the analyses include capturing of the localized deformation and the accompanying collapse mechanism using a recently developed embedded strong discontinuity model. The case study analysed, reported in previous publications, consists of a 3.5-m tall, full-scale reinforced wall model deforming in plane strain and loaded by surcharge at the surface to failure. Results of the analysis suggest strain localization developing from the toe of the wall and propagating upward to the ground surface, forming a curved failure surface. This is in agreement with a well-documented failure mechanism experienced by the physical wall model showing internal failure surfaces developing behind the wall as a result of the surface loading. Important features of the analyses include mesh sensitivity studies and a comparison of the localization properties predicted by different pre-localization constitutive models, including a family of three-invariant elastoplastic constitutive models appropriate for frictional/dilatant materials. Results of the analysis demonstrate the potential of the enhanced FE method for capturing a collapse mechanism characterized by the presence of a failure, or slip, surface through earthen materials.

A Generalized Orthotropic Elasto-Plastic Material Model for Impact Analysis

NASA Astrophysics Data System (ADS)

Hoffarth, Canio

Composite materials are now beginning to provide uses hitherto reserved for metals in structural systems such as airframes and engine containment systems, wraps for repair and rehabilitation, and ballistic/blast mitigation systems. These structural systems are often subjected to impact loads and there is a pressing need for accurate prediction of deformation, damage and failure. There are numerous material models that have been developed to analyze the dynamic impact response of polymer matrix composites. However, there are key features that are missing in those models that prevent them from providing accurate predictive capabilities. In this dissertation, a general purpose orthotropic elasto-plastic computational constitutive material model has been developed to predict the response of composites subjected to high velocity impacts. The constitutive model is divided into three components - deformation model, damage model and failure model, with failure to be added at a later date. The deformation model generalizes the Tsai-Wu failure criteria and extends it using a strain-hardening-based orthotropic yield function with a non-associative flow rule. A strain equivalent formulation is utilized in the damage model that permits plastic and damage calculations to be uncoupled and capture the nonlinear unloading and local softening of the stress-strain response. A diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The overall framework is driven by experimental tabulated temperature and rate-dependent stress-strain data as well as data that characterizes the damage matrix and failure. The developed theory has been implemented in a commercial explicit finite element analysis code, LS-DYNARTM, as MAT213. Several verification and validation tests using a commonly available carbon-fiber composite, Toyobo's T800/F3900, have been carried and the results show that the theory and implementation are efficient, robust and accurate.

Evolution of interphase and intergranular strain in zirconium-niobium alloys during deformation at room temperature

NASA Astrophysics Data System (ADS)

Cai, Song

Zr-2.5Nb is currently used for pressure tubes in the CANDU (CANada Deuterium Uranium) reactor. A complete understanding of the deformation mechanism of Zr-2.5Nb is important if we are to accurately predict the in-reactor performance of pressure tubes and guarantee normal operation of the reactors. This thesis is a first step in gaining such an understanding; the deformation mechanism of ZrNb alloys at room temperature has been evaluated through studying the effect of texture and microstructure on deformation. In-situ neutron diffraction was used to monitor the evolution of the lattice strain of individual grain families along both the loading and Poisson's directions and to track the development of interphase and intergranular strains during deformation. The following experiments were carried out with data interpreted using elasto-plastic modeling techniques: (1) Compression tests of a 100%betaZr material at room temperature. (2) Tension and compression tests of hot rolled Zr-2.5Nb plate material. (3) Compression of annealed Zr-2.5Nb. (4) Cyclic loading of the hot rolled Zr-2.5Nb. (5) Compression tests of ZrNb alloys with different Nb and oxygen contents. The experimental results were interpreted using a combination of finite element (FE) and elasto-plastic self-consistent (EPSC) models. The phase properties and phase interactions well represented by the FE model, the EPSC model successfully captured the evolution of intergranular constraint during deformation and provided reasonable estimates of the critical resolved shear stress and hardening parameters of different slip systems under different conditions. The consistency of the material parameters obtained by the EPSC model allows the deformation mechanism at room temperature and the effect of textures and microstructures of ZrNb alloys to be understood. This work provides useful information towards manufacturing of Zr-2.5Nb components and helps in producing ideal microstructures and material properties for pressure tubes. Also it is helpful in guiding the development of new materials for the next generation of nuclear reactors. Furthermore, the large data set obtained from this study can be used in evaluation and improving current and future polycrystalline deformation models.

The role of elasticity in simulating long-term tectonic extension

NASA Astrophysics Data System (ADS)

Olive, Jean-Arthur; Behn, Mark D.; Mittelstaedt, Eric; Ito, Garrett; Klein, Benjamin Z.

2016-05-01

While elasticity is a defining characteristic of the Earth's lithosphere, it is often ignored in numerical models of long-term tectonic processes in favour of a simpler viscoplastic description. Here we assess the consequences of this assumption on a well-studied geodynamic problem: the growth of normal faults at an extensional plate boundary. We conduct 2-D numerical simulations of extension in elastoplastic and viscoplastic layers using a finite difference, particle-in-cell numerical approach. Our models simulate a range of faulted layer thicknesses and extension rates, allowing us to quantify the role of elasticity on three key observables: fault-induced topography, fault rotation, and fault life span. In agreement with earlier studies, simulations carried out in elastoplastic layers produce rate-independent lithospheric flexure accompanied by rapid fault rotation and an inverse relationship between fault life span and faulted layer thickness. By contrast, models carried out with a viscoplastic lithosphere produce results that may qualitatively resemble the elastoplastic case, but depend strongly on the product of extension rate and layer viscosity U × ηL. When this product is high, fault growth initially generates little deformation of the footwall and hanging wall blocks, resulting in unrealistic, rigid block-offset in topography across the fault. This configuration progressively transitions into a regime where topographic decay associated with flexure is fully accommodated within the numerical domain. In addition, high U × ηL favours the sequential growth of multiple short-offset faults as opposed to a large-offset detachment. We interpret these results by comparing them to an analytical model for the fault-induced flexure of a thin viscous plate. The key to understanding the viscoplastic model results lies in the rate-dependence of the flexural wavelength of a viscous plate, and the strain rate dependence of the force increase associated with footwall and hanging wall bending. This behaviour produces unrealistic deformation patterns that can hinder the geological relevance of long-term rifting models that assume a viscoplastic rheology.

Microscopic processes controlling the Herschel-Bulkley exponent

NASA Astrophysics Data System (ADS)

Lin, Jie; Wyart, Matthieu

2018-01-01

The flow curve of various yield stress materials is singular as the strain rate vanishes and can be characterized by the so-called Herschel-Bulkley exponent n =1 /β . A mean-field approximation due to Hebraud and Lequeux (HL) assumes mechanical noise to be Gaussian and leads to β =2 in rather good agreement with observations. Here we prove that the improved mean-field model where the mechanical noise has fat tails instead leads to β =1 with logarithmic correction. This result supports that HL is not a suitable explanation for the value of β , which is instead significantly affected by finite-dimensional effects. From considerations on elastoplastic models and on the limitation of speed at which avalanches of plasticity can propagate, we argue that β =1 +1 /(d -df) , where df is the fractal dimension of avalanches and d the spatial dimension. Measurements of df then supports that β ≈2.1 and β ≈1.7 in two and three dimensions, respectively. We discuss theoretical arguments leading to approximations of β in finite dimensions.

The Relation of Finite Element and Finite Difference Methods

NASA Technical Reports Server (NTRS)

Vinokur, M.

1976-01-01

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Computational mechanics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Raboin, P J

1998-01-01

The Computational Mechanics thrust area is a vital and growing facet of the Mechanical Engineering Department at Lawrence Livermore National Laboratory (LLNL). This work supports the development of computational analysis tools in the areas of structural mechanics and heat transfer. Over 75 analysts depend on thrust area-supported software running on a variety of computing platforms to meet the demands of LLNL programs. Interactions with the Department of Defense (DOD) High Performance Computing and Modernization Program and the Defense Special Weapons Agency are of special importance as they support our ParaDyn project in its development of new parallel capabilities for DYNA3D.more » Working with DOD customers has been invaluable to driving this technology in directions mutually beneficial to the Department of Energy. Other projects associated with the Computational Mechanics thrust area include work with the Partnership for a New Generation Vehicle (PNGV) for ''Springback Predictability'' and with the Federal Aviation Administration (FAA) for the ''Development of Methodologies for Evaluating Containment and Mitigation of Uncontained Engine Debris.'' In this report for FY-97, there are five articles detailing three code development activities and two projects that synthesized new code capabilities with new analytic research in damage/failure and biomechanics. The article this year are: (1) Energy- and Momentum-Conserving Rigid-Body Contact for NIKE3D and DYNA3D; (2) Computational Modeling of Prosthetics: A New Approach to Implant Design; (3) Characterization of Laser-Induced Mechanical Failure Damage of Optical Components; (4) Parallel Algorithm Research for Solid Mechanics Applications Using Finite Element Analysis; and (5) An Accurate One-Step Elasto-Plasticity Algorithm for Shell Elements in DYNA3D.« less

Real time bolt preload monitoring using piezoceramic transducers and time reversal technique—a numerical study with experimental verification

NASA Astrophysics Data System (ADS)

Parvasi, Seyed Mohammad; Ho, Siu Chun Michael; Kong, Qingzhao; Mousavi, Reza; Song, Gangbing

2016-08-01

Bolted joints are ubiquitous structural elements, and form critical connections in mechanical and civil structures. As such, loosened bolted joints may lead to catastrophic failures of these structures, thus inspiring a growing interest in monitoring of bolted joints. A novel energy based wave method is proposed in this study to monitor the axial load of bolted joint connections. In this method, the time reversal technique was used to focus the energy of a piezoelectric (PZT)-generated ultrasound wave from one side of the interface to be measured as a signal peak by another PZT transducer on the other side of the interface. A tightness index (TI) was defined and used to correlate the peak amplitude to the bolt axial load. The TI bypasses the need for more complex signal processing required in other energy-based methods. A coupled, electro-mechanical analysis with elasto-plastic finite element method was used to simulate and analyze the PZT based ultrasonic wave propagation through the interface of two steel plates connected by a single nut and bolt connection. Numerical results, backed by experimental results from testing on a bolted connection between two steel plates, revealed that the peak amplitude of the focused signal increases as the bolt preload (torque level) increases due to the enlarging true contact area of the steel plates. The amplitude of the focused peak saturates and the TI reaches unity as the bolt axial load reaches a threshold value. These conditions are associated with the maximum possible true contact area between the surfaces of the bolted connection.

ANSYS duplicate finite-element checker routine

NASA Technical Reports Server (NTRS)

Ortega, R.

1995-01-01

An ANSYS finite-element code routine to check for duplicated elements within the volume of a three-dimensional (3D) finite-element mesh was developed. The routine developed is used for checking floating elements within a mesh, identically duplicated elements, and intersecting elements with a common face. A space shuttle main engine alternate turbopump development high pressure oxidizer turbopump finite-element model check using the developed subroutine is discussed. Finally, recommendations are provided for duplicate element checking of 3D finite-element models.

The Natural Neighbour Radial Point Interpolation Meshless Method Applied to the Non-Linear Analysis

NASA Astrophysics Data System (ADS)

Dinis, L. M. J. S.; Jorge, R. M. Natal; Belinha, J.

2011-05-01

In this work the Natural Neighbour Radial Point Interpolation Method (NNRPIM), is extended to large deformation analysis of elastic and elasto-plastic structures. The NNPRIM uses the Natural Neighbour concept in order to enforce the nodal connectivity and to create a node-depending background mesh, used in the numerical integration of the NNRPIM interpolation functions. Unlike the FEM, where geometrical restrictions on elements are imposed for the convergence of the method, in the NNRPIM there are no such restrictions, which permits a random node distribution for the discretized problem. The NNRPIM interpolation functions, used in the Galerkin weak form, are constructed using the Radial Point Interpolators, with some differences that modify the method performance. In the construction of the NNRPIM interpolation functions no polynomial base is required and the used Radial Basis Function (RBF) is the Multiquadric RBF. The NNRPIM interpolation functions posses the delta Kronecker property, which simplify the imposition of the natural and essential boundary conditions. One of the scopes of this work is to present the validation the NNRPIM in the large-deformation elasto-plastic analysis, thus the used non-linear solution algorithm is the Newton-Rapson initial stiffness method and the efficient "forward-Euler" procedure is used in order to return the stress state to the yield surface. Several non-linear examples, exhibiting elastic and elasto-plastic material properties, are studied to demonstrate the effectiveness of the method. The numerical results indicated that NNRPIM handles large material distortion effectively and provides an accurate solution under large deformation.

Numerical computation of transonic flows by finite-element and finite-difference methods

NASA Technical Reports Server (NTRS)

Hafez, M. M.; Wellford, L. C.; Merkle, C. L.; Murman, E. M.

1978-01-01

Studies on applications of the finite element approach to transonic flow calculations are reported. Different discretization techniques of the differential equations and boundary conditions are compared. Finite element analogs of Murman's mixed type finite difference operators for small disturbance formulations were constructed and the time dependent approach (using finite differences in time and finite elements in space) was examined.

Stress–strain state in an elastoplastic pipe taking into account the temperature and compressibility of the material

NASA Astrophysics Data System (ADS)

Gornostaev, K. K.; Kovalev, A. V.; Malygina, Y. V.

2018-03-01

In the article the authors have considered the problem of determining the stress-strain state of the elastoplastic pipe with the Mises’ condition in case of plane strain for the compressible material taking into account the temperature. The task was solved using the method of the small parameter. The expressions for the fields of stresses and displacements were received as well as the ratio of the radius of the elastoplastic boundary in the zero and first approximations.

Finite Element Modeling of In-Situ Stresses near Salt Bodies

NASA Astrophysics Data System (ADS)

Sanz, P.; Gray, G.; Albertz, M.

2011-12-01

The in-situ stress field is modified around salt bodies because salt rock has no ability to sustain shear stresses. A reliable prediction of stresses near salt is important for planning safe and economic drilling programs. A better understanding of in-situ stresses before drilling can be achieved using finite element models that account for the creeping salt behavior and the elastoplastic response of the surrounding sediments. Two different geomechanical modeling techniques can be distinguished: "dynamic" modeling and "static" modeling. "Dynamic" models, also known as forward models, simulate the development of structural processes in geologic time. This technique provides the evolution of stresses and so it is used to simulate the initiation and development of structural features, such as, faults, folds, fractures, and salt diapers. The original or initial configuration and the unknown final configuration of forward models are usually significantly different therefore geometric non-linearities need to be considered. These models may be difficult to constrain when different tectonic, deposition, and erosion events, and the timing among them, needs to be accounted for. While dynamic models provide insight into the stress evolution, in many cases is very challenging, if not impossible, to forward model a configuration to its known present-day geometry; particularly in the case of salt layers that evolve into highly irregular and complex geometries. Alternatively, "static" models use the present-day geometry and present-day far-field stresses to estimate the present-day in-situ stress field inside a domain. In this case, it is appropriate to use a small deformation approach because initial and final configurations should be very similar, and more important, because the equilibrium of stresses should be stated in the present-day initial configuration. The initial stresses and the applied boundary conditions are constrained by the geologic setting and available data. This modeling technique does not predict the evolution of structural elements or stresses with time; therefore it does not provide any insight into the formation of fractures that were previously developed under a different stress condition or the development of overpressure generated by a high sedimentation rate. This work provides a validation for predicting in-situ stresses near salt using "static" models. We compare synthetic examples using both modeling techniques and show that stresses near salt predicted with "static" models are comparable to the ones generated by "dynamic" models.

Studies of finite element analysis of composite material structures

NASA Technical Reports Server (NTRS)

Douglas, D. O.; Holzmacher, D. E.; Lane, Z. C.; Thornton, E. A.

1975-01-01

Research in the area of finite element analysis is summarized. Topics discussed include finite element analysis of a picture frame shear test, BANSAP (a bandwidth reduction program for SAP IV), FEMESH (a finite element mesh generation program based on isoparametric zones), and finite element analysis of a composite bolted joint specimens.

An 8-node tetrahedral finite element suitable for explicit transient dynamic simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Key, S.W.; Heinstein, M.W.; Stone, C.M.

1997-12-31

Considerable effort has been expended in perfecting the algorithmic properties of 8-node hexahedral finite elements. Today the element is well understood and performs exceptionally well when used in modeling three-dimensional explicit transient dynamic events. However, the automatic generation of all-hexahedral meshes remains an elusive achievement. The alternative of automatic generation for all-tetrahedral finite element is a notoriously poor performer, and the 10-node quadratic tetrahedral finite element while a better performer numerically is computationally expensive. To use the all-tetrahedral mesh generation extant today, the authors have explored the creation of a quality 8-node tetrahedral finite element (a four-node tetrahedral finite elementmore » enriched with four midface nodal points). The derivation of the element`s gradient operator, studies in obtaining a suitable mass lumping and the element`s performance in applications are presented. In particular, they examine the 80node tetrahedral finite element`s behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element only samples constant strain states and, therefore, has 12 hourglass modes. In this regard, it bears similarities to the 8-node, mean-quadrature hexahedral finite element. Given automatic all-tetrahedral meshing, the 8-node, constant-strain tetrahedral finite element is a suitable replacement for the 8-node hexahedral finite element and handbuilt meshes.« less

Patient-specific finite element modeling of bones.

PubMed

Poelert, Sander; Valstar, Edward; Weinans, Harrie; Zadpoor, Amir A

2013-04-01

Finite element modeling is an engineering tool for structural analysis that has been used for many years to assess the relationship between load transfer and bone morphology and to optimize the design and fixation of orthopedic implants. Due to recent developments in finite element model generation, for example, improved computed tomography imaging quality, improved segmentation algorithms, and faster computers, the accuracy of finite element modeling has increased vastly and finite element models simulating the anatomy and properties of an individual patient can be constructed. Such so-called patient-specific finite element models are potentially valuable tools for orthopedic surgeons in fracture risk assessment or pre- and intraoperative planning of implant placement. The aim of this article is to provide a critical overview of current themes in patient-specific finite element modeling of bones. In addition, the state-of-the-art in patient-specific modeling of bones is compared with the requirements for a clinically applicable patient-specific finite element method, and judgment is passed on the feasibility of application of patient-specific finite element modeling as a part of clinical orthopedic routine. It is concluded that further development in certain aspects of patient-specific finite element modeling are needed before finite element modeling can be used as a routine clinical tool.

Analytical modeling of structure-soil systems for lunar bases

NASA Technical Reports Server (NTRS)

Macari-Pasqualino, Jose Emir

1989-01-01

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Thermo-magneto-elastoplastic coupling model of metal magnetic memory testing method for ferromagnetic materials

NASA Astrophysics Data System (ADS)

Shi, Pengpeng; Zhang, Pengcheng; Jin, Ke; Chen, Zhenmao; Zheng, Xiaojing

2018-04-01

Metal magnetic memory (MMM) testing (also known as micro-magnetic testing) is a new non-destructive electromagnetic testing method that can diagnose ferromagnetic materials at an early stage by measuring the MMM signal directly on the material surface. Previous experiments have shown that many factors affect MMM signals, in particular, the temperature, the elastoplastic state, and the complex environmental magnetic field. However, the fact that there have been only a few studies of either how these factors affect the signals or the physical coupling mechanisms among them seriously limits the industrial applications of MMM testing. In this paper, a nonlinear constitutive relation for a ferromagnetic material considering the influences of temperature and elastoplastic state is established under a weak magnetic field and is used to establish a nonlinear thermo-magneto-elastoplastic coupling model of MMM testing. Comparing with experimental data verifies that the proposed theoretical model can accurately describe the thermo-magneto-elastoplastic coupling influence on MMM signals. The proposed theoretical model can predict the MMM signals in a complex environment and so is expected to provide a theoretical basis for improving the degree of quantification in MMM testing.

Improved finite element methodology for integrated thermal structural analysis

NASA Technical Reports Server (NTRS)

Dechaumphai, P.; Thornton, E. A.

1982-01-01

An integrated thermal-structural finite element approach for efficient coupling of thermal and structural analysis is presented. New thermal finite elements which yield exact nodal and element temperatures for one dimensional linear steady state heat transfer problems are developed. A nodeless variable formulation is used to establish improved thermal finite elements for one dimensional nonlinear transient and two dimensional linear transient heat transfer problems. The thermal finite elements provide detailed temperature distributions without using additional element nodes and permit a common discretization with lower order congruent structural finite elements. The accuracy of the integrated approach is evaluated by comparisons with analytical solutions and conventional finite element thermal structural analyses for a number of academic and more realistic problems. Results indicate that the approach provides a significant improvement in the accuracy and efficiency of thermal stress analysis for structures with complex temperature distributions.

A general algorithm using finite element method for aerodynamic configurations at low speeds

NASA Technical Reports Server (NTRS)

Balasubramanian, R.

1975-01-01

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Effects of implant drilling parameters for pilot and twist drills on temperature rise in bone analog and alveolar bones.

PubMed

Chen, Yung-Chuan; Hsiao, Chih-Kun; Ciou, Ji-Sih; Tsai, Yi-Jung; Tu, Yuan-Kun

2016-11-01

This study concerns the effects of different drilling parameters of pilot drills and twist drills on the temperature rise of alveolar bones during dental implant procedures. The drilling parameters studied here include the feed rate and rotation speed of the drill. The bone temperature distribution was analyzed through experiments and numerical simulations of the drilling process. In this study, a three dimensional (3D) elasto-plastic dynamic finite element model (DFEM) was proposed to investigate the effects of drilling parameters on the bone temperature rise. In addition, the FE model is validated with drilling experiments on artificial human bones and porcine alveolar bones. The results indicate that 3D DFEM can effectively simulate the bone temperature rise during the drilling process. During the drilling process with pilot drills or twist drills, the maximum bone temperature occurred in the region of the cancellous bones close to the cortical bones. The feed rate was one of the important factors affecting the time when the maximum bone temperature occurred. Our results also demonstrate that the elevation of bone temperature was reduced as the feed rate increased and the drill speed decreased, which also effectively reduced the risk region of osteonecrosis. These findings can serve as a reference for dentists in choosing drilling parameters for dental implant surgeries. Copyright © 2016 IPEM. Published by Elsevier Ltd. All rights reserved.

Three-dimensional numerical simulations of crustal-scale wrenching using a non-linear failure criterion

NASA Astrophysics Data System (ADS)

Braun, Jean

1994-08-01

We have developed a three-dimensional finite element model to study wrench deformation of the crust regarded as an elasto-plastic material obeying Murrell's extension of Griffith's failure criterion. Numerical experiments using this model predict that the imposed basal wrenching is accommodated by an array of oblique Riedel-like shears and Y-shears (parallel to the direction of wrenching). The partitioning of deformation between the two types of structure depends on the width of the zone of imposed basal wrenching and the existence of a component of deformation in the x-direction (normal to the direction of wrenching). The Riedel shears are arranged in spiral-like structures that root into the basal wrench zone. In cross-section, the Riedel shears resemble wedge-shaped flower structures similar to those often observed in seismic cross-sections. The 'polarity' of the flower structures is positive (or palm-tree-like) in transpression experiments and negative (or tulip-like) in transtension experiments. The orientation of the Riedel shears throughout the crust obeys Mohr's hypothesis for incipient faulting combined with Murrell's failure criterion. The model also predicts plastic dilatancy inversely proportional to the square root of the confining pressure; this result agrees qualitatively with field observations and the results of sand-box experiments and quantitatively with direct measurement of dilatancy during high-pressure rock-deformation experiments.

Modeling interface shear behavior of granular materials using micro-polar continuum approach

NASA Astrophysics Data System (ADS)

Ebrahimian, Babak; Noorzad, Ali; Alsaleh, Mustafa I.

2018-01-01

Recently, the authors have focused on the shear behavior of interface between granular soil body and very rough surface of moving bounding structure. For this purpose, they have used finite element method and a micro-polar elasto-plastic continuum model. They have shown that the boundary conditions assumed along the interface have strong influences on the soil behavior. While in the previous studies, only very rough bounding interfaces have been taken into account, the present investigation focuses on the rough, medium rough and relatively smooth interfaces. In this regard, plane monotonic shearing of an infinite extended narrow granular soil layer is simulated under constant vertical pressure and free dilatancy. The soil layer is located between two parallel rigid boundaries of different surface roughness values. Particular attention is paid to the effect of surface roughness of top and bottom boundaries on the shear behavior of granular soil layer. It is shown that the interaction between roughness of bounding structure surface and the rotation resistance of bounding grains can be modeled in a reasonable manner through considered Cosserat boundary conditions. The influence of surface roughness is investigated on the soil shear strength mobilized along the interface as well as on the location and evolution of shear localization formed within the layer. The obtained numerical results have been qualitatively compared with experimental observations as well as DEM simulations, and acceptable agreement is shown.

Improving a complex finite-difference ground water flow model through the use of an analytic element screening model

USGS Publications Warehouse

Hunt, R.J.; Anderson, M.P.; Kelson, V.A.

1998-01-01

This paper demonstrates that analytic element models have potential as powerful screening tools that can facilitate or improve calibration of more complicated finite-difference and finite-element models. We demonstrate how a two-dimensional analytic element model was used to identify errors in a complex three-dimensional finite-difference model caused by incorrect specification of boundary conditions. An improved finite-difference model was developed using boundary conditions developed from a far-field analytic element model. Calibration of a revised finite-difference model was achieved using fewer zones of hydraulic conductivity and lake bed conductance than the original finite-difference model. Calibration statistics were also improved in that simulated base-flows were much closer to measured values. The improved calibration is due mainly to improved specification of the boundary conditions made possible by first solving the far-field problem with an analytic element model.This paper demonstrates that analytic element models have potential as powerful screening tools that can facilitate or improve calibration of more complicated finite-difference and finite-element models. We demonstrate how a two-dimensional analytic element model was used to identify errors in a complex three-dimensional finite-difference model caused by incorrect specification of boundary conditions. An improved finite-difference model was developed using boundary conditions developed from a far-field analytic element model. Calibration of a revised finite-difference model was achieved using fewer zones of hydraulic conductivity and lake bed conductance than the original finite-difference model. Calibration statistics were also improved in that simulated base-flows were much closer to measured values. The improved calibration is due mainly to improved specification of the boundary conditions made possible by first solving the far-field problem with an analytic element model.

Elasto-plastic properties of Cu-Nb nanolaminate

NASA Astrophysics Data System (ADS)

Betekhtin, V. I.; Kolobov, Yu. R.; Kardashev, B. K.; Golosov, E. V.; Narykova, M. V.; Kadomtsev, A. G.; Klimenko, D. N.; Karpov, M. I.

2012-02-01

The Young's modulus, internal friction, and microplastic flow stress in Cu-Nb nanolaminate has been determined by an acoustic technique. The influence of high hydrostatic compression (1 GPa) on these elasto-plastic properties of the nanolaminate has been studied.

Establishing the 3-D finite element solid model of femurs in partial by volume rendering.

PubMed

Zhang, Yinwang; Zhong, Wuxue; Zhu, Haibo; Chen, Yun; Xu, Lingjun; Zhu, Jianmin

2013-01-01

It remains rare to report three-dimensional (3-D) finite element solid model of femurs in partial by volume rendering method, though several methods of femoral 3-D finite element modeling are already available. We aim to analyze the advantages of the modeling method by establishing the 3-D finite element solid model of femurs in partial by volume rendering. A 3-D finite element model of the normal human femurs, made up of three anatomic structures: cortical bone, cancellous bone and pulp cavity, was constructed followed by pretreatment of the CT original image. Moreover, the finite-element analysis was carried on different material properties, three types of materials given for cortical bone, six assigned for cancellous bone, and single for pulp cavity. The established 3-D finite element of femurs contains three anatomical structures: cortical bone, cancellous bone, and pulp cavity. The compressive stress primarily concentrated in the medial surfaces of femur, especially in the calcar femorale. Compared with whole modeling by volume rendering method, the 3-D finite element solid model created in partial is more real and fit for finite element analysis. Copyright © 2013 Surgical Associates Ltd. Published by Elsevier Ltd. All rights reserved.

A Novel Approach to Model Earth Fissure Caused by Extensive Aquifer Exploitation and its Application to the Wuxi Case, China

NASA Astrophysics Data System (ADS)

Ye, Shujun; Franceschini, Andrea; Zhang, Yan; Janna, Carlo; Gong, Xulong; Yu, Jun; Teatini, Pietro

2018-03-01

Initially observed in the semiarid basins of southwestern USA, earth fissures due to aquifer over-exploitation are presently threatening a large number of subsiding basins in various countries worldwide. Different mechanics have been proposed to explain this process, such as differential compaction, horizontal movements, and fault reactivation. Numerical modeling and prediction of this major geohazard caused by overuse of groundwater resources are challenging because of two main requirements: shifting from the classical continuous to discontinuous geomechanics and incorporating two-dimensional features (the earth fissures) into large three-dimensional (3-D) modeling domain (the subsiding basin). In this work, we proposed a novel modeling approach to simulate earth fissure generation and propagation in 3-D complex geological settings. A nested two-scale approach associated with an original nonlinear elastoplastic finite element/interface element simulator allows modeling the mechanics of earth discontinuities, in terms of both sliding and opening. The model is applied on a case study in Wuxi, China, where groundwater pumping between 1985 and 2004 has caused land subsidence larger than 2 m. The model outcomes highlight that the presence of a shallow (˜80 m deep) bedrock ridge crossing the Yangtze River delta is the key factor triggering the earth fissure development in this area. Bending of the alluvial deposits around the ridge tip and shear stress due to the uneven piezometric change and asymmetrical shape of the bedrock have caused the earth fissure to onset at the land surface and propagate downward to a maximum depth of about 20-30 m. Maximum sliding and opening are computed in the range of 10-40 cm, in agreement with the order of magnitude estimated in the field.

Modeling Progressive Failure of Bonded Joints Using a Single Joint Finite Element

NASA Technical Reports Server (NTRS)

Stapleton, Scott E.; Waas, Anthony M.; Bednarcyk, Brett A.

2010-01-01

Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.

The Blended Finite Element Method for Multi-fluid Plasma Modeling

DTIC Science & Technology

2016-07-01

Briefing Charts 3. DATES COVERED (From - To) 07 June 2016 - 01 July 2016 4. TITLE AND SUBTITLE The Blended Finite Element Method for Multi-fluid Plasma...BLENDED FINITE ELEMENT METHOD FOR MULTI-FLUID PLASMA MODELING Éder M. Sousa1, Uri Shumlak2 1ERC INC., IN-SPACE PROPULSION BRANCH (RQRS) AIR FORCE RESEARCH...MULTI-FLUID PLASMA MODEL 2 BLENDED FINITE ELEMENT METHOD Blended Finite Element Method Nodal Continuous Galerkin Modal Discontinuous Galerkin Model

Spatial Convergence of Three Dimensional Turbulent Flows

NASA Technical Reports Server (NTRS)

Park, Michael A.; Anderson, W. Kyle

2016-01-01

Finite-volume and finite-element schemes, both implemented within the FUN3D flow solver, are evaluated for several test cases described on the Turbulence-Modeling Resource (TMR) web site. The cases include subsonic flow over a hemisphere cylinder, subsonic flow over a swept bump configuration, and supersonic flow in a square duct. The finite- volume and finite-element schemes are both used to obtain solutions for the first two cases, whereas only the finite-volume scheme is used for the supersonic duct. For the hemisphere cylinder, finite-element solutions obtained on tetrahedral meshes are compared with finite- volume solutions on mixed-element meshes. For the swept bump, finite-volume solutions have been obtained for both hexahedral and tetrahedral meshes and are compared with finite-element solutions obtained on tetrahedral meshes. For the hemisphere cylinder and the swept bump, solutions are obtained on a series of meshes with varying grid density and comparisons are made between drag coefficients, pressure distributions, velocity profiles, and profiles of the turbulence working variable. The square duct shows small variation due to element type or the spatial accuracy of turbulence model convection. It is demonstrated that the finite-element scheme on tetrahedral meshes yields similar accuracy as the finite- volume scheme on mixed-element and hexahedral grids, and demonstrates less sensitivity to the mesh topology (biased tetrahedral grids) than the finite-volume scheme.

Dental application of novel finite element analysis software for three-dimensional finite element modeling of a dentulous mandible from its computed tomography images.

PubMed

Nakamura, Keiko; Tajima, Kiyoshi; Chen, Ker-Kong; Nagamatsu, Yuki; Kakigawa, Hiroshi; Masumi, Shin-ich

2013-12-01

This study focused on the application of novel finite-element analysis software for constructing a finite-element model from the computed tomography data of a human dentulous mandible. The finite-element model is necessary for evaluating the mechanical response of the alveolar part of the mandible, resulting from occlusal force applied to the teeth during biting. Commercially available patient-specific general computed tomography-based finite-element analysis software was solely applied to the finite-element analysis for the extraction of computed tomography data. The mandibular bone with teeth was extracted from the original images. Both the enamel and the dentin were extracted after image processing, and the periodontal ligament was created from the segmented dentin. The constructed finite-element model was reasonably accurate using a total of 234,644 nodes and 1,268,784 tetrahedral and 40,665 shell elements. The elastic moduli of the heterogeneous mandibular bone were determined from the bone density data of the computed tomography images. The results suggested that the software applied in this study is both useful and powerful for creating a more accurate three-dimensional finite-element model of a dentulous mandible from the computed tomography data without the need for any other software.

Methods for High-Order Multi-Scale and Stochastic Problems Analysis, Algorithms, and Applications

DTIC Science & Technology

2016-10-17

finite volume schemes, discontinuous Galerkin finite element method, and related methods, for solving computational fluid dynamics (CFD) problems and...approximation for finite element methods. (3) The development of methods of simulation and analysis for the study of large scale stochastic systems of...laws, finite element method, Bernstein-Bezier finite elements , weakly interacting particle systems, accelerated Monte Carlo, stochastic networks 16

A Simulation of Low and High Cycle Fatigue Failure Effects for Metal Matrix Composites Based on Innovative J2-Flow Elastoplasticity Model

PubMed Central

Wang, Zhaoling; Xiao, Heng

2017-01-01

New elastoplastic J2-flow constitutive equations at finite deformations are proposed for the purpose of simulating the fatigue failure behavior for metal matrix composites. A new, direct approach is established in a two-fold sense of unification. Namely, both low and high cycle fatigue failure effects of metal matrix composites may be simultaneously simulated for various cases of the weight percentage of reinforcing particles. Novel results are presented in four respects. First, both the yield condition and the loading–unloading conditions in a usual sense need not be involved but may be automatically incorporated into inherent features of the proposed constitutive equations; second, low-to-high cycle fatigue failure effects may be directly represented by a simple condition for asymptotic loss of the material strength, without involving any additional damage-like variables; third, both high and low cycle fatigue failure effects need not be separately treated but may be automatically derived as model predictions with a unified criterion for critical failure states, without assuming any ad hoc failure criteria; and, finally, explicit expressions for each incorporated model parameter changing with the weight percentage of reinforcing particles may be obtainable directly from appropriate test data. Numerical examples are presented for medium-to-high cycle fatigue failure effects and for complicated duplex effects from low to high cycle fatigue failure effects. Simulation results are in good agreement with experimental data. PMID:28946637

Finite element implementation of Robinson's unified viscoplastic model and its application to some uniaxial and multiaxial problems

NASA Technical Reports Server (NTRS)

Arya, V. K.; Kaufman, A.

1989-01-01

A description of the finite element implementation of Robinson's unified viscoplastic model into the General Purpose Finite Element Program (MARC) is presented. To demonstrate its application, the implementation is applied to some uniaxial and multiaxial problems. A comparison of the results for the multiaxial problem of a thick internally pressurized cylinder, obtained using the finite element implementation and an analytical solution, is also presented. The excellent agreement obtained confirms the correct finite element implementation of Robinson's model.

Finite element implementation of Robinson's unified viscoplastic model and its application to some uniaxial and multiaxial problems

NASA Technical Reports Server (NTRS)

Arya, V. K.; Kaufman, A.

1987-01-01

A description of the finite element implementation of Robinson's unified viscoplastic model into the General Purpose Finite Element Program (MARC) is presented. To demonstrate its application, the implementation is applied to some uniaxial and multiaxial problems. A comparison of the results for the multiaxial problem of a thick internally pressurized cylinder, obtained using the finite element implementation and an analytical solution, is also presented. The excellent agreement obtained confirms the correct finite element implementation of Robinson's model.

Stress-intensity factors for small surface and corner cracks in plates

NASA Technical Reports Server (NTRS)

Raju, I. S.; Atluri, S. N.; Newman, J. C., Jr.

1988-01-01

Three-dimensional finite-element and finite-alternating methods were used to obtain the stress-intensity factors for small surface and corner cracked plates subjected to remote tension and bending loads. The crack-depth-to-crack-length ratios (a/c) ranged from 0.2 to 1 and the crack-depth-to-plate-thickness ratios (a/t) ranged from 0.05 to 0.2. The performance of the finite-element alternating method was studied on these crack configurations. A study of the computational effort involved in the finite-element alternating method showed that several crack configurations could be analyzed with a single rectangular mesh idealization, whereas the conventional finite-element method requires a different mesh for each configuration. The stress-intensity factors obtained with the finite-element-alternating method agreed well (within 5 percent) with those calculated from the finite-element method with singularity elements.

An interactive graphics system to facilitate finite element structural analysis

NASA Technical Reports Server (NTRS)

Burk, R. C.; Held, F. H.

1973-01-01

The characteristics of an interactive graphics systems to facilitate the finite element method of structural analysis are described. The finite element model analysis consists of three phases: (1) preprocessing (model generation), (2) problem solution, and (3) postprocessing (interpretation of results). The advantages of interactive graphics to finite element structural analysis are defined.

Simulation of high-temperature superlocalization of plastic deformation in single-crystals of alloys with an L12 superstructure

NASA Astrophysics Data System (ADS)

Solov'eva, Yu. V.; Fakhrutdinova, Ya. D.; Starenchenko, V. A.

2015-01-01

The processes of the superlocalization of plastic deformation in L12 alloys have been studied numerically based on a combination of the model of the dislocation kinetics of the deformation-induced and heat-treatment-induced strengthening of an element of a deformable medium with the model of the mechanics of microplastic deformation described in terms of elastoplastic medium. It has been shown that the superlocalization of plastic deformation is determined by the presence of stress concentrators and by the nonmonotonic strengthening of the elements of the deformable medium. The multiple nonmonotonicity of the process of strengthening of the elementary volume of the medium can be responsible for the multiplicity of bands of microplastic localization of deformation.

Integrated transient thermal-structural finite element analysis

NASA Technical Reports Server (NTRS)

Thornton, E. A.; Dechaumphai, P.; Wieting, A. R.; Tamma, K. K.

1981-01-01

An integrated thermal structural finite element approach for efficient coupling of transient thermal and structural analysis is presented. Integrated thermal structural rod and one dimensional axisymmetric elements considering conduction and convection are developed and used in transient thermal structural applications. The improved accuracy of the integrated approach is illustrated by comparisons with exact transient heat conduction elasticity solutions and conventional finite element thermal finite element structural analyses.

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

NASA Technical Reports Server (NTRS)

Gould, Dana C.

2000-01-01

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Elasto-plastic bond mechanics of embedded fiber optic sensors in concrete under uniaxial tension with strain localization

NASA Astrophysics Data System (ADS)

Li, Qingbin; Li, Guang; Wang, Guanglun

2003-12-01

Brittleness of the glass core inside fiber optic sensors limits their practical usage, and therefore they are coated with low-modulus softer protective materials. Protective coatings absorb a portion of the strain, and hence part of the structural strain is sensed. The study reported here corrects for this error through development of a theoretical model to account for the loss of strain in the protective coating of optical fibers. The model considers the coating as an elasto-plastic material and formulates strain transfer coefficients for elastic, elasto-plastic and strain localization phases of coating deformations in strain localization in concrete. The theoretical findings were verified through laboratory experimentation. The experimental program involved fabrication of interferometric optical fiber sensors, embedding within mortar samples and tensile tests in a closed-loop servo-hydraulic testing machine. The elasto-plastic strain transfer coefficients were employed for correction of optical fiber sensor data and results were compared with those of conventional extensometers.

SEACAS Theory Manuals: Part III. Finite Element Analysis in Nonlinear Solid Mechanics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Laursen, T.A.; Attaway, S.W.; Zadoks, R.I.

1999-03-01

This report outlines the application of finite element methodology to large deformation solid mechanics problems, detailing also some of the key technological issues that effective finite element formulations must address. The presentation is organized into three major portions: first, a discussion of finite element discretization from the global point of view, emphasizing the relationship between a virtual work principle and the associated fully discrete system, second, a discussion of finite element technology, emphasizing the important theoretical and practical features associated with an individual finite element; and third, detailed description of specific elements that enjoy widespread use, providing some examples ofmore » the theoretical ideas already described. Descriptions of problem formulation in nonlinear solid mechanics, nonlinear continuum mechanics, and constitutive modeling are given in three companion reports.« less

Finite element analysis of structural engineering problems using a viscoplastic model incorporating two back stresses

NASA Technical Reports Server (NTRS)

Arya, Vinod K.; Halford, Gary R.

1993-01-01

The feasibility of a viscoplastic model incorporating two back stresses and a drag strength is investigated for performing nonlinear finite element analyses of structural engineering problems. To demonstrate suitability for nonlinear structural analyses, the model is implemented into a finite element program and analyses for several uniaxial and multiaxial problems are performed. Good agreement is shown between the results obtained using the finite element implementation and those obtained experimentally. The advantages of using advanced viscoplastic models for performing nonlinear finite element analyses of structural components are indicated.

Traction free finite elements with the assumed stress hybrid model. M.S. Thesis, 1981

NASA Technical Reports Server (NTRS)

Kafie, Kurosh

1991-01-01

An effective approach in the finite element analysis of the stress field at the traction free boundary of a solid continuum was studied. Conventional displacement and assumed stress finite elements were used in the determination of stress concentrations around circular and elliptical holes. Specialized hybrid elements were then developed to improve the satisfaction of prescribed traction boundary conditions. Results of the stress analysis indicated that finite elements which exactly satisfy the free stress boundary conditions are the most accurate and efficient in such problems. A general approach for hybrid finite elements which incorporate traction free boundaries of arbitrary geometry was formulated.

Improving finite element results in modeling heart valve mechanics.

PubMed

Earl, Emily; Mohammadi, Hadi

2018-06-01

Finite element analysis is a well-established computational tool which can be used for the analysis of soft tissue mechanics. Due to the structural complexity of the leaflet tissue of the heart valve, the currently available finite element models do not adequately represent the leaflet tissue. A method of addressing this issue is to implement computationally expensive finite element models, characterized by precise constitutive models including high-order and high-density mesh techniques. In this study, we introduce a novel numerical technique that enhances the results obtained from coarse mesh finite element models to provide accuracy comparable to that of fine mesh finite element models while maintaining a relatively low computational cost. Introduced in this study is a method by which the computational expense required to solve linear and nonlinear constitutive models, commonly used in heart valve mechanics simulations, is reduced while continuing to account for large and infinitesimal deformations. This continuum model is developed based on the least square algorithm procedure coupled with the finite difference method adhering to the assumption that the components of the strain tensor are available at all nodes of the finite element mesh model. The suggested numerical technique is easy to implement, practically efficient, and requires less computational time compared to currently available commercial finite element packages such as ANSYS and/or ABAQUS.

Phase-space finite elements in a least-squares solution of the transport equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Drumm, C.; Fan, W.; Pautz, S.

2013-07-01

The linear Boltzmann transport equation is solved using a least-squares finite element approximation in the space, angular and energy phase-space variables. The method is applied to both neutral particle transport and also to charged particle transport in the presence of an electric field, where the angular and energy derivative terms are handled with the energy/angular finite elements approximation, in a manner analogous to the way the spatial streaming term is handled. For multi-dimensional problems, a novel approach is used for the angular finite elements: mapping the surface of a unit sphere to a two-dimensional planar region and using a meshingmore » tool to generate a mesh. In this manner, much of the spatial finite-elements machinery can be easily adapted to handle the angular variable. The energy variable and the angular variable for one-dimensional problems make use of edge/beam elements, also building upon the spatial finite elements capabilities. The methods described here can make use of either continuous or discontinuous finite elements in space, angle and/or energy, with the use of continuous finite elements resulting in a smaller problem size and the use of discontinuous finite elements resulting in more accurate solutions for certain types of problems. The work described in this paper makes use of continuous finite elements, so that the resulting linear system is symmetric positive definite and can be solved with a highly efficient parallel preconditioned conjugate gradients algorithm. The phase-space finite elements capability has been built into the Sceptre code and applied to several test problems, including a simple one-dimensional problem with an analytic solution available, a two-dimensional problem with an isolated source term, showing how the method essentially eliminates ray effects encountered with discrete ordinates, and a simple one-dimensional charged-particle transport problem in the presence of an electric field. (authors)« less

Quality-assurance study of the special - purpose finite-element program - SPECTROM: I. Thermal, thermoelastic, and viscoelastic problems. [Comparison with MARC-CDC

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wagner, R.A.

1980-12-01

This comparison study involves a preliminary verification of finite element calculations. The methodology of the comparison study consists of solving four example problems with both the SPECTROM finite element program and the MARC-CDC general purpose finite element program. The results show close agreement for all example problems.

Finite Element Analysis of Particle Ionization within Carbon Nanotube Ion Micro Thruster

DTIC Science & Technology

2017-12-01

NAVAL POSTGRADUATE SCHOOL MONTEREY, CALIFORNIA THESIS Approved for public release. Distribution is unlimited. FINITE ELEMENT ...AND DATES COVERED Master’s thesis 4. TITLE AND SUBTITLE FINITE ELEMENT ANALYSIS OF PARTICLE IONIZATION WITHIN CARBON NANOTUBE ION MICRO THRUSTER 5...simulation, carbon nanotube simulation, microsatellite, finite element analysis, electric field, particle tracing 15. NUMBER OF PAGES 55 16. PRICE

The aggregated unfitted finite element method for elliptic problems

NASA Astrophysics Data System (ADS)

Badia, Santiago; Verdugo, Francesc; Martín, Alberto F.

2018-07-01

Unfitted finite element techniques are valuable tools in different applications where the generation of body-fitted meshes is difficult. However, these techniques are prone to severe ill conditioning problems that obstruct the efficient use of iterative Krylov methods and, in consequence, hinders the practical usage of unfitted methods for realistic large scale applications. In this work, we present a technique that addresses such conditioning problems by constructing enhanced finite element spaces based on a cell aggregation technique. The presented method, called aggregated unfitted finite element method, is easy to implement, and can be used, in contrast to previous works, in Galerkin approximations of coercive problems with conforming Lagrangian finite element spaces. The mathematical analysis of the new method states that the condition number of the resulting linear system matrix scales as in standard finite elements for body-fitted meshes, without being affected by small cut cells, and that the method leads to the optimal finite element convergence order. These theoretical results are confirmed with 2D and 3D numerical experiments.

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

NASA Technical Reports Server (NTRS)

Chung, T. J. (Editor); Karr, Gerald R. (Editor)

1989-01-01

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Reliability analysis of dispersion nuclear fuel elements

NASA Astrophysics Data System (ADS)

Ding, Shurong; Jiang, Xin; Huo, Yongzhong; Li, Lin an

2008-03-01

Taking a dispersion fuel element as a special particle composite, the representative volume element is chosen to act as the research object. The fuel swelling is simulated through temperature increase. The large strain elastoplastic analysis is carried out for the mechanical behaviors using FEM. The results indicate that the fission swelling is simulated successfully; the thickness increments grow linearly with burnup; with increasing of burnup: (1) the first principal stresses at fuel particles change from tensile ones to compression ones, (2) the maximum Mises stresses at the particles transfer from the centers of fuel particles to the location close to the interfaces between the matrix and the particles, their values increase with burnup; the maximum Mises stresses at the matrix exist in the middle location between the two particles near the mid-plane along the length (or width) direction, and the maximum plastic strains are also at the above region.

Use of edge-based finite elements for solving three dimensional scattering problems

NASA Technical Reports Server (NTRS)

Chatterjee, A.; Jin, J. M.; Volakis, John L.

1991-01-01

Edge based finite elements are free from drawbacks associated with node based vectorial finite elements and are, therefore, ideal for solving 3-D scattering problems. The finite element discretization using edge elements is checked by solving for the resonant frequencies of a closed inhomogeneously filled metallic cavity. Great improvements in accuracy are observed when compared to the classical node based approach with no penalty in terms of computational time and with the expected absence of spurious modes. A performance comparison between the edge based tetrahedra and rectangular brick elements is carried out and tetrahedral elements are found to be more accurate than rectangular bricks for a given storage intensity. A detailed formulation for the scattering problem with various approaches for terminating the finite element mesh is also presented.

Stabilized Finite Elements in FUN3D

NASA Technical Reports Server (NTRS)

Anderson, W. Kyle; Newman, James C.; Karman, Steve L.

2017-01-01

A Streamlined Upwind Petrov-Galerkin (SUPG) stabilized finite-element discretization has been implemented as a library into the FUN3D unstructured-grid flow solver. Motivation for the selection of this methodology is given, details of the implementation are provided, and the discretization for the interior scheme is verified for linear and quadratic elements by using the method of manufactured solutions. A methodology is also described for capturing shocks, and simulation results are compared to the finite-volume formulation that is currently the primary method employed for routine engineering applications. The finite-element methodology is demonstrated to be more accurate than the finite-volume technology, particularly on tetrahedral meshes where the solutions obtained using the finite-volume scheme can suffer from adverse effects caused by bias in the grid. Although no effort has been made to date to optimize computational efficiency, the finite-element scheme is competitive with the finite-volume scheme in terms of computer time to reach convergence.

SUPG Finite Element Simulations of Compressible Flows

NASA Technical Reports Server (NTRS)

Kirk, Brnjamin, S.

2006-01-01

The Streamline-Upwind Petrov-Galerkin (SUPG) finite element simulations of compressible flows is presented. The topics include: 1) Introduction; 2) SUPG Galerkin Finite Element Methods; 3) Applications; and 4) Bibliography.

On the effects of grid ill-conditioning in three dimensional finite element vector potential magnetostatic field computations

NASA Technical Reports Server (NTRS)

Wang, R.; Demerdash, N. A.

1990-01-01

The effects of finite element grid geometries and associated ill-conditioning were studied in single medium and multi-media (air-iron) three dimensional magnetostatic field computation problems. The sensitivities of these 3D field computations to finite element grid geometries were investigated. It was found that in single medium applications the unconstrained magnetic vector potential curl-curl formulation in conjunction with first order finite elements produce global results which are almost totally insensitive to grid geometries. However, it was found that in multi-media (air-iron) applications first order finite element results are sensitive to grid geometries and consequent elemental shape ill-conditioning. These sensitivities were almost totally eliminated by means of the use of second order finite elements in the field computation algorithms. Practical examples are given in this paper to demonstrate these aspects mentioned above.

Finite Element Modeling of Scattering from Underwater Proud and Buried Military Munitions

DTIC Science & Technology

2017-02-28

FINAL REPORT Finite Element Modeling of Scattering from Underwater Proud and Buried Military Munitions SERDP Project MR-2408 JULY 2017...solution and the red dash-dot line repre- sents the coupled finite -boundary element solution. . . . . . . . . . . . . . . . . . 11 3 The scattering...dot line represents the coupled finite -boundary element solution. . . . . . . . 11 i 4 The scattering amplitude as a function of the receiver angle for

Finite elements and finite differences for transonic flow calculations

NASA Technical Reports Server (NTRS)

Hafez, M. M.; Murman, E. M.; Wellford, L. C.

1978-01-01

The paper reviews the chief finite difference and finite element techniques used for numerical solution of nonlinear mixed elliptic-hyperbolic equations governing transonic flow. The forms of the governing equations for unsteady two-dimensional transonic flow considered are the Euler equation, the full potential equation in both conservative and nonconservative form, the transonic small-disturbance equation in both conservative and nonconservative form, and the hodograph equations for the small-disturbance case and the full-potential case. Finite difference methods considered include time-dependent methods, relaxation methods, semidirect methods, and hybrid methods. Finite element methods include finite element Lax-Wendroff schemes, implicit Galerkin method, mixed variational principles, dual iterative procedures, optimal control methods and least squares.

Cognitive-graphic method for constructing of hierarchical forms of basic functions of biquadratic finite element

NASA Astrophysics Data System (ADS)

Astionenko, I. O.; Litvinenko, O. I.; Osipova, N. V.; Tuluchenko, G. Ya.; Khomchenko, A. N.

2016-10-01

Recently the interpolation bases of the hierarchical type have been used for the problem solving of the approximation of multiple arguments functions (such as in the finite-element method). In this work the cognitive graphical method of constructing of the hierarchical form bases on the serendipity finite elements is suggested, which allowed to get the alternative bases on a biquadratic finite element from the serendipity family without internal knots' inclusion. The cognitive-graphic method allowed to improve the known interpolation procedure of Taylor and to get the modified elements with irregular arrangement of knots. The proposed procedures are universal and are spread in the area of finite-elements.

Engine dynamic analysis with general nonlinear finite element codes. Part 2: Bearing element implementation overall numerical characteristics and benchmaking

NASA Technical Reports Server (NTRS)

Padovan, J.; Adams, M.; Fertis, J.; Zeid, I.; Lam, P.

1982-01-01

Finite element codes are used in modelling rotor-bearing-stator structure common to the turbine industry. Engine dynamic simulation is used by developing strategies which enable the use of available finite element codes. benchmarking the elements developed are benchmarked by incorporation into a general purpose code (ADINA); the numerical characteristics of finite element type rotor-bearing-stator simulations are evaluated through the use of various types of explicit/implicit numerical integration operators. Improving the overall numerical efficiency of the procedure is improved.

Extension to linear dynamics for hybrid stress finite element formulation based on additional displacements

NASA Astrophysics Data System (ADS)

Sumihara, K.

Based upon legitimate variational principles, one microscopic-macroscopic finite element formulation for linear dynamics is presented by Hybrid Stress Finite Element Method. The microscopic application of Geometric Perturbation introduced by Pian and the introduction of infinitesimal limit core element (Baby Element) have been consistently combined according to the flexible and inherent interpretation of the legitimate variational principles initially originated by Pian and Tong. The conceptual development based upon Hybrid Finite Element Method is extended to linear dynamics with the introduction of physically meaningful higher modes.

Finite Macro-Element Mesh Deformation in a Structured Multi-Block Navier-Stokes Code

NASA Technical Reports Server (NTRS)

Bartels, Robert E.

2005-01-01

A mesh deformation scheme is developed for a structured multi-block Navier-Stokes code consisting of two steps. The first step is a finite element solution of either user defined or automatically generated macro-elements. Macro-elements are hexagonal finite elements created from a subset of points from the full mesh. When assembled, the finite element system spans the complete flow domain. Macro-element moduli vary according to the distance to the nearest surface, resulting in extremely stiff elements near a moving surface and very pliable elements away from boundaries. Solution of the finite element system for the imposed boundary deflections generally produces smoothly varying nodal deflections. The manner in which distance to the nearest surface has been found to critically influence the quality of the element deformation. The second step is a transfinite interpolation which distributes the macro-element nodal deflections to the remaining fluid mesh points. The scheme is demonstrated for several two-dimensional applications.

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

NASA Technical Reports Server (NTRS)

Wu, R. W.; Witmer, E. A.

1972-01-01

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Investigations into the tensile failure of doubly-convex cylindrical tablets under diametral loading using finite element methodology.

PubMed

Podczeck, Fridrun; Drake, Kevin R; Newton, J Michael

2013-09-15

In the literature various solutions exist for the calculation of the diametral compression tensile strength of doubly-convex tablets and each approach is based on experimental data obtained from single materials (gypsum, microcrystalline cellulose) only. The solutions are represented by complex equations and further differ for elastic and elasto-plastic behaviour of the compacts. The aim of this work was to develop a general equation that is applicable independently of deformation behaviour and which is based on simple tablet dimensions such as diameter and total tablet thickness only. With the help of 3D-FEM analysis the tensile failure stress of doubly-convex tables with central cylinder to total tablet thickness ratios W/D between 0.06 and 0.50 and face-curvature ratios D/R between 0.25 and 1.85 were evaluated. Both elastic and elasto-plastic deformation behaviour were considered. The results of 80 individual simulations were combined and showed that the tensile failure stress σt of doubly-convex tablets can be calculated from σt=(2P/πDW)(W/T)=2P/πDT with P being the failure load, D the diameter, W the central cylinder thickness, and T the total thickness of the tablet. This equation converts into the standard Brazilian equation (σt=2P/πDW) when W equals T, i.e. is equally valid for flat cylindrical tablets. In practice, the use of this new equation removes the need for complex measurements of tablet dimensions, because it only requires values for diameter and total tablet thickness. It also allows setting of standards for the mechanical strength of doubly-convex tablets. The new equation holds both for elastic and elasto-plastic deformation behaviour of the tablets under load. It is valid for all combinations of W/D-ratios between 0.06 and 0.50 with D/R-ratios between 0.00 and 1.85 except for W/D=0.50 in combination with D/R-ratios of 1.85 and 1.43 and for W/D-ratios of 0.40 and 0.30 in combination with D/R=1.85. FEM-analysis indicated a tendency to failure by capping or even more complex failure patterns in these exceptional cases. The FEM-results further indicated that in general W/D-ratios between 0.15 and 0.20 are favourable when the overall size and shape of the tablets is modified to give maximum tablet tensile strength. However, the maximum tensile stress of doubly-convex tablets will never exceed that of a flat-face cylindrical tablet of similar W/D-ratio. The lowest tensile stress depends on the W/D-ratio. For the thinnest central cylinder thickness, this minimum stress occurs at D/R=0.50; for W/D-ratios between 0.10 and 0.20 the D/R-ratio for the minimum tensile stress increases to 0.67, and for all other central cylinder thicknesses the minimum tensile stress is found at D/R=1.00. Copyright © 2013 Elsevier B.V. All rights reserved.

Wave Scattering in Heterogeneous Media using the Finite Element Method

DTIC Science & Technology

2016-10-21

AFRL-AFOSR-JP-TR-2016-0086 Wave Scattering in Heterogeneous Media using the Finite Element Method Chiruvai Vendhan INDIAN INSTITUTE OF TECHNOLOGY...Scattering in Heterogeneous Media using the Finite Element Method 5a.  CONTRACT NUMBER 5b.  GRANT NUMBER FA2386-12-1-4026 5c.  PROGRAM ELEMENT NUMBER 61102F 6...14.  ABSTRACT The primary aim of this study is to develop a finite element model for elastic scattering by axisymmetric bodies submerged in a

Experimental and Modelling Investigations of the Coupled Elastoplastic Damage of a Quasi-brittle Rock

NASA Astrophysics Data System (ADS)

Zhang, Jiu-Chang

2018-02-01

Triaxial compression tests are conducted on a quasi-brittle rock, limestone. The analyses show that elastoplastic deformation is coupled with damage. Based on the experimental investigation, a coupled elastoplastic damage model is developed within the framework of irreversible thermodynamics. The coupling effects between the plastic and damage dissipations are described by introducing an isotropic damage variable into the elastic stiffness and yield criterion. The novelty of the model is in the description of the thermodynamic force associated with damage, which is formulated as a state function of both elastic and plastic strain energies. The latter gives a full consideration on the comprehensive effects of plastic strain and stress changing processes in rock material on the development of damage. The damage criterion and potential are constructed to determine the onset and evolution of damage variable. The return mapping algorithms of the coupled model are deduced for three different inelastic corrections. Comparisons between test data and numerical simulations show that the coupled elastoplastic damage model is capable of describing the main mechanical behaviours of the quasi-brittle rock.

Micromechanics and effective elastoplastic behavior of two-phase metal matrix composites

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ju, J.W.; Chen, T.M.

A micromechanical framework is presented to predict effective (overall) elasto-(visco-)plastic behavior of two-phase particle-reinforced metal matrix composites (PRMMC). In particular, the inclusion phase (particle) is assumed to be elastic and the matrix material is elasto-(visco-)plastic. Emanating from Ju and Chen's (1994a,b) work on effective elastic properties of composites containing many randomly dispersed inhomogeneities, effective elastoplastic deformations and responses of PRMMC are estimated by means of the effective yield criterion'' derived micromechanically by considering effects due to elastic particles embedded in the elastoplastic matrix. The matrix material is elastic or plastic, depending on local stress and deformation, and obeys general plasticmore » flow rule and hardening law. Arbitrary (general) loadings and unloadings are permitted in the framework through the elastic predictor-plastic corrector two-step operator splitting methodology. The proposed combined micromechanical and computational approach allows one to estimate overall elastoplastic responses of PRMMCs by accounting for the microstructural information (such as the spatial distribution and micro-geometry of particles), elastic properties of constituent phases, and the plastic behavior of the matrix-only materials.« less

On the determination of stress fields and displacements in a thin elastoplastic plate containing elastic inclusion - a shim

NASA Astrophysics Data System (ADS)

Kovalev, A. V.; Rusina, E. Y.; Yakovlev, A. Y.

2018-03-01

The paper is devoted to the determination of the stress-strain state of a mechanical structure, which consists of a thin infinite elastoplastic plate with a hole and a continuous fine elastic inclusion. The complexity of this problem lies in the fact that the shape of the boundary between the elastic and plastic zones in the plate is not known in advance. The small parameter method is used as the solution method, while the small parameter characterizes the deviation of the shape of the contour from the circle and the perturbation of external static boundary conditions. As the zero solution, the axisymmetric elastoplastic state of the plate with inclusion is chosen. Two variants of inclusion fixation in a plate are considered: inclusion is enclosed with tension or soldered. As a result of solving the problem within the framework of ideal plasticity, the distribution of the stress and displacement fields and the shape of the elastoplastic boundary are obtained. To illustrate the case of a plane-stressed state, the results of a numerical experiment on the mathematical model obtained are presented.

3D viscous time dependent analysis of a slow moving landslide by finite elements

NASA Astrophysics Data System (ADS)

Bru, Guadalupe; Fernández-Merodo, Josè Antonio; García-Davalillo, Juan Carlos; Herrera, Gerardo

2017-04-01

A methodology to study the cinematic behavior in time and 3D effects of slow moving landslides is presented by its application to the Portalet landslide, located in the Spanish Central Pyrenees. The area is characterize by the presence of several complex landslides triggered by glacial retreat and developed in weathered Devonian and Carboniferous slate materials. In summer 2004, the toe of two paleolanslides was excavated in order to construct a parking area, which reactivated the preexistent rupture surfaces and generated a new and smaller roto-translational landslide of about 5·105m3. Works were paralyzed and the road A-136 connecting Spain and France was temporary closed. In 2006 stabilizing solutions were performed, although surface and deep monitoring data reveal that the landslide is currently active. In this work the cinematic behavior of the Portalet parking landslide has been reproduced since the excavation in 2004 until 2016 by an advanced 3D finite element model (FEM) analysis. A detailed 3D litho-stratigraphic geometry has been defined using the information available from previous works(1-3): real topography, geological profiles, ground water level and material properties. The first step has been a stability analysis using the shear stress reduction (SSR) technique to calibrate the value of the friction angle of the soil layer where the rupture surface develops. In this step the sensibility of the mesh size has been studied, as is a critical parameter. Secondly, stability analysis results have been verified by simulating the parking excavation with a static analysis using Mohr-Coulomb elastoplastic failure criteria. In the last step the behavior of the landslide has been recreated using a hydromechanic coupled formulation for displacements and interstitial water pressure (u-pw), a simple elevation ground water model calculated from daily rainfall and a Perzyna viscous constitutive model of the solid skeleton which represent the creep detected by monitoring. The influence of the fluidity factor of the material has been explored, obtaining the best fitting value. Results show the advantages of performing a 3D analysis, mainly because they reflect the spatial extension of the landslide without the need of assuming a critical profile and allow to study the interaction of differential movements that occur within the same landslide. This methodology can be applied to other slow moving landslides where geological structure and geotechnical data are available. REFERENCES 1. Fernández-Merodo, J. et al. Modelling the Portalet landslide mobility (Formigal, Spain). iEMSs 2008: International Congress on Environmental Modelling and Software. International Environmental Modelling and Software Society (iEMSs) (2008). 2. Herrera, G. et al. A landslide forecasting model using ground based SAR data: The Portalet case study. Engineering Geology 105, 220-230 (2009). 3. Fernández-Merodo, J., García-Davalillo, J., Herrera, G., Mira, P. & Pastor, M. 2D viscoplastic finite element modelling of slow landslides: the Portalet case study (Spain). Landslides 11, 29-42 (2014).

Carbon-dioxide storage in the subsurface: towards an understanding of crack development in the cap-rock including phase transition processes

NASA Astrophysics Data System (ADS)

Häberle, K.; Ehlers, W.

2012-04-01

Supercritical CO2 can be injected into deep saline aquifers to reduce the amount of CO2 in the atmosphere and thus, lessen the impact on the global warming. Qualified reservoirs should be in a sufficient depth to guarantee the thermodynamical environment for the supercritical state of CO2. Furthermore, an impermeable cap-rock layer must confine the reservoir layer, in order to collect the CO2 in the desired region. In CO2 storage it is crucial to guarantee the safety of the storage site and to eliminate possibilities of leakage. Therefore, deformation processes of the rock matrix and the cap-rock layer, which might be induced by the high pressure injection of CO2, must be investigated. The increase in stress may also cause crack development in the cap-rock layer. These could either be new developing fractures or the break-up of already existing but cemented cracks or faults. If such cracks occur, CO2 could migrate to shallower regions where the temperature and pressure cannot support the supercritical condition of the CO2 anymore. Thus, it is important to describe the phase transition process between supercritical, liquid and gaseous CO2. This requires a proper understanding of the thermodynamical behaviour of CO2 within the reservoir. The Theory of Porous Media (TPM) provides a useful continuum-mechanical basis to describe real natural systems in a thermodynamically consistent way. Hence, the TPM is applied to model multiphasic flow of CO2 and water and include elasto-plastic solid deformations of the porous matrix. The Peng-Robinson equation is implemented as a cubic equation of state to describe the phase behaviour of CO2 in the liquid, gaseous and supercritical region. However, in the two-phase region the isotherms show a horizontal section and kinks at the boiling and vapour curve. This cannot be represented by a continuously differentiable function such as the Peng-Robinson equation. To circumvent this problem, the Antoine equation provides additional information by defining the saturation pressure for a given temperature. The injection of CO2 will increase the reservoir pressure which then will cause solid deformations. The extended Finite Element Method (XFEM) will be used to account for the discontinuities arising from crack development due to these solid deformations. The XFEM bears the advantage that the finite element mesh must not be adapted to the crack. Instead, to describe the discontinuity of the crack, the field quantities are locally enriched by defining additional degrees of freedom at the intersected finite elements. Herein, special attention has to be paid to the matrix-fracture interaction of the fluid phases. Numerical examples are performed to investigate the injection of CO2 into a saline aquifer. These are computed with the FEM program PANDAS, which allows solutions of strongly coupled multiphasic problems in deformable porous media.

Improved finite-element methods for rotorcraft structures

NASA Technical Reports Server (NTRS)

Hinnant, Howard E.

1991-01-01

An overview of the research directed at improving finite-element methods for rotorcraft airframes is presented. The development of a modification to the finite element method which eliminates interelement discontinuities is covered. The following subject areas are discussed: geometric entities, interelement continuity, dependent rotational degrees of freedom, and adaptive numerical integration. This new methodology is being implemented as an anisotropic, curvilinear, p-version, beam, shell, and brick finite element program.

Refinement of Out of Circularity and Thickness Measurements of a Cylinder for Finite Element Analysis

DTIC Science & Technology

2016-09-01

UNCLASSIFIED UNCLASSIFIED Refinement of Out of Circularity and Thickness Measurements of a Cylinder for Finite Element Analysis...significant effect on the collapse strength and must be accurately represented in finite element analysis to obtain accurate results. Often it is necessary...to interpolate measurements from a relatively coarse grid to a refined finite element model and methods that have wide general acceptance are

Using Finite Element Method to Estimate the Material Properties of a Bearing Cage

DTIC Science & Technology

2018-02-01

UNCLASSIFIED UNCLASSIFIED AD-E403 988 Technical Report ARMET-TR-17035 USING FINITE ELEMENT METHOD TO ESTIMATE THE MATERIAL...TITLE AND SUBTITLE USING FINITE ELEMENT METHOD TO ESTIMATE THE MATERIAL PROPERTIES OF A BEARING CAGE 5a. CONTRACT NUMBER 5b. GRANT...specifications of non-metallic bearing cages are typically not supplied by the manufacturer. In order to setup a finite element analysis of a

Artificial Boundary Conditions for Finite Element Model Update and Damage Detection

DTIC Science & Technology

2017-03-01

BOUNDARY CONDITIONS FOR FINITE ELEMENT MODEL UPDATE AND DAMAGE DETECTION by Emmanouil Damanakis March 2017 Thesis Advisor: Joshua H. Gordis...REPORT TYPE AND DATES COVERED Master’s thesis 4. TITLE AND SUBTITLE ARTIFICIAL BOUNDARY CONDITIONS FOR FINITE ELEMENT MODEL UPDATE AND DAMAGE DETECTION...release. Distribution is unlimited. 12b. DISTRIBUTION CODE 13. ABSTRACT (maximum 200 words) In structural engineering, a finite element model is often

Angular Random Walk Estimation of a Time-Domain Switching Micromachined Gyroscope

DTIC Science & Technology

2016-10-19

1 2. PARAMETRIC SYSTEM IDENTIFICATION BASED ON TIME-DOMAIN SWITCHING ........ 2 3. FINITE ELEMENT MODELING OF RESONATOR...8 3. FINITE ELEMENT MODELING OF RESONATOR This section details basic finite element modeling of the resonator used with the TDSMG. While it...Based on finite element simulations of the employed resonator, it is found that the effects of thermomechanical noise is on par with 10 ps of timing

Verification of finite element analysis of fixed partial denture with in vitro electronic strain measurement.

PubMed

Wang, Gaoqi; Zhang, Song; Bian, Cuirong; Kong, Hui

2016-01-01

The purpose of the study was to verify the finite element analysis model of three-unite fixed partial denture with in vitro electronic strain analysis and analyze clinical situation with the verified model. First, strain gauges were attached to the critical areas of a three-unit fixed partial denture. Strain values were measured under 300 N load perpendicular to the occlusal plane. Secondly, a three-dimensional finite element model in accordance with the electronic strain analysis experiment was constructed from the scanning data. And the strain values obtained by finite element analysis and in vitro measurements were compared. Finally, the clinical destruction of the fixed partial denture was evaluated with the verified finite element analysis model. There was a mutual agreement and consistency between the finite element analysis results and experimental data. The finite element analysis revealed that failure will occur in the veneer layer on buccal surface of the connector under occlusal force of 570 N. The results indicate that the electronic strain analysis is an appropriate and cost saving method to verify the finite element model. The veneer layer on buccal surface of the connector is the weakest area in the fixed partial denture. Copyright © 2015 Japan Prosthodontic Society. Published by Elsevier Ltd. All rights reserved.

On conforming mixed finite element methods for incompressible viscous flow problems

NASA Technical Reports Server (NTRS)

Gunzburger, M. D; Nicolaides, R. A.; Peterson, J. S.

1982-01-01

The application of conforming mixed finite element methods to obtain approximate solutions of linearized Navier-Stokes equations is examined. Attention is given to the convergence rates of various finite element approximations of the pressure and the velocity field. The optimality of the convergence rates are addressed in terms of comparisons of the approximation convergence to a smooth solution in relation to the best approximation available for the finite element space used. Consideration is also devoted to techniques for efficient use of a Gaussian elimination algorithm to obtain a solution to a system of linear algebraic equations derived by finite element discretizations of linear partial differential equations.

Ablative Thermal Response Analysis Using the Finite Element Method

NASA Technical Reports Server (NTRS)

Dec John A.; Braun, Robert D.

2009-01-01

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Infinite Possibilities for the Finite Element.

ERIC Educational Resources Information Center

Finlayson, Bruce A.

1981-01-01

Describes the uses of finite element methods in solving problems of heat transfer, fluid flow, etc. Suggests that engineers should know the general concepts and be able to apply the principles of finite element methods. (Author/WB)

High Performance Computing Technologies for Modeling the Dynamics and Dispersion of Ice Chunks in the Arctic Ocean

DTIC Science & Technology

2016-08-23

SECURITY CLASSIFICATION OF: Hybrid finite element / finite volume based CaMEL shallow water flow solvers have been successfully extended to study wave...effects on ice floes in a simplified 10 sq-km ocean domain. Our solver combines the merits of both the finite element and finite volume methods and...ES) U.S. Army Research Office P.O. Box 12211 Research Triangle Park, NC 27709-2211 sea ice dynamics, shallow water, finite element , finite volume

Vectorial finite elements for solving the radiative transfer equation

NASA Astrophysics Data System (ADS)

Badri, M. A.; Jolivet, P.; Rousseau, B.; Le Corre, S.; Digonnet, H.; Favennec, Y.

2018-06-01

The discrete ordinate method coupled with the finite element method is often used for the spatio-angular discretization of the radiative transfer equation. In this paper we attempt to improve upon such a discretization technique. Instead of using standard finite elements, we reformulate the radiative transfer equation using vectorial finite elements. In comparison to standard finite elements, this reformulation yields faster timings for the linear system assemblies, as well as for the solution phase when using scattering media. The proposed vectorial finite element discretization for solving the radiative transfer equation is cross-validated against a benchmark problem available in literature. In addition, we have used the method of manufactured solutions to verify the order of accuracy for our discretization technique within different absorbing, scattering, and emitting media. For solving large problems of radiation on parallel computers, the vectorial finite element method is parallelized using domain decomposition. The proposed domain decomposition method scales on large number of processes, and its performance is unaffected by the changes in optical thickness of the medium. Our parallel solver is used to solve a large scale radiative transfer problem of the Kelvin-cell radiation.

A combined registration and finite element analysis method for fast estimation of intraoperative brain shift; phantom and animal model study.

PubMed

Mohammadi, Amrollah; Ahmadian, Alireza; Rabbani, Shahram; Fattahi, Ehsan; Shirani, Shapour

2017-12-01

Finite element models for estimation of intraoperative brain shift suffer from huge computational cost. In these models, image registration and finite element analysis are two time-consuming processes. The proposed method is an improved version of our previously developed Finite Element Drift (FED) registration algorithm. In this work the registration process is combined with the finite element analysis. In the Combined FED (CFED), the deformation of whole brain mesh is iteratively calculated by geometrical extension of a local load vector which is computed by FED. While the processing time of the FED-based method including registration and finite element analysis was about 70 s, the computation time of the CFED was about 3.2 s. The computational cost of CFED is almost 50% less than similar state of the art brain shift estimators based on finite element models. The proposed combination of registration and structural analysis can make the calculation of brain deformation much faster. Copyright © 2016 John Wiley & Sons, Ltd.

Global-Local Finite Element Analysis of Bonded Single-Lap Joints

NASA Technical Reports Server (NTRS)

Kilic, Bahattin; Madenci, Erdogan; Ambur, Damodar R.

2004-01-01

Adhesively bonded lap joints involve dissimilar material junctions and sharp changes in geometry, possibly leading to premature failure. Although the finite element method is well suited to model the bonded lap joints, traditional finite elements are incapable of correctly resolving the stress state at junctions of dissimilar materials because of the unbounded nature of the stresses. In order to facilitate the use of bonded lap joints in future structures, this study presents a finite element technique utilizing a global (special) element coupled with traditional elements. The global element includes the singular behavior at the junction of dissimilar materials with or without traction-free surfaces.

The constraint method: A new finite element technique. [applied to static and dynamic loads on plates

NASA Technical Reports Server (NTRS)

Tsai, C.; Szabo, B. A.

1973-01-01

An approch to the finite element method which utilizes families of conforming finite elements based on complete polynomials is presented. Finite element approximations based on this method converge with respect to progressively reduced element sizes as well as with respect to progressively increasing orders of approximation. Numerical results of static and dynamic applications of plates are presented to demonstrate the efficiency of the method. Comparisons are made with plate elements in NASTRAN and the high-precision plate element developed by Cowper and his co-workers. Some considerations are given to implementation of the constraint method into general purpose computer programs such as NASTRAN.

Element-topology-independent preconditioners for parallel finite element computations

NASA Technical Reports Server (NTRS)

Park, K. C.; Alexander, Scott

1992-01-01

A family of preconditioners for the solution of finite element equations are presented, which are element-topology independent and thus can be applicable to element order-free parallel computations. A key feature of the present preconditioners is the repeated use of element connectivity matrices and their left and right inverses. The properties and performance of the present preconditioners are demonstrated via beam and two-dimensional finite element matrices for implicit time integration computations.

Flexible parallel implicit modelling of coupled thermal-hydraulic-mechanical processes in fractured rocks

NASA Astrophysics Data System (ADS)

Cacace, Mauro; Jacquey, Antoine B.

2017-09-01

Theory and numerical implementation describing groundwater flow and the transport of heat and solute mass in fully saturated fractured rocks with elasto-plastic mechanical feedbacks are developed. In our formulation, fractures are considered as being of lower dimension than the hosting deformable porous rock and we consider their hydraulic and mechanical apertures as scaling parameters to ensure continuous exchange of fluid mass and energy within the fracture-solid matrix system. The coupled system of equations is implemented in a new simulator code that makes use of a Galerkin finite-element technique. The code builds on a flexible, object-oriented numerical framework (MOOSE, Multiphysics Object Oriented Simulation Environment) which provides an extensive scalable parallel and implicit coupling to solve for the multiphysics problem. The governing equations of groundwater flow, heat and mass transport, and rock deformation are solved in a weak sense (either by classical Newton-Raphson or by free Jacobian inexact Newton-Krylow schemes) on an underlying unstructured mesh. Nonlinear feedbacks among the active processes are enforced by considering evolving fluid and rock properties depending on the thermo-hydro-mechanical state of the system and the local structure, i.e. degree of connectivity, of the fracture system. A suite of applications is presented to illustrate the flexibility and capability of the new simulator to address problems of increasing complexity and occurring at different spatial (from centimetres to tens of kilometres) and temporal scales (from minutes to hundreds of years).

Triple shape memory polymers by 4D printing

NASA Astrophysics Data System (ADS)

Bodaghi, M.; Damanpack, A. R.; Liao, W. H.

2018-06-01

This article aims at introducing triple shape memory polymers (SMPs) by four-dimensional (4D) printing technology and shaping adaptive structures for mechanical/bio-medical devices. The main approach is based on arranging hot–cold programming of SMPs with fused decomposition modeling technology to engineer adaptive structures with triple shape memory effect (SME). Experiments are conducted to characterize elasto-plastic and hyper-elastic thermo-mechanical material properties of SMPs in low and high temperatures at large deformation regime. The feasibility of the dual and triple SMPs with self-bending features is demonstrated experimentally. It is advantageous in situations either where it is desired to perform mechanical manipulations on the 4D printed objects for specific purposes or when they experience cold programming inevitably before activation. A phenomenological 3D constitutive model is developed for quantitative understanding of dual/triple SME of SMPs fabricated by 4D printing in the large deformation range. Governing equations of equilibrium are established for adaptive structures on the basis of the nonlinear Green–Lagrange strains. They are then solved by developing a finite element approach along with an elastic-predictor plastic-corrector return map procedure accomplished by the Newton–Raphson method. The computational tool is applied to simulate dual/triple SMP structures enabled by 4D printing and explore hot–cold programming mechanisms behind material tailoring. It is shown that the 4D printed dual/triple SMPs have great potential in mechanical/bio-medical applications such as self-bending gripers/stents and self-shrinking/tightening staples.

Stress and Friction Distribution around Slab Corner in Continuous Casting Mold with Different Corner Structures

NASA Astrophysics Data System (ADS)

Yu, Sheng; Long, Mujun; Chen, Huabiao; Chen, Dengfu; Liu, Tao; Duan, Huamei; Cao, Junsheng

2018-06-01

The non-uniform friction and thermal stress in the mold are important as causes of the transverse cracks around strand corner. To analyze the stress distribution features around strand corner, a three-dimensional thermo-elastoplastic finite-element mold model with different corner structures (right-angle, big-chamfer, multi-chamfer, and fillet) was established. The temperature field in the mold was indirectly coupled through a three-dimensional fluid flow and heat transfer model. In addition, the non-uniform mold friction stress loaded on the strand surface was calculated through a friction model. The results show that the stress distribution on the shell is similar to the temperature distribution. The stress concentration appears in the strand corner and the lower part of wide face. The friction stress enhances the corner stress around the edge of the air-gap. For chamfered molds, the stress around the corner between the wide face and chamfer face is larger than that between the narrow face and chamfer face. Around the corner region, both the stress peak and the area of the large stress zone of the right-angle strand are the largest, while those of big-chamfered, multi-chamfered, and fillet strands decrease in that order. The stress peak position of the chamfered strands is closer to the mold exit than that of the right-angle strand. Compared with the use of the right-angle mold, the application of chamfered molds is able to reduce the stress concentration around the strand corner.

Stress and Friction Distribution around Slab Corner in Continuous Casting Mold with Different Corner Structures

NASA Astrophysics Data System (ADS)

Yu, Sheng; Long, Mujun; Chen, Huabiao; Chen, Dengfu; Liu, Tao; Duan, Huamei; Cao, Junsheng

2018-02-01

The non-uniform friction and thermal stress in the mold are important as causes of the transverse cracks around strand corner. To analyze the stress distribution features around strand corner, a three-dimensional thermo-elastoplastic finite-element mold model with different corner structures (right-angle, big-chamfer, multi-chamfer, and fillet) was established. The temperature field in the mold was indirectly coupled through a three-dimensional fluid flow and heat transfer model. In addition, the non-uniform mold friction stress loaded on the strand surface was calculated through a friction model. The results show that the stress distribution on the shell is similar to the temperature distribution. The stress concentration appears in the strand corner and the lower part of wide face. The friction stress enhances the corner stress around the edge of the air-gap. For chamfered molds, the stress around the corner between the wide face and chamfer face is larger than that between the narrow face and chamfer face. Around the corner region, both the stress peak and the area of the large stress zone of the right-angle strand are the largest, while those of big-chamfered, multi-chamfered, and fillet strands decrease in that order. The stress peak position of the chamfered strands is closer to the mold exit than that of the right-angle strand. Compared with the use of the right-angle mold, the application of chamfered molds is able to reduce the stress concentration around the strand corner.

idRHa+ProMod - Rail Hardening Control System

NASA Astrophysics Data System (ADS)

Ferro, L.

2016-03-01

idRHa+ProMod is the process control system developed by Primetals Technologies to foresee the thermo-mechanical evolution and micro-structural composition of rail steels subjected to slack quenching into idRHa+ Rail Hardening equipments in a simulation environment. This tool can be used both off-line or in-line, giving the user the chance to test and study the best cooling strategies or letting the automatic control system free to adjust the proper cooling recipe. Optimization criteria have been tailored in order to determine the best cooling conditions according to the metallurgical requirements imposed by the main rail standards and also taking into account the elastoplastic bending phenomena occurring during all stages of the head hardening process. The computational core of idRHa+ProMod is a thermal finite element procedure coupled with special algorithms developed to work out the main thermo-physical properties of steel, to predict the non-isothermal austenite decomposition into all the relevant phases and subsequently to evaluate the amount of latent heat of transformation released, the compound thermal expansion coefficient and the amount of plastic deformation in the material. Air mist and air blades boundary conditions have been carefully investigated by means of pilot plant tests aimed to study the jet impingement on rail surfaces and the cooling efficiency at all working conditions. Heat transfer coefficients have been further checked and adjusted directly on field during commissioning. idRHa+ is a trademark of Primetals Technologies Italy Srl

A progress report on estuary modeling by the finite-element method

USGS Publications Warehouse

Gray, William G.

1978-01-01

Various schemes are investigated for finite-element modeling of two-dimensional surface-water flows. The first schemes investigated combine finite-element spatial discretization with split-step time stepping schemes that have been found useful in finite-difference computations. Because of the large number of numerical integrations performed in space and the large sparse matrices solved, these finite-element schemes were found to be economically uncompetitive with finite-difference schemes. A very promising leapfrog scheme is proposed which, when combined with a novel very fast spatial integration procedure, eliminates the need to solve any matrices at all. Additional problems attacked included proper propagation of waves and proper specification of the normal flow-boundary condition. This report indicates work in progress and does not come to a definitive conclusion as to the best approach for finite-element modeling of surface-water problems. The results presented represent findings obtained between September 1973 and July 1976. (Woodard-USGS)

Finite Element Models and Properties of a Stiffened Floor-Equipped Composite Cylinder

NASA Technical Reports Server (NTRS)

Grosveld, Ferdinand W.; Schiller, Noah H.; Cabell, Randolph H.

2010-01-01

Finite element models were developed of a floor-equipped, frame and stringer stiffened composite cylinder including a coarse finite element model of the structural components, a coarse finite element model of the acoustic cavities above and below the beam-supported plywood floor, and two dense models consisting of only the structural components. The report summarizes the geometry, the element properties, the material and mechanical properties, the beam cross-section characteristics, the beam element representations and the boundary conditions of the composite cylinder models. The expressions used to calculate the group speeds for the cylinder components are presented.

IFEMS, an Interactive Finite Element Modeling System Using a CAD/CAM System

NASA Technical Reports Server (NTRS)

Mckellip, S.; Schuman, T.; Lauer, S.

1980-01-01

A method of coupling a CAD/CAM system with a general purpose finite element mesh generator is described. The three computer programs which make up the interactive finite element graphics system are discussed.

Recent Development of Multigrid Algorithms for Mixed and Noncomforming Methods for Second Order Elliptical Problems

NASA Technical Reports Server (NTRS)

Chen, Zhangxin; Ewing, Richard E.

1996-01-01

Multigrid algorithms for nonconforming and mixed finite element methods for second order elliptic problems on triangular and rectangular finite elements are considered. The construction of several coarse-to-fine intergrid transfer operators for nonconforming multigrid algorithms is discussed. The equivalence between the nonconforming and mixed finite element methods with and without projection of the coefficient of the differential problems into finite element spaces is described.

Development of an Anatomically Accurate Finite Element Human Ocular Globe Model for Blast-Related Fluid-Structure Interaction Studies

DTIC Science & Technology

2017-02-01

ARL-TR-7945 ● FEB 2017 US Army Research Laboratory Development of an Anatomically Accurate Finite Element Human Ocular Globe...ARL-TR-7945 ● FEB 2017 US Army Research Laboratory Development of an Anatomically Accurate Finite Element Human Ocular Globe Model... Finite Element Human Ocular Globe Model for Blast-Related Fluid-Structure Interaction Studies 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM

Finite Element Barotropic Model for the Indian and Western Pacific OceanBasin: Tidal Model Data Comparisons and Sensitivities

DTIC Science & Technology

2018-01-11

From - To) 01/11/2018 Final Technical Report June 01 2016 - Dec 30 2017 4. TITLE AND SUBTITLE Sa. CONTRACT NUMBER Finite - Element Barotropic Model...grid finite - element barotropic fully hydrodynamic model in order to understand the shallow-water dynamics of the Indian Ocean and Western Pacific Ocean...dissipative dissipative processes are explored. 15. SUBJECTTERMS finite - element , unstructured grid, barotropic tides, bathymetry, internal tide

MSC products for the simulation of tire behavior

NASA Technical Reports Server (NTRS)

Muskivitch, John C.

1995-01-01

The modeling of tires and the simulation of tire behavior are complex problems. The MacNeal-Schwendler Corporation (MSC) has a number of finite element analysis products that can be used to address the complexities of tire modeling and simulation. While there are many similarities between the products, each product has a number of capabilities that uniquely enable it to be used for a specific aspect of tire behavior. This paper discusses the following programs: (1) MSC/NASTRAN - general purpose finite element program for linear and nonlinear static and dynamic analysis; (2) MSC/ADAQUS - nonlinear statics and dynamics finite element program; (3) MSC/PATRAN AFEA (Advanced Finite Element Analysis) - general purpose finite element program with a subset of linear and nonlinear static and dynamic analysis capabilities with an integrated version of MSC/PATRAN for pre- and post-processing; and (4) MSC/DYTRAN - nonlinear explicit transient dynamics finite element program.

[Research Progress and Prospect of Applications of Finite Element Method in Lumbar Spine Biomechanics].

PubMed

Zhang, Zhenjun; Li, Yang; Liao, Zhenhua; Liu, Weiqiang

2016-12-01

Based on the application of finite element analysis in spine biomechanics,the research progress of finite element method applied in lumbar spine mechanics is reviewed and the prospect is forecasted.The related works,including lumbar ontology modeling,clinical application research,and occupational injury and protection,are summarized.The main research areas of finite element method are as follows:new accurate modeling process,the optimized simulation method,diversified clinical effect evaluation,and the clinical application of artificial lumbar disc.According to the recent research progress,the application prospects of finite element method,such as automation and individuation of modeling process,evaluation and analysis of new operation methods and simulation of mechanical damage and dynamic response,are discussed.The purpose of this paper is to provide the theoretical reference and practical guidance for the clinical lumbar problems by reviewing the application of finite element method in the field of the lumbar spine biomechanics.

Numerical modelling of bifurcation and localisation in cohesive-frictional materials

NASA Astrophysics Data System (ADS)

de Borst, René

1991-12-01

Methods are reviewed for analysing highly localised failure and bifurcation modes in discretised mechanical systems as typically arise in numerical simulations of failure in soils, rocks, metals and concrete. By the example of a plane-strain biaxial test it is shown that strain softening and lack of normality in elasto-plastic constitutive equations and the ensuing loss of ellipticity of the governing field equations cause a pathological mesh dependence of numerical solutions for such problems, thus rendering the results effectively meaningless. The need for introduction of higher-order continuum models is emphasised to remedy this shortcoming of the conventional approach. For one such a continuum model, namely the unconstrained Cosserat continuum, it is demonstrated that meaningful and convergent solutions (in the sense that a finite width of the localisation zone is computed upon mesh refinement) can be obtained.

Higher-order adaptive finite-element methods for Kohn–Sham density functional theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Motamarri, P.; Nowak, M.R.; Leiter, K.

2013-11-15

We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn–Sham density-functional theory (DFT). To this end, we develop an a priori mesh-adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss–Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100–200-fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposedmore » solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretizations of the Kohn–Sham DFT problem. Our studies suggest that staggering computational savings—of the order of 1000-fold—relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations, by using higher-order finite-element discretizations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy, suggesting that the hexic spectral-element may be an optimal choice for the finite-element discretization of the Kohn–Sham DFT problem. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations to within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688 atoms using modest computational resources, and good scalability of the present implementation up to 192 processors.« less

Comparison of Gap Elements and Contact Algorithm for 3D Contact Analysis of Spiral Bevel Gears

NASA Technical Reports Server (NTRS)

Bibel, G. D.; Tiku, K.; Kumar, A.; Handschuh, R.

1994-01-01

Three dimensional stress analysis of spiral bevel gears in mesh using the finite element method is presented. A finite element model is generated by solving equations that identify tooth surface coordinates. Contact is simulated by the automatic generation of nonpenetration constraints. This method is compared to a finite element contact analysis conducted with gap elements.

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

NASA Astrophysics Data System (ADS)

Yang, Zailin; Wang, Yao; Hei, Baoping

2013-12-01

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of onedimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.

A new conformal absorbing boundary condition for finite element meshes and parallelization of FEMATS

NASA Technical Reports Server (NTRS)

Chatterjee, A.; Volakis, J. L.; Nguyen, J.; Nurnberger, M.; Ross, D.

1993-01-01

Some of the progress toward the development and parallelization of an improved version of the finite element code FEMATS is described. This is a finite element code for computing the scattering by arbitrarily shaped three dimensional surfaces composite scatterers. The following tasks were worked on during the report period: (1) new absorbing boundary conditions (ABC's) for truncating the finite element mesh; (2) mixed mesh termination schemes; (3) hierarchical elements and multigridding; (4) parallelization; and (5) various modeling enhancements (antenna feeds, anisotropy, and higher order GIBC).

Books and monographs on finite element technology

NASA Technical Reports Server (NTRS)

Noor, A. K.

1985-01-01

The present paper proviees a listing of all of the English books and some of the foreign books on finite element technology, taking into account also a list of the conference proceedings devoted solely to finite elements. The references are divided into categories. Attention is given to fundamentals, mathematical foundations, structural and solid mechanics applications, fluid mechanics applications, other applied science and engineering applications, computer implementation and software systems, computational and modeling aspects, special topics, boundary element methods, proceedings of symmposia and conferences on finite element technology, bibliographies, handbooks, and historical accounts.

An evolution infinity Laplace equation modelling dynamic elasto-plastic torsion

NASA Astrophysics Data System (ADS)

Messelmi, Farid

2017-12-01

We consider in this paper a parabolic partial differential equation involving the infinity Laplace operator and a Leray-Lions operator with no coercitive assumption. We prove the existence and uniqueness of the corresponding approached problem and we show that at the limit the solution solves the parabolic variational inequality arising in the elasto-plastic torsion problem.

Reformulation of Nonlinear Anisotropic Crystal Elastoplasticity for Impact Physics

DTIC Science & Technology

2015-03-01

interest include metals, ceramics , minerals, and energetic materials . Accurate, efficient, stable, and thermodynamically consistent models for...Clayton JD. Phase field theory and analysis of pressure-shear induced amorphization and failure in boron carbide ceramic . AIMS Materials Science. 2014;1...of Nonlinear Anisotropic Crystal Elastoplasticity for Impact Physics by JD Clayton Weapons and Materials Research Directorate, ARL

An Automated Method for Landmark Identification and Finite-Element Modeling of the Lumbar Spine.

PubMed

Campbell, Julius Quinn; Petrella, Anthony J

2015-11-01

The purpose of this study was to develop a method for the automated creation of finite-element models of the lumbar spine. Custom scripts were written to extract bone landmarks of lumbar vertebrae and assemble L1-L5 finite-element models. End-plate borders, ligament attachment points, and facet surfaces were identified. Landmarks were identified to maintain mesh correspondence between meshes for later use in statistical shape modeling. 90 lumbar vertebrae were processed creating 18 subject-specific finite-element models. Finite-element model surfaces and ligament attachment points were reproduced within 1e-5 mm of the bone surface, including the critical contact surfaces of the facets. Element quality exceeded specifications in 97% of elements for the 18 models created. The current method is capable of producing subject-specific finite-element models of the lumbar spine with good accuracy, quality, and robustness. The automated methods developed represent advancement in the state of the art of subject-specific lumbar spine modeling to a scale not possible with prior manual and semiautomated methods.

Discrete maximum principle for the P1 - P0 weak Galerkin finite element approximations

NASA Astrophysics Data System (ADS)

Wang, Junping; Ye, Xiu; Zhai, Qilong; Zhang, Ran

2018-06-01

This paper presents two discrete maximum principles (DMP) for the numerical solution of second order elliptic equations arising from the weak Galerkin finite element method. The results are established by assuming an h-acute angle condition for the underlying finite element triangulations. The mathematical theory is based on the well-known De Giorgi technique adapted in the finite element context. Some numerical results are reported to validate the theory of DMP.

A Computational Approach for Automated Posturing of a Human Finite Element Model

DTIC Science & Technology

2016-07-01

Std. Z39.18 July 2016 Memorandum Report A Computational Approach for Automated Posturing of a Human Finite Element Model Justin McKee and Adam...protection by influencing the path that loading will be transferred into the body and is a major source of variability. The development of a finite element ...posture, human body, finite element , leg, spine 42 Adam Sokolow 410-306-2985Unclassified Unclassified Unclassified UU ii Approved for public release

Evaluation of Acoustic Propagation Paths into the Human Head

DTIC Science & Technology

2005-07-25

paths. A 3D finite-element solid mesh was constructed using a digital image database of an adult male head. Finite-element analysis was used to model the...air-borne sound pressure amplitude) via the alternate propagation paths. A 3D finite-element solid mesh was constructed using a digital image database ... database of an adult male head Coupled acoustic-mechanical finite-element analysis (FEA) was used to model the wave propagation through the fluid-solid

Semi-automatic sparse preconditioners for high-order finite element methods on non-uniform meshes

NASA Astrophysics Data System (ADS)

Austin, Travis M.; Brezina, Marian; Jamroz, Ben; Jhurani, Chetan; Manteuffel, Thomas A.; Ruge, John

2012-05-01

High-order finite elements often have a higher accuracy per degree of freedom than the classical low-order finite elements. However, in the context of implicit time-stepping methods, high-order finite elements present challenges to the construction of efficient simulations due to the high cost of inverting the denser finite element matrix. There are many cases where simulations are limited by the memory required to store the matrix and/or the algorithmic components of the linear solver. We are particularly interested in preconditioned Krylov methods for linear systems generated by discretization of elliptic partial differential equations with high-order finite elements. Using a preconditioner like Algebraic Multigrid can be costly in terms of memory due to the need to store matrix information at the various levels. We present a novel method for defining a preconditioner for systems generated by high-order finite elements that is based on a much sparser system than the original high-order finite element system. We investigate the performance for non-uniform meshes on a cube and a cubed sphere mesh, showing that the sparser preconditioner is more efficient and uses significantly less memory. Finally, we explore new methods to construct the sparse preconditioner and examine their effectiveness for non-uniform meshes. We compare results to a direct use of Algebraic Multigrid as a preconditioner and to a two-level additive Schwarz method.

Application of the Finite Element Method to Rotary Wing Aeroelasticity

NASA Technical Reports Server (NTRS)

Straub, F. K.; Friedmann, P. P.

1982-01-01

A finite element method for the spatial discretization of the dynamic equations of equilibrium governing rotary-wing aeroelastic problems is presented. Formulation of the finite element equations is based on weighted Galerkin residuals. This Galerkin finite element method reduces algebraic manipulative labor significantly, when compared to the application of the global Galerkin method in similar problems. The coupled flap-lag aeroelastic stability boundaries of hingeless helicopter rotor blades in hover are calculated. The linearized dynamic equations are reduced to the standard eigenvalue problem from which the aeroelastic stability boundaries are obtained. The convergence properties of the Galerkin finite element method are studied numerically by refining the discretization process. Results indicate that four or five elements suffice to capture the dynamics of the blade with the same accuracy as the global Galerkin method.

FINITE-ELEMENT ANALYSIS OF MULTIPHASE IMMISCIBLE FLOW THROUGH SOILS

EPA Science Inventory

A finite-element model is developed for multiphase flow through soil involving three immiscible fluids: namely, air, water, and a nonaqueous phase liquid (NAPL). A variational method is employed for the finite-element formulation corresponding to the coupled differential equation...

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Cooke, C. H.

1976-01-01

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Skeletal assessment with finite element analysis: relevance, pitfalls and interpretation.

PubMed

Campbell, Graeme Michael; Glüer, Claus-C

2017-07-01

Finite element models simulate the mechanical response of bone under load, enabling noninvasive assessment of strength. Models generated from quantitative computed tomography (QCT) incorporate the geometry and spatial distribution of bone mineral density (BMD) to simulate physiological and traumatic loads as well as orthopaedic implant behaviour. The present review discusses the current strengths and weakness of finite element models for application to skeletal biomechanics. In cadaver studies, finite element models provide better estimations of strength compared to BMD. Data from clinical studies are encouraging; however, the superiority of finite element models over BMD measures for fracture prediction has not been shown conclusively, and may be sex and site dependent. Therapeutic effects on bone strength are larger than for BMD; however, model validation has only been performed on untreated bone. High-resolution modalities and novel image processing methods may enhance the structural representation and predictive ability. Despite extensive use of finite element models to study orthopaedic implant stability, accurate simulation of the bone-implant interface and fracture progression remains a significant challenge. Skeletal finite element models provide noninvasive assessments of strength and implant stability. Improved structural representation and implant surface interaction may enable more accurate models of fragility in the future.

Efficient finite element simulation of slot spirals, slot radomes and microwave structures

NASA Technical Reports Server (NTRS)

Gong, J.; Volakis, J. L.

1995-01-01

This progress report contains the following two documents: (1) 'Efficient Finite Element Simulation of Slot Antennas using Prismatic Elements' - A hybrid finite element-boundary integral (FE-BI) simulation technique is discussed to treat narrow slot antennas etched on a planar platform. Specifically, the prismatic elements are used to reduce the redundant sampling rates and ease the mesh generation process. Numerical results for an antenna slot and frequency selective surfaces are presented to demonstrate the validity and capability of the technique; and (2) 'Application and Design Guidelines of the PML Absorber for Finite Element Simulations of Microwave Packages' - The recently introduced perfectly matched layer (PML) uniaxial absorber for frequency domain finite element simulations has several advantages. In this paper we present the application of PML for microwave circuit simulations along with design guidelines to obtain a desired level of absorption. Different feeding techniques are also investigated for improved accuracy.

Application of variational and Galerkin equations to linear and nonlinear finite element analysis

NASA Technical Reports Server (NTRS)

Yu, Y.-Y.

1974-01-01

The paper discusses the application of the variational equation to nonlinear finite element analysis. The problem of beam vibration with large deflection is considered. The variational equation is shown to be flexible in both the solution of a general problem and in the finite element formulation. Difficulties are shown to arise when Galerkin's equations are used in the consideration of the finite element formulation of two-dimensional linear elasticity and of the linear classical beam.

Plan, formulate, discuss and correlate a NASTRAN finite element vibrations model of the Boeing Model 360 helicopter airframe

NASA Technical Reports Server (NTRS)

Gabel, R.; Lang, P. F.; Smith, L. A.; Reed, D. A.

1989-01-01

Boeing Helicopter, together with other United States helicopter manufacturers, participated in a finite element applications program to emplace in the United States a superior capability to utilize finite element analysis models in support of helicopter airframe design. The activities relating to planning and creating a finite element vibrations model of the Boeing Model 36-0 composite airframe are summarized, along with the subsequent analytical correlation with ground shake test data.

Development and Application of the p-version of the Finite Element Method.

DTIC Science & Technology

1985-11-21

this property hierarchic families of finite elements. The h-version of the finite element method has been the subject of inten- sive study since the...early 1950’s and perhaps even earlier. Study of the p-version of the finite element method, on the other hand, began at Washington University in St...Louis in the early 1970’s and led to a more recent study of * .the h-p version. Research in the p-version (formerly called The Constraint Method) has

Self Diagnostic Adhesive for Bonded Joints in Aircraft Structures

DTIC Science & Technology

2016-10-04

validated under the fatigue/dynamic loading condition. 3) Both SEM (Spectral Element Modeling) and FEM ( Finite Element Modeling) simulation of the...Sensors ..................................................................... 22 Parametric Study of Sensor Performance via Finite Element Simulation...The frequency range that we are interested is around 800 kHz. Conventional linear finite element method (FEM) requires a very fine spatial

An elasto-plastic solution for channel cracking of brittle coating on polymer substrate

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Chao; Chen, Fangliang; Gray, Matthew H.

In this study, an elasto-plastic channel-cracking model is presented to study the open-mode fracture of a thin layer brittle coating grown on a polymer substrate. A linear elastic shear interlayer is introduced to describe the stress transfer from the elasto-plastic substrate to the brittle coating, on basis of the shear-lag principle. The channel cracking behavior involves three stages: elastic, elasto-plastic and plastic stages, which are solved in a continuous manner based on the deformation status of the substrate. Explicit solutions are derived for the mutli-stage cracking process. Corresponding experimental tests for a titanium oxide (TiO 2) coating on a polymore » (ethylene terephthalate) substrate are conducted. The fracture toughness of the coating layer is estimated based on the crack spacing versus layer thickness relationship at certain strain levels. This method is found to be more reliable than the traditional methods using crack onset strain. Parametric studies of the fracture energy release rate for the coating and interfacial compliance of the thin film system are conducted, through which the effect of plastic deformation on the channel cracking behavior is studied extensively. The results indicate that the tangent modulus of the substrate controls the evolution curvature of crack spacing where a smaller tangent modulus corresponds to a slower saturation of crack spacing. The energy release rate also varies significantly with the properties of the interlayer. The study highlights the necessity of an elasto-plastic model for the thin film systems of brittle coating on a plastic substrate.« less

An elasto-plastic solution for channel cracking of brittle coating on polymer substrate

DOE PAGES

Zhang, Chao; Chen, Fangliang; Gray, Matthew H.; ...

2017-04-25

In this study, an elasto-plastic channel-cracking model is presented to study the open-mode fracture of a thin layer brittle coating grown on a polymer substrate. A linear elastic shear interlayer is introduced to describe the stress transfer from the elasto-plastic substrate to the brittle coating, on basis of the shear-lag principle. The channel cracking behavior involves three stages: elastic, elasto-plastic and plastic stages, which are solved in a continuous manner based on the deformation status of the substrate. Explicit solutions are derived for the mutli-stage cracking process. Corresponding experimental tests for a titanium oxide (TiO 2) coating on a polymore » (ethylene terephthalate) substrate are conducted. The fracture toughness of the coating layer is estimated based on the crack spacing versus layer thickness relationship at certain strain levels. This method is found to be more reliable than the traditional methods using crack onset strain. Parametric studies of the fracture energy release rate for the coating and interfacial compliance of the thin film system are conducted, through which the effect of plastic deformation on the channel cracking behavior is studied extensively. The results indicate that the tangent modulus of the substrate controls the evolution curvature of crack spacing where a smaller tangent modulus corresponds to a slower saturation of crack spacing. The energy release rate also varies significantly with the properties of the interlayer. The study highlights the necessity of an elasto-plastic model for the thin film systems of brittle coating on a plastic substrate.« less

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 2: Derivation of finite-element equations and comparisons with analytical solutions

USGS Publications Warehouse

Cooley, Richard L.

1992-01-01

MODFE, a modular finite-element model for simulating steady- or unsteady-state, area1 or axisymmetric flow of ground water in a heterogeneous anisotropic aquifer is documented in a three-part series of reports. In this report, part 2, the finite-element equations are derived by minimizing a functional of the difference between the true and approximate hydraulic head, which produces equations that are equivalent to those obtained by either classical variational or Galerkin techniques. Spatial finite elements are triangular with linear basis functions, and temporal finite elements are one dimensional with linear basis functions. Physical processes that can be represented by the model include (1) confined flow, unconfined flow (using the Dupuit approximation), or a combination of both; (2) leakage through either rigid or elastic confining units; (3) specified recharge or discharge at points, along lines, or areally; (4) flow across specified-flow, specified-head, or head-dependent boundaries; (5) decrease of aquifer thickness to zero under extreme water-table decline and increase of aquifer thickness from zero as the water table rises; and (6) head-dependent fluxes from springs, drainage wells, leakage across riverbeds or confining units combined with aquifer dewatering, and evapotranspiration. The matrix equations produced by the finite-element method are solved by the direct symmetric-Doolittle method or the iterative modified incomplete-Cholesky conjugate-gradient method. The direct method can be efficient for small- to medium-sized problems (less than about 500 nodes), and the iterative method is generally more efficient for larger-sized problems. Comparison of finite-element solutions with analytical solutions for five example problems demonstrates that the finite-element model can yield accurate solutions to ground-water flow problems.

Error analysis and correction of discrete solutions from finite element codes

NASA Technical Reports Server (NTRS)

Thurston, G. A.; Stein, P. A.; Knight, N. F., Jr.; Reissner, J. E.

1984-01-01

Many structures are an assembly of individual shell components. Therefore, results for stresses and deflections from finite element solutions for each shell component should agree with the equations of shell theory. This paper examines the problem of applying shell theory to the error analysis and the correction of finite element results. The general approach to error analysis and correction is discussed first. Relaxation methods are suggested as one approach to correcting finite element results for all or parts of shell structures. Next, the problem of error analysis of plate structures is examined in more detail. The method of successive approximations is adapted to take discrete finite element solutions and to generate continuous approximate solutions for postbuckled plates. Preliminary numerical results are included.

Scalable Implementation of Finite Elements by NASA _ Implicit (ScIFEi)

NASA Technical Reports Server (NTRS)

Warner, James E.; Bomarito, Geoffrey F.; Heber, Gerd; Hochhalter, Jacob D.

2016-01-01

Scalable Implementation of Finite Elements by NASA (ScIFEN) is a parallel finite element analysis code written in C++. ScIFEN is designed to provide scalable solutions to computational mechanics problems. It supports a variety of finite element types, nonlinear material models, and boundary conditions. This report provides an overview of ScIFEi (\\Sci-Fi"), the implicit solid mechanics driver within ScIFEN. A description of ScIFEi's capabilities is provided, including an overview of the tools and features that accompany the software as well as a description of the input and output le formats. Results from several problems are included, demonstrating the efficiency and scalability of ScIFEi by comparing to finite element analysis using a commercial code.

Optimal least-squares finite element method for elliptic problems

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Povinelli, Louis A.

1991-01-01

An optimal least squares finite element method is proposed for two dimensional and three dimensional elliptic problems and its advantages are discussed over the mixed Galerkin method and the usual least squares finite element method. In the usual least squares finite element method, the second order equation (-Delta x (Delta u) + u = f) is recast as a first order system (-Delta x p + u = f, Delta u - p = 0). The error analysis and numerical experiment show that, in this usual least squares finite element method, the rate of convergence for flux p is one order lower than optimal. In order to get an optimal least squares method, the irrotationality Delta x p = 0 should be included in the first order system.

The effectiveness of element downsizing on a three-dimensional finite element model of bone trabeculae in implant biomechanics.

PubMed

Sato, Y; Wadamoto, M; Tsuga, K; Teixeira, E R

1999-04-01

More validity of finite element analysis in implant biomechanics requires element downsizing. However, excess downsizing needs computer memory and calculation time. To investigate the effectiveness of element downsizing on the construction of a three-dimensional finite element bone trabeculae model, with different element sizes (600, 300, 150 and 75 microm) models were constructed and stress induced by vertical 10 N loading was analysed. The difference in von Mises stress values between the models with 600 and 300 microm element sizes was larger than that between 300 and 150 microm. On the other hand, no clear difference of stress values was detected among the models with 300, 150 and 75 microm element sizes. Downsizing of elements from 600 to 300 microm is suggested to be effective in the construction of a three-dimensional finite element bone trabeculae model for possible saving of computer memory and calculation time in the laboratory.

Elasto-plastic analysis of interface layers for fiber reinforced metal matrix composites

NASA Technical Reports Server (NTRS)

Doghri, I.; Leckie, F. A.

1991-01-01

The mismatch in coefficients of thermal expansion (CTE) of fiber and matrix in metal matrix composites reinforced with ceramic fibers induces high thermal stresses in the matrix. Elasto-plastic analyses - with different degrees of simplification and modelization - show that an interface layer with a sufficiently high CTE can reduce the tensile hoop stress in the matrix substantially.

Using a multifrontal sparse solver in a high performance, finite element code

NASA Technical Reports Server (NTRS)

King, Scott D.; Lucas, Robert; Raefsky, Arthur

1990-01-01

We consider the performance of the finite element method on a vector supercomputer. The computationally intensive parts of the finite element method are typically the individual element forms and the solution of the global stiffness matrix both of which are vectorized in high performance codes. To further increase throughput, new algorithms are needed. We compare a multifrontal sparse solver to a traditional skyline solver in a finite element code on a vector supercomputer. The multifrontal solver uses the Multiple-Minimum Degree reordering heuristic to reduce the number of operations required to factor a sparse matrix and full matrix computational kernels (e.g., BLAS3) to enhance vector performance. The net result in an order-of-magnitude reduction in run time for a finite element application on one processor of a Cray X-MP.

Structural Acoustic Physics Based Modeling of Curved Composite Shells

DTIC Science & Technology

2017-09-19

Results show that the finite element computational models accurately match analytical calculations, and that the composite material studied in this...products. 15. SUBJECT TERMS Finite Element Analysis, Structural Acoustics, Fiber-Reinforced Composites, Physics-Based Modeling 16. SECURITY...2 4 FINITE ELEMENT MODEL DESCRIPTION

Finite element analysis of thrust angle contact ball slewing bearing

NASA Astrophysics Data System (ADS)

Deng, Biao; Guo, Yuan; Zhang, An; Tang, Shengjin

2017-12-01

In view of the large heavy slewing bearing no longer follows the rigid ring hupothesis under the load condition, the entity finite element model of thrust angular contact ball bearing was established by using finite element analysis software ANSYS. The boundary conditions of the model were set according to the actual condition of slewing bearing, the internal stress state of the slewing bearing was obtained by solving and calculation, and the calculated results were compared with the numerical results based on the rigid ring assumption. The results show that more balls are loaded in the result of finite element method, and the maximum contact stresses between the ball and raceway have some reductions. This is because the finite element method considers the ferrule as an elastic body. The ring will produce structure deformation in the radial plane when the heavy load slewing bearings are subjected to external loads. The results of the finite element method are more in line with the actual situation of the slewing bearing in the engineering.

Challenges in Integrating Nondestructive Evaluation and Finite Element Methods for Realistic Structural Analysis

NASA Technical Reports Server (NTRS)

Abdul-Aziz, Ali; Baaklini, George Y.; Zagidulin, Dmitri; Rauser, Richard W.

2000-01-01

Capabilities and expertise related to the development of links between nondestructive evaluation (NDE) and finite element analysis (FEA) at Glenn Research Center (GRC) are demonstrated. Current tools to analyze data produced by computed tomography (CT) scans are exercised to help assess the damage state in high temperature structural composite materials. A utility translator was written to convert velocity (an image processing software) STL data file to a suitable CAD-FEA type file. Finite element analyses are carried out with MARC, a commercial nonlinear finite element code, and the analytical results are discussed. Modeling was established by building MSC/Patran (a pre and post processing finite element package) generated model and comparing it to a model generated by Velocity in conjunction with MSC/Patran Graphics. Modeling issues and results are discussed in this paper. The entire process that outlines the tie between the data extracted via NDE and the finite element modeling and analysis is fully described.

An enriched finite element method to fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam

2017-08-01

In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.

On the existence and stability conditions for mixed-hybrid finite element solutions based on Reissner's variational principle

NASA Technical Reports Server (NTRS)

Karlovitz, L. A.; Atluri, S. N.; Xue, W.-M.

1985-01-01

The extensions of Reissner's two-field (stress and displacement) principle to the cases wherein the displacement field is discontinuous and/or the stress field results in unreciprocated tractions, at a finite number of surfaces ('interelement boundaries') in a domain (as, for instance, when the domain is discretized into finite elements), is considered. The conditions for the existence, uniqueness, and stability of mixed-hybrid finite element solutions based on such discontinuous fields, are summarized. The reduction of these global conditions to local ('element') level, and the attendant conditions on the ranks of element matrices, are discussed. Two examples of stable, invariant, least-order elements - a four-node square planar element and an eight-node cubic element - are discussed in detail.

GENSURF: A mesh generator for 3D finite element analysis of surface and corner cracks in finite thickness plates subjected to mode-1 loadings

NASA Technical Reports Server (NTRS)

Raju, I. S.

1992-01-01

A computer program that generates three-dimensional (3D) finite element models for cracked 3D solids was written. This computer program, gensurf, uses minimal input data to generate 3D finite element models for isotropic solids with elliptic or part-elliptic cracks. These models can be used with a 3D finite element program called surf3d. This report documents this mesh generator. In this manual the capabilities, limitations, and organization of gensurf are described. The procedures used to develop 3D finite element models and the input for and the output of gensurf are explained. Several examples are included to illustrate the use of this program. Several input data files are included with this manual so that the users can edit these files to conform to their crack configuration and use them with gensurf.

A comparative study of finite element and finite difference methods for Cauchy-Riemann type equations

NASA Technical Reports Server (NTRS)

Fix, G. J.; Rose, M. E.

1983-01-01

A least squares formulation of the system divu = rho, curlu = zeta is surveyed from the viewpoint of both finite element and finite difference methods. Closely related arguments are shown to establish convergence estimates.

Robustness analysis of elastoplastic structure subjected to double impulse

NASA Astrophysics Data System (ADS)

Kanno, Yoshihiro; Takewaki, Izuru

2016-11-01

The double impulse has extensively been used to evaluate the critical response of an elastoplastic structure against a pulse-type input, including near-fault earthquake ground motions. In this paper, we propose a robustness assessment method for elastoplastic single-degree-of-freedom structures subjected to the double impulse input. Uncertainties in the initial velocity of the input, as well as the natural frequency and the strength of the structure, are considered. As fundamental properties of the structural robustness, we show monotonicity of the robustness measure with respect to the natural frequency. In contrast, we show that robustness is not necessarily improved even if the structural strength is increased. Moreover, the robustness preference between two structures with different values of structural strength can possibly reverse when the performance requirement is changed.

CFD Analysis of the SBXC Glider Airframe

DTIC Science & Technology

2016-06-01

mathematically on finite element methods. To validate and verify the methodology developed, a mathematical comparison was made with the previous research data...greater than 15 m/s. 14. SUBJECT TERMS finite element method, computational fluid dynamics, Y Plus, mesh element quality, aerodynamic data, fluid...based mathematically on finite element methods. To validate and verify the methodology developed, a mathematical comparison was made with the

Functional Data Approximation on Bounded Domains using Polygonal Finite Elements.

PubMed

Cao, Juan; Xiao, Yanyang; Chen, Zhonggui; Wang, Wenping; Bajaj, Chandrajit

2018-07-01

We construct and analyze piecewise approximations of functional data on arbitrary 2D bounded domains using generalized barycentric finite elements, and particularly quadratic serendipity elements for planar polygons. We compare approximation qualities (precision/convergence) of these partition-of-unity finite elements through numerical experiments, using Wachspress coordinates, natural neighbor coordinates, Poisson coordinates, mean value coordinates, and quadratic serendipity bases over polygonal meshes on the domain. For a convex n -sided polygon, the quadratic serendipity elements have 2 n basis functions, associated in a Lagrange-like fashion to each vertex and each edge midpoint, rather than the usual n ( n + 1)/2 basis functions to achieve quadratic convergence. Two greedy algorithms are proposed to generate Voronoi meshes for adaptive functional/scattered data approximations. Experimental results show space/accuracy advantages for these quadratic serendipity finite elements on polygonal domains versus traditional finite elements over simplicial meshes. Polygonal meshes and parameter coefficients of the quadratic serendipity finite elements obtained by our greedy algorithms can be further refined using an L 2 -optimization to improve the piecewise functional approximation. We conduct several experiments to demonstrate the efficacy of our algorithm for modeling features/discontinuities in functional data/image approximation.

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 2: User's manual

NASA Technical Reports Server (NTRS)

Nakazawa, Shohei

1989-01-01

The user options available for running the MHOST finite element analysis package is described. MHOST is a solid and structural analysis program based on the mixed finite element technology, and is specifically designed for 3-D inelastic analysis. A family of 2- and 3-D continuum elements along with beam and shell structural elements can be utilized, many options are available in the constitutive equation library, the solution algorithms and the analysis capabilities. The outline of solution algorithms is discussed along with the data input and output, analysis options including the user subroutines and the definition of the finite elements implemented in the program package.

Influence of Finite Element Software on Energy Release Rates Computed Using the Virtual Crack Closure Technique

NASA Technical Reports Server (NTRS)

Krueger, Ronald; Goetze, Dirk; Ransom, Jonathon (Technical Monitor)

2006-01-01

Strain energy release rates were computed along straight delamination fronts of Double Cantilever Beam, End-Notched Flexure and Single Leg Bending specimens using the Virtual Crack Closure Technique (VCCT). Th e results were based on finite element analyses using ABAQUS# and ANSYS# and were calculated from the finite element results using the same post-processing routine to assure a consistent procedure. Mixed-mode strain energy release rates obtained from post-processing finite elem ent results were in good agreement for all element types used and all specimens modeled. Compared to previous studies, the models made of s olid twenty-node hexahedral elements and solid eight-node incompatible mode elements yielded excellent results. For both codes, models made of standard brick elements and elements with reduced integration did not correctly capture the distribution of the energy release rate acr oss the width of the specimens for the models chosen. The results suggested that element types with similar formulation yield matching results independent of the finite element software used. For comparison, m ixed-mode strain energy release rates were also calculated within ABAQUS#/Standard using the VCCT for ABAQUS# add on. For all specimens mod eled, mixed-mode strain energy release rates obtained from ABAQUS# finite element results using post-processing were almost identical to re sults calculated using the VCCT for ABAQUS# add on.

Advanced development of BEM for elastic and inelastic dynamic analysis of solids

NASA Technical Reports Server (NTRS)

Banerjee, P. K.; Ahmad, S.; Wang, H. C.

1989-01-01

Direct Boundary Element formulations and their numerical implementation for periodic and transient elastic as well as inelastic transient dynamic analyses of two-dimensional, axisymmetric and three-dimensional solids are presented. The inelastic formulation is based on an initial stress approach and is the first of its kind in the field of Boundary Element Methods. This formulation employs the Navier-Cauchy equation of motion, Graffi's dynamic reciprocal theorem, Stokes' fundamental solution, and the divergence theorem, together with kinematical and constitutive equations to obtain the pertinent integral equations of the problem in the time domain within the context of the small displacement theory of elastoplasticity. The dynamic (periodic, transient as well as nonlinear transient) formulations have been applied to a range of problems. The numerical formulations presented here are included in the BEST3D and GPBEST systems.

Application of finite element approach to transonic flow problems

NASA Technical Reports Server (NTRS)

Hafez, M. M.; Murman, E. M.; Wellford, L. C., Jr.

1976-01-01

A variational finite element model for transonic small disturbance calculations is described. Different strategy is adopted in subsonic and supersonic regions, and blending elements are introduced between different regions. In the supersonic region, no upstream effect is allowed. If rectangular elements with linear shape functions are used, the model is similar to Murman's finite difference operators. Higher order shape functions, nonrectangular elements, and discontinuous approximation of shock waves are also discussed.

Finite element meshing of ANSYS (trademark) solid models

NASA Technical Reports Server (NTRS)

Kelley, F. S.

1987-01-01

A large scale, general purpose finite element computer program, ANSYS, developed and marketed by Swanson Analysis Systems, Inc. is discussed. ANSYS was perhaps the first commercially available program to offer truly interactive finite element model generation. ANSYS's purpose is for solid modeling. This application is briefly discussed and illustrated.

Nonlinear finite element modeling of corrugated board

Treesearch

A. C. Gilchrist; J. C. Suhling; T. J. Urbanik

1999-01-01

In this research, an investigation on the mechanical behavior of corrugated board has been performed using finite element analysis. Numerical finite element models for corrugated board geometries have been created and executed. Both geometric (large deformation) and material nonlinearities were included in the models. The analyses were performed using the commercial...

A Version of the Graphics-Oriented Interactive Finite Element Time-Sharing System (GIFTS) for an IBM with CP/CMS.

DTIC Science & Technology

1982-03-01

POSTGRADUATE SCHOOL fMonterey, California THESIS A VERSION OF THE GRAPHICS-ORIENTED INTERACTIVE FINITE ELEMENT TIME-SHARING SYSTEM ( GIFTS ) FOR AN IBM...Master’s & Engineer’s active Finite Element Time-sharing System Thesis - March 1982 ( GIFTS ) for an IBM with CP/CMS 6. penromm.oOn. REPoRT MUlmiR 1. AUTHOIee...ss0in D dinuf 5W M memisi) ’A version of the Graphics-oriented, Interactive, Finite element, Time-sharing System ( GIFTS ) has been developed for, and

An Error Analysis for the Finite Element Method Applied to Convection Diffusion Problems.

DTIC Science & Technology

1981-03-01

D TFhG-]NOLOGY k 4b 00 \\" ) ’b Technical Note BN-962 AN ERROR ANALYSIS FOR THE FINITE ELEMENT METHOD APPLIED TO CONVECTION DIFFUSION PROBLEM by I...Babu~ka and W. G. Szym’czak March 1981 V.. UNVI I Of- ’i -S AN ERROR ANALYSIS FOR THE FINITE ELEMENT METHOD P. - 0 w APPLIED TO CONVECTION DIFFUSION ...AOAO98 895 MARYLAND UNIVYCOLLEGE PARK INST FOR PHYSICAL SCIENCE--ETC F/G 12/I AN ERROR ANALYIS FOR THE FINITE ELEMENT METHOD APPLIED TO CONV..ETC (U

Quality assessment and control of finite element solutions

NASA Technical Reports Server (NTRS)

Noor, Ahmed K.; Babuska, Ivo

1987-01-01

Status and some recent developments in the techniques for assessing the reliability of finite element solutions are summarized. Discussion focuses on a number of aspects including: the major types of errors in the finite element solutions; techniques used for a posteriori error estimation and the reliability of these estimators; the feedback and adaptive strategies for improving the finite element solutions; and postprocessing approaches used for improving the accuracy of stresses and other important engineering data. Also, future directions for research needed to make error estimation and adaptive movement practical are identified.

Use of system identification techniques for improving airframe finite element models using test data

NASA Technical Reports Server (NTRS)

Hanagud, Sathya V.; Zhou, Weiyu; Craig, James I.; Weston, Neil J.

1991-01-01

A method for using system identification techniques to improve airframe finite element models was developed and demonstrated. The method uses linear sensitivity matrices to relate changes in selected physical parameters to changes in total system matrices. The values for these physical parameters were determined using constrained optimization with singular value decomposition. The method was confirmed using both simple and complex finite element models for which pseudo-experimental data was synthesized directly from the finite element model. The method was then applied to a real airframe model which incorporated all the complexities and details of a large finite element model and for which extensive test data was available. The method was shown to work, and the differences between the identified model and the measured results were considered satisfactory.

Weak form implementation of the semi-analytical finite element (SAFE) method for a variety of elastodynamic waveguides

NASA Astrophysics Data System (ADS)

Hakoda, Christopher; Lissenden, Clifford; Rose, Joseph L.

2018-04-01

Dispersion curves are essential to any guided wave NDE project. The Semi-Analytical Finite Element (SAFE) method has significantly increased the ease by which these curves can be calculated. However, due to misconceptions regarding theory and fragmentation based on different finite-element software, the theory has stagnated, and adoption by researchers who are new to the field has been slow. This paper focuses on the relationship between the SAFE formulation and finite element theory, and the implementation of the SAFE method in a weak form for plates, pipes, layered waveguides/composites, curved waveguides, and arbitrary cross-sections is shown. The benefits of the weak form are briefly described, as is implementation in open-source and commercial finite element software.

Nonlinear Legendre Spectral Finite Elements for Wind Turbine Blade Dynamics: Preprint

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Q.; Sprague, M. A.; Jonkman, J.

2014-01-01

This paper presents a numerical implementation and examination of new wind turbine blade finite element model based on Geometrically Exact Beam Theory (GEBT) and a high-order spectral finite element method. The displacement-based GEBT is presented, which includes the coupling effects that exist in composite structures and geometric nonlinearity. Legendre spectral finite elements (LSFEs) are high-order finite elements with nodes located at the Gauss-Legendre-Lobatto points. LSFEs can be an order of magnitude more efficient that low-order finite elements for a given accuracy level. Interpolation of the three-dimensional rotation, a major technical barrier in large-deformation simulation, is discussed in the context ofmore » LSFEs. It is shown, by numerical example, that the high-order LSFEs, where weak forms are evaluated with nodal quadrature, do not suffer from a drawback that exists in low-order finite elements where the tangent-stiffness matrix is calculated at the Gauss points. Finally, the new LSFE code is implemented in the new FAST Modularization Framework for dynamic simulation of highly flexible composite-material wind turbine blades. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples showing validation and LSFE performance will be provided in the final paper.« less

[Progression on finite element modeling method in scoliosis].

PubMed

Fan, Ning; Zang, Lei; Hai, Yong; Du, Peng; Yuan, Shuo

2018-04-25

Scoliosis is a complex spinal three-dimensional malformation with complicated pathogenesis, often associated with complications as thoracic deformity and shoulder imbalance. Because the acquisition of specimen or animal models are difficult, the biomechanical study of scoliosis is limited. In recent years, along with the development of the computer technology, software and image, the technology of establishing a finite element model of human spine is maturing and it has been providing strong support for the research of pathogenesis of scoliosis, the design and application of brace, and the selection of surgical methods. The finite element model method is gradually becoming an important tool in the biomechanical study of scoliosis. Establishing a high quality finite element model is the basis of analysis and future study. However, the finite element modeling process can be complex and modeling methods are greatly varied. Choosing the appropriate modeling method according to research objectives has become researchers' primary task. In this paper, the author reviews the national and international literature in recent years and concludes the finite element modeling methods in scoliosis, including data acquisition, establishment of the geometric model, the material properties, parameters setting, the validity of the finite element model validation and so on. Copyright© 2018 by the China Journal of Orthopaedics and Traumatology Press.

Finite elements of nonlinear continua.

NASA Technical Reports Server (NTRS)

Oden, J. T.

1972-01-01

The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Finite element solution for energy conservation using a highly stable explicit integration algorithm

NASA Technical Reports Server (NTRS)

Baker, A. J.; Manhardt, P. D.

1972-01-01

Theoretical derivation of a finite element solution algorithm for the transient energy conservation equation in multidimensional, stationary multi-media continua with irregular solution domain closure is considered. The complete finite element matrix forms for arbitrarily irregular discretizations are established, using natural coordinate function representations. The algorithm is embodied into a user-oriented computer program (COMOC) which obtains transient temperature distributions at the node points of the finite element discretization using a highly stable explicit integration procedure with automatic error control features. The finite element algorithm is shown to posses convergence with discretization for a transient sample problem. The condensed form for the specific heat element matrix is shown to be preferable to the consistent form. Computed results for diverse problems illustrate the versatility of COMOC, and easily prepared output subroutines are shown to allow quick engineering assessment of solution behavior.

Periodic trim solutions with hp-version finite elements in time

NASA Technical Reports Server (NTRS)

Peters, David A.; Hou, Lin-Jun

1990-01-01

Finite elements in time as an alternative strategy for rotorcraft trim problems are studied. The research treats linear flap and linearized flap-lag response both for quasi-trim and trim cases. The connection between Fourier series analysis and hp-finite elements for periodic a problem is also examined. It is proved that Fourier series is a special case of space-time finite elements in which one element is used with a strong displacement formulation. Comparisons are made with respect to accuracy among Fourier analysis, displacement methods, and mixed methods over a variety parameters. The hp trade-off is studied for the periodic trim problem to provide an optimum step size and order of polynomial for a given error criteria. It is found that finite elements in time can outperform Fourier analysis for periodic problems, and for some given error criteria. The mixed method provides better results than does the displacement method.

Finite strain anisotropic elasto-plastic model for the simulation of the forming and testing of metal/short fiber reinforced polymer clinch joints at room temperature

NASA Astrophysics Data System (ADS)

Dean, A.; Rolfes, R.; Behrens, A.; Bouguecha, A.; Hübner, S.; Bonk, C.; Grbic, N.

2017-10-01

There is a strong trend in the automotive industry to reduce car body-, chassis- and power-train mass in order to lower carbon emissions. More wide spread use of lightweight short fiber reinforced polymer (SFRP) is a promising approach to attain this goal. This poses the challenge of how to integrate new SFRP components by joining them to traditional sheet metal structures. Recently (1), the clinching technique has been successfully applied as a suitable joining method for dissimilar material such as SFRP and Aluminum. The material pairing PA6GF30 and EN AW 5754 is chosen for this purpose due to their common application in industry. The current contribution presents a verification and validation of a finite strain anisotropic material model for SFRP developed in (2) for the FE simulation of the hybrid clinching process. The finite fiber rotation during forming and separation, and thus the change of the preferential material direction, is represented in this model. Plastic deformations in SFRP are considered in this model via an invariant based non-associated plasticity formulation following the multiplicative decomposition approach of the deformation gradient where the stress-free intermediate configuration is introduced. The model allows for six independent characterization curves. The aforementioned material model allows for a detailed simulation of the forming process as well as a simulative prediction of the shear test strength of the produced joint at room temperature.

Contact Stress Analysis of Spiral Bevel Gears Using Finite Element Analysis

NASA Technical Reports Server (NTRS)

Bibel, G. D.; Kumar, A; Reddy, S.; Handschuh, R.

1995-01-01

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Nonlinear solid finite element analysis of mitral valves with heterogeneous leaflet layers

NASA Astrophysics Data System (ADS)

Prot, V.; Skallerud, B.

2009-02-01

An incompressible transversely isotropic hyperelastic material for solid finite element analysis of a porcine mitral valve response is described. The material model implementation is checked in single element tests and compared with a membrane implementation in an out-of-plane loading test to study how the layered structures modify the stress response for a simple geometry. Three different collagen layer arrangements are used in finite element analysis of the mitral valve. When the leaflets are arranged in two layers with the collagen on the ventricular side, the stress in the fibre direction through the thickness in the central part of the anterior leaflet is homogenized and the peak stress is reduced. A simulation using membrane elements is also carried out for comparison with the solid finite element results. Compared to echocardiographic measurements, the finite element models bulge too much in the left atrium. This may be due to evidence of active muscle fibres in some parts of the anterior leaflet, whereas our constitutive modelling is based on passive material.

Coupled porohyperelastic mass transport (PHEXPT) finite element models for soft tissues using ABAQUS.

PubMed

Vande Geest, Jonathan P; Simon, B R; Rigby, Paul H; Newberg, Tyler P

2011-04-01

Finite element models (FEMs) including characteristic large deformations in highly nonlinear materials (hyperelasticity and coupled diffusive/convective transport of neutral mobile species) will allow quantitative study of in vivo tissues. Such FEMs will provide basic understanding of normal and pathological tissue responses and lead to optimization of local drug delivery strategies. We present a coupled porohyperelastic mass transport (PHEXPT) finite element approach developed using a commercially available ABAQUS finite element software. The PHEXPT transient simulations are based on sequential solution of the porohyperelastic (PHE) and mass transport (XPT) problems where an Eulerian PHE FEM is coupled to a Lagrangian XPT FEM using a custom-written FORTRAN program. The PHEXPT theoretical background is derived in the context of porous media transport theory and extended to ABAQUS finite element formulations. The essential assumptions needed in order to use ABAQUS are clearly identified in the derivation. Representative benchmark finite element simulations are provided along with analytical solutions (when appropriate). These simulations demonstrate the differences in transient and steady state responses including finite deformations, total stress, fluid pressure, relative fluid, and mobile species flux. A detailed description of important model considerations (e.g., material property functions and jump discontinuities at material interfaces) is also presented in the context of finite deformations. The ABAQUS-based PHEXPT approach enables the use of the available ABAQUS capabilities (interactive FEM mesh generation, finite element libraries, nonlinear material laws, pre- and postprocessing, etc.). PHEXPT FEMs can be used to simulate the transport of a relatively large neutral species (negligible osmotic fluid flux) in highly deformable hydrated soft tissues and tissue-engineered materials.

Electromagnetic finite elements based on a four-potential variational principle

NASA Technical Reports Server (NTRS)

Schuler, James J.; Felippa, Carlos A.

1991-01-01

Electromagnetic finite elements based on a variational principle that uses the electromagnetic four-potential as a primary variable are derived. This choice is used to construct elements suitable for downstream coupling with mechanical and thermal finite elements for the analysis of electromagnetic/mechanical systems that involve superconductors. The main advantages of the four-potential as a basis for finite element formulation are that the number of degrees of freedom per node remains modest as the problem dimensionally increases, that jump discontinuities on interfaces are naturally accommodated, and that statics as well as dynamics may be treated without any a priori approximations. The new elements are tested on an axisymmetric problem under steady state forcing conditions. The results are in excellent agreement with analytical solutions.

Application of the control volume mixed finite element method to a triangular discretization

USGS Publications Warehouse

Naff, R.L.

2012-01-01

A two-dimensional control volume mixed finite element method is applied to the elliptic equation. Discretization of the computational domain is based in triangular elements. Shape functions and test functions are formulated on the basis of an equilateral reference triangle with unit edges. A pressure support based on the linear interpolation of elemental edge pressures is used in this formulation. Comparisons are made between results from the standard mixed finite element method and this control volume mixed finite element method. Published 2011. This article is a US Government work and is in the public domain in the USA. ?? 2012 John Wiley & Sons, Ltd. This article is a US Government work and is in the public domain in the USA.

Second order tensor finite element

NASA Technical Reports Server (NTRS)

Oden, J. Tinsley; Fly, J.; Berry, C.; Tworzydlo, W.; Vadaketh, S.; Bass, J.

1990-01-01

The results of a research and software development effort are presented for the finite element modeling of the static and dynamic behavior of anisotropic materials, with emphasis on single crystal alloys. Various versions of two dimensional and three dimensional hybrid finite elements were implemented and compared with displacement-based elements. Both static and dynamic cases are considered. The hybrid elements developed in the project were incorporated into the SPAR finite element code. In an extension of the first phase of the project, optimization of experimental tests for anisotropic materials was addressed. In particular, the problem of calculating material properties from tensile tests and of calculating stresses from strain measurements were considered. For both cases, numerical procedures and software for the optimization of strain gauge and material axes orientation were developed.

Adaptive finite element methods for two-dimensional problems in computational fracture mechanics

NASA Technical Reports Server (NTRS)

Min, J. B.; Bass, J. M.; Spradley, L. W.

1994-01-01

Some recent results obtained using solution-adaptive finite element methods in two-dimensional problems in linear elastic fracture mechanics are presented. The focus is on the basic issue of adaptive finite element methods for validating the new methodology by computing demonstration problems and comparing the stress intensity factors to analytical results.

N%-Superconvergence of Finite Element Approximations in the Interior of General Meshes of Triangles

DTIC Science & Technology

1993-12-01

RODiGuEz, On the asymptotic exactness of error estimators for linear triangular finite elements, Numer. Math., 59 (1991), pp. 107-127. 27. R. DURAN ...WAHLDIN, Interior maxmum norma estimates for finite element methods, Part H, unpublished manuscript. 38. I. BABUfKA, T. STROUBOULIS, A. MATHU. AND C.S

75 FR 70623 - Airworthiness Directives; DORNIER LUFTFAHRT GmbH Models Dornier 228-100, Dornier 228-101, Dornier...

Federal Register 2010, 2011, 2012, 2013, 2014

2010-11-18

... measurements as well as finite element modelling and fatigue analyses to better understand the stress... include strain measurements as well as finite element modeling and fatigue analyses to better understand... finite element modelling and fatigue analyses to better understand the stress distribution onto the frame...

Evaluation of an improved finite-element thermal stress calculation technique

NASA Technical Reports Server (NTRS)

Camarda, C. J.

1982-01-01

A procedure for generating accurate thermal stresses with coarse finite element grids (Ojalvo's method) is described. The procedure is based on the observation that for linear thermoelastic problems, the thermal stresses may be envisioned as being composed of two contributions; the first due to the strains in the structure which depend on the integral of the temperature distribution over the finite element and the second due to the local variation of the temperature in the element. The first contribution can be accurately predicted with a coarse finite-element mesh. The resulting strain distribution can then be combined via the constitutive relations with detailed temperatures from a separate thermal analysis. The result is accurate thermal stresses from coarse finite element structural models even where the temperature distributions have sharp variations. The range of applicability of the method for various classes of thermostructural problems such as in-plane or bending type problems and the effect of the nature of the temperature distribution and edge constraints are addressed. Ojalvo's method is used in conjunction with the SPAR finite element program. Results are obtained for rods, membranes, a box beam and a stiffened panel.

Deformation analysis of rotary combustion engine housings

NASA Technical Reports Server (NTRS)

Vilmann, Carl

1991-01-01

This analysis of the deformation of rotary combustion engine housings targeted the following objectives: (1) the development and verification of a finite element model of the trochoid housing, (2) the prediction of the stress and deformation fields present within the trochoid housing during operating conditions, and (3) the development of a specialized preprocessor which would shorten the time necessary for mesh generation of a trochoid housing's FEM model from roughly one month to approximately two man hours. Executable finite element models were developed for both the Mazda and the Outboard Marine Corporation trochoid housings. It was also demonstrated that a preprocessor which would hasten the generation of finite element models of a rotary engine was possible to develop. The above objectives are treated in detail in the attached appendices. The first deals with finite element modeling of a Wankel engine center housing, and the second with the development of a preprocessor that generates finite element models of rotary combustion engine center housings. A computer program, designed to generate finite element models of user defined rotary combustion engine center housing geometries, is also included.

Finite-element simulation of blood perfusion in muscle tissue during compression and sustained contraction.

PubMed

Vankan, W J; Huyghe, J M; Slaaf, D W; van Donkelaar, C C; Drost, M R; Janssen, J D; Huson, A

1997-09-01

Mechanical interaction between tissue stress and blood perfusion in skeletal muscles plays an important role in blood flow impediment during sustained contraction. The exact mechanism of this interaction is not clear, and experimental investigation of this mechanism is difficult. We developed a finite-element model of the mechanical behavior of blood-perfused muscle tissue, which accounts for mechanical blood-tissue interaction in maximally vasodilated vasculature. Verification of the model was performed by comparing finite-element results of blood pressure and flow with experimental measurements in a muscle that is subject to well-controlled mechanical loading conditions. In addition, we performed simulations of blood perfusion during tetanic, isometric contraction and maximal vasodilation in a simplified, two-dimensional finite-element model of a rat calf muscle. A vascular waterfall in the venous compartment was identified as the main cause for blood flow impediment both in the experiment and in the finite-element simulations. The validated finite-element model offers possibilities for detailed analysis of blood perfusion in three-dimensional muscle models under complicated loading conditions.

Finite element modelling of aluminum alloy 2024-T3 under transverse impact loading

NASA Astrophysics Data System (ADS)

Abdullah, Ahmad Sufian; Kuntjoro, Wahyu; Yamin, A. F. M.

2017-12-01

Fiber metal laminate named GLARE is a new aerospace material which has great potential to be widely used in future lightweight aircraft. It consists of aluminum alloy 2024-T3 and glass-fiber reinforced laminate. In order to produce reliable finite element model of impact response or crashworthiness of structure made of GLARE, one can initially model and validate the finite element model of the impact response of its constituents separately. The objective of this study was to develop a reliable finite element model of aluminum alloy 2024-T3 under low velocity transverse impact loading using commercial software ABAQUS. Johnson-Cook plasticity and damage models were used to predict the alloy's material properties and impact behavior. The results of the finite element analysis were compared to the experiment that has similar material and impact conditions. Results showed good correlations in terms of impact forces, deformation and failure progressions which concluded that the finite element model of 2024-T3 aluminum alloy under low velocity transverse impact condition using Johnson-Cook plastic and damage models was reliable.

Transport Modeling of Hydrogen in Metals for Application to Hydrogen Assisted Cracking of Metals.

DTIC Science & Technology

1995-04-04

34 consists of a Fortran "user element" subroutine for use with the ABAQUS 2 finite element program. Documentation of the 1-D user element subroutine is...trapping theory. The use of the ABAQUS finite element "User Element" subroutines for solving 1-D problems is then outlined in full detail. This is followed...reflect the new ordering given by Eq. (57). ABAOUS User Element Subroutines ABAQUS executes a Fortran subroutine named UEL for each "user defined" finite

REVERSING CYCLIC ELASTO-PLASTIC DEMANDS ON STRUCTURES DURING STRONG MOTION EARTHQUAKE EXCITATION.

USGS Publications Warehouse

Perez, V.; Brady, A.G.; Safak, E.

1986-01-01

Using the horizontal components from El Centro 1940, Taft 1952, and 4 accelerograms from the San Fernando earthquake of 2/9/71, the time history of the elasto-plastic displacement response was calculated for oscillators having periods within the range of 1 to 6 s and ductility factors within the range of 3 to 6. The Nth largest peak of the elasto-plastic response (N equals 2,4,8,16), when expressed as a percentage of maximum response (that is, N equals 1), is fairly independent of period within our period range. When considering only plastic peaks occurring, sometimes in a one-directional group of peaks, in the reverse direction from the preceding plastic peak, the amplitude of the Nth reversing plastic peak is similar to the Nth elastic peak, regardless of the ductility factor.

Validation of High Displacement Piezoelectric Actuator Finite Element Models

NASA Technical Reports Server (NTRS)

Taleghani, B. K.

2000-01-01

The paper presents the results obtained by using NASTRAN(Registered Trademark) and ANSYS(Regitered Trademark) finite element codes to predict doming of the THUNDER piezoelectric actuators during the manufacturing process and subsequent straining due to an applied input voltage. To effectively use such devices in engineering applications, modeling and characterization are essential. Length, width, dome height, and thickness are important parameters for users of such devices. Therefore, finite element models were used to assess the effects of these parameters. NASTRAN(Registered Trademark) and ANSYS(Registered Trademark) used different methods for modeling piezoelectric effects. In NASTRAN(Registered Trademark), a thermal analogy was used to represent voltage at nodes as equivalent temperatures, while ANSYS(Registered Trademark) processed the voltage directly using piezoelectric finite elements. The results of finite element models were validated by using the experimental results.

Use of system identification techniques for improving airframe finite element models using test data

NASA Technical Reports Server (NTRS)

Hanagud, Sathya V.; Zhou, Weiyu; Craig, James I.; Weston, Neil J.

1993-01-01

A method for using system identification techniques to improve airframe finite element models using test data was developed and demonstrated. The method uses linear sensitivity matrices to relate changes in selected physical parameters to changes in the total system matrices. The values for these physical parameters were determined using constrained optimization with singular value decomposition. The method was confirmed using both simple and complex finite element models for which pseudo-experimental data was synthesized directly from the finite element model. The method was then applied to a real airframe model which incorporated all of the complexities and details of a large finite element model and for which extensive test data was available. The method was shown to work, and the differences between the identified model and the measured results were considered satisfactory.

The NASA/Industry Design Analysis Methods for Vibrations (DAMVIBS) Program - A government overview. [of rotorcraft technology development using finite element method

NASA Technical Reports Server (NTRS)

Kvaternik, Raymond G.

1992-01-01

An overview is presented of government contributions to the program called Design Analysis Methods for Vibrations (DAMV) which attempted to develop finite-element-based analyses of rotorcraft vibrations. NASA initiated the program with a finite-element modeling program for the CH-47D tandem-rotor helicopter. The DAMV program emphasized four areas including: airframe finite-element modeling, difficult components studies, coupled rotor-airframe vibrations, and airframe structural optimization. Key accomplishments of the program include industrywide standards for modeling metal and composite airframes, improved industrial designs for vibrations, and the identification of critical structural contributors to airframe vibratory responses. The program also demonstrated the value of incorporating secondary modeling details to improving correlation, and the findings provide the basis for an improved finite-element-based dynamics design-analysis capability.

Probabilistic finite elements for fracture mechanics

NASA Technical Reports Server (NTRS)

Besterfield, Glen

1988-01-01

The probabilistic finite element method (PFEM) is developed for probabilistic fracture mechanics (PFM). A finite element which has the near crack-tip singular strain embedded in the element is used. Probabilistic distributions, such as expectation, covariance and correlation stress intensity factors, are calculated for random load, random material and random crack length. The method is computationally quite efficient and can be expected to determine the probability of fracture or reliability.

Reliability of Next Generation Power Electronics Packaging Under Concurrent Vibration, Thermal and High Power Loads

DTIC Science & Technology

2008-02-01

combined thermal g effect and initial current field. The model is implemented using Abaqus user element subroutine and verified against the experimental...Finite Element Formulation The proposed model is implemented with ABAQUS general purpose finite element program using thermal -displacement analysis...option. ABAQUS and other commercially available finite element codes do not have the capability to solve general electromigration problem directly. Thermal

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

NASA Technical Reports Server (NTRS)

Bibel, G. D.; Kumar, A.; Reddy, S.; Handschuh, R.

1993-01-01

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Error analysis of finite element method for Poisson–Nernst–Planck equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sun, Yuzhou; Sun, Pengtao; Zheng, Bin

A priori error estimates of finite element method for time-dependent Poisson-Nernst-Planck equations are studied in this work. We obtain the optimal error estimates in L∞(H1) and L2(H1) norms, and suboptimal error estimates in L∞(L2) norm, with linear element, and optimal error estimates in L∞(L2) norm with quadratic or higher-order element, for both semi- and fully discrete finite element approximations. Numerical experiments are also given to validate the theoretical results.

An anisotropic elastoplastic constitutive formulation generalised for orthotropic materials

NASA Astrophysics Data System (ADS)

Mohd Nor, M. K.; Ma'at, N.; Ho, C. S.

2018-03-01

This paper presents a finite strain constitutive model to predict a complex elastoplastic deformation behaviour that involves very high pressures and shockwaves in orthotropic materials using an anisotropic Hill's yield criterion by means of the evolving structural tensors. The yield surface of this hyperelastic-plastic constitutive model is aligned uniquely within the principal stress space due to the combination of Mandel stress tensor and a new generalised orthotropic pressure. The formulation is developed in the isoclinic configuration and allows for a unique treatment for elastic and plastic orthotropy. An isotropic hardening is adopted to define the evolution of plastic orthotropy. The important feature of the proposed hyperelastic-plastic constitutive model is the introduction of anisotropic effect in the Mie-Gruneisen equation of state (EOS). The formulation is further combined with Grady spall failure model to predict spall failure in the materials. The proposed constitutive model is implemented as a new material model in the Lawrence Livermore National Laboratory (LLNL)-DYNA3D code of UTHM's version, named Material Type 92 (Mat92). The combination of the proposed stress tensor decomposition and the Mie-Gruneisen EOS requires some modifications in the code to reflect the formulation of the generalised orthotropic pressure. The validation approach is also presented in this paper for guidance purpose. The \\varvec{ψ} tensor used to define the alignment of the adopted yield surface is first validated. This is continued with an internal validation related to elastic isotropic, elastic orthotropic and elastic-plastic orthotropic of the proposed formulation before a comparison against range of plate impact test data at 234, 450 and {895 ms}^{-1} impact velocities is performed. A good agreement is obtained in each test.

Finite element analysis on the bending condition of truck frame before and after opening

NASA Astrophysics Data System (ADS)

Cai, Kaiwu; Cheng, Wei; Lu, Jifu

2018-05-01

Based on the design parameters of a truck frame, the structure design and model of the truck frame are built. Based on the finite element theory, the load, the type of fatigue and the material parameters of the frame are combined with the semi-trailer. Using finite element analysis software, after a truck frame hole in bending condition for the finite element analysis of comparison, through the analysis found that the truck frame hole under bending condition can meet the strength requirements are very helpful for improving the design of the truck frame.

Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators

NASA Technical Reports Server (NTRS)

Taleghani, Barmac K.; Campbell, Joel F.

1999-01-01

A NASTRAN non-linear finite element model has been developed for predicting the dome heights of THUNDER (THin Layer UNimorph Ferroelectric DrivER) piezoelectric actuators. To analytically validate the finite element model, a comparison was made with a non-linear plate solution using Von Karmen's approximation. A 500 volt input was used to examine the actuator deformation. The NASTRAN finite element model was also compared with experimental results. Four groups of specimens were fabricated and tested. Four different input voltages, which included 120, 160, 200, and 240 Vp-p with a 0 volts offset, were used for this comparison.

Biomechanical investigation of naso-orbitoethmoid trauma by finite element analysis.

PubMed

Huempfner-Hierl, Heike; Schaller, Andreas; Hemprich, Alexander; Hierl, Thomas

2014-11-01

Naso-orbitoethmoid fractures account for 5% of all facial fractures. We used data derived from a white 34-year-old man to make a transient dynamic finite element model, which consisted of about 740 000 elements, to simulate fist-like impacts to this anatomically complex area. Finite element analysis showed a pattern of von Mises stresses beyond the yield criterion of bone that corresponded with fractures commonly seen clinically. Finite element models can be used to simulate injuries to the human skull, and provide information about the pathogenesis of different types of fracture. Copyright © 2014 The British Association of Oral and Maxillofacial Surgeons. Published by Elsevier Ltd. All rights reserved.

A unified approach to the analysis and design of elasto-plastic structures with mechanical contact

NASA Technical Reports Server (NTRS)

Bendsoe, Martin P.; Olhoff, Niels; Taylor, John E.

1990-01-01

With structural design in mind, a new unified variational model has been developed which represents the mechanics of deformation elasto-plasticity with unilateral contact conditions. For a design problem formulated as maximization of the load carrying capacity of a structure under certain constraints, the unified model allows for a simultaneous analysis and design synthesis for a whole range of mechanical behavior.

On the uniqueness of measuring elastoplastic properties from indentation: The indistinguishable mystical materials

NASA Astrophysics Data System (ADS)

Chen, Xi; Ogasawara, Nagahisa; Zhao, Manhong; Chiba, Norimasa

2007-08-01

Indentation is widely used to extract material elastoplastic properties from the measured force-displacement curves. One of the most well-established indentation techniques utilizes dual (or plural) sharp indenters (which have different apex angles) to deduce key parameters such as the elastic modulus, yield stress, and work-hardening exponent for materials that obey the power-law constitutive relationship. However, the uniqueness of such analysis is not yet systematically studied or challenged. Here we show the existence of "mystical materials", which have distinct elastoplastic properties yet they yield almost identical indentation behaviors, even when the indenter angle is varied in a large range. These mystical materials are, therefore, indistinguishable by many existing indentation analyses unless extreme (and often impractical) indenter angles are used. Explicit procedures of deriving these mystical materials are established, and the general characteristics of the mystical materials are discussed. In many cases, for a given indenter angle range, a material would have infinite numbers of mystical siblings, and the existence maps of the mystical materials are also obtained. Furthermore, we propose two alternative techniques to effectively distinguish these mystical materials. The study in this paper addresses the important question of the uniqueness of indentation test, as well as providing useful guidelines to properly use the indentation technique to measure material elastoplastic properties.

Computer aided stress analysis of long bones utilizing computer tomography

DOE Office of Scientific and Technical Information (OSTI.GOV)

Marom, S.A.

1986-01-01

A computer aided analysis method, utilizing computed tomography (CT) has been developed, which together with a finite element program determines the stress-displacement pattern in a long bone section. The CT data file provides the geometry, the density and the material properties for the generated finite element model. A three-dimensional finite element model of a tibial shaft is automatically generated from the CT file by a pre-processing procedure for a finite element program. The developed pre-processor includes an edge detection algorithm which determines the boundaries of the reconstructed cross-sectional images of the scanned bone. A mesh generation procedure than automatically generatesmore » a three-dimensional mesh of a user-selected refinement. The elastic properties needed for the stress analysis are individually determined for each model element using the radiographic density (CT number) of each pixel with the elemental borders. The elastic modulus is determined from the CT radiographic density by using an empirical relationship from the literature. The generated finite element model, together with applied loads, determined from existing gait analysis and initial displacements, comprise a formatted input for the SAP IV finite element program. The output of this program, stresses and displacements at the model elements and nodes, are sorted and displayed by a developed post-processor to provide maximum and minimum values at selected locations in the model.« less

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

NASA Technical Reports Server (NTRS)

Rismantab-Sany, J.; Chang, B.; Shabana, A. A.

1989-01-01

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Finite element probabilistic risk assessment of transmission line insulation flashovers caused by lightning strokes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bacvarov, D.C.

1981-01-01

A new method for probabilistic risk assessment of transmission line insulation flashovers caused by lightning strokes is presented. The utilized approach of applying the finite element method for probabilistic risk assessment is demonstrated to be very powerful. The reasons for this are two. First, the finite element method is inherently suitable for analysis of three dimensional spaces where the parameters, such as three variate probability densities of the lightning currents, are non-uniformly distributed. Second, the finite element method permits non-uniform discretization of the three dimensional probability spaces thus yielding high accuracy in critical regions, such as the area of themore » low probability events, while at the same time maintaining coarse discretization in the non-critical areas to keep the number of grid points and the size of the problem to a manageable low level. The finite element probabilistic risk assessment method presented here is based on a new multidimensional search algorithm. It utilizes an efficient iterative technique for finite element interpolation of the transmission line insulation flashover criteria computed with an electro-magnetic transients program. Compared to other available methods the new finite element probabilistic risk assessment method is significantly more accurate and approximately two orders of magnitude computationally more efficient. The method is especially suited for accurate assessment of rare, very low probability events.« less

Consequences of the presence of a weak fault on the stress and strain within an active margin

NASA Astrophysics Data System (ADS)

Conin, M.; Henry, P.; Godard, V.; Bourlange, S.

2009-12-01

Accreting margins often display an outer thrust and fold belt and an inner forearc domain overlying the subduction plate. Assuming that this overlying material behaves as Coulomb material, the outer wedge and the inner wedge are classically approximated as a critical state and a stable state Coulomb wedge, respectively. Critical Coulomb wedge theory can account for the transition from wedge to forearc. However, it cannot be used to determine the state of stress in the transition zone, nor the consequences of a discontinuity within the margin. The presence of a discontinuity such as a splay fault having a low effective friction coefficient should affect the stress state within the wedge, at least locally around the splay fault. Moreover, the effective friction coefficient of the seismogenic zone is expected to vary during the seismic cycle, and this may influence the stability of the Coulomb wedges. We use the ADELI finite element code (Chery and Hassani, 2000) to model the quasi-static stress and strain of a decollement and splay fault system, within a two dimensional elasto-plastic wedge with Drucker-Prager rheology. The subduction plane, the basal decollement of the accretionary wedge and the splay fault are modeled with contact elements. The modeled margin comprises an inner and an outer domain with distinct tapers and basal friction coefficients. For a given splay fault geometry, we evaluate the friction coefficient threshold for splay fault activation as a function of the basal friction coefficients, and examine the consequences of motion along the splay fault on stress and strain within the wedge and on the surface slope at equilibrium. Friction coefficients are varied in time to mimic the consequence of the seismic cycle on the static stress state and strain distribution. Results show the possibility of coexistence of localized extensional regime above the splay fault within a regional compressional regime. Such coexistence is consistent with stress orientation estimation made from breakouts in the Nankai accretionary prim (Kinoshita et al, 2009).

A computer graphics program for general finite element analyses

NASA Technical Reports Server (NTRS)

Thornton, E. A.; Sawyer, L. M.

1978-01-01

Documentation for a computer graphics program for displays from general finite element analyses is presented. A general description of display options and detailed user instructions are given. Several plots made in structural, thermal and fluid finite element analyses are included to illustrate program options. Sample data files are given to illustrate use of the program.

Solution-adaptive finite element method in computational fracture mechanics

NASA Technical Reports Server (NTRS)

Min, J. B.; Bass, J. M.; Spradley, L. W.

1993-01-01

Some recent results obtained using solution-adaptive finite element method in linear elastic two-dimensional fracture mechanics problems are presented. The focus is on the basic issue of adaptive finite element method for validating the applications of new methodology to fracture mechanics problems by computing demonstration problems and comparing the stress intensity factors to analytical results.

2D and 3D Multiscale/Multicomponent Modeling of Impact Response of Heterogeneous Energetic Composites

DTIC Science & Technology

2016-06-01

7 Development of Cohesive Finite Element Method (CFEM) Capability ................................7 3D...Cohesive Finite Element Method (CFEM) framework A new scientific framework and technical capability is developed for the computational analyses of...this section should shift from reporting activities to reporting accomplishments. Development of Cohesive Finite Element Method (CFEM) Capability

Plane stress analysis of wood members using isoparametric finite elements, a computer program

Treesearch

Gary D. Gerhardt

1983-01-01

A finite element program is presented which computes displacements, strains, and stresses in wood members of arbitrary shape which are subjected to plane strain/stressloading conditions. This report extends a program developed by R. L. Taylor in 1977, by adding both the cubic isoparametric finite element and the capability to analyze nonisotropic materials. The...

Finite element analysis of helicopter structures

NASA Technical Reports Server (NTRS)

Rich, M. J.

1978-01-01

Application of the finite element analysis is now being expanded to three dimensional analysis of mechanical components. Examples are presented for airframe, mechanical components, and composite structure calculations. Data are detailed on the increase of model size, computer usage, and the effect on reducing stress analysis costs. Future applications for use of finite element analysis for helicopter structures are projected.

Benchmark model correction of monitoring system based on Dynamic Load Test of Bridge

NASA Astrophysics Data System (ADS)

Shi, Jing-xian; Fan, Jiang

2018-03-01

Structural health monitoring (SHM) is a field of research in the area, and it’s designed to achieve bridge safety and reliability assessment, which needs to be carried out on the basis of the accurate simulation of the finite element model. Bridge finite element model is simplified of the structural section form, support conditions, material properties and boundary condition, which is based on the design and construction drawings, and it gets the calculation models and the results.But according to the design and specification requirements established finite element model due to its cannot fully reflect the true state of the bridge, so need to modify the finite element model to obtain the more accurate finite element model. Based on Da-guan river crossing of Ma - Zhao highway in Yunnan province as the background to do the dynamic load test test, we find that the impact coefficient of the theoretical model of the bridge is very different from the coefficient of the actual test, and the change is different; according to the actual situation, the calculation model is adjusted to get the correct frequency of the bridge, the revised impact coefficient found that the modified finite element model is closer to the real state, and provides the basis for the correction of the finite model.

Adaptive mesh refinement for time-domain electromagnetics using vector finite elements :a feasibility study.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Turner, C. David; Kotulski, Joseph Daniel; Pasik, Michael Francis

This report investigates the feasibility of applying Adaptive Mesh Refinement (AMR) techniques to a vector finite element formulation for the wave equation in three dimensions. Possible error estimators are considered first. Next, approaches for refining tetrahedral elements are reviewed. AMR capabilities within the Nevada framework are then evaluated. We summarize our conclusions on the feasibility of AMR for time-domain vector finite elements and identify a path forward.

Nanoengineering Testbed for Nanosolar Cell and Piezoelectric Compounds

DTIC Science & Technology

2012-02-29

element mesh. The third model was a 3D finite element mesh that included complete geometric representation of Berkovich tip. This model allows for a...height of the specimen. These simulations suggest the proper specimen size to approximate a body of semi-infinite extent for a given indentation depth...tip nanoindentation model was the third and final finite element mesh created for analysis and comparison. The material model and the finite element

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

NASA Astrophysics Data System (ADS)

Bause, Markus

2008-02-01

In this work we study mixed finite element approximations of Richards' equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61-78]), the Raviart- Thomas mixed finite element method ( RT0) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565-581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373-394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869-1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148-167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701-724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi- Douglas- Marini element ( BDM1) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM1 approach to the RT0 one is observed, which however is less significant for the accompanying solute transport.

Prediction of local proximal tibial subchondral bone structural stiffness using subject-specific finite element modeling: Effect of selected density-modulus relationship.

PubMed

Nazemi, S Majid; Amini, Morteza; Kontulainen, Saija A; Milner, Jaques S; Holdsworth, David W; Masri, Bassam A; Wilson, David R; Johnston, James D

2015-08-01

Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain initiation. Calculation of bone elastic moduli from image data is a basic step when constructing finite element models. However, different relationships between elastic moduli and imaged density (known as density-modulus relationships) have been reported in the literature. The objective of this study was to apply seven different trabecular-specific and two cortical-specific density-modulus relationships from the literature to finite element models of proximal tibia subchondral bone, and identify the relationship(s) that best predicted experimentally measured local subchondral structural stiffness with highest explained variance and least error. Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using published density-modulus relationships and mapped to corresponding finite element models. Proximal tibial structural stiffness values were compared to experimentally measured stiffness values from in-situ macro-indentation testing directly on the subchondral bone surface (47 indentation points). Regression lines between experimentally measured and finite element calculated stiffness had R(2) values ranging from 0.56 to 0.77. Normalized root mean squared error varied from 16.6% to 337.6%. Of the 21 evaluated density-modulus relationships in this study, Goulet combined with Snyder and Schneider or Rho appeared most appropriate for finite element modeling of local subchondral bone structural stiffness. Though, further studies are needed to optimize density-modulus relationships and improve finite element estimates of local subchondral bone structural stiffness. Copyright © 2015 Elsevier Ltd. All rights reserved.

Accuracy of specimen-specific nonlinear finite element analysis for evaluation of distal radius strength in cadaver material.

PubMed

Matsuura, Yusuke; Kuniyoshi, Kazuki; Suzuki, Takane; Ogawa, Yasufumi; Sukegawa, Koji; Rokkaku, Tomoyuki; Takahashi, Kazuhisa

2014-11-01

Distal radius fracture, which often occurs in the setting of osteoporosis, can lead to permanent deformity and disability. Great effort has been directed toward developing noninvasive methods for evaluating the distal radius strength, with the goal of assessing fracture risk. The aim of this study was to evaluate distal radius strength using a finite element model and to gauge the accuracy of finite element model measurement using cadaver material. Ten wrists were obtained from cadavers with a mean age of 89.5 years at death. CT images of each wrist in an extended position were obtained. CT-based finite element models were prepared with Mechanical Finder software. Fracture on the models was simulated by applying a mechanical load to the palm in a direction parallel to the forearm axis, after which the fracture load and the site at which the fracture began were identified. For comparison, the wrists were fractured using a universal testing machine and the fracture load and the site of fracture were identified. The fracture load was 970.9 N in the finite element model group and 990.0 N in the actual measurement group. The site of the initial fracture was extra-articular to the distal radius in both groups. The finite element model was predictive for distal radius fracture when compared to the actual measurement. In this study, a finite element model for evaluation of distal radius strength was validated and can be used to predict fracture risk. We conclude that a finite element model is useful for the evaluation of distal radius strength. Knowing distal radius strength might avoid distal radius fracture because appropriate antiosteoporotic treatment can be initiated.

Experimental validation of finite element modelling of a modular metal-on-polyethylene total hip replacement.

PubMed

Hua, Xijin; Wang, Ling; Al-Hajjar, Mazen; Jin, Zhongmin; Wilcox, Ruth K; Fisher, John

2014-07-01

Finite element models are becoming increasingly useful tools to conduct parametric analysis, design optimisation and pre-clinical testing for hip joint replacements. However, the verification of the finite element model is critically important. The purposes of this study were to develop a three-dimensional anatomic finite element model for a modular metal-on-polyethylene total hip replacement for predicting its contact mechanics and to conduct experimental validation for a simple finite element model which was simplified from the anatomic finite element model. An anatomic modular metal-on-polyethylene total hip replacement model (anatomic model) was first developed and then simplified with reasonable accuracy to a simple modular total hip replacement model (simplified model) for validation. The contact areas on the articulating surface of three polyethylene liners of modular metal-on-polyethylene total hip replacement bearings with different clearances were measured experimentally in the Leeds ProSim hip joint simulator under a series of loading conditions and different cup inclination angles. The contact areas predicted from the simplified model were then compared with that measured experimentally under the same conditions. The results showed that the simplification made for the anatomic model did not change the predictions of contact mechanics of the modular metal-on-polyethylene total hip replacement substantially (less than 12% for contact stresses and contact areas). Good agreements of contact areas between the finite element predictions from the simplified model and experimental measurements were obtained, with maximum difference of 14% across all conditions considered. This indicated that the simplification and assumptions made in the anatomic model were reasonable and the finite element predictions from the simplified model were valid. © IMechE 2014.

Methods for analysis of cracks in three-dimensional solids

NASA Technical Reports Server (NTRS)

Raju, I. S.; Newman, J. C., Jr.

1984-01-01

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Evaluation of the finite element fuel rod analysis code (FRANCO)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, K.; Feltus, M.A.

1994-12-31

Knowledge of temperature distribution in a nuclear fuel rod is required to predict the behavior of fuel elements during operating conditions. The thermal and mechanical properties and performance characteristics are strongly dependent on the temperature, which can vary greatly inside the fuel rod. A detailed model of fuel rod behavior can be described by various numerical methods, including the finite element approach. The finite element method has been successfully used in many engineering applications, including nuclear piping and reactor component analysis. However, fuel pin analysis has traditionally been carried out with finite difference codes, with the exception of Electric Powermore » Research Institute`s FREY code, which was developed for mainframe execution. This report describes FRANCO, a finite element fuel rod analysis code capable of computing temperature disrtibution and mechanical deformation of a single light water reactor fuel rod.« less

Dynamic responses of graphite/epoxy laminated beam to impact of elastic spheres

NASA Technical Reports Server (NTRS)

Sun, C. T.; Wang, T.

1982-01-01

Wave propagation in 90/45/90/-45/902s and 0/45/0/-45/02s laminates of a graphite/epoxy composite due to impact of a steel ball was investigated experimentally and also by using a high order beam finite element. Dynamic strain responses at several locations were obtained using strain gages. The finite element program which incorporated statically determined contact laws was employed to calculate the contact force history as well as the target beam dynamic deformation. The comparison of the finite element solutions with the experimental data indicated that the static contact laws for loading and unloading (developed under this grant) are adequate for the dynamic impact analysis. It was found that for the 0/45/0/-45/02s laminate which has a much larger longitudinal bending rigidity, the use of beam finite elements is not suitable and plate finite element should be used instead.

A new parallel-vector finite element analysis software on distributed-memory computers

NASA Technical Reports Server (NTRS)

Qin, Jiangning; Nguyen, Duc T.

1993-01-01

A new parallel-vector finite element analysis software package MPFEA (Massively Parallel-vector Finite Element Analysis) is developed for large-scale structural analysis on massively parallel computers with distributed-memory. MPFEA is designed for parallel generation and assembly of the global finite element stiffness matrices as well as parallel solution of the simultaneous linear equations, since these are often the major time-consuming parts of a finite element analysis. Block-skyline storage scheme along with vector-unrolling techniques are used to enhance the vector performance. Communications among processors are carried out concurrently with arithmetic operations to reduce the total execution time. Numerical results on the Intel iPSC/860 computers (such as the Intel Gamma with 128 processors and the Intel Touchstone Delta with 512 processors) are presented, including an aircraft structure and some very large truss structures, to demonstrate the efficiency and accuracy of MPFEA.

Finite element analysis (FEA) analysis of the preflex beam

NASA Astrophysics Data System (ADS)

Wan, Lijuan; Gao, Qilang

2017-10-01

The development of finite element analysis (FEA) has been relatively mature, and is one of the important means of structural analysis. This method changes the problem that the research of complex structure in the past needs to be done by a large number of experiments. Through the finite element method, the numerical simulation of the structure can be used to achieve a variety of static and dynamic simulation analysis of the mechanical problems, it is also convenient to study the parameters of the structural parameters. Combined with a certain number of experiments to verify the simulation model can be completed in the past all the needs of experimental research. The nonlinear finite element method is used to simulate the flexural behavior of the prestressed composite beams with corrugated steel webs. The finite element analysis is used to understand the mechanical properties of the structure under the action of bending load.

Research on Finite Element Model Generating Method of General Gear Based on Parametric Modelling

NASA Astrophysics Data System (ADS)

Lei, Yulong; Yan, Bo; Fu, Yao; Chen, Wei; Hou, Liguo

2017-06-01

Aiming at the problems of low efficiency and poor quality of gear meshing in the current mainstream finite element software, through the establishment of universal gear three-dimensional model, and explore the rules of unit and node arrangement. In this paper, a finite element model generation method of universal gear based on parameterization is proposed. Visual Basic program is used to realize the finite element meshing, give the material properties, and set the boundary / load conditions and other pre-processing work. The dynamic meshing analysis of the gears is carried out with the method proposed in this pape, and compared with the calculated values to verify the correctness of the method. The method greatly shortens the workload of gear finite element pre-processing, improves the quality of gear mesh, and provides a new idea for the FEM pre-processing.

Stable finite element approximations of two-phase flow with soluble surfactant

NASA Astrophysics Data System (ADS)

Barrett, John W.; Garcke, Harald; Nürnberg, Robert

2015-09-01

A parametric finite element approximation of incompressible two-phase flow with soluble surfactants is presented. The Navier-Stokes equations are coupled to bulk and surfaces PDEs for the surfactant concentrations. At the interface adsorption, desorption and stress balances involving curvature effects and Marangoni forces have to be considered. A parametric finite element approximation for the advection of the interface, which maintains good mesh properties, is coupled to the evolving surface finite element method, which is used to discretize the surface PDE for the interface surfactant concentration. The resulting system is solved together with standard finite element approximations of the Navier-Stokes equations and of the bulk parabolic PDE for the surfactant concentration. Semidiscrete and fully discrete approximations are analyzed with respect to stability, conservation and existence/uniqueness issues. The approach is validated for simple test cases and for complex scenarios, including colliding drops in a shear flow, which are computed in two and three space dimensions.

Heat transfer model and finite element formulation for simulation of selective laser melting

NASA Astrophysics Data System (ADS)

Roy, Souvik; Juha, Mario; Shephard, Mark S.; Maniatty, Antoinette M.

2017-10-01

A novel approach and finite element formulation for modeling the melting, consolidation, and re-solidification process that occurs in selective laser melting additive manufacturing is presented. Two state variables are introduced to track the phase (melt/solid) and the degree of consolidation (powder/fully dense). The effect of the consolidation on the absorption of the laser energy into the material as it transforms from a porous powder to a dense melt is considered. A Lagrangian finite element formulation, which solves the governing equations on the unconsolidated reference configuration is derived, which naturally considers the effect of the changing geometry as the powder melts without needing to update the simulation domain. The finite element model is implemented into a general-purpose parallel finite element solver. Results are presented comparing to experimental results in the literature for a single laser track with good agreement. Predictions for a spiral laser pattern are also shown.

The application of super wavelet finite element on temperature-pressure coupled field simulation of LPG tank under jet fire

NASA Astrophysics Data System (ADS)

Zhao, Bin

2015-02-01

Temperature-pressure coupled field analysis of liquefied petroleum gas (LPG) tank under jet fire can offer theoretical guidance for preventing the fire accidents of LPG tank, the application of super wavelet finite element on it is studied in depth. First, review of related researches on heat transfer analysis of LPG tank under fire and super wavelet are carried out. Second, basic theory of super wavelet transform is studied. Third, the temperature-pressure coupled model of gas phase and liquid LPG under jet fire is established based on the equation of state, the VOF model and the RNG k-ɛ model. Then the super wavelet finite element formulation is constructed using the super wavelet scale function as interpolating function. Finally, the simulation is carried out, and results show that the super wavelet finite element method has higher computing precision than wavelet finite element method.

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Cooke, C. H.; Blanchard, D. K.

1975-01-01

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

Dislocation dynamics in non-convex domains using finite elements with embedded discontinuities

NASA Astrophysics Data System (ADS)

Romero, Ignacio; Segurado, Javier; LLorca, Javier

2008-04-01

The standard strategy developed by Van der Giessen and Needleman (1995 Modelling Simul. Mater. Sci. Eng. 3 689) to simulate dislocation dynamics in two-dimensional finite domains was modified to account for the effect of dislocations leaving the crystal through a free surface in the case of arbitrary non-convex domains. The new approach incorporates the displacement jumps across the slip segments of the dislocations that have exited the crystal within the finite element analysis carried out to compute the image stresses on the dislocations due to the finite boundaries. This is done in a simple computationally efficient way by embedding the discontinuities in the finite element solution, a strategy often used in the numerical simulation of crack propagation in solids. Two academic examples are presented to validate and demonstrate the extended model and its implementation within a finite element program is detailed in the appendix.

The use of Galerkin finite-element methods to solve mass-transport equations

USGS Publications Warehouse

Grove, David B.

1977-01-01

The partial differential equation that describes the transport and reaction of chemical solutes in porous media was solved using the Galerkin finite-element technique. These finite elements were superimposed over finite-difference cells used to solve the flow equation. Both convection and flow due to hydraulic dispersion were considered. Linear and Hermite cubic approximations (basis functions) provided satisfactory results: however, the linear functions were computationally more efficient for two-dimensional problems. Successive over relaxation (SOR) and iteration techniques using Tchebyschef polynomials were used to solve the sparce matrices generated using the linear and Hermite cubic functions, respectively. Comparisons of the finite-element methods to the finite-difference methods, and to analytical results, indicated that a high degree of accuracy may be obtained using the method outlined. The technique was applied to a field problem involving an aquifer contaminated with chloride, tritium, and strontium-90. (Woodard-USGS)

Computational Challenges in the Analysis of Petrophysics Using Microtomography and Upscaling

NASA Astrophysics Data System (ADS)

Liu, J.; Pereira, G.; Freij-Ayoub, R.; Regenauer-Lieb, K.

2014-12-01

Microtomography provides detailed 3D internal structures of rocks in micro- to tens of nano-meter resolution and is quickly turning into a new technology for studying petrophysical properties of materials. An important step is the upscaling of these properties as micron or sub-micron resolution can only be done on the sample-scale of millimeters or even less than a millimeter. We present here a recently developed computational workflow for the analysis of microstructures including the upscaling of material properties. Computations of properties are first performed using conventional material science simulations at micro to nano-scale. The subsequent upscaling of these properties is done by a novel renormalization procedure based on percolation theory. We have tested the workflow using different rock samples, biological and food science materials. We have also applied the technique on high-resolution time-lapse synchrotron CT scans. In this contribution we focus on the computational challenges that arise from the big data problem of analyzing petrophysical properties and its subsequent upscaling. We discuss the following challenges: 1) Characterization of microtomography for extremely large data sets - our current capability. 2) Computational fluid dynamics simulations at pore-scale for permeability estimation - methods, computing cost and accuracy. 3) Solid mechanical computations at pore-scale for estimating elasto-plastic properties - computational stability, cost, and efficiency. 4) Extracting critical exponents from derivative models for scaling laws - models, finite element meshing, and accuracy. Significant progress in each of these challenges is necessary to transform microtomography from the current research problem into a robust computational big data tool for multi-scale scientific and engineering problems.

The structural stability of lunar lava tubes

NASA Astrophysics Data System (ADS)

Blair, David M.; Chappaz, Loic; Sood, Rohan; Milbury, Colleen; Bobet, Antonio; Melosh, H. Jay; Howell, Kathleen C.; Freed, Andrew M.

2017-01-01

Mounting evidence from the SELENE, LRO, and GRAIL spacecraft suggests the presence of vacant lava tubes under the surface of the Moon. GRAIL evidence, in particular, suggests that some may be more than a kilometer in width. Such large sublunarean structures would be of great benefit to future human exploration of the Moon, providing shelter from the harsh environment at the surface-but could empty lava tubes of this size be stable under lunar conditions? And what is the largest size at which they could remain structurally sound? We address these questions by creating elasto-plastic finite element models of lava tubes using the Abaqus modeling software and examining where there is local material failure in the tube's roof. We assess the strength of the rock body using the Geological Strength Index method with values appropriate to the Moon, assign it a basaltic density derived from a modern re-analysis of lunar samples, and assume a 3:1 width-to-height ratio for the lava tube. Our results show that the stability of a lava tube depends on its width, its roof thickness, and whether the rock comprising the structure begins in a lithostatic or Poisson stress state. With a roof 2 m thick, lava tubes a kilometer or more in width can remain stable, supporting inferences from GRAIL observations. The theoretical maximum size of a lunar lava tube depends on a variety of factors, but given sufficient burial depth (500 m) and an initial lithostatic stress state, our results show that lava tubes up to 5 km wide may be able to remain structurally stable.

A Multiscale Computational Model Combining a Single Crystal Plasticity Constitutive Model with the Generalized Method of Cells (GMC) for Metallic Polycrystals.

PubMed

Ghorbani Moghaddam, Masoud; Achuthan, Ajit; Bednarcyk, Brett A; Arnold, Steven M; Pineda, Evan J

2016-05-04

A multiscale computational model is developed for determining the elasto-plastic behavior of polycrystal metals by employing a single crystal plasticity constitutive model that can capture the microstructural scale stress field on a finite element analysis (FEA) framework. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, the stand-alone GMC is applied for studying simple material microstructures such as a repeating unit cell (RUC) containing single grain or two grains under uniaxial loading conditions. For verification, the results obtained by the stand-alone GMC are compared to those from an analogous FEA model incorporating the same single crystal plasticity constitutive model. This verification is then extended to samples containing tens to hundreds of grains. The results demonstrate that the GMC homogenization combined with the crystal plasticity constitutive framework is a promising approach for failure analysis of structures as it allows for properly predicting the von Mises stress in the entire RUC, in an average sense, as well as in the local microstructural level, i.e. , each individual grain. Two-three orders of saving in computational cost, at the expense of some accuracy in prediction, especially in the prediction of the components of local tensor field quantities and the quantities near the grain boundaries, was obtained with GMC. Finally, the capability of the developed multiscale model linking FEA and GMC to solve real-life-sized structures is demonstrated by successfully analyzing an engine disc component and determining the microstructural scale details of the field quantities.

Simulation of balloon angioplasty in residually stressed blood vessels-Application of a gradient-enhanced fibre damage model.

PubMed

Polindara, César; Waffenschmidt, Tobias; Menzel, Andreas

2016-08-16

In this contribution we study the balloon angioplasty in a residually stressed artery by means of a non-local gradient-enhanced fibre damage model. The balloon angioplasty is a common surgical intervention used to extend or reopen narrowed blood vessels in order to restore the continuous blood flow in, for instance, atherosclerotic arteries. Inelastic, i.e. predominantly damage-related and elastoplastic processes are induced in the artery during its inflation resulting in an irreversible deformation. As a beneficial consequence, provided that the inelastic deformations do not exceed a specific limit, higher deformations can be obtained within the same pressure level and a continuous blood flow can be guaranteed. In order to study the mechanical response of the artery in this scenario, we make use of the non-local gradient-enhanced model proposed in Waffenschmidt et al. (2014). In this contribution, we extend this model to make use of an incompressible format in connection with a Q1Q1P0 finite element implementation. The residual stresses in the artery are also taken into account following the framework presented in Waffenschmidt (2015). From the results it becomes apparent that, when the artery is subjected to radial stresses beyond the physiological range, damage evolution is triggered in the collagen fibres. The impact of the residual stresses on the structural response and on the circumferential stress distribution along the thickness of the arterial wall is also studied. It is observed that the residual stresses have a beneficial effect on the mechanical response of the arterial wall. Copyright © 2016 Elsevier Ltd. All rights reserved.

A Multiscale Computational Model Combining a Single Crystal Plasticity Constitutive Model with the Generalized Method of Cells (GMC) for Metallic Polycrystals

PubMed Central

Ghorbani Moghaddam, Masoud; Achuthan, Ajit; Bednarcyk, Brett A.; Arnold, Steven M.; Pineda, Evan J.

2016-01-01

A multiscale computational model is developed for determining the elasto-plastic behavior of polycrystal metals by employing a single crystal plasticity constitutive model that can capture the microstructural scale stress field on a finite element analysis (FEA) framework. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, the stand-alone GMC is applied for studying simple material microstructures such as a repeating unit cell (RUC) containing single grain or two grains under uniaxial loading conditions. For verification, the results obtained by the stand-alone GMC are compared to those from an analogous FEA model incorporating the same single crystal plasticity constitutive model. This verification is then extended to samples containing tens to hundreds of grains. The results demonstrate that the GMC homogenization combined with the crystal plasticity constitutive framework is a promising approach for failure analysis of structures as it allows for properly predicting the von Mises stress in the entire RUC, in an average sense, as well as in the local microstructural level, i.e., each individual grain. Two–three orders of saving in computational cost, at the expense of some accuracy in prediction, especially in the prediction of the components of local tensor field quantities and the quantities near the grain boundaries, was obtained with GMC. Finally, the capability of the developed multiscale model linking FEA and GMC to solve real-life-sized structures is demonstrated by successfully analyzing an engine disc component and determining the microstructural scale details of the field quantities. PMID:28773458

Finite element, modal co-ordinate analysis of structures subjected to moving loads

NASA Astrophysics Data System (ADS)

Olsson, M.

1985-03-01

Some of the possibilities of the finite element method in the moving load problem are demonstrated. The bridge-vehicle interaction phenomenon is considered by deriving a general bridge-vehicle element which is believed to be novel. This element may be regarded as a finite element with time-dependent and unsymmetric element matrices. The bridge response is formulated in modal co-ordinates thereby reducing the number of equations to be solved within each time step. Illustrative examples are shown for the special case of a beam bridge model and a one-axle vehicle model.

Finite Element A Posteriori Error Estimation for Heat Conduction. Degree awarded by George Washington Univ.

NASA Technical Reports Server (NTRS)

Lang, Christapher G.; Bey, Kim S. (Technical Monitor)

2002-01-01

This research investigates residual-based a posteriori error estimates for finite element approximations of heat conduction in single-layer and multi-layered materials. The finite element approximation, based upon hierarchical modelling combined with p-version finite elements, is described with specific application to a two-dimensional, steady state, heat-conduction problem. Element error indicators are determined by solving an element equation for the error with the element residual as a source, and a global error estimate in the energy norm is computed by collecting the element contributions. Numerical results of the performance of the error estimate are presented by comparisons to the actual error. Two methods are discussed and compared for approximating the element boundary flux. The equilibrated flux method provides more accurate results for estimating the error than the average flux method. The error estimation is applied to multi-layered materials with a modification to the equilibrated flux method to approximate the discontinuous flux along a boundary at the material interfaces. A directional error indicator is developed which distinguishes between the hierarchical modeling error and the finite element error. Numerical results are presented for single-layered materials which show that the directional indicators accurately determine which contribution to the total error dominates.

Dissipation consistent fabric tensor definition from DEM to continuum for granular media

NASA Astrophysics Data System (ADS)

Li, X. S.; Dafalias, Y. F.

2015-05-01

In elastoplastic soil models aimed at capturing the impact of fabric anisotropy, a necessary ingredient is a measure of anisotropic fabric in the form of an evolving tensor. While it is possible to formulate such a fabric tensor based on indirect phenomenological observations at the continuum level, it is more effective and insightful to have the tensor defined first based on direct particle level microstructural observations and subsequently deduce a corresponding continuum definition. A practical means able to provide such observations, at least in the context of fabric evolution mechanisms, is the discrete element method (DEM). Some DEM defined fabric tensors such as the one based on the statistics of interparticle contact normals have already gained widespread acceptance as a quantitative measure of fabric anisotropy among researchers of granular material behavior. On the other hand, a fabric tensor in continuum elastoplastic modeling has been treated as a tensor-valued internal variable whose evolution must be properly linked to physical dissipation. Accordingly, the adaptation of a DEM fabric tensor definition to a continuum constitutive modeling theory must be thermodynamically consistent in regards to dissipation mechanisms. The present paper addresses this issue in detail, brings up possible pitfalls if such consistency is violated and proposes remedies and guidelines for such adaptation within a recently developed Anisotropic Critical State Theory (ACST) for granular materials.

A cell-centered Lagrangian finite volume approach for computing elasto-plastic response of solids in cylindrical axisymmetric geometries

NASA Astrophysics Data System (ADS)

Sambasivan, Shiv Kumar; Shashkov, Mikhail J.; Burton, Donald E.

2013-03-01

A finite volume cell-centered Lagrangian formulation is presented for solving large deformation problems in cylindrical axisymmetric geometries. Since solid materials can sustain significant shear deformation, evolution equations for stress and strain fields are solved in addition to mass, momentum and energy conservation laws. The total strain-rate realized in the material is split into an elastic and plastic response. The elastic and plastic components in turn are modeled using hypo-elastic theory. In accordance with the hypo-elastic model, a predictor-corrector algorithm is employed for evolving the deviatoric component of the stress tensor. A trial elastic deviatoric stress state is obtained by integrating a rate equation, cast in the form of an objective (Jaumann) derivative, based on Hooke's law. The dilatational response of the material is modeled using an equation of state of the Mie-Grüneisen form. The plastic deformation is accounted for via an iterative radial return algorithm constructed from the J2 von Mises yield condition. Several benchmark example problems with non-linear strain hardening and thermal softening yield models are presented. Extensive comparisons with representative Eulerian and Lagrangian hydrocodes in addition to analytical and experimental results are made to validate the current approach.

Simulation of hot spots formation and evolution in HMX

NASA Astrophysics Data System (ADS)

Wang, Cheng; Yang, Tonghui

2017-06-01

In order to study the formation and evolution of hot spots under shock loading, HMX explosives were selected as the object of study for the two-dimensional finite difference numerical simulation. A fifth order finite difference weighted essentially non-oscillatory (WENO) scheme and a third order TVD Runge-Kutta method are utilized for the spatial discretization and the time advance, respectively. The governing equations are based on the fluid elasto-plastic control equations. The Mie-Gruneisen equation of state and the ideal gas equation of state are selected to use in the state equation of the solid explosives and gas material. In order to simplify the calculation of the model, the reaction can be considered to complete in one step. The calculated area is [ 3.0 ×10-5 m ] × [ 3.0 ×10-5 m ] . The radius is 0.6 ×10-5 m, and the internal gas is not involved in the reaction. The calculation area is divided into 300×300 grids and 10 grids are selected from the bottom of each column to give the particle velocity u as the initial condition. In the selected grid, different initial velocity 100m/s and 200m/s are loaded respectively to study the influence of hot spot formation and evolution in different impact intensity.

Modeling crack propagation in polycrystalline microstructure using variational multiscale method

DOE PAGES

Sun, Shang; Sundararaghavan, Veera

2016-01-01

Crack propagation in a polycrystalline microstructure is analyzed using a novel multiscale model. The model includes an explicit microstructural representation at critical regions (stress concentrators such as notches and cracks) and a reduced order model that statistically captures the microstructure at regions far away from stress concentrations. Crack propagation is modeled in these critical regions using the variational multiscale method. In this approach, a discontinuous displacement field is added to elements that exceed the critical values of normal or tangential tractions during loading. Compared to traditional cohesive zone modeling approaches, the method does not require the use of any specialmore » interface elements in the microstructure and thus can model arbitrary crack paths. As a result, the capability of the method in predicting both intergranular and transgranular failure modes in an elastoplastic polycrystal is demonstrated under tensile and three-point bending loads.« less

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

NASA Technical Reports Server (NTRS)

Thornton, E. A.

1976-01-01

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Supercomputer implementation of finite element algorithms for high speed compressible flows

NASA Technical Reports Server (NTRS)

Thornton, E. A.; Ramakrishnan, R.

1986-01-01

Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes.

Boundary-integral methods in elasticity and plasticity. [solutions of boundary value problems

NASA Technical Reports Server (NTRS)

Mendelson, A.

1973-01-01

Recently developed methods that use boundary-integral equations applied to elastic and elastoplastic boundary value problems are reviewed. Direct, indirect, and semidirect methods using potential functions, stress functions, and displacement functions are described. Examples of the use of these methods for torsion problems, plane problems, and three-dimensional problems are given. It is concluded that the boundary-integral methods represent a powerful tool for the solution of elastic and elastoplastic problems.

Development of an hp-version finite element method for computational optimal control

NASA Technical Reports Server (NTRS)

Hodges, Dewey H.; Warner, Michael S.

1993-01-01

The purpose of this research effort is to develop a means to use, and to ultimately implement, hp-version finite elements in the numerical solution of optimal control problems. The hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element.

Plane elasto-plastic analysis of v-notched plate under bending by boundary integral equation method. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Rzasnicki, W.

1973-01-01

A method of solution is presented, which, when applied to the elasto-plastic analysis of plates having a v-notch on one edge and subjected to pure bending, will produce stress and strain fields in much greater detail than presently available. Application of the boundary integral equation method results in two coupled Fredholm-type integral equations, subject to prescribed boundary conditions. These equations are replaced by a system of simultaneous algebraic equations and solved by a successive approximation method employing Prandtl-Reuss incremental plasticity relations. The method is first applied to number of elasto-static problems and the results compared with available solutions. Good agreement is obtained in all cases. The elasto-plastic analysis provides detailed stress and strain distributions for several cases of plates with various notch angles and notch depths. A strain hardening material is assumed and both plane strain and plane stress conditions are considered.

Study of the Anisotropic Elastoplastic Properties of β-Ga2O3 Films Synthesized on SiC/Si Substrates

NASA Astrophysics Data System (ADS)

Grashchenko, A. S.; Kukushkin, S. A.; Nikolaev, V. I.; Osipov, A. V.; Osipova, E. V.; Soshnikov, I. P.

2018-05-01

The structural and mechanical properties of gallium oxide films grown on silicon crystallographic planes (001), (011), and (111) with a buffer layer of silicon carbide are investigated. Nanoindentation was used to study the elastoplastic properties of gallium oxide and also to determine the elastic recovery parameter of the films under study. The tensile strength, hardness, elasticity tensor, compliance tensor, Young's modulus, Poisson's ratio, and other characteristics of gallium oxide were calculated using quantum chemistry methods. It was found that the gallium oxide crystal is auxetic because, for some stretching directions, the Poisson's ratio takes on negative values. The calculated values correspond quantitatively to the experimental data. It is concluded that the elastoplastic properties of gallium oxide films approximately correspond to the properties of bulk crystals and that a change in the orientation of the silicon surface leads to a significant change in the orientation of gallium oxide.

Elasto-plastic bending of cracked plates, including the effects of crack closure. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Jones, D. P.

1972-01-01

A capability for solving elasto-plastic plate bending problems is developed using assumptions consistent with Kirchhoff plate theory. Both bending and extensional modes of deformation are admitted with the two modes becoming coupled as yielding proceeds. Equilibrium solutions are obtained numerically by determination of the stationary point of a functional which is analogous to the potential strain energy. The stationary value of the functional for each load increment is efficiently obtained through use of the conjugate gradient. This technique is applied to the problem of a large centrally through cracked plate subject to remote circular bending. Comparison is drawn between two cases of the bending problem. The first neglects the possibility of crack face interference with bending, and the second includes a kinematic prohibition against the crack face from passing through the symmetry plane. Results are reported which isolate the effects of elastoplastic flow and crack closure.

A mixed shear flexible finite element for the analysis of laminated plates

NASA Technical Reports Server (NTRS)

Putcha, N. S.; Reddy, J. N.

1984-01-01

A mixed shear flexible finite element based on the Hencky-Mindlin type shear deformation theory of laminated plates is presented and their behavior in bending is investigated. The element consists of three displacements, two rotations, and three moments as the generalized degrees of freedom per node. The numerical convergence and accuracy characteristics of the element are investigated by comparing the finite element solutions with the exact solutions. The present study shows that reduced-order integration of the stiffness coefficients due to shear is necessary to obtain accurate results for thin plates.

Expansion and improvement of the FORMA system for response and load analysis. Volume 2C: Listings, finite element FORMA subroutines

NASA Technical Reports Server (NTRS)

Wohlen, R. L.

1976-01-01

A listing of the source deck of each finite element FORMA subroutine is given to remove the 'black-box' aura of the subroutines so that the analyst may better understand the detailed operations of each subroutine. The FORTRAN 4 programming language is used in all finite element FORMA subroutines.

Constitutive Model Calibration via Autonomous Multiaxial Experimentation (Postprint)

DTIC Science & Technology

2016-09-17

test machine. Experimental data is reduced and finite element simulations are conducted in parallel with the test based on experimental strain...data is reduced and finite element simulations are conducted in parallel with the test based on experimental strain conditions. Optimization methods...be used directly in finite element simulations of more complex geometries. Keywords Axial/torsional experimentation • Plasticity • Constitutive model

Toward automatic finite element analysis

NASA Technical Reports Server (NTRS)

Kela, Ajay; Perucchio, Renato; Voelcker, Herbert

1987-01-01

Two problems must be solved if the finite element method is to become a reliable and affordable blackbox engineering tool. Finite element meshes must be generated automatically from computer aided design databases and mesh analysis must be made self-adaptive. The experimental system described solves both problems in 2-D through spatial and analytical substructuring techniques that are now being extended into 3-D.

Numerical Analysis of Solids at Failure

DTIC Science & Technology

2011-08-20

failure analyses include the formulation of invariant finite elements for thin Kirchhoff rods, and preliminary initial studies of growth in...analysis of the failure of other structural/mechanical systems, including the finite element modeling of thin Kirchhoff rods and the constitutive...algorithm based on the connectivity graph of the underlying finite element mesh. In this setting, the discontinuities are defined by fronts propagating

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wiley, J.C.

The author describes a general `hp` finite element method with adaptive grids. The code was based on the work of Oden, et al. The term `hp` refers to the method of spatial refinement (h), in conjunction with the order of polynomials used as a part of the finite element discretization (p). This finite element code seems to handle well the different mesh grid sizes occuring between abuted grids with different resolutions.

A finite element conjugate gradient FFT method for scattering

NASA Technical Reports Server (NTRS)

Collins, Jeffery D.; Ross, Dan; Jin, J.-M.; Chatterjee, A.; Volakis, John L.

1991-01-01

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

Computer animation of modal and transient vibrations

NASA Technical Reports Server (NTRS)

Lipman, Robert R.

1987-01-01

An interactive computer graphics processor is described that is capable of generating input to animate modal and transient vibrations of finite element models on an interactive graphics system. The results from NASTRAN can be postprocessed such that a three dimensional wire-frame picture, in perspective, of the finite element mesh is drawn on the graphics display. Modal vibrations of any mode shape or transient motions over any range of steps can be animated. The finite element mesh can be color-coded by any component of displacement. Viewing parameters and the rate of vibration of the finite element model can be interactively updated while the structure is vibrating.

Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms

NASA Technical Reports Server (NTRS)

Kurdila, Andrew J.; Sharpley, Robert C.

1999-01-01

This paper presents a final report on Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms. The focus of this research is to derive and implement: 1) Wavelet based methodologies for the compression, transmission, decoding, and visualization of three dimensional finite element geometry and simulation data in a network environment; 2) methodologies for interactive algorithm monitoring and tracking in computational mechanics; and 3) Methodologies for interactive algorithm steering for the acceleration of large scale finite element simulations. Also included in this report are appendices describing the derivation of wavelet based Particle Image Velocity algorithms and reduced order input-output models for nonlinear systems by utilizing wavelet approximations.

Finite Element Analysis Generates an Increasing Interest in Dental Research: A Bibliometric Study.

PubMed

Diarra, Abdoulaziz; Mushegyan, Vagan; Naveau, Adrien

2016-01-01

The purpose was to provide a longitudinal overview of published studies that use finite element analysis in dental research, by using the SCI-expanded database of Web of Science(®) (Thomson Reuters). Eighty publications from 1999-2000 and 473 from 2009-2010 were retrieved. This literature grew faster than the overall dental literature. The number of publishing countries doubled. The main journals were American or English, and dealt with implantology. For the top 10 journals publishing dental finite element papers, the mean impact factor increased by 75% during the decade. Finite elements generate an increasing interest from dental authors and publishers worldwide.

On the Finite Element Implementation of the Generalized Method of Cells Micromechanics Constitutive Model

NASA Technical Reports Server (NTRS)

Wilt, T. E.

1995-01-01

The Generalized Method of Cells (GMC), a micromechanics based constitutive model, is implemented into the finite element code MARC using the user subroutine HYPELA. Comparisons in terms of transverse deformation response, micro stress and strain distributions, and required CPU time are presented for GMC and finite element models of fiber/matrix unit cell. GMC is shown to provide comparable predictions of the composite behavior and requires significantly less CPU time as compared to a finite element analysis of the unit cell. Details as to the organization of the HYPELA code are provided with the actual HYPELA code included in the appendix.

Finite element methods on supercomputers - The scatter-problem

NASA Technical Reports Server (NTRS)

Loehner, R.; Morgan, K.

1985-01-01

Certain problems arise in connection with the use of supercomputers for the implementation of finite-element methods. These problems are related to the desirability of utilizing the power of the supercomputer as fully as possible for the rapid execution of the required computations, taking into account the gain in speed possible with the aid of pipelining operations. For the finite-element method, the time-consuming operations may be divided into three categories. The first two present no problems, while the third type of operation can be a reason for the inefficient performance of finite-element programs. Two possibilities for overcoming certain difficulties are proposed, giving attention to a scatter-process.

Finite elements: Theory and application

NASA Technical Reports Server (NTRS)

Dwoyer, D. L. (Editor); Hussaini, M. Y. (Editor); Voigt, R. G. (Editor)

1988-01-01

Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

NASA Technical Reports Server (NTRS)

Atluri, S. N.; Nakagaki, M.; Kathiresan, K.

1980-01-01

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Shear-flexible finite-element models of laminated composite plates and shells

NASA Technical Reports Server (NTRS)

Noor, A. K.; Mathers, M. D.

1975-01-01

Several finite-element models are applied to the linear static, stability, and vibration analysis of laminated composite plates and shells. The study is based on linear shallow-shell theory, with the effects of shear deformation, anisotropic material behavior, and bending-extensional coupling included. Both stiffness (displacement) and mixed finite-element models are considered. Discussion is focused on the effects of shear deformation and anisotropic material behavior on the accuracy and convergence of different finite-element models. Numerical studies are presented which show the effects of increasing the order of the approximating polynomials, adding internal degrees of freedom, and using derivatives of generalized displacements as nodal parameters.

Application of Finite Element Method in Traffic Injury and Its Prospect in Forensic Science.

PubMed

Liu, C G; Lu, Y J; Gao, J; Liu, Q

2016-06-01

The finite element method （FEM） is a numerical computation method based on computer technology, and has been gradually applied in the fields of medicine and biomechanics. The finite element analysis can be used to explore the loading process and injury mechanism of human body in traffic injury. FEM is also helpful for the forensic investigation in traffic injury. This paper reviews the development of the finite element models and analysis of brain, cervical spine, chest and abdomen, pelvis, limbs at home and aboard in traffic injury in recent years. Copyright© by the Editorial Department of Journal of Forensic Medicine.

Tools for Modeling & Simulation of Molecular and Nanoelectronics Devices

DTIC Science & Technology

2012-06-14

implemented a prototype DFT simulation software using two different open source Finite Element (FE) libraries: DEALII and FENICS . These two libraries have been...ATK. In the first part of this Phase I project we investigated two different candidate finite element libraries, DEAL II and FENICS . Although both...element libraries, Deal.II and FEniCS /dolfin, for use as back-ends to a finite element DFT in ATK, Quantum Insight and QuantumWise A/S, October 2011.

Review of literature on the finite-element solution of the equations of two-dimensional surface-water flow in the horizontal plane

USGS Publications Warehouse

Lee, Jonathan K.; Froehlich, David C.

1987-01-01

Published literature on the application of the finite-element method to solving the equations of two-dimensional surface-water flow in the horizontal plane is reviewed in this report. The finite-element method is ideally suited to modeling two-dimensional flow over complex topography with spatially variable resistance. A two-dimensional finite-element surface-water flow model with depth and vertically averaged velocity components as dependent variables allows the user great flexibility in defining geometric features such as the boundaries of a water body, channels, islands, dikes, and embankments. The following topics are reviewed in this report: alternative formulations of the equations of two-dimensional surface-water flow in the horizontal plane; basic concepts of the finite-element method; discretization of the flow domain and representation of the dependent flow variables; treatment of boundary conditions; discretization of the time domain; methods for modeling bottom, surface, and lateral stresses; approaches to solving systems of nonlinear equations; techniques for solving systems of linear equations; finite-element alternatives to Galerkin's method of weighted residuals; techniques of model validation; and preparation of model input data. References are listed in the final chapter.

[Construction of platform on the three-dimensional finite element model of the dentulous mandibular body of a normal person].

PubMed

Gong, Lu-Lu; Zhu, Jing; Ding, Zu-Quan; Li, Guo-Qiang; Wang, Li-Ming; Yan, Bo-Yong

2008-04-01

To develop a method to construct a three-dimensional finite element model of the dentulous mandibular body of a normal person. A series of pictures with the interval of 0.1 mm were taken by CT scanning. After extracting the coordinates of key points of some pictures by the procedure, we used a C program to process the useful data, and constructed a platform of the three-dimensional finite element model of the dentulous mandibular body with the Ansys software for finite element analysis. The experimental results showed that the platform of the three-dimensional finite element model of the dentulous mandibular body was more accurate and applicable. The exact three-dimensional shape of model was well constructed, and each part of this model, such as one single tooth, can be deleted, which can be used to emulate various tooth-loss clinical cases. The three-dimensional finite element model is constructed with life-like shapes of dental cusps. Each part of this model can be easily removed. In conclusion, this experiment provides a good platform of biomechanical analysis on various tooth-loss clinical cases.

The Applications of Finite Element Analysis in Proximal Humeral Fractures.

PubMed

Ye, Yongyu; You, Wei; Zhu, Weimin; Cui, Jiaming; Chen, Kang; Wang, Daping

2017-01-01

Proximal humeral fractures are common and most challenging, due to the complexity of the glenohumeral joint, especially in the geriatric population with impacted fractures, that the development of implants continues because currently the problems with their fixation are not solved. Pre-, intra-, and postoperative assessments are crucial in management of those patients. Finite element analysis, as one of the valuable tools, has been implemented as an effective and noninvasive method to analyze proximal humeral fractures, providing solid evidence for management of troublesome patients. However, no review article about the applications and effects of finite element analysis in assessing proximal humeral fractures has been reported yet. This review article summarized the applications, contribution, and clinical significance of finite element analysis in assessing proximal humeral fractures. Furthermore, the limitations of finite element analysis, the difficulties of more realistic simulation, and the validation and also the creation of validated FE models were discussed. We concluded that although some advancements in proximal humeral fractures researches have been made by using finite element analysis, utility of this powerful tool for routine clinical management and adequate simulation requires more state-of-the-art studies to provide evidence and bases.

The Elastic Behaviour of Sintered Metallic Fibre Networks: A Finite Element Study by Beam Theory

PubMed Central

Bosbach, Wolfram A.

2015-01-01

Background The finite element method has complimented research in the field of network mechanics in the past years in numerous studies about various materials. Numerical predictions and the planning efficiency of experimental procedures are two of the motivational aspects for these numerical studies. The widespread availability of high performance computing facilities has been the enabler for the simulation of sufficiently large systems. Objectives and Motivation In the present study, finite element models were built for sintered, metallic fibre networks and validated by previously published experimental stiffness measurements. The validated models were the basis for predictions about so far unknown properties. Materials and Methods The finite element models were built by transferring previously published skeletons of fibre networks into finite element models. Beam theory was applied as simplification method. Results and Conclusions The obtained material stiffness isn’t a constant but rather a function of variables such as sample size and boundary conditions. Beam theory offers an efficient finite element method for the simulated fibre networks. The experimental results can be approximated by the simulated systems. Two worthwhile aspects for future work will be the influence of size and shape and the mechanical interaction with matrix materials. PMID:26569603

An accurate nonlinear finite element analysis and test correlation of a stiffened composite wing panel

NASA Astrophysics Data System (ADS)

Davis, D. D., Jr.; Krishnamurthy, T.; Stroud, W. J.; McCleary, S. L.

1991-05-01

State-of-the-art nonlinear finite element analysis techniques are evaluated by applying them to a realistic aircraft structural component. A wing panel from the V-22 tiltrotor aircraft is chosen because it is a typical modern aircraft structural component for which there is experimental data for comparison of results. From blueprints and drawings, a very detailed finite element model containing 2284 9-node Assumed Natural-Coordinate Strain elements was generated. A novel solution strategy which accounts for geometric nonlinearity through the use of corotating element reference frames and nonlinear strain-displacement relations is used to analyze this detailed model. Results from linear analyses using the same finite element model are presented in order to illustrate the advantages and costs of the nonlinear analysis as compared with the more traditional linear analysis.

An accurate nonlinear finite element analysis and test correlation of a stiffened composite wing panel

NASA Technical Reports Server (NTRS)

Davis, D. D., Jr.; Krishnamurthy, T.; Stroud, W. J.; Mccleary, S. L.

1991-01-01

State-of-the-art nonlinear finite element analysis techniques are evaluated by applying them to a realistic aircraft structural component. A wing panel from the V-22 tiltrotor aircraft is chosen because it is a typical modern aircraft structural component for which there is experimental data for comparison of results. From blueprints and drawings, a very detailed finite element model containing 2284 9-node Assumed Natural-Coordinate Strain elements was generated. A novel solution strategy which accounts for geometric nonlinearity through the use of corotating element reference frames and nonlinear strain-displacement relations is used to analyze this detailed model. Results from linear analyses using the same finite element model are presented in order to illustrate the advantages and costs of the nonlinear analysis as compared with the more traditional linear analysis.

Application of 'steady' state finite element and transient finite difference theory to sound propagation in a variable duct - A comparison with experiment

NASA Technical Reports Server (NTRS)

Baumeister, K. J.; Eversman, W.; Astley, R. J.; White, J. W.

1981-01-01

Experimental data are presented for sound propagation in a simulated infinite hard wall duct with a large change in duct cross sectional area. The data are conveniently tabulated for further use. The 'steady' state finite element theory of Astley and Eversman (1981) and the transient finite difference theory of White (1981) are in good agreement with the data for both the axial and transverse pressure profiles and the axial phase angle. Therefore, numerical finite difference and finite element theories appear to be ideally suited for handling duct propagation problems which encounter large axial gradients in acoustic parameters. The measured energy reflection coefficient agrees with the values from the Astley-Eversman modal coupling model.

The Use of Sparse Direct Solver in Vector Finite Element Modeling for Calculating Two Dimensional (2-D) Magnetotelluric Responses in Transverse Electric (TE) Mode

NASA Astrophysics Data System (ADS)

Yihaa Roodhiyah, Lisa’; Tjong, Tiffany; Nurhasan; Sutarno, D.

2018-04-01

The late research, linear matrices of vector finite element in two dimensional(2-D) magnetotelluric (MT) responses modeling was solved by non-sparse direct solver in TE mode. Nevertheless, there is some weakness which have to be improved especially accuracy in the low frequency (10-3 Hz-10-5 Hz) which is not achieved yet and high cost computation in dense mesh. In this work, the solver which is used is sparse direct solver instead of non-sparse direct solverto overcome the weaknesses of solving linear matrices of vector finite element metod using non-sparse direct solver. Sparse direct solver will be advantageous in solving linear matrices of vector finite element method because of the matrix properties which is symmetrical and sparse. The validation of sparse direct solver in solving linear matrices of vector finite element has been done for a homogen half-space model and vertical contact model by analytical solution. Thevalidation result of sparse direct solver in solving linear matrices of vector finite element shows that sparse direct solver is more stable than non-sparse direct solver in computing linear problem of vector finite element method especially in low frequency. In the end, the accuracy of 2D MT responses modelling in low frequency (10-3 Hz-10-5 Hz) has been reached out under the efficient allocation memory of array and less computational time consuming.

A meta-model analysis of a finite element simulation for defining poroelastic properties of intervertebral discs.

PubMed

Nikkhoo, Mohammad; Hsu, Yu-Chun; Haghpanahi, Mohammad; Parnianpour, Mohamad; Wang, Jaw-Lin

2013-06-01

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 meta-model 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 anulus fibrosus and nucleus pulposus is critical for the quality of the simulation. We developed a material property updating protocol, which is basically a fitting algorithm consisted of finite element simulations and a quadratic response surface regression. This protocol was used to find the material properties, such as the hydraulic permeability, elastic modulus, and Poisson's ratio, of intact and degenerated porcine discs. The results showed that the in vitro disc experimental deformations were well fitted with limited finite element simulations and a quadratic response surface regression. The comparison of material properties of intact and degenerated discs showed that the hydraulic permeability significantly decreased but Poisson's ratio significantly increased for the degenerated discs. This study shows that the developed protocol is efficient and effective in defining material properties of a complex structure such as the intervertebral disc.

Micro-scale finite element modeling of ultrasound propagation in aluminum trabecular bone-mimicking phantoms: A comparison between numerical simulation and experimental results.

PubMed

Vafaeian, B; Le, L H; Tran, T N H T; El-Rich, M; El-Bialy, T; Adeeb, S

2016-05-01

The present study investigated the accuracy of micro-scale finite element modeling for simulating broadband ultrasound propagation in water-saturated trabecular bone-mimicking phantoms. To this end, five commercially manufactured aluminum foam samples as trabecular bone-mimicking phantoms were utilized for ultrasonic immersion through-transmission experiments. Based on micro-computed tomography images of the same physical samples, three-dimensional high-resolution computational samples were generated to be implemented in the micro-scale finite element models. The finite element models employed the standard Galerkin finite element method (FEM) in time domain to simulate the ultrasonic experiments. The numerical simulations did not include energy dissipative mechanisms of ultrasonic attenuation; however, they expectedly simulated reflection, refraction, scattering, and wave mode conversion. The accuracy of the finite element simulations were evaluated by comparing the simulated ultrasonic attenuation and velocity with the experimental data. The maximum and the average relative errors between the experimental and simulated attenuation coefficients in the frequency range of 0.6-1.4 MHz were 17% and 6% respectively. Moreover, the simulations closely predicted the time-of-flight based velocities and the phase velocities of ultrasound with maximum relative errors of 20 m/s and 11 m/s respectively. The results of this study strongly suggest that micro-scale finite element modeling can effectively simulate broadband ultrasound propagation in water-saturated trabecular bone-mimicking structures. Copyright © 2016 Elsevier B.V. All rights reserved.

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

NASA Technical Reports Server (NTRS)

Kandil, Osama A.

1998-01-01

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Eigenvalues of Rectangular Waveguide Using FEM With Hybrid Elements

NASA Technical Reports Server (NTRS)

Deshpande, Manohar D.; Hall, John M.

2002-01-01

A finite element analysis using hybrid triangular-rectangular elements is developed to estimate eigenvalues of a rectangular waveguide. Use of rectangular vector-edge finite elements in the vicinity of the PEC boundary and triangular elements in the interior region more accurately models the physical nature of the electromagnetic field, and consequently quicken the convergence.

Simulation of thin slot spirals and dual circular patch antennas using the finite element method with mixed elements

NASA Technical Reports Server (NTRS)

Gong, Jian; Volakis, John L.; Nurnberger, Michael W.

1995-01-01

This semi-annual report describes progress up to mid-January 1995. The report contains five sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. Of significance is the development of decomposition schemes which separate the different regions of the antenna volume. Substantial effort was devoted to improving the feed model in the context of the finite element method (FEM). Finally, an innovative scheme for truncating finite element meshes is presented.

An Analysis of a Finite Element Method for Convection-Diffusion Problems. Part II. A Posteriori Error Estimates and Adaptivity.

DTIC Science & Technology

1983-03-01

AN ANALYSIS OF A FINITE ELEMENT METHOD FOR CONVECTION- DIFFUSION PROBLEMS PART II: A POSTERIORI ERROR ESTIMATES AND ADAPTIVITY by W. G. Szymczak Y 6a...PERIOD COVERED AN ANALYSIS OF A FINITE ELEMENT METHOD FOR final life of the contract CONVECTION- DIFFUSION PROBLEM S. Part II: A POSTERIORI ERROR ...Element Method for Convection- Diffusion Problems. Part II: A Posteriori Error Estimates and Adaptivity W. G. Szvmczak and I. Babu~ka# Laboratory for

Vertical discretization with finite elements for a global hydrostatic model on the cubed sphere

NASA Astrophysics Data System (ADS)

Yi, Tae-Hyeong; Park, Ja-Rin

2017-06-01

A formulation of Galerkin finite element with basis-spline functions on a hybrid sigma-pressure coordinate is presented to discretize the vertical terms of global Eulerian hydrostatic equations employed in a numerical weather prediction system, which is horizontally discretized with high-order spectral elements on a cubed sphere grid. This replaces the vertical discretization of conventional central finite difference that is first-order accurate in non-uniform grids and causes numerical instability in advection-dominant flows. Therefore, a model remains in the framework of Galerkin finite elements for both the horizontal and vertical spatial terms. The basis-spline functions, obtained from the de-Boor algorithm, are employed to derive both the vertical derivative and integral operators, since Eulerian advection terms are involved. These operators are used to discretize the vertical terms of the prognostic and diagnostic equations. To verify the vertical discretization schemes and compare their performance, various two- and three-dimensional idealized cases and a hindcast case with full physics are performed in terms of accuracy and stability. It was shown that the vertical finite element with the cubic basis-spline function is more accurate and stable than that of the vertical finite difference, as indicated by faster residual convergence, fewer statistical errors, and reduction in computational mode. This leads to the general conclusion that the overall performance of a global hydrostatic model might be significantly improved with the vertical finite element.

Finite Element Analysis (FEA) in Design and Production.

ERIC Educational Resources Information Center

Waggoner, Todd C.; And Others

1995-01-01

Finite element analysis (FEA) enables industrial designers to analyze complex components by dividing them into smaller elements, then assessing stress and strain characteristics. Traditionally mainframe based, FEA is being increasingly used in microcomputers. (SK)

High-Accuracy Finite Element Method: Benchmark Calculations

NASA Astrophysics Data System (ADS)

Gusev, Alexander; Vinitsky, Sergue; Chuluunbaatar, Ochbadrakh; Chuluunbaatar, Galmandakh; Gerdt, Vladimir; Derbov, Vladimir; Góźdź, Andrzej; Krassovitskiy, Pavel

2018-02-01

We describe a new high-accuracy finite element scheme with simplex elements for solving the elliptic boundary-value problems and show its efficiency on benchmark solutions of the Helmholtz equation for the triangle membrane and hypercube.

Workshop on the Integration of Finite Element Modeling with Geometric Modeling

NASA Technical Reports Server (NTRS)

Wozny, Michael J.

1987-01-01

The workshop on the Integration of Finite Element Modeling with Geometric Modeling was held on 12 May 1987. It was held to discuss the geometric modeling requirements of the finite element modeling process and to better understand the technical aspects of the integration of these two areas. The 11 papers are presented except for one for which only the abstract is given.

Effects of Crimped Fiber Paths on Mixed Mode Delamination Behaviors in Woven Fabric Composites

DTIC Science & Technology

2016-09-01

continuum finite - element models. Three variations of a plain-woven fabric architecture—each of which had different crimped fiber paths—were considered... Finite - Element Analysis Fracture Mechanics Fracture Toughness Mixed Modes Strain Energy Release Rate 16. SECURITY...polymer FB Fully balanced laminate FEA Finite - element analysis FTCM Fracture toughness conversion mechanism G Shear modulus GI, GII, GIII Mode

A simple finite element method for linear hyperbolic problems

DOE PAGES

Mu, Lin; Ye, Xiu

2017-09-14

Here, we introduce a simple finite element method for solving first order hyperbolic equations with easy implementation and analysis. Our new method, with a symmetric, positive definite system, is designed to use discontinuous approximations on finite element partitions consisting of arbitrary shape of polygons/polyhedra. Error estimate is established. Extensive numerical examples are tested that demonstrate the robustness and flexibility of the method.

Reference Models for Multi-Layer Tissue Structures

DTIC Science & Technology

2016-09-01

simulation,  finite   element  analysis 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT 18. NUMBER OF PAGES 19a. NAME OF RESPONSIBLE PERSON USAMRMC...Physiologically realistic, fully specimen-specific, nonlinear reference models. Tasks. Finite element analysis of non-linear mechanics of cadaver...models. Tasks. Finite element analysis of non-linear mechanics of multi-layer tissue regions of human subjects. Deliverables. Partially subject- and

A simple finite element method for linear hyperbolic problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Ye, Xiu

Here, we introduce a simple finite element method for solving first order hyperbolic equations with easy implementation and analysis. Our new method, with a symmetric, positive definite system, is designed to use discontinuous approximations on finite element partitions consisting of arbitrary shape of polygons/polyhedra. Error estimate is established. Extensive numerical examples are tested that demonstrate the robustness and flexibility of the method.

Component mode synthesis and large deflection vibration of complex structures. Volume 3: Multiple-mode nonlinear free and forced vibrations of beams using finite element method

NASA Technical Reports Server (NTRS)

Mei, Chuh; Shen, Mo-How

1987-01-01

Multiple-mode nonlinear forced vibration of a beam was analyzed by the finite element method. Inplane (longitudinal) displacement and inertia (IDI) are considered in the formulation. By combining the finite element method and nonlinear theory, more realistic models of structural response are obtained more easily and faster.

Examples of finite element mesh generation using SDRC IDEAS

NASA Technical Reports Server (NTRS)

Zapp, John; Volakis, John L.

1990-01-01

IDEAS (Integrated Design Engineering Analysis Software) offers a comprehensive package for mechanical design engineers. Due to its multifaceted capabilities, however, it can be manipulated to serve the needs of electrical engineers, also. IDEAS can be used to perform the following tasks: system modeling, system assembly, kinematics, finite element pre/post processing, finite element solution, system dynamics, drafting, test data analysis, and project relational database.

Finite element analysis of transonic flows in cascades: Importance of computational grids in improving accuracy and convergence

NASA Technical Reports Server (NTRS)

Ecer, A.; Akay, H. U.

1981-01-01

The finite element method is applied for the solution of transonic potential flows through a cascade of airfoils. Convergence characteristics of the solution scheme are discussed. Accuracy of the numerical solutions is investigated for various flow regions in the transonic flow configuration. The design of an efficient finite element computational grid is discussed for improving accuracy and convergence.

AutoCAD-To-NASTRAN Translator Program

NASA Technical Reports Server (NTRS)

Jones, A.

1989-01-01

Program facilitates creation of finite-element mathematical models from geometric entities. AutoCAD to NASTRAN translator (ACTON) computer program developed to facilitate quick generation of small finite-element mathematical models for use with NASTRAN finite-element modeling program. Reads geometric data of drawing from Data Exchange File (DXF) used in AutoCAD and other PC-based drafting programs. Written in Microsoft Quick-Basic (Version 2.0).

The Constraint Method for Solid Finite Elements.

DTIC Science & Technology

1980-09-30

9. ’Hierarchical Approximation in Finite Element Analysis", by I. Norman Katz, International Symposium on Innovative Numerical Analysis In Applied ... Engineering Science, Versailles, France, May 23-27, 1977. 10. "Efficient Generation of Hierarchal Finite Elamnts Through the Use of Precomputed Arrays

Adaptive Shape Functions and Internal Mesh Adaptation for Modelling Progressive Failure in Adhesively Bonded Joints

NASA Technical Reports Server (NTRS)

Stapleton, Scott; Gries, Thomas; Waas, Anthony M.; Pineda, Evan J.

2014-01-01

Enhanced finite elements are elements with an embedded analytical solution that can capture detailed local fields, enabling more efficient, mesh independent finite element analysis. The shape functions are determined based on the analytical model rather than prescribed. This method was applied to adhesively bonded joints to model joint behavior with one element through the thickness. This study demonstrates two methods of maintaining the fidelity of such elements during adhesive non-linearity and cracking without increasing the mesh needed for an accurate solution. The first method uses adaptive shape functions, where the shape functions are recalculated at each load step based on the softening of the adhesive. The second method is internal mesh adaption, where cracking of the adhesive within an element is captured by further discretizing the element internally to represent the partially cracked geometry. By keeping mesh adaptations within an element, a finer mesh can be used during the analysis without affecting the global finite element model mesh. Examples are shown which highlight when each method is most effective in reducing the number of elements needed to capture adhesive nonlinearity and cracking. These methods are validated against analogous finite element models utilizing cohesive zone elements.

An anisotropic elastoplasticity model implemented in FLAG

DOE Office of Scientific and Technical Information (OSTI.GOV)

Buechler, Miles Allen; Canfield, Thomas R.

2017-10-12

Many metals, including Tantalum and Zirconium, exhibit anisotropic elastoplastic behavior at the single crystal level, and if components are manufactured from these metals through forming processes the polycrystal (component) may also exhibit anisotropic elastoplastic behavior. This is because the forming can induce a preferential orientation of the crystals in the polycrystal. One example is a rolled plate of Uranium where the sti /strong orientation of the crystal (c-axis) tends to align itself perpendicular to the rolling direction. If loads are applied to this plate in di erent orientations the sti ness as well as the ow strength of the materialmore » will be greater in the through thickness direction than in other directions. To better accommodate simulations of such materials, an anisotropic elastoplasticity model has been implemented in FLAG. The model includes an anisotropic elastic stress model as well as an anisotropic plasticity model. The model could represent single crystals of any symmetry, though it should not be confused with a high- delity crystal plasticity model with multiple slip planes and evolutions. The model is most appropriate for homogenized polycrystalline materials. Elastic rotation of the material due to deformation is captured, so the anisotropic models are appropriate for arbitrary large rotations, but currently they do not account for signi cant change in material texture beyond the elastic rotation of the entire polycrystal.« less

A comparison of the two approaches of the theory of critical distances based on linear-elastic and elasto-plastic analyses

NASA Astrophysics Data System (ADS)

Terekhina, A. I.; Plekhov, O. A.; Kostina, A. A.; Susmel, L.

2017-06-01

The problem of determining the strength of engineering structures, considering the effects of the non-local fracture in the area of stress concentrators is a great scientific and industrial interest. This work is aimed on modification of the classical theory of critical distance that is known as a method of failure prediction based on linear-elastic analysis in case of elasto-plastic material behaviour to improve the accuracy of estimation of lifetime of notched components. Accounting plasticity has been implemented with the use of the Simplified Johnson-Cook model. Mechanical tests were carried out using a 300 kN electromechanical testing machine Shimadzu AG-X Plus. The cylindrical un-notched specimens and specimens with stress concentrators of titanium alloy Grade2 were tested under tensile loading with different grippers travel speed, which ensured several orders of strain rate. The results of elasto-plastic analyses of stress distributions near a wide variety of notches are presented. The results showed that the use of the modification of the TCD based on elasto-plastic analysis gives us estimates falling within an error interval of ±5-10%, that more accurate predictions than the linear elastic TCD solution. The use of an improved description of the stress-strain state at the notch tip allows introducing the critical distances as a material parameter.

CONSTRUCTION OF SCALAR AND VECTOR FINITE ELEMENT FAMILIES ON POLYGONAL AND POLYHEDRAL MESHES

PubMed Central

GILLETTE, ANDREW; RAND, ALEXANDER; BAJAJ, CHANDRAJIT

2016-01-01

We combine theoretical results from polytope domain meshing, generalized barycentric coordinates, and finite element exterior calculus to construct scalar- and vector-valued basis functions for conforming finite element methods on generic convex polytope meshes in dimensions 2 and 3. Our construction recovers well-known bases for the lowest order Nédélec, Raviart-Thomas, and Brezzi-Douglas-Marini elements on simplicial meshes and generalizes the notion of Whitney forms to non-simplicial convex polygons and polyhedra. We show that our basis functions lie in the correct function space with regards to global continuity and that they reproduce the requisite polynomial differential forms described by finite element exterior calculus. We present a method to count the number of basis functions required to ensure these two key properties. PMID:28077939

CONSTRUCTION OF SCALAR AND VECTOR FINITE ELEMENT FAMILIES ON POLYGONAL AND POLYHEDRAL MESHES.

PubMed

Gillette, Andrew; Rand, Alexander; Bajaj, Chandrajit

2016-10-01

We combine theoretical results from polytope domain meshing, generalized barycentric coordinates, and finite element exterior calculus to construct scalar- and vector-valued basis functions for conforming finite element methods on generic convex polytope meshes in dimensions 2 and 3. Our construction recovers well-known bases for the lowest order Nédélec, Raviart-Thomas, and Brezzi-Douglas-Marini elements on simplicial meshes and generalizes the notion of Whitney forms to non-simplicial convex polygons and polyhedra. We show that our basis functions lie in the correct function space with regards to global continuity and that they reproduce the requisite polynomial differential forms described by finite element exterior calculus. We present a method to count the number of basis functions required to ensure these two key properties.

Design sensitivity analysis with Applicon IFAD using the adjoint variable method

NASA Technical Reports Server (NTRS)

Frederick, Marjorie C.; Choi, Kyung K.

1984-01-01

A numerical method is presented to implement structural design sensitivity analysis using the versatility and convenience of existing finite element structural analysis program and the theoretical foundation in structural design sensitivity analysis. Conventional design variables, such as thickness and cross-sectional areas, are considered. Structural performance functionals considered include compliance, displacement, and stress. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. That is, design sensitivity analysis software does not have to be imbedded in an existing finite element code. The finite element structural analysis program used in the implementation presented is IFAD. Feasibility of the method is shown through analysis of several problems, including built-up structures. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of a finite difference perturbation.

Parallel processing in finite element structural analysis

NASA Technical Reports Server (NTRS)

Noor, Ahmed K.

1987-01-01

A brief review is made of the fundamental concepts and basic issues of parallel processing. Discussion focuses on parallel numerical algorithms, performance evaluation of machines and algorithms, and parallelism in finite element computations. A computational strategy is proposed for maximizing the degree of parallelism at different levels of the finite element analysis process including: 1) formulation level (through the use of mixed finite element models); 2) analysis level (through additive decomposition of the different arrays in the governing equations into the contributions to a symmetrized response plus correction terms); 3) numerical algorithm level (through the use of operator splitting techniques and application of iterative processes); and 4) implementation level (through the effective combination of vectorization, multitasking and microtasking, whenever available).

A Simplified Finite Element Simulation for Straightening Process of Thin-Walled Tube

NASA Astrophysics Data System (ADS)

Zhang, Ziqian; Yang, Huilin

2017-12-01

The finite element simulation is an effective way for the study of thin-walled tube in the two cross rolls straightening process. To determine the accurate radius of curvature of the roll profile more efficiently, a simplified finite element model based on the technical parameters of an actual two cross roll straightening machine, was developed to simulate the complex straightening process. Then a dynamic simulation was carried out using ANSYS LS-DYNA program. The result implied that the simplified finite element model was reasonable for simulate the two cross rolls straightening process, and can be obtained the radius of curvature of the roll profile with the tube’s straightness 2 mm/m.

Finite element modeling and analysis of tires

NASA Technical Reports Server (NTRS)

Noor, A. K.; Andersen, C. M.

1983-01-01

Predicting the response of tires under various loading conditions using finite element technology is addressed. Some of the recent advances in finite element technology which have high potential for application to tire modeling problems are reviewed. The analysis and modeling needs for tires are identified. Reduction methods for large-scale nonlinear analysis, with particular emphasis on treatment of combined loads, displacement-dependent and nonconservative loadings; development of simple and efficient mixed finite element models for shell analysis, identification of equivalent mixed and purely displacement models, and determination of the advantages of using mixed models; and effective computational models for large-rotation nonlinear problems, based on a total Lagrangian description of the deformation are included.

Application of finite element substructuring to composite micromechanics. M.S. Thesis - Akron Univ., May 1984

NASA Technical Reports Server (NTRS)

Caruso, J. J.

1984-01-01

Finite element substructuring is used to predict unidirectional fiber composite hygral (moisture), thermal, and mechanical properties. COSMIC NASTRAN and MSC/NASTRAN are used to perform the finite element analysis. The results obtained from the finite element model are compared with those obtained from the simplified composite micromechanics equations. A unidirectional composite structure made of boron/HM-epoxy, S-glass/IMHS-epoxy and AS/IMHS-epoxy are studied. The finite element analysis is performed using three dimensional isoparametric brick elements and two distinct models. The first model consists of a single cell (one fiber surrounded by matrix) to form a square. The second model uses the single cell and substructuring to form a nine cell square array. To compare computer time and results with the nine cell superelement model, another nine cell model is constructed using conventional mesh generation techniques. An independent computer program consisting of the simplified micromechanics equation is developed to predict the hygral, thermal, and mechanical properties for this comparison. The results indicate that advanced techniques can be used advantageously for fiber composite micromechanics.

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Wang, Junping; Ye, Xiu

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for bothmore » the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.« less

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

DOE PAGES

Mu, Lin; Wang, Junping; Ye, Xiu

2017-08-17

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for bothmore » the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.« less

Finite Element Analysis of Tube Hydroforming in Non-Symmetrical Dies

NASA Astrophysics Data System (ADS)

Nulkar, Abhishek V.; Gu, Randy; Murty, Pilaka

2011-08-01

Tube hydroforming has been studied intensively using commercial finite element programs. A great deal of the investigations dealt with models with symmetric cross-sections. It is known that additional constraints due to symmetry may be imposed on the model so that it is properly supported. For a non-symmetric model, these constraints become invalid and the model does not have sufficient support resulting in a singular finite element system. Majority of commercial codes have a limited capability in solving models with insufficient supports. Recently, new algorithms using penalty variable and air-like contact element (ALCE) have been developed to solve positive semi-definite finite element systems such as those in contact mechanics. In this study the ALCE algorithm is first validated by comparing its result against a commercial code using a symmetric model in which a circular tube is formed to polygonal dies with symmetric shapes. Then, the study investigates the accuracy and efficiency of using ALCE in analyzing hydroforming of tubes with various cross-sections in non-symmetrical dies in 2-D finite element settings.

Interpolation Hermite Polynomials For Finite Element Method

NASA Astrophysics Data System (ADS)

Gusev, Alexander; Vinitsky, Sergue; Chuluunbaatar, Ochbadrakh; Chuluunbaatar, Galmandakh; Gerdt, Vladimir; Derbov, Vladimir; Góźdź, Andrzej; Krassovitskiy, Pavel

2018-02-01

We describe a new algorithm for analytic calculation of high-order Hermite interpolation polynomials of the simplex and give their classification. A typical example of triangle element, to be built in high accuracy finite element schemes, is given.

Dynamic characterization, monitoring and control of rotating flexible beam-mass structures via piezo-embedded techniques

NASA Technical Reports Server (NTRS)

Lai, Steven H.-Y.

1992-01-01

A variational principle and a finite element discretization technique were used to derive the dynamic equations for a high speed rotating flexible beam-mass system embedded with piezo-electric materials. The dynamic equation thus obtained allows the development of finite element models which accommodate both the original structural element and the piezoelectric element. The solutions of finite element models provide system dynamics needed to design a sensing system. The characterization of gyroscopic effect and damping capacity of smart rotating devices are addressed. Several simulation examples are presented to validate the analytical solution.

Finite element stress analysis of the human left ventricle whose irregular shape is developed from single plane cineangiocardiogram

NASA Technical Reports Server (NTRS)

Ghista, D. N.; Hamid, M. S.

1977-01-01

The three-dimensional left ventricular chamber geometrical model is developed from single plane cineangiocardiogram. This left ventricular model is loaded by an internal pressure monitored by cardiac catheterization. The resulting stresses in the left ventricular model chamber's wall are determined by computerized finite element procedure. For the discretization of this left ventricular model structure, a 20-node, isoparametric finite element is employed. The analysis and formulation of the computerised procedure is presented in the paper, along with the detailed algorithms and computer programs. The procedure is applied to determine the stresses in a left ventricle at an instant, during systole. Next, a portion (represented by a finite element) of this left ventricular chamber is simulated as being infarcted by making its active-state modulus value equal to its passive-state value; the neighbouring elements are shown to relieve the 'infarcted' element of stress by themselves taking on more stress.

Least-squares finite element methods for compressible Euler equations

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Carey, G. F.

1990-01-01

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L-sq-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L-sq method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H1-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.

A triangular prism solid and shell interactive mapping element for electromagnetic sheet metal forming process

NASA Astrophysics Data System (ADS)

Cui, Xiangyang; Li, She; Feng, Hui; Li, Guangyao

2017-05-01

In this paper, a novel triangular prism solid and shell interactive mapping element is proposed to solve the coupled magnetic-mechanical formulation in electromagnetic sheet metal forming process. A linear six-node "Triprism" element is firstly proposed for transient eddy current analysis in electromagnetic field. In present "Triprism" element, shape functions are given explicitly, and a cell-wise gradient smoothing operation is used to obtain the gradient matrices without evaluating derivatives of shape functions. In mechanical field analysis, a shear locking free triangular shell element is employed in internal force computation, and a data mapping method is developed to transfer the Lorentz force on solid into the external forces suffered by shell structure for dynamic elasto-plasticity deformation analysis. Based on the deformed triangular shell structure, a "Triprism" element generation rule is established for updated electromagnetic analysis, which means inter-transformation of meshes between the coupled fields can be performed automatically. In addition, the dynamic moving mesh is adopted for air mesh updating based on the deformation of sheet metal. A benchmark problem is carried out for confirming the accuracy of the proposed "Triprism" element in predicting flux density in electromagnetic field. Solutions of several EMF problems obtained by present work are compared with experiment results and those of traditional method, which are showing excellent performances of present interactive mapping element.

Finite element area and line integral transforms for generalization of aperture function and geometry in Kirchhoff scalar diffraction theory

NASA Astrophysics Data System (ADS)

Kraus, Hal G.

1993-02-01

Two finite element-based methods for calculating Fresnel region and near-field region intensities resulting from diffraction of light by two-dimensional apertures are presented. The first is derived using the Kirchhoff area diffraction integral and the second is derived using a displaced vector potential to achieve a line integral transformation. The specific form of each of these formulations is presented for incident spherical waves and for Gaussian laser beams. The geometry of the two-dimensional diffracting aperture(s) is based on biquadratic isoparametric elements, which are used to define apertures of complex geometry. These elements are also used to build complex amplitude and phase functions across the aperture(s), which may be of continuous or discontinuous form. The finite element transform integrals are accurately and efficiently integrated numerically using Gaussian quadrature. The power of these methods is illustrated in several examples which include secondary obstructions, secondary spider supports, multiple mirror arrays, synthetic aperture arrays, apertures covered by screens, apodization, phase plates, and off-axis apertures. Typically, the finite element line integral transform results in significant gains in computational efficiency over the finite element Kirchhoff transform method, but is also subject to some loss in generality.

Dynamic and thermal response finite element models of multi-body space structural configurations

NASA Technical Reports Server (NTRS)

Edighoffer, Harold H.

1987-01-01

Presented is structural dynamics modeling of two multibody space structural configurations. The first configuration is a generic space station model of a cylindrical habitation module, two solar array panels, radiator panel, and central connecting tube. The second is a 15-m hoop-column antenna. Discussed is the special joint elimination sequence used for these large finite element models, so that eigenvalues could be extracted. The generic space station model aided test configuration design and analysis/test data correlation. The model consisted of six finite element models, one of each substructure and one of all substructures as a system. Static analysis and tests at the substructure level fine-tuned the finite element models. The 15-m hoop-column antenna is a truss column and structural ring interconnected with tension stabilizing cables. To the cables, pretensioned mesh membrane elements were attached to form four parabolic shaped antennae, one per quadrant. Imposing thermal preloads in the cables and mesh elements produced pretension in the finite element model. Thermal preload variation in the 96 control cables was adjusted to maintain antenna shape within the required tolerance and to give pointing accuracy.

The Model Experiments and Finite Element Analysis on Deformation and Failure by Excavation of Grounds in Foregoing-roof Method

NASA Astrophysics Data System (ADS)

Sotokoba, Yasumasa; Okajima, Kenji; Iida, Toshiaki; Tanaka, Tadatsugu

We propose the trenchless box culvert construction method to construct box culverts in small covering soil layers while keeping roads or tracks open. When we use this construction method, it is necessary to clarify deformation and shear failure by excavation of grounds. In order to investigate the soil behavior, model experiments and elasto-plactic finite element analysis were performed. In the model experiments, it was shown that the shear failure was developed from the end of the roof to the toe of the boundary surface. In the finite element analysis, a shear band effect was introduced. Comparing the observed shear bands in model experiments with computed maximum shear strain contours, it was found that the observed direction of the shear band could be simulated reasonably by the finite element analysis. We may say that the finite element method used in this study is useful tool for this construction method.

Models and finite element approximations for interacting nanosized piezoelectric bodies and acoustic medium

NASA Astrophysics Data System (ADS)

Nasedkin, A. V.

2017-01-01

This research presents the new size-dependent models of piezoelectric materials oriented to finite element applications. The proposed models include the facilities of taking into account different mechanisms of damping for mechanical and electric fields. The coupled models also incorporate the equations of the theory of acoustics for viscous fluids. In particular cases, these models permit to use the mode superposition method with full separation of the finite element systems into independent equations for the independent modes for transient and harmonic problems. The main boundary conditions were supplemented with the facilities of taking into account the coupled surface effects, allowing to explore the nanoscale piezoelectric materials in the framework of theories of continuous media with surface stresses and their generalizations. For the considered problems we have implemented the finite element technologies and various numerical algorithms to maintain a symmetrical structure of the finite element quasi-definite matrices (matrix structure for the problems with a saddle point).

Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mohamed, M. Shadi, E-mail: m.s.mohamed@durham.ac.uk; Seaid, Mohammed; Trevelyan, Jon

2013-10-15

We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach canmore » be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.« less

Evaluation of stress changes in the mandible with a fixed functional appliance: a finite element study.

PubMed

Chaudhry, Anshul; Sidhu, Maninder S; Chaudhary, Girish; Grover, Seema; Chaudhry, Nimisha; Kaushik, Ashutosh

2015-02-01

The aim of this study was to evaluate the effects of a fixed functional appliance (Forsus Fatigue Resistant Device; 3M Unitek, Monrovia, Calif) on the mandible with 3-dimensional finite element stress analysis. A 3-dimensional finite element model of the mandible was constructed from the images generated by cone-beam computed tomography of a patient undergoing fixed orthodontic treatment. The changes were studied with the finite element method, in the form of highest von Mises stress and maximum principal stress regions. More areas of stress were seen in the model of the mandible with the Forsus compared with the model of the mandible in the resting stage. This fixed functional appliance studied by finite element model analysis caused increases in the maximum principal stress and the von Mises stress in both the cortical bone and the condylar region of the mandible by more than 2 times. Copyright © 2015 American Association of Orthodontists. Published by Elsevier Inc. All rights reserved.

Determination of ankle external fixation stiffness by expedited interactive finite element analysis.

PubMed

Nielsen, Jonathan K; Saltzman, Charles L; Brown, Thomas D

2005-11-01

Interactive finite element analysis holds the potential to quickly and accurately determine the mechanical stiffness of alternative external fixator frame configurations. Using as an example Ilizarov distraction of the ankle, a finite element model and graphical user interface were developed that provided rapid, construct-specific information on fixation rigidity. After input of specific construct variables, the finite element software determined the resulting tibial displacement for a given configuration in typically 15s. The formulation was employed to investigate constructs used to treat end-stage arthritis, both in a parametric series and for five specific clinical distraction cases. Parametric testing of 15 individual variables revealed that tibial half-pins were much more effective than transfixion wires in limiting axial tibial displacement. Factors most strongly contributing to stiffening the construct included placing the tibia closer to the fixator rings, and mounting the pins to the rings at the nearest circumferential location to the bone. Benchtop mechanical validation results differed inappreciably from the finite element computations.

Multiphysics elastodynamic finite element analysis of space debris deorbit stability and efficiency by electrodynamic tethers

NASA Astrophysics Data System (ADS)

Li, Gangqiang; Zhu, Zheng H.; Ruel, Stephane; Meguid, S. A.

2017-08-01

This paper developed a new multiphysics finite element method for the elastodynamic analysis of space debris deorbit by a bare flexible electrodynamic tether. Orbital motion limited theory and dynamics of flexible electrodynamic tethers are discretized by the finite element method, where the motional electric field is variant along the tether and coupled with tether deflection and motion. Accordingly, the electrical current and potential bias profiles of tether are solved together with the tether dynamics by the nodal position finite element method. The newly proposed multiphysics finite element method is applied to analyze the deorbit dynamics of space debris by electrodynamic tethers with a two-stage energy control strategy to ensure an efficient and stable deorbit process. Numerical simulations are conducted to study the coupled effect between the motional electric field and the tether dynamics. The results reveal that the coupling effect has a significant influence on the tether stability and the deorbit performance. It cannot be ignored when the libration and deflection of the tether are significant.

OpenSeesPy: Python library for the OpenSees finite element framework

NASA Astrophysics Data System (ADS)

Zhu, Minjie; McKenna, Frank; Scott, Michael H.

2018-01-01

OpenSees, an open source finite element software framework, has been used broadly in the earthquake engineering community for simulating the seismic response of structural and geotechnical systems. The framework allows users to perform finite element analysis with a scripting language and for developers to create both serial and parallel finite element computer applications as interpreters. For the last 15 years, Tcl has been the primary scripting language to which the model building and analysis modules of OpenSees are linked. To provide users with different scripting language options, particularly Python, the OpenSees interpreter interface was refactored to provide multi-interpreter capabilities. This refactoring, resulting in the creation of OpenSeesPy as a Python module, is accomplished through an abstract interface for interpreter calls with concrete implementations for different scripting languages. Through this approach, users are able to develop applications that utilize the unique features of several scripting languages while taking advantage of advanced finite element analysis models and algorithms.

Extended Finite Element Method with Simplified Spherical Harmonics Approximation for the Forward Model of Optical Molecular Imaging

PubMed Central

Li, Wei; Yi, Huangjian; Zhang, Qitan; Chen, Duofang; Liang, Jimin

2012-01-01

An extended finite element method (XFEM) for the forward model of 3D optical molecular imaging is developed with simplified spherical harmonics approximation (SPN). In XFEM scheme of SPN equations, the signed distance function is employed to accurately represent the internal tissue boundary, and then it is used to construct the enriched basis function of the finite element scheme. Therefore, the finite element calculation can be carried out without the time-consuming internal boundary mesh generation. Moreover, the required overly fine mesh conforming to the complex tissue boundary which leads to excess time cost can be avoided. XFEM conveniences its application to tissues with complex internal structure and improves the computational efficiency. Phantom and digital mouse experiments were carried out to validate the efficiency of the proposed method. Compared with standard finite element method and classical Monte Carlo (MC) method, the validation results show the merits and potential of the XFEM for optical imaging. PMID:23227108

Extended finite element method with simplified spherical harmonics approximation for the forward model of optical molecular imaging.

PubMed

Li, Wei; Yi, Huangjian; Zhang, Qitan; Chen, Duofang; Liang, Jimin

2012-01-01

An extended finite element method (XFEM) for the forward model of 3D optical molecular imaging is developed with simplified spherical harmonics approximation (SP(N)). In XFEM scheme of SP(N) equations, the signed distance function is employed to accurately represent the internal tissue boundary, and then it is used to construct the enriched basis function of the finite element scheme. Therefore, the finite element calculation can be carried out without the time-consuming internal boundary mesh generation. Moreover, the required overly fine mesh conforming to the complex tissue boundary which leads to excess time cost can be avoided. XFEM conveniences its application to tissues with complex internal structure and improves the computational efficiency. Phantom and digital mouse experiments were carried out to validate the efficiency of the proposed method. Compared with standard finite element method and classical Monte Carlo (MC) method, the validation results show the merits and potential of the XFEM for optical imaging.

Optimum element density studies for finite-element thermal analysis of hypersonic aircraft structures

NASA Technical Reports Server (NTRS)

Ko, William L.; Olona, Timothy; Muramoto, Kyle M.

1990-01-01

Different finite element models previously set up for thermal analysis of the space shuttle orbiter structure are discussed and their shortcomings identified. Element density criteria are established for the finite element thermal modelings of space shuttle orbiter-type large, hypersonic aircraft structures. These criteria are based on rigorous studies on solution accuracies using different finite element models having different element densities set up for one cell of the orbiter wing. Also, a method for optimization of the transient thermal analysis computer central processing unit (CPU) time is discussed. Based on the newly established element density criteria, the orbiter wing midspan segment was modeled for the examination of thermal analysis solution accuracies and the extent of computation CPU time requirements. The results showed that the distributions of the structural temperatures and the thermal stresses obtained from this wing segment model were satisfactory and the computation CPU time was at the acceptable level. The studies offered the hope that modeling the large, hypersonic aircraft structures using high-density elements for transient thermal analysis is possible if a CPU optimization technique was used.

Asymptotic analysis of stability for prismatic solids under axial loads

NASA Astrophysics Data System (ADS)

Scherzinger, W.; Triantafyllidis, N.

1998-06-01

This work addresses the stability of axially loaded prismatic beams with any simply connected crosssection. The solids obey a general class of rate-independent constitutive laws, and can sustain finite strains in either compression or tension. The proposed method is based on multiple scale asymptotic analysis, and starts with the full Lagrangian formulation for the three-dimensional stability problem, where the boundary conditions are chosen to avoid the formation of boundary layers. The calculations proceed by taking the limit of the beam's slenderness parameter, ɛ (ɛ 2 ≡ area/length 2), going to zero, thus resulting in asymptotic expressions for the critical loads and modes. The analysis presents a consistent and unified treatment for both compressive (buckling) and tensile (necking) instabilities, and is carried out explicitly up to o( ɛ4) in each case. The present method circumvents the standard structural mechanics approach for the stability problem of beams which requires the choice of displacement and stress field approximations in order to construct a nonlinear beam theory. Moreover, this work provides a consistent way to calculate the effect of the beam's slenderness on the critical load and mode to any order of accuracy required. In contrast, engineering theories give accurately the lowest order terms ( O( ɛ2)—Euler load—in compression or O(1)—maximum load—in tension) but give only approximately the next higher order terms, with the exception of simple section geometries where exact stability results are available. The proposed method is used to calculate the critical loads and eigenmodes for bars of several different cross-sections (circular, square, cruciform and L-shaped). Elastic beams are considered in compression and elastoplastic beams are considered in tension. The O( ɛ2) and O( ɛ4) asymptotic results are compared to the exact finite element calculations for the corresponding three-dimensional prismatic solids. The O( ɛ4) results give significant improvement over the O( ɛ2) results, even for extremely stubby beams, and in particular for the case of cross-sections with commensurate dimensions.

Slices: A Scalable Partitioner for Finite Element Meshes

NASA Technical Reports Server (NTRS)

Ding, H. Q.; Ferraro, R. D.

1995-01-01

A parallel partitioner for partitioning unstructured finite element meshes on distributed memory architectures is developed. The element based partitioner can handle mixtures of different element types. All algorithms adopted in the partitioner are scalable, including a communication template for unpredictable incoming messages, as shown in actual timing measurements.

Development and validation of a subject-specific finite element model of the functional spinal unit to predict vertebral strength.

PubMed

Lee, Chu-Hee; Landham, Priyan R; Eastell, Richard; Adams, Michael A; Dolan, Patricia; Yang, Lang

2017-09-01

Finite element models of an isolated vertebral body cannot accurately predict compressive strength of the spinal column because, in life, compressive load is variably distributed across the vertebral body and neural arch. The purpose of this study was to develop and validate a patient-specific finite element model of a functional spinal unit, and then use the model to predict vertebral strength from medical images. A total of 16 cadaveric functional spinal units were scanned and then tested mechanically in bending and compression to generate a vertebral wedge fracture. Before testing, an image processing and finite element analysis framework (SpineVox-Pro), developed previously in MATLAB using ANSYS APDL, was used to generate a subject-specific finite element model with eight-node hexahedral elements. Transversely isotropic linear-elastic material properties were assigned to vertebrae, and simple homogeneous linear-elastic properties were assigned to the intervertebral disc. Forward bending loading conditions were applied to simulate manual handling. Results showed that vertebral strengths measured by experiment were positively correlated with strengths predicted by the functional spinal unit finite element model with von Mises or Drucker-Prager failure criteria ( R 2  = 0.80-0.87), with areal bone mineral density measured by dual-energy X-ray absorptiometry ( R 2  = 0.54) and with volumetric bone mineral density from quantitative computed tomography ( R 2  = 0.79). Large-displacement non-linear analyses on all specimens did not improve predictions. We conclude that subject-specific finite element models of a functional spinal unit have potential to estimate the vertebral strength better than bone mineral density alone.

Numerical Methods for 2-Dimensional Modeling

DTIC Science & Technology

1980-12-01

high-order finite element methods, and a multidimensional version of the method of lines, both utilizing an optimized stiff integrator for the time...integration. The finite element methods have proved disappointing, but the method of lines has provided an unexpectedly large gain in speed. Two...diffusion problems with the same number of unknowns (a 21 x 41 grid), solved by second-order finite element methods, took over seven minutes on the Cray-i

A Moving Discontinuous Galerkin Finite Element Method for Flows with Interfaces

DTIC Science & Technology

2017-12-07

Naval Research Laboratory Washington, DC 20375-5320 NRL/MR/6040--17-9765 A Moving Discontinuous Galerkin Finite Element Method for Flows with...guidance to revise the method to ensure such properties. Acknowledgements This work was sponsored by the Office of Naval Research through the Naval...18. NUMBER OF PAGES 17. LIMITATION OF ABSTRACT A Moving Discontinuous Galerkin Finite Element Method for Flows with Interfaces Andrew Corrigan, Andrew

Structure and conformational dynamics of scaffolded DNA origami nanoparticles

DTIC Science & Technology

2017-05-08

all-atom molecular dynamics and coarse-grained finite element modeling to DX-based nanoparticles to elucidate their fine-scale and global conforma... finite element (FE) modeling approach CanDo is also routinely used to predict the 3D equilibrium conformation of programmed DNA assemblies based on a...model with both experimental cryo-electron microscopy (cryo-EM) data and all-atom modeling. MATERIALS AND METHODS Lattice-free finite element model

Object Toolkit Version 4.3 User’s Manual

DTIC Science & Technology

2016-12-31

unlimited. (OPS-17-12855 dtd 19 Jan 2017) 13. SUPPLEMENTARY NOTES 14. ABSTRACT Object Toolkit is a finite - element model builder specifically designed for...INTRODUCTION 1 What Is Object Toolkit? Object Toolkit is a finite - element model builder specifically designed for creating representations of spacecraft...Nascap-2k and EPIC, the user is not required to purchase or learn expensive finite element generators to create system models. Second, Object Toolkit

Dynamic Fracture of Concrete. Part 1

DTIC Science & Technology

1990-02-14

unnotched) by Mindess and the Charpy type impact tests by Shah. In both cases, dynamic finite element modeling with the adjusted constitutive equavm for the...Mindess and the Charpy type impact tests by Shah. In both cases, dynamic finite element modeling with the adjusted constitutive equations for the...Modeling Shah’s Charpy Impact Tests ................ 190 Figure 7.20 Specimen Configuration and Finite Element Model for Concrete and Mortar Beam Impact

Numerical model of glulam beam delamination in dependence on cohesive strength

NASA Astrophysics Data System (ADS)

Kawecki, Bartosz; Podgórski, Jerzy

2018-01-01

This paper presents an attempt of using a finite element method for predicting delamination of a glue laminated timber beam through a cohesive layer. There were used cohesive finite elements, quadratic stress damage initiation criterion and mixed mode energy release rate failure model. Finite element damage was equal to its complete stiffness degradation. Timber material was considered to be an orthotropic with plastic behaviour after reaching bending limit.

United States Air Force Graduate Student Research Program. Program Management Report

DTIC Science & Technology

1988-12-01

PRELIMINARY STRUCTURAL DESIGN/OPTIMIZATION by Richard A. Swift ABSTRACT Finite element analysis for use in structural design has advanced to the point where...Plates Subjected Gregory Schoeppner to Low Velocity Impact *** Same Report as Prof. William Wolfe * 57 Finite Element Analysis for Preliminary Richard...and dynamic load conditions using both radial and bias- ply tires. A detailed three-dimensional finite - element model of the wheel was generated for

Biomechanical Evaluation of Different Fixation Methods for Mandibular Anterior Segmental Osteotomy Using Finite Element Analysis, Part Two: Superior Repositioning Surgery With Bone Allograft.

PubMed

Kilinç, Yeliz; Erkmen, Erkan; Kurt, Ahmet

2016-01-01

In this study, the biomechanical behavior of different fixation methods used to fix the mandibular anterior segment following various amounts of superior repositioning was evaluated by using Finite Element Analysis (FEA). The three-dimensional finite element models representing 3 and 5 mm superior repositioning were generated. The gap in between segments was assumed to be filled by block bone allograft and resignated to be in perfect contact with the mandible and segmented bone. Six different finite element models with 2 distinct mobilization rate including 3 different fixation configurations, double right L (DRL), double left L (DLL), or double I (DI) miniplates with monocortical screws, correspondingly were created. A comparative evaluation has been made under vertical, horizontal and oblique loads. The von Mises and principal maximum stress (Pmax) values were calculated by finite element solver programme. The first part of our ongoing Finite Element Analysis research has been addressed to the mechanical behavior of the same fixation configurations in nongrafted models. In comparison with the findings of the first part of the study, it was concluded that bone graft offers superior mechanical stability without any limitation of mobilization and less stress on the fixative appliances as well as in the bone.

A combined experimental and finite element approach to analyse the fretting mechanism of the head-stem taper junction in total hip replacement.

PubMed

Bitter, Thom; Khan, Imran; Marriott, Tim; Lovelady, Elaine; Verdonschot, Nico; Janssen, Dennis

2017-09-01

Fretting corrosion at the taper interface of modular hip implants has been implicated as a possible cause of implant failure. This study was set up to gain more insight in the taper mechanics that lead to fretting corrosion. The objectives of this study therefore were (1) to select experimental loading conditions to reproduce clinically relevant fretting corrosion features observed in retrieved components, (2) to develop a finite element model consistent with the fretting experiments and (3) to apply more complicated loading conditions of activities of daily living to the finite element model to study the taper mechanics. The experiments showed similar wear patterns on the taper surface as observed in retrievals. The finite element wear score based on Archard's law did not correlate well with the amount of material loss measured in the experiments. However, similar patterns were observed between the simulated micromotions and the experimental wear measurements. Although the finite element model could not be validated, the loading conditions based on activities of daily living demonstrate the importance of assembly load on the wear potential. These findings suggest that finite element models that do not incorporate geometry updates to account for wear loss may not be appropriate to predict wear volumes of taper connections.

Numerical evaluation of implantable hearing devices using a finite element model of human ear considering viscoelastic properties.

PubMed

Zhang, Jing; Tian, Jiabin; Ta, Na; Huang, Xinsheng; Rao, Zhushi

2016-08-01

Finite element method was employed in this study to analyze the change in performance of implantable hearing devices due to the consideration of soft tissues' viscoelasticity. An integrated finite element model of human ear including the external ear, middle ear and inner ear was first developed via reverse engineering and analyzed by acoustic-structure-fluid coupling. Viscoelastic properties of soft tissues in the middle ear were taken into consideration in this model. The model-derived dynamic responses including middle ear and cochlea functions showed a better agreement with experimental data at high frequencies above 3000 Hz than the Rayleigh-type damping. On this basis, a coupled finite element model consisting of the human ear and a piezoelectric actuator attached to the long process of incus was further constructed. Based on the electromechanical coupling analysis, equivalent sound pressure and power consumption of the actuator corresponding to viscoelasticity and Rayleigh damping were calculated using this model. The analytical results showed that the implant performance of the actuator evaluated using a finite element model considering viscoelastic properties gives a lower output above about 3 kHz than does Rayleigh damping model. Finite element model considering viscoelastic properties was more accurate to numerically evaluate implantable hearing devices. © IMechE 2016.

Improved Finite Element Modeling of the Turbofan Engine Inlet Radiation Problem

NASA Technical Reports Server (NTRS)

Roy, Indranil Danda; Eversman, Walter; Meyer, H. D.

1993-01-01

Improvements have been made in the finite element model of the acoustic radiated field from a turbofan engine inlet in the presence of a mean flow. The problem of acoustic radiation from a turbofan engine inlet is difficult to model numerically because of the large domain and high frequencies involved. A numerical model with conventional finite elements in the near field and wave envelope elements in the far field has been constructed. By employing an irrotational mean flow assumption, both the mean flow and the acoustic perturbation problem have been posed in an axisymmetric formulation in terms of the velocity potential; thereby minimizing computer storage and time requirements. The finite element mesh has been altered in search of an improved solution. The mean flow problem has been reformulated with new boundary conditions to make it theoretically rigorous. The sound source at the fan face has been modeled as a combination of positive and negative propagating duct eigenfunctions. Therefore, a finite element duct eigenvalue problem has been solved on the fan face and the resulting modal matrix has been used to implement a source boundary condition on the fan face in the acoustic radiation problem. In the post processing of the solution, the acoustic pressure has been evaluated at Gauss points inside the elements and the nodal pressure values have been interpolated from them. This has significantly improved the results. The effect of the geometric position of the transition circle between conventional finite elements and wave envelope elements has been studied and it has been found that the transition can be made nearer to the inlet than previously assumed.

Life assessment of structural components using inelastic finite element analyses

NASA Technical Reports Server (NTRS)

Arya, Vinod K.; Halford, Gary R.

1993-01-01

The need for enhanced and improved performance of structural components subject to severe cyclic thermal/mechanical loadings, such as in the aerospace industry, requires development of appropriate solution technologies involving time-dependent inelastic analyses. Such analyses are mandatory to predict local stress-strain response and to assess more accurately the cyclic life time of structural components. The NASA-Lewis Research Center is cognizant of this need. As a result of concerted efforts at Lewis during the last few years, several such finite element solution technologies (in conjunction with the finite element program MARC) were developed and successfully applied to numerous uniaxial and multiaxial problems. These solution technologies, although developed for use with MARC program, are general in nature and can easily be extended for adaptation with other finite element programs such as ABAQUS, ANSYS, etc. The description and results obtained from two such inelastic finite element solution technologies are presented. The first employs a classical (non-unified) creep-plasticity model. An application of this technology is presented for a hypersonic inlet cowl-lip problem. The second of these technologies uses a unified creep-plasticity model put forth by Freed. The structural component for which this finite element solution technology is illustrated, is a cylindrical rocket engine thrust chamber. The advantages of employing a viscoplastic model for nonlinear time-dependent structural analyses are demonstrated. The life analyses for cowl-lip and cylindrical thrust chambers are presented. These analyses are conducted by using the stress-strain response of these components obtained from the corresponding finite element analyses.

Slave finite elements for nonlinear analysis of engine structures, volume 1

NASA Technical Reports Server (NTRS)

Gellin, S.

1991-01-01

A 336 degrees of freedom slave finite element processing capability to analyze engine structures under severe thermomechanical loading is presented. Description of the theoretical development and demonstration of that element is presented in this volume.

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

NASA Technical Reports Server (NTRS)

Kim, Sang-Wook

1988-01-01

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

Integral finite element analysis of turntable bearing with flexible rings

NASA Astrophysics Data System (ADS)

Deng, Biao; Liu, Yunfei; Guo, Yuan; Tang, Shengjin; Su, Wenbin; Lei, Zhufeng; Wang, Pengcheng

2018-03-01

This paper suggests a method to calculate the internal load distribution and contact stress of the thrust angular contact ball turntable bearing by FEA. The influence of the stiffness of the bearing structure and the plastic deformation of contact area on the internal load distribution and contact stress of the bearing is considered. In this method, the load-deformation relationship of the rolling elements is determined by the finite element contact analysis of a single rolling element and the raceway. Based on this, the nonlinear contact between the rolling elements and the inner and outer ring raceways is same as a nonlinear compression spring and bearing integral finite element analysis model including support structure was established. The effects of structural deformation and plastic deformation on the built-in stress distribution of slewing bearing are investigated on basis of comparing the consequences of load distribution, inner and outer ring stress, contact stress and other finite element analysis results with the traditional bearing theory, which has guiding function for improving the design of slewing bearing.

Applications of discrete element method in modeling of grain postharvest operations

USDA-ARS?s Scientific Manuscript database

Grain kernels are finite and discrete materials. Although flowing grain can behave like a continuum fluid at times, the discontinuous behavior exhibited by grain kernels cannot be simulated solely with conventional continuum-based computer modeling such as finite-element or finite-difference methods...

Finite element methods for the biomechanics of soft hydrated tissues: nonlinear analysis and adaptive control of meshes.

PubMed

Spilker, R L; de Almeida, E S; Donzelli, P S

1992-01-01

This chapter addresses computationally demanding numerical formulations in the biomechanics of soft tissues. The theory of mixtures can be used to represent soft hydrated tissues in the human musculoskeletal system as a two-phase continuum consisting of an incompressible solid phase (collagen and proteoglycan) and an incompressible fluid phase (interstitial water). We first consider the finite deformation of soft hydrated tissues in which the solid phase is represented as hyperelastic. A finite element formulation of the governing nonlinear biphasic equations is presented based on a mixed-penalty approach and derived using the weighted residual method. Fluid and solid phase deformation, velocity, and pressure are interpolated within each element, and the pressure variables within each element are eliminated at the element level. A system of nonlinear, first-order differential equations in the fluid and solid phase deformation and velocity is obtained. In order to solve these equations, the contributions of the hyperelastic solid phase are incrementally linearized, a finite difference rule is introduced for temporal discretization, and an iterative scheme is adopted to achieve equilibrium at the end of each time increment. We demonstrate the accuracy and adequacy of the procedure using a six-node, isoparametric axisymmetric element, and we present an example problem for which independent numerical solution is available. Next, we present an automated, adaptive environment for the simulation of soft tissue continua in which the finite element analysis is coupled with automatic mesh generation, error indicators, and projection methods. Mesh generation and updating, including both refinement and coarsening, for the two-dimensional examples examined in this study are performed using the finite quadtree approach. The adaptive analysis is based on an error indicator which is the L2 norm of the difference between the finite element solution and a projected finite element solution. Total stress, calculated as the sum of the solid and fluid phase stresses, is used in the error indicator. To allow the finite difference algorithm to proceed in time using an updated mesh, solution values must be transferred to the new nodal locations. This rezoning is accomplished using a projected field for the primary variables. The accuracy and effectiveness of this adaptive finite element analysis is demonstrated using a linear, two-dimensional, axisymmetric problem corresponding to the indentation of a thin sheet of soft tissue. The method is shown to effectively capture the steep gradients and to produce solutions in good agreement with independent, converged, numerical solutions.

Physics-based animation of large-scale splashing liquids, elastoplastic solids, and model-reduced flow

NASA Astrophysics Data System (ADS)

Gerszewski, Daniel James

Physical simulation has become an essential tool in computer animation. As the use of visual effects increases, the need for simulating real-world materials increases. In this dissertation, we consider three problems in physics-based animation: large-scale splashing liquids, elastoplastic material simulation, and dimensionality reduction techniques for fluid simulation. Fluid simulation has been one of the greatest successes of physics-based animation, generating hundreds of research papers and a great many special effects over the last fifteen years. However, the animation of large-scale, splashing liquids remains challenging. We show that a novel combination of unilateral incompressibility, mass-full FLIP, and blurred boundaries is extremely well-suited to the animation of large-scale, violent, splashing liquids. Materials that incorporate both plastic and elastic deformations, also referred to as elastioplastic materials, are frequently encountered in everyday life. Methods for animating such common real-world materials are useful for effects practitioners and have been successfully employed in films. We describe a point-based method for animating elastoplastic materials. Our primary contribution is a simple method for computing the deformation gradient for each particle in the simulation. Given the deformation gradient, we can apply arbitrary constitutive models and compute the resulting elastic forces. Our method has two primary advantages: we do not store or compare to an initial rest configuration and we work directly with the deformation gradient. The first advantage avoids poor numerical conditioning and the second naturally leads to a multiplicative model of deformation appropriate for finite deformations. One of the most significant drawbacks of physics-based animation is that ever-higher fidelity leads to an explosion in the number of degrees of freedom. This problem leads us to the consideration of dimensionality reduction techniques. We present several enhancements to model-reduced fluid simulation that allow improved simulation bases and two-way solid-fluid coupling. Specifically, we present a basis enrichment scheme that allows us to combine data-driven or artistically derived bases with more general analytic bases derived from Laplacian Eigenfunctions. Additionally, we handle two-way solid-fluid coupling in a time-splitting fashion---we alternately timestep the fluid and rigid body simulators, while taking into account the effects of the fluid on the rigid bodies and vice versa. We employ the vortex panel method to handle solid-fluid coupling and use dynamic pressure to compute the effect of the fluid on rigid bodies. Taken together, these contributions have advanced the state-of-the art in physics-based animation and are practical enough to be used in production pipelines.

Contact mechanics at nanometric scale using nanoindentation technique for brittle and ductile materials.

PubMed

Roa, J J; Rayon, E; Morales, M; Segarra, M

2012-06-01

In the last years, Nanoindentation or Instrumented Indentation Technique has become a powerful tool to study the mechanical properties at micro/nanometric scale (commonly known as hardness, elastic modulus and the stress-strain curve). In this review, the different contact mechanisms (elastic and elasto-plastic) are discussed, the recent patents for each mechanism (elastic and elasto-plastic) are summarized in detail, and the basic equations employed to know the mechanical behaviour for brittle and ductile materials are described.

Analysis of random structure-acoustic interaction problems using coupled boundary element and finite element methods

NASA Technical Reports Server (NTRS)

Mei, Chuh; Pates, Carl S., III

1994-01-01

A coupled boundary element (BEM)-finite element (FEM) approach is presented to accurately model structure-acoustic interaction systems. The boundary element method is first applied to interior, two and three-dimensional acoustic domains with complex geometry configurations. Boundary element results are very accurate when compared with limited exact solutions. Structure-interaction problems are then analyzed with the coupled FEM-BEM method, where the finite element method models the structure and the boundary element method models the interior acoustic domain. The coupled analysis is compared with exact and experimental results for a simplistic model. Composite panels are analyzed and compared with isotropic results. The coupled method is then extended for random excitation. Random excitation results are compared with uncoupled results for isotropic and composite panels.

The GPRIME approach to finite element modeling

NASA Technical Reports Server (NTRS)

Wallace, D. R.; Mckee, J. H.; Hurwitz, M. M.

1983-01-01

GPRIME, an interactive modeling system, runs on the CDC 6000 computers and the DEC VAX 11/780 minicomputer. This system includes three components: (1) GPRIME, a user friendly geometric language and a processor to translate that language into geometric entities, (2) GGEN, an interactive data generator for 2-D models; and (3) SOLIDGEN, a 3-D solid modeling program. Each component has a computer user interface of an extensive command set. All of these programs make use of a comprehensive B-spline mathematics subroutine library, which can be used for a wide variety of interpolation problems and other geometric calculations. Many other user aids, such as automatic saving of the geometric and finite element data bases and hidden line removal, are available. This interactive finite element modeling capability can produce a complete finite element model, producing an output file of grid and element data.

Modelling of thick composites using a layerwise laminate theory

NASA Technical Reports Server (NTRS)

Robbins, D. H., Jr.; Reddy, J. N.

1993-01-01

The layerwise laminate theory of Reddy (1987) is used to develop a layerwise, two-dimensional, displacement-based, finite element model of laminated composite plates that assumes a piecewise continuous distribution of the tranverse strains through the laminate thickness. The resulting layerwise finite element model is capable of computing interlaminar stresses and other localized effects with the same level of accuracy as a conventional 3D finite element model. Although the total number of degrees of freedom are comparable in both models, the layerwise model maintains a 2D-type data structure that provides several advantages over a conventional 3D finite element model, e.g. simplified input data, ease of mesh alteration, and faster element stiffness matrix formulation. Two sample problems are provided to illustrate the accuracy of the present model in computing interlaminar stresses for laminates in bending and extension.

Finite element method for optimal guidance of an advanced launch vehicle

NASA Technical Reports Server (NTRS)

Hodges, Dewey H.; Bless, Robert R.; Calise, Anthony J.; Leung, Martin

1992-01-01

A temporal finite element based on a mixed form of Hamilton's weak principle is summarized for optimal control problems. The resulting weak Hamiltonian finite element method is extended to allow for discontinuities in the states and/or discontinuities in the system equations. An extension of the formulation to allow for control inequality constraints is also presented. The formulation does not require element quadrature, and it produces a sparse system of nonlinear algebraic equations. To evaluate its feasibility for real-time guidance applications, this approach is applied to the trajectory optimization of a four-state, two-stage model with inequality constraints for an advanced launch vehicle. Numerical results for this model are presented and compared to results from a multiple-shooting code. The results show the accuracy and computational efficiency of the finite element method.

On finite element implementation and computational techniques for constitutive modeling of high temperature composites

NASA Technical Reports Server (NTRS)

Saleeb, A. F.; Chang, T. Y. P.; Wilt, T.; Iskovitz, I.

1989-01-01

The research work performed during the past year on finite element implementation and computational techniques pertaining to high temperature composites is outlined. In the present research, two main issues are addressed: efficient geometric modeling of composite structures and expedient numerical integration techniques dealing with constitutive rate equations. In the first issue, mixed finite elements for modeling laminated plates and shells were examined in terms of numerical accuracy, locking property and computational efficiency. Element applications include (currently available) linearly elastic analysis and future extension to material nonlinearity for damage predictions and large deformations. On the material level, various integration methods to integrate nonlinear constitutive rate equations for finite element implementation were studied. These include explicit, implicit and automatic subincrementing schemes. In all cases, examples are included to illustrate the numerical characteristics of various methods that were considered.

A simple finite element method for non-divergence form elliptic equation

DOE PAGES

Mu, Lin; Ye, Xiu

2017-03-01

Here, we develop a simple finite element method for solving second order elliptic equations in non-divergence form by combining least squares concept with discontinuous approximations. This simple method has a symmetric and positive definite system and can be easily analyzed and implemented. We could have also used general meshes with polytopal element and hanging node in the method. We prove that our finite element solution approaches to the true solution when the mesh size approaches to zero. Numerical examples are tested that demonstrate the robustness and flexibility of the method.

Convergence rates for finite element problems with singularities. Part 1: Antiplane shear. [crack

NASA Technical Reports Server (NTRS)

Plunkett, R.

1980-01-01

The problem of a finite crack in an infinite medium under antiplane shear load is considered. It is shown that the nodal forces at the tip of the crack accurately gives the order of singularity, that n energy release methods can give the strength to better than 1 percent with element size 1/10 the crack length, and that nodal forces give a much better estimate of the stress field than do the elements themselves. The finite element formulation and the factoring of tridiagonal matrices are discussed.

SUTRA (Saturated-Unsaturated Transport). A Finite-Element Simulation Model for Saturated-Unsaturated, Fluid-Density-Dependent Ground-Water Flow with Energy Transport or Chemically-Reactive Single-Species Solute Transport.

DTIC Science & Technology

1984-12-30

as three dimensional, when the assumption is made that all SUTRA parameters and coefficients have a constant value in the third space direction. A...finite element. The type of element employed by SUTRA for two-dimensional simulation is a quadrilateral which has a finite thickness in the third ... space dimension. This type of a quad- rilateral element and a typical two-dimensional mesh is shown in Figure 3.1. - All twelve edges of the two

A simple finite element method for non-divergence form elliptic equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Ye, Xiu

Here, we develop a simple finite element method for solving second order elliptic equations in non-divergence form by combining least squares concept with discontinuous approximations. This simple method has a symmetric and positive definite system and can be easily analyzed and implemented. We could have also used general meshes with polytopal element and hanging node in the method. We prove that our finite element solution approaches to the true solution when the mesh size approaches to zero. Numerical examples are tested that demonstrate the robustness and flexibility of the method.

Domain decomposition methods for nonconforming finite element spaces of Lagrange-type

NASA Technical Reports Server (NTRS)

Cowsar, Lawrence C.

1993-01-01

In this article, we consider the application of three popular domain decomposition methods to Lagrange-type nonconforming finite element discretizations of scalar, self-adjoint, second order elliptic equations. The additive Schwarz method of Dryja and Widlund, the vertex space method of Smith, and the balancing method of Mandel applied to nonconforming elements are shown to converge at a rate no worse than their applications to the standard conforming piecewise linear Galerkin discretization. Essentially, the theory for the nonconforming elements is inherited from the existing theory for the conforming elements with only modest modification by constructing an isomorphism between the nonconforming finite element space and a space of continuous piecewise linear functions.

Development library of finite elements for computer-aided design system of reed sensors

NASA Astrophysics Data System (ADS)

Kozlov, A. S.; Shmakov, N. A.; Tkalich, V. L.; Labkovskaia, R. I.; Kalinkina, M. E.; Pirozhnikova, O. I.

2018-05-01

The article is devoted to the development of a modern highly reliable element base of devices for security and fire alarm systems, in particular, to the improvement of the quality of contact cores (reed and membrane) of reed sensors. Modeling of elastic sensitive elements uses quadrangular elements of plates and shells, considered in the system of curvilinear orthogonal coordinates. The developed mathematical models and the formed finite element library are designed for systems of automated design of reed switch detectors to create competitive devices alarms. The finite element library is used for the automated system production of reed switch detectors both in series production and in the implementation of individual orders.

A class of hybrid finite element methods for electromagnetics: A review

NASA Technical Reports Server (NTRS)

Volakis, J. L.; Chatterjee, A.; Gong, J.

1993-01-01

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

A comparative study of an ABC and an artificial absorber for truncating finite element meshes

NASA Technical Reports Server (NTRS)

Oezdemir, T.; Volakis, John L.

1993-01-01

The type of mesh termination used in the context of finite element formulations plays a major role on the efficiency and accuracy of the field solution. The performance of an absorbing boundary condition (ABC) and an artificial absorber (a new concept) for terminating the finite element mesh was evaluated. This analysis is done in connection with the problem of scattering by a finite slot array in a thick ground plane. The two approximate mesh truncation schemes are compared with the exact finite element-boundary integral (FEM-BI) method in terms of accuracy and efficiency. It is demonstrated that both approximate truncation schemes yield reasonably accurate results even when the mesh is extended only 0.3 wavelengths away from the array aperture. However, the artificial absorber termination method leads to a substantially more efficient solution. Moreover, it is shown that the FEM-BI method remains quite competitive with the FEM-artificial absorber method when the FFT is used for computing the matrix-vector products in the iterative solution algorithm. These conclusions are indeed surprising and of major importance in electromagnetic simulations based on the finite element method.

[Three-dimensional finite element study on the change of glossopharyngeum in patient with obstructive sleep apnea hypopnea syndrome during titrated mandible advancement].

PubMed

Yang, Suixing; Feng, Jing; Zhang, Zuo; Qu, Aili; Gong, Miao; Tang, Jie; Fan, Junheng; Li, Songqing; Zhao, Yanling

2013-04-01

To construct a three-dimensional finite element model of the upper airway and adjacent structure of an obstructive sleep apnea hypopnea syndrome (OSAHS) patient for biomechanical analysis. And to study the influence of glossopharyngeum of an OSAHS patient with three-dimensional finite element model during titrated mandible advancement. DICOM format image information of an OSAHS patient's upper airway was obtained by thin-section CT scanning and digital image processing were utilized to construct a three-dimensional finite element model by Mimics 10.0, Imageware 10.0 and Ansys software. The changes and the law of glossopharyngeum were observed by biomechanics and morphology after loading with titrated mandible advancement. A three-dimensional finite element model of the adjacent upper airway structure of OSAHS was established successfully. After loading, the transverse diameter of epiglottis tip of glossopharyngeum increased significantly, although the sagittal diameter decreased correspondingly. The principal stress was mainly distributed in anterior wall of the upper airway. The location of principal stress concentration did not change significantly with the increasing of distance. The stress of glossopharyngeum increased during titrated mandible advancement. A more precise three-dimensional finite model of upper airway and adjacent structure of an OSAHS patient is established and improved efficiency by Mimics, Imageware and Ansys software. The glossopharyngeum of finite element model of OSAHS is analyzed by titrated mandible advancement and can effectively show the relationship between mandible advancement and the glossopharyngeum.

Application of numerical methods to heat transfer and thermal stress analysis of aerospace vehicles

NASA Technical Reports Server (NTRS)

Wieting, A. R.

1979-01-01

The paper describes a thermal-structural design analysis study of a fuel-injection strut for a hydrogen-cooled scramjet engine for a supersonic transport, utilizing finite-element methodology. Applications of finite-element and finite-difference codes to the thermal-structural design-analysis of space transports and structures are discussed. The interaction between the thermal and structural analyses has led to development of finite-element thermal methodology to improve the integration between these two disciplines. The integrated thermal-structural analysis capability developed within the framework of a computer code is outlined.

Generalizing the TRAPRG and TRAPAX finite elements

NASA Technical Reports Server (NTRS)

Hurwitz, M. M.

1983-01-01

The NASTRAN TRAPRG and TRAPAX finite elements are very restrictive as to shape and grid point numbering. The elements must be trapezoidal with two sides parallel to the radial axis. In addition, the ordering of the grid points on the element connection card must follow strict rules. The paper describes the generalization of these elements so that these restrictions no longer apply.

Finite elements and fluid dynamics. [instability effects on solution of nonlinear equations

NASA Technical Reports Server (NTRS)

Fix, G.

1975-01-01

Difficulties concerning a use of the finite element method in the solution of the nonlinear equations of fluid dynamics are partly related to various 'hidden' instabilities which often arise in fluid calculations. The instabilities are typically due to boundary effects or nonlinearities. It is shown that in certain cases these instabilities can be avoided if certain conservation laws are satisfied, and that the latter are often intimately related to finite elements.

A weak Galerkin least-squares finite element method for div-curl systems

NASA Astrophysics Data System (ADS)

Li, Jichun; Ye, Xiu; Zhang, Shangyou

2018-06-01

In this paper, we introduce a weak Galerkin least-squares method for solving div-curl problem. This finite element method leads to a symmetric positive definite system and has the flexibility to work with general meshes such as hybrid mesh, polytopal mesh and mesh with hanging nodes. Error estimates of the finite element solution are derived. The numerical examples demonstrate the robustness and flexibility of the proposed method.

BUCKY instruction manual, version 3.3

NASA Technical Reports Server (NTRS)

Smith, James P.

1994-01-01

The computer program BUCKY is a p-version finite element package for the solution of structural problems. The current version of BUCKY solves the 2-D plane stress, 3-D plane stress plasticity, 3-D axisymmetric, Mindlin and Kirchoff plate bending, and buckling problems. The p-version of the finite element method is a highly accurate version of the traditional finite element method. Example cases are presented to show the accuracy and application of BUCKY.

Decomposition of the Seismic Source Using Numerical Simulations and Observations of Nuclear Explosions

DTIC Science & Technology

2017-05-31

SUBJECT TERMS nonlinear finite element calculations, nuclear explosion monitoring, topography 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT 18...3D North Korea calculations........ Figure 6. The CRAM 3D finite element outer grid (left) is rectangular......................... Figure 7. Stress...Figure 6. The CRAM 3D finite element outer grid (left) is rectangular. The inner grid (center) is shaped to match the shape of the explosion shock wave