Search results

The great magnitude differences between the integral tunnel and its structure details make it impossible to numerically model and analyze the global and local seismic behavior of large-scale shield tunnels using a unified spatial scale, even with the help of supercomputers. The paper aims to present a combined equivalent & multi-scale simulation method, by which the tunnel's major mechanical properties under seismic loads can be represented by the equivalent model, and the seismic responses of the interested details can be studied efficiently by the coupled multi-scale model.

The nominal orthotropic material constants of the equivalent tunnel model are inversely determined by fitting the modal characteristics of the equivalent model with the corresponding segmental lining model. The critical sections are selected by comprehensive analyzing of the integral compression/extension and bending loads in the equivalent lining under the seismic shaking and the coupled multi-scale model containing the details of interest is solved by the mixed time explicit integration algorithm.

The combined equivalent & multi-scale simulation method is an effective and efficient way for seismic analyses of large-scale tunnels. The response of each flexible joint is related to its polar location on the lining ring, and the mixed time integration method can speed-up the calculation process for hybrid FE model with great differences in element sizes.

The orthotropic equivalent assumption is, to the best of the authors’ knowledge, for the first time, used in the 3D simulation of the shield tunnel lining, representing the rigidity discrepancies caused by the structural property.

The multi-scale numerical simulation method, able to represent the complexity of the random structures and capture phase degradation, is an effective way to investigate the long-term behavior of concrete in service and bridges the gap between research on the material and on the structural level. However, the combined chemical-physical deterioration mechanisms of concrete remain a challenging task. The purpose of this paper is to investigate the degradation mechanism of concrete at the waterline in cold regions induced by combined calcium leaching and frost damage.

With the help of the NIST’s three-dimensional (3D) hydration model and the random aggregate model, realistic 3D representative volume elements (RVEs) of concrete at the micro-, the meso-, and the macro-scales can be reconstructed. The boundary problem method is introduced to compute the homogenized mechanical properties for both sound and damaged RVEs. According to the damage characteristics, the staggering method including a random dissolution model and a thermo-mechanical coupling model is developed to simulate the synergy deterioration effects of interacted calcium leaching and frost attacks. The coupled damage procedure for the frost damage process is based on the hydraulic pressure theory and the ice lens growth theory considering the relationship between the frozen temperature and the radius of the capillary pore. Finally, regarding calcium leaching as the leading role in actual engineering, the numerical methodology for combined leaching and frost damage on concrete property is proposed using a successive multi-scale method.

On the basis of available experimental data, this methodology is employed to explore the deterioration process. The results agree with the experimental ones to some extent, chemical leaching leads to the nucleation of some micro-cracks (i.e. damage), and consequently, to the decrease of the frost resistance.

It is demonstrated that the multi-scale numerical methodology can capture potential aging and deterioration evolution processes, and can give an insight into the macroscopic property degradation of concrete under long-term aggressive conditions.

This paper aims to introduce a multiscale computational method for structural failure analysis with inheriting simulation of moving trans-scale boundary (MTB). This method is motivated from the error in domain bridging caused by cross-scale damage evolution, which is common in structural failure induced by damage accumulation.

Within the method, vulnerable regions with high stress level are described by continuum damage mechanics, while elastic structural theory is sufficient for the rest, dividing the structural model into two scale domains. The two domains are bridged to generate mixed dimensional finite element equation of the whole system. Inheriting simulation is developed to make the computation of MTB sustainable.

Numerical tests of a notched three-point bending beam and a steel frame show that this MTB method can improve efficiency and ensure accuracy while capturing the effect of material damage on deterioration of components and structure.

The proposed MTB method with inheriting simulation is an extension of multiscale simulation to structural failure analysis. Most importantly, it can deal with cross-scale damage evolution and improve computation efficiency significantly.

Continuous fiber reinforced thermoplastic composites (CFRTPCs) with great mechanical properties and green recyclability have been widely used in aerospace, transportation, sports and leisure products, etc. However, the conventional molding technologies of CFRTPCs, with high cost and low efficiency, limit the property design and broad application of composite materials. The purpose of this paper is to study the effect of the 3D printing process on the integrated rapid manufacturing of CFRTPCs.

Tensile and flexural simulations and tests were performed on CFRTPCs. The effect of key process parameters on mechanical properties and molding qualities was evaluated individually and mutually to optimize the printing process. The micro morphologies of tensile and flexural breakages of the printed CFRTPCs were observed and analyzed to study the failure mechanism.

The results proved that the suitable process parameters for great printing qualities and mechanical properties included the glass hot bed with the microporous and solid glue coatings at 60°C and the nozzle temperature at 295°C. The best parameters of the nozzle temperature, layer thickness, feed rate and printing speed for the best elastic modulus and tensile strength were 285°C, 0.5 mm, 6.5r/min and 500 mm/min, respectively, whereas those for the smallest sectional porosity were 305°C, 0.6 mm, 5.5r/min and 550 mm/min, respectively.

This work promises a significant contribution to the improvement of the printing quality and mechanical properties of 3D printed CFRTPCs parts by the optimization of 3D printing processes.

The purpose of this research is to establish a robust numerical framework for the calibration of macroscopic constitutive parameters, based on the analysis of polycrystalline RVEs with computational homogenisation.

This framework is composed of four building-blocks: (1) the multi-scale model, consisting of polycrystalline RVEs, where the grains are modelled with anisotropic crystal plasticity, and computational homogenisation to link the scales, (2) a set of loading cases to generate the reference responses, (3) the von Mises elasto-plastic model to be calibrated, and (4) the optimisation algorithms to solve the inverse identification problem. Several optimisation algorithms are assessed through a reference identification problem. Thereafter, different calibration strategies are tested. The accuracy of the calibrated models is evaluated by comparing their results against an FE2 model and experimental data.

In the initial tests, the LIPO optimiser performs the best. Good results accuracy is obtained with the calibrated constitutive models. The computing time needed by the FE2 simulations is 5 orders of magnitude larger, compared to the standard macroscopic simulations, demonstrating how this framework is suitable to obtain efficient micro-mechanics-informed constitutive models.

This contribution proposes a numerical framework, based on FE2 and macro-scale single element simulations, where the calibration of constitutive laws is informed by multi-scale analysis. The most efficient combination of optimisation algorithm and definition of the objective function is studied, and the robustness of the proposed approach is demonstrated by validation with both numerical and experimental data.

Hybrid simulation, which is a general technique for obtaining the seismic response of an entire structure, is an improvement of the traditional seismic test technique. In order to improve the analysis accuracy of the numerical substructure in hybrid simulation, the purpose of this paper is to propose an innovative hybrid simulation technique. The technique combines the multi-scale finite element (MFE) analysis method and hybrid simulation method with the objective of achieving the balance between the accuracy and efficiency for the numerical substructure simulation.

To achieve this goal, a hybrid simulation system is established based on the MTS servo control system to develop a hybrid analysis model using an MFE model. Moreover, in order to verify the efficiency of the technique, the hybrid simulation of a three-storey benchmark structure is conducted. In this simulation, a ductile column—represented by a half-scale scale specimen—is selected as the experimental element, meanwhile the rest of the frame is modelled as microscopic and macroscopic elements in the Abaqus software simultaneously. Finally, to demonstrate the stability and accuracy of the proposed technique, the seismic response of the target structure obtained via hybrid simulation using the MFE model is compared with that of the numerical simulation.

First, the use of the hybrid simulation with the MFE model yields results similar to those obtained by the fine finite element (FE) model using solid elements without adding excessive computing burden, thus advancing the application of the hybrid simulation in large complex structures. Moreover, the proposed hybrid simulation is found to be more versatile in structural seismic analysis than other techniques. Second, the hybrid simulation system developed in this paper can perform hybrid simulation with the MFE model as well as handle the integration and coupling of the experimental elements with the numerical substructure, which consists of the macro- and micro-level elements. Third, conducting the hybrid simulation by applying earthquake motion to simulate seismic structural behaviour is feasible by using Abaqus to model the numerical substructure and harmonise the boundary connections between three different scale elements.

In terms of the implementation of the hybrid simulation with the MFE model, this work is helpful to advance the hybrid simulation method in the structural experiment field. Nevertheless, there is still a need to refine and enhance the current technique, especially when the hybrid simulation is used in real complex engineering structures, having numerous micro-level elements. A large number of these elements may render the relevant hybrid simulations unattainable because the time consumed in the numeral calculations can become excessive, making the testing of the loading system almost difficult to run smoothly.

The MFE model is implemented in hybrid simulation, enabling to overcome the problems related to the testing accuracy caused by the numerical substructure simplifications using only macro-level elements.

This paper is the first to recognise the advantage of the MFE analysis method in hybrid simulation and propose an innovative hybrid simulation technique, combining the MFE analysis method with hybrid simulation method to strike a delicate balance between the accuracy and efficiency of the numerical substructure simulation in hybrid simulation. With the help of the coordinated analysis of FEs at different scales, not only the accuracy and reliability of the overall seismic analysis of the structure is improved, but the computational cost can be restrained to ensure the efficiency of hybrid simulation.

Depending on an experimental approach to find optimal parameters for producing fully dense (relative density > 99%) Inconel 718 (IN718) components in the selective laser melting (SLM) process is expensive and offers no guarantee of success. Accordingly, this study aims to propose a multi-scale simulation framework to guide the choice of processing parameters in a more pragmatic manner.

In the proposed approach, a powder layer, ray tracing and heat transfer simulation models are used to calculate the melt pool dimensions and evaporation volume corresponding to a small number of laser power and scanning speed conditions within the input design space. A layer-heating model is then used to determine the inter-layer idle time required to maximize the temperature convergence rate of the solidified layer beneath the power bed. The simulation results are used to train surrogate models to construct SLM process maps for 3,600 pairs of the laser power and scanning speed within the input design space given three different values of the underlying solidified layer temperature (i.e., 353 K, 673 K and 873 K). The ideal selection of laser power and scanning speed of each process map is chosen based on four quality-related criteria listed as follows: without the appearance of key-hole melting; an evaporation volume less than the volume of the d90 powder particles; ensuring the stability of single scan tracks; and avoiding a weak contact between the melt pool and substrate. Finally, the optimal laser power and scanning speed parameters for the SLM process are determined by superimposing the optimal regions of the individual process maps.

The feasibility of the proposed approach is demonstrated by fabricating IN718 test specimens using the optimal processing conditions identified by the simulation framework. It is shown that the maximum density of the fabricated parts is 99.94%, while the average density is 99.88% and the standard deviation is less than 0.05%.

The present study proposed a multi-scale simulation model which can efficiently predict the optimal processing conditions for producing fully dense components in the SLM process. If the geometry of the three-dimensional printed part is changed or the machine and powder material is altered, users can use the proposed method for predicting the processing conditions that can produce the high-density part.

The flow phenomenon of particle image velocimetry has revealed the transition process of the complex multi-scale vortex between the boundary layer and mainstream region. Nonetheless, present computational fluid dynamics methods inadequately distinguish the discernable flows in detail. A multi-physical field coupling model, which was applied in rotor-stator fluid machinery (Umavathi, 2015; Syawitri et al., 2020), was put forward to ensure the identification of multi-scale vortexes and the improvement of performance prediction in torque converter.

A newly-developed multi-physical field simulation framework that coupled the scale-resolving simulation method with a dynamic modified viscosity coefficient was proposed to comparatively investigate the influence of energy exchange on thermal and flow characteristics and the description of the flow field in detail.

Regardless of whether quantitative or qualitative, its description ability on turbulence statistics, pressure-streamline, vortex structure and eddy viscosity ratio were visually experimentally and numerically analyzed. The results revealed that the modification of transmission medium viscous can identify flows more exactly between the viscous sublayer and outer boundary layer. Compared with RANS and large eddy simulation, a stress-blended eddy simulation model with a dynamic modified viscosity coefficient, which was further used to achieve blending on the stress level, can effectively solve the calculating problem of the transition region between the near-wall boundary layer and mainstream region.

This indeed provides an excellent description of the transient flow field and vortex structure in different physical flow states. Furthermore, the experimental data has proven that the maximum error of the external performance prediction was less than 4%.

An improved model was applied to simulate and analyze the flow mechanism through the evolution of vortex structures in a working chamber, to deepen the designer with a fundamental understanding on how to reduce flow losses and flow non-uniformity in manufacturing.

The purpose of this paper is to develop an empiricism free, first principle‐based model to simulate fluid flow and heat transfer through porous media.

Conventional approaches to the problem are reviewed. A multi‐scale approach that makes use of the sample simulations at the individual pore levels is employed. The effect of porous structures on the global fluid flow is accounted for via local volume averaged governing equations, while the closure terms are accounted for via averaging flow characteristics around the pores.

The performance of the model has been tested for an isothermal flow case. Good agreement with experimental data were achieved. Both the permeability and Ergun coefficient are shown to be flow properties as opposed to the empirical approach which typically results in constant values of these parameters independent of the flow conditions. Hence, the present multi‐scale approach is more versatile and can account for the possible changes in flow characteristics.

Further validation including non‐isothermal cases is necessary. Current scope of the model is limited to incompressible flows. The methodology can accommodate extension to compressible flows.

This paper proposes a method that eliminates the dependence of the numerical porous media simulations on empirical data. Although the model increases the fidelity of the simulations, it is still computationally affordable due to the use of a multi‐scale methodology.

The purpose of this paper is to capture the actual structural behavior of the longest timber footbridge in Spain by means of a multi-scale model updating approach in conjunction with ambient vibration tests.

In a first stage, a numerical pre-test analysis of the full bridge is performed, using standard beam-type finite elements with isotropic material properties. This approach offers a first structural model in which optimal sensor placement (OSP) methodologies are applied to improve the system identification process. In particular, the effective independence (EFI) method is used to determine the optimal locations of a set of sensors. Ambient vibration tests are conducted to determine experimentally the modal characteristics of the structure. The identified modal parameters are compared with those values obtained from this preliminary model. To improve the accuracy of the numerical predictions, the material response is modeled by means of a homogenization-based multi-scale computational approach. In a second stage, the structure is modeled by means of three-dimensional solid elements with the above material definition, capturing realistically the full orthotropic mechanical properties of wood. A genetic algorithm (GA) technique is adopted to calibrate the micromechanical parameters which are either not well-known or susceptible to considerable variations when measured experimentally.

An overall good agreement is found between the results of the updated numerical simulations and the corresponding experimental measurements. The longitudinal and transverse Young's moduli, sliding and rolling shear moduli, density and natural frequencies are computed by the present approach. The obtained results reveal the potential predictive capabilities of the present GA/multi-scale/experimental approach to capture accurately the actual behavior of complex materials and structures.

The uniqueness and importance of this structure leads to an intensive study of its structural behavior. Ambient vibration tests are carried out under environmental excitation. Extraction of modal parameters is obtained from output-only experimental data. The EFI methodology is applied for the OSP on a large-scale structure. Information coming from several length scales, from sub-micrometer dimensions to macroscopic scales, is included in the material definition. The strong differences found between the stiffness along the longitudinal and transverse directions of wood lumbers are incorporated in the structural model. A multi-scale model updating approach is carried out by means of a GA technique to calibrate the micromechanical parameters which are either not well-known or susceptible to considerable variations when measured experimentally.