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 pr...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.展开更多
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 betwe...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.展开更多
A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale met...A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.展开更多
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 distribu...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.展开更多
In order to study the dynamic responses in the microstructures of the pavement structure, the multi-scale modeling subjected to moving load is analyzed using the discrete element method (DEM). The macro-scale discre...In order to study the dynamic responses in the microstructures of the pavement structure, the multi-scale modeling subjected to moving load is analyzed using the discrete element method (DEM). The macro-scale discrete element model of the flexible pavement structure is established. The stress and strain at the bottom of the asphalt concrete layer under moving load are calculated. The DEM model is validated through comparison between DEM predictions and the results from the classical program. Based on the validated macro-scale DEM model, the distribution and the volumetric fraction of coarse aggregate, mastics and air voids at the bottom of the asphalt layer are modeled, and then the multi-scale model is constructed. The dynamic response in the microstructures of the multi-scale model are calculated and compared with the results from the macro model. The influence of mastic stiffness on the distribution of dynamic response in the microstructures is also analyzed. Results show that the average values and the variation coefficient of the tensile stress at the aggregate-mastic interface are far more than those within the mastics. The dynamic response including stress and strain distributes non-uniformly in both mastics and the interface. An increase in mastic stiffness tends to a uniform distribution of tensile stress in asphalt concrete.展开更多
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 ...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.展开更多
This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material ...This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material for the specimens consists of a single hole model(25% void/cell,16% void/cell and 10% void/cell)and a four-hole model(25%void/cell).Using a representative volume element(RVE),we try to produce the equivalent homogenized properties and work on a homogeneous specimen for the study of fretting fatigue.Next,the fretting fatigue contact problem is performed for 3 new cases of models that consist of a homogeneous and a heterogeneous part(single hole cell)in the contact area.The aim is to analyze the normal and shear stresses of these models and compare them with the results of the corresponding heterogeneous models based on the Direct Numerical Simulation(DNS)method.Finally,by comparing the computational time and%deviations,we draw conclusions about the reliability and effectiveness of the proposed method.展开更多
Watershed segmentation is sensitive to noises and irregular details within the image,which frequently leads to a serious over-segmentation Linear filtering before watershed segmentation can reduce over-segmentation to...Watershed segmentation is sensitive to noises and irregular details within the image,which frequently leads to a serious over-segmentation Linear filtering before watershed segmentation can reduce over-segmentation to some extent,however,it often causes the position offset of object contours.For the purpose of reducing over-segmentation to preserve the location of object contours,the watershed segmentation based on the hierarchical multi-scale modification of morphological gradient is proposed.Firstly,multi-scale morphological filtering was employed to smooth the original image.Then,the gradient image was divided into multi-levels by the volume of three-dimension topographic relief,where the lower gradient layers were further modifiedby morphological closing with larger-sized structuring-elements,and the higher layers with the smaller one.In this way,most local minimums caused by irregular details and noises can be removed,while region contour positions corresponding to the target area were largely preserved.Finally,morphological watershed algorithm was employed to implement segmentation on the modified gradient image.The experimental results show that the proposed method can greatly reduce the over-segmentation of the watershed and avoid the position offset of the object contours.展开更多
Working conditions of rolling bearings of wind turbine generators are complicated, and their vibration signals often show non-linear and non-stationary characteristics. In order to improve the efficiency of feature ex...Working conditions of rolling bearings of wind turbine generators are complicated, and their vibration signals often show non-linear and non-stationary characteristics. In order to improve the efficiency of feature extraction of wind turbine rolling bearings and to strengthen the feature information, a new structural element and an adaptive algorithm based on the peak energy are proposed,which are combined with spectral correlation analysis to form a fault diagnosis algorithm for wind turbine rolling bearings. The proposed method firstly addresses the problem of impulsive signal omissions that are prone to occur in the process of fault feature extraction of traditional structural elements and proposes a "W" structural element to capture more characteristic information. Then, the proposed method selects the scale of multi-scale mathematical morphology, aiming at the problem of multi-scale mathematical morphology scale selection and structural element expansion law. An adaptive algorithm based on peak energy is proposed to carry out morphological scale selection and structural element expansion by improving the computing efficiency and enhancing the feature extraction effect.Finally, the proposed method performs spectral correlation analysis in the frequency domain for an unknown signal of the extracted feature and identifies the fault based on the correlation coefficient. The method is verified by numerical examples using experimental rig bearing data and actual wind field acquisition data and compared with traditional triangular and flat structural elements. The experimental results show that the new structural elements can more effectively extract the pulses in the signal and reduce noise interference,and the fault-diagnosis algorithm can accurately identify the fault category and improve the reliability of the results.展开更多
Multi-scale finite element method is adopted to simulate wood compression behavior under axial and transverse loading. Representative volume elements (RVE) of wood microfibril and cell are proposed to analyze orthotro...Multi-scale finite element method is adopted to simulate wood compression behavior under axial and transverse loading. Representative volume elements (RVE) of wood microfibril and cell are proposed to analyze orthotropic mechanical behavior. Lignin, hemicellulose and crystalline-amorphous cellulose core of spruce are concerned in spruce nanoscale model. The equivalent elastic modulus and yield strength of the microfibril are gained by the RVE simulation. The anisotropism of the crystalline-amorphous cellulose core brings the microfibril buckling deformation during compression loading. The failure mechanism of the cell-wall under axial compression is related to the distribution of amorphous cellulose and crystalline cellulose. According to the spruce cell observation by scanning electron microscope, numerical model of spruce cell is established using simplified circular hole and regular hexagon arrangement respectively. Axial and transverse compression loadings are taken into account in the numerical simulations. It indicates that the compression stress-strain curves of the numerical simulation are consistent with the experimental results. The wood microstructure arrangement has an important effect on the stress plateau during compression process. Cell-wall buckling in axial compression induces the stress value drops rapidly. The wide stress plateau duration means wood is with large energy dissipation under a low stress level. The numerical results show that loading velocity affects greatly wood microstructure failure modes in axial loading. For low velocity axial compression, shear sliding is the main failure mode. For high velocity axial compression, wood occur fold and collapse. In transverse compression, wood deformation is gradual and uniform, which brings stable stress plateau.展开更多
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 subsequ...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.展开更多
A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by com...A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by combining the proposed multi-scale crystal plasticity finite element method with experimental measurements.The complex local micro region load was progressively extracted from the simulation results of a macro model and applied to an established three-dimensional multi-grain microstructure model.Subsequently,the evolution histories of the grain shape,size,and orientation inside the adiabatic shear band were quantitatively simulated.The results corresponded closely to the experimental results obtained via transmission electron microscopy and precession electron diffraction.Furthermore,by calculating the grain rotation and temperature rise inside the adiabatic shear band,the microstructural softening and thermal softening effects of typical heavily-deformed α grains were successfully decoupled.The results revealed that the microstructural softening stress was triggered and then stabilized(in general)at a relatively high value.This indicated that the mechanical strength was lowered mainly by the grain orientation evolution or dynamic recrystallization occurring during early plastic deformation.Subsequently,thermal softening increased linearly and became the main softening mechanism.Noticeably,in the final stage,the thermal softening stress accounted for 78.4% of the total softening stress due to the sharp temperature increase,which inevitably leads to the stress collapse and potential failure of the alloy.展开更多
Ceramic Matrix Composite (CMC) turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress flu...Ceramic Matrix Composite (CMC) turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation, Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration (CVI) pro- cessed SiC/SiC-CMC are predicted by a Representative Volume Element (RVE) model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.展开更多
The challenge of estimating mine-wide subsidence and linkages to seismicity over tabular deposits is addressed by a special finite element technique(dual node–dual mesh). Subsidence and mine-induced seismicity begins...The challenge of estimating mine-wide subsidence and linkages to seismicity over tabular deposits is addressed by a special finite element technique(dual node–dual mesh). Subsidence and mine-induced seismicity begins near the face when caving occurs and propagates to the surface as extraction reaches a critical extent. Thus, the challenge is to obtain details at the face at the meter scale and also at the surface over the whole mine at the kilometer scale. Interactions between old and new sections of a mine are automatically taken into account with this technique. The finite element method is well established technology based on fundamentals of physical laws, kinematics and material laws. With this technique, no empirical ‘‘scaling" or fitting computer output by input data ‘‘adjustment" to mine measurements is necessary. Capability is demonstrated for doing practical whole-mine subsidence analysis from first principles. Mine-induced seismicity is shown to correlate well with face advance and element failure.展开更多
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 ...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.展开更多
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 ...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.展开更多
The void evolution of large-section plastic mold steel during multi-directional forging(MDF)was investigated using multiscale analysis.To simulate the forging process of the plastic mold steel(SDP1 steel)and realize m...The void evolution of large-section plastic mold steel during multi-directional forging(MDF)was investigated using multiscale analysis.To simulate the forging process of the plastic mold steel(SDP1 steel)and realize micro-void reconstruction in a representative volume element(RVE),MDF experiment and void-characteristic evaluation of the SDP1 steel were carried out.Traditional upsetting and stretching forging(TUSF)and MDF were simulated to comparatively analyze the evolution of temperature,effective stress,and effective strain.By embedding RVE with a micro-void and using boundary condition by point tracking into the forging process,the single-void evolution in TUSF and MDF was studied.The effect of void orientation on single-void evolution was also investigated.The multi-scale analysis revealed the following results.(1)Compared with TUSF,MDF achieved a higher efficiency in void closure.(2)The closing efficiency of the void increased with the increase in angle h(the angle between the Z and long axes of the void).(3)The closing efficiency increased with the increase in the orientation angle during the forging process.On the basis of the important role of the main stress in each forging step on the void closure,an integral formula of the main stress was proposed.When main compressive-stress integration reached-0.4,the closed state of the void could be accurately determined.展开更多
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 yar...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.展开更多
Deals with a study which discussed the mechanical behaviour for subdivided periodic elastic structures of composite materials, from the viewpoint of macro- and meso-scale coupling. Multiscale asymptotic expansion and ...Deals with a study which discussed the mechanical behaviour for subdivided periodic elastic structures of composite materials, from the viewpoint of macro- and meso-scale coupling. Multiscale asymptotic expansion and truncation error estimates; Discussion on multiscale finite element method; Details of higher order difference quotients and total error estimates; Numerical experiments.展开更多
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 c...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.展开更多
基金Supported by Science Center for Gas Turbine Project of China (Grant No.P2022-B-IV-014-001)Frontier Leading Technology Basic Research Special Project of Jiangsu Province of China (Grant No.BK20212007)the BIT Research and Innovation Promoting Project of China (Grant No.2022YCXZ019)。
文摘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.
基金supported by National Basic Research Program of China (973 Program, Grant No. 2013CB035400)Science Fund for Creative Research Groups of NSFC of China (Grant No. 51221004)National Natural Science Foundation of China (Grant No. 51075357)
文摘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.
基金Supported by the National Natural Science Foundation of China(51105195,51075204)the Aeronautical Science Foundation of China(2011ZB52024)
文摘A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.
基金This research has been supported by the National Natural Science Foundation of China(Grant No.51778135)the National Key R&D Program Foundation of China(Grant No.2017YFC0806001)+2 种基金the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20160207)Aeronautical Science Foundation of China(Grant No.20130969010)the Fundamental Research Funds for the Central Universities and Postgraduate Research&Practice Innovation Program of Jiangsu Province,China(Grant Nos.KYCX18_0113 and KYLX16_0253).
文摘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.
基金The National Natural Science Foundation of China (No.51208178,51108157)China Postdoctoral Science Foundation (No.2012M520991)
文摘In order to study the dynamic responses in the microstructures of the pavement structure, the multi-scale modeling subjected to moving load is analyzed using the discrete element method (DEM). The macro-scale discrete element model of the flexible pavement structure is established. The stress and strain at the bottom of the asphalt concrete layer under moving load are calculated. The DEM model is validated through comparison between DEM predictions and the results from the classical program. Based on the validated macro-scale DEM model, the distribution and the volumetric fraction of coarse aggregate, mastics and air voids at the bottom of the asphalt layer are modeled, and then the multi-scale model is constructed. The dynamic response in the microstructures of the multi-scale model are calculated and compared with the results from the macro model. The influence of mastic stiffness on the distribution of dynamic response in the microstructures is also analyzed. Results show that the average values and the variation coefficient of the tensile stress at the aggregate-mastic interface are far more than those within the mastics. The dynamic response including stress and strain distributes non-uniformly in both mastics and the interface. An increase in mastic stiffness tends to a uniform distribution of tensile stress in asphalt concrete.
基金National Program on Key Basic Research Project of China(973) under Grant No.2011CB013603the National Natural Science Foundation of China under Grant Nos.51427901,91315301 and 51408410the Natural Science Foundation of Tianjin,China under Grant No.15JCQNJC07200
文摘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.
文摘This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material for the specimens consists of a single hole model(25% void/cell,16% void/cell and 10% void/cell)and a four-hole model(25%void/cell).Using a representative volume element(RVE),we try to produce the equivalent homogenized properties and work on a homogeneous specimen for the study of fretting fatigue.Next,the fretting fatigue contact problem is performed for 3 new cases of models that consist of a homogeneous and a heterogeneous part(single hole cell)in the contact area.The aim is to analyze the normal and shear stresses of these models and compare them with the results of the corresponding heterogeneous models based on the Direct Numerical Simulation(DNS)method.Finally,by comparing the computational time and%deviations,we draw conclusions about the reliability and effectiveness of the proposed method.
基金National Natural Science Foundation of China(No.61261029)
文摘Watershed segmentation is sensitive to noises and irregular details within the image,which frequently leads to a serious over-segmentation Linear filtering before watershed segmentation can reduce over-segmentation to some extent,however,it often causes the position offset of object contours.For the purpose of reducing over-segmentation to preserve the location of object contours,the watershed segmentation based on the hierarchical multi-scale modification of morphological gradient is proposed.Firstly,multi-scale morphological filtering was employed to smooth the original image.Then,the gradient image was divided into multi-levels by the volume of three-dimension topographic relief,where the lower gradient layers were further modifiedby morphological closing with larger-sized structuring-elements,and the higher layers with the smaller one.In this way,most local minimums caused by irregular details and noises can be removed,while region contour positions corresponding to the target area were largely preserved.Finally,morphological watershed algorithm was employed to implement segmentation on the modified gradient image.The experimental results show that the proposed method can greatly reduce the over-segmentation of the watershed and avoid the position offset of the object contours.
基金supported by National Natural Science Foundation of China (No. 61763037)Inner Mongolia Autonomous Region Natural Science Foundation of China(No. 2019LH06007)Science and Technology Plan Project of Inner Mongolia (No. 2019,2020GG028)。
文摘Working conditions of rolling bearings of wind turbine generators are complicated, and their vibration signals often show non-linear and non-stationary characteristics. In order to improve the efficiency of feature extraction of wind turbine rolling bearings and to strengthen the feature information, a new structural element and an adaptive algorithm based on the peak energy are proposed,which are combined with spectral correlation analysis to form a fault diagnosis algorithm for wind turbine rolling bearings. The proposed method firstly addresses the problem of impulsive signal omissions that are prone to occur in the process of fault feature extraction of traditional structural elements and proposes a "W" structural element to capture more characteristic information. Then, the proposed method selects the scale of multi-scale mathematical morphology, aiming at the problem of multi-scale mathematical morphology scale selection and structural element expansion law. An adaptive algorithm based on peak energy is proposed to carry out morphological scale selection and structural element expansion by improving the computing efficiency and enhancing the feature extraction effect.Finally, the proposed method performs spectral correlation analysis in the frequency domain for an unknown signal of the extracted feature and identifies the fault based on the correlation coefficient. The method is verified by numerical examples using experimental rig bearing data and actual wind field acquisition data and compared with traditional triangular and flat structural elements. The experimental results show that the new structural elements can more effectively extract the pulses in the signal and reduce noise interference,and the fault-diagnosis algorithm can accurately identify the fault category and improve the reliability of the results.
基金supported by the National Natural Science Foundation of China(Grants Nos 11302211,11390361,and 11572299).
文摘Multi-scale finite element method is adopted to simulate wood compression behavior under axial and transverse loading. Representative volume elements (RVE) of wood microfibril and cell are proposed to analyze orthotropic mechanical behavior. Lignin, hemicellulose and crystalline-amorphous cellulose core of spruce are concerned in spruce nanoscale model. The equivalent elastic modulus and yield strength of the microfibril are gained by the RVE simulation. The anisotropism of the crystalline-amorphous cellulose core brings the microfibril buckling deformation during compression loading. The failure mechanism of the cell-wall under axial compression is related to the distribution of amorphous cellulose and crystalline cellulose. According to the spruce cell observation by scanning electron microscope, numerical model of spruce cell is established using simplified circular hole and regular hexagon arrangement respectively. Axial and transverse compression loadings are taken into account in the numerical simulations. It indicates that the compression stress-strain curves of the numerical simulation are consistent with the experimental results. The wood microstructure arrangement has an important effect on the stress plateau during compression process. Cell-wall buckling in axial compression induces the stress value drops rapidly. The wide stress plateau duration means wood is with large energy dissipation under a low stress level. The numerical results show that loading velocity affects greatly wood microstructure failure modes in axial loading. For low velocity axial compression, shear sliding is the main failure mode. For high velocity axial compression, wood occur fold and collapse. In transverse compression, wood deformation is gradual and uniform, which brings stable stress plateau.
基金supported by the National Basic Research Program of China (Grant No. 2009CB623202)the Specialized Research Fund for the Doctoral Program of Higher Education,China (Grant No.20100092110049)the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘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.
基金financially supported by the National Natural Science Foundation of China(No.51571031)。
文摘A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by combining the proposed multi-scale crystal plasticity finite element method with experimental measurements.The complex local micro region load was progressively extracted from the simulation results of a macro model and applied to an established three-dimensional multi-grain microstructure model.Subsequently,the evolution histories of the grain shape,size,and orientation inside the adiabatic shear band were quantitatively simulated.The results corresponded closely to the experimental results obtained via transmission electron microscopy and precession electron diffraction.Furthermore,by calculating the grain rotation and temperature rise inside the adiabatic shear band,the microstructural softening and thermal softening effects of typical heavily-deformed α grains were successfully decoupled.The results revealed that the microstructural softening stress was triggered and then stabilized(in general)at a relatively high value.This indicated that the mechanical strength was lowered mainly by the grain orientation evolution or dynamic recrystallization occurring during early plastic deformation.Subsequently,thermal softening increased linearly and became the main softening mechanism.Noticeably,in the final stage,the thermal softening stress accounted for 78.4% of the total softening stress due to the sharp temperature increase,which inevitably leads to the stress collapse and potential failure of the alloy.
文摘Ceramic Matrix Composite (CMC) turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation, Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration (CVI) pro- cessed SiC/SiC-CMC are predicted by a Representative Volume Element (RVE) model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.
基金Financial support for‘‘Analysis of Mine Seismicity and Geotechnical Modeling for Improved Safety in Underground Coal Mines”by NIOSH/CDC
文摘The challenge of estimating mine-wide subsidence and linkages to seismicity over tabular deposits is addressed by a special finite element technique(dual node–dual mesh). Subsidence and mine-induced seismicity begins near the face when caving occurs and propagates to the surface as extraction reaches a critical extent. Thus, the challenge is to obtain details at the face at the meter scale and also at the surface over the whole mine at the kilometer scale. Interactions between old and new sections of a mine are automatically taken into account with this technique. The finite element method is well established technology based on fundamentals of physical laws, kinematics and material laws. With this technique, no empirical ‘‘scaling" or fitting computer output by input data ‘‘adjustment" to mine measurements is necessary. Capability is demonstrated for doing practical whole-mine subsidence analysis from first principles. Mine-induced seismicity is shown to correlate well with face advance and element failure.
文摘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.
基金This project is supported by Motorola (China) Electronics Ltd. and the work of Pingwen Zhang is also partially supported by Special Funds for Major State Basic Research Projects of China G1999032804.
文摘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.
基金This work is supported by National KeyR&D Program of China(Gran Nt oS.2016YFB0300400 and 2016YFB0300404).
文摘The void evolution of large-section plastic mold steel during multi-directional forging(MDF)was investigated using multiscale analysis.To simulate the forging process of the plastic mold steel(SDP1 steel)and realize micro-void reconstruction in a representative volume element(RVE),MDF experiment and void-characteristic evaluation of the SDP1 steel were carried out.Traditional upsetting and stretching forging(TUSF)and MDF were simulated to comparatively analyze the evolution of temperature,effective stress,and effective strain.By embedding RVE with a micro-void and using boundary condition by point tracking into the forging process,the single-void evolution in TUSF and MDF was studied.The effect of void orientation on single-void evolution was also investigated.The multi-scale analysis revealed the following results.(1)Compared with TUSF,MDF achieved a higher efficiency in void closure.(2)The closing efficiency of the void increased with the increase in angle h(the angle between the Z and long axes of the void).(3)The closing efficiency increased with the increase in the orientation angle during the forging process.On the basis of the important role of the main stress in each forging step on the void closure,an integral formula of the main stress was proposed.When main compressive-stress integration reached-0.4,the closed state of the void could be accurately determined.
基金BASTRI Subtopic Research about Digital Sampler Technology of Body Structure Performance Study Based on Big Data Calculation Model,China(No.MIIT Civil aircraft special purpose MJ-2017-F-20)
文摘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.
基金The Project Supported by National Natural Science Foundation of China (No.19801006)and SpecialFunds for Major State Basic Rese
文摘Deals with a study which discussed the mechanical behaviour for subdivided periodic elastic structures of composite materials, from the viewpoint of macro- and meso-scale coupling. Multiscale asymptotic expansion and truncation error estimates; Discussion on multiscale finite element method; Details of higher order difference quotients and total error estimates; Numerical experiments.
基金financially supported by the Key Research Program of Frontier Sciences,CAS(No.QYZDJ-SSWJSC015)the National Natural Science Foundation of China(Nos.51931009,51871214 and 51871215)the National Key R&D Program of China(No.2017YFB0703104)。
文摘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.