Concurrent multi-scale modeling of civil infrastructures for analyses on structural deterioration—Part I : Modeling methodology and strategy

This paper aims to develop the methodology and strategy for concurrent finite element modeling of civil infrastructures at the different scale levels for the purposes of analyses of structural deteriorating. The modeling strategy and method were investigated to develop the concurrent multi-scale model of structural behavior (CMSM-of-SB) in which the global structural behavior and nonlinear damage features of local details in a large complicated structure could be concurrently analyzed in order to meet the needs of structural-state evaluation as well as structural deteriorating. In the proposed method, the “large-scale” modeling is adopted for the global structure with linear responses between stress and strain and the “small-scale” modeling is available for nonlinear damage analyses of the local welded details. A longitudinal truss in steel bridge decks was selected as a case to study how a CMSM-of-SB was developed. The reduced-scale specimen of the longitudinal truss was studied in the laboratory to measure its dynamic and static behavior in global truss and local welded details, while the multi-scale models using constraint equations and substructuring were developed for numerical simulation. The comparison of dynamic and static response between the calculated results by different models indicated that the proposed multi-scale model was found to be the most efficient and accurate. The verification of the model with results from the tested truss under the specific loading showed that, responses at the material scale in the vicinity of local details as well as structural global behaviors could be obtained and fit well with the measured results. The proposed concurrent multi-scale modeling strategy and implementation procedures were applied to Runyang cable-stayed bridge (RYCB) and the CMSM-of-SB of the bridge deck system was accordingly constructed as a practical application.

Concurrent multi-scale modeling of civil infrastructures for analyses on structural deteriorating—Part II : Model updating and verification

This paper is a continuation of the paper titled “Concurrent multi-scale modeling of civil infrastructure for analyses on structural deteriorating—Part I: Modeling methodology and strategy” with the emphasis on model updating and verification for the developed concurrent multi-scale model. The sensitivity-based parameter updating method was applied and some important issues such as selection of reference data and model parameters, and model updating procedures on the multi-scale model were investigated based on the sensitivity analysis of the selected model parameters. The experimental modal data as well as static response in terms of component nominal stresses and hot-spot stresses at the concerned locations were used for dynamic response- and static response-oriented model updating, respectively. The updated multi-scale model was further verified to act as the baseline model which is assumed to be finite-element model closest to the real situation of the structure available for the subsequent arbitrary numerical simulation. The comparison of dynamic and static responses between the calculated results by the final model and measured data indicated the updating and verification methods applied in this paper are reliable and accurate for the multi-scale model of frame-like structure. The general procedures of multi-scale model updating and verification were finally proposed for nonlinear physical-based modeling of large civil infrastructure, and it was applied to the model verification of a long-span bridge as an actual engineering practice of the proposed procedures.

Mixed-dimensional consistent coupling by multi-point constraint equations for efficient multi-scale modeling

Finding an appropriate linking method to connect different dimensional element types in a single finite element model is a key issue in the multi-scale modeling. This paper presents a mixed dimensional coupling method using multi-point constraint equations derived by equating the work done on either side of interface connecting beam elements and shell elements for constructing a finite element multiscale model. A typical steel truss frame structure is selected as case example and the reduced scale specimen of this truss section is then studied in the laboratory to measure its dynamic and static behavior in global truss and local welded details while the different analytical models are developed for numerical simulation. Comparison of dynamic and static response of the calculated results among different numerical models as well as the good agreement with those from experimental results indicates that the proposed multi-scale model is efficient and accurate.

Accommodating curvature in a highly ordered functionalized metal oxide nanofiber: synthesis, characterization and multi-scale modeling of layered nanosheets

A key element in the rational design of hybrid organic-inorganic nanostructures, is control of surfactant packing and adsorption onto the inorganic phase in crystal growth and assembly. In layered single crystal nanofibers and bilayered 2D nanosheets of vanadium oxide, we show how the chemisorption of preferred densities of surfactant molecules can direct formation of ordered, curved layers. The atom-scale features of the structures are described using molecular dynamics simulations that quantify surfactant packing effects and confirm the preference for a density of 5 dodecanethiol molecules per 8 vanadium attachment sites in the synthesised structures. This assembly maintains a remarkably well ordered interlayer spacing, even when curved. The assemblies of interdigitated organic bilayers on V2O5 are shown to be sufficiently flexible to tolerate curvature while maintaining a constant interlayer distance without rupture, delamination or cleavage. The accommodation of curvature and invariant structural integrity points to a beneficial role for oxide-directed organic film packing effects in layered architectures such as stacked nanofibers and hybrid 2D nanosheet systems.

Fracture of DP590 steel : a multi-scale modeling approach

Advanced high strength steel sheets are one of the higher strength advance material developed by the steel industry for automotive bodies. One of the categories of this advanced high strength steel is Dual Phase (DP) steel. This steel consists of a two phase microstructure where soft and hard phase acts together to offer a high strength composite effect. The combination of high strength and ductility exhibited by these sheets allows the design and manufacture of complex parts. However, during forming certain grades of DP steel sudden cracking can occur without any intimation of necking. This abnormal forming behavior is difficult to accurately predict because most classical modelling approaches are not designed for such micro-structurally heterogeneous materials. These modelling approaches are generally based on an average representation of the material behaviour in a continuum mechanics formulation. This works for materials that are homogenous, or at least could be assumed to be homogenous at scales lower than the naked eye can see. However, for a material like advanced high strength steel, the microstructure plays a significant role in dictating the mechanical behavior at the macro-scale. This paper studies the forming and fracture behavior through multi-scale modeling of DPO590 steel. It is found that the sufficient accuracy can be achieved from multi-scale modeling when comparing with experiments.

Multi-Scale Modeling for Image Analysis of Brain Tumor Studies

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

Updating and verification for multi-scale finite element model of structure behavior

The method on concurrent multi-scale model of structural behavior (CMSM-of-SB) for the purpose of structural health monitoring including model updating and validating has been studied. The detailed process of model updating and validating is discussed in terms of reduced scale specimen of the steel box girder in longitudinal stiffening truss of a long span bridge. Firstly, some influence factors affecting the accuracy of the CMSM-of-SB including the boundary restraint regidity, the geometry and material parameters on the toe of the weld and its neighbor are analyzed using sensitivity method. Then, sensitivity-based model updating technology is adopted to update the developed CMSM-of-SB and model verification is carried out through calculating and comparing stresses on different locations under various loading from dynamic characteristic and static response. It can be concluded that the CMSM-of-SB based on the substructure method is valid.

Multi-scale rheological model for discontinuous phenomena in materials under deformation conditions

A study of possibilities given by the developed Cellular Automata–Finite Element (CAFE) multi-scale model for prediction of the initiation and propagation of micro-shear bands and shear bands in metallic materials subjected to plastic deformation is described in the paper. Particular emphasis in defining the criterion for initiation of micro-shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of those phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi-scale model of strain localization. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the paper. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is presented in the paper as well. In this approach remeshing becomes possible and mesh distortion, which limits application of the CAFE method to simple deformation processes, is eliminated.

Development of the multi-scale analysis model to simulate strain localization occurring during material processing

Abstract A detailed description of possibilities given by the developed Cellular Automata—Finite Element (CAFE) multi scale model for prediction of the initiation and propagation of micro shear bands and shear bands in metallic materials subjected to plastic deformation is presented in the work. Particular emphasis in defining the criterion for initiation of micro shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of these phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi scale model of strain localization phenomena. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the work. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is a subject of this work as well. In the developed model remeshing becomes possible and difficulties limiting application of the CAFE method to simple deformation processes are solved. Obtained results of numerical simulaA detailed description of possibilities given by the developed Cellular Automata—Finite Element (CAFE) multi scale model for prediction of the initiation and propagation of micro shear bands and shear bands in metallic materials subjected to plastic deformation is presented in the work. Particular emphasis in defining the criterion for initiation of micro shear and shear bands, as well as in defining the transition rules for the cellular automata, is put on accounting for the physical aspects of these phenomena occurring in two different scales in the material. The proposed approach led to the creation of the real multi scale model of strain localization phenomena. This model predicts material behavior in various thermo-mechanical processes. Selected examples of applications of the developed model to simulations of metal forming processes, which involve strain localization, are presented in the work. An approach based on the Smoothed Particle Hydrodynamic, which allows to overcome difficulties with remeshing in the traditional CAFE method, is a subject of this work as well. In the developed model remeshing becomes possible and difficulties limiting application of the CAFE method to simple deformation processes are solved. Obtained results of numerical simulations are compared with the experimental results of cold rolling process to show good predicative capabilities of the developed model.tions are compared with the experimental results of cold rolling process to show good predicative capabilities of the developed model.

First steps towards modeling a multi-scale Earth system

Recent advances in computational geodynamics are applied to explore the link between Earth’s heat, its chemistry and its mechanical behavior. Computational thermal-mechanical solutions are now allowing us to understand Earth patterns by solving the basic physics of heat transfer. This approach is currently used to solve basic convection patterns of terrestrial planets. Applying the same methodology to smaller scales delivers promising similarities between observed and predicted structures which are often the site of mineral deposits. The new approach involves a fully coupled solution to the energy, momentum and continuity equations of the system at all scales, allowing the prediction of fractures, shear zones and other typical geological patterns out of a randomly perturbed initial state. The results of this approach are linking a global geodynamic mechanical framework over regional-scale mineral deposits down to the underlying micro-scale processes. Ongoing work includes the challenge of incorporating chemistry into the formulation.

Probabilistic downscaling of remote sensing data with applications for multi-scale biogeochemical flux modeling

Upscaling ecological information to larger scales in space and downscaling remote sensing observations or model simulations to finer scales remain grand challenges in Earth system science. Downscaling often involves inferring subgrid information from coarse-scale data, and such ill-posed problems are classically addressed using regularization. Here, we apply two-dimensional Tikhonov Regularization (2DTR) to simulate subgrid surface patterns for ecological applications. Specifically, we test the ability of 2DTR to simulate the spatial statistics of high-resolution (4 m) remote sensing observations of the normalized difference vegetation index (NDVI) in a tundra landscape. We find that the 2DTR approach as applied here can capture the major mode of spatial variability of the high-resolution information, but not multiple modes of spatial variability, and that the Lagrange multiplier (γ) used to impose the condition of smoothness across space is related to the range of the experimental semivariogram. We used observed and 2DTR-simulated maps of NDVI to estimate landscape-level leaf area index (LAI) and gross primary productivity (GPP). NDVI maps simulated using a γ value that approximates the range of observed NDVI result in a landscape-level GPP estimate that differs by ca 2% from those created using observed NDVI. Following findings that GPP per unit LAI is lower near vegetation patch edges, we simulated vegetation patch edges using multiple approaches and found that simulated GPP declined by up to 12% as a result. 2DTR can generate random landscapes rapidly and can be applied to disaggregate ecological information and compare of spatial observations against simulated landscapes.

Predicting active material utilisation in LiFePO4 electrodes using a multi-scale mathematical model

A mathematical model is developed to simulate the discharge of a LiFePO4 cathode. This model contains 3 size scales, which match with experimental observations present in the literature on the multi-scale nature of LiFePO4 material. A shrinking-core is used on the smallest scale to represent the phase-transition of LiFePO4 during discharge. The model is then validated against existing experimental data and this validated model is then used to investigate parameters that influence active material utilisation. Specifically the size and composition of agglomerates of LiFePO4 crystals is discussed, and we investigate and quantify the relative effects that the ionic and electronic conductivities within the oxide have on oxide utilisation. We find that agglomerates of crystals can be tolerated under low discharge rates. The role of the electrolyte in limiting (cathodic) discharge is also discussed, and we show that electrolyte transport does limit performance at high discharge rates, confirming the conclusions of recent literature.

Time dependent deformations in concrete: A multi-scale approach

Estimation of creep and shrinkage are critical in order to compute loss of prestress with time in order to compute leak tightness and assess safety margins available in containment structures of nuclear power plants. Short-term creep and shrinkage experiments have been conducted using in-house test facilities developed specifically for the present research program on 35 and 45 MPa normal concrete and 25 MPa heavy density concrete. The extensive experimental program for creep, has cylinders subject to sustained levels of load typically for several days duration (till negligible strain increase with time is observed in the creep specimen), to provide the total creep strain versus time curves for the two normal density concrete grades and one heavy density concrete grade at different load levels, different ages at loading, and at different relative humidity’s. Shrinkage studies on prism specimen for concrete of the same mix grades are also being studied. In the first instance, creep and shrinkage prediction models reported in the literature has been used to predict the creep and shrinkage levels in subsequent experimental data with acceptable accuracy. While macro-scale short experiments and analytical model development to estimate time dependent deformation under sustained loads over long term, accounting for the composite rheology through the influence of parameters such as the characteristic strength, age of concrete at loading, relative humidity, temperature, mix proportion (cement: fine aggregate: coarse aggregate: water) and volume to surface ratio and the associated uncertainties in these variables form one part of the study, it is widely believed that strength, early age rheology, creep and shrinkage are affected by the material properties at the nano-scale that are not well established. In order to understand and improve cement and concrete properties, investigation of the nanostructure of the composite and how it relates to the local mechanical properties is being undertaken. While results of creep and shrinkage obtained at macro-scale and their predictions through rheological modeling are satisfactory, the nano and micro indenting experimental and analytical studies are presently underway. Computational mechanics based models for creep and shrinkage in concrete must necessarily account for numerous parameters that impact their short and long term response. A Kelvin type model with several elements representing the influence of various factors that impact the behaviour is under development. The immediate short term deformation (elastic response), effects of relative humidity and temperature, volume to surface ratio, water cement ratio and aggregate cement ratio, load levels and age of concrete at loading are parameters accounted for in this model. Inputs to this model, such as the pore structure and mechanical properties at micro/nano scale have been taken from scanning electron microscopy and micro/nano-indenting of the sample specimen.