Downloads & Free Reading Options - Results

Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.

A number of finite element methods for the analysis of nonlinear problems in rubber elasticity are outlined. Several different finite element schemes are discussed. These include the augmented Lagrangian method, continuation or incremental loading methods, and associated Riks-type methods which have the capability of incorporating limit point behavior and bifurcations. Algorithms for the analysis of limit point behavior and bifurcations are described and the results of several numerical experiments are presented. In addition, a brief survey of some recent work on modelling contact and friction in elasticity problems is given. These results pertain to the use of new nonlocal and nonlinear friction laws.

In this article we study adaptive finite element methods (AFEM) with inexact solvers for a class of semilinear elliptic interface problems. We are particularly interested in nonlinear problems with discontinuous diffusion coefficients, such as the nonlinear Poisson-Boltzmann equation and its regularizations. The algorithm we study consists of the standard SOLVE-ESTIMATE-MARK-REFINE procedure common to many adaptive finite element algorithms, but where the SOLVE step involves only a full solve on the coarsest level, and the remaining levels involve only single Newton updates to the previous approximate solution. We summarize a recently developed AFEM convergence theory for inexact solvers, and present a sequence of numerical experiments that give evidence that the theory does in fact predict the contraction properties of AFEM with inexact solvers. The various routines used are all designed to maintain a linear-time computational complexity.

In this paper, we propose a multiscale empirical interpolation method for solving nonlinear multiscale partial differential equations. The proposed method combines empirical interpolation techniques and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM). To solve nonlinear equations, the GMsFEM is used to represent the solution on a coarse grid with multiscale basis functions computed offline. Computing the GMsFEM solution involves calculating the residuals on the fine grid. We use empirical interpolation concepts to evaluate the residuals and the Jacobians of the multiscale system with a computational cost which is proportional to the coarse scale problem rather than the fully-resolved fine scale one. Empirical interpolation methods use basis functions and an inexpensive inversion which are computed in the offline stage for finding the coefficients in the expansion based on a limited number of nonlinear function evaluations. The proposed multiscale empirical interpolation techniques: (1) divide computing the nonlinear function into coarse regions; (2) evaluate contributions of nonlinear functions in each coarse region taking advantage of a reduced-order representation of the solution; and (3) introduce multiscale proper-orthogonal-decomposition techniques to find appropriate interpolation vectors. We demonstrate the effectiveness of the proposed methods on several examples of nonlinear multiscale PDEs that are solved with Newton's methods and fully-implicit time marching schemes. Our numerical results show that the proposed methods provide a robust framework for solving nonlinear multiscale PDEs on a coarse grid with bounded error.

Book Source: Digital Library of India Item 2015.193316 dc.contributor.author: Amiya Kumar Pani dc.date.accessioned: 2015-07-08T02:30:55Z dc.date.available: 2015-07-08T02:30:55Z dc.date.digitalpublicationdate: 2005-08-27 dc.identifier.barcode: 1990010086728 dc.identifier.origpath: /rawdataupload/upload/0086/728 dc.identifier.copyno: 1 dc.identifier.uri: http://www.new.dli.ernet.in/handle/2015/193316 dc.description.scannerno: 14 dc.description.scanningcentre: IIIT, Allahabad dc.description.main: 1 dc.description.tagged: 0 dc.description.totalpages: 180 dc.format.mimetype: application/pdf dc.language.iso: English dc.publisher: Indian Institute Of Technology Kanpur dc.rights: Out_of_copyright dc.source.library: Indian Institute Of Technology Kanpur dc.subject.classification: Technology dc.subject.classification: Engineering. Technology In General dc.subject.classification: Mechanical Engineering In General. Nuclear Technology. Electrical Engineering. Machinery dc.title: Finite Element Methods For Nonlinear Free Boundary Problems

Book Source: Digital Library of India Item 2015.193316 dc.contributor.author: Amiya Kumar Pani dc.date.accessioned: 2015-07-08T02:30:55Z dc.date.available: 2015-07-08T02:30:55Z dc.date.digitalpublicationdate: 2005-08-27 dc.identifier.barcode: 1990010086728 dc.identifier.origpath: /rawdataupload/upload/0086/728 dc.identifier.copyno: 1 dc.identifier.uri: http://www.new.dli.ernet.in/handle/2015/193316 dc.description.scannerno: 14 dc.description.scanningcentre: IIIT, Allahabad dc.description.main: 1 dc.description.tagged: 0 dc.description.totalpages: 180 dc.format.mimetype: application/pdf dc.language.iso: English dc.publisher: Indian Institute Of Technology Kanpur dc.rights: Out_of_copyright dc.source.library: Indian Institute Of Technology Kanpur dc.subject.classification: Technology dc.subject.classification: Engineering. Technology In General dc.subject.classification: Mechanical Engineering In General. Nuclear Technology. Electrical Engineering. Machinery dc.title: Finite Element Methods For Nonlinear Free Boundary Problems

No Abstract Available

This paper studies mixed finite element approximations of the viscosity solution to the Dirichlet problem for the fully nonlinear Monge-Amp\`ere equation $\det(D^2u^0)=f$ based on the vanishing moment method which was proposed recently by the authors in \cite{Feng2}. In this approach, the second order fully nonlinear Monge-Amp\`ere equation is approximated by the fourth order quasilinear equation $-\epsilon\Delta^2 u^\epsilon + \det{D^2u^\epsilon} =f$. It was proved in \cite{Feng1} that the solution $u^\epsilon$ converges to the unique convex viscosity solution $u^0$ of the Dirichlet problem for the Monge-Amp\`ere equation. This result then opens a door for constructing convergent finite element methods for the fully nonlinear second order equations, a task which has been impracticable before. The goal of this paper is threefold. First, we develop a family of Hermann-Miyoshi type mixed finite element methods for approximating the solution $u^\epsilon$ of the regularized fourth order problem, which computes simultaneously $u^\vepsi$ and the moment tensor $\sigma^\vepsi:=D^2u^\epsilon$. Second, we derive error estimates, which track explicitly the dependence of the error constants on the parameter $\vepsi$, for the errors $u^\epsilon-u^\epsilon_h$ and $\sigma^\vepsi-\sigma_h^\vepsi$. Finally, we present a detailed numerical study on the rates of convergence in terms of powers of $\vepsi$ for the error $u^0-u_h^\vepsi$ and $\sigma^\vepsi-\sigma_h^\vepsi$, and numerically examine what is the "best" mesh size $h$ in relation to $\vepsi$ in order to achieve these rates.

In this paper we use the GeneralizedMultiscale Finite ElementMethod (GMsFEM) framework, introduced in [20], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation. With this convention, we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin (CG) or discontinuous Galerkin (DG) global formulations. Here, we use Symmetric Interior Penalty Discontinuous Galerkin approach. Both methods yield a predictable error decline that depends on the respective coarse space dimension, and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples.

A finite element method, developed for static analysis of composite plates, was enhanced to be used with sandwich plates. The theory considers geometric nonlinearity and transverse shear. Furthermore, a new postprocessor was written to check for initial failure using the maximum stress criteria. It also includes a procedure for evaluating transverse normal stresses by enforcing equilibrium. The code modifications for sandwiches were verified by by comparing finite element solutions to closed-form sandwich theories. Both methods showed good correlation. In addition, comparisons between one type of composite plate and a similar sandwich plate found that the sandwich had better specific stiffness for various geometries. Finally, quasi-static models of low-velocity impacts on sandwich plates succeeded in predicting incipient damage due to core crushing in the impact region, but they were inconsistent in estimating the actual load levels for initial failure, based on experiments from other research.

A finite element formulation and algorithm for the nonlinear analysis of the large deflection, materially nonlinear response of impulsively loaded shells is presented. A unique feature of this algorithm is the use of a bilinear four node quadrilateral element with single point quadrature and a simple hourglass control which is orthogonal to rigid body modes on an element level and does not compromise the consistency of the equations. The geometric nonlinearities are treated by using a corotational description wherein a coordinate system that rotates with the material is embedded at the integration point; thus the algorithm is directly applicable to anisotropic materials without any corrections for frame invariance of material property tensors. This algorithm can treat about 200 element-time-steps per CPU second on a CYBER 170/730 computer in the explicit time integration mode. Numerous results are presented for both elastic and elastic-plastic problems with large strains that show that the method in most cases is comparable in accuracy with an earlier version of this algorithm employing a cubic triangular plate-shell element, but substantially faster. (Author)

We consider the mixed formulation of the linear elasticity problem. We identify stable families of mixed finite element spaces for both curvilinear quadrilateral and triangular elements. We give conditions under which these elements converge at asymptotically optimal (or near-optimal) rates, both in terms of the mesh width h and polynomial degree p. We validate these findings computationally. Finally, we identify a mixed formulation for the geometrically nonlinear elasticity problem and experimentally investigate the stability of our mixed finite element spaces for this problem.