Iterative methods for mixed finite element equations
NASA Technical Reports Server (NTRS)
Nakazawa, S.; Nagtegaal, J. C.; Zienkiewicz, O. C.
1985-01-01
Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.
Material nonlinear analysis via mixed-iterative finite element method
NASA Technical Reports Server (NTRS)
Sutjahjo, Edhi; Chamis, Christos C.
1992-01-01
The performance of elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors are tested using 4-node quadrilateral finite elements. The membrane result is excellent, which indicates the implementation of elastic-plastic mixed-iterative analysis is appropriate. On the other hand, further research to improve bending performance of the method seems to be warranted.
Kim, S.
1994-12-31
Parallel iterative procedures based on domain decomposition techniques are defined and analyzed for the numerical solution of wave propagation by finite element and finite difference methods. For finite element methods, in a Lagrangian framework, an efficient way for choosing the algorithm parameter as well as the algorithm convergence are indicated. Some heuristic arguments for finding the algorithm parameter for finite difference schemes are addressed. Numerical results are presented to indicate the effectiveness of the methods.
Iterative methods for elliptic finite element equations on general meshes
NASA Technical Reports Server (NTRS)
Nicolaides, R. A.; Choudhury, Shenaz
1986-01-01
Iterative methods for arbitrary mesh discretizations of elliptic partial differential equations are surveyed. The methods discussed are preconditioned conjugate gradients, algebraic multigrid, deflated conjugate gradients, an element-by-element techniques, and domain decomposition. Computational results are included.
Elastic-plastic mixed-iterative finite element analysis: Implementation and performance assessment
NASA Technical Reports Server (NTRS)
Sutjahjo, Edhi; Chamis, Christos C.
1993-01-01
An elastic-plastic algorithm based on Von Mises and associative flow criteria is implemented in MHOST-a mixed iterative finite element analysis computer program developed by NASA Lewis Research Center. The performance of the resulting elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors of 4-node quadrilateral shell finite elements are tested for elastic-plastic performance. Generally, the membrane results are excellent, indicating the implementation of elastic-plastic mixed-iterative analysis is appropriate.
An approach to probabilistic finite element analysis using a mixed-iterative formulation
NASA Technical Reports Server (NTRS)
Dias, J. B.; Nakazawa, S.
1988-01-01
An efficient algorithm for computing the response sensitivity of finite element problems based on a mixed-iterative formulation is proposed. This method does not involve explicit differentiation of the tangent stiffness array and can be used with formulations for which a consistent tangent stiffness is not readily available. The method has been successfully applied to probabilistic finite element analysis of problems using the proposed mixed formulation, and this exercise has provided valuable insights regarding the extension of the method to a more general class of problems to include material and geometric nonlinearities.
NASA Technical Reports Server (NTRS)
Achar, N. S.; Gaonkar, G. H.
1993-01-01
Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used, and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.
NASA Astrophysics Data System (ADS)
Kergrene, Kenan; Babuška, Ivo; Banerjee, Uday
2016-06-01
The Generalized Finite Element Method (GFEM) is an extension of the Finite Element Method (FEM), where the standard finite element space is augmented with a space of non-polynomial functions, called the enrichment space. The functions in the enrichment space mimic the local behavior of the unknown solution of the underlying variational problem. GFEM has been successfully applied to a wide range of problems. However, it often suffers from bad conditioning, i.e., its conditioning may not be robust with respect to the mesh and in fact, the conditioning could be much worse than that of the standard FEM. In this paper, we present a numerical study that shows that if the "angle" between the finite element space and the enrichment space is bounded away from 0, uniformly with respect to the mesh, then the GFEM is stable, i.e., the conditioning of GFEM is not worse than that of the standard FEM. A GFEM with this property is called a Stable GFEM (SGFEM). The last part of the paper is devoted to the derivation of a robust iterative solver exploiting this angle condition. It is shown that the required "wall-clock" time is greatly reduced compared to popular GFEMs used in the literature.
NASA Astrophysics Data System (ADS)
Lavery, N.; Taylor, C.
1999-07-01
Multigrid and iterative methods are used to reduce the solution time of the matrix equations which arise from the finite element (FE) discretisation of the time-independent equations of motion of the incompressible fluid in turbulent motion. Incompressible flow is solved by using the method of reduce interpolation for the pressure to satisfy the Brezzi-Babuska condition. The k-l model is used to complete the turbulence closure problem. The non-symmetric iterative matrix methods examined are the methods of least squares conjugate gradient (LSCG), biconjugate gradient (BCG), conjugate gradient squared (CGS), and the biconjugate gradient squared stabilised (BCGSTAB). The multigrid algorithm applied is based on the FAS algorithm of Brandt, and uses two and three levels of grids with a V-cycling schedule. These methods are all compared to the non-symmetric frontal solver. Copyright
An iterative finite-element collocation method for parabolic problems using domain decomposition
Curran, M.C.
1992-01-01
Advection-dominated flows occur widely in the transport of groundwater contaminants, the movements of fluids in enhanced oil recovery projects, and many other contexts. In numerical models of such flows, adaptive local grid refinement is a conceptually attractive approach for resolving the sharp fronts or layers that tend to characterize the solutions. However, this approach can be difficult to implement in practice. A domain decomposition method developed by Bramble, Ewing, Pasciak, and Schatz, known as the BEPS method, overcomes many of the difficulties. We demonstrate the applicability of the iterative BEPS ideas to finite-element collocation on trial spaces of piecewise Hermite bicubics. The resulting scheme allows one to refine selected parts of a spatial grid without destroying algebraic efficiencies associated with the original coarse grid. We apply the method to two dimensional time-dependent advection-diffusion problems.
An iterative finite-element collocation method for parabolic problems using domain decomposition
Curran, M.C.
1992-11-01
Advection-dominated flows occur widely in the transport of groundwater contaminants, the movements of fluids in enhanced oil recovery projects, and many other contexts. In numerical models of such flows, adaptive local grid refinement is a conceptually attractive approach for resolving the sharp fronts or layers that tend to characterize the solutions. However, this approach can be difficult to implement in practice. A domain decomposition method developed by Bramble, Ewing, Pasciak, and Schatz, known as the BEPS method, overcomes many of the difficulties. We demonstrate the applicability of the iterative BEPS ideas to finite-element collocation on trial spaces of piecewise Hermite bicubics. The resulting scheme allows one to refine selected parts of a spatial grid without destroying algebraic efficiencies associated with the original coarse grid. We apply the method to two dimensional time-dependent advection-diffusion problems.
A comparison of direct and iterative finite element inversion techniques in dynamic elastography
NASA Astrophysics Data System (ADS)
Honarvar, M.; Rohling, R.; Salcudean, S. E.
2016-04-01
As part of tissue elasticity imaging or elastography, an inverse problem needs to be solved to find the elasticity distribution from the measured displacements. The finite element method (FEM) is a common method for solving the inverse problem in dynamic elastography. This problem has been solved with both direct and iterative FEM schemes. Each of these methods has its own advantages and disadvantages which are examined in this paper. Choosing the data resolution and the excitation frequency are critical for achieving the best estimation of the tissue elasticity in FEM methods. In this paper we investigate the performance of both direct and iterative FEMs for different ranges of excitation frequency. A new form of iterative method is suggested here which requires a lower mesh density compared to the original form. Also two forms of the direct method are compared in this paper: one using the exact fit for derivatives calculation and the other using the least squares fit. We also perform a study on the spatial resolution of these methods using simulations. The comparison is also validated using a phantom experiment. The results suggest that the direct method with least squares fit is more robust to noise compared to other methods but has slightly lower resolution results. For example, for the homogenous region with 20 dB noise added to the data, the RMS error for the direct method with least squares fit is approximately half of the iterative method. It was observed that the ratio of voxel size to the wavelength should be within a specific range for the results to be reliable. For example for the direct method with least squares fit, for the case of 20 dB noise level, this ratio should be between 0.1 to 0.2. On balance, considering the much higher computational cost of the iterative method, the dependency of the iterative method on the initial guess, and the greater robustness of the direct method to noise, we suggest using the direct method with least squares fit for
Koteras, J.R.
1993-07-01
Tunnels buried deep within the earth constitute an important class geomechanics problems. Two numerical techniques used for the analysis of geomechanics problems, the finite element method and the boundary element method, have complementary characteristics for applications to problems of this type. The usefulness of combining these two methods for use as a geomechanics analysis tool has been recognized for some time, and a number of coupling techniques have been proposed. However, not all of them lend themselves to efficient computational implementations for large-scale problems. This report examines a coupling technique that can form the basis for an efficient analysis tool for large scale geomechanics problems through the use of an iterative equation solver.
NASA Technical Reports Server (NTRS)
Horowitz, S. J.; Sigman, R. K.; Zinn, B. T.
1982-01-01
A new iterative solution technique for predicting the sound field radiated from a turbofan inlet in steady flight is presented. The sound field is divided into two regions: the sound field within and near the inlet which is computed using the finite element method and the radiation field beyond the inlet which is calculated using an integral solution technique. A continuous solution is obtained by matching the finite element and integral solutions at the interface between the two regions. The applicability of the iterative technique is demonstrated by comparison of experimental results with the theoretical results for several different inlet configurations with and without flow. These examples show that good agreement between experiment and theory is obtained within five iterations.
Domain decomposition based iterative methods for nonlinear elliptic finite element problems
Cai, X.C.
1994-12-31
The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.
Cwik, T.; Jamnejad, V.; Zuffada, C.
1994-12-31
The usefulness of finite element modeling follows from the ability to accurately simulate the geometry and three-dimensional fields on the scale of a fraction of a wavelength. To make this modeling practical for engineering design, it is necessary to integrate the stages of geometry modeling and mesh generation, numerical solution of the fields-a stage heavily dependent on the efficient use of a sparse matrix equation solver, and display of field information. The stages of geometry modeling, mesh generation, and field display are commonly completed using commercially available software packages. Algorithms for the numerical solution of the fields need to be written for the specific class of problems considered. Interior problems, i.e. simulating fields in waveguides and cavities, have been successfully solved using finite element methods. Exterior problems, i.e. simulating fields scattered or radiated from structures, are more difficult to model because of the need to numerically truncate the finite element mesh. To practically compute a solution to exterior problems, the domain must be truncated at some finite surface where the Sommerfeld radiation condition is enforced, either approximately or exactly. Approximate methods attempt to truncate the mesh using only local field information at each grid point, whereas exact methods are global, needing information from the entire mesh boundary. In this work, a method that couples three-dimensional finite element (FE) solutions interior to the bounding surface, with an efficient integral equation (IE) solution that exactly enforces the Sommerfeld radiation condition is developed. The bounding surface is taken to be a surface of revolution (SOR) to greatly reduce computational expense in the IE portion of the modeling.
NASA Astrophysics Data System (ADS)
Cao, Yong; Chu, Yuchuan; He, Xiaoming; Lin, Tao
2015-01-01
Interface problems involving the non-homogeneous flux jump condition are critical for engineering designs in the magnetostatic/electrostatic field. In applications, such as plasma simulation, we often only know the total electric quantity on the surface of the object, not the charge density distribution on the surface which appears as the non-homogeneous flux jump condition in the usual interface problems considered in the literature for the magnetostatic/electrostatic field. Based on structured meshes independent of the interface, this article proposes an iterative method that employs both the immersed finite element (IFE) method with non-homogeneous flux jump conditions and the regular finite element method with ghost nodes introduced in the object to solve the 2D interface problem for the potential field according to the given total electric quantity on the surface of the object. Numerical experiments are provided to illustrate the accuracy and efficiency of the proposed method.
Ozgun, Ozlem Mittra, Raj; Kuzuoglu, Mustafa
2009-04-01
In this paper, we introduce a parallelized version of a novel, non-iterative domain decomposition algorithm, called Characteristic Basis Finite Element Method (CBFEM-MPI), for efficient solution of large-scale electromagnetic scattering problems, by utilizing a set of specially defined characteristic basis functions (CBFs). This approach is based on the decomposition of the computational domain into a number of non-overlapping subdomains wherein the CBFs are generated by employing a novel procedure, which differs from all those that have been used in the past. Clearly, the CBFs are obtained by calculating the fields radiated by a finite number of dipole-type sources, which are placed hypothetically along the boundary of the conducting object. The major advantages of the proposed technique are twofold: (i) it provides a substantial reduction in the matrix size, and thus, makes use of direct solvers efficiently and (ii) it enables the utilization of parallel processing techniques that considerably decrease the overall computation time. We illustrate the application of the proposed approach via several 3D electromagnetic scattering problems.
NASA Astrophysics Data System (ADS)
Sharafat, S.; El-Awady, J.; Liu, S.; Diegele, E.; Ghoniem, N. M.
2007-08-01
Large-scale finite element modeling (FEM) of the US Dual Coolant Lead Lithium ITER Test Blanket Module including damage evolution is under development. A comprehensive rate-theory based radiation damage creep deformation code was integrated with the ABACUS FEM code. The advantage of this approach is that time-dependent in-reactor deformations and radiation damage can now be directly coupled with 'material properties' of FEM analyses. The coupled FEM-Creep damage model successfully simulated the simultaneous microstructure and stress evolution in small tensile test-bar structures. Applying the integrated Creep/FEM code to large structures is still computationally prohibitive. Instead, for thermo-structural analysis of the DCLL TBM structure the integrated FEM-creep damage model was used to develop true stress-strain behavior of F82H ferritic steel. Based on this integrated damage evolution-FEM approach it is proposed to use large-scale FEM analysis to identify and isolate critical stress areas for follow up analysis using detailed and fully integrated creep-FEM approach.
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.
Automatic finite element generators
NASA Technical Reports Server (NTRS)
Wang, P. S.
1984-01-01
The design and implementation of a software system for generating finite elements and related computations are described. Exact symbolic computational techniques are employed to derive strain-displacement matrices and element stiffness matrices. Methods for dealing with the excessive growth of symbolic expressions are discussed. Automatic FORTRAN code generation is described with emphasis on improving the efficiency of the resultant code.
Energy Science and Technology Software Center (ESTSC)
2006-03-08
MAPVAR-KD is designed to transfer solution results from one finite element mesh to another. MAPVAR-KD draws heavily from the structure and coding of MERLIN II, but it employs a new finite element data base, EXODUS II, and offers enhanced speed and new capabilities not available in MERLIN II. In keeping with the MERLIN II documentation, the computational algorithms used in MAPVAR-KD are described. User instructions are presented. Example problems are included to demonstrate the operationmore » of the code and the effects of various input options. MAPVAR-KD is a modification of MAPVAR in which the search algorithm was replaced by a kd-tree-based search for better performance on large problems.« less
NASA Astrophysics Data System (ADS)
Leibs, Christopher A.
Efforts are currently being directed towards a fully implicit, electromagnetic, JFNK-based solver, motivating the necessity of developing a fluid-based, electromag- netic, preconditioning strategy. The two-fluid plasma (TFP) model is an ideal approximation to the kinetic Jacobian. The TFP model couples both an ion and an electron fluid with Maxwell's equations. The fluid equations consist of the conservation of momentum and number density. A Darwin approximation of Maxwell is used to eliminate light waves from the model in order to facilitate coupling to non-relativistic particle models. We analyze the TFP-Darwin system in the context of a stand-alone solver with consideration of preconditioning a kinetic-JFNK approach. The TFP-Darwin system is addressed numerically by use of nested iteration (NI) and a First-Order Systems Least Squares (FOSLS) discretization. An important goal of NI is to produce an approximation that is within the basis of attraction for Newton's method on a relatively coarse mesh and, thus, on all subsequent meshes. After scaling and modification, the TFP-Darwin model yields a nonlinear, first-order system of equa- tions whose Frechet derivative is shown to be uniformly H1-elliptic in a neighborhood of the exact solution. H1 ellipticity yields optimal finite element performance and lin- ear systems amenable to solution with Algebraic Multigrid (AMG). To efficiently focus computational resources, an adaptive mesh refinement scheme, based on the accuracy per computational cost, is leveraged. Numerical tests demonstrate the efficacy of the approach, yielding an approximate solution within discretization error in a relatively small number of computational work units.
Energy Science and Technology Software Center (ESTSC)
2005-05-07
CONEX is a code for joining sequentially in time multiple exodusll database files which all represent the same base mesh topology and geometry. It is used to create a single results or restart file from multiple results or restart files which typically arise as the result of multiple restarted analyses. CONEX is used to postprocess the results from a series of finite element analyses. It can join sequentially the data from multiple results databases intomore » a single database which makes it easier to postprocess the results data.« less
Energy Science and Technology Software Center (ESTSC)
2005-06-26
Exotxt is an analysis code that reads finite element results data stored in an exodusII file and generates a file in a structured text format. The text file can be edited or modified via a number of text formatting tools. Exotxt is used by analysis to translate data from the binary exodusII format into a structured text format which can then be edited or modified and then either translated back to exodusII format or tomore » another format.« less
Probabilistic fracture finite elements
NASA Technical Reports Server (NTRS)
Liu, W. K.; Belytschko, T.; Lua, Y. J.
1991-01-01
The Probabilistic Fracture Mechanics (PFM) is a promising method for estimating the fatigue life and inspection cycles for mechanical and structural components. The Probability Finite Element Method (PFEM), which is based on second moment analysis, has proved to be a promising, practical approach to handle problems with uncertainties. As the PFEM provides a powerful computational tool to determine first and second moment of random parameters, the second moment reliability method can be easily combined with PFEM to obtain measures of the reliability of the structural system. The method is also being applied to fatigue crack growth. Uncertainties in the material properties of advanced materials such as polycrystalline alloys, ceramics, and composites are commonly observed from experimental tests. This is mainly attributed to intrinsic microcracks, which are randomly distributed as a result of the applied load and the residual stress.
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan
1993-01-01
A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.
Finite element analysis of wrinkling membranes
NASA Technical Reports Server (NTRS)
Miller, R. K.; Hedgepeth, J. M.; Weingarten, V. I.; Das, P.; Kahyai, S.
1984-01-01
The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.
Finite element computational fluid mechanics
NASA Technical Reports Server (NTRS)
Baker, A. J.
1983-01-01
Finite element analysis as applied to the broad spectrum of computational fluid mechanics is analyzed. The finite element solution methodology is derived, developed, and applied directly to the differential equation systems governing classes of problems in fluid mechanics. The heat conduction equation is used to reveal the essence and elegance of finite element theory, including higher order accuracy and convergence. The algorithm is extended to the pervasive nonlinearity of the Navier-Stokes equations. A specific fluid mechanics problem class is analyzed with an even mix of theory and applications, including turbulence closure and the solution of turbulent flows.
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.
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.
Probabilistic Finite Element: Variational Theory
NASA Technical Reports Server (NTRS)
Belytschko, T.; Liu, W. K.
1985-01-01
The goal of this research is to provide techniques which are cost-effective and enable the engineer to evaluate the effect of uncertainties in complex finite element models. Embedding the probabilistic aspects in a variational formulation is a natural approach. In addition, a variational approach to probabilistic finite elements enables it to be incorporated within standard finite element methodologies. Therefore, once the procedures are developed, they can easily be adapted to existing general purpose programs. Furthermore, the variational basis for these methods enables them to be adapted to a wide variety of structural elements and to provide a consistent basis for incorporating probabilistic features in many aspects of the structural problem. Tasks concluded include the theoretical development of probabilistic variational equations for structural dynamics, the development of efficient numerical algorithms for probabilistic sensitivity displacement and stress analysis, and integration of methodologies into a pilot computer code.
The NESSUS finite element code
NASA Technical Reports Server (NTRS)
Dias, J. B.; Nagiegaal, J. C.; Nakazawa, S.
1987-01-01
The objective of this development is to provide a new analysis tool which integrates the structural modeling versatility of a modern finite element code with the latest advances in the area of probabilistic modeling and structural reliability. Version 2.0 of the NESSUS finite element code was released last February, and is currently being exercised on a set of problems which are representative of typical Space Shuttle Main Engine (SSME) applications. NESSUS 2.0 allows linear elastostatic and eigenvalue analysis of structures with uncertain geometry, material properties and boundary conditions, which are subjected to a random mechanical and thermal loading environment. The NESSUS finite element code is a key component in a broader software system consisting of five major modules. NESSUS/EXPERT is an expert system under development at Southwest Research Institute, with the objective of centralizing all component-specific knowledge useful for conducting probabilistic analysis of typical Space Shuttle Main Engine (SSME) components. NESSUS/FEM contains the finite element code used for the structural analysis and parameter sensitivity evaluation of these components. The task of parametrizing a finite element mesh in terms of the random variables present is facilitated with the use of the probabilistic data preprocessor in NESSUS/PRE. An external database file is used for managing the bulk of the data generated by NESSUS/FEM.
Mixed finite elements for the Richards' equation: linearization procedure
NASA Astrophysics Data System (ADS)
Pop, I. S.; Radu, F.; Knabner, P.
2004-07-01
We consider mixed finite element discretization for a class of degenerate parabolic problems including the Richards' equation. After regularization, time discretization is achieved by an Euler implicit scheme, while mixed finite elements are employed for the discretization in space. Based on the results obtained in (Radu et al. RANA Preprint 02-06, Eindhoven University of Technology, 2002), this paper considers a simple iterative scheme to solve the emerging nonlinear elliptic problems.
Updating finite element dynamic models using an element-by-element sensitivity methodology
NASA Astrophysics Data System (ADS)
Farhat, Charbel; Hemez, Francois M.
1993-09-01
A sensitivity-based methodology for improving the finite element model of a given structure using test modal data and a few sensors is presented. The proposed method searches for both the location and sources of the mass and stiffness errors and does not interfere with the theory behind the finite element model while correcting these errors. The updating algorithm is derived from the unconstrained minimization of the squared L sub 2 norms of the modal dynamic residuals via an iterative two-step staggered procedure. At each iteration, the measured mode shapes are first expanded assuming that the model is error free, then the model parameters are corrected assuming that the expanded mode shapes are exact. The numerical algorithm is implemented in an element-by-element fashion and is capable of 'zooming' on the detected error locations. Several simulation examples which demonstate the potential of the proposed methodology are discussed.
Nonlinear, finite deformation, finite element analysis
NASA Astrophysics Data System (ADS)
Nguyen, Nhung; Waas, Anthony M.
2016-06-01
The roles of the consistent Jacobian matrix and the material tangent moduli, which are used in nonlinear incremental finite deformation mechanics problems solved using the finite element method, are emphasized in this paper, and demonstrated using the commercial software ABAQUS standard. In doing so, the necessity for correctly employing user material subroutines to solve nonlinear problems involving large deformation and/or large rotation is clarified. Starting with the rate form of the principle of virtual work, the derivations of the material tangent moduli, the consistent Jacobian matrix, the stress/strain measures, and the objective stress rates are discussed and clarified. The difference between the consistent Jacobian matrix (which, in the ABAQUS UMAT user material subroutine is referred to as DDSDDE) and the material tangent moduli ( C e ) needed for the stress update is pointed out and emphasized in this paper. While the former is derived based on the Jaumann rate of the Kirchhoff stress, the latter is derived using the Jaumann rate of the Cauchy stress. Understanding the difference between these two objective stress rates is crucial for correctly implementing a constitutive model, especially a rate form constitutive relation, and for ensuring fast convergence. Specifically, the implementation requires the stresses to be updated correctly. For this, the strains must be computed directly from the deformation gradient and corresponding strain measure (for a total form model). Alternatively, the material tangent moduli derived from the corresponding Jaumann rate of the Cauchy stress of the constitutive relation (for a rate form model) should be used. Given that this requirement is satisfied, the consistent Jacobian matrix only influences the rate of convergence. Its derivation should be based on the Jaumann rate of the Kirchhoff stress to ensure fast convergence; however, the use of a different objective stress rate may also be possible. The error associated
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.
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)
Peridynamic Multiscale Finite Element Methods
Costa, Timothy; Bond, Stephen D.; Littlewood, David John; Moore, Stan Gerald
2015-12-01
The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the
Finite-element solutions for geothermal systems
NASA Technical Reports Server (NTRS)
Chen, J. C.; Conel, J. E.
1977-01-01
Vector potential and scalar potential are used to formulate the governing equations for a single-component and single-phase geothermal system. By assuming an initial temperature field, the fluid velocity can be determined which, in turn, is used to calculate the convective heat transfer. The energy equation is then solved by considering convected heat as a distributed source. Using the resulting temperature to compute new source terms, the final results are obtained by iterations of the procedure. Finite-element methods are proposed for modeling of realistic geothermal systems; the advantages of such methods are discussed. The developed methodology is then applied to a sample problem. Favorable agreement is obtained by comparisons with a previous study.
Finite element methods in numerical relativity.
NASA Astrophysics Data System (ADS)
Mann, P. J.
The finite element method is very successful in Newtonian fluid simulations, and can be extended to relativitstic fluid flows. This paper describes the general method, and then outlines some preliminary results for spherically symmetric geometries. The mixed finite element - finite difference scheme is introduced, and used for the description of spherically symmetric collapse. Baker's (Newtonian) shock modelling method and Miller's moving finite element method are also mentioned. Collapse in double-null coordinates requires non-constant time slicing, so the full finite element method in space and time is described.
Finite element and finite difference methods in electromagnetic scattering
NASA Astrophysics Data System (ADS)
Morgan, Michael A.
Finite-difference and finite-element methods for the computational analysis of EM scattering phenomena are examined in chapters contributed by leading experts. Topics addressed include an FEM for composite scatterers, coupled finite- and boundary-element methods for EM scattering, absorbing boundary conditions for the direct solution PDEs arising in EM scattering problems, application of the control-region approximation to two-dimensional EM scattering, coupled potentials for EM fields in inhomogeneous media, the method of conforming boundary elements for transient electromagnetics, and the finite-difference time-domain method for numerical modeling of EM wave interactions with arbitrary structures. Extensive diagrams and graphs of typical results are provided.
Ferencz, R.M.
1988-10-01
Past work with element-by-element (EBE) preconditioned conjugate gradient iterative solution strategies has shown these techniques can be effective for large-scale, three-dimensional calculations in solid and structural mechanics. Significant gains over the profile storage direct solution method traditionally used in implicit finite element codes have been observed for a variety of real engineering analyses, especially in solid mechanics. Structural mechanics applications have proved less successful due to the ill-conditioned linear systems engendered by standard structural discretizations. This lack of robustness has recently motivated reconsideration of Lanczos-based algorithms as alternative iterative drivers. In this paper we compare the relative strengths of the conjugate gradient and Lanczos drivers when coupled with EBE preconditioning. The performance of the two methods is characterized, and compared with direct solution, using a model problem and a number of real engineering meshes. 26 refs., 14 figs., 2 tabs.
Fourier analysis of finite element preconditioned collocation schemes
NASA Technical Reports Server (NTRS)
Deville, Michel O.; Mund, Ernest H.
1990-01-01
The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.
User's Guide for ENSAERO_FE Parallel Finite Element Solver
NASA Technical Reports Server (NTRS)
Eldred, Lloyd B.; Guruswamy, Guru P.
1999-01-01
A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaces ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multizonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRANKOSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft model. It is used to create the stiffness matrices for each substructure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the substructure boundary nodes.
Probabilistic finite element analysis of a craniofacial finite element model.
Berthaume, Michael A; Dechow, Paul C; Iriarte-Diaz, Jose; Ross, Callum F; Strait, David S; Wang, Qian; Grosse, Ian R
2012-05-01
We employed a probabilistic finite element analysis (FEA) method to determine how variability in material property values affects stress and strain values in a finite model of a Macaca fascicularis cranium. The material behavior of cortical bone varied in three ways: isotropic homogeneous, isotropic non-homogeneous, and orthotropic non-homogeneous. The material behavior of the trabecular bone and teeth was always treated as isotropic and homogeneous. All material property values for the cranium were randomized with a Gaussian distribution with either coefficients of variation (CVs) of 0.2 or with CVs calculated from empirical data. Latin hypercube sampling was used to determine the values of the material properties used in the finite element models. In total, four hundred and twenty six separate deterministic FE simulations were executed. We tested four hypotheses in this study: (1) uncertainty in material property values will have an insignificant effect on high stresses and a significant effect on high strains for homogeneous isotropic models; (2) the effect of variability in material property values on the stress state will increase as non-homogeneity and anisotropy increase; (3) variation in the in vivo shear strain values reported by Strait et al. (2005) and Ross et al. (2011) is not only due to variations in muscle forces and cranial morphology, but also due to variation in material property values; (4) the assumption of a uniform coefficient of variation for the material property values will result in the same trend in how moderate-to-high stresses and moderate-to-high strains vary with respect to the degree of non-homogeneity and anisotropy as the trend found when the coefficients of variation for material property values are calculated from empirical data. Our results supported the first three hypotheses and falsified the fourth. When material properties were varied with a constant CV, as non-homogeneity and anisotropy increased the level of variability in
Domain decomposition methods for mortar finite elements
Widlund, O.
1996-12-31
In the last few years, domain decomposition methods, previously developed and tested for standard finite element methods and elliptic problems, have been extended and modified to work for mortar and other nonconforming finite element methods. A survey will be given of work carried out jointly with Yves Achdou, Mario Casarin, Maksymilian Dryja and Yvon Maday. Results on the p- and h-p-version finite elements will also be discussed.
FEBio: finite elements for biomechanics.
Maas, Steve A; Ellis, Benjamin J; Ateshian, Gerard A; Weiss, Jeffrey A
2012-01-01
In the field of computational biomechanics, investigators have primarily used commercial software that is neither geared toward biological applications nor sufficiently flexible to follow the latest developments in the field. This lack of a tailored software environment has hampered research progress, as well as dissemination of models and results. To address these issues, we developed the FEBio software suite (http://mrl.sci.utah.edu/software/febio), a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. This paper provides an overview of the theoretical basis of FEBio and its main features. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures. Software verification is a large part of the development and maintenance of FEBio, and to demonstrate the general approach, the description and results of several problems from the FEBio Verification Suite are presented and compared to analytical solutions or results from other established and verified FE codes. An additional simulation is described that illustrates the application of FEBio to a research problem in biomechanics. Together with the pre- and postprocessing software PREVIEW and POSTVIEW, FEBio provides a tailored solution for research and development in computational biomechanics. PMID:22482660
Finite element coiled cochlea model
NASA Astrophysics Data System (ADS)
Isailovic, Velibor; Nikolic, Milica; Milosevic, Zarko; Saveljic, Igor; Nikolic, Dalibor; Radovic, Milos; Filipović, Nenad
2015-12-01
Cochlea is important part of the hearing system, and thanks to special structure converts external sound waves into neural impulses which go to the brain. Shape of the cochlea is like snail, so geometry of the cochlea model is complex. The simplified cochlea coiled model was developed using finite element method inside SIFEM FP7 project. Software application is created on the way that user can prescribe set of the parameters for spiral cochlea, as well as material properties and boundary conditions to the model. Several mathematical models were tested. The acoustic wave equation for describing fluid in the cochlea chambers - scala vestibuli and scala timpani, and Newtonian dynamics for describing vibrations of the basilar membrane are used. The mechanical behavior of the coiled cochlea was analyzed and the third chamber, scala media, was not modeled because it does not have a significant impact on the mechanical vibrations of the basilar membrane. The obtained results are in good agreement with experimental measurements. Future work is needed for more realistic geometry model. Coiled model of the cochlea was created and results are compared with initial simplified coiled model of the cochlea.
Finite element neural networks for electromagnetic inverse problems
NASA Astrophysics Data System (ADS)
Ramuhalli, P.; Udpa, L.; Udpa, S.
2002-05-01
Iterative approaches using numerical forward models are commonly used for solving inverse problems in nondestructive evaluation. The drawbacks of these approaches include their high computational cost and the difficulty in computing gradients for updating defect profiles. This paper proposes a finite element neural network (FENN) that embeds finite element models into a neural network format. This approach enables fast and accurate solution of the forward problem. The FENN can then be used as the forward model in an iterative approach to solve the inverse problem. Gradient-based optimization methods are easily applied since the FENN provides an explicit functional mapping between the defect profile and the measured signal. Results of applying the FENN to several simple electromagnetic forward and inverse problems are presented.
Feng, Xiaobing
1996-12-31
A non-overlapping domain decomposition iterative method is proposed and analyzed for mixed finite element methods for a sequence of noncoercive elliptic systems with radiation boundary conditions. These differential systems describe the motion of a nearly elastic solid in the frequency domain. The convergence of the iterative procedure is demonstrated and the rate of convergence is derived for the case when the domain is decomposed into subdomains in which each subdomain consists of an individual element associated with the mixed finite elements. The hybridization of mixed finite element methods plays a important role in the construction of the discrete procedure.
Implicit extrapolation methods for multilevel finite element computations
Jung, M.; Ruede, U.
1994-12-31
The finite element package FEMGP has been developed to solve elliptic and parabolic problems arising in the computation of magnetic and thermomechanical fields. FEMGP implements various methods for the construction of hierarchical finite element meshes, a variety of efficient multilevel solvers, including multigrid and preconditioned conjugate gradient iterations, as well as pre- and post-processing software. Within FEMGP, multigrid {tau}-extrapolation can be employed to improve the finite element solution iteratively to higher order. This algorithm is based on an implicit extrapolation, so that the algorithm differs from a regular multigrid algorithm only by a slightly modified computation of the residuals on the finest mesh. Another advantage of this technique is, that in contrast to explicit extrapolation methods, it does not rely on the existence of global error expansions, and therefore neither requires uniform meshes nor global regularity assumptions. In the paper the authors will analyse the {tau}-extrapolation algorithm and present experimental results in the context of the FEMGP package. Furthermore, the {tau}-extrapolation results will be compared to higher order finite element solutions.
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.
Graphics for Finite-Element Analysis
NASA Technical Reports Server (NTRS)
Thornton, E. A.; Sawyer, L. M.
1982-01-01
ELPLOT program is a passive computer graphics system that could be utilized for display of models and responses of general finite-element analyses. Program includes: Wide range of view-orientation selections, number of alternative data-input formats, extensive family of finite-element types, and capabilities for both static and dynamic-response displays.
3-D Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
TAURUS is an interactive post-processing application supporting visualization of finite element analysis results on unstructured grids. TAURUS provides the ability to display deformed geometries and contours or fringes of a large number of derived results on meshes consisting of beam, plate, shell, and solid type finite elements. Time history plotting is also available.
Finite-Element Composite-Analysis Program
NASA Technical Reports Server (NTRS)
Bowles, David E.
1990-01-01
Finite Element Composite Analysis Program, FECAP, special-purpose finite-element program for analyzing behavior of composite material with microcomputer. Procedure leads to set of linear simultaneous equations relating unknown nodal displacement to applied loads. Written in HP BASIC 3.0.
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.
A finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.; Nayani, S.
1990-01-01
Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.
FDIPS: Finite Difference Iterative Potential-field Solver
NASA Astrophysics Data System (ADS)
Toth, Gabor; van der Holst, Bartholomeus; Huang, Zhenguang
2016-06-01
FDIPS is a finite difference iterative potential-field solver that can generate the 3D potential magnetic field solution based on a magnetogram. It is offered as an alternative to the spherical harmonics approach, as when the number of spherical harmonics is increased, using the raw magnetogram data given on a grid that is uniform in the sine of the latitude coordinate can result in inaccurate and unreliable results, especially in the polar regions close to the Sun. FDIPS is written in Fortran 90 and uses the MPI library for parallel execution.
Will Finite Elements Replace Structural Mechanics?
NASA Astrophysics Data System (ADS)
Ojalvo, I. U.
1984-01-01
This paper presents a personal view regarding the need for a continued interest and activity in structural methods in general, while viewing finite elements and the computer as simply two specific tools for assisting in this endeavor. An attempt is made to provide some insight as to why finite element methods seem to have "won the war," and to give examples of their more (and less) intelligent use. Items addressed include a highlight of unnecessary limitations of many existing standard finite element codes and where it is felt that further development work is needed.
The finite element method in thermomechanics
Hsu, T.
1986-01-01
Thermal stress analysis is critical in the design and operation of energy-efficient power plant components and engines as well as in nuclear and aerospace systems. The Finite Element Method in Thermomechanics attempts to embrace a wide range of topics in the nonlinear thermomechanical analysis. The book covers the basic principles of the finite element method: the formulations for the base thermomechanical analysis, including thermoelastic-plastic-creep stress analysis; the use of Fourier series for nonaxisymmetric loadings, and stress waves in solids in thermal environments; and the base finite element code called TEPSAC.
Assignment Of Finite Elements To Parallel Processors
NASA Technical Reports Server (NTRS)
Salama, Moktar A.; Flower, Jon W.; Otto, Steve W.
1990-01-01
Elements assigned approximately optimally to subdomains. Mapping algorithm based on simulated-annealing concept used to minimize approximate time required to perform finite-element computation on hypercube computer or other network of parallel data processors. Mapping algorithm needed when shape of domain complicated or otherwise not obvious what allocation of elements to subdomains minimizes cost of computation.
A sweeping preconditioner for time-harmonic Maxwell's equations with finite elements
NASA Astrophysics Data System (ADS)
Tsuji, Paul; Engquist, Bjorn; Ying, Lexing
2012-05-01
This paper is concerned with preconditioning the stiffness matrix resulting from finite element discretizations of Maxwell's equations in the high frequency regime. The moving PML sweeping preconditioner, first introduced for the Helmholtz equation on a Cartesian finite difference grid, is generalized to an unstructured mesh with finite elements. The method dramatically reduces the number of GMRES iterations necessary for convergence, resulting in an almost linear complexity solver. Numerical examples including electromagnetic cloaking simulations are presented to demonstrate the efficiency of the proposed method.
Visualization of higher order finite elements.
Thompson, David C.; Pebay, Philippe Pierre; Crawford, Richard H.; Khardekar, Rahul Vinay
2004-04-01
Finite element meshes are used to approximate the solution to some differential equation when no exact solution exists. A finite element mesh consists of many small (but finite, not infinitesimal or differential) regions of space that partition the problem domain, {Omega}. Each region, or element, or cell has an associated polynomial map, {Phi}, that converts the coordinates of any point, x = ( x y z ), in the element into another value, f(x), that is an approximate solution to the differential equation, as in Figure 1(a). This representation works quite well for axis-aligned regions of space, but when there are curved boundaries on the problem domain, {Omega}, it becomes algorithmically much more difficult to define {Phi} in terms of x. Rather, we define an archetypal element in a new coordinate space, r = ( r s t ), which has a simple, axis-aligned boundary (see Figure 1(b)) and place two maps onto our archetypal element:
A survey of mixed finite element methods
NASA Technical Reports Server (NTRS)
Brezzi, F.
1987-01-01
This paper is an introduction to and an overview of mixed finite element methods. It discusses the mixed formulation of certain basic problems in elasticity and hydrodynamics. It also discusses special techniques for solving the discrete problem.
Hybrid finite element-finite difference method for thermal analysis of blood vessels.
Blanchard, C H; Gutierrez, G; White, J A; Roemer, R B
2000-01-01
A hybrid finite-difference/finite-element technique for the thermal analysis of blood vessels embedded in perfused tissue has been developed and evaluated. This method provides efficient and accurate solutions to the conjugated heat transfer problem of convection by blood coupled to conduction in the tissue. The technique uses a previously developed 3D automatic meshing method for creating a finite element mesh in the tissue surrounding the vessels, coupled iteratively with a 1-D marching finite difference method for the interior of the vessels. This hybrid technique retains the flexibility and ease of automated finite-element meshing techniques for modelling the complex geometry of blood vessels and irregularly shaped tissues, and speeds the solution time by using a simple finite-difference method to calculate the bulk mean temperatures within all blood vessels. The use of the 1D finite-difference technique in the blood vessels also eliminates the large computer memory requirements needed to accurately solve large vessel network problems when fine FE meshes are used in the interior of vessels. The accuracy of the hybrid technique has been verified against previously verified numerical solutions. In summary, the hybrid technique combines the accuracy and flexibility found in automated finite-element techniques, with the speed and reduction of computational memory requirements associated with the 1D finite-difference technique, something which has not been done before. This method, thus, has the potential to provide accurate, flexible and relatively fast solutions for the thermal analysis of coupled perfusion/blood vessel problems, and large vessel network problems. PMID:10949130
Finite element modeling of the human pelvis
Carlson, B.
1995-11-01
A finite element model of the human pelvis was created using a commercial wire frame image as a template. To test the final mesh, the model`s mechanical behavior was analyzed through finite element analysis and the results were displayed graphically as stress concentrations. In the future, this grid of the pelvis will be integrated with a full leg model and used in side-impact car collision simulations.
Finite-Element Modeling For Structural Analysis
NASA Technical Reports Server (NTRS)
Min, J. B.; Androlake, S. G.
1995-01-01
Report presents study of finite-element mathematical modeling as used in analyzing stresses and strains at joints between thin, shell-like components (e.g., ducts) and thicker components (e.g., flanges or engine blocks). First approach uses global/local model to evaluate system. Provides correct total response and correct representation of stresses away from any discontinuities. Second approach involves development of special transition finite elements to model transitions between shells and thicker structural components.
Multigrid waveform relaxation on spatial finite element meshes
Janssen, J.; Vandewalle, S.
1994-12-31
The authors shall discuss the numerical solution of a parabolic partial differential equation {partial_derivative}u/{partial_derivative}t(x,t) = Lu(x,t) + f(x,t), x{element_of}{Omega}, t>0, (1) supplied with a boundary condition and given initial values. The spatial finite element discretization of (1) on a discrete grid {Omega}{sub h} leads to an initial value problem of the form B{dot u} + Au = f, u(0) = u{sub o}, t > 0, (2) with B a non-singular matrix. The waveform relaxation method is a method for solving ordinary differential equations. It differs from most standard iterative techniques in that it is a continuous-time method, iterating with functions in time, and thereby well-suited for parallel computation.
Finite element analysis of flexible, rotating blades
NASA Technical Reports Server (NTRS)
Mcgee, Oliver G.
1987-01-01
A reference guide that can be used when using the finite element method to approximate the static and dynamic behavior of flexible, rotating blades is given. Important parameters such as twist, sweep, camber, co-planar shell elements, centrifugal loads, and inertia properties are studied. Comparisons are made between NASTRAN elements through published benchmark tests. The main purpose is to summarize blade modeling strategies and to document capabilities and limitations (for flexible, rotating blades) of various NASTRAN elements.
Finite Element Interface to Linear Solvers
Williams, Alan
2005-03-18
Sparse systems of linear equations arise in many engineering applications, including finite elements, finite volumes, and others. The solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver capable of solving all of the linear systems that arise. This motivates the desire to switch an application from one solver librwy to another, depending on the problem being solved. The interfaces provided by solver libraries differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application Can be greatly reduced by having an abstraction layer between solver libraries and the application, putting a common "face" on various solver libraries. One such abstraction layer is the Finite Element Interface to Linear Solvers (EEl), which has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory.
Finite-Element Analysis of Multiphase Immiscible Flow Through Soils
NASA Astrophysics Data System (ADS)
Kuppusamy, T.; Sheng, J.; Parker, J. C.; Lenhard, R. J.
1987-04-01
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 equations governing flow in a three-fluid phase porous medium system with constant air phase pressure. Constitutive relationships for fluid conductivities and saturations as functions of fluid pressures, which are derived in a companion paper by J. C. Parker et al. (this issue) and which may be calibrated from two-phase laboratory measurements, are employed in the finite-element program. The solution procedure uses backward time integration with iteration by a modified Picard method to handle the nonlinear properties. Laboratory experiments involving water displacement from soil columns by p cymene (a benzene-derivative hydrocarbon) under constant pressure were simulated by the finite-element program to validate the numerical model and formulation for constitutive properties. Transient water outflow predicted using independently measured saturation-capillary head data agreed with observed outflow data within the limits of precision of the predictions as estimated by a first-order Taylor series approximation considering parameter uncertainty due to experimental reproducability and constitutive model accuracy. Two-dimensional simulations are presented for a hypothetical field case involving introduction of NAPL near the soil surface due to leakage from an underground storage tank. Subsequent transport of NAPL in the variably saturated vadose and groundwater zones is analyzed.
Quadrilateral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Benzley, Steven E
2012-10-16
Techniques for coarsening a quadrilateral mesh are described. These techniques include identifying a coarsening region within the quadrilateral mesh to be coarsened. Quadrilateral elements along a path through the coarsening region are removed. Node pairs along opposite sides of the path are identified. The node pairs along the path are then merged to collapse the path.
Finite-element models of continental extension
NASA Technical Reports Server (NTRS)
Lynch, H. David; Morgan, Paul
1990-01-01
Numerical models of the initial deformation of extending continental lithosphere, computed to investigate the control of preexisting thermal and mechanical heterogeneities on the style of deformation, are presented. The finite element method is used to calculate deformation with a viscoelastic-plastic model for the lithosphere. Comparisons of the results of analytic models and finite-element models using this method show that good results may be obtained by the numerical technique, even with elements containing both brittle and viscoelastic sampling points. It is shown that the gross style of initial extensional deformation is controlled by the depth and width of the initial heterogeneity which localizes deformation.
Finite element simulation of temperature dependent free surface flows
NASA Technical Reports Server (NTRS)
Engelman, M. S.; Sani, R. L.
1985-01-01
The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.
Surface subsidence prediction by nonlinear finite-element analysis
Najjar, Y. . Dept. of Civil Engineering); Zaman, M. . School of Civil Engineering and Environmental Science)
1993-11-01
An improved two-dimensional plane-strain numerical procedure based on the incremental-iterative nonlinear finite-element is developed to predict ground subsidence caused by underground mining. The procedure emphasizes the use of the following features: (1) an appropriate constitutive model that can accurately describe the nonlinear behavior of geological strata; and (2) an accurate algorithm for simulation of excavation sequences consistent with the actual underground mining process. The computer code is used to analyze a collapse that occurred in the Blue Goose Lease [number sign]1 Mine in northeastern Oklahoma. A parametric study is conducted to investigate the effects of some selected factors on the shape and extent of subsidence profiles. Analyses of the numerical results indicate that the nonlinear finite-element technique can be employed to meaningfully predict and characterize the potential for ground subsidence due to underground mining.
NASA Astrophysics Data System (ADS)
1990-10-01
Topics presented include the finite element analysis of confined turbulent swirling flows, compressible viscous flow calculations using compatible finite element approximations, the equilibrium and stability of Tokamaks, and a coupled finite element solution of biharmonic problems for vector potentials. Also presented are the Godunov-mixed methods for immiscible displacement, the iterative adaptive implicit-explicit methods for flow problems, finite element methods for one-dimensional combustion problems, and a technique for analyzing finite element methods for viscous incompressible flow.
Finite element model calibration using frequency responses with damping equalization
NASA Astrophysics Data System (ADS)
Abrahamsson, T. J. S.; Kammer, D. C.
2015-10-01
Model calibration is a cornerstone of the finite element verification and validation procedure, in which the credibility of the model is substantiated by positive comparison with test data. The calibration problem, in which the minimum deviation between finite element model data and experimental data is searched for, is normally characterized as being a large scale optimization problem with many model parameters to solve for and with deviation metrics that are nonlinear in these parameters. The calibrated parameters need to be found by iterative procedures, starting from initial estimates. Sometimes these procedures get trapped in local deviation function minima and do not converge to the globally optimal calibration solution that is searched for. The reason for such traps is often the multi-modality of the problem which causes eigenmode crossover problems in the iterative variation of parameter settings. This work presents a calibration formulation which gives a smooth deviation metric with a large radius of convergence to the global minimum. A damping equalization method is suggested to avoid the mode correlation and mode pairing problems that need to be solved in many other model updating procedures. By this method, the modal damping of a test data model and the finite element model is set to be the same fraction of critical modal damping. Mode pairing for mapping of experimentally found damping to the finite element model is thus not needed. The method is combined with model reduction for efficiency and employs the Levenberg-Marquardt minimizer with randomized starts to achieve the calibration solution. The performance of the calibration procedure, including a study of parameter bias and variance under noisy data conditions, is demonstrated by two numerical examples.
Massively parallel computation of RCS with finite elements
NASA Technical Reports Server (NTRS)
Parker, Jay
1993-01-01
One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary
Finite-element calculations on alliant FX/80
Watson, B.C.; Kamat, M.P.
1994-10-01
The finite-element method has proven to be an invaluable tool for analysis and design of complex, high-performance systems, such as those typically encountered in the aerospace or automotive industries. However, as the size of the finite-element models of such systems increases, analysis computation time using conventional computers can become prohibitively high. Parallel processing computers provide the means to overcome these computation-time limits, provided the algorithms used in the analysis can take advantage of multiple processors. The writers have examined several algorithms for linear and nonlinear static analysis, as well as dynamic finite-element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been investigated. For single load-case linear static analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear static analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms. For eigenvalue/eigenvector analysis, the subspace iteration algorithm with a parallel decomposition is shown to achieve a relatively high parallel efficiency. 12 refs.
Efficient finite element method for grating profile reconstruction
NASA Astrophysics Data System (ADS)
Zhang, Ruming; Sun, Jiguang
2015-12-01
This paper concerns the reconstruction of grating profiles from scattering data. The inverse problem is formulated as an optimization problem with a regularization term. We devise an efficient finite element method (FEM) and employ a quasi-Newton method to solve it. For the direct problems, the FEM stiff and mass matrices are assembled once at the beginning of the numerical procedure. Then only minor changes are made to the mass matrix at each iteration, which significantly saves the computation cost. Numerical examples show that the method is effective and robust.
Finite element analysis of periodic transonic flow problems
NASA Technical Reports Server (NTRS)
Fix, G. J.
1978-01-01
Flow about an oscillating thin airfoil in a transonic stream was considered. It was assumed that the flow field can be decomposed into a mean flow plus a periodic perturbation. On the surface of the airfoil the usual Neumman conditions are imposed. Two computer programs were written, both using linear basis functions over triangles for the finite element space. The first program uses a banded Gaussian elimination solver to solve the matrix problem, while the second uses an iterative technique, namely SOR. The only results obtained are for an oscillating flat plate.
Verification of Orthogrid Finite Element Modeling Techniques
NASA Technical Reports Server (NTRS)
Steeve, B. E.
1996-01-01
The stress analysis of orthogrid structures, specifically with I-beam sections, is regularly performed using finite elements. Various modeling techniques are often used to simplify the modeling process but still adequately capture the actual hardware behavior. The accuracy of such 'Oshort cutso' is sometimes in question. This report compares three modeling techniques to actual test results from a loaded orthogrid panel. The finite element models include a beam, shell, and mixed beam and shell element model. Results show that the shell element model performs the best, but that the simpler beam and beam and shell element models provide reasonable to conservative results for a stress analysis. When deflection and stiffness is critical, it is important to capture the effect of the orthogrid nodes in the model.
Visualizing higher order finite elements. Final report
Thompson, David C; Pebay, Philippe Pierre
2005-11-01
This report contains an algorithm for decomposing higher-order finite elements into regions appropriate for isosurfacing and proves the conditions under which the algorithm will terminate. Finite elements are used to create piecewise polynomial approximants to the solution of partial differential equations for which no analytical solution exists. These polynomials represent fields such as pressure, stress, and momentum. In the past, these polynomials have been linear in each parametric coordinate. Each polynomial coefficient must be uniquely determined by a simulation, and these coefficients are called degrees of freedom. When there are not enough degrees of freedom, simulations will typically fail to produce a valid approximation to the solution. Recent work has shown that increasing the number of degrees of freedom by increasing the order of the polynomial approximation (instead of increasing the number of finite elements, each of which has its own set of coefficients) can allow some types of simulations to produce a valid approximation with many fewer degrees of freedom than increasing the number of finite elements alone. However, once the simulation has determined the values of all the coefficients in a higher-order approximant, tools do not exist for visual inspection of the solution. This report focuses on a technique for the visual inspection of higher-order finite element simulation results based on decomposing each finite element into simplicial regions where existing visualization algorithms such as isosurfacing will work. The requirements of the isosurfacing algorithm are enumerated and related to the places where the partial derivatives of the polynomial become zero. The original isosurfacing algorithm is then applied to each of these regions in turn.
Finite element radiation transport in one dimension
Painter, J.F.
1997-05-09
A new physics package solves radiation transport equations in one space dimension, multiple energy groups and directions. A discontinuous finite element method discretizes radiation intensity with respect to space and angle, and a continuous finite element method discretizes electron temperature `in space. A splitting method solves the resulting linear equations. This is a one-dimensional analog of Kershaw and Harte`s two-dimensional package. This package has been installed in a two-dimensional inertial confinement fusion code, and has given excellent results for both thermal waves and highly directional radiation. In contrast, the traditional discrete ordinate and spherical harmonic methods show less accurate results in both cases.
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.
Slave finite elements: The temporal element approach to nonlinear analysis
NASA Technical Reports Server (NTRS)
Gellin, S.
1984-01-01
A formulation method for finite elements in space and time incorporating nonlinear geometric and material behavior is presented. The method uses interpolation polynomials for approximating the behavior of various quantities over the element domain, and only explicit integration over space and time. While applications are general, the plate and shell elements that are currently being programmed are appropriate to model turbine blades, vanes, and combustor liners.
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.
Quadrilateral/hexahedral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Scott, Michael A; Benzley, Steven E
2012-10-16
A technique for coarsening a finite element mesh ("FEM") is described. This technique includes identifying a coarsening region within the FEM to be coarsened. Perimeter chords running along perimeter boundaries of the coarsening region are identified. The perimeter chords are redirected to create an adaptive chord separating the coarsening region from a remainder of the FEM. The adaptive chord runs through mesh elements residing along the perimeter boundaries of the coarsening region. The adaptive chord is then extracted to coarsen the FEM.
Evolution of assumed stress hybrid finite element
NASA Technical Reports Server (NTRS)
Pian, T. H. H.
1984-01-01
Early versions of the assumed stress hybrid finite elements were based on the a priori satisifaction of stress equilibrium conditions. In the new version such conditions are relaxed but are introduced through additional internal displacement functions as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are now of complete polynomials up to the same degree as that of the assumed stresses. Several example problems indicate that optimal element properties are resulted by this method.
Finite element modeling of nonisothermal polymer flows
NASA Technical Reports Server (NTRS)
Roylance, D.
1981-01-01
A finite element formulation designed to simulate polymer melt flows in which both conductive and convective heat transfer are important is described, and the numerical model is illustrated by means of computer experiments using extruder drag flow and entry flow as trial problems. Fluid incompressibility is enforced by a penalty treatment of the element pressures, and the thermal convective transport is modeled by conventional Galerkin and optimal upwind treatments.
A multidimensional finite element method for CFD
NASA Technical Reports Server (NTRS)
Pepper, Darrell W.; Humphrey, Joseph W.
1991-01-01
A finite element method is used to solve the equations of motion for 2- and 3-D fluid flow. The time-dependent equations are solved explicitly using quadrilateral (2-D) and hexahedral (3-D) elements, mass lumping, and reduced integration. A Petrov-Galerkin technique is applied to the advection terms. The method requires a minimum of computational storage, executes quickly, and is scalable for execution on computer systems ranging from PCs to supercomputers.
Finite Element Simulation of Smart Structures
NASA Technical Reports Server (NTRS)
Cui, Y. Lawrence; Panahandeh, M.
1996-01-01
Finite element equations representing the behavior of piezoelectric materials when bounded to a typical structure and used as sensors and actuators were developed. Emphasis was placed on generating sensor output equations of piezoelectric sensors and responses of a typical structure bonded with piezoelectric sensors and actuators on the basis of finite element formulation. The model can predict not only structural responses due to both mechanical and electrical loading but also electrical potential due to mechanical or thermal effects. The resulted finite element equations were then used for simple control design and performance evaluation. In the control algorithm, voltages coming out from piezoelectric sensors, which are proportional to strains at sensing locations, are taken as input. The voltages applied to the piezoelectric actuators are used as output. The feasibility of integrating control algorithm with the element routine developed herein and FEAP was demonstrated. In particular, optimal independent modal space control was implemented in a software package on the basis of finite element formulation. A rudimentary finite element-control algorithm package was also developed to evaluate the performance of candidate control laws. A few numerical simulations using the software package developed herein were given. The integrated software package will provide a design tool to address issues such as how adaptive smart systems will scale to a full size aircraft, the amount of piezoelectric materials and the powers needed to actuate it for desired performance. It will also provide a viable new structural control design concept for practical applications in large flexible structures such as aerospace vehicles and aircraft.
Finite-approximation-error-based discrete-time iterative adaptive dynamic programming.
Wei, Qinglai; Wang, Fei-Yue; Liu, Derong; Yang, Xiong
2014-12-01
In this paper, a new iterative adaptive dynamic programming (ADP) algorithm is developed to solve optimal control problems for infinite horizon discrete-time nonlinear systems with finite approximation errors. First, a new generalized value iteration algorithm of ADP is developed to make the iterative performance index function converge to the solution of the Hamilton-Jacobi-Bellman equation. The generalized value iteration algorithm permits an arbitrary positive semi-definite function to initialize it, which overcomes the disadvantage of traditional value iteration algorithms. When the iterative control law and iterative performance index function in each iteration cannot accurately be obtained, for the first time a new "design method of the convergence criteria" for the finite-approximation-error-based generalized value iteration algorithm is established. A suitable approximation error can be designed adaptively to make the iterative performance index function converge to a finite neighborhood of the optimal performance index function. Neural networks are used to implement the iterative ADP algorithm. Finally, two simulation examples are given to illustrate the performance of the developed method. PMID:25265640
Finite element displacement analysis of a lung.
NASA Technical Reports Server (NTRS)
Matthews, F. L.; West, J. B.
1972-01-01
A method is given based on the technique of finite elements which determines theoretically the mechanical behavior of a lung-shaped body loaded by its own weight. The results of this theoretical analysis have been compared with actual measurements of alveolar size and pleural pressures in animal lungs.
Animation of finite element models and results
NASA Technical Reports Server (NTRS)
Lipman, Robert R.
1992-01-01
This is not intended as a complete review of computer hardware and software that can be used for animation of finite element models and results, but is instead a demonstration of the benefits of visualization using selected hardware and software. The role of raw computational power, graphics speed, and the use of videotape are discussed.
A triangular thin shell finite element: Nonlinear analysis. [structural analysis
NASA Technical Reports Server (NTRS)
Thomas, G. R.; Gallagher, R. H.
1975-01-01
Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.
Least-squares finite element methods for quantum chromodynamics
Ketelsen, Christian; Brannick, J; Manteuffel, T; Mccormick, S
2008-01-01
A significant amount of the computational time in large Monte Carlo simulations of lattice quantum chromodynamics (QCD) is spent inverting the discrete Dirac operator. Unfortunately, traditional covariant finite difference discretizations of the Dirac operator present serious challenges for standard iterative methods. For interesting physical parameters, the discretized operator is large and ill-conditioned, and has random coefficients. More recently, adaptive algebraic multigrid (AMG) methods have been shown to be effective preconditioners for Wilson's discretization of the Dirac equation. This paper presents an alternate discretization of the Dirac operator based on least-squares finite elements. The discretization is systematically developed and physical properties of the resulting matrix system are discussed. Finally, numerical experiments are presented that demonstrate the effectiveness of adaptive smoothed aggregation ({alpha}SA ) multigrid as a preconditioner for the discrete field equations resulting from applying the proposed least-squares FE formulation to a simplified test problem, the 2d Schwinger model of quantum electrodynamics.
A combined finite element-boundary element formulation for solution of axially symmetric bodies
NASA Technical Reports Server (NTRS)
Collins, Jeffrey D.; Volakis, John L.
1991-01-01
A new method is presented for the computation of electromagnetic scattering from axially symmetric bodies. To allow the simulation of inhomogeneous cross sections, the method combines the finite element and boundary element techniques. Interior to a fictitious surface enclosing the scattering body, the finite element method is used which results in a sparce submatrix, whereas along the enclosure the Stratton-Chu integral equation is enforced. By choosing the fictitious enclosure to be a right circular cylinder, most of the resulting boundary integrals are convolutional and may therefore be evaluated via the FFT with which the system is iteratively solved. In view of the sparce matrix associated with the interior fields, this reduces the storage requirement of the entire system to O(N) making the method attractive for large scale computations. The details of the corresponding formulation and its numerical implementation are described.
The mixed finite element multigrid method for stokes equations.
Muzhinji, K; Shateyi, S; Motsa, S S
2015-01-01
The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q2-Q1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results. PMID:25945361
The Mixed Finite Element Multigrid Method for Stokes Equations
Muzhinji, K.; Shateyi, S.; Motsa, S. S.
2015-01-01
The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q2-Q1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results. PMID:25945361
On Hybrid and mixed finite element methods
NASA Technical Reports Server (NTRS)
Pian, T. H. H.
1981-01-01
Three versions of the assumed stress hybrid model in finite element methods and the corresponding variational principles for the formulation are presented. Examples of rank deficiency for stiffness matrices by the hybrid stress model are given and their corresponding kinematic deformation modes are identified. A discussion of the derivation of general semi-Loof elements for plates and shells by the hybrid stress method is given. It is shown that the equilibrium model by Fraeijs de Veubeke can be derived by the approach of the hybrid stress model as a special case of semi-Loof elements.
Revolution in Orthodontics: Finite element analysis
Singh, Johar Rajvinder; Kambalyal, Prabhuraj; Jain, Megha; Khandelwal, Piyush
2016-01-01
Engineering has not only developed in the field of medicine but has also become quite established in the field of dentistry, especially Orthodontics. Finite element analysis (FEA) is a computational procedure to calculate the stress in an element, which performs a model solution. This structural analysis allows the determination of stress resulting from external force, pressure, thermal change, and other factors. This method is extremely useful for indicating mechanical aspects of biomaterials and human tissues that can hardly be measured in vivo. The results obtained can then be studied using visualization software within the finite element method (FEM) to view a variety of parameters, and to fully identify implications of the analysis. This is a review to show the applications of FEM in Orthodontics. It is extremely important to verify what the purpose of the study is in order to correctly apply FEM. PMID:27114948
Finite element computation with parallel VLSI
NASA Technical Reports Server (NTRS)
Mcgregor, J.; Salama, M.
1983-01-01
This paper describes a parallel processing computer consisting of a 16-bit microcomputer as a master processor which controls and coordinates the activities of 8086/8087 VLSI chip set slave processors working in parallel. The hardware is inexpensive and can be flexibly configured and programmed to perform various functions. This makes it a useful research tool for the development of, and experimentation with parallel mathematical algorithms. Application of the hardware to computational tasks involved in the finite element analysis method is demonstrated by the generation and assembly of beam finite element stiffness matrices. A number of possible schemes for the implementation of N-elements on N- or n-processors (N is greater than n) are described, and the speedup factors of their time consumption are determined as a function of the number of available parallel processors.
A finite element thermal analysis procedure for several temperature-dependent parameters
NASA Technical Reports Server (NTRS)
Thornton, E. A.; Wieting, A. R.
1978-01-01
A finite-element thermal analysis procedure for elements with several temperature-dependent thermal parameters is presented. The procedure, based on an application of the Newton-Raphson iteration technique, is formulated by resolving element matrices into component matrices, one component for each thermal parameter. Component conductance matrices are evaluated by assuming constant thermal parameters within an element and are computed once per unit thermal parameter. Significant savings in computer time result from the unit thermal parameter concept. The solution procedure applied to a convectively cooled structure with significantly varying thermal parameters converged in four iterations.
Finite Element Interface to Linear Solvers
Energy Science and Technology Software Center (ESTSC)
2005-03-18
Sparse systems of linear equations arise in many engineering applications, including finite elements, finite volumes, and others. The solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver capable of solving all of the linear systems that arise. This motivates the desire to switch an application from one solver librwy to another, depending on themore » problem being solved. The interfaces provided by solver libraries differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application Can be greatly reduced by having an abstraction layer between solver libraries and the application, putting a common "face" on various solver libraries. One such abstraction layer is the Finite Element Interface to Linear Solvers (EEl), which has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory.« less
Finite Element Heat & Mass Transfer Code
Energy Science and Technology Software Center (ESTSC)
1996-10-10
FEHM is a numerical simulation code for subsurface transport processes. It models 3-D, time-dependent, multiphase, multicomponent, non-isothermal, reactive flow through porous and fractured media. It can accurately represent complex 3-D geologic media and structures and their effects on subsurface flow and transport. Its capabilities include flow of gas, water, and heat; flow of air, water, and heat; multiple chemically reactive and sorbing tracers; finite element/finite volume formulation; coupled stress module; saturated and unsaturated media; andmore » double porosity and double porosity/double permeability capabilities.« less
Diagonal multisoliton matrix elements in finite volume
NASA Astrophysics Data System (ADS)
Pálmai, T.; Takács, G.
2013-02-01
We consider diagonal matrix elements of local operators between multisoliton states in finite volume in the sine-Gordon model and formulate a conjecture regarding their finite size dependence which is valid up to corrections exponential in the volume. This conjecture extends the results of Pozsgay and Takács which were only valid for diagonal scattering. In order to test the conjecture, we implement a numerical renormalization group improved truncated conformal space approach. The numerical comparisons confirm the conjecture, which is expected to be valid for general integrable field theories. The conjectured formula can be used to evaluate finite temperature one-point and two-point functions using recently developed methods.
Analysis of aircraft tires via semianalytic finite elements
NASA Technical Reports Server (NTRS)
Noor, Ahmed K.; Kim, Kyun O.; Tanner, John A.
1990-01-01
A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.
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.
Finite Element Analysis of Honeycomb Impact Attenuator
NASA Astrophysics Data System (ADS)
Yang, Seung-Yong; Choi, Seung-Kyu; Kim, Nohyu
To participate in Student Formula Society of Automotive Engineers (SAE) competitions, it is necessary to build an impact attenuator that would give an average deceleration not to exceed 20g when it runs into a rigid wall. Students can use numerical simulations or experimental test data to show that their car satisfies this safety requirement. A student group to study formula cars at the Korea University of Technology and Education has designed a vehicle to take part in a SAE competition, and a honeycomb structure was adopted as the impact attenuator. In this paper, finite element calculations were carried out to investigate the dynamic behavior of the honeycomb attenuator. Deceleration and deformation behaviors were studied. Effect of the yield strength was checked by comparing the numerical results. ABAQUS/Explicit finite element code was used.
Plasticity - Theory and finite element applications.
NASA Technical Reports Server (NTRS)
Armen, H., Jr.; Levine, H. S.
1972-01-01
A unified presentation is given of the development and distinctions associated with various incremental solution procedures used to solve the equations governing the nonlinear behavior of structures, and this is discussed within the framework of the finite-element method. Although the primary emphasis here is on material nonlinearities, consideration is also given to geometric nonlinearities acting separately or in combination with nonlinear material behavior. The methods discussed here are applicable to a broad spectrum of structures, ranging from simple beams to general three-dimensional bodies. The finite-element analysis methods for material nonlinearity are general in the sense that any of the available plasticity theories can be incorporated to treat strain hardening or ideally plastic behavior.
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmore » along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.« less
Finite element analysis of human joints
Bossart, P.L.; Hollerbach, K.
1996-09-01
Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described.
Finite element based electric motor design optimization
NASA Technical Reports Server (NTRS)
Campbell, C. Warren
1993-01-01
The purpose of this effort was to develop a finite element code for the analysis and design of permanent magnet electric motors. These motors would drive electromechanical actuators in advanced rocket engines. The actuators would control fuel valves and thrust vector control systems. Refurbishing the hydraulic systems of the Space Shuttle after each flight is costly and time consuming. Electromechanical actuators could replace hydraulics, improve system reliability, and reduce down time.
Finite Element Analysis of Reverberation Chambers
NASA Technical Reports Server (NTRS)
Bunting, Charles F.; Nguyen, Duc T.
2000-01-01
The primary motivating factor behind the initiation of this work was to provide a deterministic means of establishing the validity of the statistical methods that are recommended for the determination of fields that interact in -an avionics system. The application of finite element analysis to reverberation chambers is the initial step required to establish a reasonable course of inquiry in this particularly data-intensive study. The use of computational electromagnetics provides a high degree of control of the "experimental" parameters that can be utilized in a simulation of reverberating structures. As the work evolved there were four primary focus areas they are: 1. The eigenvalue problem for the source free problem. 2. The development of a complex efficient eigensolver. 3. The application of a source for the TE and TM fields for statistical characterization. 4. The examination of shielding effectiveness in a reverberating environment. One early purpose of this work was to establish the utility of finite element techniques in the development of an extended low frequency statistical model for reverberation phenomena. By employing finite element techniques, structures of arbitrary complexity can be analyzed due to the use of triangular shape functions in the spatial discretization. The effects of both frequency stirring and mechanical stirring are presented. It is suggested that for the low frequency operation the typical tuner size is inadequate to provide a sufficiently random field and that frequency stirring should be used. The results of the finite element analysis of the reverberation chamber illustrate io-W the potential utility of a 2D representation for enhancing the basic statistical characteristics of the chamber when operating in a low frequency regime. The basic field statistics are verified for frequency stirring over a wide range of frequencies. Mechanical stirring is shown to provide an effective frequency deviation.
ExodusII Finite Element Data Model
Energy Science and Technology Software Center (ESTSC)
2005-05-14
EXODUS II is a model developed to store and retrieve data for finite element analyses. It is used for preprocessing (problem definition), postprocessing (results visualization), as well as code to code data transfer. An EXODUS II data file is a random access, machine independent, binary file that is written and read via C, C++, or Fortran library routines which comprise the Application Programming Interface. (exodus II is based on netcdf)
Finite Element Results Visualization for Unstructured Grids
Speck, Douglas E.; Dovey, Donald J.
1996-07-15
GRIZ is a general-purpose post-processing application supporting interactive visualization of finite element analysis results on unstructured grids. In addition to basic pseudocolor renderings of state variables over the mesh surface, GRIZ provides modern visualization techniques such as isocontours and isosurfaces, cutting planes, vector field display, and particle traces. GRIZ accepts both command-line and mouse-driven input, and is portable to virtually any UNIX platform which provides Motif and OpenGl libraries.
Finite element model of needle electrode sensitivity
NASA Astrophysics Data System (ADS)
Høyum, P.; Kalvøy, H.; Martinsen, Ø. G.; Grimnes, S.
2010-04-01
We used the Finite Element (FE) Method to estimate the sensitivity of a needle electrode for bioimpedance measurement. This current conducting needle with insulated shaft was inserted in a saline solution and current was measured at the neutral electrode. FE model resistance and reactance were calculated and successfully compared with measurements on a laboratory model. The sensitivity field was described graphically based on these FE simulations.
FESDIF -- Finite Element Scalar Diffraction theory code
Kraus, H.G.
1992-09-01
This document describes the theory and use of a powerful scalar diffraction theory based computer code for calculation of intensity fields due to diffraction of optical waves by two-dimensional planar apertures and lenses. This code is called FESDIF (Finite Element Scalar Diffraction). It is based upon both Fraunhofer and Kirchhoff scalar diffraction theories. Simplified routines for circular apertures are included. However, the real power of the code comes from its basis in finite element methods. These methods allow the diffracting aperture to be virtually any geometric shape, including the various secondary aperture obstructions present in telescope systems. Aperture functions, with virtually any phase and amplitude variations, are allowed in the aperture openings. Step change aperture functions are accommodated. The incident waves are considered to be monochromatic. Plane waves, spherical waves, or Gaussian laser beams may be incident upon the apertures. Both area and line integral transformations were developed for the finite element based diffraction transformations. There is some loss of aperture function generality in the line integral transformations which are typically many times more computationally efficient than the area integral transformations when applicable to a particular problem.
Variational approach to probabilistic finite elements
NASA Technical Reports Server (NTRS)
Belytschko, T.; Liu, W. K.; Mani, A.; Besterfield, G.
1987-01-01
Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.
Enhancements to modal testing using finite elements
NASA Astrophysics Data System (ADS)
Jarvis, Brian
In calculating the natural frequencies and mode shapes from a finite element analysis, there are generally many more degrees of freedom than can be handled for the eigensolution. A reduction process is employed to reduce the number to a master set and chosen so that the modes of interest are well defined. By choosing those freedoms where the inertia terms are high or the stiffness terms are low then an automatic procedure for selecting the best freedoms can be defined. For modal testing, these master freedoms also indicate the best transducer locations for optimum low order mode identification. Having carried out the modal test, the mode shapes obtained can be forced onto the finite element model giving greatly enhanced results. By examining terms in all mode shapes from the finite element model in the frequency range of interest, the best reference or excitation position can be found. An example of the use of this technique to study the modal properties of an aero-engine compressor blade is given.
Variational approach to probabilistic finite elements
NASA Technical Reports Server (NTRS)
Belytschko, T.; Liu, W. K.; Mani, A.; Besterfield, G.
1991-01-01
Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.
Finite element-wavelet hybrid algorithm for atmospheric tomography.
Yudytskiy, Mykhaylo; Helin, Tapio; Ramlau, Ronny
2014-03-01
Reconstruction of the refractive index fluctuations in the atmosphere, or atmospheric tomography, is an underlying problem of many next generation adaptive optics (AO) systems, such as the multiconjugate adaptive optics or multiobject adaptive optics (MOAO). The dimension of the problem for the extremely large telescopes, such as the European Extremely Large Telescope (E-ELT), suggests the use of iterative schemes as an alternative to the matrix-vector multiply (MVM) methods. Recently, an algorithm based on the wavelet representation of the turbulence has been introduced in [Inverse Probl.29, 085003 (2013)] by the authors to solve the atmospheric tomography using the conjugate gradient iteration. The authors also developed an efficient frequency-dependent preconditioner for the wavelet method in a later work. In this paper we study the computational aspects of the wavelet algorithm. We introduce three new techniques, the dual domain discretization strategy, a scale-dependent preconditioner, and a ground layer multiscale method, to derive a method that is globally O(n), parallelizable, and compact with respect to memory. We present the computational cost estimates and compare the theoretical numerical performance of the resulting finite element-wavelet hybrid algorithm with the MVM. The quality of the method is evaluated in terms of an MOAO simulation for the E-ELT on the European Southern Observatory (ESO) end-to-end simulation system OCTOPUS. The method is compared to the ESO version of the Fractal Iterative Method [Proc. SPIE7736, 77360X (2010)] in terms of quality. PMID:24690653
Numerical performance of projection methods in finite element consolidation models
NASA Astrophysics Data System (ADS)
Gambolati, Giuseppe; Pini, Giorgio; Ferronato, Massimiliano
2001-12-01
Projection, or conjugate gradient like, methods are becoming increasingly popular for the efficient solution of large sparse sets of unsymmetric indefinite equations arising from the numerical integration of (initial) boundary value problems. One such problem is soil consolidation coupling a flow and a structural model, typically solved by finite elements (FE) in space and a marching scheme in time (e.g. the Crank-Nicolson scheme). The attraction of a projection method stems from a number of factors, including the ease of implementation, the requirement of limited core memory and the low computational cost if a cheap and effective matrix preconditioner is available. In the present paper, biconjugate gradient stabilized (Bi- CGSTAB) is used to solve FE consolidation equations in 2-D and 3-D settings with variable time integration steps. Three different nodal orderings are selected along with the preconditioner ILUT based on incomplete triangular factorization and variable fill-in. The overall cost of the solver is made up of the preconditioning cost plus the cost to converge which is in turn related to the number of iterations and the elementary operations required by each iteration. The results show that nodal ordering affects the perfor mance of Bi-CGSTAB. For normally conditioned consolidation problems Bi-CGSTAB with the best ILUT preconditioner may converge in a number of iterations up to two order of magnitude smaller than the size of the FE model and proves an accurate, cost-effective and robust alternative to direct methods.
Finite Elements approach for Density Functional Theory calculations on locally refined meshes
Fattebert, J; Hornung, R D; Wissink, A M
2006-03-27
We present a quadratic Finite Elements approach to discretize the Kohn-Sham equations on structured non-uniform meshes. A multigrid FAC preconditioner is proposed to iteratively solve the equations by an accelerated steepest descent scheme. The method was implemented using SAMRAI, a parallel software infrastructure for general AMR applications. Examples of applications to small nanoclusters calculations are presented.
Finite element procedures for coupled linear analysis of heat transfer, fluid and solid mechanics
NASA Technical Reports Server (NTRS)
Sutjahjo, Edhi; Chamis, Christos C.
1993-01-01
Coupled finite element formulations for fluid mechanics, heat transfer, and solid mechanics are derived from the conservation laws for energy, mass, and momentum. To model the physics of interactions among the participating disciplines, the linearized equations are coupled by combining domain and boundary coupling procedures. Iterative numerical solution strategy is presented to solve the equations, with the partitioning of temporal discretization implemented.
Finite Element approach for Density Functional Theory calculations on locally refined meshes
Fattebert, J; Hornung, R D; Wissink, A M
2007-02-23
We present a quadratic Finite Element approach to discretize the Kohn-Sham equations on structured non-uniform meshes. A multigrid FAC preconditioner is proposed to iteratively solve the equations by an accelerated steepest descent scheme. The method was implemented using SAMRAI, a parallel software infrastructure for general AMR applications. Examples of applications to small nanoclusters calculations are presented.
Regularised finite element model updating using measured incomplete modal data
NASA Astrophysics Data System (ADS)
Chen, Hua-Peng; Maung, Than Soe
2014-10-01
This paper presents an effective approach for directly updating finite element model from measured incomplete vibration modal data with regularised algorithms. The proposed method is based on the relationship between the perturbation of structural parameters such as stiffness change and the modal data measurements of the tested structure such as measured mode shape readings. In order to adjust structural parameters at detailed locations, structural updating parameters will be selected at critical point level to reflect the modelling errors at the connections of structural elements. These updating parameters are then evaluated by an iterative or a direct solution procedure, which gives optimised solutions in the least squares sense without requiring an optimisation technique. In order to reduce the influence of modal measurement uncertainty, the Tikhonov regularisation method incorporating the L-curve criterion is employed to produce reliable solutions for the chosen updating parameters. Numerical simulation investigations and experimental studies for the laboratory tested space steel frame structure are undertaken to verify the accuracy and effectiveness of the proposed methods for adjusting the stiffness at the joints of structural members. The results demonstrate that the proposed methods provide reliable estimates of finite element model updating using the measured incomplete modal data.
Mixed Finite Element Method for Melt Migration
NASA Astrophysics Data System (ADS)
Taicher, A. L.; Hesse, M. A.; Arbogast, T.
2012-12-01
Multi-phase flow arises during partial melting in the earth mantle, where the porosity is small and material has the characteristics of a compacting porous medium. The equations governing multi-phase flow have been specialized to partially molten materials by McKenzie and Fowler. Their model, also called a Darcy-Stokes system, is highly coupled and non-linear. Melt flow is governed by Darcy's Law while the high temperature, ductile creep of the solid matrix is modeled using viscous non-Newtonian Stokes rheology. In addition, the melt and solid pressures are related through a compaction relation. This nearly elliptic mechanical problem is then coupled with both solute transport and thermal evolution according to the enthalpy method developed by Katz. A suitable numerical method must solve the Darcy-Stokes problem in a manner compatible with the transport problem. Moreover, unlike most porous media problems, partially molten materials transition dynamically from non-porous solid to porous medium. Therefore, a numerical method must also carefully account for the limit of zero porosity. The Darcy-Stokes system for modeling partial melting in the mantle is a novel problem. As far as we know, there currently does not exist a finite element solution in the literature solving these coupled equations. The finite element framework provides support for additional analysis of error and convergence. Moreover, both mesh refinement and anisotropy are naturally incorporated into finite elements. In particular, the mixed finite element method presents a good candidate because it works in both limiting cases: Darcy and incompressible Stokes flow. Mixed methods also produce discretely conservative fluxes that are required for the transport problem to remains stable without violating conservation of mass. Based preliminary investigations in 1D and derived energy estimates, we present a mixed formulation for the Darcy-Stokes system. Next, using novel elements of lowest order and
Modelling bucket excavation by finite element
NASA Astrophysics Data System (ADS)
Pecingina, O. M.
2015-11-01
Changes in geological components of the layers from lignite pits have an impact on the sustainability of the cup path elements and under the action of excavation force appear efforts leading to deformation of the entire assembly. Application of finite element method in the optimization of components leads to economic growth, to increase the reliability and durability of the studied machine parts thus the machine. It is obvious usefulness of knowledge the state of mechanical tensions that the designed piece or the assembly not to break under the action of tensions that must cope during operation. In the course of excavation work on all bucket cutting force components, the first coming into contact with the material being excavated cutting edge. Therefore in the study with finite element analysis is retained only cutting edge. To study the field of stress and strain on the cutting edge will be created geometric patterns for each type of cup this will be subject to static analysis. The geometric design retains the cutting edge shape and on this on the tooth cassette location will apply an areal force on the abutment tooth. The cutting edge real pattern is subjected to finite element study for the worst case of rock cutting by symmetrical and asymmetrical cups whose profile is different. The purpose of this paper is to determine the displacement and tensions field for both profiles considering the maximum force applied on the cutting edge and the depth of the cutting is equal with the width of the cutting edge of the tooth. It will consider the worst case when on the structure will act both the tangential force and radial force on the bucket profile. For determination of stress and strain field on the form design of cutting edge profile will apply maximum force assuming uniform distribution and on the edge surface force will apply a radial force. After geometric patterns discretization on the cutting knives and determining stress field, can be seen that at the
The use of Galerkin finite-element methods to solve mass-transport equations
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)
Finite element modeling of pulsed eddy current NDT phenomena
Allen, B.; Ida, N.; Lord, W.
1985-05-15
Transient fields for nondestructive testing (pulsed eddy current methods) have been used experimentally for such applications as coating thickness measurements and the inspection of reactor fuel tubing. The lack of suitable models to facilitate understanding of the interaction of the pulsed field with the test specimen has hindered a wider acceptance of the method as a tool in NDT. Two models, based on the finite element technique are described. The first model, used for repetitive pulse train sources makes use of the Fourier series of the source current to solve a steady state problem for each significant harmonic. The harmonic solutions are then summed to produce the total EMF in the pickup coil. The second model is used for single pulse application. The response is calculated using an iterative time stepping solution. In both cases axisymmetric geometries are studied using a magnetic vector potential formulation. Solutions are compared with experimental results. 3 refs., 3 figs.
Efficient algorithms for use in probabilistic finite element analysis
NASA Technical Reports Server (NTRS)
Dias, J. B.; Nagtegaal, J. C.
1985-01-01
This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.
Finite element implementation of state variable-based viscoplasticity models
NASA Technical Reports Server (NTRS)
Iskovitz, I.; Chang, T. Y. P.; Saleeb, A. F.
1991-01-01
The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.
Chemorheology of reactive systems: Finite element analysis
NASA Technical Reports Server (NTRS)
Douglas, C.; Roylance, D.
1982-01-01
The equations which govern the nonisothermal flow of reactive fluids are outlined, and the means by which finite element analysis is used to solve these equations for the sort of arbitrary boundary conditions encountered in industrial practice are described. The performance of the computer code is illustrated by several trial problems, selected more for their value in providing insight to polymer processing flows than as practical production problems. Although a good deal remains to be learned as to the performance and proper use of this numerical technique, it is undeniably useful in providing better understanding of today's complicated polymer processing problems.
Finite element solutions of free surface flows
NASA Technical Reports Server (NTRS)
Zarda, P. R.; Marcus, M. S.
1977-01-01
A procedure is presented for using NASTRAN to determine the flow field about arbitrarily shaped bodies in the presence of a free surface. The fundamental unknown of the problem is the velocity potential which must satisfy Laplace's equation in the fluid region. Boundary conditions on the free surface may involve second order derivatives in space and time. In cases involving infinite domains either a tractable radiation condition is applied at a truncated boundary or a series expansion is used and matched to the local finite elements. Solutions are presented for harmonic, transient, and steady state problems and compared to either exact solutions or other numerical solutions.
Algebraic surface design and finite element meshes
NASA Technical Reports Server (NTRS)
Bajaj, Chandrajit L.
1992-01-01
Some of the techniques are summarized which are used in constructing C sup 0 and C sup 1 continuous meshes of low degree, implicitly defined, algebraic surface patches in three dimensional space. These meshes of low degree algebraic surface patches are used to construct accurate computer models of physical objects. These meshes are also used in the finite element simulation of physical phenomena (e.g., heat dissipation, stress/strain distributions, fluid flow characteristics) required in the computer prototyping of both the manufacturability and functionality of the geometric design.
Finite element methods in probabilistic mechanics
NASA Technical Reports Server (NTRS)
Liu, Wing Kam; Mani, A.; Belytschko, Ted
1987-01-01
Probabilistic methods, synthesizing the power of finite element methods with second-order perturbation techniques, are formulated for linear and nonlinear problems. Random material, geometric properties and loads can be incorporated in these methods, in terms of their fundamental statistics. By construction, these methods are applicable when the scale of randomness is not too large and when the probabilistic density functions have decaying tails. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. Applications showing the effects of combined random fields and cyclic loading/stress reversal are studied and compared with Monte Carlo simulation results.
Shape optimization including finite element grid adaptation
NASA Technical Reports Server (NTRS)
Kikuchi, N.; Taylor, J. E.
1984-01-01
The prediction of optimal shape design for structures depends on having a sufficient level of precision in the computation of structural response. These requirements become critical in situations where the region to be designed includes stress concentrations or unilateral contact surfaces, for example. In the approach to shape optimization discussed here, a means to obtain grid adaptation is incorporated into the finite element procedures. This facility makes it possible to maintain a level of quality in the computational estimate of response that is surely adequate for the shape design problem.
2-D Finite Element Heat Conduction
Energy Science and Technology Software Center (ESTSC)
1989-10-30
AYER is a finite element program which implicitly solves the general two-dimensional equation of thermal conduction for plane or axisymmetric bodies. AYER takes into account the effects of time (transient problems), in-plane anisotropic thermal conductivity, a three-dimensional velocity distribution, and interface thermal contact resistance. Geometry and material distributions are arbitrary, and input is via subroutines provided by the user. As a result, boundary conditions, material properties, velocity distributions, and internal power generation may be mademore » functions of, e.g., time, temperature, location, and heat flux.« less
Dynamic analysis of mechanisms by finite elements
Botsali, F.M.; Uenuevar, A.
1996-11-01
The need to increase productivity in order to decrease manufacturing costs lead to an increase in the working speeds of machines and mechanical systems used in manufacturing. A method is presented for investigating the dynamics of mechanisms with elastic links. Finite element method is used in the formulation of the dynamic problem. Modal transformation is used in order to reduce the number of equations of motion. Using the presented technique, elastic and rigid body motions of mechanism links are solved simultaneously. The presented method may be applied to spatial and open loop mechanisms including robot manipulators as well.
Adaptive Finite Element Methods in Geodynamics
NASA Astrophysics Data System (ADS)
Davies, R.; Davies, H.; Hassan, O.; Morgan, K.; Nithiarasu, P.
2006-12-01
Adaptive finite element methods are presented for improving the quality of solutions to two-dimensional (2D) and three-dimensional (3D) convection dominated problems in geodynamics. The methods demonstrate the application of existing technology in the engineering community to problems within the `solid' Earth sciences. Two-Dimensional `Adaptive Remeshing': The `remeshing' strategy introduced in 2D adapts the mesh automatically around regions of high solution gradient, yielding enhanced resolution of the associated flow features. The approach requires the coupling of an automatic mesh generator, a finite element flow solver and an error estimator. In this study, the procedure is implemented in conjunction with the well-known geodynamical finite element code `ConMan'. An unstructured quadrilateral mesh generator is utilised, with mesh adaptation accomplished through regeneration. This regeneration employs information provided by an interpolation based local error estimator, obtained from the computed solution on an existing mesh. The technique is validated by solving thermal and thermo-chemical problems with known benchmark solutions. In a purely thermal context, results illustrate that the method is highly successful, improving solution accuracy whilst increasing computational efficiency. For thermo-chemical simulations the same conclusions can be drawn. However, results also demonstrate that the grid based methods employed for simulating the compositional field are not competitive with the other methods (tracer particle and marker chain) currently employed in this field, even at the higher spatial resolutions allowed by the adaptive grid strategies. Three-Dimensional Adaptive Multigrid: We extend the ideas from our 2D work into the 3D realm in the context of a pre-existing 3D-spherical mantle dynamics code, `TERRA'. In its original format, `TERRA' is computationally highly efficient since it employs a multigrid solver that depends upon a grid utilizing a clever
System software for the finite element machine
NASA Technical Reports Server (NTRS)
Crockett, T. W.; Knott, J. D.
1985-01-01
The Finite Element Machine is an experimental parallel computer developed at Langley Research Center to investigate the application of concurrent processing to structural engineering analysis. This report describes system-level software which has been developed to facilitate use of the machine by applications researchers. The overall software design is outlined, and several important parallel processing issues are discussed in detail, including processor management, communication, synchronization, and input/output. Based on experience using the system, the hardware architecture and software design are critiqued, and areas for further work are suggested.
NASA Technical Reports Server (NTRS)
Gelinas, R. J.; Doss, S. K.; Vajk, J. P.; Djomehri, J.; Miller, K.
1983-01-01
The mathematical background regarding the moving finite element (MFE) method of Miller and Miller (1981) is discussed, taking into account a general system of partial differential equations (PDE) and the amenability of the MFE method in two dimensions to code modularization and to semiautomatic user-construction of numerous PDE systems for both Dirichlet and zero-Neumann boundary conditions. A description of test problem results is presented, giving attention to aspects of single square wave propagation, and a solution of the heat equation.
Schwarze, M.; Reese, S.
2007-05-17
Finite element simulations of sheet metal forming processes are highly non-linear problems. The non-linearity arises not only from the kinematical relations and the material formulation, furthermore the contact between workpiece and the forming tools leads to an increased number of iterations within the Newton-Raphson scheme. This fact puts high demands on the robustness of finite element formulations. For this reason we study the enhanced assumed strain (EAS) concept as proposed in [1]. The goal is to improve the robustness of the solid-shell formulation in deep drawing simulations.
A nonlinear dynamic finite element approach for simulating muscular hydrostats.
Vavourakis, V; Kazakidi, A; Tsakiris, D P; Ekaterinaris, J A
2014-01-01
An implicit nonlinear finite element model for simulating biological muscle mechanics is developed. The numerical method is suitable for dynamic simulations of three-dimensional, nonlinear, nearly incompressible, hyperelastic materials that undergo large deformations. These features characterise biological muscles, which consist of fibres and connective tissues. It can be assumed that the stress distribution inside the muscles is the superposition of stresses along the fibres and the connective tissues. The mechanical behaviour of the surrounding tissues is determined by adopting a Mooney-Rivlin constitutive model, while the mechanical description of fibres is considered to be the sum of active and passive stresses. Due to the nonlinear nature of the problem, evaluation of the Jacobian matrix is carried out in order to subsequently utilise the standard Newton-Raphson iterative procedure and to carry out time integration with an implicit scheme. The proposed methodology is implemented into our in-house, open source, finite element software, which is validated by comparing numerical results with experimental measurements and other numerical results. Finally, the numerical procedure is utilised to simulate primitive octopus arm manoeuvres, such as bending and reaching. PMID:23025686
Aeroelastic Stability of Rotor Blades Using Finite Element Analysis
NASA Technical Reports Server (NTRS)
Chopra, I.; Sivaneri, N.
1982-01-01
The flutter stability of flap bending, lead-lag bending, and torsion of helicopter rotor blades in hover is investigated using a finite element formulation based on Hamilton's principle. The blade is divided into a number of finite elements. Quasi-steady strip theory is used to evaluate the aerodynamic loads. The nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The equations of motion are linearized assuming blade motion to be a small perturbation about the steady deflected shape. The normal mode method based on the coupled rotating natural modes is used to reduce the number of equations in the flutter analysis. First the formulation is applied to single-load-path blades (articulated and hingeless blades). Numerical results show very good agreement with existing results obtained using the modal approach. The second part of the application concerns multiple-load-path blades, i.e. bearingless blades. Numerical results are presented for several analytical models of the bearingless blade. Results are also obtained using an equivalent beam approach wherein a bearingless blade is modelled as a single beam with equivalent properties. Results show the equivalent beam model.
A finite element conjugate gradient FFT method for scattering
NASA Technical Reports Server (NTRS)
Collins, Jeffery D.; Zapp, John; Hsa, Chang-Yu; Volakis, John L.
1990-01-01
An extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation (FFT) is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. By virtue of the finite element method, the algorithm is applicable to structures of arbitrary material composition. Several improvements to the two dimensional algorithm are also described. These include: (1) modifications for terminating the mesh at circular boundaries without distorting the convolutionality of the boundary integrals; (2) the development of nonproprietary mesh generation routines for two dimensional applications; (3) the development of preprocessors for interfacing SDRC IDEAS with the main algorithm; and (4) the development of post-processing algorithms based on the public domain package GRAFIC to generate two and three dimensional gray level and color field maps.
Finite element modeling of piezoelectric elements with complex electrode configuration
NASA Astrophysics Data System (ADS)
Paradies, R.; Schläpfer, B.
2009-02-01
It is well known that the material properties of piezoelectric materials strongly depend on the state of polarization of the individual element. While an unpolarized material exhibits mechanically isotropic material properties in the absence of global piezoelectric capabilities, the piezoelectric material properties become transversally isotropic with respect to the polarization direction after polarization. Therefore, for evaluating piezoelectric elements the material properties, including the coupling between the mechanical and the electromechanical behavior, should be addressed correctly. This is of special importance for the micromechanical description of piezoelectric elements with interdigitated electrodes (IDEs). The best known representatives of this group are active fiber composites (AFCs), macro fiber composites (MFCs) and the radial field diaphragm (RFD), respectively. While the material properties are available for a piezoelectric wafer with a homogeneous polarization perpendicular to its plane as postulated in the so-called uniform field model (UFM), the same information is missing for piezoelectric elements with more complex electrode configurations like the above-mentioned ones with IDEs. This is due to the inhomogeneous field distribution which does not automatically allow for the correct assignment of the material, i.e. orientation and property. A variation of the material orientation as well as the material properties can be accomplished by including the polarization process of the piezoelectric transducer in the finite element (FE) simulation prior to the actual load case to be investigated. A corresponding procedure is presented which automatically assigns the piezoelectric material properties, e.g. elasticity matrix, permittivity, and charge vector, for finite element models (FEMs) describing piezoelectric transducers according to the electric field distribution (field orientation and strength) in the structure. A corresponding code has been
Quantum algorithms and the finite element method
NASA Astrophysics Data System (ADS)
Montanaro, Ashley; Pallister, Sam
2016-03-01
The finite element method is used to approximately solve boundary value problems for differential equations. The method discretizes the parameter space and finds an approximate solution by solving a large system of linear equations. Here we investigate the extent to which the finite element method can be accelerated using an efficient quantum algorithm for solving linear equations. We consider the representative general question of approximately computing a linear functional of the solution to a boundary value problem and compare the quantum algorithm's theoretical performance with that of a standard classical algorithm—the conjugate gradient method. Prior work claimed that the quantum algorithm could be exponentially faster but did not determine the overall classical and quantum run times required to achieve a predetermined solution accuracy. Taking this into account, we find that the quantum algorithm can achieve a polynomial speedup, the extent of which grows with the dimension of the partial differential equation. In addition, we give evidence that no improvement of the quantum algorithm can lead to a superpolynomial speedup when the dimension is fixed and the solution satisfies certain smoothness properties.
Impeller deflection and modal finite element analysis.
Spencer, Nathan A.
2013-10-01
Deflections of an impeller due to centripetal forces are calculated using finite element analysis. The lateral, or out of plane, deflections are an important design consideration for this particular impeller because it incorporates an air bearing with critical gap tolerances. The target gap distance is approximately 10 microns at a rotational velocity of 2500 rpm. The centripetal forces acting on the impeller cause it deflect in a concave fashion, decreasing the initial gap distance as a function of radial position. This deflection is characterized for a previous and updated impeller design for comparative purposes. The impact of design options such as material selection, geometry dimensions, and operating rotational velocity are also explored, followed by a sensitivity study with these parameters bounded by specific design values. A modal analysis is also performed to calculate the impeller's natural frequencies which are desired to be avoided during operation. The finite element modeling techniques continue to be exercised by the impeller design team to address specific questions and evaluate conceptual designs, some of which are included in the Appendix.
Finite element analysis of multilayer coextrusion.
Hopkins, Matthew Morgan; Schunk, Peter Randall; Baer, Thomas A.; Mrozek, Randy A.; Lenhart, Joseph Ludlow; Rao, Rekha Ranjana; Collins, Robert; Mondy, Lisa Ann
2011-09-01
Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.
A finite element model for ultrasonic cutting.
Lucas, Margaret; MacBeath, Alan; McCulloch, Euan; Cardoni, Andrea
2006-12-22
Using a single-blade ultrasonic cutting device, a study of ultrasonic cutting of three very different materials is conducted using specimens of cheese, polyurethane foam and epoxy resin. Initial finite element models are created, based on the assumption that the ultrasonic blade causes a crack to propagate in a controlled mode 1 opening, and these are validated against experimental data from three point bend fracture tests and ultrasonic cutting experiments on the materials. Subsequently, the finite element model is developed to represent ultrasonic cutting of a multi-layered material. Materials are chosen whose properties allow a model to be developed that could represent a multi-layer food product or biological structure, to enable ultrasonic cutting systems to be designed for applications both in the field of food processing and surgical procedures. The model incorporates an estimation of the friction condition between the cutting blade and the material to be cut and allows adjustment of the frequency, cutting amplitude and cutting speed. PMID:16814351
Overcoming element erosion limitations within Lagrangian finite element codes
NASA Astrophysics Data System (ADS)
Vignjevic, Rade; Hughes, Kevin; Walker, Andrew; Taylor, Emma A.
2001-10-01
Lagrangian finite element methods have been used extensively in the past to study the non-linear transient behaviour of materials, ranging from crash test of cars to simulating bird strikes on planes.... However, as this type of space discretization does not allow for motion of the material through the mesh when modelling extremely large deformations, the mesh becomes highly distorted. This paper describes some limitations and applicability of this type of analysis for high velocity impacts. A method for dealing with this problem is by the erosion of elements is proposed where the main issue is the deformation of element failure strains. Results were compared with empirical perforation results and were found to be in good agreement. The results were then used to simulate high velocity impacts upon a multi-layered aluminium target, in order to predict a ballistic limit curve. LS-DYNA3D was used as the FE solver for all simulations. Meshes were generated with Truegrid.
A multigrid solution method for mixed hybrid finite elements
Schmid, W.
1996-12-31
We consider the multigrid solution of linear equations arising within the discretization of elliptic second order boundary value problems of the form by mixed hybrid finite elements. Using the equivalence of mixed hybrid finite elements and non-conforming nodal finite elements, we construct a multigrid scheme for the corresponding non-conforming finite elements, and, by this equivalence, for the mixed hybrid finite elements, following guidelines from Arbogast/Chen. For a rectangular triangulation of the computational domain, this non-conforming schemes are the so-called nodal finite elements. We explicitly construct prolongation and restriction operators for this type of non-conforming finite elements. We discuss the use of plain multigrid and the multilevel-preconditioned cg-method and compare their efficiency in numerical tests.
Mixed Finite Element Methods for Melt Migration
NASA Astrophysics Data System (ADS)
Taicher, A. L.
2013-12-01
Multi-phase flow arises during partial melting in the earth mantle, where the porosity is small and material has the characteristics of a compacting porous medium. The equations governing multi-phase flow have been specialized to partially molten materials by McKenzie and Fowler. Their model, also called a Darcy-Stokes system, is highly coupled and non-linear. Melt flow is governed by Darcy's Law while the high temperature, ductile creep of the solid matrix is modeled using viscous non-Newtonian Stokes rheology. In addition, the melt and solid pressures are related through a compaction relation. This nearly elliptic mechanical problem is then coupled with both solute transport and thermal evolution according to the enthalpy method developed by Katz. A suitable numerical method must solve the Darcy-Stokes problem in a manner compatible with the transport problem. Moreover, unlike most porous media problems, partially molten materials transition dynamically from non-porous solid to porous medium so must carefully account for the limit of zero porosity. The Darcy-Stokes system for modeling partial melting in the mantle is a novel problem. As far as we know, there currently does not exist a finite element solution in the literature solving these coupled equations. In particular, the mixed finite element method presents a good candidate because it works in both limiting cases: Darcy and incompressible Stokes flow. We present a mixed formulation for the Darcy-Stokes system. Next, we present novel elements of lowest order and compatible with both Darcy and Stokes flow Finally, we present our 2D mixed FEM code result for solving Stokes and Darcy flow as well as the coupled Darcy-Stokes system the mid-ocean ridge or corner flow problem.
Finite element analysis enhancement of cryogenic testing
NASA Astrophysics Data System (ADS)
Thiem, Clare D.; Norton, Douglas A.
1991-12-01
Finite element analysis (FEA) of large space optics enhances cryogenic testing by providing an analytical method by which to ensure that a test article survives proposed testing. The analyses presented in this paper were concerned with determining the reliability of a half meter mirror in an environment where the exact environmental profile was unknown. FEA allows the interaction between the test object and the environment to be simulated to detect potential problems prior to actual testing. These analyses examined worse case scenerios related to cooling the mirror, its structural integrity for the proposed test environment, and deformation of the reflective surface. The FEA was conducted in-house on the System's Reliability Division's VAX 11-750 and Decstation 3100 using Engineering Mechanics Research Corporation's numerically integrated elements for systems analysis finite element software. The results of the analyses showed that it would take at least 48 hours to cool the mirror to its desired testing temperature. It was also determined that the proposed mirror mount would not cause critical concentrated thermal stresses that would fracture the mirror. FEA and actual measurements of the front reflective face were compared and good agreement between computer simulation and physical tests were seen. Space deployment of large optics requires lightweight mirrors which can perform under the harsh conditions of space. The physical characteristics of these mirrors must be well understood in order that their deployment and operation are successful. Evaluating design approaches by analytical simulation, like FEA, verifies the reliability and structural integrity of a space optic during design prior to prototyping and testing. Eliminating an optic's poor design early in its life saves money, materials, and human resources while ensuring performance.
Elbow stress indices using finite element analysis
NASA Astrophysics Data System (ADS)
Yu, Lixin
Section III of the ASME Boiler and Pressure Vessel Code (the Code) specifies rules for the design of nuclear power plant components. NB-3600 of the Code presents a simplified design method using stress indices---Scalar Coefficients used the modify straight pipe stress equations so that they can be applied to elbows, tees and other piping components. The stress indices of piping components are allowed to be determined both analytically and experimentally. This study concentrates on the determination of B2 stress indices for elbow components using finite element analysis (FEA). First, the previous theoretical, numerical and experimental investigations on elbow behavior were comprehensively reviewed, as was the philosophy behind the use of stress indices. The areas of further research was defined. Then, a comprehensive investigation was carried out to determine how the finite element method should be used to correctly simulate an elbow's structural behavior. This investigation included choice of element type, convergence of mesh density, use of boundary restraint and a reconciliation study between FEA and laboratory experiments or other theoretical formulations in both elastic and elasto-plastic domain. Results from different computer programs were also compared. Reasonably good reconciliation was obtained. Appendix II of the Code describes the experimental method to determine B2 stress indices based on load-deflection curves. This procedure was used to compute the B2 stress indices for various loading modes on one particular elbow configuration. The B2 stress indices thus determined were found to be about half of the value calculated from the Code equation. Then the effect on B2 stress indices of those factors such as internal pressure and flange attachments were studied. Finally, the investigation was extended to other configurations of elbow components. A parametric study was conducted on different elbow sizes and schedules. Regression analysis was then used to
NASA Astrophysics Data System (ADS)
Beilina, Larisa
2016-08-01
We present domain decomposition finite element/finite difference method for the solution of hyperbolic equation. The domain decomposition is performed such that finite elements and finite differences are used in different subdomains of the computational domain: finite difference method is used on the structured part of the computational domain and finite elements on the unstructured part of the domain. Explicit discretizations for both methods are constructed such that the finite element and the finite difference schemes coincide on the common structured overlapping layer between computational subdomains. Then the resulting approach can be considered as a pure finite element scheme which avoids instabilities at the interfaces. We derive an energy estimate for the underlying hyperbolic equation with absorbing boundary conditions and illustrate efficiency of the domain decomposition method on the reconstruction of the conductivity function in three dimensions.
Parallel Finite Element Domain Decomposition for Structural/Acoustic Analysis
NASA Technical Reports Server (NTRS)
Nguyen, Duc T.; Tungkahotara, Siroj; Watson, Willie R.; Rajan, Subramaniam D.
2005-01-01
A domain decomposition (DD) formulation for solving sparse linear systems of equations resulting from finite element analysis is presented. The formulation incorporates mixed direct and iterative equation solving strategics and other novel algorithmic ideas that are optimized to take advantage of sparsity and exploit modern computer architecture, such as memory and parallel computing. The most time consuming part of the formulation is identified and the critical roles of direct sparse and iterative solvers within the framework of the formulation are discussed. Experiments on several computer platforms using several complex test matrices are conducted using software based on the formulation. Small-scale structural examples are used to validate thc steps in the formulation and large-scale (l,000,000+ unknowns) duct acoustic examples are used to evaluate the ORIGIN 2000 processors, and a duster of 6 PCs (running under the Windows environment). Statistics show that the formulation is efficient in both sequential and parallel computing environmental and that the formulation is significantly faster and consumes less memory than that based on one of the best available commercialized parallel sparse solvers.
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 analyses is presented. New thermal finite elements which yield exact nodal and element temperature 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 Viscoelastic Hybrid Shell Finite Element
NASA Technical Reports Server (NTRS)
Johnson, Arthur
1999-01-01
An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.
Asymmetric quadrilateral shell elements for finite strains
NASA Astrophysics Data System (ADS)
Areias, P.; Dias-da-Costa, D.; Pires, E. B.; Van Goethem, N.
2013-07-01
Very good results in infinitesimal and finite strain analysis of shells are achieved by combining either the enhanced-metric technique or the selective-reduced integration for the in-plane shear energy and an assumed natural strain technique (ANS) in a non-symmetric Petrov-Galerkin arrangement which complies with the patch-test. A recovery of the original Wilson incompatible mode element is shown for the trial functions in the in-plane components. As a beneficial side-effect, Newton-Raphson convergence behavior for non-linear problems is improved with respect to symmetric formulations. Transverse-shear and in-plane patch tests are satisfied while distorted-mesh accuracy is higher than with symmetric formulations. Classical test functions with assumed-metric components are required for compatibility reasons. Verification tests are performed with advantageous comparisons being observed in all of them. Applications to large displacement elasticity and finite strain plasticity are shown with both low sensitivity to mesh distortion and (relatively) high accuracy. A equilibrium-consistent (and consistently linearized) updated-Lagrangian algorithm is proposed and tested. Concerning the time-step dependency, it was found that the consistent updated-Lagrangian algorithm is nearly time-step independent and can replace the multiplicative plasticity approach if only moderate elastic strains are present, as is the case of most metals.
Finite element or Galerkin type semidiscrete schemes
NASA Technical Reports Server (NTRS)
Durgun, K.
1983-01-01
A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.
Finite-element modeling of nanoindentation
Knapp, J.A.; Follstaedt, D.M.; Myers, S.M.; Barbour, J.C.; Friedmann, T.A.
1999-02-01
Procedures have been developed based on finite-element modeling of nanoindentation data to obtain the mechanical properties of thin films and ion-beam-modified layers independently of the properties of the underlying substrates. These procedures accurately deduce the yield strength, Young{close_quote}s elastic modulus, and layer hardness from indentations as deep as 50{percent} of the layer thickness or more. We have used these procedures to evaluate materials ranging from ion implanted metals to deposited, diamond-like carbon layers. The technique increases the applicability of indentation testing to very thin layers, composite layers, and modulated compositions. This article presents an overview of the procedures involved and illustrates them with selected examples. {copyright} {ital 1999 American Institute of Physics.}
Fuzzy finite element analysis of smart structures
NASA Astrophysics Data System (ADS)
Akpan, Unyime O.; Koko, Tamunoiyala S.; Orisamolu, Irewole R.; Gallant, B. Keith
2000-06-01
A fuzzy finite element based approach is developed for modelling smart structures with vague or imprecise uncertainties. Fuzzy sets are used to represent the uncertainties present in the piezoelectric, mechanical, thermal, and physical properties of the smart structure. In order to facilitate efficient computation, a sensitivity analysis procedure is used to streamline the number of input fuzzy variables, and the vertex fuzzy analysis technique is then used to compute the possibility distributions of the responses of the smart structural system. The methodology has been developed within the framework of the SMARTCOM computational tool for the design/analysis of smart composite structures. The methodology developed is found to be accurate and computationally efficient for solution of practical problems.
Continuation finite element analysis of viscoelastic fluids
NASA Astrophysics Data System (ADS)
Chow, Tai-Whang
A finite element procedure using a mixed formulation and a predictor-corrector type continuation algorithm for the analysis of two dimensional steady state flows of viscoelastic fluids is described. As a simple but nontrivial test example, radial flow immenating from a line by the numerical discretization and believed to be the cause for previous numerical failures, are shown and branch solution paths are followed by step length adjustment and by convergent tolerance relaxation. A technique for jumping over bifurcation points is presented and used to increase the Weissenberg number with no apparent limit for the radial flow problem. A second example related to extrusion of viscoelastic material is also analyzed. Steady state velocity fields, deviatoric stress distributions and pressure distributions for several different Weissenberg numbers are presented with bifurcation points and turning points noted.
Quality management of finite element analysis
NASA Astrophysics Data System (ADS)
Barlow, John
1991-09-01
A quality management system covering the use of finite element analysis is described. The main topics are as follows: acquisition, development and verification of software (including the software suppliers software quality control system), support, documentation, error control, internal software, software acceptance and release; development and qualification of analysis methods, including software evaluation, analysis procedure qualification and documentation, procedure quality checks, control of analysis procedure errors; product design and integrity analysis, including project quality assurance and analysis planning, task specification and allocation, analysis, execution, results checking and analysis records. Other issues include the commercial and business advantages of quality systems, project and technical management and the training and experience of personnel. The items are correlated with the requirements of International Standard Organization 9001.
Optimizing electroslag cladding with finite element modeling
Li, M.V.; Atteridge, D.G.; Meekisho, L.
1996-12-31
Electroslag cladding of nickel alloys onto carbon steel propeller shafts was optimized in terms of interpass temperatures. A two dimensional finite element model was used in this study to analyze the heat transfer induced by multipass electroslag cladding. Changes of interpass temperatures during a cladding experiment with uniform initial temperature distribution on a section of shaft were first simulated. It was concluded that uniform initial temperature distribution would lead to interpass temperatures out of the optimal range if continuous cladding is expected. The difference in the cooling conditions among experimental and full size shafts and its impact on interpass temperatures during the cladding were discussed. Electroslag cladding onto a much longer shaft, virtually an semi infinite long shaft, was analyzed with specific reference to the practical applications of electroslag cladding. Optimal initial preheating temperature distribution was obtained for continuous cladding on full size shafts which would keep the interpass temperatures within the required range.
Finite element simulation of pipe dynamic response
Slagis, G.C.; Litton, R.W.
1996-12-01
Nonlinear finite element dynamic analyses of the response of a pipe span to controlled-displacement, sinusoidal vibration have been performed. The objective of this preliminary study is to compare strain and acceleration response data to those generated by Beaney in the Berkeley Nuclear Laboratories experiments. Results for an unpressurized, 5 Hz, carbon steel pipe are in good agreement with the experiments. Hence, it appears that analytical simulation will be useful to assess seismic margins. Recommendations for additional studies are provided. The analyses confirm the test results--dynamic response is greatly attenuated by material plasticity. Analytical strains and accelerations are about 30% higher than test data. There are several possible explanations for the differences. To assess the effect of frequency on response, the length of the pipe span was increased. Analysis of the longer, 2 Hz, pipe span shows significantly greater cyclic strains than the 5 Hz span at the same input excitation levels.
3-D Finite Element Heat Transfer
Energy Science and Technology Software Center (ESTSC)
1992-02-01
TOPAZ3D is a three-dimensional implicit finite element computer code for heat transfer analysis. TOPAZ3D can be used to solve for the steady-state or transient temperature field on three-dimensional geometries. Material properties may be temperature-dependent and either isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functionalmore » representation of boundary conditions and internal heat generation. TOPAZ3D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
Finite element analysis: A boon to dentistry
Trivedi, Shilpa
2014-01-01
The finite element analysis (FEA) is an upcoming and significant research tool for biomechanical analyses in biological research. It is an ultimate method for modeling complex structures and analyzing their mechanical properties. In Implantology, FEA has been used to study the stress patterns in various implant components and also in the peri-implant bone. It is also useful for studying the biomechanical properties of implants as well as for predicting the success of implants in clinical condition. FEA of simulated traumatic loads can be used to understand the biomechanics of fracture. FEA has various advantages compared with studies on real models. The experiments are repeatable, there are no ethical considerations and the study designs may be modified and changed as per the requirement. There are certain limitations of FEA too. It is a computerized in vitro study in which clinical condition may not be completely replicated. So, further FEA research should be supplemented with clinical evaluation. PMID:25737944
Boundary element and finite element coupling for aeroacoustics simulations
NASA Astrophysics Data System (ADS)
Balin, Nolwenn; Casenave, Fabien; Dubois, François; Duceau, Eric; Duprey, Stefan; Terrasse, Isabelle
2015-08-01
We consider the scattering of acoustic perturbations in the presence of a flow. We suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential. In the first zone, we apply a Prandtl-Glauert transformation to recover the Helmholtz equation. The well-known setting of boundary element method for the Helmholtz equation is available. In the second zone, the flow quantities are space dependent, we have to consider a local resolution, namely the finite element method. Herein, we carry out the coupling of these two methods and present various applications and validation test cases. The source term is given through the decomposition of an incident acoustic field on a section of the computational domain's boundary. Validations against analytic, another numerical method and measurements on different test cases are presented.
Structural optimization of thin shells using finite element method
NASA Technical Reports Server (NTRS)
Gotsis, Pascal K.
1992-01-01
The objective of the present work was the structural optimization of thin shell structures that are subjected to stress and displacement constraints. In order to accomplish this, the structural optimization computer program DESAP1 was modified and improved. In the static analysis part of the DESAP1 computer program the torsional spring elements, which are used to analyze thin, shallow shell structures, were eliminated by modifying the membrane stiffness matrix of the triangular elements in the local coordinate system and adding a fictitious rotational stiffness matrix. This simplified the DESAP1 program input, improved the accuracy of the analysis, and saved computation time. In the optimization part of the DESAP1 program the stress ratio formula, which redesigns the thickness of each finite element of the structure, was solved by an analytical method. This scheme replaced the iterative solution that was previously used in the DESAP1 program, thus increasing the accuracy and speed of the design. The modified program was used to design a thin, cylindrical shell structure with optimum weight, and the results are reported in this paper.
Finite element analyses of CCAT preliminary design
NASA Astrophysics Data System (ADS)
Sarawit, Andrew T.; Kan, Frank W.
2014-07-01
This paper describes the development of the CCAT telescope finite element model (FEM) and the analyses performed to support the preliminary design work. CCAT will be a 25 m diameter telescope operating in the 0.2 to 2 mm wavelength range. It will be located at an elevation of 5600 m on Cerro Chajnantor in Northern Chile, near ALMA. The telescope will be equipped with wide-field cameras and spectrometers mounted at the two Nasmyth foci. The telescope will be inside an enclosure to protect it from wind buffeting, direct solar heating, and bad weather. The main structures of the telescope include a steel Mount and a carbon-fiber-reinforced-plastic (CFRP) primary truss. The finite element model developed in this study was used to perform modal, frequency response, seismic response spectrum, stress, and deflection analyses of telescope. Modal analyses of telescope were performed to compute the structure natural frequencies and mode shapes and to obtain reduced order modal output at selected locations in the telescope structure to support the design of the Mount control system. Modal frequency response analyses were also performed to compute transfer functions at these selected locations. Seismic response spectrum analyses of the telescope subject to the Maximum Likely Earthquake were performed to compute peak accelerations and seismic demand stresses. Stress analyses were performed for gravity load to obtain gravity demand stresses. Deflection analyses for gravity load, thermal load, and differential elevation drive torque were performed so that the CCAT Observatory can verify that the structures meet the stringent telescope surface and pointing error requirements.
Finite-element-based photoacoustic tomography: phantom and chicken bone experiments
NASA Astrophysics Data System (ADS)
Yuan, Zhen; Zhao, Hongzhi; Wu, Changfeng; Zhang, Qizhi; Jiang, Huabei
2006-05-01
We describe a photoacoustic image reconstruction algorithm that is based on the finite-element solution to the photoacoustic wave equation in the frequency domain. Our reconstruction approach is an iterative Newton method coupled with combined Marquardt and Tikhonov regularizations that can extract the spatial distribution of optical-absorption property in heterogeneous media. We demonstrate this algorithm by using phantom and chicken bone measurements from a circular scanning photoacoustic tomography system. The results obtained show that millimeter-sized phantom objects and chicken bones and/or joints can be clearly detected using our finite-element-based photoacoustic tomography method.
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.
NASA Astrophysics Data System (ADS)
Zhu, Yu; Cangellaris, Andreas C.
2002-05-01
A new finite element methodology is presented for fast and robust numerical simulation of three-dimensional electromagnetic wave phenomena. The new methodology combines nested multigrid techniques with the ungauged vector and scalar potential formulation of the finite element method. The finite element modeling is performed on nested meshes over the computational domain of interest. The iterative solution of the finite element matrix on the finest mesh is performed using the conjugate gradient method, while the nested multigrid vector and scalar potential algorithm acts as the preconditioner for the iterative solver. Numerical experiments from the application of the new methodology to three-dimensional electromagnetic scattering are used to demonstrate its superior numerical convergence and efficient memory usage.
NASA Technical Reports Server (NTRS)
1976-01-01
The development of two new shell finite elements for applications to large deflection problems is considered. The elements in question are doubly curved and of triangular and quadrilateral planform. They are restricted to small strains of elastic materials, and can accommodate large rotations. The elements described, which are based on relatively simple linear elements, make use of a new displacement function approach specifically designed for strongly nonlinear problems. The displacement function development for nonlinear applications is based on certain beam element formulations, and the strain-displacement equations are of a shallow shell type. Additional terms were included in these equations in an attempt to avoid the large errors characteristic of shallow shell elements in certain types of problems. An incremental nonlinear solution procedure specifically adopted to the element formulation was developed. The solution procedure is of combined incremental and total Lagrangian type, and uses a new updating scheme. A computer program was written to evaluate the developed formulations. This program can accommodate small element groups in arbitrary arrangements. Two simple programs were successfully solved. The results indicate that this new type of element has definite promise and should be a fruitful area for further research.
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.
Impact of new computing systems on finite element computations
NASA Technical Reports Server (NTRS)
Noor, A. K.; Storassili, O. O.; Fulton, R. E.
1983-01-01
Recent advances in computer technology that are likely to impact finite element computations are reviewed. The characteristics of supersystems, highly parallel systems, and small systems (mini and microcomputers) are summarized. The interrelations of numerical algorithms and software with parallel architectures are discussed. A scenario is presented for future hardware/software environment and finite element systems. A number of research areas which have high potential for improving the effectiveness of finite element analysis in the new environment are identified.
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.
NASA Astrophysics Data System (ADS)
Kouhi, Mohammad; Oñate, Eugenio
2015-07-01
A new implicit stabilized formulation for the numerical solution of the compressible Navier-Stokes equations is presented. The method is based on the finite calculus (FIC) scheme using the Galerkin finite element method (FEM) on triangular grids. Via the FIC formulation, two stabilization terms, called streamline term and transverse term, are added to the original conservation equations in the space-time domain. The non-linear system of equations resulting from the spatial discretization is solved implicitly using a damped Newton method benefiting from the exact Jacobian matrix. The matrix system is solved at each iteration with a preconditioned GMRES method. The efficiency of the proposed stabilization technique is checked out in the solution of 2D inviscid and laminar viscous flow problems where appropriate solutions are obtained especially near the boundary layer and shock waves. The work presented here can be considered as a follow up of a previous work of the authors Kouhi, Oñate (Int J Numer Methods Fluids 74:872-897, 2014). In that paper, the stabilized Galerkin FEM based on the FIC formulation was derived for the Euler equations together with an explicit scheme. In the present paper, the extension of this work to the Navier-Stokes equations using an implicit scheme is presented.
Finite element analysis in a minicomputer/mainframe environment
NASA Technical Reports Server (NTRS)
Storaasli, O. O.; Murphy, R. C.
1978-01-01
Design considerations were evaluated for general purpose finite element systems to maximize performance when installed on distributed computer hardware/software systems. It is shown how the features of current minicomputers complement those of a modular implementation of the finite element method for increasing the control, speed, and visibility (interactive graphics) in solving structural problems at reduced cost. The approach used is to implement a finite element system in a distributed computer environment to solve structural problems and to explore alternatives in distributing finite element computations.
A multi-microprocessor system for finite element structural analysis
NASA Technical Reports Server (NTRS)
Jordan, H. F.; Sawyer, P. L.
1978-01-01
During the last few years, advances in microprocessor technology have spurred a renewed interest in special-purpose computers. The microprocessor has become small, inexpensive, and powerful enough to be considered as a building block for special-purpose hardware. A description is presented of the architecture of a prototype 'finite element machine' currently being built. Attention is given to details regarding the finite element analysis problem, the arrangement of the processors as finite element nodes in the structural model, the influence of the architecture on the solution algorithm, interprocessor communication primitives, and the performance of the finite element machine.
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.
NASA Technical Reports Server (NTRS)
Nakazawa, Shohei
1989-01-01
The internal structure is discussed of the MHOST finite element program designed for 3-D inelastic analysis of gas turbine hot section components. The computer code is the first implementation of the mixed iterative solution strategy for improved efficiency and accuracy over the conventional finite element method. The control structure of the program is covered along with the data storage scheme and the memory allocation procedure and the file handling facilities including the read and/or write sequences.
A fast algorithm for nonlinear finite element analysis using equivalent magnetization current
NASA Astrophysics Data System (ADS)
Lee, Joon-Ho; Park, Il-Han; Kim, Dong-Hun; Lee, Ki-Sik
2002-05-01
A fast algorithm for iterative nonlinear finite element analysis is presented in this paper. The algorithm replaces updated permeability by an equivalent magnetization current and moves it to the source current term. Once the initial system matrix is decomposed in the LU form, the iterative procedure involves the trivial step of back-substitution from the LU form. Consequently, the computation time for the nonlinear analysis is greatly reduced. A numerical model of a cylindrical conductor enclosed with saturable iron is tested to validate the proposed algorithm. Numerical results are compared with those obtained using conventional Newton-Raphson algorithm in respect to accuracy and computational time.
Integrated finite element model of composite materials
NASA Astrophysics Data System (ADS)
Teply, Jan L.; Herbein, William C.
1989-05-01
Two problems traditionally addressed in the area of micromechanics of composite materials can be briefly summarized as follows: (1) for a macroscopically uniform volume of composite material, which is subjected to macroscopically uniform boundary tractions, displacements or heat influx, find overall thermomechanical properties in terms of the thermomechanical properties of the individual constituents; and (2) for the same material volume and boundary conditions as above, find the local stress, strain, and temperature fields in the constituents and on the interfaces. Two different types of micromechanical models are usually applied to the solutions of these two types of problems. For linear elastic materials, the micromechanical models to solve problem (1) offer simple solutions of overall thermomechanical properties either in terms of bound which are derived from periodic or random microstructures, or in terms of single estimates, which are derived from a solution of an isolated inclusion. The finite element variational approaches are applied to integrate the solutions of problems (1) and (2) into one model. The application of displacement and equilibrium variational approaches to the calculation of overall elastic-plastic properties, are extended to the solution of the second problem. The integrated model is then applied to calculate the overall properties and local stress and strain fields of boron-aluminum composites subjected to transverse tension, in-plane shear and bending.
Laterally displaced pipelines: Finite element analysis
Altaee, A.; Boivin, R.
1995-12-31
The rate effect of lateral soil movement against buried pipes in clay soils is investigated in finite element analyzes using two different computer programs, AGAC and CRISP. Rapid and slow ground movements are considered in ideal undrained and ideal drained analysis, respectively, which represent the two extreme boundaries with respect to rate of loading (rate of ground movement). The analyses address a typical full-scale buried pipe as described by Rizkalla et al. (1992). The pipe considered for the analysis has a diameter of 0.914 m and is placed in a backfilled 2.0 m wide and 1.8 m deep excavation. Results from both AGAC and CRISP analyzes are similar in terms of total lateral force versus lateral pipe movement. For example, both programs indicate the same clear difference in the resulting pipe movement for cases of rapid and slow ground movement, especially at large movement. When the ground movement is rapid, the pipe moves both laterally and upward. One the other hand, when the ground movement is slow, the pipe experiences only lateral movement and no noticeable vertical movement. The total force acting on the pipe (and stresses and strains within the pipe) is larger for the slow rate of loading. The results of analyzes presented herein agree with results of tests on a 5.5 m beam centrifuge performed by the Center for Cold Oceans Resources Engineering.
Finite Element Modeling of Human Placental Tissue
Yu, Mao; Manoogian, Sarah; Duma, Stefan M.; Stitzel, Joel D.
2009-01-01
Motor vehicle crashes account for a large portion of placental abruption and fetal losses. To better understand the material properties of the human placenta, a Finite Element (FE) model of human placenta tissue was created and verified using data from uniaxial tension tests. Sixty-four tensile tests at three different strain rates of 7% strain/s, 70% strain/s, and 700% strain/s from six whole human placentas were used for model development. Nominal stresses were calculated by dividing forces at the grips by the original cross-sectional area. Nominal strains were calculated by dividing cross-head displacement by the original gauge length. A detailed methodology for interpreting experimental data for application to material model development is presented. A model of the tension coupon was created in LS-DYNA and stretched in the same manner as the uniaxial tension tests. The behavior of the material was optimized to the uniaxial tension test using a multi-island genetic algorithm. The results demonstrate good correlation between experiments and the model, with an average difference of 2% between the optimized FE and experimental first principal stress at the termination state. The material parameters found in this study can be utilized in FE models of placental tissues for behavior under dynamic loading. PMID:20184849
Finite element modeling of retinal prosthesis mechanics
NASA Astrophysics Data System (ADS)
Basinger, B. C.; Rowley, A. P.; Chen, K.; Humayun, M. S.; Weiland, J. D.
2009-10-01
Epiretinal prostheses used to treat degenerative retina diseases apply stimulus via an electrode array fixed to the ganglion cell side of the retina. Mechanical pressure applied by these arrays to the retina, both during initial insertion and throughout chronic use, could cause sufficient retinal damage to reduce the device's effectiveness. In order to understand and minimize potential mechanical damage, we have used finite element analysis to model mechanical interactions between an electrode array and the retina in both acute and chronic loading configurations. Modeling indicates that an acute tacking force distributes stress primarily underneath the tack site and heel edge of the array, while more moderate chronic stresses are distributed more evenly underneath the array. Retinal damage in a canine model chronically implanted with a similar array occurred in correlating locations, and model predictions correlate well with benchtop eyewall compression tests. This model provides retinal prosthesis researchers with a tool to optimize the mechanical electrode array design, but the techniques used here represent a unique effort to combine a modifiable device and soft biological tissues in the same model and those techniques could be extended to other devices that come into mechanical contact with soft neural tissues.
TACO: a finite element heat transfer code
Mason, W.E. Jr.
1980-02-01
TACO is a two-dimensional implicit finite element code for heat transfer analysis. It can perform both linear and nonlinear analyses and can be used to solve either transient or steady state problems. Either plane or axisymmetric geometries can be analyzed. TACO has the capability to handle time or temperature dependent material properties and materials may be either isotropic or orthotropic. A variety of time and temperature dependent loadings and boundary conditions are available including temperature, flux, convection, and radiation boundary conditions and internal heat generation. Additionally, TACO has some specialized features such as internal surface conditions (e.g., contact resistance), bulk nodes, enclosure radiation with view factor calculations, and chemical reactive kinetics. A user subprogram feature allows for any type of functional representation of any independent variable. A bandwidth and profile minimization option is also available in the code. Graphical representation of data generated by TACO is provided by a companion post-processor named POSTACO. The theory on which TACO is based is outlined, the capabilities of the code are explained, the input data required to perform an analysis with TACO are described. Some simple examples are provided to illustrate the use of the code.
Updated lagrangian mixed finite element formulation for quasi and fully incompressible fluids
NASA Astrophysics Data System (ADS)
Oñate, Eugenio; Carbonell, Josep M.
2014-12-01
We present a mixed velocity-pressure finite element formulation for solving the updated Lagrangian equations for quasi and fully incompressible fluids. Details of the governing equations for the conservation of momentum and mass are given in both differential and variational form. The finite element interpolation uses an equal order approximation for the velocity and pressure unknowns. The procedure for obtaining stabilized FEM solutions is outlined. The solution in time of the discretized governing conservation equations using an incremental iterative segregated scheme is described. The linearization of these equations and the derivation of the corresponding tangent stiffness matrices is detailed. Other iterative schemes for the direct computation of the nodal velocities and pressures at the updated configuration are presented. The advantages and disadvantages of choosing the current or the updated configuration as the reference configuration in the Lagrangian formulation are discussed.
Cwik, T.; Katz, D.S.
1996-12-31
Finite element modeling has proven useful for accurately simulating scattered or radiated electromagnetic fields from complex three-dimensional objects whose geometry varies on the scale of a fraction of an electrical wavelength. An unstructured finite element model of realistic objects leads to a large, sparse, system of equations that needs to be solved efficiently with regard to machine memory and execution time. Both factorization and iterative solvers can be used to produce solutions to these systems of equations. Factorization leads to high memory requirements that limit the electrical problem size of three-dimensional objects that can be modeled. An iterative solver can be used to efficiently solve the system without excessive memory use and in a minimal amount of time if the convergence rate is controlled.
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)
Large spatial, temporal, and algorithmic adaptivity for implicit nonlinear finite element analysis
Engelmann, B.E.; Whirley, R.G.
1992-07-30
The development of effective solution strategies to solve the global nonlinear equations which arise in implicit finite element analysis has been the subject of much research in recent years. Robust algorithms are needed to handle the complex nonlinearities that arise in many implicit finite element applications such as metalforming process simulation. The authors experience indicates that robustness can best be achieved through adaptive solution strategies. In the course of their research, this adaptivity and flexibility has been refined into a production tool through the development of a solution control language called ISLAND. This paper discusses aspects of adaptive solution strategies including iterative procedures to solve the global equations and remeshing techniques to extend the domain of Lagrangian methods. Examples using the newly developed ISLAND language are presented to illustrate the advantages of embedding temporal, algorithmic, and spatial adaptivity in a modem implicit nonlinear finite element analysis code.
A finite element model updating technique for adjustment of parameters near boundaries
NASA Astrophysics Data System (ADS)
Gwinn, Allen Fort, Jr.
Even though there have been many advances in research related to methods of updating finite element models based on measured normal mode vibration characteristics, there is yet to be a widely accepted method that works reliably with a wide range of problems. This dissertation focuses on the specific class of problems having to do with changes in stiffness near the clamped boundary of plate structures. This class of problems is especially important as it relates to the performance of turbine engine blades, where a change in stiffness at the base of the blade can be indicative of structural damage. The method that is presented herein is a new technique for resolving the differences between the physical structure and the finite element model. It is a semi-iterative technique that incorporates a "physical expansion" of the measured eigenvectors along with appropriate scaling of these expanded eigenvectors into an iterative loop that uses the Engel's model modification method to then calculate adjusted stiffness parameters for the finite element model. Three example problems are presented that use eigenvalues and mass normalized eigenvectors that have been calculated from experimentally obtained accelerometer readings. The test articles that were used were all thin plates with one edge fully clamped. They each had a cantilevered length of 8.5 inches and a width of 4 inches. The three plates differed from one another in thickness from 0.100 inches to 0.188 inches. These dimensions were selected in order to approximate a gas turbine engine blade. The semi-iterative modification technique is shown to do an excellent job of calculating the necessary adjustments to the finite element model so that the analytically determined eigenvalues and eigenvectors for the adjusted model match the corresponding values from the experimental data with good agreement. Furthermore, the semi-iterative method is quite robust. For the examples presented here, the method consistently converged
Application of edge-based finite elements and vector ABCs in 3D scattering
NASA Technical Reports Server (NTRS)
Chatterjee, A.; Jin, J. M.; Volakis, John L.
1992-01-01
A finite element absorbing boundary condition (FE-ABC) solution of the scattering by arbitrary 3-D structures is considered. The computational domain is discretized using edge-based tetrahedral elements. In contrast to the node-based elements, edge elements can treat geometries with sharp edges, are divergence-less, and easily satisfy the field continuity condition across dielectric interfaces. They do, however, lead to a higher unknown count but this is balanced by the greater sparsity of the resulting finite element matrix. Thus, the computation time required to solve such a system iteratively with a given degree of accuracy is less than the traditional node-based approach. The purpose is to examine the derivation and performance of the ABC's when applied to 2-D and 3-D problems and to discuss the specifics of our FE-ABC implementation.
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.
FINITE-ELEMENT ANALYSIS OF MULTIPHASE IMMISCIBLE FLOW THROUGH SOILS
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 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.
Modular Finite Element Methods Library Version: 1.0
Energy Science and Technology Software Center (ESTSC)
2010-06-22
MFEM is a general, modular library for finite element methods. It provides a variety of finite element spaces and bilinear/linear forms in 2D and 3D. MFEM also includes classes for dealing with various types of meshes and their refinement.
Finite Element analyses of soil bioengineered slopes
NASA Astrophysics Data System (ADS)
Tamagnini, Roberto; Switala, Barbara Maria; Sudan Acharya, Madhu; Wu, Wei; Graf, Frank; Auer, Michael; te Kamp, Lothar
2014-05-01
Soil Bioengineering methods are not only effective from an economical point of view, but they are also interesting as fully ecological solutions. The presented project is aimed to define a numerical model which includes the impact of vegetation on slope stability, considering both mechanical and hydrological effects. In this project, a constitutive model has been developed that accounts for the multi-phase nature of the soil, namely the partly saturated condition and it also includes the effects of a biological component. The constitutive equation is implemented in the Finite Element (FE) software Comes-Geo with an implicit integration scheme that accounts for the collapse of the soils structure due to wetting. The mathematical formulation of the constitutive equations is introduced by means of thermodynamics and it simulates the growth of the biological system during the time. The numerical code is then applied in the analysis of an ideal rainfall induced landslide. The slope is analyzed for vegetated and non-vegetated conditions. The final results allow to quantitatively assessing the impact of vegetation on slope stability. This allows drawing conclusions and choosing whenever it is worthful to use soil bioengineering methods in slope stabilization instead of traditional approaches. The application of the FE methods show some advantages with respect to the commonly used limit equilibrium analyses, because it can account for the real coupled strain-diffusion nature of the problem. The mechanical strength of roots is in fact influenced by the stress evolution into the slope. Moreover, FE method does not need a pre-definition of any failure surface. FE method can also be used in monitoring the progressive failure of the soil bio-engineered system as it calculates the amount of displacements and strains of the model slope. The preliminary study results show that the formulated equations can be useful for analysis and evaluation of different soil bio
Nondestructive Evaluation Correlated with Finite Element Analysis
NASA Technical Reports Server (NTRS)
Abdul-Azid, Ali; Baaklini, George Y.
1999-01-01
Advanced materials are being developed for use in high-temperature gas turbine applications. For these new materials to be fully utilized, their deformation properties, their nondestructive evaluation (NDE) quality and material durability, and their creep and fatigue fracture characteristics need to be determined by suitable experiments. The experimental findings must be analyzed, characterized, modeled and translated into constitutive equations for stress analysis and life prediction. Only when these ingredients - together with the appropriate computational tools - are available, can durability analysis be performed in the design stage, long before the component is built. One of the many structural components being evaluated by the NDE group at the NASA Lewis Research Center is the flywheel system. It is being considered as an energy storage device for advanced space vehicles. Such devices offer advantages over electrochemical batteries in situations demanding high power delivery and high energy storage per unit weight. In addition, flywheels have potentially higher efficiency and longer lifetimes with proper motor-generator and rotor design. Flywheels made of fiber-reinforced polymer composite material show great promise for energy applications because of the high energy and power densities that they can achieve along with a burst failure mode that is relatively benign in comparison to those of flywheels made of metallic materials Therefore, to help improve durability and reduce structural uncertainties, we are developing a comprehensive analytical approach to predict the reliability and life of these components under these harsh loading conditions. The combination of NDE and two- and three-dimensional finite element analyses (e.g., stress analyses and fracture mechanics) is expected to set a standardized procedure to accurately assess the applicability of using various composite materials to design a suitable rotor/flywheel assembly.
FEWA: a Finite Element model of Water flow through Aquifers
Yeh, G.T.; Huff, D.D.
1983-11-01
This report documents the implementation and demonstration of a Finite Element model of Water flow through Aquifers (FEWA). The particular features of FEWA are its versatility and flexibility to deal with as many real-world problems as possible. Point as well as distributed sources/sinks are included to represent recharges/pumpings and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Source/sink strength over each element and node, hydraulic head at each Dirichlet boundary node, and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution method for the matrix equation approximating the partial differential equation of groundwater flow. FEWA also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. It is then demonstrated by two examples of how the model can be applied to heterogeneous and anisotropic aquifers with transient boundary conditions, time-dependent sources/sinks, and confining aquitards for a confined aquifer of variable thickness and for a free surface problem in an unconfined aquifer, respectively. 20 references, 25 figures, 8 tables.
Advances in 3D electromagnetic finite element modeling
Nelson, E.M.
1997-08-01
Numerous advances in electromagnetic finite element analysis (FEA) have been made in recent years. The maturity of frequency domain and eigenmode calculations, and the growth of time domain applications is briefly reviewed. A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will also be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis is also discussed.
A hybrid finite-difference and analytic element groundwater model.
Haitjema, H M; Feinstein, D T; Hunt, R J; Gusyev, M A
2010-01-01
Regional finite-difference models tend to have large cell sizes, often on the order of 1-2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW-MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models. PMID:20132324
ELLIPT2D: A Flexible Finite Element Code Written Python
Pletzer, A.; Mollis, J.C.
2001-03-22
The use of the Python scripting language for scientific applications and in particular to solve partial differential equations is explored. It is shown that Python's rich data structure and object-oriented features can be exploited to write programs that are not only significantly more concise than their counter parts written in Fortran, C or C++, but are also numerically efficient. To illustrate this, a two-dimensional finite element code (ELLIPT2D) has been written. ELLIPT2D provides a flexible and easy-to-use framework for solving a large class of second-order elliptic problems. The program allows for structured or unstructured meshes. All functions defining the elliptic operator are user supplied and so are the boundary conditions, which can be of Dirichlet, Neumann or Robbins type. ELLIPT2D makes extensive use of dictionaries (hash tables) as a way to represent sparse matrices.Other key features of the Python language that have been widely used include: operator over loading, error handling, array slicing, and the Tkinter module for building graphical use interfaces. As an example of the utility of ELLIPT2D, a nonlinear solution of the Grad-Shafranov equation is computed using a Newton iterative scheme. A second application focuses on a solution of the toroidal Laplace equation coupled to a magnetohydrodynamic stability code, a problem arising in the context of magnetic fusion research.
Three dimensional inelastic finite element analysis of laminated composites
NASA Technical Reports Server (NTRS)
Griffin, O. H., Jr.; Kamat, M. P.
1980-01-01
Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.
The L sub 1 finite element method for pure convection problems
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan
1991-01-01
The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.
Interpolation functions in the immersed boundary and finite element methods
NASA Astrophysics Data System (ADS)
Wang, Xingshi; Zhang, Lucy T.
2010-03-01
In this paper, we review the existing interpolation functions and introduce a finite element interpolation function to be used in the immersed boundary and finite element methods. This straightforward finite element interpolation function for unstructured grids enables us to obtain a sharper interface that yields more accurate interfacial solutions. The solution accuracy is compared with the existing interpolation functions such as the discretized Dirac delta function and the reproducing kernel interpolation function. The finite element shape function is easy to implement and it naturally satisfies the reproducing condition. They are interpolated through only one element layer instead of smearing to several elements. A pressure jump is clearly captured at the fluid-solid interface. Two example problems are studied and results are compared with other numerical methods. A convergence test is thoroughly conducted for the independent fluid and solid meshes in a fluid-structure interaction system. The required mesh size ratio between the fluid and solid domains is obtained.
A modified finite element procedure for underwater shock analysis
Chan, S.K.
1990-12-31
Using the regular finite element method for analyzing wave propagation problems presents difficulties: (a) The finite element mesh gives spurious reflection of the traveling wave and (b) Since a finite element model has to have a finite boundary, the wave is reflected by the outside boundary. However, for underwater shock problems, only the response of the structure is of major interest, not the behavior of the wave itself, and the shock wave can be assumed to be spherical. By taking advantage of the limited scope of the underwater shock problem, a finite element procedure can be developed that eliminates the above difficulties. This procedure not only can give very accurate solutions but it may also include structural nonlinearities and effect of cavitation.
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.
Nonlinear finite element modeling of THUNDER piezoelectric actuators
NASA Astrophysics Data System (ADS)
Taleghani, Barmac K.; Campbell, Joel F.
1999-06-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.
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.
P-Finite-Element Program For Analysis Of Plates
NASA Technical Reports Server (NTRS)
Smith, James P.
1995-01-01
BUCKY is p-finite-element computer program for highly accurate analysis of structures. Used to analyze buckling, bending, and in-plane stress-and-strain behaviors of plates. Provides elastic-plastic solutions for isotropic plates in states of plane stress, and axisymmetric solution sequence used to treat three-dimensional problems. Computes response of plate to variety of loading and boundary conditions by use of higher-order displacement function in p-finite-element method. Enables user to obtain results more accurate than obtained by use of traditional h-finite elements. Written in FORTRAN 77.
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.
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.
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.
Higher-Order Finite Elements for Computing Thermal Radiation
NASA Technical Reports Server (NTRS)
Gould, Dana C.
2004-01-01
Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other
NASA Astrophysics Data System (ADS)
Lin, Yong; Hu, Jiasheng; Wu, Kenan
2009-08-01
The vector fuzzy control iterative algorithm (VFCIA) is proposed for the design of phase-only sub-wavelength diffractive optical elements (SWDOEs) for beam shaping. The vector diffraction model put forward by Mansuripur is applied to relate the field distributions between the SWDOE plane and the output plane. Fuzzy control theory is used to decide the constraint method for each iterative process of the algorithm. We have designed a SWDOE that transforms a circular flat-top beam to a square irradiance pattern. Computer design results show that the SWDOE designed by the VFCIA can produce better results than the vector iterative algorithm (VIA). And the finite difference time-domain method (FDTD), a rigorous electromagnetic analysis technique, is used to analyze the designed SWDOE for further confirming the validity of the proposed method.
North Atlantic Finite Element Ocean Modeling
NASA Astrophysics Data System (ADS)
Veluthedathekuzhiyil, Praveen
This thesis presents a modified version of the Finite Element Ocean Model (FEOM) developed at Alfred Wegener Institute for Polar and Marine Research (AWI) for the North Atlantic Ocean. A reasonable North Atlantic Ocean simulation is obtained against the observational data sets in a Control simulation (CS) where the surface boundary conditions are relaxed to a climatology. The vertical mixing in the model was tuned to represent convection in the model, also the horizontal mixing and diffusion coefficients to represent the changes in the resolution of the model’s unstructured grid. In addition, the open boundaries in the model are treated with a sponge layer where tracers are relaxed to climatology. The model is then further modified to accept the atmospheric flux forcing at the surface boundary with an added net heat flux correction and freshwater forcing from major rivers that are flowing into the North Atlantic Ocean. The impact of this boundary condition on the simulation results is then analyzed and shows many improvements albeit the drift in tracer properties around the Gulf Stream region remains as that of the CS case. However a comparison of the vertical sections at Cape Desolation and Cape Farewell with the available observational data sets shows many improvements in this simulation compared to that of the CS case. But the freshwater content in the Labrador Sea interior shows a continued drift as that of the CS case with an improvement towards the 10th model year. A detailed analysis of the boundary currents around the Labrador Sea shows the weak offshore transport of freshwater from the West Greenland Current (WGC) as one of the causes. To further improve the model and reasonably represent the boundary currents and associated sub-grid scale eddies in the model, a modified sub-grid scale parameterization based on Gent and McWilliams, (1990) is adopted. The sensitivity of using various approaches in the thickness diffusion parameter ( Kgm) for this
Turcksin, B.; Ragusa, J. C.
2013-07-01
A DSA technique to accelerate the iterative convergence of S{sub n} transport solves is derived for bilinear discontinuous (BLD) finite elements on rectangular grids. The diffusion synthetic acceleration equations are discretized using BLD elements by adapting the Modified Interior Penalty technique, introduced in [4] for triangular grids. The MIP-DSA equations are SPD and thus are solved using a preconditioned CG technique. Fourier analyses and implementation of the technique in a BLD S{sub n} transport code show that the technique is stable is effective. (authors)
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.
Application of Mass Lumped Higher Order Finite Elements
Chen, J.; Strauss, H. R.; Jardin, S. C.; Park, W.; Sugiyama, L. E.; G. Fu; Breslau, J.
2005-11-01
There are many interesting phenomena in extended-MHD such as anisotropic transport, mhd, 2-fluid effects stellarator and hot particles. Any one of them challenges numerical analysts, and researchers are seeking for higher order methods, such as higher order finite difference, higher order finite elements and hp/spectral elements. It is true that these methods give more accurate solution than their linear counterparts. However, numerically they are prohibitively expensive. Here we give a successful solution of this conflict by applying mass lumped higher order finite elements. This type of elements not only keep second/third order accuracy but also scale closely to linear elements by doing mass lumping. This is especially true for second order lump elements. Full M3D and anisotropic transport models are studied.
Error analysis of finite element solutions for postbuckled cylinders
NASA Technical Reports Server (NTRS)
Sistla, Rajaram; Thurston, Gaylen A.
1989-01-01
A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.
Generalized multiscale finite element method. Symmetric interior penalty coupling
NASA Astrophysics Data System (ADS)
Efendiev, Y.; Galvis, J.; Lazarov, R.; Moon, M.; Sarkis, M.
2013-12-01
Motivated by applications to numerical simulations of flows in highly heterogeneous porous media, we develop multiscale finite element methods for second order elliptic equations. We discuss a multiscale model reduction technique in the framework of the discontinuous Galerkin finite element method. We propose two different finite element spaces on the coarse mesh. The first space is based on a local eigenvalue problem that uses an interior weighted L2-norm and a boundary weighted L2-norm for computing the “mass” matrix. The second choice is based on generation of a snapshot space and subsequent selection of a subspace of a reduced dimension. The approximation with these multiscale spaces is based on the discontinuous Galerkin finite element method framework. We investigate the stability and derive error estimates for the methods and further experimentally study their performance on a representative number of numerical examples.
The finite element machine: An experiment in parallel processing
NASA Technical Reports Server (NTRS)
Storaasli, O. O.; Peebles, S. W.; Crockett, T. W.; Knott, J. D.; Adams, L.
1982-01-01
The finite element machine is a prototype computer designed to support parallel solutions to structural analysis problems. The hardware architecture and support software for the machine, initial solution algorithms and test applications, and preliminary results are described.
An enhanced finite element technique for diffuse phase transition
NASA Astrophysics Data System (ADS)
Münch, I.; Krauß, M.
2015-10-01
We propose a finite element technique to enhance phase-field simulations. As adaptive p-method it and can be generally applied to finite element formulations. However, diffuse interfaces have non-linear gradients within regions typically smaller compared to the size of the overall model. Thus, enhanced field interpolation with higher polynomial functions on demand allows for coarser meshing or lower regularization length for the phase transition. Our method preserves continuity of finite elements and is particularly advantageous in the context of parallelized computing. An analytical solution for the evolution of a phase-field variable governed by the Allen-Cahn equation is used to define an error measure and to investigate the proposed method. Several examples demonstrate the capability of this finite element technique.
Validation of high displacement piezoelectric actuator finite element models
NASA Astrophysics Data System (ADS)
Taleghani, Barmac K.
2000-08-01
The paper presents the results obtained by using NASTRAN and ANSYS 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 and important parameters for users of such devices. Therefore, finite element models were used to assess the effects of these parameters. NASTRAN and ANSYS used different methods for modeling piezoelectric effects. In NASTRAN, a thermal analogy was used to represent voltage at nodes as equivalent temperatures, while ANSYS processed the voltage directly using piezoelectric finite elements. The results of finite element models were validated by using the experimental results.
Adaptive Finite-Element Computation In Fracture Mechanics
NASA Technical Reports Server (NTRS)
Min, J. B.; Bass, J. M.; Spradley, L. W.
1995-01-01
Report discusses recent progress in use of solution-adaptive finite-element computational methods to solve two-dimensional problems in linear elastic fracture mechanics. Method also shown extensible to three-dimensional problems.
Scalable, Finite Element Analysis of Electromagnetic Scattering and Radiation
NASA Technical Reports Server (NTRS)
Cwik, T.; Lou, J.; Katz, D.
1997-01-01
In this paper a method for simulating electromagnetic fields scattered from complex objects is reviewed; namely, an unstructured finite element code that does not use traditional mesh partitioning algorithms.
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.
Comparison of different precondtioners for nonsymmtric finite volume element methods
Mishev, I.D.
1996-12-31
We consider a few different preconditioners for the linear systems arising from the discretization of 3-D convection-diffusion problems with the finite volume element method. Their theoretical and computational convergence rates are compared and discussed.
NASA Astrophysics Data System (ADS)
Bing, Zhou; Greenhalgh, S. A.
2001-06-01
The finite element method is a powerful tool for 3-D DC resistivity modelling and inversion. The solution accuracy and computational efficiency are critical factors in using the method in 3-D resistivity imaging. This paper investigates the solution accuracy and the computational efficiency of two common element-type schemes: trilinear interpolation within a regular 8-node solid parallelepiped, and linear interpolations within six tetrahedral bricks within the same 8-node solid block. Four iterative solvers based on the pre-conditioned conjugate gradient method (SCG, TRIDCG, SORCG and ICCG), and one elimination solver called the banded Choleski factorization are employed for the solutions. The comparisons of the element schemes and solvers were made by means of numerical experiments using three synthetic models. The results show that the tetrahedron element scheme is far superior to the parallelepiped element scheme, both in accuracy and computational efficiency. The tetrahedron element scheme may save 43 per cent storage for an iterative solver, and achieve an accuracy of the maximum relative error of <1 per cent with an appropriate element size. The two iterative solvers, SORCG and ICCG, are suitable options for 3-D resistivity computations on a PC, and both perform comparably in terms of convergence speed in the two element schemes. ICCG achieves the best convergence rate, but nearly doubles the total storage size of the computation. Simple programming codes for the two iterative solvers are presented. We also show that a fine grid, which doubles the density of a coarse grid, will require at least 27=128 times as much computing time when using the banded Choleski factorization. Such an increase, especially for 3-D resistivity inversion, should be compared with SORCG and ICCG solvers in order to find the computationally most efficient method when dealing with a large number of electrodes.
Finite element analysis of vibration and damping of laminated composites
NASA Astrophysics Data System (ADS)
Rikards, Rolands
Simple finite elements are used to form a special laminated beam and plate superelements excluding all degrees of freedom in the nodes of the middle layer, and the finite element analysis of this structure is performed. To estimate damping of structures, modal loss factors are calculated, using two methods: the 'exact' method of complex eigenvalues and the approximate energy method. It was found that both methods give satisfactory results. However, the energy method needs less computer time than the exact method.
Mathematical aspects of finite element methods for incompressible viscous flows
NASA Technical Reports Server (NTRS)
Gunzburger, M. D.
1986-01-01
Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.
Evaluation of a hybrid, anisotropic, multilayered, quadrilateral finite element
NASA Technical Reports Server (NTRS)
Robinson, J. C.; Blackburn, C. L.
1978-01-01
A multilayered finite element with bending-extensional coupling is evaluated for: (1) buckling of general laminated plates; (2) thermal stresses of laminated plates cured at elevated temperatures; (3) displacements of a bimetallic beam; and (4) displacement and stresses of a single-cell box beam with warped cover panels. Also, displacements and stresses for flat and spherical orthotropic and anisotropic segments are compared with results from higher order plate and shell finite-element analyses.
Finite element analysis of a composite wheelchair wheel design
NASA Technical Reports Server (NTRS)
Ortega, Rene
1994-01-01
The finite element analysis of a composite wheelchair wheel design is presented. The design is the result of a technology utilization request. The designer's intent is to soften the riding feeling by incorporating a mechanism attaching the wheel rim to the spokes that would allow considerable deflection upon compressive loads. A finite element analysis was conducted to verify proper structural function. Displacement and stress results are presented and conclusions are provided.
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.
Simulation of two-dimensional waterflooding using mixed finite elements
Chavent, G.; Jaffre, J.; Cohen, G.; Dupuy, M.; Dieste, I.
1982-01-01
A new method for the simulation of incompressible diphasic flows in two dimensions is presented, the distinctive features of which are: (1) reformation of the basic equation and specific choices of the finite element approximation of the same; (11) use of a mixed finite elements method, approximating both scalar and vector functions. Several test examples are shown, including gravity and capillary effects. The use of discontinuous basis functions proved successful for an accurate representation of sharp fronts. 16 refs.
Integration of geometric modeling and advanced finite element preprocessing
NASA Technical Reports Server (NTRS)
Shephard, Mark S.; Finnigan, Peter M.
1987-01-01
The structure to a geometry based finite element preprocessing system is presented. The key features of the system are the use of geometric operators to support all geometric calculations required for analysis model generation, and the use of a hierarchic boundary based data structure for the major data sets within the system. The approach presented can support the finite element modeling procedures used today as well as the fully automated procedures under development.
Azimuthally-dependent Finite Element Solution to the Cylindrical Resonator
NASA Technical Reports Server (NTRS)
Osegueda, R.; Pierluissi, J.; Gil, L.; Revilla, A.; Villalva, G.; Dick, G.; Wang, D. SantiagoR.
1994-01-01
The cylindrical cavity resonator loaded with an anisotropic dielectric is analyzed as a two-dimensional problem using a finite element approach that assumes sinusoidal dependence in azimuth. This methodology allows the first finite element treatment of the technically important case of a resonator containing a sapphire element with a cylindrically aligned c axis. Second order trial functions together with quadrilateral elements are adopted in the calculations. The method was validated through comparisons with the analytical solutions for the hollow metal cavity and a coaxial cavity, as well as through measurements on a shielded sapphire resonator.
Finite element thermal analysis of convectively-cooled aircraft structures
NASA Technical Reports Server (NTRS)
Wieting, A. R.; Thornton, E. A.
1981-01-01
The design complexity and size of convectively-cooled engine and airframe structures for hypersonic transports necessitate the use of large general purpose computer programs for both thermal and structural analyses. Generally thermal analyses are based on the lumped-parameter finite difference technique, and structural analyses are based on the finite element technique. Differences in these techniques make it difficult to achieve an efficient interface. It appears, therefore, desirable to conduct an integrated analysis based on a common technique. A summary is provided of efforts by NASA concerned with the development of an integrated thermal structural analysis capability using the finite element method. Particular attention is given to the development of conduction/forced-convection finite element methodology and applications which illustrate the capabilities of the developed concepts.
Recent developments in finite element analysis for transonic airfoils
NASA Technical Reports Server (NTRS)
Hafez, M. M.; Murman, E. M.
1979-01-01
The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.
Dynamical observer for a flexible beam via finite element approximations
NASA Technical Reports Server (NTRS)
Manitius, Andre; Xia, Hong-Xing
1994-01-01
The purpose of this view-graph presentation is a computational investigation of the closed-loop output feedback control of a Euler-Bernoulli beam based on finite element approximation. The observer is part of the classical observer plus state feedback control, but it is finite-dimensional. In the theoretical work on the subject it is assumed (and sometimes proved) that increasing the number of finite elements will improve accuracy of the control. In applications, this may be difficult to achieve because of numerical problems. The main difficulty in computing the observer and simulating its work is the presence of high frequency eigenvalues in the finite-element model and poor numerical conditioning of some of the system matrices (e.g. poor observability properties) when the dimension of the approximating system increases. This work dealt with some of these difficulties.
Two-dimensional finite element neutron diffusion analysis using hierarchic shape functions
Carpenter, D.C.
1997-04-01
Recent advances have been made in the use of p-type finite element method (FEM) for structural and fluid dynamics problems that hold promise for reactor physics problems. These advances include using hierarchic shape functions, element-by-element iterative solvers and more powerful mapping techniques. Use of the hierarchic shape functions allows greater flexibility and efficiency in implementing energy-dependent flux expansions and incorporating localized refinement of the solution space. The irregular matrices generated by the p-type FEM can be solved efficiently using element-by-element conjugate gradient iterative solvers. These solvers do not require storage of either the global or local stiffness matrices and can be highly vectorized. Mapping techniques based on blending function interpolation allow exact representation of curved boundaries using coarse element grids. These features were implemented in a developmental two-dimensional neutron diffusion program based on the use of hierarchic shape functions (FEM2DH). Several aspects in the effective use of p-type analysis were explored. Two choices of elemental preconditioning were examined--the proper selection of the polynomial shape functions and the proper number of functions to use. Of the five shape function polynomials tested, the integral Legendre functions were the most effective. The serendipity set of functions is preferable over the full tensor product set. Two global preconditioners were also examined--simple diagonal and incomplete Cholesky. The full effectiveness of the finite element methodology was demonstrated on a two-region, two-group cylindrical problem but solved in the x-y coordinate space, using a non-structured element grid. The exact, analytic eigenvalue solution was achieved with FEM2DH using various combinations of element grids and flux expansions.