Search results

This paper aims to develop a new isogeometric boundary element formulation based on non-uniform rational basis splines (NURBS) curves for solving Reissner’s shear-deformable plates.

The generalized displacements and tractions along the problem boundary are approximated as NURBS curves having the same rational B-spline basis functions used to describe the geometrical boundary of the problem. The source points positions are determined over the problem boundary by the well-known Greville abscissae definition. The singular integrals are accurately evaluated using the singularity subtraction technique.

Numerical examples are solved to demonstrate the validity and the accuracy of the developed formulation.

This formulation is considered to preserve the exact geometry of the problem and to reduce or cancel mesh generation time by using NURBS curves employed in computer aided designs as a tool for isogeometric analysis. The present formulation extends such curves to be implemented as a stress analysis tool.

The purpose of this paper is to further simplify the use of NURBS in industrial environnements. Although isogeometric analysis (IGA) has been the object of intensive studies over the past decade, its massive deployment in industrial analysis still appears quite marginal. This is partly due to its implementation, which is not straightforward with respect to the elementary structure of finite element (FE) codes. This often discourages industrial engineers from adopting isogeometric capabilities in their well-established simulation environment.

Based on the concept of Bézier and Lagrange extractions, a novel method is proposed to implement IGA from an existing industrial FE code with the aim of bringing human implementation effort to the minimal possible level (only using standard input-output of finite element analysis (FEA) codes, avoid code-dependent subroutines implementation). An approximate global link to go from Lagrange polynomials to non-uniform-rational-B-splines functions is formulated, which enables the whole FE routines to be untouched during the implementation.

As a result, only the linear system resolution step is bypassed: the resolution is performed in an external script after projecting the FE system onto the reduced, more regular and isogeometric basis. The novel procedure is successfully validated through different numerical experiments involving linear and nonlinear isogeometric analyses using the standard input/output of the industrial FE software Code_Aster.

A non-invasive implementation of IGA into FEA software is proposed. The whole FE routines are untouched during the novel implementation procedure; a focus is made on the IGA solution of nonlinear problems from existing FEA software; technical details on the approach are provided by means of illustrative examples and step-by-step implementation; the methodology is evaluated on a range of two- and three-dimensional elasticity and elastoplasticity benchmarks solved using the commercial software Code_Aster.

As one of the earliest high-level programming languages, Fortran with easy accessibility and computational efficiency is widely used in the engineering field. The purpose of this paper is to present a Fortran implementation of isogeometric analysis (IGA) for thin plate problems.

IGA based on non-uniform rational B-splines (NURBS) offers exact geometries and is more accurate than finite element analysis (FEA). Unlike the basis functions in FEA, NURBS basis functions are non-interpolated. Hence, the penalty method is used to enforce boundary conditions.

Several thin plate examples based on the Kirchhoff-Love theory were illustrated to demonstrate the accuracy of the implementation in contrast with analytical solutions, and the efficiency was validated in comparison with another open method.

A Fortran implementation of NURBS-based IGA was developed to solve Kirchhoff-Love plate problems. It easily obtained high-continuity basis functions, which are necessary for Kirchhoff formulation. In comparison with theoretical solutions, the numerical examples demonstrated higher accuracy and faster convergence of the Fortran implementation. The Fortran implementation can well solve the time-consuming problem, and it was validated by the time-consumption comparison with the Matlab implementation. Due to the non-interpolation of NURBS, the penalty method was used to impose boundary conditions. A suggestion of the selection of penalty coefficients was given.

The purpose of this paper is to propose an efficient and accurate numerical model that employs isogeometric analysis (IGA) for the geometrically nonlinear analysis of functionally graded plates (FGPs). This model is utilized to investigate the effects of boundary conditions, gradient index, and geometric shape on the nonlinear responses of FGPs.

A geometrically nonlinear analysis of thin and moderately thick functionally graded ceramic-metal plates based on IGA in conjunction with first-order shear deformation theory and von Kármán strains is presented. The displacement fields and geometric description are approximated with nonuniform rational B-splines (NURBS) basis functions. The Newton-Raphson iterative scheme is employed to solve the nonlinear equation system. Material properties are assumed to vary along the thickness direction with a power law distribution of the volume fraction of the constituents.

The present model for analysis of the geometrically nonlinear behavior of thin and moderately thick FGPs exhibited high accuracy. The shear locking phenomenon is avoided without extra numerical efforts when cubic or high-order NURBS basis functions are utilized.

This paper shows that IGA is particularly well suited for the geometrically nonlinear analysis of plates because of its exact geometrical modelling and high-order continuity.

This paper aims to propose a gradient-based shape optimization framework in which traditional time-consuming conversions between computer-aided design and computer-aided engineering and the mesh update procedure are avoided/eliminated. The scheme is general so that it can be used in all cases as a black box, no matter what the objective and/or design variables are, whilst the efficiency and accuracy are guaranteed.

The authors integrated CAD and CAE by using isogeometric analysis (IGA), enabling the present methodology to be robust and accurate. To overcome the difficulty in evaluating the sensitivities of objective and/or constraint functions by analytic method in some cases, the authors adopt the finite difference method to calculate these sensitivities, thereby providing a universal approach. Moreover, to further eliminate the inefficiency caused by the finite difference method, the authors advance the exact reanalysis method, the indirect factorization updating (IFU), to exactly and efficiently calculate functions and their sensitivities, which guarantees its generality and efficiency at the same time.

The proposed isogeometric gradient-based shape optimization using our IFU approach is reliable and accurate, as well as general and efficient.

The authors proposed a gradient-based shape optimization framework in which they first integrate IGA and the proposed exact reanalysis method for applicability to structural response and sensitivity analysis.

This study aims to develop a new analysis approach devised by incorporating a gradient-enhanced microplane damage model (GeMpDM) into isogeometric analysis (IGA), which shows computational stability and capability in accurately predicting crack propagations in structures with complex geometries.

For the non-local microplane damage modeling, the maximum modified von-Mises equivalent strain among all microplanes is regularized as a representative quantity. This characterization implies that only one additional governing equation is considered, which improves computational efficiency dramatically. By combined use of GeMpDM and IGA, quasi-static and dynamic numerical analyses are conducted to demonstrate the capability in predicting crack paths of complex geometries in comparison to FEM and experimental results.

The implicit scheme with the adopted damage model shows favorable numerical stability and the numerical results exhibit appropriate convergence characteristics concerning the mesh size. The damage evolution is successfully controlled by a tension-compression damage factor. Thanks to the advanced geometric design capability of IGA, the details of crack patterns can be predicted reliably, which are somewhat difficult to be acquired by FEM. Additionally, the damage distribution obtained in the dynamic analysis is in close agreement with experimental results.

The paper originally incorporates GeMpDM into IGA. Especially, only one non-local variable is considered besides the displacement field, which improves the computational efficiency and favorable convergence characteristics within the IGA framework. Also, enjoying the geometric design ability of IGA, the proposed analysis method is capable of accurately predicting crack paths reflecting the complex geometries of target structures.

This paper proposes a novel approach to impose the Neumann boundary condition for isogeometric analysis (IGA) of Euler–Bernoulli beam with 1-D formulation. The proposed method is for only IGA in which it is difficult to handle the Neumann boundary conditions. The control points of B-spline are equivalent to nodes in finite element method. With 1-D formulation, it is not possible to accommodate multiple degrees of freedom in IGA. This case arises in the analysis of beams. The paper aims to propose a way to work around this issue in a simple way.

Neumann boundary conditions, which are even-order derivatives (example: double derivative) of the primary variable, are inherently satisfied in the weak form. Boundary conditions with an odd number of derivatives (example: slope) are imposed with the introduction of a new penalty matrix.

The proposed method can impose a slope boundary condition for IGA of a beam using 1-D formulation.

From the literature, it can be observed that the beam is formulated in 1-D by considering it as either a rotation-free element or a 2-D formulation by considering shear strain along with the normal strain. The work represents 1-D formulation of a beam while considering the slope boundary condition, which is easy and effective to formulate, compared with the slope boundary conditions reported in previous works.

This paper aims to propose a new boundary condition and a web-spline basis of finite element space approximation to remedy the problems of constraints due to homogeneous and non-homogeneous; Dirichlet boundary conditions. This paper considered the two-dimensional linear elasticity equation of Navier–Lamé with the condition CAB. The latter allows to have a total insertion of the essential boundary condition in the linear system obtained; without using a numerical method as Lagrange multiplier. This study have developed mixed finite element; method using the B-splines Web-spline space. These provide an exact implementation of the homogeneous; Dirichlet boundary conditions, which removes the constraints caused by the standard; conditions. This paper showed the existence and the uniqueness of the weak solution, as well as the convergence of the numerical solution for the quadratic case are proved. The weighted extended B-spline; approach have become a much more workmanlike solution.

In this paper, this study used the implementation of weighted finite element methods to solve the Navier–Lamé system with a new boundary condition CA, B (Koubaiti et al., 2020), that generalises the well-known basis, especially the Dirichlet and the Neumann conditions. The novel proposed boundary condition permits to use a single Matlab code, which summarises all kind of boundary conditions encountered in the system. By using this model is possible to save time and programming recourses while reap several programs in a single directory.

The results have shown that the Web-spline-based quadratic-linear finite elements satisfy the inf–sup condition, which is necessary for existence and uniqueness of the solution. It was demonstrated by the existence of the discrete solution. A full convergence was established using the numerical solution for the quadratic case. Due to limited regularity of the Navier–Lamé problem, it will not change by increasing the degree of the Web-spline. The computed relative errors and their rates indicate that they are of order 1/H. Thus, it was provided their theoretical validity for the numerical solution stability. The advantage of this problem that uses the CA, B boundary condition is associated to reduce Matlab programming complexity.

The mixed finite element method is a robust technique to solve difficult challenges from engineering and physical sciences using the partial differential equations. Some of the important applications include structural mechanics, fluid flow, thermodynamics and electromagnetic fields (Zienkiewicz and Taylor, 2000) that are mainly based on the approximation of Lagrange. However, this type of approximation has experienced a great restriction in the level of domain modelling, especially in the case of complicated boundaries such as that in the form of curvilinear graphs. Recently, the research community tried to develop a new way of approximation based on the so-called B-spline that seems to have superior results in solving the engineering problems.

This study aims to present a new hybrid formulation based on non-uniform rational b-splines functions and enrichment strategies applied to free and forced vibration of straight bars and trusses.

Based on the idea of enrichment from generalized finite element method (GFEM)/extended finite element method (XFEM), an extended isogeometric formulation (partition of unity isogeometric analysis [PUIGA]) is conceived. By numerical examples the methods are tested and compared with isogeometric analysis, finite element method and GFEM in terms of convergence, error spectrum, conditioning and adaptivity capacity.

The results show a high convergence rate and accuracy for PUIGA and the advantage of input enrichment functions and material parameters on parametric space.

The enrichment strategies demonstrated considerable improvements in numerical solutions. The applications of computer-aided design mapped enrichments applied to structural dynamics are not known in the literature.

This paper aims to describe the application of the virtual element method (VEM) to contact problems between elastic bodies.

Polygonal elements with arbitrary shape allow a stable node-to-node contact enforcement. By adaptively adjusting the polygonal mesh, this methodology is extended to problems undergoing large frictional sliding.

The virtual element is well suited for large deformation contact problems. The issue of element stability for this specific application is discussed, and the capability of the method is demonstrated by means of numerical examples.

This work is completely new as this is the first time, as per the authors’ knowledge, the VEM is applied to large deformation contact.