Search

Purpose: Hexahedral meshing is one of the most important steps in performing an accurate simulation using the finite element analysis (FEA). However, the current hexahedral meshing method in the industry is a nonautomatic and inefficient method, i.e. manually decomposing the model into suitable blocks and obtaining the hexahedral mesh from these blocks by mapping or sweeping algorithms. The purpose of this paper is to propose an almost automatic decomposition algorithm based on the 3D frame field and model features to replace the traditional time-consuming and laborious manual decomposition method. Design/methodology/approach: The proposed algorithm is based on the 3D frame field and features, where features are used to construct feature-cutting surfaces and the 3D frame field is used to construct singular-cutting surfaces. The feature-cutting surfaces constructed from concave features first reduce the complexity of the model and decompose it into some coarse blocks. Then, an improved 3D frame field algorithm is performed on these coarse blocks to extract the singular structure and construct singular-cutting surfaces to further decompose the coarse blocks. In most modeling examples, the proposed algorithm uses both types of cutting surfaces to decompose models fully automatically. In a few examples with special requirements for hexahedral meshes, the algorithm requires manual input of some user-defined cutting surfaces and constructs different singular-cutting surfaces to ensure the effectiveness of the decomposition. Findings: Benefiting from the feature decomposition and the 3D frame field algorithm, the output blocks of the proposed algorithm have no inner singular structure and are suitable for the mapping or sweeping algorithm. The introduction of internal constraints makes 3D frame field generation more robust in this paper, and it can automatically correct some invalid 3–5 singular structures. In a few examples with special requirements, the proposed algorithm successfully generates valid blocks even though the singular structure of the model is modified by user-defined cutting surfaces. Originality/value: The proposed algorithm takes the advantage of feature decomposition and the 3D frame field to generate suitable blocks for a mapping or sweeping algorithm, which saves a lot of simulation time and requires less experience. The user-defined cutting surfaces enable the creation of special hexahedral meshes, which was difficult with previous algorithms. An improved 3D frame field generation method is proposed to correct some invalid singular structures and improve the robustness of the previous methods. [ABSTRACT FROM AUTHOR]

Purpose: Pores and shrink holes are unavoidable defects in the die-casting mass production process which may significantly influence the strength, fatigue and fracture behaviour as well as the life span of structures, especially if they are subjected to high static and dynamic loads. Such defects should be considered during the design process or after production, where the defects could be detected with the help of computed tomography (CT) measurements. However, this is usually not done in today's mass production environments. This paper deals with the stress analysis of die-cast structural parts with pores found from CT measurements or that are artificially placed within a structure. Design/methodology/approach: In this paper the authors illustrate two general methodologies to take into account the porosity of die-cast components in the stress analysis. The detailed geometry of a die-cast part including all discontinuities such as pores and shrink holes can be included via STL data provided by CT measurements. The first approach is a combination of the finite element method (FEM) and the finite cell method (FCM), which extends the FEM if the real geometry cuts finite elements. The FCM is only applied in regions with pores. This procedure has the advantage that all simulations with different pore distributions, real or artificial, can be calculated without changing the base finite element mesh. The second approach includes the pore information as STL data into the original CAD model and creates a new adapted finite element mesh for the simulation. Both methods are compared and evaluated for an industrial problem. Findings: The STL data of defects which the authors received from CT measurements could not be directly applied without repairing them. Therefore, for FEM applications an appropriate repair procedure is proposed. The first approach, which combines the FEM with the FCM, the authors have realized within the commercial software tool Abaqus. This combination performs well, which is demonstrated for test examples, and is also applied for a complex industrial project. The developed in-house code still has some limitations which restrict broader application in industry. The second pure FEM-based approach works well without limitations but requires increasing computational effort if many different pore distributions are to be investigated. Originality/value: A new simulation approach which combines the FEM with the FCM has been developed and implemented into the commercial Abaqus FEM software. This approach the authors have applied to simulate a real engineering die-cast structure with pores. This approach could become a preferred way to consider pores in practical applications, where the porosity can be derived either from CT measurements or are artificially adopted for design purposes. The authors have also shown how pores can be considered in the standard FEM analysis as well. [ABSTRACT FROM AUTHOR]

Purpose: The purpose is to predict the distribution of the residual pretightening force of the bolt group under the action of any initial pretightening force, and to achieve the final residual pretightening force as the target to solve the initial pretightening force value to be applied. Design/methodology/approach: Based on the finite element method and the elastic interaction theory between bolt group, this paper establishes a prediction model for the residual pretightening force distribution of bolt group for one-step pretightening and multi-step pretightening of gasketless flange connection systems. In addition, using the general modeling method given in this paper, the prediction model of residual pretightening force of long plate bolt connection system is established, and compared with reference, which fully proves the effectiveness and universality of the general prediction model of residual pretightening force of bolt group. Findings: The appropriate pretightening sequence, increasing the number of pretightening steps and variable amplitude loading can effectively reduce the influence of elastic interaction and improve the uniformity of residual pretightening force of the bolt group. And the selection of material, number of bolts and connected thickness of bolt connection system also has a great influence on the distribution of residual pretightening force of bolt groups. Originality/value: The general prediction model for the residual pretightening force of bolt group of connecting structural components considering elastic interaction given in this paper can provide a reference for the design and optimization of the bolt assembly process of the rotor system and the casing system in aero-engine and the prediction of the performance of the connecting system. [ABSTRACT FROM AUTHOR]

Purpose: Realistic composite vehicles with 2, 3, 5 and 9 axles, consisting of a truck with one or two trailers, are addressed in this paper by computational models for vehicle–bridge interaction analysis. Design/methodology/approach: The vehicle–bridge interaction (VBI) models are formed by sets of 2-D rigid blocks interconnected by mass, damping and stiffness elements to simulate their suspension system. The passage of the vehicles is performed at different speeds. Several rolling surface profiles are admitted, considering the maintenance grade of the pavement. The spectral density functions are generated from an experimental database to form the longitudinal surface irregularity profiles. A computational code written in Phyton based on the finite element method was developed considering the Euler–Bernoulli beam model. Findings: Several models of composite heavy vehicles are presented as manufactured and currently travel on major roads. Dynamic amplification factors are presented for each type of composite vehicle. Research limitations/implications: The VBI models for compound heavy vehicles are 2-D. Social implications: This work contributes to improving the safety and lifetime of the bridges, as well as the stability and comfort of the vehicles when passing over a bridge. Originality/value: The structural response of the bridge is affected by the type and size of the compound vehicles, their speed and the conservative grade of the pavement. Moreover, one axle produces vibrations that can be superposed by the vibrations of the other axles. This effect can generate not usual dynamic responses. [ABSTRACT FROM AUTHOR]

Purpose: In this paper, based on the density functional theory (DFT) and finite element method (FEM), the elastic, buckling and vibrational behaviors of the monolayer bismuthene are studied. Design/methodology/approach: The computed elastic properties based on DFT are used to develop a finite element (FE) model for the monolayer bismuthene in which the Bi-Bi bonds are simulated by beam elements. Furthermore, mass elements are used to model the Bi atoms. The developed FE model is used to compute Young's modulus of monolayer bismuthene. The model is then used to evaluate the buckling force and fundamental natural frequency of the monolayer bismuthene with different geometrical parameters. Findings: Comparing the results of the FEM and DFT, it is shown that the proposed model can predict Young's modulus of the monolayer bismuthene with an acceptable accuracy. It is also shown that the influence of the vertical side length on the fundamental natural frequency of the monolayer bismuthene is not significant. However, vibrational characteristics of the bismuthene are significantly affected by the horizontal side length. Originality/value: DFT and FEM are used to study the elastic, vibrational and buckling properties of the monolayer bismuthene. The developed model can be used to predict Young's modulus of the monolayer bismuthene accurately. Effect of the vertical side length on the fundamental natural frequency is negligible. However, vibrational characteristics are significantly affected by the horizontal side length. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle. Design/methodology/approach: The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods. Findings: The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements. Originality/value: This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures. [ABSTRACT FROM AUTHOR]

Purpose: Natural gas hydrate (NGH) has been regarded as one of the most important resources due to NGH's large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Therefore, developing a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior is very important and necessary. The purpose of this article is to develop and solve a multi-phase, strong nonlinearity and large-scale fully coupled thermal-hydro-mechanical–chemical (THMC) model for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition. Design/methodology/approach: In this paper, a multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH dissociation. The fully coupled THMC model is solved by using a fully implicit finite element method, in which the gas pressure, water pressure, temperature and displacement are taken as basic unknown variables. The proposed model is validated against with the experimental data, showing high accuracy and reliability. Findings: A multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition. The proposed model is validated against with the experimental data, showing high accuracy and reliability. Research limitations/implications: Some assumptions are made to make the model tractable, including (1) the composition gas of hydrate is pure methane; (2) the gas-liquid multi-phase flow in the pore obeys Darcy's law; (3) hydrate occurs on the surface of soil particles, both of them form the composite consolidation material; (4) the small-strain assumption is applied to composite solid materials, which are treated as skeletons and cannot be moved; (5) momentum change caused by phase change is not considered. Practical implications: NGH has been regarded as one of the most important resources due to its large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Most of the existing studies decouple the process with solid deformation and seepage behavior, but the accuracy of the numerical results will be sacrificed to certain extent. Therefore, it is very important and necessary to develop a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior. Social implications: NGH, widely distributed in shallow seabed or permanent frozen region, has the characteristics of high energy density and high combustion efficiency (Yan et al., 2020). A total of around 7.5 × 1,018 m3 has been proved to exist around the world and 1 m3 of NGH can release about 160–180 m3 of natural gas (Kvenvolden and Lorenson) under normal conditions. Safely and sustainably extracting NGH commercially can effectively relieve global energy pressure and contribute to achieving carbon reduction goals. Originality/value: The novelty of the present work lies in mainly two aspects. First, a fully coupled THMC model is developed for studying the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and solid deformation during NGH dissociation. Second, the numerical solution is obtained by using a fully implicit finite element method (FEM) and is validated against experimental data. [ABSTRACT FROM AUTHOR]

Purpose: In this paper, a superconvergent patch recovery method is proposed for superconvergent solutions of modes in the finite element post-processing stage of variable geometrical Timoshenko beams. The proposed superconvergent patch recovery method improves the solution speed and accuracy of the finite element analysis of a curved beam. The free vibration and natural frequency of the beam were considered for studying forced vibrations and structural resonance. Beam vibration mode analysis was performed for high-precision vibration mode solutions and frequency values. The proposed method can be used to compute beam vibration modes of beams with different shapes and boundary conditions as well as variable cross sections and curvatures. The purpose of this paper is to address these issues. Design/methodology/approach: An adaptive method was proposed to analyse the in-plane and out-of-plane free vibrations of the variable geometrical Timoshenko beams. In the post-processing stage of the displacement-based finite element method, the superconvergent patch recovery method and high-order shape function interpolation technique were used to obtain the superconvergent solution of mode (displacement). The superconvergent solution of mode was used to estimate the error of the finite element solution of mode in the energy form under the current mesh. Furthermore, an adaptive mesh refinement was proposed by mesh subdivision to derive an optimised mesh and accurate finite element solution to meet the preset error tolerance. Findings: The results computed using the proposed algorithm were in good agreement with those computed using other high-precision algorithms, thus validating the accuracy of the proposed algorithm for beam analysis. The numerical analysis of parabolic curved beams, beams with variable cross sections and curvatures, elliptically curved beams and circularly curved beams helped verify that the solutions of frequencies were consistent with the results obtained using other specially developed methods. The proposed method is well suited for the mesh refinement analysis of a curved beam structure for analysing the changes in high-order vibration mode. The parts where the vibration mode changed significantly were locally densified; a relatively fine mesh division was adopted that validated the reliability of the mesh optimisation processing of the proposed algorithm. Originality/value: The proposed algorithm can obtain high-precision vibration solutions of variable geometrical Timoshenko beams based on more optimized and reasonable meshes than the conventional finite element method. Furthermore, it can be used for vibration problems of parabolic curved beams, beams with variable cross sections and curvatures, elliptically curved beams and circularly curved beams. The proposed algorithm can be extended for application in superconvergent computation and adaptive analysis of finite element solutions of general structures and solid deformation fields and used for adaptive analysis of more complex plates, shells and three-dimensional structures. Additionally, this method can analyse the vibration and stability of curved members with crack damage to obtain high-precision vibration modes and instability modes under damage defects. [ABSTRACT FROM AUTHOR]

Purpose: The present paper aims to explore a computational homogenisation procedure to investigate the full geometric representation of yield surfaces for isotropic porous ductile media. The effects of cell morphology and imposed boundary conditions are assessed. The sensitivity of the yield surfaces to the Lode angle is also investigated in detail. Design/methodology/approach: The microscale of the material is modelled by the concept of Representative Volume Element (RVE) or unit cell, which is numerically simulated through three-dimensional finite element analyses. Numerous loading conditions are considered to create complete yield surfaces encompassing high, intermediate and low triaxialities. The influence of cell morphology on the yield surfaces is assessed considering a spherical cell with spherical void and a cubic RVE with spherical void, both under uniform strain boundary condition. The use of spherical cell is interesting as preferential directions in the effective behaviour are avoided. The periodic boundary condition, which favours strain localization, is imposed on the cubic RVE to compare the results. Small strains are assumed and the cell matrix is considered as a perfect elasto-plastic material following the von Mises yield criterion. Findings: Different morphologies for the cell imply in different yield conditions for the same load situations. The yield surfaces in correspondence to periodic boundary condition show significant differences compared to those obtained by imposing uniform strain boundary condition. The stress Lode angle has a strong influence on the geometry of the yield surfaces considering low and intermediate triaxialities. Originality/value: The exhaustive computational study of the effects of cell morphologies and imposed boundary conditions fills a gap in the full representation of the flow surfaces. The homogenisation-based strategy allows us to further investigate the influence of the Lode angle on the yield surfaces. [ABSTRACT FROM AUTHOR]

Purpose: The objective of this paper is to establish a solution strategy for obtaining dual solutions, namely trivial (conventional) and nontrivial (unconventional) solutions, of coupled pore-fluid flow and chemical dissolution problems in heterogeneous porous media. Design/methodology/approach: Through applying a perturbation to the pore-fluid velocity, original governing partial differential equations of a coupled pore-fluid flow and chemical dissolution problem in heterogeneous porous media are transformed into perturbed ones, which are then solved by using the semi-analytical finite element method. Through switching off and on the applied perturbation terms in the resulting perturbed governing partial differential equations, both the trivial and nontrivial solutions can be obtained for the original governing partial differential equations of the coupled pore-fluid flow and chemical dissolution problem in fluid-saturated heterogeneous porous media. Findings: When a coupled pore-fluid flow and chemical dissolution system is in a stable state, the trivial and nontrivial solutions of the system are identical. However, if a coupled pore-fluid flow and chemical dissolution system is in an unstable state, then the trivial and nontrivial solutions of the system are totally different. This recognition can be equally used to judge whether a coupled pore-fluid flow and chemical dissolution system involving heterogeneous porous media is in a stable state or in an unstable state. The proposed solution strategy can produce dual solutions for simulating coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media. Originality/value: A solution strategy is proposed to obtain the nontrivial solution, which is often overlooked in the computational simulation of coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media. The proposed solution strategy provides a useful way for understanding the underlying dynamic mechanisms of the chemical damage effect associated with the stability of structures that are built on soil foundations. [ABSTRACT FROM AUTHOR]

Purpose: When bending a large diameter thin-walled tube, the thickn ess of outer side wall will reduce greatly, which leads to a decrease of structural strength of the tube. To solve this problem, this paper investigated the deformation principles of an eccentric tube in the rotary draw bending process, trying to find a way to reduce the wall thickness difference between inner and outer diameters. Design/methodology/approach: An finite element model is established for analyzing the deformation of an eccentric tube in rotary draw bending process. The wall thickness distribution of the formed pipe was analyzed along the axis and diameter, respectively. Findings: It is found that there exists an optimal eccentricity between the inner and outer circle center of the tube cross-section. If the eccentricity of the tube is chosen properly, it is possible to get a bent tube with equal thickness of inner and outer side walls. In addition, it is also found the optimal eccentricity on the cross-section can be influenced by bending radius, wall thickness, diameter and bending angle. The optimal eccentricity increases greatly with the decreasing of bending radius, the increase of outer diameter and the increase of wall thickness. The influence of bending angle on the optimal eccentricity can be divided into two situations. When the bending angle is small, the optimal eccentricity increases with the increase of bending angle. When the bending angle exceeds a certain value, the pipe enters a stable forming state. The optimal eccentricity of the stable forming region does not change with the bending angle. Originality/value: Such a research is beneficial for reducing the thickness difference between inner and outer side walls in the rotary draw bending process. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to contribute to the solution of the fatigue damage problem of reinforced concrete frames in bending. Design/methodology/approach: The problem of fatigue damage is formulated based on the rheological–dynamical analogy, including a scalar damage variable to address the reduction of stiffness in strain softening. The modal analysis is used by the finite element method for the determination of modal parameters and resonance stability of the selected frame cross-section. The objectivity of the presented method is verified by numerical examples, predicting the ductility in bending of the frame whose basic mechanical properties were obtained by non-destructive testing systems. Findings: The modal analysis in the frame of the finite element method is suitable for the determination of modal parameters and resonance stability of the selected frame cross-section. It is recommended that the modulus of elasticity be determined by non-destructive methods, e.g. from the acoustic response. Originality/value: The paper presents a novel method of solving the ductility in bending taking into account both the creep coefficient and the aging coefficient. The rheological-dynamical analogy (RDA) method uses the resonant method to find material properties. The characterization of the structural damping via the damping ratio is original and effective. [ABSTRACT FROM AUTHOR]

Purpose: The response of an idealized triangular concrete gravity dam is studied due to horizontal and vertical ground motions for both fully reflective and absorptive reservoir bottom conditions. For each combination, in this paper different orders of Givoli-Neta (G-N) high-order truncation condition are aimed to be evaluated from accuracy point of view by comparing the results against corresponding exact solutions which relies on utilizing a two-dimensional fluid hyper-element. Design/methodology/approach: In present study, the dynamic analysis of concrete gravity dam-reservoir systems is formulated by Finite Element (FE)-(FE-TE) approach. In this technique, dam and reservoir are discretized by plane solid and fluid finite elements. Moreover, the G-N high-order condition imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model. Findings: It was observed that trend in gaining accuracy with increase in the order of G-N condition were basically the same for both horizontal and vertical ground motions under full reflective reservoir bottom condition. Moreover, convergence rate increases for absorptive reservoir bottom condition cases in comparison with fully reflective cases. It is also noticed that in certain cases, the responses are hardly distinguishable from corresponding exact responses. This reveals that proposed FE-(FE-TE) analysis technique based on G-N condition is quite successful, and one may fully rely on that for accurate and efficient analysis of concrete gravity dam-reservoir systems. Originality/value: Dynamic analysis of concrete gravity dam-reservoir systems are formulated by a new method. The salient aspect of the technique is that it utilizes G-N high-order condition at the truncation boundary. This is achieved by developing a special truncation element which its generalized matrices are derived for Finite Element Method (FEM) programmers. The method is discussed for all types of excitation and reservoir bottom conditions. It must be emphasized that although time harmonic analysis is considered in the present study, the main part of formulation is explained in the context of time domain. Therefore, the approach can easily be extended for transient type of analysis. [ABSTRACT FROM AUTHOR]

Purpose: The objective of this paper is to provide a better understanding of the effect of irregular columnar jointed structure on the permeability and flow characteristics of rock masses. Design/methodology/approach: An efficient numerical procedure is proposed to investigate the permeability and fluid flow in columnar jointed rock masses (CJRMs), of which the columnar jointed networks are generated by a modified constrained centroid Voronoi algorithm according to the field statistical results. The fractures are represented explicitly by using the lower-dimensional zero thickness elements. And the modeling scheme is validated by a benchmark test for flow in fractured porous media. The effective permeability and representative elementary volume (REV) size of CJRMs are estimated using finite element method (FEM). The influences of joint density and variation coefficient of columnar joint structure on the permeability of the rock mass are discussed. Findings: The simulation results indicate that the permeability is scale-dependent and tends to be stable with increase of model size. The hydraulic REV size is determined as 3.5 m for CJRMs in the present study. Moreover, the joint density is a dominant factor affecting the permeability of CJRMs. The average permeability of columnar jointed structures increases linearly with the joint density under the same REV size, while the influence from the coefficient of variation can be neglected. Originality/value: The present paper investigates the REV size of the CJRMs and the effect of joint parameters on the permeability. The proposed method and the results obtained are useful on understanding the hydraulic characteristic of the irregular CJRMs in engineering projects. [ABSTRACT FROM AUTHOR]

Purpose: An improved adaptive finite element analysis based on local error estimate is proposed via the element energy projection (EEP) technique. This paper aims to discuss the aforementioned idea. Design/methodology/approach: The computational region for a posteriori error estimation based on EEP method is further confined to a critical set of local elements generated in the previous adaptive step, enhancing efficiency while maintaining accuracy. The adaptive procedure incorporated with hierarchical mesh refinement is then developed. Findings: The effectiveness of the improved error estimation of the overall adaptive analysis is confirmed by several benchmark examples. The results show that the shrinkage of the local computational region has little negative influence on the accuracy of a posteriori error estimation, thus yielding an improved adaptive procedure with simplified logic and reduced cost. Originality/value: By localizing the computational region for error estimation, two crucial but cumbersome tricks, i.e. treatments of virtual elements and hanging nodes, are removed, giving the proposed approach full clarity and flexibility. The improved adaptive procedure characterizes simpler and faster computational algorithm and can produce results with required accuracy measured in maximum norm. [ABSTRACT FROM AUTHOR]

Purpose: The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect. Design/methodology/approach: In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method. Findings: The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised. Originality/value: A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to propose a new boundary condition and a web-spline basis of finite element space approximation to remedy the problems of constraints due to homogeneous and non-homogeneous; Dirichlet boundary conditions. This paper considered the two-dimensional linear elasticity equation of Navier–Lamé with the condition CAB. The latter allows to have a total insertion of the essential boundary condition in the linear system obtained; without using a numerical method as Lagrange multiplier. This study have developed mixed finite element; method using the B-splines Web-spline space. These provide an exact implementation of the homogeneous; Dirichlet boundary conditions, which removes the constraints caused by the standard; conditions. This paper showed the existence and the uniqueness of the weak solution, as well as the convergence of the numerical solution for the quadratic case are proved. The weighted extended B-spline; approach have become a much more workmanlike solution. Design/methodology/approach: In this paper, this study used the implementation of weighted finite element methods to solve the Navier–Lamé system with a new boundary condition CA, B (Koubaiti et al., 2020), that generalises the well-known basis, especially the Dirichlet and the Neumann conditions. The novel proposed boundary condition permits to use a single Matlab code, which summarises all kind of boundary conditions encountered in the system. By using this model is possible to save time and programming recourses while reap several programs in a single directory. Findings: The results have shown that the Web-spline-based quadratic-linear finite elements satisfy the inf–sup condition, which is necessary for existence and uniqueness of the solution. It was demonstrated by the existence of the discrete solution. A full convergence was established using the numerical solution for the quadratic case. Due to limited regularity of the Navier–Lamé problem, it will not change by increasing the degree of the Web-spline. The computed relative errors and their rates indicate that they are of order 1/H. Thus, it was provided their theoretical validity for the numerical solution stability. The advantage of this problem that uses the CA, B boundary condition is associated to reduce Matlab programming complexity. Originality/value: The mixed finite element method is a robust technique to solve difficult challenges from engineering and physical sciences using the partial differential equations. Some of the important applications include structural mechanics, fluid flow, thermodynamics and electromagnetic fields (Zienkiewicz and Taylor, 2000) that are mainly based on the approximation of Lagrange. However, this type of approximation has experienced a great restriction in the level of domain modelling, especially in the case of complicated boundaries such as that in the form of curvilinear graphs. Recently, the research community tried to develop a new way of approximation based on the so-called B-spline that seems to have superior results in solving the engineering problems. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to propose the five-phase sphere equivalent model of recycled concrete, which can be used to deduce the theoretical formulas for the Poisson's ratio and effective elastic modulus. Design/methodology/approach: At a mesoscopic level, the equivalent model converts the interfacial layer, which consists of the new interfacial transition zone (ITZ), the old mortar and the old (ITZ), into a uniform equivalent medium. This paper deduces a strength expression for the interfacial transition zone at the microscopic level using the equivalent model and elastic theory. In addition, a new finite element method called the base force element method was used in this research. Findings: Through numerical simulation, it was found that the mechanical property results from the five-phase sphere equivalent model were in good agreement with those of the random aggregate model. Furthermore, the proposed model agree on quite well with the available experimental data. Originality/value: The equivalent model can eliminate the influence of the interfacial layer on the macroscopic mechanical properties, thereby improving the calculation accuracy and computational efficiency. The proposed model can also provide a suitable model for multi-scale calculations. [ABSTRACT FROM AUTHOR]

Purpose: For gravity dams built on foundations with directional joint sets, the seepage in the foundation possesses anisotropic characteristics and may have adverse effects on the foundation stability. A methodology for system reliability analysis of gravity dam foundations considering anisotropic seepage and multiple sliding surfaces is proposed in this paper. Design/methodology/approach: Anisotropic seepages in dam foundations are simulated using finite element method (FEM) with the equivalent continuum model (ECM), and their effect on dam foundation stability is involved by uplift pressures acting on the potential sliding surfaces. The system failure probability of the dam foundation is efficiently estimated using Monte Carlo method (MCM) combined with response surface method (RSM). Findings: The case study shows that it is necessary to consider the possibly adverse effect of anisotropic seepage on foundation stability of gravity dams and the deterministic analysis of the foundation stability may be misleading. The system reliability analysis of the dam foundation is justified, as the uncertainties in shear strength parameters of the foundation rocks and joint sets as well as aperture, connectivity and spacing of the joint sets are quantified and the system effect of the multiple potential sliding surfaces on the foundation reliability is reasonably considered. Originality/value: (1) A methodology is proposed for efficient system reliability analysis of foundation stability of gravity dams considering anisotropic seepage and multiple sliding surfaces (2) The influence of anisotropic seepage on the stability of gravity dam foundation is revealed (3) The influence of estimation errors of RSMs on the system reliability assessment of dam foundation is investigated. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to propose an effective and efficient numerical method that can consider natural frequency in multi-material topology optimization (MMTO) and which is scalable for complex three-dimensional (3D) problems. Design/methodology/approach: The optimization algorithm is developed by combining custom FORTRAN code for MMTO with the open-source software Mystran, which is used as a finite element analysis (FEA) solver. The proposed algorithm allows the designer to shift the fundamental frequency of the design beyond a defined frequency spectrum from the initial designing phase. The methodology is formulated in a smooth and differentiable manner, with the sensitivity expressions, required by gradient-based optimization solvers, presented. Findings: Natural frequency constraint has been successfully implemented into MMTO. The use of open-source software Mystran as an FEA solver in the algorithm provides ability to solve complex problems. Mystran offers powerful built-in functions for eigenvalue extraction using methods like Givens, modified Givens, inverse power and the Lanczos method, which provide the ability to solve complex models. The algorithm is successfully able to solve both two- and three-material MMTO jobs for two-dimensional and 3D geometries. Originality/value: Natural frequency constraint consideration into topology optimization is very challenging due to three common issues: localized eigenmodes, mode switching and high computational cost. The proposed algorithm addresses these inherent issues, implements natural frequency constraint to MMTO and solves for complex models, which is hardly possible using conventional methods. [ABSTRACT FROM AUTHOR]

Purpose: Integrality of surface mesh is requisite for computational engineering. Nonwatertight meshes with holes can bring inconvenience to applications. Unlike simple modeling or visualization, the downstream industrial application scenarios put forward higher requirements for hole-filling, although many related algorithms have been developed. This study aims at the hole-filling issue in industrial application scenarios. Design/methodology/approach: This algorithm overcomes some inherent weakness of general methods and generates a high-level resulting mesh. Initially, the primitive hole boundary is filled with a more appropriate triangulation which introduces fewer geometric errors. And in order for better performances on shape approximation of the background mesh, the algorithm also refines the initial triangulation with topology optimization. When obtaining the background mesh defining the geometry and size field, spheres on it are packed to determine the vertex configuration and then the resulting high-level mesh is generated. Findings: Through emphasizing geometry recovery and mesh quality, the proposed algorithm works well in hole-filling in industrial application scenarios. Many experimental results demonstrate the reliability and the performance of the algorithm. And the processed meshes are capable of being used for industrial simulation computations directly. Originality/value: This paper makes input meshes more adaptable for solving programs through local modifications on meshes and perfects the preprocessing technology of finite element analysis (FEA). [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this study is to develop a homogenization approach that ensures both high accuracy and time-efficient solution for elastic-plastic functionally graded composites. Design/methodology/approach: The paper presents a novel two-stage hybrid homogenization approach that combines advantages of the mean field homogenization and homogenization based on the finite element method (FEM). The groundbreaking nature of the developed approach is associated with division of the hybrid homogenization procedure into two stages, which allows to very efficiently determine the solution for arbitrary volume fraction of the reinforcement. This paper concerns also on modelling of composites with randomly distributed prolate and oblate particles. For this purpose, the hybrid homogenization was implemented in the framework of the discrete orientation averaging procedure involving pseudo-grain discretization method. Findings: Agreement between the results obtained using the proposed approach and the standard FEM-based homogenization is very good (up to the volume fraction of 0.3). Originality/value: The proposed two-stage homogenization approach allows to obtain the solution for materials with arbitrary volume fraction of the reinforcement very efficiently; therefore, it is highly beneficial for the two-scale modeling of nonlinear functionally graded materials and structures. [ABSTRACT FROM AUTHOR]

Purpose: The paper aims to presents a numerical analysis of free vibration of micromorphic structures subjected to various boundary conditions. Design/methodology/approach: To accomplish this objective, first, a two-dimensional (2D) micromorphic formulation is presented and the matrix representation of this formulation is given. Then, two size-dependent quadrilateral and triangular elements are developed within the commercial finite element software ABAQUS. User element subroutine (UEL) is used to implement the micromorphic elements. These non-classical elements are capable of capturing the micro-structure effects by considering the micro-motion of materials. The effects of the side length-to-length scale parameter ratio and boundary conditions on the vibration behavior of 2D micro-structures are discussed in detail. The reliability of the present finite element method (FEM) is confirmed by the convergence studies and the obtained results are validated with the results available in the literature. Also, the results of micromorphic theory (MMT) are compared with those of micropolar and classical elasticity theories. Findings: The study found that the size effect becomes very significant when the side length of micro-structures is close to the length scale parameter. Originality/value: The study is to analyze the free vibrations of 2D micro-structures based on MMT; to develop a 2D formulation for micromorphic continua within ABAQUS; to propose quadrilateral and triangular micromorphic elements using UEL and to investigate size effects on the vibrational behavior of micro-structures with various geometries. [ABSTRACT FROM AUTHOR]

Purpose: The objective of this paper is to develop a semi-analytical finite element method for solving chemical dissolution-front instability problems in fluid-saturated porous media. Design/methodology/approach: The porosity, horizontal and vertical components of the pore-fluid velocity and solute concentration are selected as four fundamental unknown variables for describing chemical dissolution-front instability problems in fluid-saturated porous media. To avoid the use of numerical integration, analytical solutions for the property matrices of a rectangular element are precisely derived in a purely mathematical manner. This means that the proposed finite element method is a kind of semi-analytical method. The column pivot element solver is used to solve the resulting finite element equations of the chemical dissolution-front instability problem. Findings: The direct use of horizontal and vertical components of the pore-fluid velocity as fundamental unknown variables can improve the accuracy of the related numerical solution. The column pivot element solver is useful for solving the finite element equations of a chemical dissolution-front instability problem. The proposed semi-analytical finite element method can produce highly accurate numerical solutions for simulating chemical dissolution-front instability problems in fluid-saturated porous media. Originality/value: Analytical solutions for the property matrices of a rectangular element are precisely derived for solving chemical dissolution-front instability problems in fluid-saturated porous media. The proposed semi-analytical finite element method provides a useful way for understanding the underlying dynamic mechanisms of the washing land method involved in the contaminated land remediation. [ABSTRACT FROM AUTHOR]

Purpose: In this paper the use of classical strain measures in analysis of trusses at finite deformations will be discussed. The results will be compared to the ones acquired using a novel strain measure based on the Hyperbolic Sine function. Through the evaluation of results, algebraic development and graph analysis, the properties of the Hyperbolic Sine strain measure will be examined. Design/methodology/approach: Through graph plotting, comparisons between the novel strain measure and the classic ones will be made. The formulae for the implementation of the Hyperbolic Sine strain measure into a positional finite element method are developed. Four engineering applications are presented and comparisons between results obtained using all strain measures studied are made. Findings: The proposed strain measure, Hyperbolic Sine, has objectivity and symmetry. The linear constitutive model formed by the Hyperbolic Sine strain and its conjugated stress presents an increasing stiffness, both in compression and tension, a behavior that can be useful in the modeling of several materials. Research limitations/implications: The structural analysis performed on the four examples of trusses in this article did not consider the variation of the cross-sectional area of the elements or the buckling phenomenon, moreover, only elastic behavior is considered. Originality/value: The present article proposes the use of a novel strain measure family, based on the Hyperbolic Sine function and suitable for structural applications. Mathematical expressions for the use of the Hyperbolic Sine strain measure are established following the energetic concepts of the positional formulation of the finite element method. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to present a skin-based approach able to generate the variability model for a component in composite material due to its manufacturing process. It generates a skin-based model of the manufactured part. The skin model discretizes the part surfaces by points to take into account the geometric deviations, those points are the nodes of finite element analysis used for tolerance analysis of compliant assemblies. Design/methodology/approach: The paper presents a general and systematic simulation model for generating a variability meta-model for a component in composite material due to its manufacturing process. The model is constituted by three steps: definition and pre-processing of the nominal model, generation of the manufacturing process model and evaluation of the part variability. Findings: The advantage of this approach is related to the fact that it is designed as a part of a digital process that establishes a continuous and unambiguous flow of variation information from the part design to manufacturing and assembly and that takes into account the manufacturing signature. This is its uniqueness compared to other simulation approaches focused only on manufacturing. Research limitations/implications: Considering the variability around the nominal value of all the process parameters and parts with more complex geometries are not taken into account now, which will be modelled in practical applications. Practical implications: To properly manage uncertainty since conceptual design of complex product, next generation geometry assurance requires simulation models to realistically consider process signatures due to the manufacturing process. This work focusses on this next generation tool for geometry assurance. Originality/value: The literature is focused on metal sheets joined by welding or riveting. There are other materials widely used and typically compliant: the composite materials that typically used mechanical fixing elements (bolting, riveting) and structural adhesives to joint parts. No software tools exist in the literature to deal with uncertainty from manufacturing to assembly processes in products made by composite. This is the reason of the present work. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to propose a new scheme for obtaining acceptable solutions for problems of continuum topology optimization of structures, regarding the distribution and limitation of discretization errors by considering h-adaptivity. Design/methodology/approach: The new scheme encompasses, simultaneously, the solution of the optimization problem considering a solid isotropic microstructure with penalization (SIMP) and the application of the h-adaptive finite element method. An analysis of discretization errors is carried out using an a posteriori error estimator based on both the recovery and the abrupt variation of material properties. The estimate of new element sizes is computed by a new h-adaptive technique named "Isotropic Error Density Recovery", which is based on the construction of the strain energy error density function together with the analytical solution of an optimization problem at the element level. Findings: Two-dimensional numerical examples, regarding minimization of the structure compliance and constraint over the material volume, demonstrate the capacity of the methodology in controlling and equidistributing discretization errors, as well as obtaining a great definition of the void–material interface, thanks to the h-adaptivity, when compared with results obtained by other methods based on microstructure. Originality/value: This paper presents a new technique to design a mesh made with isotropic triangular finite elements. Furthermore, this technique is applied to continuum topology optimization problems using a new iterative scheme to obtain solutions with controlled discretization errors, measured in terms of the energy norm, and a great resolution of the material boundary. Regarding the computational cost in terms of degrees of freedom, the present scheme provides approximations with considerable less error if compared to the optimization process on fixed meshes. [ABSTRACT FROM AUTHOR]

Purpose: The paper aims to propose two new 8-node quadrilateral membrane elements with good distortion tolerance for the modified couple stress elasticity based on the unsymmetric finite element method (FEM). Design/methodology/approach: The nodal rotation degrees of freedom (DOFs) are introduced into the virtual work principle and constrained by the penalty function for approximating the test functions of the physical rotation and curvature. Therefore, only the C0 continuity instead of C1 continuity is required for the displacement during the element construction. The first unsymmetric element assumes the test functions of the displacement and strain using the standard 8-node isoparametric interpolations, while these test functions in the second model are further enhanced by the nodal rotation DOFs. Besides, the trial functions in these two elements are constructed based on the stress functions that can a priori satisfy related governing equations. Findings: The benchmark tests show that both the two elements can efficiently simulate the size-dependent plane problems, exhibiting good numerical accuracies and high distortion tolerances. In particular, they can still exactly reproduce the constant couple stress state when the element shape deteriorates severely into the degenerated triangle. Moreover, it can also be observed that the second element model, in which the linked interpolation technique is used, has better performance than the first one, especially in capturing the steep gradients of the physical rotations. Originality/value: As the proposed new elements use only three DOFs per node, they can be readily incorporated into the existing finite element (FE) programs. Thus, they are of great benefit to analysis of size-dependent membrane behaviors of micro/nano structures. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to analyze the stability of laminated shells subjected to axial loads or external pressure with considering various geometries and boundary conditions. The main aim of the present study is developing an efficient combined method which uses the advantages of different methods, such as finite element method (FEM) and isogeometric analysis (IGA), to achieve multipurpose targets. Two types of material including laminated composite and sandwich functionally graded material are considered. Design/methodology/approach: A novel type of finite strip method called isogeometric B3-spline finite strip method (IG-SFSM) is used to solve the eigenvalue buckling problem. IG-SFSM uses B3-spline basis functions to interpolate the buckling displacements and mapping operations in the longitudinal direction of the strips, whereas the Lagrangian functions are used in transverse direction. The current presented IG-SFSM is formulated based on the degenerated shell method. Findings: The buckling behavior of laminated shells is discussed by solving several examples corresponding to shells with various geometries, boundary conditions and material properties. The effects of mechanical and geometrical properties on critical loads of shells are investigated using the related results obtained by IG-SFSM. Originality/value: This paper shows that the proposed IG-SFSM leads to the critical loads with an approved accuracy comparing with the same examples extracted from the literature. Moreover, it leads to a high level of convergence rate and low cost of solving the stability problems in comparison to the FEM. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to present an improved finite element method used for achieving faster convergence in simulations of incompressible fluid flows. For stable computations of incompressible fluid flows, it is important to ensure that the flow field satisfies the equation of continuity in each element of a generally distorted mesh. The study aims to develop a numerical approach that satisfies this requirement based on the highly simplified marker-and-cell (HSMAC) method and increases computational speed by introducing a new algorithm into the simultaneous relaxation of velocity and pressure. Design/methodology/approach: First, the paper shows that the classical HSMAC method is equivalent to a Jacobi-type method in terms of the simultaneous relaxation of velocity and pressure. Then, a Gauss–Seidel or successive over-relaxation (SOR)-type method is introduced in the Newton–Raphson iterations to take into account all the derivative terms in the first-order Taylor series expansion of a nodal-averaged error explicitly. Here, the nine-node quadrilateral (Q2–Q1) elements are used. Findings: The new finite element approach based on the improved HSMAC algorithm is tested on fluid flow problems including the lid-driven square cavity flow and the flow past a circular cylinder. The results show significant improvement of the convergence property with the accuracy of the numerical solutions kept unchanged even on a highly distorted mesh. Originality/value: To the best of the author's knowledge, the idea of using the Gauss–Seidel or SOR method in the simultaneous relaxation procedure of the HSMAC method has not been proposed elsewhere. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to propose a new strong-form numerical method, called the free element method, for solving general boundary value problems governed by partial differential equations. The main idea of the method is to use a locally formed element for each point to set up the system of equations. The proposed method is used to solve the fluid mechanics problems. Design/methodology/approach: The proposed free element method adopts the isoparametric elements as used in the finite element method (FEM) to represent the variation of coordinates and physical variables and collocates equations node-by-node based on the newly derived element differential formulations by the authors. The distinct feature of the method is that only one independently formed individual element is used at each point. The final system of equations is directly formed by collocating the governing equations at internal points and the boundary conditions at boundary points. The method can effectively capture phenomena of sharply jumped variables and discontinuities (e.g. the shock waves). Findings: a) A new numerical method called the FEM is proposed; b) the proposed method is used to solve the compressible fluid mechanics problems for the first time, in which the shock wave can be naturally captured; and c) the method can directly set up the system of equations from the governing equations. Originality/value: This paper presents a completely new numerical method for solving compressible fluid mechanics problems, which has not been submitted anywhere else for publication. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to present a novel scheme for imposing periodic boundary conditions with downscaled macroscopic strain measures of gradient Cosserat continuum on the representative volume element (RVE) of discrete particle assembly in the frame of the second-order computational homogenization methods for granular materials. Design/methodology/approach: The proposed scheme is based on the generalized Hill's lemma of gradient Cosserat continuum and the incremental non-linear constitutive relation condensed to the peripheral particles of the RVE of discrete particle assembly. The generalized Hill's lemma conducts to downscale the macroscopic strain or stress measures and to impose the periodic boundary conditions on the RVE boundary so that the Hill-Mandel energy equivalence condition is ensured. Because of the incremental non-linear constitutive relation condensed to the peripheral particles of the RVE, the periodic boundary displacement and traction constraints together with the downscaled macroscopic strains and strain gradients, micro-rotations and curvatures are imposed in the point-wise sense without the need of introducing the Lagrange multipliers for enforcing the periodic boundary displacement and traction constraints in a weak sense. Findings: Numerical results demonstrate that the applicability and effectiveness of the proposed scheme in imposing the periodic boundary conditions on the RVE. The results of the RVE subjected to the periodic boundary conditions together with the displacement boundary conditions in the second-order computational homogenization for granular materials provide the desired estimations, which lie between the upper and the lower bounds provided by the displacement and the traction boundary conditions imposed on the RVE respectively. Research limitations/implications: Each grain in the particulate system under consideration is assumed to be rigid and circular. Practical implications: The proposed scheme for imposing periodic boundary conditions on the RVE can be adopted solely for estimating the effective mechanical properties of granular materials and/or integrated into the frame of the second-order computational homogenization method with a nested finite element method-discrete element method solution procedure for granular materials. It will tend to provide, at least theoretically, more reasonable results for effective material properties and solutions of a macroscopic boundary value problem simulated by the computational homogenization method. Originality/value: This paper presents a novel scheme for imposing periodic boundary conditions with downscaled macroscopic strain measures of gradient Cosserat continuum on the RVE of discrete particle assembly for granular materials without need of introducing Lagrange multipliers for enforcing periodic boundary conditions in a weak (integration) sense. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to study the energy absorption properties of the thin-walled square tube with lateral piecewise variable thickness under axial crashing and the influence of the tube parameters on energy absorption. Design/methodology/approach: In this work, the energy absorption properties of the thin-walled square tube were analyzed by theoretical, numerical and experimental approach. The numerical results are obtained based on the finite element method. The explicit formulation for predicting the mean crushing force of the tube with lateral piecewise variable thickness was derived based on Super Folding Element method. The limitation of the prediction formulation was analyzed by numerical calculation. The numerical calculation was also used to compare the energy absorption between the tube with lateral piecewise variable thickness and other tubes, and to carry out the parametric analysis. Findings: Results indicate that the thin-walled tube with lateral piecewise variable thickness has higher energy absorption properties than the uniform thickness tubes and the tubes with lateral linear variable thickness. The thickness of the corner is the key factor for the energy absorption of the tubes. The thickness of the non-corner region is the secondary factor. Increasing the corner thickness and decreasing the non-corner thickness can make the energy absorption improved. It is also found that the prediction formulation of the mean crushing force given in this paper can quickly and accurately predict the energy absorption of the square tube. Originality/value: The outcome of the present research provides a design idea to improve the energy absorption of thin-walled tube by designing cross-section thickness and gives an explicit formulation for predicting the mean crushing force quickly and accurately. [ABSTRACT FROM AUTHOR]

Purpose: This meshless collocation method is applicable not only to the Helmholtz equation with Dirichlet boundary condition but also mixed boundary conditions. It can calculate not only the high wavenumber problems, but also the variable wave number problems. Design/methodology/approach: In this paper, the authors developed a meshless collocation method by using barycentric Lagrange interpolation basis function based on the Chebyshev nodes to deduce the scheme for solving the three-dimensional Helmholtz equation. First, the spatial variables and their partial derivatives are treated by interpolation basis functions, and the collocation method is established for solving second order differential equations. Then the differential matrix is employed to simplify the differential equations which is on a given test node. Finally, numerical experiments show the accuracy and effectiveness of the proposed method. Findings: The numerical experiments show the advantages of the present method, such as less number of collocation nodes needed, shorter calculation time, higher precision, smaller error and higher efficiency. What is more, the numerical solutions agree well with the exact solutions. Research limitations/implications: Compared with finite element method, finite difference method and other traditional numerical methods based on grid solution, meshless method can reduce or eliminate the dependence on grid and make the numerical implementation more flexible. Practical implications: The Helmholtz equation has a wide application background in many fields, such as physics, mechanics, engineering and so on. Originality/value: This meshless method is first time applied for solving the 3D Helmholtz equation. What is more the present work not only gives the relationship of interpolation nodes but also the test nodes. [ABSTRACT FROM AUTHOR]

Purpose: Various time integration methods and time finite element methods have been developed to obtain the responses of structural dynamic problems, but the accuracy and computational efficiency of them are sometimes not satisfactory. The purpose of this paper is to present a more accurate and efficient formulation on the basis of the weak form quadrature element method to solve linear structural dynamic problems. Design/methodology/approach: A variational principle for linear structural dynamics, which is inspired by Noble's work, is proposed to develop the weak form temporal quadrature element formulation. With Lobatto quadrature rule and the differential quadrature analog, a system of linear equations is obtained to solve the responses at sampling time points simultaneously. Computation for multi-elements can be carried out by a time-marching technique, using the end point results of the last element as the initial conditions for the next. Findings: The weak form temporal quadrature element formulation is conditionally stable. The relation between the normalized length of element and the suggested number of integration points in one element is given by a simple formula. Results show that the present formulation is much more accurate than other time integration methods and its dissipative property is also illustrated. Originality/value: The weak form temporal quadrature element formulation provides a choice with high accuracy and efficiency for solution of linear structural dynamic problems. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to extend a recent unsymmetric 8-node, 24-DOF hexahedral solid element US-ATFH8 with high distortion tolerance, which uses the analytical solutions of linear elasticity governing equations as the trial functions (analytical trial function) to geometrically nonlinear analysis. Design/methodology/approach: Based on the assumption that these analytical trial functions can still properly work in each increment step during the nonlinear analysis, the present work concentrates on the construction of incremental nonlinear formulations of the unsymmetric element US-ATFH8 through two different ways: the general updated Lagrangian (UL) approach and the incremental co-rotational (CR) approach. The key innovation is how to update the stresses containing the linear analytical trial functions. Findings: Several numerical examples for 3D structures show that both resulting nonlinear elements, US-ATFH8-UL and US-ATFH8-CR, perform very well, no matter whether regular or distorted coarse mesh is used, and exhibit much better performances than those conventional symmetric nonlinear solid elements. Originality/value: The success of the extension of element US-ATFH8 to geometrically nonlinear analysis again shows the merits of the unsymmetric finite element method with analytical trial functions, although these functions are the analytical solutions of linear elasticity governing equations. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to analyze the influence of rotated anisotropy on the stability of slope, the random finite element method is used in this study. Design/methodology/approach: The random field is generated by the discrete cosine transform (DCT) method, which can generate random field with different rotated angles conveniently. Findings: Two idealized slopes are analyzed; it is observed that the rotated angle significantly affects the slope failure risk. The two examples support the conclusion that when the orientation of the layers is nearly perpendicular to the slip surface, the slope is in a relative stable condition. The results of heterogeneous slope with two clay layers demonstrate that the rotated angle of lower layer mainly controls the failure mechanism of the slope, and the rotated angle of upper layer exhibits a significant influence on the probability of slope failure. Originality/value: The method for rotated anisotropy random field generation based on the DCT has a simple expression with few parameters and is convenient for implementation and practical application. The proposed method and the results obtained are useful for analyzing the stability of the heterogeneous slopes in engineering projects. [ABSTRACT FROM AUTHOR]

Purpose: The finite element method (FEM) is used to calculate the two-dimensional anti-plane dynamic response of structure embedded in D'Alembert viscoelastic multilayered soil on the rigid bedrock. This paper aims to research a time-domain absorbing boundary condition (ABC), which should be imposed on the truncation boundary of the finite domain to represent the dynamic interaction between the truncated infinite domain and the finite domain. Design/methodology/approach: A high-order ABC for scalar wave propagation in the D'Alembert viscoelastic multilayered media is proposed. A new operator separation method and the mode reduction are adopted to construct the time-domain ABC. Findings: The derivation of the ABC is accurate for the single layer but less accurate for the multilayer. To achieve high accuracy, therefore, the distance from the truncation boundary to the region of interest can be zero for the single layer but need to be about 0.5 times of the total layer height of the infinite domain for the multilayer. Both single-layered and multilayered numerical examples verify that the accuracy of the ABC is almost the same for both cases of only using the modal number excited by dynamic load and using the full modal number of infinite domain. Using the ABC with reduced modes can not only reduce the computation cost but also be more friendly to the stability. Numerical examples demonstrate the superior properties of the proposed ABC with stability, high accuracy and remarkable coupling with the FEM. Originality/value: A high-order time-domain ABC for scalar wave propagation in the D'Alembert viscoelastic multilayered media is proposed. The proposed ABC is suitable for both linear elastic and D'Alembert viscoelastic media, and it can be coupled seamlessly with the FEM. A new operator separation method combining mode reduction is presented with better stability than the existing methods. [ABSTRACT FROM AUTHOR]

Purpose: A general first-invariant constitutive model has been derived in literature for incompressible, isotropic hyperelastic materials, known as Marlow model, which reproduces test data exactly without the need of curve-fitting procedures. This paper aims to describe how to extend Marlow's constitutive model to the more general case of compressible hyperelastic materials. Design/methodology/approach: The isotropic constitutive model is based on a strain energy function, whose isochoric part is solely dependent on the first modified strain invariant. Based on Marlow's idea, a principle of energetically equivalent deformation states is derived for the compressible case, which is used to determine the underlying strain energy function directly from measured test data. No particular functional of the strain energy function is assumed. It is shown how to calibrate the volumetric and isochoric strain energy functions uniquely with uniaxial or biaxial test data only. The constitutive model is implemented into a finite element program to demonstrate its applicability. Findings: The model is well suited for use in finite element analysis. Only one set of test data is required for calibration without any need for curve-fitting procedures. These test data are reproduced exactly, and the model prediction is reasonable for other deformation modes. Originality/value: Marlow's basic concept is extended to the compressible case and applied to both the volumetric and isochoric part of the compressible strain energy function. Moreover, a novel approach is described on how both compressive and tensile test data can be used simultaneously to calibrate the model. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to analyze the liquid sloshing behaviors in two-dimensional tanks with various porous baffles under the external excitation. Design/methodology/approach: Adopting the finite element method (FEM) and control variable method to study the impacts of the height, length, number, location, shape, porous-effect parameter of the porous baffle, the external load frequency and the shape of the tank on the liquid sloshing response. Findings: The amplitude of the free surface can be reduced effectively when the baffle opening is appropriate. The anti-sway ability of the system increases in pace with the baffle's height growing. Under the same conditions, the shapes of the baffles have an important effect on improving the anti-sway ability of the system. Originality/value: As there exist the differences of the velocity potential between each side of the porous baffle, which means that there are two different velocity potentials at a point on the porous baffle, the conventional finite element modeling technologies are not suitable to be applied here. To deal with this problem, the points on the porous baffle are regarded as two nodes with the same coordinate to model and calculate. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to quantitatively investigate the effect of process parameters (including welding current, voltage and speed) and plate thickness on in-plane inherent deformations in typical fillet welded joint; meanwhile, the plastic strains remaining in the weld zone are also analyzed under different influencing factors. Design/methodology/approach: To achieve the purpose of this study, a thermal-elastic-plastic finite element (TEP FE) model is developed to analyze the thermal-mechanical behavior of the T-welded joint during the welding process. Experimental measurements have verified the validity of the established TEP FE model. Using the effective model, a series of numerical experiments are performed to obtain the inherent deformations under the conditions of different influencing factors, and then the calculation results are discussed based on the relevant data obtained. Findings: Through numerical simulation analysis, it is found that the longitudinal and transverse inherent deformations decrease with the increase of welding speed and plate thickness, whereas as the nominal heat input increases, the inherent deformations increase significantly. The longitudinal shrinkage presents a quasi-linear and nonlinear distribution in the middle and end of the weld, respectively. The plastic strains in the cross section of the T-joint also vary greatly because of the process parameters and plate thickness, but the maximum value always appears near the location of the welding toe, which means that this point faces a relatively large risk of fatigue cracking. The inherent deformations are closely related to the plastic strains remaining in the weld zone and are also affected by many influencing factors such as process parameters and plate thickness. Research limitations/implications: In this study, relatively few influencing factors such as welding current, voltage, speed and plate thickness are considered to analyze the inherent deformations in the T-welded joint. Also, these influencing factors are all within a certain range of parameters, which shows that only limited applicability can be provided. In addition, only in-plane inherent deformations are considered in this study, without considering the other two out-of-plane components of inherent deformations. Originality/value: This study can help to expand the understanding of the relationship between the inherent deformations and its influencing factors for a specific form of the welded joint, and can also provide basic data to supplement the inherent deformation database, thereby facilitating further researches on welding deformations for stiffened-panel structures in shipbuilding or steel bridges. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to synthetically investigate the impact damage responses of carbon fiber reinforced polymer (CFRP) and its influence on the compression mechanical responses of CFRP laminates, including damage distribution, residual compressive strength and fracture morphology. Design/methodology/approach: A progressive damage simulation model is developed to analyze the complicated damage responses of CFRP laminates that are manufactured by resin transfer method (RTM) technology. Based on the ABAQUS/explicit finite element analysis solver, a VUMAT code is proposed to descript the composite materials' damage behaviors under both impact and compression load. Adopting this proposed model, the primary mechanical indicators of four groups' 5284RTM/U3160 CFRP laminates with different stacking sequences are predicted. Moreover, impact and compression after impact tests are conducted to verify the accuracy of simulation results. Findings: Both simulation and experimental results show that the impact damage with low visible detectability can significantly reduce composites' compressive strength. For all four groups' composite laminates, the residual strength ratio is around 35% or even lower. The kernel impact damage near the plates' geometric center promotes the degradation process of local materials and finally leads to the early occurrence of mechanical fracture. In addition, the impact damage projection area is not sensitive to the parameters of stacking sequences, while the residual compression strength is proportional to the number of 0-degree layers within whole laminates. Originality/value: This study helps to understand the effect of an impact event on CFRP laminates' compressive bearing capacity and provides a numerical method in simulating the damage responses under both impact and compression load. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to obtain the single-tooth stiffness, single-tooth time-varying meshing stiffness and comprehensive meshing stiffness of the internal and external face gears and to analyze the influence of the modulus, pressure angle and tooth width of each face gear on the single-tooth stiffness of the gear in nutation face gear transmission. Design/methodology/approach: From the point of view of material mechanics, the gear teeth of nutation face gear are simplified as spacial variable cross-section beams. The shear deformation of gear teeth, the bending deformation of tooth root and the additional elastic deformation caused by the base deformation are gotten by simplified trapezoidal section method, thus the stiffness of nutation face gear teeth can be obtained. The comparison with finite element method results verifies the rationality of simplified trapezoidal section method for calculating the tooth stiffness of nutation face gear. Findings: The variation of stiffness of internal and external face gears along the meshing line and tooth height in nutation face gear transmission is studied, and the variation laws of single tooth stiffness, single-tooth-pair mesh stiffness and single tooth time-varying meshing stiffness of nutation face gear teeth are obtained. Originality/value: Nutation face gear transmission is a new type of transmission. The stiffness of face gear teeth is analyzed, and the variation rules of single tooth stiffness, single-tooth-pair mesh stiffness and single tooth time-varying meshing stiffness of nutation face gear teeth are obtained, which not only enriches the research of nutation face gear transmission but also has important guiding significance for the application of nutation face gear in engineering practice. [ABSTRACT FROM AUTHOR]

Purpose: This study aims to overcome the involved challenging issues and provide high-precision eigensolutions. General eigenproblems in the system of ordinary differential equations (ODEs) serve as mathematical models for vector Sturm-Liouville (SL) and free vibration problems. High-precision eigenvalue and eigenfunction solutions are crucial bases for the reliable dynamic analysis of structures. However, solutions that meet the error tolerances specified are difficult to obtain for issues such as coefficients of variable matrices, coincident and adjacent approximate eigenvalues, continuous orders of eigenpairs and varying boundary conditions. Design/methodology/approach: This study presents an h-version adaptive finite element method based on the superconvergent patch recovery displacement method for eigenproblems in system of second-order ODEs. The high-order shape function interpolation technique is further introduced to acquire superconvergent solution of eigenfunction, and superconvergent solution of eigenvalue is obtained by computing the Rayleigh quotient. Superconvergent solution of eigenfunction is used to estimate the error of finite element solution in the energy norm. The mesh is then, subdivided to generate an improved mesh, based on the error. Findings: Representative eigenproblems examples, containing typical vector SL and free vibration of beams problems involved the aforementioned challenging issues, are selected to evaluate the accuracy and reliability of the proposed method. Non-uniform refined meshes are established to suit eigenfunctions change, and numerical solutions satisfy the pre-specified error tolerance. Originality/value: The proposed combination of methodologies described in the paper, leads to a powerful h-version mesh refinement algorithm for eigenproblems in system of second-order ODEs, that can be extended to other classes of applications in damage detection of multiple cracks in structures based on the high-precision eigensolutions. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to distinguish the relationship between the morphology characteristics of different scales and the contact performance of the mating surfaces. Also, an integrated method of the spectrum analysis and the wavelet transform is used to separate the morphology characteristics of the actual machined parts. Design/methodology/approach: First, a three-dimensional (3D) surface profilometer is used to obtain the surface morphology data of the actual machined parts. Second, the morphology characteristics of different scales are realized by the wavelet analysis and the power spectral density. Third, the reverse modeling engineering is used to construct the 3D contact models for the macroscopic characteristics. Finally, the finite element method is used to analyze the contact stiffness and the contact area of the 3D contact model. Findings: The contact area and the nominal contact pressure Pn have a nonlinear relationship in the whole compression process for the 3D contact model. The percentage of the total contact area of the macro-scale mating surface is about 70 per cent when the contact pressure Pn is in the range of 0-100 MPa, and the elastic contact area accounts for the vast majority. Meanwhile, when the contact pressure Pn is less than 10MPa, the influence factor (the relative error of contact stiffness) is larger than 50 per cent, so the surface macro-scale morphology has a weakening effect on the normal contact stiffness of the mating surfaces. Originality/value: This paper provides an effective method for the multi-scale separation of the surface morphology and then lays a certain theoretical foundation for improving the surface quality of parts and the morphology design. [ABSTRACT FROM AUTHOR]

Purpose: This paper aims to evaluate the dynamic response of surrounding foundation and study the vibration characteristics of track system. Design/methodology/approach: A double-line underground coupling analysis model was established, which included two moving train, track, liner and the ground field. Findings: Based on the 2.5D (D is diameter) finite element analysis, the influence of the important factors such as the depth of the subway tunnel, the nature of the foundation soil, the relative position relation of the double tunnel, the subway driving speed on the foundation and the orbital vibration are analyzed in this article. Originality/value: The results in paper may have reference value for the prediction of train induced vibrations and for the research of dynamic response of ground field. [ABSTRACT FROM AUTHOR]

Purpose The purpose of this paper is to discuss about evaluating the integrals involving B-spline wavelet on the interval (BSWI), in wavelet finite element formulations, using Gauss Quadrature.Design/methodology/approach In the proposed scheme, background cells are placed over each BSWI element and Gauss quadrature rule is defined for each of these cells. The nodal discretization used for BSWI WFEM element is independent to the selection of number of background cells used for the integration process. During the analysis, background cells of various lengths are used for evaluating the integrals for various combination of order and resolution of BSWI scaling functions. Numerical examples based on one-dimensional (1D) and two-dimensional (2D) plane elasto-statics are solved. Problems on beams based on Euler Bernoulli and Timoshenko beam theory under different boundary conditions are also examined. The condition number and sparseness of the formulated stiffness matrices are analyzed.Findings It is found that to form a well-conditioned stiffness matrix, the support domain of every wavelet scaling function should possess sufficient number of integration points. The results are analyzed and validated against the existing analytical solutions. Numerical examples demonstrate that the accuracy of displacements and stresses is dependent on the size of the background cell and number of Gauss points considered per background cell during the analysis.Originality/value The current paper gives the details on implementation of Gauss Quadrature scheme, using a background cell-based approach, for evaluating the integrals involved in BSWI-based wavelet finite element method, which is missing in the existing literature. [ABSTRACT FROM AUTHOR]

Purpose The traditional precision design only takes the influence of geometric tolerance of the parts and does not involve the load deformation in the assembly process. This paper aims to analyze the influence mechanism of flexible parts deformation on the geometric precision, and then to ensure the reliability and stability of the mechanical system.Design/methodology/approach Firstly, this paper adopts the N-GPS to analyze the influence mechanism of flexible parts deformation on the geometric precision and constructs a coupling 3D tolerance mathematical model of the geometric tolerance and the load deformation deviation based on the SDT theory, homogeneous coordinate transformation theory and surface authentication idea. Secondly, the least square method is used to fit the deformation surface of the mating surface under load so as to complete the conversion from the non-ideal element to the ideal element.Findings This paper takes the horizontal machining center as a case to obtain the deformation information of the mating surface under the self-weight load. The results show that the deformation deviation of the parts has the trend of transmission and accumulation under the load. The terminal deformation cumulative amount of the system is up to –0.0249 mm, which indicated that the influence of parts deformation on the mechanical system precision cannot be ignored.Originality/value This paper establishes a comprehensive 3D tolerance mathematical model, which comprehensively considers the effect of the dimensional tolerance, geometric tolerance and load deformation deviation. By this way, the assembly precision of mechanical system can be accurately predicted. [ABSTRACT FROM AUTHOR]

Purpose This paper aims to provide a comprehensive review of uncertainty quantification methods supported by evidence-based comparison studies. Uncertainties are widely encountered in engineering practice, arising from such diverse sources as heterogeneity of materials, variability in measurement, lack of data and ambiguity in knowledge. Academia and industries have long been researching for uncertainty quantification (UQ) methods to quantitatively account for the effects of various input uncertainties on the system response. Despite the rich literature of relevant research, UQ is not an easy subject for novice researchers/practitioners, where many different methods and techniques coexist with inconsistent input/output requirements and analysis schemes.Design/methodology/approach This confusing status significantly hampers the research progress and practical application of UQ methods in engineering. In the context of engineering analysis, the research efforts of UQ are most focused in two largely separate research fields: structural reliability analysis (SRA) and stochastic finite element method (SFEM). This paper provides a state-of-the-art review of SRA and SFEM, covering both technology and application aspects. Moreover, unlike standard survey papers that focus primarily on description and explanation, a thorough and rigorous comparative study is performed to test all UQ methods reviewed in the paper on a common set of reprehensive examples.Findings Over 20 uncertainty quantification methods in the fields of structural reliability analysis and stochastic finite element methods are reviewed and rigorously tested on carefully designed numerical examples. They include FORM/SORM, importance sampling, subset simulation, response surface method, surrogate methods, polynomial chaos expansion, perturbation method, stochastic collocation method, etc. The review and comparison tests comment and conclude not only on accuracy and efficiency of each method but also their applicability in different types of uncertainty propagation problems.Originality/value The research fields of structural reliability analysis and stochastic finite element methods have largely been developed separately, although both tackle uncertainty quantification in engineering problems. For the first time, all major uncertainty quantification methods in both fields are reviewed and rigorously tested on a common set of examples. Critical opinions and concluding remarks are drawn from the rigorous comparative study, providing objective evidence-based information for further research and practical applications. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this paper is to analyze the repair of a straight and angular crack in the structure using a piezoelectric material under thermo-mechanical loading by the extended finite element method (XFEM) approach. This provides a general and simple solution for the modeling of crack in the structure to analyze the repair. Design/methodology/approach: The extended finite element method is used to model crack geometry. The crack surface is modeled by Heaviside enrichment function while the crack front is modeled by branch enrichment functions. Findings: The effectiveness of the repair is measured in terms of stress intensity factor and J-integral. The critical voltage at which patch repair is most effective is evaluated and presented. Optimal patch shape, location of patch, adhesive thickness and adhesive modulus are obtained for effective repair under thermo-mechanical loading environment. Originality/value: The presented numerical modeling and simulation by the XFEM approach are of great benefit to analyze crack repair in two-dimensional and three-dimensional structures using piezoelectric patch material under thermo-mechanical loading. [ABSTRACT FROM AUTHOR]