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 mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
Geared-rotor systems are critical components in mechanical applications,and their performance can be severely affected by faults,such as profile errors,wear,pitting,spalling,flaking,and cracks.Profile errors in gear t...Geared-rotor systems are critical components in mechanical applications,and their performance can be severely affected by faults,such as profile errors,wear,pitting,spalling,flaking,and cracks.Profile errors in gear teeth are inevitable in manufacturing and subsequently accumulate during operations.This work aims to predict the status of gear profile deviations based on gear dynamics response using the digital model of an experimental rig setup.The digital model comprises detailed CAD models and has been validated against the expected physical behavior using commercial finite element analysis software.The different profile deviations are then modeled using gear charts,and the dynamic response is captured through simulations.The various features are then obtained by signal processing,and various ML models are then evaluated to predict the fault/no-fault condition for the gear.The best performance is achieved by an artificial neural network with a prediction accuracy of 97.5%,which concludes a strong influence on the dynamics of the gear rotor system due to profile deviations.展开更多
In this study,an alternative modelling approach for absorbed hydrogen stress corrosion cracking(SCC)is proposed,with hydrogen-enhanced decohesion(HEDE)identified as the key failure mechanism.All analyses have been per...In this study,an alternative modelling approach for absorbed hydrogen stress corrosion cracking(SCC)is proposed,with hydrogen-enhanced decohesion(HEDE)identified as the key failure mechanism.All analyses have been performed by utilising only ABAQUS standard elements,COH2D4T and CPE4T,already available within the software and without the need to develop external subroutines.The study also tends to highlight the criticality of implementing a correct Traction Separation Law(TSL)curve to simulate the hydrogen diffusion within the specimen and using the concept of dynamic hydrogen penetration by continuously updating the hydrogen concentration boundary conditions as the crack propagates.In conclusion,this study successfully demonstrated that standard software elements(COH2D4T and CPE4T)can effectively model physical problems and crack velocity propagation without custom subroutines.It emphasized that while the specific shape of the Traction-Separation Law(TSL)is less critical,its correct implementation is vital for simulating dynamic hydrogen coverage.Crucially,excluding this dynamic coverage—a common practice—risks significantly underestimating crack propagation speed.Although results incorporating dynamic coverage aligned well with experimental data,minor discrepancies are likely due to unmodeled factors like material property variations,hydrogen trapping,temperature,and granular microstructure,which are proposed for future research.展开更多
In this work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on...In this work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on the actual 3D microstructure achieved by synchrotron tomography.The results show that the average grain size of composite increases from 0.57μm of 8μm-SiC/ZA63 to 8.73μm of 50μm-SiC/ZA63.The type of texture transforms from the typicalfiber texture in 8μm-SiC/ZA63 to intense basal texture in 50μm-SiC/ZA63 composite and the intensity of texture increases sharply with increase of SiC particle size.The dynamic recrystallization(DRX)mechanism is also changed with increasing SiC particle size.Experimental and simulation results verify that the strength and elongation both decrease with increase of SiC particle size.The 8μm-SiC/ZA63 composite possesses the optimal mechanical property with yield strength(YS)of 383 MPa,ultimate tensile strength(UTS)of 424 MPa and elongation of 6.3%.The outstanding mechanical property is attributed to the ultrafine grain size,high-density precipitates and dislocation,good loading transfer effect and the interface bonding between SiC and matrix,as well as the weakened basal texture.The simulation results reveal that the micro-cracks tend to initiate at the interface between SiC and matrix,and then propagate along the interface between particle and Mg matrix or at the high strain and stress regions,and further connect with other micro-cracks.The main fracture mechanism in 8μm-SiC/ZA63 composite is ductile damage of matrix and interfacial debonding.With the increase of particle size,interface strength and particle strength decrease,and interface debonding and particle rupture become the main fracture mechanism in the 30μm-and 50μm-SiC/ZA63 composites.展开更多
In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a ge...In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.展开更多
Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient applicat...Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.展开更多
Convective heat transfer associated with the circulation of porefluid in porous rocks and fractures within the upper crust of the Earth is substantial when the temperature gradient is sufficiently high. In order to un...Convective heat transfer associated with the circulation of porefluid in porous rocks and fractures within the upper crust of the Earth is substantial when the temperature gradient is sufficiently high. In order to understand the process of Snpolymetallic mineralization in the Dachang ore district of Guangxi, a finite element method has been used in this study to simulate both pore-fluid flow and heat transfer in this district. On the basis of related geological, tectonic and geophysical constraints, a computational model was established. It enables a computational simulation and sensitivity analysis to be carried out for investigating ore-forming pore-fluid flow and other key factors that may affect hydrothermal ore genesis in the district. The related simulation results have indicated that： （1） permeable fault zones in the Dacbang ore district can serve as preferential pathways for pore-fluid flow on a regional-scale; and （2） the pore-fluid flow can affect the salinity distribution. This latter factor is part of the reason why Sn-polymetallic mineralization has taken place in this district.展开更多
Recent seismic events have raised concerns over the safety and vulnerability of reinforced concrete moment resisting frame ＂RC-MRF＂ buildings. The seismic response of such buildings is greatly dependent on the compu...Recent seismic events have raised concerns over the safety and vulnerability of reinforced concrete moment resisting frame ＂RC-MRF＂ buildings. The seismic response of such buildings is greatly dependent on the computational tools used and the inherent assumptions in the modelling process. Thus, it is essential to investigate the sensitivity of the response demands to the corresponding modelling assumption. Many parameters and assumptions are justified to generate effective structural finite element（FE） models of buildings to simulate lateral behaviour and evaluate seismic design demands. As such, the present study focuses on the development of reliable FE models with various levels of refinement. The effects of the FE modelling assumptions on the seismic response demands on the design of buildings are investigated. the predictive ability of a FE model is tied to the accuracy of numerical analysis; a numerical analysis is performed for a series of symmetric buildings in active seismic zones. The results of the seismic response demands are presented in a comparative format to confirm drift and strength limits requirements. A proposed model is formulated based on a simplified modeling approach, where the most refined model is used to calibrate the simplified model.展开更多
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.展开更多
The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are comple...The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are complex and inter-connected.Finite element method(FEM)is considered as an effective method to predict process variables and reveal microscopic physical phenomena in the cutting process.Therefore,the finite element(FE)simulation is used to research the conventional and micro scale cutting process,and the differences in the establishment of process variable FE simulation models are distinguished,thereby improving the accuracy of FE simulation.The reliability and effectiveness of FE simulation model largely depend on the accuracy of the simulation method,constitutive model,friction model,damage model in describing mesh element,the dynamic mechanical behavior of materials,the tool-chip-workpiece contact process and the chip formation mechanism.In this paper,the FE models of conventional and micro process variables are comprehensively and up-to-date reviewed for different materials and machining methods.The purpose is to establish a FE model that is more in line with the real cutting conditions,and to provide the possibility for optimizing the cutting process variables.The development direction of FE simulation of metal cutting process is discussed,which provides guidance for future cutting process modeling.展开更多
The microstructure-based finite element modeling was conducted to study the mechanical properties of Q&P 980 steel at the microscopic level.The two-dimensional representative volume elements of real microstructure...The microstructure-based finite element modeling was conducted to study the mechanical properties of Q&P 980 steel at the microscopic level.The two-dimensional representative volume elements of real microstructure were obtained from electron backscattered diffraction mapping.Mecking-Kocks equation was used to predict the constitutive strain-stress relationships of individual phases.Mechanical-induced martensitic transformation takes place when the driving force exceeds the critical driving force according to a stress-invariant-based model.The macroscopic stress-strain curves and the work-hardening rate curves obtained from modeling fit well with the experimental results.The simulation results also indicate that the local distributions of stress and strain in constituent phases are dependent on their strength.Soft ferrite carries the highest strain,while hard mechanical-induced martensite carries the highest stress.By comparing the modeling results of the microstructures with and without austenite,it shows that the transformation of retained austenite to hard martensite can increase the .work-hardening ability and hence improve the strength and ductility of the steel.The detailed finite element modeling methods and results are presented and discussed.展开更多
A linear viscoelastic finite element model was built to investigate factors that influenced the intraocular pressure (IOP) elevations due to micro-volumetric changes in the eye at three different rates. The viscoelast...A linear viscoelastic finite element model was built to investigate factors that influenced the intraocular pressure (IOP) elevations due to micro-volumetric changes in the eye at three different rates. The viscoelastic properties of the cornea and the sclera, including the instantaneous modulus, equilibrium modulus, and relaxation time constants, parametrically varied to examine their effects on IOP elevations at different rates of volumetric changes. The simulated responses were in good agreement with the previously reported experimental results obtained from porcine globes, showing the general trend of higher IOP elevations at faster rates. The simulations showed that all viscoelastic properties influenced the profile of the dynamic IOP due to volumetric changes, and the relative significance of a specific parameter was highly dependent on the rate of change.展开更多
In finite element modeling of impact,it is necessary to define appropriate values of the normal contact stiffness,Kn,and the Integration Time Step(ITS).Because impacts are usually of very short duration,very small ITS...In finite element modeling of impact,it is necessary to define appropriate values of the normal contact stiffness,Kn,and the Integration Time Step(ITS).Because impacts are usually of very short duration,very small ITSs are required.Moreover,the selection of a suitable value of Kn is a critical issue,as the impact behavior depends dramatically on this parameter.In this work,a number of experimental tests and finite element analyses have been performed in order to obtain an appropriate value of Kn for the interaction between a bristle of a gutter brush for road sweeping and a concrete surface.Furthermore,a suitable ITS is determined.The experiments consist of releasing a steel bristle that is placed vertically at a certain distance from a concrete surface and tracking the impact.Similarly,in the finite element analyses,a beam is modeled in free fall and impacting a surface;contact and target elements are attached to the beam and the surface,respectively.The results of the experiments and the modeling are integrated through the principle of conservation of energy,the principle of linear impulse and momentum,and Newton’s second law.The results demonstrate that,for the case studied,Kn and the impact time tend to be independent of the velocity just before impact and that Kn has a very large variation,as concrete is a composite material with a rough surface.Also,the ratio between the largest height of the bristle after impact and the initial height tends to be constant.展开更多
The consolidation process of SiC<sub>f</sub>/Ti-6Al-4V composites by matrix-coated fiber (MCF) method via hot pressing was investigated using finite element modeling (FEM). By analyzing the elastic–plasti...The consolidation process of SiC<sub>f</sub>/Ti-6Al-4V composites by matrix-coated fiber (MCF) method via hot pressing was investigated using finite element modeling (FEM). By analyzing the elastic–plastic contact deformation of the representative aligned coated fibers, the consolidation maps delineating the time–temperature–pressure relationship for full densification were constructed. Both the flow coefficient and the contact area coefficient used to describe the contact deformation were calculated according to the model. In addition, the effect of fiber content on matrix stress distribution was analyzed. The results show that fiber content is a significant factor that influences the densification process. Higher fiber content will lower the consolidation rate.展开更多
Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By rep...Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By repeated tests and revisions, the boundary conditions of the model are determined. And then the background stress field, the stress field caused by fault creep and the stress field triggered by strong earthquake in Sichuan-Yunnan region, as well as their dynamic variations are calculated. The results indicate that the latter earthquake often occurs in the area with positive Coulomb rupture stress change associated with the former one, the former earthquake has a triggering effect on the latter one to a certain extent, and strong earthquake often occur in groups under the background of high stress, which is of great significance for distinguishing seismic anomalies, as well as for improving the level of earthquake prediction.展开更多
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 this study,we propose a novel and simple exact semi-analytical model for superelastic Shape Memory Alloy(SMA)wire reinforced composites subjected to bending loads.In order to study the mechanical response of the co...In this study,we propose a novel and simple exact semi-analytical model for superelastic Shape Memory Alloy(SMA)wire reinforced composites subjected to bending loads.In order to study the mechanical response of the composite during loading/unloading,a Representative Volume Element(RVE)is extracted to examine the bending response of the composite.Analytical moment–curvature,and shear force-shear strain relations are derived based on a 3-Dimensional(3 D)thermomechanical SMA model and Timoshenko beam theory.The composite Simpson’s rule is adopted to numerically solve the exact analytical moment–curvature and shear force-shear strain relationships.Results including the moment–curvature response,axial stress distribution along the vertical and longitudinal directions,martensite volume fraction,and the tip deflection are reported and validated against 3 D finite element simulations.The influence of temperature,martensite volume fraction distribution,and matrix stiffness on the mechanical performance of the composite is also investigated.In particular,the composite is found to behave superelastically under certain conditions of temperature,SMA volume fraction,and elastic stiffness of the matrix.Such behavior is advantageous in applications requiring large recoverable strains or high energy dissipation density.展开更多
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.展开更多
基金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.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
文摘Geared-rotor systems are critical components in mechanical applications,and their performance can be severely affected by faults,such as profile errors,wear,pitting,spalling,flaking,and cracks.Profile errors in gear teeth are inevitable in manufacturing and subsequently accumulate during operations.This work aims to predict the status of gear profile deviations based on gear dynamics response using the digital model of an experimental rig setup.The digital model comprises detailed CAD models and has been validated against the expected physical behavior using commercial finite element analysis software.The different profile deviations are then modeled using gear charts,and the dynamic response is captured through simulations.The various features are then obtained by signal processing,and various ML models are then evaluated to predict the fault/no-fault condition for the gear.The best performance is achieved by an artificial neural network with a prediction accuracy of 97.5%,which concludes a strong influence on the dynamics of the gear rotor system due to profile deviations.
文摘In this study,an alternative modelling approach for absorbed hydrogen stress corrosion cracking(SCC)is proposed,with hydrogen-enhanced decohesion(HEDE)identified as the key failure mechanism.All analyses have been performed by utilising only ABAQUS standard elements,COH2D4T and CPE4T,already available within the software and without the need to develop external subroutines.The study also tends to highlight the criticality of implementing a correct Traction Separation Law(TSL)curve to simulate the hydrogen diffusion within the specimen and using the concept of dynamic hydrogen penetration by continuously updating the hydrogen concentration boundary conditions as the crack propagates.In conclusion,this study successfully demonstrated that standard software elements(COH2D4T and CPE4T)can effectively model physical problems and crack velocity propagation without custom subroutines.It emphasized that while the specific shape of the Traction-Separation Law(TSL)is less critical,its correct implementation is vital for simulating dynamic hydrogen coverage.Crucially,excluding this dynamic coverage—a common practice—risks significantly underestimating crack propagation speed.Although results incorporating dynamic coverage aligned well with experimental data,minor discrepancies are likely due to unmodeled factors like material property variations,hydrogen trapping,temperature,and granular microstructure,which are proposed for future research.
基金supported by the National Natural Science Foundation of China[51974058,52371005,52022017,51927801]the Fundamental Research Funds for the Central Universities(DUT23YG104).
文摘In this work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on the actual 3D microstructure achieved by synchrotron tomography.The results show that the average grain size of composite increases from 0.57μm of 8μm-SiC/ZA63 to 8.73μm of 50μm-SiC/ZA63.The type of texture transforms from the typicalfiber texture in 8μm-SiC/ZA63 to intense basal texture in 50μm-SiC/ZA63 composite and the intensity of texture increases sharply with increase of SiC particle size.The dynamic recrystallization(DRX)mechanism is also changed with increasing SiC particle size.Experimental and simulation results verify that the strength and elongation both decrease with increase of SiC particle size.The 8μm-SiC/ZA63 composite possesses the optimal mechanical property with yield strength(YS)of 383 MPa,ultimate tensile strength(UTS)of 424 MPa and elongation of 6.3%.The outstanding mechanical property is attributed to the ultrafine grain size,high-density precipitates and dislocation,good loading transfer effect and the interface bonding between SiC and matrix,as well as the weakened basal texture.The simulation results reveal that the micro-cracks tend to initiate at the interface between SiC and matrix,and then propagate along the interface between particle and Mg matrix or at the high strain and stress regions,and further connect with other micro-cracks.The main fracture mechanism in 8μm-SiC/ZA63 composite is ductile damage of matrix and interfacial debonding.With the increase of particle size,interface strength and particle strength decrease,and interface debonding and particle rupture become the main fracture mechanism in the 30μm-and 50μm-SiC/ZA63 composites.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12371393,11971150 and 11801143)Natural Science Foundation of Henan Province(Grant No.242300421047).
文摘In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.
文摘Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.
基金financially supported by the Natural Science Foundation of China(Grant No:10872219)
文摘Convective heat transfer associated with the circulation of porefluid in porous rocks and fractures within the upper crust of the Earth is substantial when the temperature gradient is sufficiently high. In order to understand the process of Snpolymetallic mineralization in the Dachang ore district of Guangxi, a finite element method has been used in this study to simulate both pore-fluid flow and heat transfer in this district. On the basis of related geological, tectonic and geophysical constraints, a computational model was established. It enables a computational simulation and sensitivity analysis to be carried out for investigating ore-forming pore-fluid flow and other key factors that may affect hydrothermal ore genesis in the district. The related simulation results have indicated that： （1） permeable fault zones in the Dacbang ore district can serve as preferential pathways for pore-fluid flow on a regional-scale; and （2） the pore-fluid flow can affect the salinity distribution. This latter factor is part of the reason why Sn-polymetallic mineralization has taken place in this district.
基金Scientific Research Deanship,Taibah University Grant No.6363/436
文摘Recent seismic events have raised concerns over the safety and vulnerability of reinforced concrete moment resisting frame ＂RC-MRF＂ buildings. The seismic response of such buildings is greatly dependent on the computational tools used and the inherent assumptions in the modelling process. Thus, it is essential to investigate the sensitivity of the response demands to the corresponding modelling assumption. Many parameters and assumptions are justified to generate effective structural finite element（FE） models of buildings to simulate lateral behaviour and evaluate seismic design demands. As such, the present study focuses on the development of reliable FE models with various levels of refinement. The effects of the FE modelling assumptions on the seismic response demands on the design of buildings are investigated. the predictive ability of a FE model is tied to the accuracy of numerical analysis; a numerical analysis is performed for a series of symmetric buildings in active seismic zones. The results of the seismic response demands are presented in a comparative format to confirm drift and strength limits requirements. A proposed model is formulated based on a simplified modeling approach, where the most refined model is used to calibrate the simplified model.
基金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 Natural Science Foundation of China(No.52175393)。
文摘The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are complex and inter-connected.Finite element method(FEM)is considered as an effective method to predict process variables and reveal microscopic physical phenomena in the cutting process.Therefore,the finite element(FE)simulation is used to research the conventional and micro scale cutting process,and the differences in the establishment of process variable FE simulation models are distinguished,thereby improving the accuracy of FE simulation.The reliability and effectiveness of FE simulation model largely depend on the accuracy of the simulation method,constitutive model,friction model,damage model in describing mesh element,the dynamic mechanical behavior of materials,the tool-chip-workpiece contact process and the chip formation mechanism.In this paper,the FE models of conventional and micro process variables are comprehensively and up-to-date reviewed for different materials and machining methods.The purpose is to establish a FE model that is more in line with the real cutting conditions,and to provide the possibility for optimizing the cutting process variables.The development direction of FE simulation of metal cutting process is discussed,which provides guidance for future cutting process modeling.
基金National Key Research and Development Program of China (No.2017YFB0304401)National Natural Science Foundation of China (Nos.U1564203,51571141 and 51201105).
文摘The microstructure-based finite element modeling was conducted to study the mechanical properties of Q&P 980 steel at the microscopic level.The two-dimensional representative volume elements of real microstructure were obtained from electron backscattered diffraction mapping.Mecking-Kocks equation was used to predict the constitutive strain-stress relationships of individual phases.Mechanical-induced martensitic transformation takes place when the driving force exceeds the critical driving force according to a stress-invariant-based model.The macroscopic stress-strain curves and the work-hardening rate curves obtained from modeling fit well with the experimental results.The simulation results also indicate that the local distributions of stress and strain in constituent phases are dependent on their strength.Soft ferrite carries the highest strain,while hard mechanical-induced martensite carries the highest stress.By comparing the modeling results of the microstructures with and without austenite,it shows that the transformation of retained austenite to hard martensite can increase the .work-hardening ability and hence improve the strength and ductility of the steel.The detailed finite element modeling methods and results are presented and discussed.
文摘A linear viscoelastic finite element model was built to investigate factors that influenced the intraocular pressure (IOP) elevations due to micro-volumetric changes in the eye at three different rates. The viscoelastic properties of the cornea and the sclera, including the instantaneous modulus, equilibrium modulus, and relaxation time constants, parametrically varied to examine their effects on IOP elevations at different rates of volumetric changes. The simulated responses were in good agreement with the previously reported experimental results obtained from porcine globes, showing the general trend of higher IOP elevations at faster rates. The simulations showed that all viscoelastic properties influenced the profile of the dynamic IOP due to volumetric changes, and the relative significance of a specific parameter was highly dependent on the rate of change.
文摘In finite element modeling of impact,it is necessary to define appropriate values of the normal contact stiffness,Kn,and the Integration Time Step(ITS).Because impacts are usually of very short duration,very small ITSs are required.Moreover,the selection of a suitable value of Kn is a critical issue,as the impact behavior depends dramatically on this parameter.In this work,a number of experimental tests and finite element analyses have been performed in order to obtain an appropriate value of Kn for the interaction between a bristle of a gutter brush for road sweeping and a concrete surface.Furthermore,a suitable ITS is determined.The experiments consist of releasing a steel bristle that is placed vertically at a certain distance from a concrete surface and tracking the impact.Similarly,in the finite element analyses,a beam is modeled in free fall and impacting a surface;contact and target elements are attached to the beam and the surface,respectively.The results of the experiments and the modeling are integrated through the principle of conservation of energy,the principle of linear impulse and momentum,and Newton’s second law.The results demonstrate that,for the case studied,Kn and the impact time tend to be independent of the velocity just before impact and that Kn has a very large variation,as concrete is a composite material with a rough surface.Also,the ratio between the largest height of the bristle after impact and the initial height tends to be constant.
基金financially supported by the National Natural Science Foundation of China(Nos.51071122 and51271147)
文摘The consolidation process of SiC<sub>f</sub>/Ti-6Al-4V composites by matrix-coated fiber (MCF) method via hot pressing was investigated using finite element modeling (FEM). By analyzing the elastic–plastic contact deformation of the representative aligned coated fibers, the consolidation maps delineating the time–temperature–pressure relationship for full densification were constructed. Both the flow coefficient and the contact area coefficient used to describe the contact deformation were calculated according to the model. In addition, the effect of fiber content on matrix stress distribution was analyzed. The results show that fiber content is a significant factor that influences the densification process. Higher fiber content will lower the consolidation rate.
文摘Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By repeated tests and revisions, the boundary conditions of the model are determined. And then the background stress field, the stress field caused by fault creep and the stress field triggered by strong earthquake in Sichuan-Yunnan region, as well as their dynamic variations are calculated. The results indicate that the latter earthquake often occurs in the area with positive Coulomb rupture stress change associated with the former one, the former earthquake has a triggering effect on the latter one to a certain extent, and strong earthquake often occur in groups under the background of high stress, which is of great significance for distinguishing seismic anomalies, as well as for improving the level of earthquake prediction.
基金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 financial support of Khalifa University through research grant No.CIRA 2019024。
文摘In this study,we propose a novel and simple exact semi-analytical model for superelastic Shape Memory Alloy(SMA)wire reinforced composites subjected to bending loads.In order to study the mechanical response of the composite during loading/unloading,a Representative Volume Element(RVE)is extracted to examine the bending response of the composite.Analytical moment–curvature,and shear force-shear strain relations are derived based on a 3-Dimensional(3 D)thermomechanical SMA model and Timoshenko beam theory.The composite Simpson’s rule is adopted to numerically solve the exact analytical moment–curvature and shear force-shear strain relationships.Results including the moment–curvature response,axial stress distribution along the vertical and longitudinal directions,martensite volume fraction,and the tip deflection are reported and validated against 3 D finite element simulations.The influence of temperature,martensite volume fraction distribution,and matrix stiffness on the mechanical performance of the composite is also investigated.In particular,the composite is found to behave superelastically under certain conditions of temperature,SMA volume fraction,and elastic stiffness of the matrix.Such behavior is advantageous in applications requiring large recoverable strains or high energy dissipation density.
基金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.
