Search

In this paper, a finite element procedure is used to study the dynamic behaviour of a paper web in a free span between two rollers, including effects of transport velocity and surrounding air. The paper web is modelled as a three-dimensional orthotropic structure. The influence of air is accounted for by utilizing fluid–solid interaction analyses based on acoustic theory. The contribution of transport velocity is included through gyroscopic matrices and forces. The structural response on harmonic excitations has been studied using linear and non-linear models. Results show that air significantly reduces eigenfrequencies of the web. So called “edge-flutter” is nothing but the result of skew tension profile. Excessive web vibration can be eliminated by adjusting the web tension. [Copyright &y& Elsevier]

This paper presents the formulation of a finite element method for nonlinear Kirchhoff rods (i.e. without transverse shear strain) in the general three-dimensional setting defined by a Cosserat director treatment of the cross sections attached to the rod's axis. The new element is based on a G 1 interpolation of the rod's geometry in terms of Hermite shape functions of the rod's axis (including its tangent defining the tangential director), while the transversal directors defining the different bending and torsional responses of the rod consider a Lagrangian interpolation of the section directors. This direct interpolation of the directors, as opposed of underlying rotation vectors, assures the objectivity of the proposed formulation. In fact, the invariance properties of the resulting finite element are analyzed in detail, assuring the correct resolution of the local fundamental equilibrium relations between forces and moments, hence avoiding the so-called "self-straining" associated to separate treatments of the rod's geometry and its kinematics. Several representative numerical simulations are presented illustrating these properties as well as the appropriateness of the proposed formulation for the analysis of thin rods undergoing large finite deformations in the three-dimensional range. • The paper develops a new finite element for Kirchhoff rods. • Such rods involve no transverse shear, leading to a high-order formulation. • The fully nonlinear geometrically exact range in the 3D setting is considered. • A detailed analysis of the invariance of the proposed finite element is presented. • Numerical tests compare with analytical solutions and other literature proposals. [ABSTRACT FROM AUTHOR]

The present paper introduces a global/local approach for the analysis of three-dimensional (3D) stress states of composite laminated structures. It consists of a two-step procedure. In particular, the first step makes use of finite element modeling based on classical 2D plate elements, whereas a refined layer-wise model based on Carrera Unified Formulation (CUF) is adopted to extract the 3D stress and strain fields in some critical regions that may have arbitrary dimensions. This approach allows dealing with large local areas, increasing the accuracy of the static solution, along with the possibility of embedding this technique in more complex procedures, such as the least-weight design of large heterogeneous complex assemblies, stiffness optimization, or localized progressive failure analysis. The numerical results shown in this paper want to assess the physical and numerical validity of the global/local approach. Particular attention is focused on the choice of the dimensions of the local area (patch) subjected to detailed refined analysis. Also, the convergence properties of the present hybrid FE model are discussed. The first examples deal with laminated composites. Then, the advantages of this methodology are further highlighted by considering free-edge problems and a complex wing structure. • Global to local 3D stress state analysis of composite structures. • Classical finite element software tool is coupled with layerwise CUF formulation. • The method can deal with failure onset mechanisms, including free-edge effects. • A full wing structure is analyzed. [ABSTRACT FROM AUTHOR]

The paper presents a numerical method for simulating flow and mechanics in fractured rock. The governing equations that couple the effects in the rock mass and in the fractures are obtained using the discrete fracture-matrix approach. The fracture flow is driven by the cubic law, and the contact conditions prevent fractures from self-penetration. A stable finite element discretization is proposed for the displacement-pressure-flux formulation. The resulting nonlinear algebraic system of equations and inequalities is decoupled using a robust iterative splitting into the linearized flow subproblem, and the quadratic programming problem for the mechanical part. The non-penetration conditions are solved by means of dualization and an optimal quadratic programming algorithm. The capability of the numerical scheme is demonstrated on a benchmark problem for tunnel excavation with hundreds of fractures in 3D. The paper's novelty consists in a combination of three crucial ingredients: (i) application of discrete fracture-matrix approach to poroelasticity, (ii) robust iterative splitting of resulting nonlinear algebraic system working for real-world 3D problems, and (iii) efficient solution of its mechanical quadratic programming part with a large number of fractures in mutual contact by means of own solvers implemented into an in-house software library. • A numerical method for accurate resolution of flow and mechanics in fractured rock is presented. • Discrete fracture-matrix approach links processes in rock mass and fracture network. • Robust iterative splitting enables solving nonlinear hydro-mechanical problems. • Efficient algorithm for quadratic programming solves non-penetration conditions on fractures. • Parallel implementation is capable of handling complex geometries with hundreds of fractures in 3D. [ABSTRACT FROM AUTHOR]

• A novel surrogate-based crack identification method is presented. • Structure frequency shift characteristics is analyzed. • The crack identification is converted to an optimization problem. • Numerical examples and experimental validation prove satisfactory precision. A novel crack identification method is presented in this paper for cantilever beam-type structures within the scope of vibration-based damage identification. Natural frequencies are more easily to be measured yet are weakly sensitive to damage in comparison to other dynamic parameters such as mode shape and damping ratio. Therefore, this paper develops a surrogate-based crack identification method in order to overcome the low sensitivity of frequency to damage. The surrogate model was constructed using radial point interpolation method (RPIM) and subsequently trained with the frequency shift database established from parametric finite element analysis (FEA) modeling of cracked/ non-cracked cantilever beams. On this basis, the crack identification problem had been converted to a standard optimization model with crack length and location as design variables. Subsequently, genetic algorithm (GA) was employed for solving the optimization model to achieve crack parameters. In the end, numerical examples and experimental test have been conducted to validate the developed method, revealing maximum relative errors within 0.5% and 2.5%, respectively. [ABSTRACT FROM AUTHOR]

In this paper, the statistical finite element method is developed further to synthesize distributed rail response data with nonlinear finite element model predictions within and outside the measured loading range. In the data-generating model of the statistical finite element method, the distributed sensing data is decomposed into a finite element model component, a model-reality mismatch component, and a noise component. Each component is considered uncertain and is represented as a Gaussian random vector with a corresponding prior density. The finite element prior density is updated using the Bayesian statistical framework in light of the distributed fiber optic sensing data. The calculated posterior density enables one to infer the true structural response. The required finite element prior density is determined by solving a conventional stochastic forward problem using a polynomial chaos expansion of random fields and a non-intrusive pseudo-spectral projection approach. In this paper, a lab test involving a section of a rail (i.e., a beam-column structural member) instrumented with distributed fiber optic sensors and subjected to axial load causing progressive lateral buckling is considered to demonstrate the use of the statistical finite element method to improve rail response prediction. The results show improved prediction of true rail strain in linear and nonlinear regimes compared with a pure finite element model-based prediction. • Statistical finite element method enables integration of nonlinear model and sensor data. • Model selection/misspecification & data uncertainty incorporated using Bayesian formalism. • A polynomial chaos expansion is utilized to derive the necessary finite element prior. • The proposed approach is validated with progressive buckling tests and strain data. [ABSTRACT FROM AUTHOR]

This paper aims for developing topology optimization methodology to design the shape of electrodes in Dielectrophoresis (DEP)-based devices. The DEP force is due to a non-uniform electric field induced by applied voltages to the electrodes. Shape of the electrodes has the principal effect on the direction and magnitude of the DEP force. In medical therapy microfluidic devices, DEP force is used for cell sorting and cell separation. While the direction and magnitude of the DEP force are desired to be determined and maximized respectively, the magnitude of the electric field should be minimized to avoid damaging cells. Approaching these goals is counter intuitive where the existing electrode designs are basic. Therefore, a detailed finite element model (FEM) is developed for DEP force and electric field to formulate an optimization problem to maximize the DEP force in a particular direction while there is a constraint on electric field's magnitude. Using the developed FEM, explicit formulations for sensitivity analysis are derived to implement a gradient-based topology optimization. The performance of developed methodology is assessed numerically to determine the direction of the DEP force and constraining the electric field and experimentally in a practical case study of particle trapping in a microfluidic channel. • Topology optimization can produce novel electrode shapes for DEP devices. • Sensitivity analysis is provided for topology optimization of electrodes in DEP devices. • Optimal electrodes can determine direction of DEP force and minimize electric field. • 2D optimal electrodes are efficient in real 3D microfluidic applications. • Optimal electrodes are efficient for particle manipulation in a micro channel. [ABSTRACT FROM AUTHOR]

In this paper, we propose a novel thermal deformation formulation for elastic bodies' thermomechanical simulation with total Lagrangian smoothed particle hydrodynamics (TLSPH), which employing fixed reference configuration. We utilize implicit corotated TLSPH as a base methodology to address tensile instability and enhance the robustness of simulation. The reference configuration of heat transfer model is updated, but the reference configuration of thermal deformation model is fixed like TLSPH. As particle's temperature is changed, the inter-particle distance in reference configuration is updated by using linear thermal expansion theory. Validation involved two basic cantilever beam loading cases and a thermomechanical analysis with heat-induced thermal deformation on a cantilever beam. The TLSPH simulation results were compared against Finite Element Method simulations and analytical solutions, demonstrating equivalent accuracy across all tests. This study indicates the potential for conducting a wide array of simulations involving thermal deformation using the TLSPH method, which efficiently supports parallel computing on an individual particle basis. • An implicit corotated total Lagrangian smoothed particle hydrodynamics (SPH) is employed for simulating elastic bodies. • A novel methodology for thermal deformation, based on the TLSPH framework, is introduced. • To verify basic structural analyses, dynamic oscillation simulations are performed. • For thermomechanical validation, the simulations analyzing heat transfer and thermal deformation are carried out. [ABSTRACT FROM AUTHOR]

• This Shape Optimization allows to achieve Uniform Strength for entire Structure. • The Uniform Strength Beam shape is derived analytically for whatever type of load. • Proposed procedure work with a generic section with one axis of symmetry. • Optimization as a search for zeros of kinematic-congruence and balance conditions. • Optimization workflow given fast results also for very large-size Frame Structures. This paper provides a procedure to obtain the uniform strength of frame and lattice structures. Uniform strength condition is achieved by performing the shape optimization of all beam elements of the structure. The beam shape which guarantees uniform strength is analytically deduced from the one-dimensional Timoshenko model. The optimization problem presents itself as the search for the zeros of the objective-functions vector, which is a non-linear system of equations representing the kinematic-congruence and forces balance at every node of the structure. The analytical formulation of the optimization problem allows to construct the objective-functions vector without the use of external structural computation, i.e. not recurring to any Finite Element Analysis to accomplish iterations. This latter feature entails a great advantage in terms of computing time required to perform optimization. The proposed analytical formulation allows to directly insert the uniform strength condition into the objective-functions vector, transforming the optimization into an unconstrained problem. Some examples are shown in which the performance of the optimization procedure is discussed in terms of robustness and rate of computational complexity while increasing the degrees of freedom of the structure. The reliability and the quality of the optimization are verified through Finite Element Analysis. [ABSTRACT FROM AUTHOR]

• Propose a neural networks-based line element (NNLE) method for large deflection frame analysis. • Significantly outperform the line finite element method in efficiency. • Perform great modeling flexibility and compatibility with the exiting numerical approaches. • Simplify the implementation of machine learning (ML) in large-scale structural analysis. • Provide a novel strategy for seamlessly integrating finite element methods and ML techniques. The line finite element method (LFEM) is the predominant simulation method in structural design due to its robustness in large-scale structural analysis. However, it sometimes suffers from the tedious computational process due to its fine-mesh requirement to ensure accuracy. The machine learning (ML) technique provides an efficient mesh-free alternative but necessitating tremendous training datasets for modeling large-scale structural systems. In this paper, a novel numerical framework, named the neural networks-based line element (NNLE) method, synergizing the unique advantages of the finite element method and ML technique, is proposed and presented within the context of large deflection frame analysis. The neural networks (NN) model is only trained for modeling single components, thereby significantly diminishing the model scale and the required training dataset. Then, the NN model is used to formulate a new NNLE and implemented within the existing LFEM framework to simulate the entire structural system. Extensive examples are performed to demonstrate the accuracy, efficiency, compatibility, and flexibility of the proposed NNLE method compared with the conventional LFEM and ML techniques. It is convinced that the proposed NNLE method will offer new insights into the combination of the traditional finite element method and the emerging ML approach. [ABSTRACT FROM AUTHOR]

This paper introduces a formulation of the robust topology optimization problem that is tailored for designing fiber-reinforced composite structures with spatially varying principal mechanical properties. Specifically, a methodology is developed that incorporates the spatial variability in the engineering constants of the composite lamina into the concurrent topology (i.e., material distribution) and morphology (i.e., fiber orientation distribution) optimization problem for the minimization of the robust compliance function. The spatial variability in the mechanical properties of the lamina is modeled as a homogeneous random field within the design domain by means of the Karhunen-Lo e ´ ve series expansion, and is thereafter intrusively propagated into the stochastic finite element analysis of the composite structure. To carry out the stochastic finite element analysis per iteration of the optimization cycle, the first-order perturbation method is utilized for approximating the current state variables of the physical system. The resulting robust topology and fiber orientation optimization problem is formulated step-by-step for the minimization of the robust compliance function. With the view of solving the optimization problem at hand by means of gradient-based solution algorithms, the first-order derivatives of the involved design functions w. r. t. the associated design variables are analytically derived. The present work concludes with a series of numerical examples, focusing on the benchmark academic case studies of the 2D cantilever and the half part of the Messerschmitt-Bölkow-Blohm beam, aiming to demonstrate the developed methodology as well as to explore the effect that different parameterization instances of the random field bear on the predicted topology and morphology of the beams. • A methodology for designing robust composites under material uncertainty is presented. • The impact of the material uncertainty level on the design's topology is examined. • Higher uncertainty levels are shown to result in refined topological designs. [ABSTRACT FROM AUTHOR]

• Quadrilateral overlapping finite elements are developed to solve incompressible fluid flows; • Novel stable pairs of elements that satisfy the inf-sup condition are proposed; • The stability of the elements is tested using the inf-sup test; • The elements are polynomial-based, insensitive to distortions and for use in AMORE; • The element interpolations can also be used in incompressible solid mechanics. The purpose of this paper is to investigate the performance of a certain family of low-order quadrilateral finite elements in the solution of the stationary Navier-Stokes equations governing incompressible fluid flows. These finite elements are derived from the 'overlapping finite elements', first developed for the solution of problems in solid mechanics [5,9,48] and incorporate features of both meshfree and traditional finite element methods. One of their most remarkable properties is the insensitivity to mesh distortions. Also, since the Shepard functions are replaced by a suitable interpolation, the resulting basis functions are entirely polynomial, which allows the numerical integration of the weak forms to be performed with few integration points [6,49]. We also discuss the theoretical reasons for the stability of the solutions, that is, the absence of spurious pressure modes, and show that the proposed discretization scheme passes the relevant inf-sup test. [ABSTRACT FROM AUTHOR]

• The collective performance of Masonry Aggregates affects the individual units. • A novel BIM to FEM framework for Masonry Aggregates is outlined. • The Level of Information Need for solid FEA of Masonry Aggregates is defined. • An automated process from IFC to tetrahedron finite elements is developed. • The developed methods follow Open BIM logic. Masonry building aggregates are complex historical structures comprising interconnected buildings. The safety of each building depends not only on its structural integrity but also on the interactions with neighbouring structures and their collective performance. The structural assessment of such complex constructions can benefit from BIM methodology. The latter eases modelling complicated solid geometries and assigning mechanical properties information essential for developing a reliable finite element model. Thus far, BIM to finite element analysis procedures applied to historic constructions have remained laborious, semi-automatic, and based on proprietary software. This paper introduces an open-source automated solid finite element analysis method using OpenBIM data (IFC). The tool implemented automatically creates tetrahedron mesh and assigns mechanical properties, self-weights, and fixities at the structure's base. The methodology is applied to an actual masonry building aggregate to validate its capability of working with complex geometries and scalability. [ABSTRACT FROM AUTHOR]

Highlights • A new and robust numerical method based on use of isoparametric elements is proposed for solving general thermal and mechanical problems. • Only one element is used for each collocation node and the element can be freely formed by itself and surround nodes. • The bandwidth of the formed final system coefficient matrix is extremely narrow. • A new 21-node quadratic element for 3D problems is constructed in this paper for the first time. Abstract In this paper, a new numerical method, named as the Free Element Collocation Method (FECM), is proposed for solving general engineering problems governed by the second order partial differential equations (PDEs). The method belongs to the group of the collocation method, but the spatial partial derivatives of physical quantities are computed based on the isoparametric elements as used in FEM. The key point of the method is that the isoparametric elements used can be freely formed by the nodes around the collocation node. To achieve a narrow bandwidth of the final system of equations, elements with a central node are recommended. For this purpose, a new 21-node quadratic element for 3D problems is constructed for the first time. Attributed to the use of isoparametric elements, FECM can result in more stable results than the traditional collocation method. In addition, the elements can be freely formed by local nodes, FECM has the advantage of mesh-free methods to fit complicated geometries of engineering problems. A number of numerical examples of 2D and 3D thermal and mechanical problems are given to demonstrate the correctness and efficiency of the proposed method. [ABSTRACT FROM AUTHOR]

• A methodology for analysis of large-scale masonry structures is presented which incorporates an inelastic macroscopic framework coupled with a series of artificial neural networks. • The neural networks serve the purpose of identification of material functions embedded in anisotropic formulation as well as the specification of localization plane. • The data required for training of neural networks is generated using 'virtual experiments' that involve a mesoscale finite element analysis of masonry wallets. • Finite element analysis of a large masonry wall with multiple openings is carried out. The results are compared with those based on a detailed mesoscale model. This paper presents a methodology for analysis of large-scale masonry structures. The approach involves development of a series of artificial neural networks which enable the identification of main variables employed in the macroscopic formulation that incorporates an inelastic constitutive law with embedded discontinuity. The data required for training of neural networks is generated using 'virtual experiments', whereby the 'equivalent' anisotropic response of masonry is obtained through a mesoscale finite element analysis of masonry wallets. The paper outlines the procedure for identification of approximation coefficients describing the orientation-dependency of strength, and other relevant parameters. A numerical example is provided involving analysis of a large masonry wall with multiple openings. The results of macroscale approach are compared with those based on a detailed mesoscale model for the same geometry and boundary conditions. [ABSTRACT FROM AUTHOR]

• An energy-based homogenization method (EBHM) for the lattice structures with generalized periodicity is developed. • The limitation of classical EBHM, which is only applicable to unit cells with a square shape is overcome. • The closed-form expressions are derived for equivalent properties in terms of the elastic strain energy of the unit cell. • The effect of temperature rise operating conditions on the accuracy of calculations is highlighted. • The computational accuracy of the EBHM, classical homogenization method and direct finite element analysis are investigated. This paper presents an energy-based homogenization method (EBHM) to calculate the equivalent elastic properties of lattice structures with generalized periodicity. Unlike the traditional implementation of the homogenization method, expressions of closed-form are derived for the equivalent elastic matrix, equivalent coefficients of thermal stress and thermal expansion in terms of the elastic strain energy of the unit cell so that the tedious numerical solution and programming are avoided. It is shown that the elastic strain energy can easily be calculated by mapping the unit cell with the imposition of specific periodic boundary conditions. The implementation can resort to any available finite element tools. Numerical examples are used to compare the EBHM with the homogenization mapping method, classical homogenization method and direct finite element analysis (FEA). The computational accuracy is investigated to show the effectiveness of the EBHM. [ABSTRACT FROM AUTHOR]

• The innovative HFT approach offers a practical solution to avoid expensive large fire tests. • Achieves quick convergence for accurate structural analysis in HFT. • TCP/IP platform ensures continuous data exchange between physical and numerical substructures. • A user-friendly software implemented in C#, enabling seamless HFT execution, validated in a 3-story steel structure test. In large-scale structural fire resistance tests, the interaction between the individual elements and the surrounding structure causes discrepancies in behaviour compared to single-element fire tests. Large-scale tests of real structures are challenging due to financial and time limitations. To bridge this gap, the concept of "Hybrid Fire Testing (HFT)" emerges, where a portion of the structural system (i.e., physical substructure) is experimentally tested while the remaining structure (i.e., numerical substructure) is analyzed numerically. The primary challenges in HFT involve ensuring stability throughout the analysis by considering the varying stiffness of the fire-exposed element during the test and establishing a versatile communication platform between the physical substructure (PS) and numerical substructure (NS) components. This paper presents a comprehensive HFT framework, implemented within a user-friendly software interface, facilitating both virtual and experimental testing. The software incorporates a new method addressing stability concerns by predicting PS stiffness during the test, achieving convergence within a limited number of iterations. Additionally, the framework includes a communication platform utilizing internet protocols (IP) and COM ports for rapid and easy connection to diverse experimental control systems and finite element software packages. The functionality of the developed software is validated through its successful application in an HFT conducted on a 3-story steel structure within a simulated environment. Both force-controlled and displacement-controlled approaches confirm the method's adaptivity to the employed test procedures. [ABSTRACT FROM AUTHOR]

The use of Reduced Dimensional Models (RDM) discretized like beams, plates and shell elements drastically decreases the computational cost of solving a full 3D elastic problem with a Finite Element Method (FEM). However, its kinematic assumptions are only applicable to bodies with regular sections or continuous layouts. For the correct analysis of irregular regions, it is necessary to rely on bi-dimensional or solid models that fully reproduce the geometry of the body and its behavior but have a much higher computational cost. The Mixing Dimensional Coupling (MDC) technique allows linking models discretized with elements of different topologies, allowing the possibility of considering the most cost-effective model in each region. This coupling takes place at the interface that delimits both models and relies on the equilibrium of work and reactions on its two faces. In this paper, the formulation is presented for coupling beams with laminar sections and 2D Plane-Stress (PS) models demonstrating its proper behavior. Finally, this coupling is used for defining a new beam element, the Beam-Like Reduced Order Model (BLROM), which is obtained from a Plane-Stress model of their longitudinal section. • The use of beam elements drastically reduces the computational cost of analogous full-dimensional models. • The MDC can be used to establish a weak coupling between a layered beam model and its plane-stress representation. • This technique can also be used to characterize beams with irregular cross section as ROM, compatible with beam elements. [ABSTRACT FROM AUTHOR]

Exploring the wave propagation in nano-waveguides enables the development of on-chip nano-resonators, leading innovative capabilities that enhance acoustics and micro-mechanics. This paper conducts a numerical investigation into the dynamic properties of periodic nano-waveguides. The Second Strain Gradient (SSG) elasticity, which captures the size effects, is used in the determination of the constitutive relations. The weak form including the element matrices is deduced by the Hamilton's principle. Furthermore, a finite element method, based on the periodic structure theory, is introduced to explore the free wave propagation by solving the eigenvalue problems. The stiffness hardening phenomenon manifested in the dispersion curves, band structures, and energy flow vector fields is discussed. The sensitivity of higher order parameters existing in the SSG theory is analyzed. The numerical investigation for wave propagation in the periodic nano-waveguides through the SSG theory is an innovative work, highlighting the significance of the proposed methodology in explaining the special dynamic characteristics of complex nano-waveguides. • Modeling of 3D nano-waveguides by second strain gradient theory. • Eigenvalue problems within 2D wave finite element method framework. • 2D wave motion properties in two different nano-waveguides. • Sensitivity of higher order parameters existing in second strain gradient elasticity. [ABSTRACT FROM AUTHOR]

A new approach is proposed to account for the impact of sprain variables for computational simulation of fracture in quasi-brittle materials. For an isotropic solid, these variables are equal to the Laplacians of the components of the displacement field, and they were introduced in a recent study by Zhang and Bazant to resolve spurious mesh effects for cases involving softening material response. This is possible by establishing localization-resisting forces which are conjugate to the displacements with respect to a sprain energy functional, defined as a function of the sprain variables. To enable the evaluation of the sprain variables in a conventional, C 0 continuum finite element approximation, this paper introduces a set of scalar fields, each being weakly equal to the Laplacian of one of the components of the displacement vector field. This leads to a straightforward implementation with standard continuum (solid) elements, eliminating the need for approximations satisfying higher continuity requirements or for ad-hoc, finite-difference algorithms which were adopted in the original pertinent work for the calculation of the sprain variables. • Use of localization-resisting forces in continuum simulations of strain softening, to eliminate spurious mesh-size effects. • These forces are governed by a set of scalar fields, weakly equivalent to the Laplacians of the displacement components. • Weak equivalence allows to obtain the Laplacians of the displacement field with conventional, C0 finite element approximations. • The method is found to eliminate spurious effects of mesh size and of mesh orientation with respect to crack propagation path. [ABSTRACT FROM AUTHOR]

• A new methodology for multi-parameter identification of concrete dams is proposed. • For the first time, polynomial chaos expansion and slime mould algorithm are integrated. • A very high computational efficiency is achieved using the proposed methodology. This paper presents a novel methodology that combines polynomial chaos expansion and slime mould algorithm for multi-parameter identification of concrete dams. This methodology not only incorporates the merits of low computational cost in the polynomial chaos expansion and fast convergence of slime mould algorithm, but also considers the priori uncertainty in the input parameters by introducing statistical probability theory. By considering two examples with different complexity, this paper verifies the effectiveness of the proposed method with a univariate simply supported beam model, followed by a complex multivariate dam model to demonstrate its practicability in real engineering problems. In addition, parameter sensitivity analysis of the dam model is conducted at an extremely low cost by polynomial chaos expansion based on Sobol' indices. Furthermore, the conventional parameter identification methods based on optimization methods directly combined with the finite element model are employed for comparison, highlighting two distinct advantages of the proposed method: (i) the proposed method improves the computational efficiency by nearly 52 times while ensuring a high accuracy, and (ii) the classical non-population optimization algorithm, Bayesian optimization, is used for comparison, revealing the outstanding performance of slime mould algorithm in terms of convergence speed and robustness. The application of the proposed algorithm is not only limited to dams, but also it can be extended to any structure. [ABSTRACT FROM AUTHOR]

• This paper presents an improved algorithm for the Moore-Penrose continuation method. • Its ability to compute solution curves with critical points is greatly improved. • Difficulties observed using the standard approach are not encountered. • Easy implementation. Using typical solution strategies to compute the solution curve of challenging problems often leads to the break down of the algorithm. To improve the solution process, numerical continuation methods have proved to be a very efficient tool. However, these methods can still lead to undesired results. In particular, near severe limit points and cusps, the solution process frequently encounters one of the following situations: divergence of the algorithm, a change in direction which makes the algorithm backtrack on a part of the solution curve that has already been obtained and omitting important regions of the solution curve by converging to a point that is much farther than the one anticipated. Detecting these situations is not an easy task when solving practical problems since the shape of the solution curve is not known in advance. This paper will therefore present a modified Moore-Penrose continuation method that will include two key aspects to solve challenging problems: detection of problematic regions during the solution process and additional steps to deal with them. The proposed approach can either be used as a basic continuation method or simply activated when difficulties occur. Numerical examples will be presented to show the efficiency of the new approach. [ABSTRACT FROM AUTHOR]

• Higher-order layer-wise plate models for composite variable angle tow plate analysis. • Passive damping using viscoelastic inter-layers with frequency-dependent properties. • Mixed interpolation of tensorial components for more efficient and accurate solutions. In this paper, the damped free-vibration and frequency response analysis of variable-angle-tow composite plates embedding viscoelastic layers with frequency-dependent properties are performed. The governing equations are derived from the Principle of Virtual Displacements and higher-order Layer-Wise models are used for the unknown variables description in the thickness direction, furthermore a nine-node finite plate element is employed to solve them. The plate formulation presented in this paper is more efficient with respect to three-dimensional approaches, typically used for this kind of structures, and it is more accurate than existing reduced order methods thanks to the use of the mixed interpolation of tensorial components technique, employed to solve the problem of the shear and membrane locking phenomena. Different multilayered structures have been considered, laminated with unidirectional cross-ply layers or with curvilinear fibers path ones. The use of viscoelastic soft sheets permits to have a passive damping of the structural vibration. The frequency dependent properties of the viscoelastic materials are described through the use of a Kelvin-Voigt model. Numerical solutions are presented for the analysis of variable-angle-tow composites with different boundary-conditions and various lamination schemes. [ABSTRACT FROM AUTHOR]

We present novel overlapping finite elements for a new paradigm of solution proposed in our previous papers, see Bathe (2016) and Bathe and Zhang (2017). We give the formulation of the new overlapping elements and the solutions of basic numerical examples to investigate the robustness and efficiency of the new finite elements. The results show that the new overlapping elements are quite distortion insensitive and the numerical integration of the element matrices is efficient. The computational effort to integrate the matrices is much less than in meshfree methods. Finally, we illustrate the complete solution scheme of the new paradigm using the overlapping finite elements in the analysis of the bracket problem already considered in Bathe and Zhang (2017). While the paper proposes and studies new overlapping elements, we conclude that further research is needed to fully harvest the potential of the new analysis approach. [ABSTRACT FROM AUTHOR]

Predicting the mechanical response of components requires simplifications and idealizations that affect the fidelity of the results and introduce errors. Some errors correspond to the limited knowledge of intrinsic physical attributes while others are introduced by the modeling framework and mathematical approximations. This paper studies the dependence of the force-displacement response of threaded fasteners on modeling attributes such as geometry, material, and friction resistance using finite element simulations. A systematic comparison of 1 D , 2.5 D or 3 D computational models demonstrates the influence of model properties and the limitations of the methodologies. Finally, the paper discusses the sources of model inputs and model form errors for threaded fasteners. [ABSTRACT FROM AUTHOR]

A major difficulty in finite element analysis is the preparation of an effective mesh leading to a good response solution. In engineering analyses of complex components, oftentimes, much more time is spent to arrive at an adequate mesh than to obtain the solution of the established finite element model. The difficulty in meshing is due to the fact that finite elements need to abut each other and cannot overlap. This can lead to highly distorted elements, e.g. sliver elements, reducing the accuracy of solution. In practice, to establish an effective mesh, frequently, significant expertise in building meshes is needed, because great care must be taken in cleaning up the CAD geometry and preparing an effective mesh. The objective in this paper is to present a new finite element solution scheme including meshing in which the elements can overlap. The property that finite elements can overlap removes many of the meshing difficulties, leads to an effective meshing procedure and an overall easy-to-use solution scheme for an analyst or a designer. We first present the meshing scheme that we propose, which combines the use of traditional finite elements and overlapping finite elements. A particular feature is that the meshing procedure can be directly embedded in CAD driven solutions. We then present the theory used for the formulation of the overlapping finite elements and the coupling with traditional finite elements. We consider spherical and brick-shaped overlapping finite elements for which the theory is largely based on the formulation of the method of finite spheres. Finally, we illustrate the complete solution scheme in the analysis of some two-dimensional problems using the CAD geometry as the starting point. While the paper presents a new paradigm for analysis in CAD environments, with much potential, we realize that much further study and research is needed on some of the important ingredients of the method to render the complete procedure effective for general practical engineering analyses of static and dynamic problems. [ABSTRACT FROM AUTHOR]

The aim of this study is an alternative approach to structure response or modal analysis. The structure consists of one-dimensional bars with continuously distributed mass and stiffness. The analysis is considered on an abstract basis as a problem of a differential system on an oriented graph. This graph is a geometric representation of the investigated mechanical system, where elements of the graph are individual bars of the system, recti- or curvilinear. The system as a whole is fixed through boundary conditions or interconnected with other sub-systems. Hence the paper can be taken as a follow up to earlier works presented at the CC2013 and CST2014 Conferences, where full mathematical background dealing with a general problem has been discussed. This paper is focused on the problem of dynamics of a system with straight prismatic bars with uniformly distributed mass. Dissipation of energy is omitted in order to keep the formulation in the real domain. The detailed assembly procedure of the Dynamic Stiffness Matrix (DSM) and transformation from local to global coordinates is outlined and demonstrated. Conventional way of eigenvalue searching by means of discrete alternative of the Newton-Raphson method is sketched out and later two possibilities based on polynomial and hyperbolic approximations of the DSM elements are pointed out. Lambda matrices as a tool are introduced together with a couple of application possibilities. The Wittrick-Williams algorithm is discussed and applied to localize and facilitate the eigenvalues searching on the whole frequency interval investigated. Finally, an illustrative example of the eigenvalue analysis of a structure is included. Strengths and shortcomings of the approach are discussed. Some open problems and orientation of further investigation are briefly outlined. The rapid development of analytical and semi-analytical methods being based at DSM idea is commented and justified in comparison with FEM analysis. Even if a certain overlap exists, the DSM approach represents an irreplaceable alternative of analysis especially in broadband dynamic problems stemming from area of the stochastic dynamics, like wind, earthquake or traffic engineering. Problems of dynamic stability are easier to be managed as well, when using DSM oriented approaches together with related algorithms of the analytical and semi-analytical method family. Moreover, many problems of physics at the nano-scale emerged recently where these procedures are preferable to be used. [ABSTRACT FROM AUTHOR]

• Strength-based structural/material layout optimization for solid reinforced concrete. • The design cost is minimized based on the individual price of concrete and reinforcement. • The optimization is strength-based, i.e., strength scaling parameters on the concrete and reinforcement ratios are used. • The material layout is defined based on the concrete strength scaling parameters and reinforcement ratios. Whether monetary, environmental, or otherwise, the cost of the structure is always an important aspect when designing reinforced concrete. A structure with a low cost is desirable, but other considerations must be made. First and foremost, the structures should be safe. The two conditions define the structural optimization problem that this paper seeks to solve as finding the structure with the lowest cost, which is also safe. Validation of reinforced concrete in the ultimate limit state is often performed using plasticity and strength-based methods. This paper combines Finite Element Limit Analysis (FELA) with a yield criterion that includes material layout variables to create a strength-based structural optimization framework for solid reinforced concrete structures. The yield criterion combines the Modified Mohr–Coulomb criterion, where variables scale the compressive strength, with the inclusion of smeared reinforcement, where variables scale the reinforcement ratios. The framework is presented and shown in examples using free optimization and examples where material layout variables in parts of the structure are coupled and optimized as a group. [ABSTRACT FROM AUTHOR]

Tangled (non-convex) elements, i.e. elements with negative Jacobian determinant, can lead to erroneous results in the standard finite element method (FEM). Constructing tangle-free, well-structured meshes for complex geometries is often impossible. Hence there is a need to explore analysis methods that can directly handle such tangled meshes. In this paper, we propose the isoparametric tangled finite element method (i-TFEM) for free and forced vibration problems over tangled meshes. By employing piece-wise invertible mapping, a variational formulation is derived, leading to a simple modification of the standard FEM stiffness and mass matrices with the incorporation of additional compatibility constraints. Moreover, i-TFEM reduces to standard FEM for non-tangled (regular) meshes. The proposed method is implemented for three types of elements: 4-node quadrilateral, 9-node quadrilateral, and 8-node hexahedral elements. The numerical results demonstrate that i-TFEM is able to consistently handle general tangled (non-convex) elements, enabling convenient meshing for complex geometries. • Isoparametric tangled finite element method (i-TFEM) is proposed for free and forced vibration problems over tangled meshes. • To handle concave elements using i-TFEM, simple modifications to the standard FEM stiffness and mass matrices are needed. • I-TFEM is implemented for three element types: 4-node quadrilateral, 9-node quadrilateral, and 8-node hexahedral elements. • The numerical results demonstrate that i-TFEM consistently handles general non-convex elements. [ABSTRACT FROM AUTHOR]

• FEM analyzes of two representative historic railroad arch masonry bridges are shown. • Careful selection of the masonry properties is crucial and always required. • Modern traffic loads may contribute to the damage of masonry railroad bridges. • Cohesive zone modeling provides efficient description of masonry cracking. • Broad analysis of numerous masonry bridges require effective approaches as shown. The problems encountered during the analyzes of structural response of historic masonry railroad arch bridges are described in this paper. The attention is mainly focused on the stiffness of the masonry arches, their strengths and appropriate estimation of railroad load intensity. Issues related to computational modelling of two, existing, almost 130 years old masonry arch railroad bridges are presented in this context. The main properties and results of detailed inspection of the structures are shown. Computational models that were created in the finite element method environment in order to efficiently describe the responses of the bridges under typical loading conditions and estimate their load carrying abilities are presented. The outcomes of several nonlinear static analyses that were conducted for this purpose are discussed. What is more the results of finite element analyses are reviewed against the inspected bridge condition and final conclusions are formulated on that basis. All the analyses allowed to find the possible causes of the deterioration of the bridges condition. [ABSTRACT FROM AUTHOR]

This paper presents a numerical model of coupled hydro-mechanical behaviour of partially saturated soils. The model is based on Bishop's effective stress concept and the Cam-Clay model, with the compression stiffness influenced by the effect of partial saturation. The mechanical model is accompanied by the simplified water transport model where the water pressure is the only variable and van Genuchten moisture retention curve was adopted. The governing equations are discretized in the framework of the finite element method, and the resulting model is implemented in the in-house open source software package SIFEL. To improve convergence of the Newton-Raphson method used at the global level of the equilibrium finding, the algorithmic stiffness matrix was derived for the used stress return algorithm. A numerical examples of foundation problem illustrates the effective application of the model and proved that the model can capture wetting collapse. • Coupled hydro-mechanical model for partially saturated soils. • Improvement of the model performance by the algorithmic stiffness matrix for the cutting plane stress algorithm. • The effect of the wetting collapse can be captured by the model, as demonstrated in the example. • Comparison of model performance on the partially and fully saturated cases of the footing strip problems. [ABSTRACT FROM AUTHOR]

This paper introduces a novel hybrid formulation of a peridynamics-based machine learning model for thermal diffusion analysis in one-dimensional and two-dimensional problems with evolving discontinuities. The hybrid model employs a multivariate linear regression approach to establish the relationship between the temperature values of material points, their neighboring points, and the applied external heat flux. The thermal modal analysis method uses the finite element method to generate training and testing data. An efficient numerical procedure is also developed to couple the peridynamics model and the peridynamics-based machine learning model. The model is analyzed under multiple configurations of micro-thermal conductivity functions for one-dimensional thermal bar problems, under both steady-state and transient loading conditions, to identify the configuration that exhibits superior convergence behavior towards the local solution. Furthermore, benchmark problems demonstrate the high accuracy of the formulated hybrid model, including the analysis of thermal plates with discontinuities, such as a plate with a hole, a plate with a pre-existing insulated crack, and plate with a bi-material interface with stationary and evolving discontinuances. The hybrid model effectively captures intricate discontinuities and boundaries while being computationally efficient, indicating its potential for thermal diffusion analysis in one- and two-dimensional problems with stationary and evolving discontinuities. • The peridynamics-based machine learning model is proposed for diffusion problems. • An efficient scheme is proposed to couple peridynamics and machine learning models. • Machine learning data is generated using modal analysis for different thermal conditions. • The efficiency and robustness of the model is demonstrated through numerical examples. • The model is straightforward to extend to diffusion problems associated with pitting corrosion. [ABSTRACT FROM AUTHOR]

• The use of the Embedded Finite Element Method (E-FEM) to simulate local material heterogeneities is analyzed in detail. • A new fully consistent weak discontinuity enhancement function is proposed. • Basic simulations are used to compare the classical approach, commonly used in literature, and the new fully consistent approach proposed in the paper. • A discussion evaluates both model performance and ease of implementation. This paper analyses in detail the use of the Embedded Finite Element Method (E-FEM) to simulate local material heterogeneities. The work starts by a short review on the evolution of weak discontinuity models within the E-FEM framework to discuss how they account for the presence of multiple materials within a single element structure. A theoretical basis is introduced through some mathematical weak discontinuity definitions and the Hu-Washizu variational principle, for then establishing a set of requirements for retaining variational and kinematic consistency for any weak discontinuity enhancement proposal. From a general definition of a displacement enhancement field, two particular enhancement functions are derived by considering different consistency requirements: one which has been typically used in previous works and other which truly possesses variational consistency. A discussion is held on enhancement stability properties and the impact to global finite element solution processes. In the end, numerical simulations are made to assess the performance of each of these enhancements on the task of modelling a classical bi-material layered 3D tension problem and a more realistic heterogeneous sample having spherical inclusions of different radii. The final discussion evaluates both model performance and ease of implementation. [ABSTRACT FROM AUTHOR]

In the paper by Ham et al. (2014), the method of finite spheres enriched for transient wave propagation problems was used and a monotonic convergence of the calculated solutions with decreasing time step size was seen. This is an important property for practical analyses and different from what is seen in traditional finite element solutions. In this paper we explicitly show and study this characteristic through a dispersion analysis of the solutions calculated by the method of finite spheres using an implicit time integration method, the Bathe method. Another important property identified is that in uniform spatial discretizations for the problems solved, the calculated solution accuracy is almost independent of the solution direction considered. Numerical solutions of some wave propagation problems are given to demonstrate these attributes and show that with the schemes discussed very accurate response predictions can be obtained. [ABSTRACT FROM AUTHOR]

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton–Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results. [ABSTRACT FROM AUTHOR]

• The paper proposed a physics-informed DNN framework for forecasting the hysteretic curves of S-shaped dampers. • The proposed physics-informed DNNs were calibrated numerically and experimentally. • The physics-informed DNNs can be used for the seismic design and mechanical characteristics assessment of S-shaped dampers Numerical simulation that combines finite element methods and experimental data has been recognized as effective in modeling hysteretic behaviors and capturing the principle mechanical trend of passive energy dissipation devices. However, the seismic design and mechanical characteristics assessment for passive energy dissipation devices require a laborious effort and massive computational resources to tune mechanical-oriented parameters. Particularly, the potential risk of departure from the actual system is rising for the desirable seismic design of dampers due to the numerical model simplification and assumption. To eliminate the potential weakness of numerical models, this paper explores a surrogate model by implementing physics-informed deep neural networks (DNNs) to approximate hysteretic behaviors of S-shaped steel dampers. The proposed physics-informed DNNs mainly consists of recurrent neural networks (RNNs) and long short-term networks (LSTMs), which can encode the Bouc-Wen model into the direct graph and incorporate the effect of design-oriented geometry parameters. To validate the generality of the network, the optimization model was calibrated numerically and experimentally, respectively, which exhibits good performance in predicting nonlinear behaviors of different dampers with reasonable accuracy. The proposed physics-informed DNNs can be an alternative to relieve the laboriousness of the seismic design and mechanical characteristics assessment of passive energy dissipation devices. [ABSTRACT FROM AUTHOR]

In this paper, a solid finite element formulation applied to the analysis of functionally graded materials (FGMs) under finite elastoplastic deformation with mixed hardening is presented. The main novelty of this paper is the use of gradually variable material coefficients in finite elastoplastic strain regime. The numerical simulations are performed with the same constitutive models but with different variations of the material coefficients. Full integration and high order tetrahedral elements are used. Results show that high order elements and mesh refinement avoids general locking problems. Finally, the differences (regarding final results) between the homogeneous and the FG cases are depicted. [ABSTRACT FROM AUTHOR]

• CS-FEM is applied to the numerical simulation of viscoelastic fluid flow for the first time. • The gradient smoothing for viscoelastic fluid flow is arbitrary within each smoothing cell. • Two smoothed CBS schemes are presented to solve the viscoelastic equations. • The cell-based smoothing procedure is used to evaluate the viscoelastic fluid drag. This paper generalizes the cell-based smoothed finite element method (CS-FEM) to viscoelastic fluid flows under the fractional-step umbrella. The characteristic-based split (CBS) scheme is employed to decouple the Navier–Stokes equations and the Oldroyd-B constitutive model. Moreover, the discrete elastic-viscous-split-stress (DEVSS) procedure and stabilized pressure gradient projection (SPGP) are introduced into a second-order CBS framework. The projections of velocity and pressure gradients therein combine to form the DEVSS-G/CBS(B)-SPGP algorithm that further stabilizes the discretized model. Both proposed schemes are suitable for use in CS-FEM because they accept equal low-order interpolations for primitive variables. The spatial discretization is based upon bilinear four-node quadrilateral elements where all gradient-related items are smoothed without any methodological difficulties. In particular, the cell-based smoothing concept enables accurate evaluation of viscoelastic fluid force. The developed methods are validated against previously published data for several benchmark problems. A satisfactory agreement exists between the present and earlier results. [ABSTRACT FROM AUTHOR]

• A new sectional analysis tool for concrete member torque-twist response is presented. • Elastic shear strain distribution assumed for post-cracked torsional behaviour. • Framework is flexible, permitting any constitutive model for concrete or steel. • The approach is validated against a large set of experiments in the literature. • Detailed results obtained with lower computational costs then nonlinear finite element analysis. This paper presents a sectional analysis tool which can compute the complete torque-twist behaviour of reinforced concrete beams and shells. The cross section is modelled as a 2-D grid of triaxial elements representing the concrete and uniaxial elements representing the longitudinal reinforcement. Fixed strain patterns based on kinematic assumptions are used instead of using the finite element method to calculate strain distributions corresponding to sectional strains. Constitutive models for the concrete are based on the Modified Compression-Field Theory, though tension stiffening is neglected in the proposed model's current implementation. A key assumption used in this model is that the shear strain distribution caused by twists obtained from a linear elastic analysis can be used to describe the nonlinear behaviour of a member following cracking. The validity of this assumption when applied to rectangular sections is confirmed based on a study of 115 tests from the literature. Excellent agreement of peak strength and overall torque-twist behaviour are observed when comparing the model predictions and experimental data from these tests. Areas of future work are to improve the capabilities of the model are identified. [ABSTRACT FROM AUTHOR]

• Two new functionals are obtained for geometrical nonlinear analysis of beams. • Mixed finite element method is used with boundary terms effectively. • The structure of incremental formulation is given to solve nonlinear equations. • Comparison of mixed formulations for the planar curved beams is presented. In the present paper, two mixed type functionals are presented to analyze the geometrical nonlinear bending behavior of slender planar curved beam by Finite Element Method (FEM). Bernoulli-Euler beam theory is taken into consideration within the context of the small strains but large displacements. Non-linear bending of curved beams in-plane has been investigated in both local and Cartesian coordinate systems. The Euler-Lagrange equations of the Functional I and Functional II which give the equilibrium and constitutive equations are in relation with the Green-Lagrange strain tensor and the first-order approximation (consistent linearization) of the Green-Lagrange strain tensor, respectively. The details of this process are given in Appendix. In the formulation of the problem in Cartesian coordinate, in addition to the boundary terms, that expression has been used, also to obtain the strain quantities which should be in the constitutive expressions in-prior to the equilibrium equations for the first time in this work. For the solution of the non-linear equations of the problems "incremental formulation" has been used. The validity, reliability and robustness of the proposed functionals have been presented by solving sample problems found in the literature. [ABSTRACT FROM AUTHOR]

• A multi-scale beam model is derived using FE2 framework and absolute nodal coordinate formulation. • It has both the efficiency of absolute nodal coordinate formulation and multi-scale framework. • Performance of the proposed method is demonstrated through numerical examples. In this paper, a multi-scale modeling approach for the dynamic response of flexible multibody system made of composite material is suggested using the absolute nodal coordinate formulation (ANCF) – multilevel finite element (FE2) method. On many occasions, the analysis of composite structures requires a multi-scale modeling approach. FE2 method is one of the famous multi-scale modeling methods. And it is a versatile computational homogenization approach that is commonly applicable to many different materials. Because the present FE2 is capable of computing the mechanical response at two different scales simultaneously, ANCF element will be used at the macroscopic scale. And the microscopic RVE at each macroscopic integration point will be obtained by the finite element method. Since ANCF element uses only the absolute nodes and gradients, the mass matrix will be constant. Due to such characteristics, when combined with FE2, the linearization process will become easier compared against the existing formulation. And Coriolis or centrifugal force will not need to be considered in the dynamic formulation. Therefore, computational cost will be reduced, and the relevant algebraic manipulation will become simplified. Several numerical examples are presented to demonstrate the improved accuracy and diminished computational cost of the present suggestion. [ABSTRACT FROM AUTHOR]

• Propose a smart component model replacement (SCMR) approach to simulate large-scale nonlinear structure. • The SCMR approach can finely simulate the local damage and failure mechanism in structures. • Present "replacement" strategies for smoothly replacing between different models. • The efficiency and accuracy of SCMR approach are verified by three application examples. This paper proposes a novel smart component model replacement (SCMR) approach that is not only efficient for large-scale nonlinear structural analysis, but also accurate in the refined simulation of components of interest. The method creates a FE model with a relatively coarse mesh to simulate the entire structure. During analysis, when a component is identified to predefined criteria defined by the users, its model will be replaced by another more refined "substructure" model that can simulate the damage and failure more accurately. A few "replacement" strategies are presented for smoothly transiting between different models: (1) A gradually transiting strategy is presented to avoid the possible jumps of responses due to the model differences to guarantee numerical stability; (2) The historical boundary conditions are applied to the substructure model to address the historical dependence of responses; (3) The boundary connection between main and substructure is considered by using a flexible client/server technology, to simplify the replacement. The SCMR approach is verified by three application examples. The results demonstrate that the proposed modeling approach can efficiently perform large-scale nonlinear structural analysis, while capturing the failure mechanism in RC structures well. [ABSTRACT FROM AUTHOR]

This paper presents a novel derivation of the direct parametrisation method for invariant manifolds able to build simulation-free reduced-order models for nonlinear piezoelectric structures, with a particular emphasis on applications to Micro Electro Mechanical Systems. The constitutive model adopted accounts for the hysteretic and electrostrictive response of the piezoelectric material by resorting to the Landau-Devonshire theory of ferroelectrics. Results are validated with full-order simulations operated with a harmonic balance finite element method to highlight the reliability of the proposed reduction procedure. Numerical results show a remarkable gain in terms of computing time as a result of the dimensionality reduction process over low dimensional invariant sets. Results are also compared with experimental data to highlight the remarkable benefits of the proposed model order reduction technique. • Model order reduction for structures with geometric and piezoelectric nonlinearity. • Direct computation of the invariant manifold from finite element problem. • Efficient technique accounting for electrostrictive response of the piezoelectric. • High accuracy of the results compared to full-order simulation. • Application to micro-electro mechanical systems with complex geometry. [ABSTRACT FROM AUTHOR]

• A modified stain energy release rate method is proposed to enable the method for solid element. • An algorithm is proposed to map the crack information obtained by vision method in the FE model. • The effectiveness and generality of the proposed framework are validated using three cracking FE configurations. Vision inspection information of damages in civil engineering structures is hard to be transformed into an analysis-available resistance model to conduct in-service condition assessment. This paper proposes an updating crack-embedded finite element modeling framework using periodical vision inspection information, which offers a convenient crack embedding approach for safety assessment of cracked reinforced concrete (RC) frame structures. For this purpose, a modified stain energy release rate (SERR) method is proposed to enable the method for solid element with high accuracy. Consequently, an algorithm is proposed to map the crack information obtained by vision method in the FE model. The combination of modified SERR method and proposed mapping algorithm could achieve flexibility matrix updating of cracked element. In addition, the effectiveness and generality of the proposed framework are validated using three cracking configurations. This research can provide a feasible framework for the rapid condition assessment of in-service civil engineering structures, showing a great prospect of intelligent disaster prevention. [ABSTRACT FROM AUTHOR]

This publication presents a novel numerical avenue for the assessment of masonry structures. To circumvent the labor-intensive finite element method modeling and pre-processing step, the isogeometric analysis is enhanced within this contribution to cope with the necessities for masonry structures. This incorporates the study of apparent kinematic shell formulations, including their respective properties. Additionally, the connection towards non-linear material laws is outlined, and an empirically evaluated approach for the element size regularization is proposed. This generic approach would be inherently valid for other size dependent constitutive material models, as e.g. concrete models. The current paper is explored in the context of a number of experimental problems. These validation examples are used to provide structural assessments that extend beyond the scope of the experiment. Finally, a potential example that provides an outlook towards anticipated applications is presented. • Advancing the isogeometric analysis to cope with shell-based masonry structures. • Proposing a damage model in the context of the isogeometric analysis. • Newly proposed element length regularization technique for NURBS. • Enhancement of parametric and CAD-integrated workflow for masonry analysis. • Comprehensive mechanical assessments of masonry structures. [ABSTRACT FROM AUTHOR]

• A novel stochastic homotopy method is proposed to solve static problems involving arbitrarily distributed random variables. • The optimal homotopy series expansion is obtained by minimizing the stochastic residual error. • The proposed method shows better convergency and efficiency than the aPC method in the case of non-Gaussian random field. • This method is successfully applied to the stochastic static analysis of the large-scale engineering structure. The modelling of realistic engineering structures with uncertainties often involves various probabilistic distribution types, which bring forward higher requirements for the generality of stochastic analysis methods. This paper proposes a novel stochastic homotopy method to evaluate the static response of the structure with random variables of arbitrary distributions. In this method, a homotopy series expansion is used to approximate the stochastic static response, and the optimal form of the expansion can be determined by minimizing the residual error about the stochastic static equilibrium equation regardless of distribution types of random variables. The numerical results of a logarithmic function and a thin plate show that the new method exhibits excellent accuracy and stability compared to the homotopy stochastic finite element method depending on the sample selection. Compared with the arbitrary polynomial chaos method (aPC), the proposed method is more efficient under equivalent accuracy. On the other hand, as the expansion order increases, this new method shows better convergence than the aPC method and the perturbation stochastic finite element method in the case of non-Gaussian distributed random variables of large fluctuation. In addition, a cable-stayed bridge example illustrates the application of the proposed method on the large-scale structure. [ABSTRACT FROM AUTHOR]

• The wave input method has a deficiency in expressing complex incident site waves. • A random non-uniform wave input method is proposed for seismic response analyses. • The new method is expressible for complex incident waves and motion differences. To improve the deficiency of the wave input method (WIM) in expressing complex incident site waves, in this paper, a random non-uniform wave input method (RNUWIM) is derived from the use of the random incident waves and the visco-elastic boundary condition. The random incident waves are generated and adjusted by the comparison of the horizontal ground motion's response spectrum and the target spectrum, so that they are representative of the actual wavefield in project sites. Furthermore, the RNUWIM is implemented and verified by the finite element method. The result shows that the RNUWIM inherits the advantages of the WIM in simulating the wave propagation effect and the radiation damping effect. In the new method, the spatial variation of ground motions such as time delays, peak variation, and time-history shape differences can be simulated synthetically. The errors of RNUWIM in ground acceleration, velocity, and displacement at the reference point are respectively 3.23%, 0.66%, and 0.07% in one case, which proves the efficiency of the RNUWIM in seismic response analyses. [ABSTRACT FROM AUTHOR]