Search

The generalized Schwarz algorithm for a class of elliptic quasi-variational inequalities related to impulse control problems is studied in this paper. The principal result is to prove the error estimate in L ∞ -norm for m subdomains with overlapping nonmatching grids. This approach combines the geometrical convergence and the uniform convergence. [ABSTRACT FROM AUTHOR]

In this paper, a second-order hyperbolic equation is solved by a two-grid algorithm combined with the expanded mixed finite element method. The error estimate of the expanded mixed finite element method with discrete-time scheme is demonstrated. Moreover, we present a two-grid method and analyze its convergence. It is shown that the algorithm can achieve asymptotically optimal approximation as long as the mesh sizes satisfy H = O (h 1 2) . Finally, some numerical experiments are provided to illustrate the efficiency and accuracy of the proposed method. [ABSTRACT FROM AUTHOR]

This paper presents a new isogeometric mortar contact formulation based on an extended finite element interpolation to capture physical pressure discontinuities at the contact boundary. The so called two-half-pass algorithm is employed, which leads to an unbiased formulation and, when applied to the mortar setting, has the additional advantage that the mortar coupling term is no longer present in the contact forces. As a result, the computationally expensive segmentation at overlapping master-slave element boundaries, usually required in mortar methods (although often simplified with loss of accuracy), is not needed from the outset. For the numerical integration of general contact problems, the so-called refined boundary quadrature is employed, which is based on adaptive partitioning of contact elements along the contact boundary. The contact patch test shows that the proposed formulation passes the test without using either segmentation or refined boundary quadrature. Several numerical examples are presented to demonstrate the robustness and accuracy of the proposed formulation. [ABSTRACT FROM AUTHOR]

Thermal deformation is one of the principal factors that influence the machining accuracy of machine tools, and it can be improved by thermal error compensation. This paper presents a method of thermal error modeling using multisource information fusion, which can further improve the thermal error compensation accuracy. To set up a fusion model with optimal performance, two or more thermal error models should be established, from which a few models should be chosen to complement each other, and then combined into a synthesis model. In this paper, a dynamic thermal error model and a finite element model are combined to build a fusion model for lathe z-direction thermal error according to a fusion algorithm. The inputs for the fusion model are the values detected by the thermal sensors and the infrared imaging. An experiment carried out on a lathe verifies the validity of this modeling method. The results show that the multisource fusion model of thermal error can not only improve the prediction accuracy of thermal error over that of a single model, but also possesses better robustness. [ABSTRACT FROM AUTHOR]

The element energy projection (EEP) method for computation of superconvergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a result, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation of second-order as the model problem, this paper describes the related fundamental idea, the implementation strategy, and the computational algorithm. Representative numerical examples are given to show the efficiency, stability, versatility, and reliability of the proposed approach. [ABSTRACT FROM AUTHOR]

Mesh cleaning is the procedure of removing duplicate nodes, sequencing the indices of remaining nodes, and then updating the mesh connectivity for a topologically invalid Finite Element mesh. To the best of our knowledge, there has been no previously reported work specifically focused on the cleaning of large Finite Element meshes. In this paper we specifically present a generic and straightforward algorithm, MeshCleaner, for cleaning large Finite Element meshes. The presented mesh cleaning algorithm is composed of (1) the stage of compacting and reordering nodes and (2) the stage of updating mesh topology. The basic ideas for performing the above two stages efficiently both in sequential and in parallel are introduced. Furthermore, one serial and two parallel implementations of the algorithm MeshCleaner are developed on multi-core CPU and/or many-core GPU. To evaluate the performance of our algorithm, three groups of experimental tests are conducted. Experimental results indicate that the algorithm MeshCleaner is capable of cleaning large meshes very efficiently, both in sequential and in parallel. The presented mesh cleaning algorithm MeshCleaner is generic, simple, and practical. [ABSTRACT FROM AUTHOR]

Shot peen forming is often used to form the aerodynamic configurations of aircraft integral wing panels. The component is traditionally regarded as perfectly constrained target in peen forming process. However, weak-stiffness of the component results in concomitant macroscopic bending and stretching deformations, and releasing of peening stresses. Subsequent peening stresses are induced by shot impacts based on the current released stresses. In this paper, a intermittently peening-releasing finite element simulation model is used to study peening stress evolution in lowrigidity plate. An algorithm is proposed to calculate the increment of subsequent peening stress by connecting current stresses to the peening stresses in perfectly constrained target.With the algorithm, the peen forming deformations of freely constrained plates are calculated by equivalently simulating the peen forming process. And strip peen forming deformations are calculated and compared with experimental data. The comparison shows the peening deformation of a low-rigidity component should be treated as an incremental process and the algorithm can be used to predict the deformation with good accuracy. [ABSTRACT FROM AUTHOR]

A review of the recent development of the machining fixture configuration/layout is presented in this paper. The literature review is mainly focused on the recently developed optimal fixture configuration under the dynamic conditions of the workpiece. In this review paper, the fixture design, fixture analysis, fixture synthesis, fixture layout design, optimization of fixture layout design, various optimization algorithms, and case studies of two- and three-dimensional workpiece geometries under dynamic conditions have been emphasized specially. The further scopes of the research are finally summarized. [ABSTRACT FROM AUTHOR]

Rapid prototyping technology makes the manufacturing of functionally graded material (FGM) objects possible. The FGM production process requires concurrent applications of digital design and manufacturing methods. For this reason, it is imperative to construct and capture information and other critical data regarding the geometry and materials about FGM objects for use in the manufacturing process. This paper proposes a simplex-clear data format for modeling FGM objects. This proposed data format captures both the geometry and material information of an FGM object. It also presents a finite element-based slicing algorithm which slices the FGM object into layers and captures each layer's information for rapid prototyping manufacturing of the object. An example is given at the end of the paper to validate the data format and demonstrate the adaptive slicing algorithm. [ABSTRACT FROM AUTHOR]

In this paper, we construct a second order algorithm based on the spectral deferred correction method for the time-dependent magnetohydrodynamics flows at a low magnetic Reynolds number. We present a complete theoretical analysis to prove that this algorithm is unconditionally stable, consistent and second order accuracy. Finally, two numerical examples are given to illustrate the convergence and effectiveness of our algorithm. [ABSTRACT FROM AUTHOR]

In the present paper, a recent deterministic continuum-strategy two-player discrete-time dynamic leader-follower game with fixed finite time duration and closed-loop information structure is studied. The types of the considered payoff functions can be widely used in different applications (mainly in conflicts of consuming a limited resource, where one player, called the leader, is a superior authority choosing a strategy choice first, and another player, called the follower, chooses after). In case of certain payoff convexity, explicit conditions are given, when it can be known in advance that an equilibrium exists and consists of only two possible choices of both players at each step. The sub-game equilibrium from a given step may depend on the former selections of the players. Thus the continuum-strategy problem has been reduced to a general finite game of two possible choices corresponding to both players. Such type of games could be solved in a standard way with dynamic programming using a computer. Nevertheless, the game can be further simplified, and then an equilibrium can be directly determined, such decreasing the computational demand to a great extent. A solution algorithm and practical examples are also given to support the real-life application of the results. [ABSTRACT FROM AUTHOR]

In this paper, a new type of hybrid finite element method (FEM), hybrid fundamental-solution-based FEM (HFS-FEM), is developed for analyzing plane piezoelectric problems by employing fundamental solutions (Green's functions) as internal interpolation functions. A modified variational functional used in the proposed model is first constructed, and then the assumed intra-element displacement fields satisfying a priori the governing equations of the problem are constructed by using a linear combination of fundamental solutions at a number of source points located outside the element domain. To ensure continuity of fields over inter-element boundaries, conventional shape functions are employed to construct the independent element frame displacement fields defined over the element boundary. The proposed methodology is assessed by several examples with different boundary conditions and is also used to investigate the phenomenon of stress concentration in infinite piezoelectric medium containing a hole under remote loading. The numerical results show that the proposed algorithm has good performance in numerical accuracy and mesh distortion insensitivity compared with analytical solutions and those from ABAQUS. In addition, some new insights on the stress concentration have been clarified and presented in the paper. [ABSTRACT FROM AUTHOR]

Many applications in computational science and engineering concern composite materials, which are characterized by large discontinuities in the material properties. Such applications require fine-scale finite-element meshes, which lead to large linear systems that are challenging to solve with current direct and iterative solutions algorithms. In this paper, we consider the simulation of asphalt concrete, which is a mixture of components with large differences in material stiffness. The discontinuities in material stiffness give rise to many small eigenvalues that negatively affect the convergence of iterative solution algorithms such as the preconditioned conjugate gradient (PCG) method. This paper considers the deflated preconditioned conjugate gradient (DPCG) method in which the rigid body modes of sets of elements with homogeneous material properties are used as deflation vectors. As preconditioner we consider several variants of the algebraic multigrid smoothed aggregation method. We evaluate the performance of the DPCG method on a parallel computer using up to 64 processors. Our test problems are derived from real asphalt core samples, obtained using CT scans. We show that the DPCG method is an efficient and robust technique for solving these challenging linear systems. [ABSTRACT FROM AUTHOR]

The immersed finite element method (IFEM) is a novel numerical approach to solve fluid-structure interaction types of problems that utilizes non-conforming meshing concept. The fluid and the solid domains are represented independently. The original algorithm of the IFEM follows the interpolation process as illustrated in the original immersed boundary method where the fluid velocity and the interaction force are explicitly coupled. However, the original approach presents many numerical difficulties when the fluid and solid physical properties have large mismatches, such as when the density difference is large and when the solid is a very stiff material. Both situations will lead to divergent or unstable solutions if not handled properly. In this paper, we develop a semi-implicit formulation of the IFEM algorithm so that several terms of the interfacial forces are implicitly evaluated without going through the force distribution process. Based on the 2-D and 3-D examples that we study in this paper, we show that the semi-implicit approach is robust and is capable of handling these highly discontinuous physical properties quite well without any numerical difficulties. [ABSTRACT FROM AUTHOR]

A sub-cycling integration algorithm (or named multi-time-steps integration algorithm), which has been successfully applied to FEM dynamical analysis, was firstly presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18:259–275, 1979). However, the problem of how to apply this type of algorithm to flexible multi-body dynamics (FMD) problems still lacks investigation up to now. Similar to the region-partitioning method used in FEM, this paper presents a central-difference-based sub-cycling integral method by decomposing the variables of an FMD equation into several groups and adopting different integral step sizes to each group of the variables. Based on the condensed form of an FMD equation, a group of common update formulae and a sub-step update formula, which constitute the sub-cycling together, are established in the paper. Furthermore, an implementation flowchart of the sub-cycling is presented. Stability of the sub-cycling will be analyzed and numerical examples will be performed to verify availability and precision of the sub-cycling in part II of the paper. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to detect damage (stiffness degradation) of laminated composite plates from noisy impact response data. The combined finite element method (FEM) with five degrees of freedom (DOF) and the advanced noise filtering algorithm described in this paper may allow us not only to detect the deteriorated elements but also to find their locations and the extents. A first order shear deformation theory (FSDT) is used to predict the structural behavior and to detect damage of laminated composite plates. The filtering procedure is designed by means of a wavelet decomposition together with a selection of the measuring points, and the optimization criterion is constructed on an estimate of the probability of detection using genetic algorithms. All these techniques are applied for the first time to composites. The effects of filtered noise associated with the uncertainty of measurements due to the complex nature of composites are considered for different layup sequences, number of layers, and length–thickness ratios. Several numerical results show that the noise filtering system is computationally efficient in identifying stiffness degradation for complex structures such as laminated composites. [ABSTRACT FROM AUTHOR]

In this paper a three-dimensional Voronoi cell finite element model is developed for analyzing heterogeneous materials containing a dispersion of ellipsoidal inclusions or voids in the matrix. The paper starts with a description of 3D tessellation of a domain with ellipsoidal heterogeneities, to yield a 3D mesh of Voronoi cells containing the heterogeneities. A surface based tessellation algorithm is developed to account for the shape and size of the ellipsoids in point based tessellation methods. The 3D Voronoi cell finite element model, using the assumed stress hybrid formulation, is developed for determining stresses and displacements in a linear elastic material domain. Special stress functions that introduce classical Lamé functions in ellipsoidal coordinates are implemented to enhance solution convergence. Numerical methods for implementation of algorithms and yielding stable solutions are discussed. Numerical examples are conducted with inclusions and voids to demonstrate the effectiveness of the model. [ABSTRACT FROM AUTHOR]

In order to generate efficient tool path with given precision requirements, scallop height should be kept under a given limit, while the tool path should be as short as possible to reduce machining time. Traditional methods generate CC curves one by one, which makes the final tool path far from being globally optimal. This paper presents an optimal tool path generation model for a ball-end tool which strives to globally optimize a tool path with various objectives and constraints. Two scalar functions are constructed over the part surface to represent the path intervals and the feedrate (with directions). Using the finite element method (FEM), the tool path length minimization model and the machining time minimization model are solved numerically. The proposed method is also suitable for tool path generation on mesh surfaces. Simulation results show that the generated tool path can be direction parallel or contour parallel with different boundary conditions. Compared to most of the conventional tool path generation methods, the proposed method is able to generate more effective tool paths due to the global optimization strategy. [ABSTRACT FROM AUTHOR]

This paper presents a method to find all intersections between curved lines such as structural line elements and finite element meshes with intentions to generate smaller, non-compatible, line cells (e.g. bar elements) between crossings. The intersection finding algorithm works for two and three-dimensional meshes constituted by linear, quadratic or higher order elements. Using the proposed algorithm, meshes can then be automatically prepared for finite element analyses with techniques for embedding elements within others or analyses that require lines within solids. The application of the method is demonstrated by a number of numerical examples illustrating its capabilities in handling complex geometries, relative speed and convenience. [ABSTRACT FROM AUTHOR]

In this paper, we study a two level method based on Newton's iteration for the nonlinear system arising from the Galerkin finite element approximation to the equations of motion described by the Kelvin-Voigt viscoelastic fluid flow model. The two-grid algorithm is based on three steps. In the first step, the nonlinear system is solved on a coarse mesh $\mathcal {T}_{H}$ to obtain an approximate solution u. In the second step, the nonlinear system is linearized around u based on Newton's iteration and the linear system is solved on a finer mesh $\mathcal {T}_{h}$. Finally, in the third step, a correction to the results obtained in the second step is achieved by solving a linear problem with a different right hand side on $\mathcal {T}_{h}$. Optimal error estimates in L( L)-norm, when $h=\mathcal {O} (H^{2-\delta })$ and in L()-norm, when $h=\mathcal {O}(H^{5-2\delta })$ for velocity and in L( L)-norm, when $h=\mathcal {O}(H^{5-2\delta })$ for pressure are established, where δ > 0 is arbitrarily small for two dimensions and $\delta =\frac {1}{2}$ for three dimensions. [ABSTRACT FROM AUTHOR]

In this paper, a coupled hydro-mechanical formulation is developed for deformable porous media subjected to crack interfaces in the framework of extended finite element method. Governing equations of the porous medium consist of the momentum balance of the bulk together with the momentum balance and continuity equations of the fluid phase, known as [InlineEquation not available: see fulltext.] formulation. The discontinuity in fractured porous medium is modeled for both opening and closing modes that results in the fluid flow within the fracture, and/or contact behavior at the crack edges. The fluid flow through the fracture is assumed to be viscous and is modeled by employing the Darcy law in which the permeability of fracture is obtained using the cubic law. The contact condition in fractured porous medium is handled by taking the advantage from two different algorithms of LATIN method and penalty algorithm. The effect of contact on fluid phase is employed by considering no leak-off from/into the porous medium. The nonlinearity of coupled equations produced due to opening and closing modes is carried out using an iterative algorithm in the Newton-Raphson procedure. Finally, several numerical examples are solved to illustrate the performance of proposed X-FEM method for hydro-mechanical behavior of fractured porous media with opening and closing modes. [ABSTRACT FROM AUTHOR]

In the current paper we present a fast, reliable technique for simulating wave propagation in complex structures made of heterogeneous materials. The proposed approach, the spectral cell method, is a combination of the finite cell method and the spectral element method that significantly lowers preprocessing and computational expenditure. The spectral cell method takes advantage of explicit time-integration schemes coupled with a diagonal mass matrix to reduce the time spent on solving the equation system. By employing a fictitious domain approach, this method also helps to eliminate some of the difficulties associated with mesh generation. Besides introducing a proper, specific mass lumping technique, we also study the performance of the low-order and high-order versions of this approach based on several numerical examples. Our results show that the high-order version of the spectral cell method together requires less memory storage and less CPU time than other possible versions, when combined simultaneously with explicit time-integration algorithms. Moreover, as the implementation of the proposed method in available finite element programs is straightforward, these properties turn the method into a viable tool for practical applications such as structural health monitoring [-], quantitative ultrasound applications [], or the active control of vibrations and noise [, ]. [ABSTRACT FROM AUTHOR]

The present paper focuses on the development of a fully implicit, incrementally objective integration algorithm for a hypoelastic formulation of $$J_{2}$$ -viscoplasticity, which employs the mechanical threshold strength model to compute the material's flow stress, taking into account its dependence on strain rate and temperature. Heat generation due to high-rate viscoplastic deformation is accounted for, assuming adiabatic conditions. The implementation of the algorithm is discussed, and its performance is assessed in the contexts of implicit and explicit dynamic finite element analysis, with the aid of example problems involving a wide range of loading rates. Computational results are compared to experimental data, showing very good agreement. [ABSTRACT FROM AUTHOR]

Compressed sensing ensures the accurate reconstruction of sparse signals from far fewer samples than required in the classical Shannon-Nyquist theorem. In this paper, a generalized hard thresholding pursuit (GHTP) algorithm is presented that can recover unknown vectors without the sparsity level information. We also analyze the convergence of the proposed algorithm. Numerical experiments are given for synthetic and real-world data to illustrate the validity and the good performance of the proposed algorithm. [ABSTRACT FROM AUTHOR]

Finite element based formulations for flexible multibody systems are becoming increasingly popular and as the complexity of the configurations to be treated increases, so does the computational cost. It seems natural to investigate the applicability of parallel processing to this type of problems; domain decomposition techniques have been used extensively for this purpose. In this approach, the computational domain is divided into non-overlapping sub-domains, and the continuity of the displacement field across sub-domain boundaries is enforced via the Lagrange multiplier technique. In the finite element literature, this approach is presented as a mathematical algorithm that enables parallel processing. In this paper, the divided system is viewed as a flexible multibody system, and the sub-domains are connected by kinematic constraints. Consequently, all the techniques applicable to the enforcement of constraints in multibody systems become applicable to the present problem. In particular, it is shown that a combination of the localized Lagrange multiplier technique with the augmented Lagrange formulation leads to interesting solution strategies. The proposed algorithm is compared with the well-known FETI approach with regards to convergence and efficiency characteristics. The present algorithm is relatively simple and leads to improved convergence and efficiency characteristics. Finally, implementation on a parallel computer was conducted for the proposed approach. [ABSTRACT FROM AUTHOR]

In this paper, heat transfer problems with sharp spatial gradients are analyzed using the Generalized Finite Element Method with global-local enrichment functions ( GFEM). With this approach, scale-bridging enrichment functions are generated on the fly, providing specially-tailored enrichment functions for the problem to be analyzed with no a-priori knowledge of the exact solution. In this work, a decomposition of the linear system of equations is formulated for both steady-state and transient heat transfer problems, allowing for a much more computationally efficient analysis of the problems of interest. With this algorithm, only a small portion of the global system of equations, i.e., the hierarchically added enrichments, need to be re-computed for each loading configuration or time-step. Numerical studies confirm that the condensation scheme does not impact the solution quality, while allowing for more computationally efficient simulations when large problems are considered. We also extend the GFEM to allow for the use of hexahedral elements in the global domain, while still using tetrahedral elements in the local domain, to allow for automatic localized mesh refinement without the use of constrained approximations. Simulations are run with the use of linear and quadratic hexahedral and tetrahedral elements in the global domain. Convergence studies indicate that the use of a different partition of unity (PoU) in the global (hexahedral elements) and local (tetrahedral elements) domains does not adversely impact the solution quality. [ABSTRACT FROM AUTHOR]

The rotational quaternions represent a unique four dimensional parametrization of rotations in the three dimensional Euclidean space. In the present paper they are used as the basic rotational parameters in formulating the finite-element approach of geometrically exact beam-like structures. The classical concept of parameterizing the rotation matrix by the rotational vector is completely abandoned so that the only rotational parameters are the rotational quaternions representing both rotations and rotational strains in the beam. The space discretization based on the collocation method is used and the adjustment of the Newmark time-integration algorithm to the quaternion parameterizations of rotation is presented. [ABSTRACT FROM AUTHOR]

In this paper we present a new approach for solving elastoplastic problems as second order cone complementarity problems (SOCCPs). Specially, two classes of elastoplastic problems, i.e. the J plasticity problems with combined linear kinematic and isotropic hardening laws and the Drucker-Prager plasticity problems with associative or non-associative flow rules, are taken as the examples to illustrate the main idea of our new approach. In the new approach, firstly, the classical elastoplastic constitutive equations are equivalently reformulated as second order cone complementarity conditions. Secondly, by employing the finite element method and treating the nodal displacements and the plasticity multiplier vectors of Gaussian integration points as the unknown variables, we obtain a standard SOCCP formulation for the elastoplasticity analysis, which enables the using of general SOCCP solvers developed in the field of mathematical programming be directly available in the field of computational plasticity. Finally, a semi-smooth Newton algorithm is suggested to solve the obtained SOCCPs. Numerical results of several classical plasticity benchmark problems confirm the effectiveness and robustness of the SOCCP approach. [ABSTRACT FROM AUTHOR]

This paper studies adaptive finite element methods (AFEMs), based on piecewise linear elements and newest vertex bisection, for solving second order elliptic equations with piecewise constant coefficients on a polygonal domain Ω⊂ℝ. The main contribution is to build algorithms that hold for a general right-hand side f∈ H( Ω). Prior work assumes almost exclusively that f∈ L( Ω). New data indicators based on local H norms are introduced, and then the AFEMs are based on a standard bulk chasing strategy (or Dörfler marking) combined with a procedure that adapts the mesh to reduce these new indicators. An analysis of our AFEM is given which establishes a contraction property and optimal convergence rates N with 0< s≤1/2. In contrast to previous work, it is shown that it is not necessary to assume a compatible decay s<1/2 of the data estimator, but rather that this is automatically guaranteed by the approximability assumptions on the solution by adaptive meshes, without further assumptions on f; the borderline case s=1/2 yields an additional factor log N. Computable surrogates for the data indicators are introduced and shown to also yield optimal convergence rates N with s≤1/2. [ABSTRACT FROM AUTHOR]

Analysis of transverse vibration of the gear found in a high-speed gearbox considered as an annular plate reflecting gear geometry is the subject of this paper. How gear angular velocity affects the deformation of normal modes of transverse vibration of the system under consideration is analysed. Models considered were discretized by the finite elements method. Numerical computations have been performed in the ANSYS environment. The algorithm to identify the proper distorted mode shapes is presented. The Campbell diagram for the system under consideration is elaborated. The problems discussed here can be useful for engineers dealing with dynamics of rotating machine systems. [ABSTRACT FROM AUTHOR]

In this paper, the constitutive relations of the finite element method are constructed and used for solving two- and three-dimensional problems of the asymmetric theory of elasticity. Different variants of finite elements are considered. The numerical experiments are carried out to evaluate the reliability and computational efficiency of the finite element algorithm based on the comparison between the numerical and analytical solutions, numerical estimation of the convergence and checking of the degree of accuracy, to which the natural boundary conditions are satisfied. The obtained solutions to the two- and three-dimensional problems are interpreted from the viewpoint of their applicability to a design of experiments capable of revealing the facts of couple-stress effects in material under elastic deformation and identification of material constants for the asymmetric theory of elasticity. The capabilities of the finite element algorithm to interpret experimental data and estimate the errors occurring in real experiments have been tested by solving several example problems. [ABSTRACT FROM AUTHOR]

In this paper we propose a finite element approach which simulates the mechanical behaviour of beam assemblies that are subject to large deformations and that develop contact-friction interactions. We focus on detecting and modeling contact-friction interactions within the assembly of beams. Contact between beams-or between parts of the same beam in the case of self-contact, is detected from intermediate geometries defined within proximity zones associating close parts of beam axes. The discretization of contact-friction interactions is performed on these intermediate geometries by means of contact elements, constituted of pairs of material particles which are predicted to enter into contact. A 3D finite strain beam model is used to represent the behaviour of each beam. This model describes the kinematics of each beam cross-section using nine degrees of freedom, and is therefore able to represent plane deformations of these cross-sections. Algorithms are proposed to solve the global nonlinear problem using an implicit scheme, under quasi-static assumptions. Simulation results of the tightening and releasing of knots made on monofilament and multifilament yarns are shown as an application. Straight fibers are first twisted together to make a yarn, before suitable conditions are applied to their ends to form and tighten the knot. Tightening forces are finally released to obtain an equilibrium configuration of the knot without external forces. [ABSTRACT FROM AUTHOR]

A contact algorithm in the context of the combined discrete element (DE) and finite element (FE) method is proposed. The algorithm, which is based on the node-to-surface method used in finite element method, treats each spherical discrete element as a slave node and the surfaces of the finite element domain as the master surfaces. The contact force on the contact interface is processed by using a penalty function method. Afterward, a modification of the combined DE/FE method is proposed. Following that, the corresponding numerical code is implemented into the in-house developed code. To test the accuracy of the proposed algorithm, the impact between two identical bars and the vibration process of a laminated glass plate under impact of elastic sphere are simulated in elastic range. By comparing the results with the analytical solution and/or that calculated by using LS-DYNA, it is found that they agree with each other very well. The accuracy of the algorithm proposed in this paper is proved. [ABSTRACT FROM AUTHOR]

This paper focuses on studying and adapting modeling techniques using the finite element method to simulate the rigid die compaction of metal powders. First, it presents the implementation of the cap constitutive model into ABAQUS FE software using the closest point projection algorithm. Then, an inverse modeling procedure was proposed to alleviate the problems raised by the interpretation of the experimental tests and to more accurately determine the material parameters. The objective function is formed, based on the discrepancy in density data between the numerical model prediction and the experiment. Minimization of the objective function with respect to the material parameters was performed using an in-house optimization software shell built on a modified Levenberg-Marquardt method. Thus, an integrated simulation module consisting of an inverse optimization method and a finite element method was developed for modeling the powder compaction process as a whole. The simulation and identification module developed was applied to simulate the compaction of some industrial parts. The results reveal that the maximum absolute error between densities is 2.3%. It corresponds to the precision of the experimental method. [ABSTRACT FROM AUTHOR]

In this paper, an implicit coupling algorithm for fluid-structure interaction problems with under-time steps for the solid is presented. Its implementation on two configurations is achieved by using the CASTEM finite-elements code. First, the free oscillations of a cylinder in an annular fluid domain where its movement is determined by the coupled fluid-solid action is considered in the case of viscous fluid. It should be noted that the implicit coupling algorithm gives the best prediction of the structure oscillations. The under-time steps for the solid are introduced in order to obtain better results. Then, an application whose final objective is to model a floating barrage is studied. The main goal of this application is to predict the displacements of a ring completely immersed and anchored by a cable to the lower boundary of the fluid domain. The finite-element discretization of the Navier-Stokes equations in the ALE formulation is used [ABSTRACT FROM AUTHOR]

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, ), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples. [ABSTRACT FROM AUTHOR]

Abrasive waterjet machining (AWJM) is a non-conventional process. The mechanism of material removing in AWJM for ductile materials and existing erosion models are reviewed in this paper. To overcome the difficulties of fluid–solid interaction and extra-large deformation problem using finite element method (FEM), the SPH-coupled FEM modeling for abrasive waterjet machining simulation is presented, in which the abrasive waterjet is modeled by SPH particles and the target material is modeled by FE. The two parts interact through contact algorithm. The creativity of this model is multi-materials SPH particles, which contain abrasive and water and mix together uniformly. To build the model, a randomized algorithm is proposed. The material model for the abrasive is first presented. Utilizing this model, abrasive waterjet penetrating the target materials with high velocity is simulated and the mechanism of erosion is depicted. The relationship between the depth of penetration and jet parameters, including water pressure and traverse speed, etc., are analyzed based on the simulation. The results agree with the experimental data well. It will be a benefit to understand the abrasive waterjet cutting mechanism and optimize the operating parameters. [ABSTRACT FROM AUTHOR]

An algorithm combining three-dimensional (3D) discrete and finite element methods is proposed. This new approach is conducted by decomposing the calculation domain into a finite element (FE) calculation domain and a discrete element (DE) calculation domain; the interaction between the two sub-domains is processed by using a penalty function method. Following the established model that combines spherical DEs and FEs, the corresponding numerical code is developed. The vibration process of two cantilever beams under dynamic force is simulated. By comparing the results calculated with different penalty factors set and also with that calculated by the finite element code LS-DYNA, it is found that the calculated results are unanimous and the precision is almost the same as LS-DYNA, as long as the penalty factor is large enough. Moreover, the vibration processes of two plates under impact of rigid spheres are simulated and the accuracy of the model proposed in this paper is further proved in the field of contact mechanics by comparing the simulating results with that calculated by using LS-DYNA. Finally, the impact fracture behavior of a laminated glass plate is simulated, with the influence of model parameters taken into consideration. And the numerical experiments show that the combined model can be used to predict some macroscopical physical quantities, such as the impact force of impactor. [ABSTRACT FROM AUTHOR]

For certain continuum problems, it is desirable and beneficial to combine two different methods together in order to exploit their advantages while evading their disadvantages. In this paper, a bridging transition algorithm is developed for the combination of the meshfree method (MM) with the finite element method (FEM). In this coupled method, the MM is used in the sub-domain where the MM is required to obtain high accuracy, and the FEM is employed in other sub-domains where FEM is required to improve the computational efficiency. The MM domain and the FEM domain are connected by a transition (bridging) region. A modified variational formulation and the Lagrange multiplier method are used to ensure the compatibility of displacements and their gradients. To improve the computational efficiency and reduce the meshing cost in the transition region, regularly distributed transition particles, which are independent of either the meshfree nodes or the FE nodes, can be inserted into the transition region. The newly developed coupled method is applied to the stress analysis of 2D solids and structures in order to investigate its’ performance and study parameters. Numerical results show that the present coupled method is convergent, accurate and stable. The coupled method has a promising potential for practical applications, because it can take advantages of both the MM and FEM when overcome their shortcomings. [ABSTRACT FROM AUTHOR]

Fluid–structure interaction computations based on interface-tracking (moving-mesh) techniques are often hindered if the structural surfaces come in contact with each other. As the distance between two structural surfaces tends to zero, the fluid mesh in between distorts severely and eventually becomes invalid. Our objective is to develop a technique for modeling problems where the contacting structural surfaces would otherwise inhibit flow modeling or even fluid-mesh update. In this paper, we present our contact tracking technique that detects impending contact and maintains a minimum distance between the contacting structural surfaces. Our Surface-Edge-Node Contact Tracking (SENCT) technique conducts a topologically hierarchical search to detect contact between each node and the elements (“surfaces”), edges and other nodes. To keep the contacting surfaces apart by a small distance, we apply to the contacted nodes penalty forces in SENCT-Force (SENCT-F) and displacement restrictions in SENCT-Displacement (SENCT-D). By keeping a minimum distance between the contacting surfaces, we are able to update the fluid mesh in between and model the flow accurately. [ABSTRACT FROM AUTHOR]

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations. [ABSTRACT FROM AUTHOR]

Imposing rigidity constraints of an immersed elastic body in a transient flow field is not trivial. It requires solution stability and accuracy. In this paper, we present an efficient and accurate algorithm implemented to enforce fluid–structure interface constraints used in the immersed finite element method (IFEM). This interface treatment is a constraint applied onto the rigid bodies based on the fluid structure interaction force evaluated from the immersed solid object. It requires no ad hoc constants or adjustments, thus providing numerical stability and avoiding unnecessary trial-and-error procedures in defining the stiffness of the elastic body. This force term can be evaluated for both uniform and nonuniform fluid grids based on the higher order interpolation function adopted in the IFEM. The ability in handling nonuniform interpolations offers the convenience in modeling arbitrary geometrical shapes and provides solution refinements around interfaces. The results we obtained from flow past a rigid cylinder demonstrate that this convenient way of constraining the interface is a reliable and robust numerical approach to solve unsteady fluid flow interacting with immersed rigid bodies. [ABSTRACT FROM AUTHOR]

Abstract  The first extension of the virtual fields method to the reconstruction of heterogeneous stiffness properties from 3D bulk full-field displacement data is presented in this paper. Data are provided by Magnetic Resonance Imaging (MRI). Two main issues are addressed: 1. the identification of the stiffness ratio between two different media in a heterogeneous solid; 2. the reconstruction of stiffness heterogeneities buried in a heterogeneous solid. The approach is based on a finite element discretization of the equilibrium equations. It is tested on experimental full-field data obtained on a phantom with the stimulated echo MRI technique. The phantom is made of a stiff spherical inclusion buried within a lower modulus material. Preliminary independent tests showed that the material of the inclusion was four times stiffer than the surrounding material. This ratio value is correctly identified by our approach directly on the phantom with the MRI data. Moreover, the modulus distribution is promisingly reconstructed across the whole investigated volume. However, the resulting modulus distribution is highly variable. This is explained by the fact that the approach relies on a second order differentiation of the data, which tends to amplify noise. Noise is significantly reduced by using appropriate filtering algorithms. [ABSTRACT FROM AUTHOR]

Abstract  This paper deals with the identification of elasto-plastic constitutive parameters from deformation fields measured over the surface of thin flat specimens with the grid method. The approach for recovering the constitutive parameters is the virtual fields method. A dedicated algorithm is used for deriving the distribution of the 2D stress components from the measured deformation fields. A state of plane stress is assumed. Guesses of the constitutive parameters are input in the algorithm and updated until the stresses satisfy the principle of virtual work in the least squares sense. The advantage of this approach is that it can handle very heterogeneous plastic flows and it is much faster than classical finite element model updating approaches. An experimental application is provided to demonstrate it. Six mild steel double-notched specimens have been tested in a configuration combining tension and in-plane bending. The identified parameters are in good agreement with their reference counterparts. Stress fields are eventually reconstructed across the specimen all along the test for analyzing the evolution of the plastic flow. [ABSTRACT FROM AUTHOR]

Abstract  In this paper, adaptive finite element method is developed for the estimation of distributed parameter in elliptic equation. Both upper and lower error bound are derived and used to improve the accuracy by appropriate mesh refinement. An efficient preconditioned project gradient algorithm is employed to solve the nonlinear least-squares problem arising in the context of parameter identification problem. The efficiency of our error estimators is demonstrated by some numerical experiments. [ABSTRACT FROM AUTHOR]

In this paper a general formulation for finding the maximum allowable dynamic load (MADL) of flexible link mobile manipulators is presented. The main constraints used for the proposed algorithm are the actuator torque capacity and the limited error bound for the end-effector during motion on a given trajectory. The accuracy constraint is taken into account with two boundary lines which are equally offset due to the given end-effector trajectory, while a speed-torque characteristics curve of a typical DC motor is used for applying the actuator torque constraint. The finite element method (FEM), which is able to consider the full nonlinear dynamics of mobile manipulators, is applied to derive the kinematic and dynamic equations. In order to verify the effectiveness of the presented algorithm, two simulation studies considering a flexible two-link planar manipulator mounted on a mobile base are presented and the results are discussed. [ABSTRACT FROM AUTHOR]

In this paper a general formula for finding the maximum allowable dynamic load (MADL) of flexible link mobile manipulators is presented. The main constraints used for the proposed algorithm are the actuator torque capacity and the limited error bound for the end-effector during motion on a given trajectory. The accuracy constraint is taken into account with two boundary lines which are equally offset due to the given end-effector trajectory, while a speed-torque characteristics curve of a typical DC motor, is used for applying the actuator torque constraint. Finite element method (FEM), which is able to consider the full nonlinear dynamic of mobile manipulator is applied to derive the kinematic and dynamic equations. In order to verify the effectiveness of the presented algorithm, two simulation studies considering a flexible two-link planar manipulator mounted on a mobile base are presented and the results are discussed. [ABSTRACT FROM AUTHOR]

This paper presents our new development of parallel finite element algorithms for elastic–plastic problems. The proposed method is based on dividing the original structure under consideration into a number of substructures which are treated as isolated finite element models via the interface conditions. Throughout the analysis, each processor stores only the information relevant to its substructure and generates the local stiffness matrix. A parallel substructure oriented preconditioned conjugate gradient method, which is combined with MR smoothing and diagonal storage scheme are employed to solve linear systems of equations. After having obtained the displacements of the problem under consideration, a substepping scheme is used to integrate elastic–plastic stress–strain relations. The procedure outlined controls the error of the computed stress by choosing each substep size automatically according to a prescribed tolerance. The combination of these algorithms shows a good speedup when increasing the number of processors and the effective solution of 3D elastic–plastic problems whose size is much too large for a single workstation becomes possible. [ABSTRACT FROM AUTHOR]

This paper presents an optimization methodology for non-rigid sheet metal assembly variation by considering part variation, fixture variation, fixture layout, and joint positions, as well as the assembly spring back. The proposed algorithm integrates the finite element analysis (FEA) with a powerful global optimization algorithm, called the mode-pursuing sampling (MPS) method to simultaneously search for the optimal fixture and joint positions in order to minimize the assembly variation. An example application study is presented to demonstrate the optimization procedure and its effectiveness. [ABSTRACT FROM AUTHOR]

This paper presents a new contact searching algorithm for large deformation mortar-based contact formulations. In this algorithm, a bounding volume hierarchy, defined in the context of a binary tree, is built for each contact surface based on the geometry of the surface. A global contact searching procedure based on these bounding volume trees is first performed to find all candidate contact element pairs, and then a local searching procedure is done to find all the mortar segments having contributions to the mortar integrals that define the contact formulation. The searching algorithm is shown to be very efficient and readily applicable to a variety of large sliding contact problems. [ABSTRACT FROM AUTHOR]