Nonlinear three–dimensional finite–element modelling of reinforced–concrete beams: Computational challenges and experimental validation

Three–dimensional nonlinear finite–element models have been developed to investigate the loading–unloading–reloading behaviour of two reinforced–concrete beams under four–point bending using explicit dynamics in ABAQUS. The damaged–plasticity model proposed by Lubliner and collaborators was employed for the plain concrete, and elastic–perfectly plastic models were employed for the steel reinforcement. A perfect bond was assumed between the steel rebars and concrete, whereby the bond–slip behaviour, as well as damage along crack patterns, were modelled through concrete damage. The influence of the shape of the tension–softening law on the numerical load–deflection response was studied by considering bi–linear, exponential and linear post–failure stress–displacement and stress–strain relationships. The effect of modelling steel rebars with truss or beam elements was also investigated. Structured meshes of linear hexahedral elements either with incompatible modes or with reduced integration, and unstructured meshes of either linear or ‘modified’ quadratic tetrahedral elements were considered. In terms of load–deflection curves, both the structured and the unstructured meshes gave results in very good agreement with test results. In terms of crack patterns, results predicted by the structured meshes exhibited some mesh bias, which was less pronounced with the unstructured meshes. In the post–yield phase, if a geometrically nonlinear model is used, discrepancies were found when truss elements were used for steel rebars, whereas good agreement was found if the bending stiffness of the rebars is included using beam elements. This is a non–obvious result that may be important to consider when studying the progressive collapse of RC structures.     

Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part II—A framework for volume mesh optimization and the condition number of the Jacobian matrix

Three‐dimensional unstructured tetrahedral and hexahedral finite element mesh optimization is studied from a theoretical perspective and by computer experiments to determine what objective functions are most effective in attaining valid, high‐quality meshes. The approach uses matrices and matrix norms to extend the work in Part I to build suitable 3D objective functions. Because certain matrix norm identities which hold for 2×2 matrices do not hold for 3×3 matrices, significant differences arise between surface and volume mesh optimization objective functions. It is shown, for example, that the equality in two dimensions of the smoothness and condition number of the Jacobian matrix objective functions does not extend to three dimensions and further, that the equality of the Oddy and condition number of the metric tensor objective functions in two dimensions also fails to extend to three dimensions. Matrix norm identities are used to systematically construct dimensionally homogeneous groups of objective functions. The concept of an ideal minimizing matrix is introduced for both hexahedral and tetrahedral elements. Non‐dimensional objective functions having barriers are emphasized as the most logical choice for mesh optimization. The performance of a number of objective functions in improving mesh quality was assessed on a suite of realistic test problems, focusing particularly on all‐hexahedral ‘whisker‐weaved’ meshes. Performance is investigated on both structured and unstructured meshes and on both hexahedral and tetrahedral meshes. Although several objective functions are competitive, the condition number objective function is particularly attractive. The objective functions are closely related to mesh quality measures. To illustrate, it is shown that the condition number metric can be viewed as a new tetrahedral element quality measure. Published in 2000 by John Wiley & Sons, Ltd.     

Three‐dimensional immersed finite element methods for electric field simulation in composite materials

This paper presents two immersed finite element (IFE) methods for solving the elliptic interface problem arising from electric field simulation in composite materials. The meshes used in these IFE methods can be independent of the interface geometry and position; therefore, if desired, a structured mesh such as a Cartesian mesh can be used in an IFE method to simulate 3‐D electric field in a domain with non‐trivial interfaces separating different materials. Numerical examples are provided to demonstrate that the accuracies of these IFE methods are comparable to the standard linear finite element method with unstructured body‐fit mesh. Copyright © 2005 John Wiley & Sons, Ltd.     

A continuum damage mechanics method for fracture simulation of quasi‐brittle materials

This paper addresses a novel continuum damage‐based method for simulating failure process of quasi‐brittle materials starting from local damage initiation to final fracture. In the developed method, the preset characteristic length field is used to evaluate damage instead of element, which is used to reduce the spurious sensitivity. In addition, damage is only updated in the most dangerous location at a time for considering stress redistribution due to damage evolution, which is used to simulate competitive fracture process. As cases study, representative numerical simulations of two benchmark tests are given to verify the performance of the developed continuum damage‐based method together with a used damage model. The simulation results of the crack paths for two concrete specimens obtained from the developed method matched well with the corresponding experimental results. The results show that the developed continuum damage‐based method is effective and can be used to simulate damage and fracture process of brittle or quasi‐brittle materials. And the simulation results based on the developed method depend only the preset characteristic length field and not grid mesh.     

Simulating the propagation of displacement discontinuities in a regularized strain‐softening medium

A numerical model is developed which allows the inclusion of displacement discontinuities in a strain‐softening medium, independent of the finite element mesh structure. Inelastic deformations develop in the continuum and, when a critical threshold of inelastic deformation is reached, a displacement discontinuity is inserted. Discontinuities are introduced using the partition of unity concept which allows discontinuous functions to be added to the standard finite element basis. It is shown that the introduction of displacement discontinuities at the later stages of the failure process can lead to a failure mode that is fundamentally different than that using a continuum model only. This combined continuum‐discontinuous model is better able to describe the entire failure process than a continuum or a discrete model alone and treats mode‐I and mode‐II failure in a unified fashion. Copyright © 2001 John Wiley & Sons, Ltd.     

A mesh objective method for modeling crack propagation using the smeared crack approach

This paper presents a mesh objective method for modeling crack propagation in brittle materials using a conventional finite element formulation. The primary shortcoming of the smeared crack approach is its pathological sensitivity to the mesh orientation, which is manifested by shear locking and stress field misalignment around the crack tip. Such undesirable characteristics preclude the ability to model arbitrary crack propagation at an angle through the mesh. Several techniques are developed to address these shortcomings. First, to preclude shear locking, a modified failure constitutive model is developed, which projects out the spurious stress increments as the crack opens. If a crack exists in an element, a crack tracking algorithm is used to identify the neighboring elements most likely to show crack continuation. This algorithm also identifies a crossover element when a crack passes through adjacent sides of an element. Then, the characteristic element length used in the constitutive equation is changed with the objective of providing the correct failure energy per unit crack length, a procedure called crossover scaling. The examples provided demonstrate that the developed methods work collectively to provide a simple and efficient method for modeling failure in brittle materials without mesh bias.     

小冲孔试样弹塑性蠕变损伤的有限元分析

综合考虑了小冲孔试样实验过程中的弹性、塑性及蠕变等因素,采用了Mises型流动法则和改进的K-R蠕变损伤本构方程,建立了ABAQUS有限元模型,模拟计算了550℃三种载荷条件下,Cr5Mo试样的弹塑性蠕变过程,分析试样的中心蠕变挠度、蠕变应变随时间的变化规律,当量应力的分布情况,以及试样的蠕变损伤与位置的关系,同时分析了挠度与应力对网格的敏感性。结果表明:模拟结果与实验结果很接近,且与单轴实验的应变曲线十分相似;整个试样在中心损伤比较严重,呈现明显局部化特征;讨论了试样网格的敏感性,试样半径方向0mm~1mm之内应力变化幅度比较大(除约束点),并且应力的最大值也出现在这个区域,在这个区域进行网格细化提高计算的精确性,而其他区域采用粗网格减少计算成本,达到精确性和经济性的平衡。径向上网格大小对挠度影响比较大,应力对网格大小不敏感。     

A new method for modelling cohesive cracks using finite elements

A model which allows the introduction of displacements jumps to conventional finite elements is developed. The path of the discontinuity is completely independent of the mesh structure. Unlike so‐called ‘embedded discontinuity’ models, which are based on incompatible strain modes, there is no restriction on the type of underlying solid finite element that can be used and displacement jumps are continuous across element boundaries. Using finite element shape functions as partitions of unity, the displacement jump across a crack is represented by extra degrees of freedom at existing nodes. To model fracture in quasi‐brittle heterogeneous materials, a cohesive crack model is used. Numerical simulations illustrate the ability of the method to objectively simulate fracture with unstructured meshes. Copyright © 2001 John Wiley & Sons, Ltd.     

Automatic generation of multiblock decompositions of surfaces

Multiblock‐structured meshes have significant advantages over fully unstructured meshes in numerical simulation, but automatically generating these meshes is considerably more difficult. A method is described herein for automatically generating high‐quality multiblock decompositions of surfaces with boundaries. Controllability and flexibility are useful capabilities of the method. Additional alignment constraints for forcing the appearance of particular features in the decomposition can be easily handled. Also, adjustments are made according to input metric tensor fields that describe target element size properties. The general solution strategy is based around using a four‐way symmetry vector‐field, called a cross‐field, to describe the local mesh orientation on a triangulation of the surface. Initialisation is performed by propagating the boundary alignment constraints to the interior in a fast marching method. This is similar in a way to an advancing‐front or paving method but much more straightforward and flexible because mesh connectivity does not have to be managed in the cross‐field. Multiblock decompositions are generated by tracing the separatrices of the cross‐field to partition the surface into quadrilateral blocks with square corners. The final task of meshing the decomposition requires solving an integer programming problem for block division numbers. Copyright © 2015 John Wiley & Sons, Ltd.     

The generation of hexahedral meshes for assembly geometry: survey and progress

The finite element method is being used today to model component assemblies in a wide variety of application areas, including structural mechanics, fluid simulations, and others. Generating hexahedral meshes for these assemblies usually requires the use of geometry decomposition, with different meshing algorithms applied to different regions. While the primary motivation for this approach remains the lack of an automatic, reliable all‐hexahedral meshing algorithm, requirements in mesh quality and mesh configuration for typical analyses are also factors. For these reasons, this approach is also sometimes required when producing other types of unstructured meshes. This paper will review progress to date in automating many parts of the hex meshing process, which has halved the time to produce all‐hex meshes for large assemblies. Particular issues which have been exposed due to this progress will also be discussed, along with their applicability to the general unstructured meshing problem. Published in 2001 by John Wiley & Sons, Ltd.     

Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part I—a framework for surface mesh optimization

Structured mesh quality optimization methods are extended to optimization of unstructured triangular, quadrilateral, and mixed finite element meshes. New interpretations of well‐known nodally based objective functions are made possible using matrices and matrix norms. The matrix perspective also suggests several new objective functions. Particularly significant is the interpretation of the Oddy metric and the smoothness objective functions in terms of the condition number of the metric tensor and Jacobian matrix, respectively. Objective functions are grouped according to dimensionality to form weighted combinations. A simple unconstrained local optimum is computed using a modified Newton iteration. The optimization approach was implemented in the CUBIT mesh generation code and tested on several problems. Results were compared against several standard element‐based quality measures to demonstrate that good mesh quality can be achieved with nodally based objective functions. Published in 2000 by John Wiley & Sons, Ltd.     

Shape design sensitivities improvement using convected unstructured meshes

This paper focuses on various forms of direct differentiation methods for design sensitivity computation in the shape optimisation of continuum structures and the role of convected meshes on the accuracy of the sensitivities. A Pseudo-Analytical Sensitivity Analysis (P-ASA) method is presented and tested. In this method the response analysis component uses unstructured finite element meshes and the sensitivity algorithm entails shape-perturbation for each design variable. A material point is convected during a change of shape and the design sensitivities are therefore intrinsically associated with the mesh-sensitivities of the finite element discretization. Such mesh sensitivities are obtained using a very efficient boundary element point-tracking analysis of an affine notional underlying elastic domain. All of the differentiation, with respect to shape variables, is done exactly except for the case of mesh-sensitivities: hence the method is almost analytical. In contrast to many other competing methods, the P-ASA method is, by definition, independent of perturbation step-size, making it particularly robust. Furthermore, the sensitivity accuracy improves with mesh refinement. The boundary element point-tracking method is also combined with two popular methods of sensitivity computation, namely the global finite difference method and the semi-analytical method. Increases in accuracy and perturbation range are observed for both methods.     

An algebraic two‐level preconditioner for asymmetric,positive‐definite systems

A two‐level, linear algebraic solver for asymmetric, positive‐definite systems is developed using matrices arising from stabilized finite element formulations to motivate the approach. Supported by an analysis of a representative smoother, the parent space is divided into oscillatory and smooth subspaces according to the eigenvectors of the associated normal system. Using a mesh‐based aggregation technique, which relies only on information contained in the matrix, a restriction/prolongation operator is constructed. Various numerical examples, on both structured and unstructured meshes, are performed using the two‐level cycle as the basis for a preconditioner. Results demonstrate the complementarity between the smoother and the coarse‐level correction as well as convergence rates that are nearly independent of the problem size. Copyright © 2001 John Wiley & Sons, Ltd.     

Development of physically based meshes for two‐dimensional models of meandering channel flow

The choice of mesh generation and numerical solution strategies for two‐dimensional finite element models of fluvial flow have previously been based chiefly on experience and rule of thumb. This paper develops a rationale for the finite element modelling of flow in river channels, based on a study of flow around an annular reach. Analytical solutions of the two‐dimensional Shallow Water (St. Venant) equations are developed in plane polar co‐ordinates, and a comparison with results obtained from the TELEMAC‐2‐D finite element model indicates that of the two numerical schemes for the advection terms tested, a flux conservative transport scheme gives better results than a streamline upwind Petrov–Galerkin technique. In terms of mesh discretization, the element angular deviation is found to be the most significant control on the accuracy of the finite element solutions. A structured channel mesh generator is therefore developed which takes local channel curvature into account in the meshing process. Results indicate that simulations using curvature‐dependent meshes offer similar levels of accuracy to finer meshes made up of elements of uniform length, with the added advantage of improved model mass conservation in regions of high channel curvature. Copyright © 2000 John Wiley & Sons, Ltd.     

A parallel implicit/explicit hybrid time domain method for computational electromagnetics

The numerical solution of Maxwell's curl equations in the time domain is achieved by combining an unstructured mesh finite element algorithm with a cartesian finite difference method. The practical problem area selected to illustrate the application of the approach is the simulation of three‐dimensional electromagnetic wave scattering. The scattering obstacle and the free space region immediately adjacent to it are discretized using an unstructured mesh of linear tetrahedral elements. The remainder of the computational domain is filled with a regular cartesian mesh. These two meshes are overlapped to create a hybrid mesh for the numerical solution. On the cartesian mesh, an explicit finite difference method is adopted and an implicit/explicit finite element formulation is employed on the unstructured mesh. This approach ensures that computational efficiency is maintained if, for any reason, the generated unstructured mesh contains elements of a size much smaller than that required for accurate wave propagation. A perfectly matched layer is added at the artificial far field boundary, created by the truncation of the physical domain prior to the numerical solution. The complete solution approach is parallelized, to enable large‐scale simulations to be effectively performed. Examples are included to demonstrate the numerical performance that can be achieved. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluation of directional mesh bias in concrete fracture simulations using continuum damage models

In the present comparative study, we investigate the influence of directional mesh bias on the results of failure simulations performed with isotropic and anisotropic damage models. Several fracture tests leading to curved crack trajectories are simulated on different meshes. The isotropic damage model with a realistic biaxial strength envelope for concrete is highly sensitive to the mesh orientation, even for fine meshes. The sensitivity is reduced if the definition of the damage-driving variable (equivalent strain) is based on the modified von Mises criterion, but the corresponding biaxial strength envelope is not realistic for concrete. The anisotropic damage models used in this study capture reasonably well arbitrary crack trajectories even if the biaxial strength envelope remains close to typical experimental data. Their superior performance can be at least partially attributed to their ability to capture dilatancy under shear, which is revealed by a comparative analysis of the behavior of individual models under shear with restricted or free volume expansion.     

CFD-based optimization of hovering rotors using radial basis functions for shape parameterization and mesh deformation

Aerodynamic shape optimization of a helicopter rotor in hover is presented, using compressible CFD as the aerodynamic model. An efficient domain element shape parameterization method is used as the surface control and deformation method, and is linked to a radial basis function global interpolation, to provide direct transfer of domain element movements into deformations of the design surface and the CFD volume mesh, and so both the geometry control and volume mesh deformation problems are solved simultaneously. This method is independent of mesh type (structured or unstructured) or size, and optimization independence from the flow solver is achieved by obtaining sensitivity information for an advanced parallel gradient-based algorithm by finite-difference, resulting in a flexible method of ‘wrap-around’ optimization. This paper presents results of the method applied to hovering rotors using local and global design parameters, allowing a large geometric design space. Results are presented for two transonic tip Mach numbers, with minimum torque as the objective, and strict constraints applied on thrust, internal volume and root moments. This is believed to be the first free form design optimization of a rotor blade using compressible CFD as the aerodynamic model, and large geometric deformations are demonstrated, resulting in significant torque reductions, with off-design performance also improved.     

Finite element solution of free‐surface ship‐wave problems

An unstructured finite element solver to evaluate the ship‐wave problem is presented. The scheme uses a non‐structured finite element algorithm for the Euler or Navier–Stokes flow as for the free‐surface boundary problem. The incompressible flow equations are solved via a fractional step method whereas the non‐linear free‐surface equation is solved via a reference surface which allows fixed and moving meshes. A new non‐structured stabilized approximation is used to eliminate spurious numerical oscillations of the free surface. Copyright © 1999 John Wiley & Sons, Ltd.