Multi-scale Modeling and Finite Element Analyses of Thermal Conductivity of 3D C/SiC Composites Fabricating by Flexible-Oriented Woven Process

Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.

Multi-scale Simulation Method with Coupled Finite/Discrete Element Model and Its Application

The existing research on continuous structure is usually analyzed with finite element method (FEM) and granular medium with discrete element method (DEM), but there are few researches on the coupling interaction between continuous structure and discrete medium. To the issue of this coupling interaction, a multi-scale simulation method with coupled finite/discrete element model is put forward, in their respective domains of discrete and finite elements, the nodes follow force law and motion law of their own method, and on the their interaction interface, the touch type between discrete and finite elements is distinguished as two types: full touch and partial touch, the interaction force between them is calculated with linear elastic model. For full touch, the contact force is proportional to the overlap distance between discrete element and finite element patch. For partial touch, first the finite element patch is extended on all sides indefinitely to be a complete plane, the full contact force can be obtained with the touch type between discrete element and plane being viewed as full touch, then the full overlap area between them and the actual overlap area between discrete element and finite element patch are computed, the actual contact force is obtained by scaling the full contact force with a factor which is determined by the ratio of the actual overlap area to the full overlap area. The contact force is equivalent to the finite element nodes and the force and displacement on the nodes can be computed, so the ideal simulation results can be got. This method has been used to simulate the cutter disk of the earth pressure balance shield machine (EPBSM) made in North Heavy Industry (NHI) with its excavation diameter of 6.28 m cutting and digging the sandy clay layer. The simulation results show that as the gradual increase of excavating depth of the cutter head, the maximum stress occurs at the roots of cutters on the cutter head, while for the soil, the largest stress is distributed at the region which directly contacted with the cutters. The proposed research can provide good solutions for correct design and installation of cutters, and it is necessary to design mounting bracket to fix cutters on cutter head.

Fatigue Performance Evaluation for Welded Details in Orthotropic Steel Deck Bridges Using Multi-Scale Finite Element Method

In order to study the fatigue properties of rib-to-deck welded connection and rib-to-rib welded connection in orthotropic steel bridge decks,a multi-scale finite element model was set up to analyze the stress distribution characteristics and the load test was conducted on the Taizhou Yangtze River Bridge.Comparing the vehicle test results with the muli-scale finite element model results to verify the accuracy of the finite element simulation for the stress response of two welded details.The results indicated that The stress at the rib-to-deck welded connection and the rib-to-rib welded connection are the bending stress and the membrane stress,respectively;the stress response of the two welded connection has strong local characteristics;the lateral stress influence line of the two welded connection is relatively short and the length of the lateral stress influence line is greatly affected by the longitudinal ribs;increasing the thickness of the roof and longitudinal ribs can reduce the stress response and improve the stress performance of the heavy lanes.For the two welded details,the fatigue damage increment of the ordinary lane is greater than the heavy lane.The thickened roof and longitudinal ribs at the position of the heavy lane still cannot balance the fatigue damage caused by the heavy truck.Therefore,it is necessary to strictly control the fatigue effect of overloaded vehicles on steel box girders.

Numerical failure analysis of a continuous reinforced concrete bridge under strong earthquakes using multi-scale models

Previous failure analyses of bridges typically focus on substructure failure or superstructure failure separately. However, in an actual bridge, the seismic induced substructure failure and superstructure failure may influence each other. Moreover, previous studies typically use simplified models to analyze the bridge failure; however, there are inherent defects in the calculation accuracy compared with using a detailed three-dimensional （3D） finite element （FE） model. Conversely, a detailed 3D FE model requires more computational costs, and a proper erosion criterion of the 3D elements is necessary. In this paper, a multi-scale FE model, including a corresponding erosion criterion, is proposed and validated that can significantly reduce computational costs with high precision by modelling a pseudo-dynamic test of an reinforced concrete （RC） pier. Numerical simulations of the seismic failures of a continuous RC bridge based on the multi-scale FE modeling method using LS-DYNA are performed. The nonlinear properties of the bridge, various connection strengths and bidirectional excitations are considered. The numerical results demonstrate that the failure of the connections will induce large pounding responses of the girders. The nonlinear deformation of the piers will aggravate the pounding damages. Furthermore, bidirectional earthquakes will induce eccentric poundingsto the girders and different failure modes to the adjacent piers.

Multi-scale finite element analysis of chloride diffusion in concrete incorporating paste/aggregate ITZs

In this paper,we propose a concurrent multi-scale finite element(FE) model coupling equations of the degree of freedoms of meso-scale model of ITZs and macroscopic model of bulk pastes.The multi-scale model is subsequently implemented and integrated into ABAQUS resulting in easy application to complex concrete structures.A few benchmark numerical examples are performed to test both the accuracy and efficiency of the developed model in analyzing chloride diffusion in concrete.These examples clearly demonstrate that high diffusivity of ITZs,primarily because of its porous microstructure,tends to accelerate chloride penetration along concentration gradient.The proposed model provides new guidelines for the durability analysis of concrete structures under adverse operating conditions.

Three-dimensional modelling of shear keys in concrete gravity dams using an advanced grillage method

Contraction joint shear keys are resilient features of gravity dams that can be considered to increase the sliding safety factors or minimise seismic residual sliding displacements,allowing costly remedial actions to be avoided.This paper presents a novel,robust,and computationally efficient three-dimensional(3D)modelling and simulation strategy of gravity dams,using a series of adjacent cantilever beam elements to represent individual monoliths.These monoliths are interconnected in the longitudinal direction by 3D no-tension link elements representing the lumped shear key stiffness contributions at a particular elevation.The objective is to assess the shear key internal force demands,including the axial force,shear,and moment demands.Shear key demand-capacity ratios can then be assessed with related multi-axial failure envelopes.The 3D link element stiffness coefficients were derived from a series of 3D finite element(FE)solid models with a detailed representation of geometrical features of multiple shear keys.The results from the proposed method based on advanced grillage analysis show strong agreement with reference solutions from 3D FE solid models,demonstrating high accuracy and performance of the proposed method.The application of the proposed advanced grillage method to a dam model with two monoliths clearly shows the advantage of the proposed method,in comparison to the classical approach used in practise.

Elastic Predictions of 3D Orthogonal Woven Composites Using Micro/meso-scale Repeated Unit Cell Models

This presentation predicts the elastic properties of three-dimensional(3D)orthogonal woven composite(3DOWC)by finite element analysis based on micro/meso repeated unit cell(RUC)models.First,the properties of fiber yarn are obtained by analysis on a micro-scale RUC model assuming fibers in a hexagonal distribution pattern in the polymer matrix.Then a full thickness meso-scale RUC model including weft yarns,warp yarns,Z-yarns and pure resin zones is established and full stiffness matrix of the 3DOWC including the in-plane and flexural constants are predicted.For thick 3DOWC with large number of weft,warp layers,an alternative analysis method is proposed in which an inner meso-RUC and a surface meso-RUC are established,respectively.Then the properties of 3DOWC are deduced based on laminate theory and properties of the inner and surface layers.The predicted results by the above two alternative methods are in good experimental agreement.

MULTI-SCALE METHODS FOR INVERSE MODELING IN 1-D MOS CAPACITOR

In this paper, we investigate multi-scale methods for the inverse modeling in 1-D Metal-Oxide-Silicon (MOS) capacitor. First, the mathematical model of the device is given and the numerical simulation for the forward problem of the model is implemented using finite element method with adaptive moving mesh. Then numerical analysis of these parameters in the model for the inverse problem is presented. Some matrix analysis tools are applied to explore the parameters' sensitivities. And third, the parameters are extracted using Levenberg-Marquardt optimization method. The essential difficulty arises from the effect of multi-scale physical difference of the parameters. We explore the relationship between the parameters' sensitivities and the sequence for optimization, which can seriously affect the final inverse modeling results. An optimal sequence can efficiently overcome the multi-scale problem of these parameters. Numerical experiments show the efficiency of the proposed methods.

双头单螺杆泵定子衬套力学行为分析

为研究双头单螺杆泵定子衬套的受力状态及内轮廓线的变形规律,分析了内压、压差、材料非线性、橡胶泊松比等多种因素对定子衬套应力、应变及内轮廓线变形规律的影响,并拟合出双头单螺杆泵定子衬套在均匀内压作用下内轮廓线位移的方程式及非均匀内压下内轮廓线最大位移公式。分析结果表明:均匀压力下,衬套应力、应变呈周期性变化,内轮廓线位移呈正弦规律变化,位移与内压成线性关系;压差作用下,衬套应力、应变最大值向压差作用处转移,内轮廓线变形规律发生改变,位移及其变形幅度均增大;衬套橡胶泊松比微小的变化即会对衬套变形产生非常大的影响。因此,实际生产中应注意注胶工艺及泊松比的变化。研究结果为定子、转子过盈量的选取及线型的优化设计提供了理论基础,对提高双头单螺杆泵工作效率和使用寿命具有重要工程实际意义。

Finite Element Prediction of the Thermal Conductivity of GNP/Al Composites

A 3 D multi-scale finite element model was developed to predict the effective thermal conductivity of graphene nanoplatelet(GNP)/Al composites.The factors influencing the effective thermal conductivity of the GNP/Al composites were investigated,including the orientation,shape,aspect ratio,configuration and volume fraction of GNPs.The results show that GNPs shape has a little influence on the thermal conductivity of GNP/Al composites,and composites with elliptic GNPs have the highest thermal conductivity.In addition,with increasing the aspect ratio of GNPs,the thermal conductivity of GNP/Al composites increases and finally tends to be stable.The GNPs configuration strongly influences the thermal conductivity of GNP/Al composites,and the thermal conductivity of the composites with layered GNPs is the highest among the five configurations.The effective thermal conductivity is sensitive to volume fraction of GNPs.Ideally,when the volume fraction of layered GNPs reaches 1.54%,the thermal conductivity of GNP/Al composites is as high as 400 W/m K.The findings of this study could provide a good theoretical basis for designing high thermal conductivity GNP/Al composites.

Back-Analysis of the Response of Shield Tunneling by 3D Finite Element Method

This paper presents a numerical back-analysis of the response of a shield tunnel during construction. An important issue in the construction of shallow tunnels, especially in soft ground conditions, is the surface settlement caused by shield tunneling. The tunnel test system with 10 m length, 7 m width and 6.7 m height, which was completed in China in 2009, is a research shield tunnel system. Using shield tunneling technique known as earth pressure balance (EPB) and slurry shield method, it could be excavated in a region consisting of original soft soils, such as silty clay, and different types of underlain soft soils. Based on the test results, the real-life tunnel response can be analyzed by back-analysis technique. The back-analysis technique is adapted to the three-dimensional finite element method (FEM). Parameter analyses are calibrated to study the behavior of the multi-scale diameter tunnel under various conditions. The suggested multi-scale model results show a well agreement between the prediction and the measurement.

Research on failure criterion of composite based on unified macro- and micro-mechanical model

A new unified macro- and micro-mechanics failure analysis method for composite structures was developed in order to take the effects of composite micro structure into consideration. In this method, the macro stress distribution of composite structure was calculated by commercial finite element analysis software. According to the macro stress distribution, the damage point was searched and the micro-stress distribution was calculated by reformulated finite-volume direct averaging micromechanics （FVDAM）, which was a multi-scale finite element method for composite. The micro structure failure modes were estimated with the failure strength of constituents. A unidirectional composite plate with a circular hole in the center under two kinds of loads was analyzed with the traditional macro-mechanical failure analysis method and the unified macro- and micro-mechanics failure analysis method. The results obtained by the two methods are consistent, which show this new method＇s accuracy and efficiency.