Sparse tensor product finite element method for nonlinear multiscale variational inequalities of monotone type

2019 ◽  
Vol 40 (3) ◽  
pp. 1875-1907
Author(s):  
Wee Chin Tan ◽  
Viet Ha Hoang

Abstract We study an essentially optimal finite element (FE) method for locally periodic nonlinear multiscale variational inequalities of monotone type in a domain $D\subset{\mathbb{R}}^d$ that depend on a macroscopic and $n$ microscopic scales. The scales are separable. Using multiscale convergence we deduce a multiscale homogenized variational inequality in a tensorized domain in the high-dimensional space ${\mathbb R}^{(n+1)d}$. Given sufficient regularity on the solution the sparse tensor product FE method is developed for this problem, which attains an essentially equal (i.e., it differs by only a logarithmic factor) level of accuracy to that of the full tensor product FE method, but requires an essentially optimal number of degrees of freedom which is equal to that for solving a problem in ${{\mathbb{R}}}^d$ apart from a logarithmic factor. For two-scale problems we deduce a new homogenization error for the nonlinear monotone variational inequality. A numerical corrector is then constructed with an explicit error in terms of the homogenization and the FE errors. For general multiscale problems we deduce a numerical corrector from the FE solution of the multiscale homogenized problem, but without an explicit error as such a homogenization error is not available.

2014 ◽  
Vol 34 (4) ◽  
pp. 323-332 ◽  
Author(s):  
Gang Liu ◽  
Wei Tang ◽  
Ying-Lin Ke ◽  
Qing-Liang Chen ◽  
Yunbo Bi

Purpose – The purpose of this paper is to propose a new model for optimizing pre-joining processes quickly and accurately, guiding workers to standardized operations. For the automatic riveting in panel assemblies, the traditional approach of determination of pre-joining processes entirely rests on the experience of workers, which leads to the improper number, location and sequence of pre-joining, the low quality stability and the high repair rate in most cases. Design/methodology/approach – The clearances computation with the complete finite element model for every process combination is time-consuming. Therefore a fast pre-joining processes optimization model (FPPOM) is proposed. This model treats both the measured initial clearances and the stiffness matrices of key points of panels as an input; considers the permissive clearances as an evaluation criterion; regards the optimal number, location and sequence as an objective; and takes the neighborhood-search-based adaptive genetic algorithm as a solution. Findings – A comparison between the FPPOM and complete finite element model with clearances (CFEMC) was made in practice. Further, the results indicate that running the FPPOM is time-saving by >90 per cent compared with the CFEMC. Practical implications – This paper provides practical insights into realizing the pre-joining processes optimization quickly. Originality/value – This paper is the first to propose the FPPOM, which could simplify the processes, reduce the degrees of freedom of nodes and conduct the manufacturers to standardized manipulations.


Author(s):  
Lassaad Ben Fekih ◽  
Georges Kouroussis ◽  
David Wattiaux ◽  
Olivier Verlinden ◽  
Christophe De Fruytier

An approach is proposed to identify the modal properties of a subsystem made up of an arbitrary chosen inner module of embedded space equipment. An experimental modal analysis was carried out along the equipment transverse direction with references taken onto its outer housing. In parallel, a numerical model using the finite element (FE) method was developed to correlate with the measured results. A static Guyan reduction has led to a set of master degrees of freedom in which the experimental mode shapes were expanded. An updating technique consisting in minimizing the dynamic residual induced by the FE model and the measurements has been investigated. A last verification has consisted in solving the numeric model composed of the new mass and stiffness matrices obtained by means of a minimization of the error in the constitutive equation method.


2021 ◽  
Vol 37 (3) ◽  
pp. 477-487
Author(s):  
MONDAY OGUDU NNAKWE ◽  
◽  
" JERRY N." EZEORA ◽  

In this paper, using a sunny generalized non-expansive retraction which is different from the metric projection and generalized metric projection in Banach spaces, we present a retractive iterative algorithm of Krasnosel’skii-type, whose sequence approximates a common solution of a mono-variational inequality of a finite family of η-strongly-pseudo-monotone-type maps and fixed points of a countable family of generalized non-expansive-type maps. Furthermore, some new results relevant to the study are also presented. Finally, the theorem proved complements, improves and extends some important related recent results in the literature.


Author(s):  
M. H. Refaat ◽  
S. A. Meguid

Abstract Current solution schemes of variational inequalities arising in frictional contact problems adopt penalty and regularization techniques. The convergence and accuracy of these schemes are governed by user-defined parameters. To overcome the difficulties associated with the ad hoc use of such parameters, the variational inequality of the general frictional contact problem is treated in this paper using mathematical programming. A new non-differential optimization (NDO) technique, in association with quadratic programming, is used to treat the resulting variational inequalities.


Author(s):  
M. Logesh ◽  
S. Palani ◽  
S. Shanmugan ◽  
M. Selvam ◽  
K.A. Harish

Finite element (FE) method is commonly used to study cracks in structures. In this paper, J-integral method is applied over FE model of a cracked body to determine stress intensity factor (SIF) in the domain of linear elastic fracture mechanics (LEFM). This paper formulates the J-integral methodology for 2D FE model using a coarse mesh with less degrees of freedom. Two cases , a finite plate with edge cracks and a normal crack growth in fiber metal laminated plate, are demonstrated. Numerical implementation and mesh refinement issues to maintain path independent J-integral values are explored.


2019 ◽  
Vol 91 (1) ◽  
pp. 465-474
Author(s):  
Chunfang Meng ◽  
Chen Gu ◽  
Bradford Hager

Abstract We present a fundamental solution‐based finite‐element (FE) method to homogenize heterogeneous elastic medium, that is, fault zone, under static, and dynamic loading. This method incorporates Eshelby’s strain perturbation into FE weak forms. The resulting numerical model implicitly considers the existence of inhomogeneity bodies within each element, without introducing additional degrees of freedom. The new method is implemented within an open‐source FE package that is applicable to alternating seismic and aseismic cycles. To demonstrate this method, we modify a dynamic fault‐slip problem, hosted at Southern California Earthquake Center (SCEC), by introducing a fault zone that contains different microstructures than the host matrix. The preliminary results suggest that the fault‐zone microstructure orientation has effects on fault slip, seismic arrivals and waveform frequency contents.


1966 ◽  
Vol 25 ◽  
pp. 373
Author(s):  
Y. Kozai

The motion of an artificial satellite around the Moon is much more complicated than that around the Earth, since the shape of the Moon is a triaxial ellipsoid and the effect of the Earth on the motion is very important even for a very close satellite.The differential equations of motion of the satellite are written in canonical form of three degrees of freedom with time depending Hamiltonian. By eliminating short-periodic terms depending on the mean longitude of the satellite and by assuming that the Earth is moving on the lunar equator, however, the equations are reduced to those of two degrees of freedom with an energy integral.Since the mean motion of the Earth around the Moon is more rapid than the secular motion of the argument of pericentre of the satellite by a factor of one order, the terms depending on the longitude of the Earth can be eliminated, and the degree of freedom is reduced to one.Then the motion can be discussed by drawing equi-energy curves in two-dimensional space. According to these figures satellites with high inclination have large possibilities of falling down to the lunar surface even if the initial eccentricities are very small.The principal properties of the motion are not changed even if plausible values ofJ3andJ4of the Moon are included.This paper has been published in Publ. astr. Soc.Japan15, 301, 1963.


2020 ◽  
Vol 1 (1) ◽  
pp. 93-102
Author(s):  
Carsten Strzalka ◽  
◽  
Manfred Zehn ◽  

For the analysis of structural components, the finite element method (FEM) has become the most widely applied tool for numerical stress- and subsequent durability analyses. In industrial application advanced FE-models result in high numbers of degrees of freedom, making dynamic analyses time-consuming and expensive. As detailed finite element models are necessary for accurate stress results, the resulting data and connected numerical effort from dynamic stress analysis can be high. For the reduction of that effort, sophisticated methods have been developed to limit numerical calculations and processing of data to only small fractions of the global model. Therefore, detailed knowledge of the position of a component’s highly stressed areas is of great advantage for any present or subsequent analysis steps. In this paper an efficient method for the a priori detection of highly stressed areas of force-excited components is presented, based on modal stress superposition. As the component’s dynamic response and corresponding stress is always a function of its excitation, special attention is paid to the influence of the loading position. Based on the frequency domain solution of the modally decoupled equations of motion, a coefficient for a priori weighted superposition of modal von Mises stress fields is developed and validated on a simply supported cantilever beam structure with variable loading positions. The proposed approach is then applied to a simplified industrial model of a twist beam rear axle.


Author(s):  
Wen Zhang ◽  
Wenliang Wang ◽  
Hao Wang ◽  
Jiong Tang

A method for dynamic analysis of flexible bladed-disk/shaft coupled systems is presented in this paper. Being independant substructures first, the rigid-disk/shaft and each of the bladed-disk assemblies are analyzed separately in a centrifugal force field by means of the finite element method. Then through a modal synthesis approach the equation of motion for the integral system is derived. In the vibration analysis of the rotating bladed-disk substructure, the geometrically nonlinear deformation is taken into account and the rotationally periodic symmetry is utilized to condense the degrees of freedom into one sector. The final equation of motion for the coupled system involves the degrees of freedom of the shaft and those of only one sector of each of the bladed-disks, thereby reducing the computer storage. Some computational and experimental results are given.


Mathematics ◽  
2021 ◽  
Vol 9 (3) ◽  
pp. 266 ◽  
Author(s):  
Savin Treanţă

A new class of differential variational inequalities (DVIs), governed by a variational inequality and an evolution equation formulated in infinite-dimensional spaces, is investigated in this paper. More precisely, based on Browder’s result, optimal control theory, measurability of set-valued mappings and the theory of semigroups, we establish that the solution set of DVI is nonempty and compact. In addition, the theoretical developments are accompanied by an application to differential Nash games.


Sign in / Sign up

Export Citation Format

Share Document