Search results

When simulating fluid-structure interaction (FSI), it is often essential that the no-slip condition is accurately enforced at the wetted boundary of the structure. This paper aims to evaluate the relative strengths and limitations of the penalty and Lagrange multiplier methods, within the context of modelling FSI, through a comparative analysis.

In the immersed boundary method, the no-slip condition is typically imposed by augmenting the governing equations of the fluid with an artificial body force. The relative accuracy and computational time of the penalty and Lagrange multiplier formulations of this body force are evaluated by using each to solve three test problems, namely, flow through a channel, the harmonic motion of a cylinder through a stationary fluid and the vortex-induced vibration (VIV) of a cylinder.

The Lagrange multiplier formulation provided an accurate solution, especially when enforcing the no-slip condition, and was robust as it did not require “tuning” of problem specific parameters. However, these benefits came at a higher computational cost relative to the penalty formulation. The penalty formulation achieved similar levels of accuracy to the Lagrange multiplier formulation, but only if the appropriate penalty factor was selected, which was difficult to determine a priori.

Both the Lagrange multiplier and penalty formulations of the immersed boundary method are prominent in the literature. A systematic quantitative comparison of these two methods is presented within the same computational environment. A novel application of the Lagrange multiplier method to the modelling of VIV is also provided.

The purpose of this paper is to propose a combined reordering scheme with a wide range of application, called Reversed Cuthill-McKee-approximate minimum degree (RCM-AMD), to improve a preconditioned general minimal residual method for solving equations using Lagrange multiplier method, and facilitates the choice of the reordering for the iterative method.

To reordering the coefficient matrix before a preconditioned iterative method will greatly impact its convergence behavior, but the effect is very problem-dependent, even performs very differently when different preconditionings applied for an identical problem or the scale of the problem varies. The proposed reordering scheme is designed based on the features of two popular ordering schemes, RCM and AMD, and benefits from each of them.

Via numerical experiments for the cases of various scales and difficulties, the effects of RCM-AMD on the preconditioner and the convergence are investigated and the comparisons of RCM, AMD and RCM-AMD are presented. The results show that the proposed reordering scheme RCM-AMD is appropriate for large-scale and difficult problems and can be used more generally and conveniently. The reason of the reordering effects is further analyzed as well.

The proposed RCM-AMD reordering scheme preferable for solving equations using Lagrange multiplier method, especially considering that the large-scale and difficult problems are very common in practical application. This combined reordering scheme is more wide-ranging and facilitates the choice of the reordering for the iterative method, and the proposed iterative method has good performance for practical cases in in-house and commercial codes on PC.

This meshless collocation method is applicable not only to the Helmholtz equation with Dirichlet boundary condition but also mixed boundary conditions. It can calculate not only the high wavenumber problems, but also the variable wave number problems.

In this paper, the authors developed a meshless collocation method by using barycentric Lagrange interpolation basis function based on the Chebyshev nodes to deduce the scheme for solving the three-dimensional Helmholtz equation. First, the spatial variables and their partial derivatives are treated by interpolation basis functions, and the collocation method is established for solving second order differential equations. Then the differential matrix is employed to simplify the differential equations which is on a given test node. Finally, numerical experiments show the accuracy and effectiveness of the proposed method.

The numerical experiments show the advantages of the present method, such as less number of collocation nodes needed, shorter calculation time, higher precision, smaller error and higher efficiency. What is more, the numerical solutions agree well with the exact solutions.

Compared with finite element method, finite difference method and other traditional numerical methods based on grid solution, meshless method can reduce or eliminate the dependence on grid and make the numerical implementation more flexible.

The Helmholtz equation has a wide application background in many fields, such as physics, mechanics, engineering and so on.

This meshless method is first time applied for solving the 3D Helmholtz equation. What is more the present work not only gives the relationship of interpolation nodes but also the test nodes.

The purpose of this paper is to propose an efficient iterative method for large-scale finite element equations of bad numerical stability arising from deformation analysis with multi-point constraint using Lagrange multiplier method.

In this paper, taking warpage analysis of polymer injection molding based on surface model as an example, the performance of several popular Krylov subspace methods, including conjugate gradient, BiCGSTAB and generalized minimal residual (GMRES), with diffident Incomplete LU (ILU)-type preconditions is investigated and compared. For controlling memory usage, GMRES(m) is also considered. And the ordering technique, commonly used in the direct method, is introduced into the presented iterative method to improve the preconditioner.

It is found that the proposed preconditioned GMRES method is robust and effective for solving problems considered in this paper, and approximate minimum degree (AMD) ordering is most beneficial for the reduction of fill-ins in the ILU preconditioner and acceleration of the convergence, especially for relatively accurate ILU-type preconditioning. And because of concerns about memory usage, GMRES(m) is a good choice if necessary.

In this paper, for overcoming difficulties of bad numerical stability resulting from Lagrange multiplier method, together with increasing scale of problems in engineering applications and limited hardware conditions of computer, a stable and efficient preconditioned iterative method is proposed for practical purpose. Before the preconditioning, AMD reordering, commonly used in the direct method, is introduced to improve the preconditioner. The numerical experiments show the good performance of the proposed iterative method for practical cases, which is implemented in in-house and commercial codes on PC.

This article aims to present a direct solution to handle linear constraints in finite element (FE) analysis without penalties or the Lagrange multipliers introduced.

First, the system of linear equations corresponding to the linear constraints is solved for the leading variables in terms of the free variables and the constants. Then, the reduced system of equilibrium equations with respect to the free variables is derived from the finite-dimensional virtual work equation. Finally, the algorithm is designed.

The proposed procedure is promising in three typical cases: (1) to enforce displacement constraints in any direction; (2) to implement local refinements by allowing hanging nodes from element subdivision and (3) to treat non-matching grids of distinct parts of the problem domain. The procedure is general and suitable for 3D non-linear analyses.

The algorithm is fitted only to the Galerkin-based numerical methods.

The proposed procedure does not need Lagrange multipliers or penalties. The tangential stiffness matrix of the reduced system of equilibrium equations reserves positive definiteness and symmetry. Besides, many contemporary Galerkin-based numerical methods need to tackle the enforcement of the essential conditions, whose weak forms reduce to linear constraints. As a result, the proposed procedure is quite promising.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.

The purpose of this paper is to describe reduced order modelling based on dynamic flexibility approximation and applied to transient analyses.

This work is based on a recently proposed flexibility‐based component modes synthesis (CMS) approach which was shown to be very efficient for solving large eigenvalue problems. The model reduction approach is based on partionning via the localized Lagrange multipliers method, which makes it very appropriate to handle coupled problems.

In particular, it is demonstrated in this paper how the utilised model reduction method can be applied only to one part of the structure and efficiently coupled to a full finite element model. The performance of the method is investigated on numerical examples of plate and 3D problems.

The proposed flexibility‐based CMS approach can be used as a very efficient tool for complex engineering structures under dynamic load where the mode superposition method applies. The efficiency of the computations is brought about by the model reduction.

This article presents a geometrically non-linear finite element formulation for the analysis of planar two-layer beam-columns taking into account the inter-layer slip and uplift.

The co-rotational method is adopted, in which the motion of the element is decomposed into a rigid body motion and a small deformational one. The geometrically linear formulation can be used in the local frame and automatically be transformed into a geometrically nonlinear one. In co-rotational frame, both layers are assumed to be discretely connected at the element ends. Slips and uplifts are assumed to be small. Consequently, the condition of non interpenetration between the layers can be treated using a node-to-node contact algorithm. The resolution methods such as penalty (PM) and augmented Lagrangian method (ALM) with Uzawa updating scheme can be used.

The non-penetration condition between the layers of composite beams can be formulated by using contact law. It is found that despite a low convergence rate of augmented Lagrangian method compared to penalty method, the former prevents the unrealistic penetration. Besides, it is shown that the buckling load of the composite beam-column is largely affected by the uplift stiffness of the connectors.

The proposed finite element model is capable of simulating accurately the geometrically non-linear behavior of planar two-layer beam-columns taking into account the inter-layer slip and uplift. Regarding uplift, the non-penetration condition is strictly enforced by considering rigorous contact conditions at the interface. The constraint problem is solved using the penalty method or the augmented Lagrangian method with the Uzawa updating scheme.

The purpose of this paper is to present an application of the simulated annealing algorithm to the redundant system reliability optimization. Its main aim is to analyze and compare this optimization method performance with those of similar application.

The methods that were used to compare results are the genetic algorithm, the Lagrange Multipliers, and the evolution strategy. A hybrid algorithm composed by simulated annealing and genetic algorithm was developed in order to achieve the general applicability of the methods. The hybrid algorithm also tries to exploit the positive aspects of each method.

The results presented by the simulated annealing and the hybrid algorithm are significant, and validate the methods as a robust tool for parameter optimization in mechanical projects development.

The main objective is to propose a method for redundancy optimization in mechanical systems, which are not as large as electric and electronic systems, but involves high costs associated to redundancy and requires a high level of safety standards like: automotive and aerospace systems.

This paper aims to present a 3D static performance analysis model for the gas foil journal bearing to provide better understanding of the gas foil journal bearing and extend the development of the calculation about the static performance.

The foil bearing can be seen as a shell structure, and the mixed interpolation of tensorial components (MITC) element was used to build the shell model. The augmented Lagrange method was used to calculate the contact involving friction between foils and between the foil and the bearing sleeve. A displacement-controlled load scheme was used to calculate the deformation of the foils. A mapping operator was used to map the film pressure from the gas to the surface of the top foil.

This method provides high precision of calculation in the prediction of the static performance. The calculation results were compared with the experimental data, and they show good agreement. Meanwhile, the model can be applied in the prediction of the bearing performance in a broad range of working conditions.

This method extends the calculation of the gas foil journal bearing to a 3D scale and shows good agreement with the experimental data. Meanwhile, the present model has a good adaptability on the revolution speed and can be applied to the predictions in varied working conditions.