An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for ...An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for the equivalent continuum. A yielding model is formulated and compared with the results of other models. The bedding-in effect is considered to include the compliance of the lattice joints. The predicted stiffness and strength are in good agreement with the experimental data, validating the present model in the prediction of the mechanical properties of stretching dominated lattice structures.展开更多
In this paper,a continuum model with dynamic earth pressure coefficient is established to describe the granular slump process by introducingμ(I)rheology.This rheology is adopted to quantify the normal stresses in our...In this paper,a continuum model with dynamic earth pressure coefficient is established to describe the granular slump process by introducingμ(I)rheology.This rheology is adopted to quantify the normal stresses in our proposed model rather than shear stresses in classical models.The constitutive laws of different depth-averaged continuum approaches including the hydrodynamic,Savage–Hutter and proposed models are comparatively investigated in terms of the rheological effects on the spread of a granular column.The simulation results indicate that the proposed dynamic model captures some significant features during granular slump on inclined planes with different inclination angles(for example,the runout distance,runout time,and final profile).The proposed model can also reproduce the inner static sided axisymmetric region observed in tests when the granular column's initial aspect ratio(ratio of height to radii)is small.展开更多
Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surfac...Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.展开更多
The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of ...The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these defo...Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these deformation patterns.This study utilized geophysical surveys and the continuum-based discrete element method(CDEM)to investigate how fault activity influences rock deformation and failure.The results demonstrate that:1)Acting in mechanically weak zones,faults exerted a pronounced barrier effect on deformation propagation and stress redistribution within the surrounding rock,leading to markedly divergent displacement patterns on either side of the fault plane.Comparative analyses between single-fault and double-fault models revealed an 18%−22%expansion of the damage zone under the latter,together with significantly intensified deformation and failure;2)The double-fault model exhibited a larger maximum cumulative vertical displacement and a spatial shift in the location of peak deformation,thereby posing a heightened threat to mine safety;3)Acting in an orebody substitute,backfill effectively constrained surrounding rock deformation,enhanced its load-bearing capacity,and delayed the overburden subsidence.Nevertheless,backfill only reduced the amplitude of deformation;it could not entirely prevent settlement.These findings provide essential theoretical insights and foundational knowledge for safer submarine mining practices.展开更多
The search for fast radio bursts(FRBs) is a hot topic in current radio astronomy studies. In this work, we carry out a single pulse search with a very long baseline interferometry(VLBI) pulsar observation data set usi...The search for fast radio bursts(FRBs) is a hot topic in current radio astronomy studies. In this work, we carry out a single pulse search with a very long baseline interferometry(VLBI) pulsar observation data set using both auto spectrum and cross spectrum search methods. The cross spectrum method,first proposed in Liu et al., maximizes the signal power by fully utilizing the fringe phase information of the baseline cross spectrum. The auto spectrum search method is based on the popular pulsar software package PRESTO, which extracts single pulses from the auto spectrum of each station. According to our comparison, the cross spectrum method is able to enhance the signal power and therefore extract single pulses from data contaminated by high levels of radio frequency interference(RFI), which makes it possible to carry out a search for FRBs in regular VLBI observations when RFI is present.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
Fatigue life and reliability of aero-engine blade are always of important significance to flight safety.The establishment of damage model is one of the key factors in blade fatigue research.Conventional linear Miner＇...Fatigue life and reliability of aero-engine blade are always of important significance to flight safety.The establishment of damage model is one of the key factors in blade fatigue research.Conventional linear Miner＇s sum method is not suitable for aero-engine because of its low accuracy.A back propagation neutral network（BPNN） based on the combination of Levenberg-Marquardt（LM） and finite element method（FEM） is used to describe process of nonlinear damage accumulation behavior in material and predict fatigue life of the blade.Fatigue tests of standard specimen made from TC4 are carried out to obtain material fatigue parameters and S-N curve.A nonlinear continuum damage model（CDM）,based on the BPNN with one hidden layer and ten neurons,is built to investigate the nonlinear damage accumulation behavior,in which the results from the tests are used as training set.Comparing with linear models and previous nonlinear models,BPNN has the lowest calculation error in full load range.It has significant accuracy when the load is below 500 MPa.Especially,when the load is 350 MPa,the calculation error of the BPNN is only 0.4%.The accurate model of the blade is built by using 3D coordinate measurement technology.The loading cycle in fatigue analysis is defined from takeoff to cruise in 10 min,and the load history is obtained from finite element analysis（FEA）.Then the fatigue life of the compressor blade is predicted by using the BPNN model.The final fatigue life of the aero-engine blade is 6.55 104 cycles（10 916 h） based on the BPNN model,which is effective for the virtual design of aero-engine blade.展开更多
This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock m...This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.展开更多
In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. ...In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.展开更多
The fine-scale heterogeneity of granular material is characterized by its polydisperse microstructure with randomness and no periodicity. To predict the mechanical response of the material as the microstructure evolve...The fine-scale heterogeneity of granular material is characterized by its polydisperse microstructure with randomness and no periodicity. To predict the mechanical response of the material as the microstructure evolves, it is demonstrated to develop computational multiscale methods using discrete particle assembly-Cosserat continuum modeling in micro- and macro- scales,respectively. The computational homogenization method and the bridge scale method along the concurrent scale linking approach are briefly introduced. Based on the weak form of the Hu-Washizu variational principle, the mixed finite element procedure of gradient Cosserat continuum in the frame of the second-order homogenization scheme is developed. The meso-mechanically informed anisotropic damage of effective Cosserat continuum is characterized and identified and the microscopic mechanisms of macroscopic damage phenomenon are revealed. c 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi： 10.1063/2.1301101]展开更多
基金Project supported by the Key Project of Chinese Ministry of Education (No.106015).
文摘An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for the equivalent continuum. A yielding model is formulated and compared with the results of other models. The bedding-in effect is considered to include the compliance of the lattice joints. The predicted stiffness and strength are in good agreement with the experimental data, validating the present model in the prediction of the mechanical properties of stretching dominated lattice structures.
基金funded by the National Natural Science Foundation of China(Grant Nos.52008261 and 52090081)the Key Research and Development Program of China Railway(Grant No.K2019G033).
文摘In this paper,a continuum model with dynamic earth pressure coefficient is established to describe the granular slump process by introducingμ(I)rheology.This rheology is adopted to quantify the normal stresses in our proposed model rather than shear stresses in classical models.The constitutive laws of different depth-averaged continuum approaches including the hydrodynamic,Savage–Hutter and proposed models are comparatively investigated in terms of the rheological effects on the spread of a granular column.The simulation results indicate that the proposed dynamic model captures some significant features during granular slump on inclined planes with different inclination angles(for example,the runout distance,runout time,and final profile).The proposed model can also reproduce the inner static sided axisymmetric region observed in tests when the granular column's initial aspect ratio(ratio of height to radii)is small.
文摘Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.
基金Sponsored by the Ministerial Level Advanced Research Foundation (010896367)
文摘The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金Project(42072305)supported by the National Natural Science Foundation of China。
文摘Seabed mining operations have been found to induce significant movement and deformation in overlying rock strata,posing serious threats to mining safety.The presence of geological faults further complicates these deformation patterns.This study utilized geophysical surveys and the continuum-based discrete element method(CDEM)to investigate how fault activity influences rock deformation and failure.The results demonstrate that:1)Acting in mechanically weak zones,faults exerted a pronounced barrier effect on deformation propagation and stress redistribution within the surrounding rock,leading to markedly divergent displacement patterns on either side of the fault plane.Comparative analyses between single-fault and double-fault models revealed an 18%−22%expansion of the damage zone under the latter,together with significantly intensified deformation and failure;2)The double-fault model exhibited a larger maximum cumulative vertical displacement and a spatial shift in the location of peak deformation,thereby posing a heightened threat to mine safety;3)Acting in an orebody substitute,backfill effectively constrained surrounding rock deformation,enhanced its load-bearing capacity,and delayed the overburden subsidence.Nevertheless,backfill only reduced the amplitude of deformation;it could not entirely prevent settlement.These findings provide essential theoretical insights and foundational knowledge for safer submarine mining practices.
基金funded by the National Natural Science Foundation of China (Grant Nos. 11373061, 11573057 and U1631122)the Key Laboratory of Radio Astronomy of Chinese Academy of Sciences (CAS)+2 种基金Shanghai Key Laboratory of Space Navigation and Positioning Techniques (ZZXT-201702)the CAS Key Technology Talent Programthe National R&D Infrastructure and Facility Development Program of China, “Fundamental Science Data Sharing Platform” (DKA2017-12-02-09)
文摘The search for fast radio bursts(FRBs) is a hot topic in current radio astronomy studies. In this work, we carry out a single pulse search with a very long baseline interferometry(VLBI) pulsar observation data set using both auto spectrum and cross spectrum search methods. The cross spectrum method,first proposed in Liu et al., maximizes the signal power by fully utilizing the fringe phase information of the baseline cross spectrum. The auto spectrum search method is based on the popular pulsar software package PRESTO, which extracts single pulses from the auto spectrum of each station. According to our comparison, the cross spectrum method is able to enhance the signal power and therefore extract single pulses from data contaminated by high levels of radio frequency interference(RFI), which makes it possible to carry out a search for FRBs in regular VLBI observations when RFI is present.
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金supported by National Natural Science Foundation of China (Grant No. 60879002)Tianjin Municipal Science and Technology Support Plan of China (Grant No. 10ZCKFGX03800)
文摘Fatigue life and reliability of aero-engine blade are always of important significance to flight safety.The establishment of damage model is one of the key factors in blade fatigue research.Conventional linear Miner＇s sum method is not suitable for aero-engine because of its low accuracy.A back propagation neutral network（BPNN） based on the combination of Levenberg-Marquardt（LM） and finite element method（FEM） is used to describe process of nonlinear damage accumulation behavior in material and predict fatigue life of the blade.Fatigue tests of standard specimen made from TC4 are carried out to obtain material fatigue parameters and S-N curve.A nonlinear continuum damage model（CDM）,based on the BPNN with one hidden layer and ten neurons,is built to investigate the nonlinear damage accumulation behavior,in which the results from the tests are used as training set.Comparing with linear models and previous nonlinear models,BPNN has the lowest calculation error in full load range.It has significant accuracy when the load is below 500 MPa.Especially,when the load is 350 MPa,the calculation error of the BPNN is only 0.4%.The accurate model of the blade is built by using 3D coordinate measurement technology.The loading cycle in fatigue analysis is defined from takeoff to cruise in 10 min,and the load history is obtained from finite element analysis（FEA）.Then the fatigue life of the compressor blade is predicted by using the BPNN model.The final fatigue life of the aero-engine blade is 6.55 104 cycles（10 916 h） based on the BPNN model,which is effective for the virtual design of aero-engine blade.
基金supported by the National Natural Science Youth Foundation of China(Grant No.51309101)the Henan Province Major Scientific and Technological Projects(Grant No.172102210372)the Cooperative Project of Production,Teaching and Research in Henan Province(Grant No.18210700031)
文摘This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.
文摘In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.
基金supported by the National Natural Science Foundation of China(11072046,10672033,90715011 and 11102036)the National Basic Research and Development Program(973Program,2010CB731502)
文摘The fine-scale heterogeneity of granular material is characterized by its polydisperse microstructure with randomness and no periodicity. To predict the mechanical response of the material as the microstructure evolves, it is demonstrated to develop computational multiscale methods using discrete particle assembly-Cosserat continuum modeling in micro- and macro- scales,respectively. The computational homogenization method and the bridge scale method along the concurrent scale linking approach are briefly introduced. Based on the weak form of the Hu-Washizu variational principle, the mixed finite element procedure of gradient Cosserat continuum in the frame of the second-order homogenization scheme is developed. The meso-mechanically informed anisotropic damage of effective Cosserat continuum is characterized and identified and the microscopic mechanisms of macroscopic damage phenomenon are revealed. c 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi： 10.1063/2.1301101]
