An extended edge-based smoothed discrete shear gap method for free vibration analysis of cracked Reissner–Mindlin plate

An extended edge-based smoothed discrete shear gap method (XES-DSG3) is proposed for free vibration analysis of cracked Reissner–Mindlin plate by implementing the edge-based strain smoothing operation into the discrete shear gap-based extended finiteelement method (XFEM-DSG3). In present method, the strain smoothing operation is implemented into the bending strain gradient matrices, in which the enriched functions are included. Then, the derivatives of element shape functions and derivatives of crack-tip singular enriched functions are not required in the computation. The calculation of element matrices is performed over the smoothing domains which are associated with edges of elements. The transverse shear locking of Reissner–Mindlin plate can be avoided by using the integration of discrete shear gap (DSG) method. Several numerical examples are investigated to illustrate the accuracy of XES-DSG3 for the free vibration analysis of cracked Reissner–Mindlin plate. Moreover, numerical results show that the present method is insensitive to mesh distortion and it is more stable than the pervious XFEM-DSG3.     

3D hp-adaptive finite element simulations of bend,step, and magic-T electromagnetic waveguide structures

Metallic rectangular waveguides are often the preferred choice on telecommunication systems and medical equipment working on the upper microwave and millimeter wave frequency bands due to the simplicity of its geometry, low losses, and the capacity to handle high powers. Waveguide translational symmetry is interrupted by the unavoidable presence of bends, transitions, and junctions, among others. This paper employs a 3D hp self-adaptive grid-refinement finite element strategy for the solution of these relevant electromagnetic waveguide problems. These structures often incorporate dielectrics, metallic screws, round corners, and so on, which may facilitate its construction or improve its design, but significantly difficult its modeling when employing semi-analytical techniques. The hp-adaptive finite element method enables accurate modeling of these structures even in the presence of complex materials and geometries. Numerical results demonstrate the suitability of the hp-adaptive method for modeling these waveguide structures, delivering errors below 0.5% with a limited number of unknowns. Solutions of waveguide problems obtained with the self-adaptive hp-FEM are of comparable accuracy to those obtained with semi-analytical techniques such as the Mode Matching method, for problems where the latest methods can be applied. At the same time, the hp-adaptive FEM enables accurate modeling of more complex waveguide structures.     

Optimal balance between mass and smoothed stiffness in simulation of acoustic problems

The Finite Element Method (FEM) is known to behave overly-stiff, which leads to an imbalance between the mass and stiffness matrices within discretized systems. In this work, for the first time, a model is developed that provides optimal balance between discretized mass and smoothed stiffness—the mass-redistributed alpha finite element method (MR-αFEM). This new method improves on the computational efficiency of the FEM and Smoothed Finite Element Methods (S-FEM). The rigorous research conducted ensures that stiffness with the parameter, α, optimally matches the mass with a flexible integration point, q. The optimal balance system significantly reduces the dispersion error of acoustic problems, including those of single and multi-fluids in both time and frequency domains. The excellent properties of the proposed MR-αFEM are validated using theoretical analyses and numerical examples.     

Analysis of laminated composite plates using higher-order shear deformation plate theory and node-based smoothed discrete shear gap method

This paper presents a novel finite element formulation for static, free vibration and buckling analyses of laminated composite plates. The idea relies on a combination of node-based smoothing discrete shear gap method with the higher-order shear deformation plate theory (HSDT) to give a so-called NS-DSG3 element. The higher-order shear deformation plate theory (HSDT) is introduced in the present method to remove the shear correction factors and improve the accuracy of transverse shear stresses. The formulation uses only linear approximations and its implementation into finite element programs is quite simple and efficient. The numerical examples demonstrated that the present element is free of shear locking and shows high reliability and accuracy compared to other published solutions in the literature.     

A simple proof for a theorem of Luxemburg and Zaanen

In this paper a simple proof for the following theorem, due to Luxemburg and Zaanen is given: an Archimedean vector lattice A is Dedekind σ-complete if and only if A has the principal projection property and A is uniformly complete. As an application, we give a new and short proof for the following version of Freudenthal's spectral theorem: let A be a uniformly complete vector lattice with the principal projection property and let 0<u∈A. For any element w in A such that 0?w?u there exists a sequence in A which satisfies , where each element sn is of the form , with real numbers α₁,…,αk such that 0?αi?1 (i=1,…,k) and mutually disjoint components p₁,…,pk of u.     

On Tauberian theorems for (A)(C, α) summability method

In this paper we introduce some Tauberian conditions for the (A)(C, α) summability method. These results extend and generalize some of the classical Tauberian theorems for the Abel summability method.     

Solution of some plate bending problems using the boundary element method

Some fundamental aspects of the boundary element method of the Kirchhoff theory of thin plate flexure are given. The direct boundary integral equation method with higher conforming properties (using first-order Hermitian interpolation for plate displacement ω, and zero-order Hermitian interpolation for angle of rotation θ, the moment M andthe equivalent shear V) are used for several computational examples. They are: square plate with simply-supported or clamped edges, the same square plate with square central opening and the cantilevered triangular plates. The results of computation as compared with some known experimental and theoritical results showed that the numerical schemes seemed to be satisfactory for the practical applications.     

无锁定的退化的等参壳元

构造了一个无锁定的8节点退化的等参壳元。在自然坐标系中,对横向剪应变项插值修正使单元满足必要的模式(平移、转动和纯弯曲)以消除"剪切锁定";在局部坐标系中,对薄膜应变插值修正以消除"薄膜锁定"。这样构造的8节点单元的刚度矩阵具有正确的秩,同时具有恰当的零特征值和相应的刚体位移模式。这种单元对大跨-厚比情形既无"剪切锁定"又无"薄膜锁定",无虚假的零能模式和机构出现,可用于厚壳和薄壳。     

Analysis of multiscale finite element method for the stationary Navier-Stokes equations

In this work, a multiscale finite element method is proposed for the stationary incompressible Navier-Stokes equations. And the inf-sup stability of the method for the P₁/P₁ triangular element is established. The optimal error estimates are obtained.     

Stabilized multiscale finite element method for the stationary Navier-Stokes equations

In the paper, a stabilized multiscale finite element method for the stationary incompressible Navier-Stokes equations is considered. The method is a Petrov-Galerkin approach based on the multiscale enrichment of the standard polynomial space enriched with the unusual bubble functions which no longer vanish on every element boundary for the velocity space. The stability of the P₁-P₀ triangular element (or the Q₁-P₀ quadrilateral element) is established. And the optimal error estimates of the stabilized multiscale finite element method for the stationary Navier-Stokes equations are obtained.     

Analysis of a 3‐node mixed‐shear‐projected triangular element method for Reissner–Mindlin plates

This article analyses an existing 3‐node hybrid triangular element, called MiSP3, for Reissner–Mindlin plates which behaves robustly in numerical benchmark tests (Ayad, Dhatt, and Batoz, Int J Numer Method Eng 42 (1998), 1149–1179). Based on Hellinger‐Reissner variational principle and the mixed shear interpolation/projection technique of MITC family, the MiSP3 element uses continuous piecewise linear polynomials for the approximations of displacements and a piecewise‐independent equilibrium mode for the approximations of bending moments/shear stresses. Due to local elimination of the parameters of moments/stresses, the element is almost of the same computational cost as the conforming linear triangular displacement element. We derive uniform stability and convergence results with respect to the plate thickness. The main tools of our analysis are the self‐equilibrium relation of the moments/stresses approximations, the properties of the mixed shear interpolation and the discrete Helmholtz decomposition of the shear stress approximation. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 241–258, 2017     

Strong convergence of an iterative method for nonexpansive and accretive operators

Let X be a Banach space and A an m-accretive operator with a zero. Consider the iterative method that generates the sequence {xn} by the algorithm xn+1=αnu+(1−αn)Jrnxn, where {αn} and {rn} are two sequences satisfying certain conditions, and Jr denotes the resolvent ⁻¹(I+rA) for r>0. Strong convergence of the algorithm {xn} is proved assuming X either has a weakly continuous duality map or is uniformly smooth.     

On degrees of freedom of certain conservative turbulence models for the Navier-Stokes equations

We study the degrees of freedom of several conservative computational turbulence models that are derived via a non-dissipative regularizations of the Navier-Stokes equations. For the Navier-Stokes-α, the Leray-α and the Navier-Stokes-ω equations we prove that the longtime behavior of their respective solutions is completely determined by a finite set of grid values and by a finite set of Fourier modes. For each turbulence model the number of determining nodes and of determining modes is estimated in terms of flow parameters, such as viscosity, smoothing length, forcing and domain size. These estimates are global as they do not depend on an individual solution.     

三角形REISSNER-MINDLIN板元

本文提出构造无自锁现象的Reissuer－Mindlin板元的一个一般性方法.此方法将剪切应变用它的适当的插值多项式代替,当板厚趋于零时这对应于插值点的Kirchhoff条件,因而单元无自锁现象.根据这种方法我们构造两个三角形元──一个3节点元和一个6节点元,并给出数值结果.     

A class of semigroups regularized in space and time

We consider α-times integrated C-regularized semigroups, which are a hybrid between semigroups regularized in space (C-semigroups) and integrated semigroups regularized in time. We study the basic properties of these objects, also in absence of exponential boundedness. We discuss their generators and establish an equivalence theorem between existence of integrated regularized semigroups and well-posedness of certain Cauchy problems. We investigate the effect of smoothing regularized semigroups by fractional integration.     

Free vibration analysis of a rigid bar supported by arbitrary elastic beams

For convenience, a two-node conventional elastic beam element (C beam element) with the displacements of its 2nd node replaced by those of center of gravity (c.g.) of the joined rigid bar is called the modified beam element (M beam element). The objective of this paper is to present a modified finite element method (modified FEM) such that the free vibration characteristics of a rigid bar supported by a number of elastic beams can be easily determined. First of all, the displacements for the 2nd node of a C beam element joined with the rigid bar are determined in terms of those for the c.g. of the joined rigid bar to establish the M beam element. Next, the mass and stiffness matrices for the M beam element are derived based on the displacements for the 1st node of the C beam element and those for the c.g. of the joined rigid bar. Then, the overall property matrices of the entire unconstrained vibrating system (i.e. a rigid bar supported by a number of elastic beams) can be determined by using the assembly technique of the conventional FEM and considering the effects of lumped mass and rotary inertia of the rigid bar. Finally, the boundary (supporting) conditions are imposed to produce the effective property matrices of the constrained vibrating system and then the free vibration characteristics are determined with the standard approach. In order to confirm the presented theory and the developed computer program, the rigid bar is modeled by a number of C beam elements with bigger Young’s modulus (E_R) and the conventional FEM is used to determine the natural frequencies and associated mode shapes of the vibrating system. It is found that the latter will converge to the corresponding ones obtained from the presented modified FEM when the magnitude of E_R increases to certain values.     

On fractional powers of generators of fractional resolvent families

We show that if −A generates a bounded α-times resolvent family for some α∈(0,2], then −Aβ generates an analytic γ-times resolvent family for and γ∈(0,2). And a generalized subordination principle is derived. In particular, if −A generates a bounded α-times resolvent family for some α∈(1,2], then −A1/α generates an analytic C₀-semigroup. Such relations are applied to study the solutions of Cauchy problems of fractional order and first order.     

Analysis of adiabatic shear bands in thermo-elasto-viscoplastic materials by using piece-wise discontinuous basis functions

An adiabatic shear band (ASB) is a narrow region of intense plastic deformation that forms when some metallic alloys and some polymers are deformed at high strain rates and there is not enough time for the heat generated by plastic deformations to diffuse away. The study of ASBs is important because an ASB is a precursor to shear/ductile fractures. Initial-boundary-value problems simulating the initiation and propagation of an ASB have been analyzed usually using the finite element method (FEM). Because of the large plastic strains involved, the FE mesh needs to be refined several times to delineate the ASB width. Each refinement requires, in turn, interpolation of data from the previous mesh to the new one which causes a smoothening of the sharp gradients of the deformation fields, and affects characteristics of the ASB. In this paper, we propose the application of the finite element method with piecewise discontinuous basis functions for studying the occurrence of ASBs in simple shearing deformations of a body composed of an isotropic and homogeneous thermo-elastoviscoplastic material. The mathematical model of the problem is defined by a system of coupled nonlinear partial differential equations and an inequality constraint associated with the plastic strain rates admissibility.