The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
In the field of calculating the attack area of air-to-air missiles in modern air combat scenarios,the limitations of existing research,including real-time calculation,accuracy efficiency trade-off,and the absence of t...In the field of calculating the attack area of air-to-air missiles in modern air combat scenarios,the limitations of existing research,including real-time calculation,accuracy efficiency trade-off,and the absence of the three-dimensional attack area model,restrict their practical applications.To address these issues,an improved backtracking algorithm is proposed to improve calculation efficiency.A significant reduction in solution time and maintenance of accuracy in the three-dimensional attack area are achieved by using the proposed algorithm.Furthermore,the age-layered population structure genetic programming(ALPS-GP)algorithm is introduced to determine an analytical polynomial model of the three-dimensional attack area,considering real-time requirements.The accuracy of the polynomial model is enhanced through the coefficient correction using an improved gradient descent algorithm.The study reveals a remarkable combination of high accuracy and efficient real-time computation,with a mean error of 91.89 m using the analytical polynomial model of the three-dimensional attack area solved in just 10^(-4)s,thus meeting the requirements of real-time combat scenarios.展开更多
This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabri...This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabricated fragments are employed to examine the damage under blast shock waves and combined blast and fragments loading on various liquid-filled cylindrical shell structures.The test results are compared to numerical calculations and theoretical analysis for the structure's deformation,the liquid medium's movement,and the pressure waves'propagation characteristics under different liquid-filling methods.The results showed that the filling method influences the blast protection and the struc-ture's energy absorption performance.The external filling method reduces the structural deformation,and the internal filling method increases the damage effect.The gapped internal filling method improves the structure's energy absorption efficiency.The pressure wave loading on the liquid-filled cylindrical shell structure differs depending on filling methods.Explosive shock waves and high-speed fragments show a damage enhancement effect on the liquid-filled cylindrical shell structure,depending on the thickness of the internal liquid container layer.The specific impulse on the inner surface of the cylindrical shell positively correlates to the radial deformation of the cylindrical shell structure,and the external liquid layer limits the radial structural deformation.展开更多
Modern warfare demands weapons capable of penetrating substantial structures,which presents sig-nificant challenges to the reliability of the electronic devices that are crucial to the weapon's perfor-mance.Due to...Modern warfare demands weapons capable of penetrating substantial structures,which presents sig-nificant challenges to the reliability of the electronic devices that are crucial to the weapon's perfor-mance.Due to miniaturization of electronic components,it is challenging to directly measure or numerically predict the mechanical response of small-sized critical interconnections in board-level packaging structures to ensure the mechanical reliability of electronic devices in projectiles under harsh working conditions.To address this issue,an indirect measurement method using the Bayesian regularization-based load identification was proposed in this study based on finite element(FE)pre-dictions to estimate the load applied on critical interconnections of board-level packaging structures during the process of projectile penetration.For predicting the high-strain-rate penetration process,an FE model was established with elasto-plastic constitutive models of the representative packaging ma-terials(that is,solder material and epoxy molding compound)in which material constitutive parameters were calibrated against the experimental results by using the split-Hopkinson pressure bar.As the impact-induced dynamic bending of the printed circuit board resulted in an alternating tensile-compressive loading on the solder joints during penetration,the corner solder joints in the edge re-gions experience the highest S11 and strain,making them more prone to failure.Based on FE predictions at different structural scales,an improved Bayesian method based on augmented Tikhonov regulariza-tion was theoretically proposed to address the issues of ill-posed matrix inversion and noise sensitivity in the load identification at the critical solder joints.By incorporating a wavelet thresholding technique,the method resolves the problem of poor load identification accuracy at high noise levels.The proposed method achieves satisfactorily small relative errors and high correlation coefficients in identifying the mechanical response of local interconnections in board-level packaging structures,while significantly balancing the smoothness of response curves with the accuracy of peak identification.At medium and low noise levels,the relative error is less than 6%,while it is less than 10%at high noise levels.The proposed method provides an effective indirect approach for the boundary conditions of localized solder joints during the projectile penetration process,and its philosophy can be readily extended to other scenarios of multiscale analysis for highly nonlinear materials and structures under extreme loading conditions.展开更多
A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂F...A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
The searching method of failure surface which consists of complex geological structures in high and steep rock slopes was studied. Based on computer simulation technology and Monte-Carlo method, three dimensional mult...The searching method of failure surface which consists of complex geological structures in high and steep rock slopes was studied. Based on computer simulation technology and Monte-Carlo method, three dimensional multi-scale geological structures such as engineering scale and statistical scale structures of the slope were simulated. The searching method of failure route which consists of joints and rock bridges was determined via simulation annealing method by considering the shear strength of joints or rock bridges in one supposed route. When shear strengths of all the supposed routes were computed, the least shear strength route was considered failure route. Then, the inclined slice of joint slices and rock bridge slices were separated according to the position of joints and rock bridges. For the rock bridge slices, by distinguishing the failure model, the force direction to the next slice was defined. Finally, the limit equilibrium equations for every slice were established, and the slope stability factor was obtained. One practical example indicates that the discussed method is more closely to the real condition.展开更多
In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite ...In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.展开更多
The static performance of inflatable structures has been well studied and the dynamic deployment simulation has received much attention. However, very few studies focus on its deflation behavior. Although there are se...The static performance of inflatable structures has been well studied and the dynamic deployment simulation has received much attention. However, very few studies focus on its deflation behavior. Although there are several dynamic finite element algorithms that can be applied to the deflation simulation, their computation costs are expensive, especially for large scale structures. In this work, a simple method based on classic thermodynamics and the analytical relationship between air and membrane was proposed to efficiently analyze the air state variables under the condition of ventilation. Combined with failure analysis of static bearing capacity, a fast incremental analytical method was presented to predict both elastic and post wrinkling deflation process of inflatable structures. Comparisons between simplified analysis, dynamic finite element simulation, and a full-scale experimental test are presented and the suitability of this simple method for solving the air state and predicting the deflation behavior of inflatable structures is proved.展开更多
In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this wo...Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.展开更多
In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation...In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.展开更多
The dynamic characteristic analysis model of antenna structures is built,in which the structural physical parameters and geometrical dimensions are all considered as unascertained variables.And a structure dynamic cha...The dynamic characteristic analysis model of antenna structures is built,in which the structural physical parameters and geometrical dimensions are all considered as unascertained variables.And a structure dynamic characteristic analysis method based on the unascertained factor method is given.The computational expression of structural characteristic is developed by the mathematics expression of unascertained factor and the principles of unascertained rational numbers arithmetic.An example is given,in which the possible values and confidence degrees of the unascertained structure characteristics are obtained.The calculated results show that the method is feasible and effective.展开更多
For the 110 mining method,it is challenging to accurately calculate the support resistance of the roadway due to the lack of understanding of the dynamic movement of the overlying strata in this method.The consequenti...For the 110 mining method,it is challenging to accurately calculate the support resistance of the roadway due to the lack of understanding of the dynamic movement of the overlying strata in this method.The consequential excessive support results in a significant increase in the cost of roadway support.The authors explored the overlying strata movement and roadway deformation of the gob-entry retaining in the 110 mining method to solve this problem.First,the typical stages of the roof-cutting gob-side entry were defined.Second,the mechanical model and calculation formula of the support resistance on the roof were explored.Then,using numerical simulation software,the starting ranges of the specific supports at different stages were verified and the feasibility of the support scheme was examined.Finally,combined with the field measurement data,the stress and the deformation of the gob roadway at different stages under the influence of two mining processes in the 110 mining method were obtained.The numerical simulation results obtained are consistent with the field test results,providing a theoretical basis for precision support at different stages by the 110 mining method.展开更多
In order to optimize the crashworthy characteristic of energy-absorbing structures, the surrogate models of specific energy absorption （SEA） and ratio of SEA to initial peak force （REAF） with respect to the design...In order to optimize the crashworthy characteristic of energy-absorbing structures, the surrogate models of specific energy absorption （SEA） and ratio of SEA to initial peak force （REAF） with respect to the design parameters were respectively constructed based on surrogate model optimization methods （polynomial response surface method （PRSM） and Kriging method （KM））. Firstly, the sample data were prepared through the design of experiment （DOE）. Then, the test data models were set up based on the theory of surrogate model, and the data samples were trained to obtain the response relationship between the SEA ＆amp; REAF and design parameters. At last, the structure optimal parameters were obtained by visual analysis and genetic algorithm （GA）. The results indicate that the KM, where the local interpolation method is used in Gauss correlation function, has the highest fitting accuracy and the structure optimal parameters are obtained as： the SEA of 29.8558 kJ/kg （corresponding toa=70 mm andt= 3.5 mm） and REAF of 0.2896 （corresponding toa=70 mm andt=1.9615 mm）. The basis function of the quartic PRSM with higher order than that of the quadratic PRSM, and the mutual influence of the design variables are considered, so the fitting accuracy of the quartic PRSM is higher than that of the quadratic PRSM.展开更多
Finite element modeling methods of steel-concrete composite structure with overlap slab were investigated. A two-step finite element method was presented. It was applied to analyze an extra long span composite bridge....Finite element modeling methods of steel-concrete composite structure with overlap slab were investigated. A two-step finite element method was presented. It was applied to analyze an extra long span composite bridge. The conversion of structure system and the mechanical behavior of the bridge were analyzed with two different construction methods. The stresses of steel beams, precast slabs and in-situ-place concrete under the total load were compared. The results show that steel-concrete composite structure with overlap slab has many advantages, the construction method that the top in-situ concrete and the concrete in construction joints are cast respectively is rather reasonable than the one that the top in-situ concrete and the concrete in construction joints are cast at the same time, and the two-step finite element method is affective to such large-scale structures.展开更多
Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing me...Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing mechanism of fluid structure interaction between circulation valve and liquid of hydraulic shock absorber. The fluid mesh distortion was controlled by the CEL language, and the fluid struc^tre interaction mathematical model was established. The finite element model was established by ANSYS CFX software and was analyzed by dynamic mesh technique. The local sensitive computational area was meshed by prismatic grid, which could reduce the negative volume problem during the simulation. The circulation valve and liquid of hydraulic shock absorber were simulated and analyzed under the condition of sinusoidal inlet velocity loads. Flow characteristic and dynamics characteristic were obtained. The pressure distribution and the displacement of circulation value were obtained, and the acceleration curve of circulation valve was simulated and analyzed. The conformity of the final simulation results with the experimental datum indicates that this method is accurate and reliable to analyze the dynamics characteristic between circulation valve and liquid of hydraulic shock absorber, which can provide a theoretical foundation for optimizing hydraulic shock absorber in the future.展开更多
The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathe...The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.展开更多
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金National Natural Science Foundation of China(62373187)Forward-looking Layout Special Projects(ILA220591A22)。
文摘In the field of calculating the attack area of air-to-air missiles in modern air combat scenarios,the limitations of existing research,including real-time calculation,accuracy efficiency trade-off,and the absence of the three-dimensional attack area model,restrict their practical applications.To address these issues,an improved backtracking algorithm is proposed to improve calculation efficiency.A significant reduction in solution time and maintenance of accuracy in the three-dimensional attack area are achieved by using the proposed algorithm.Furthermore,the age-layered population structure genetic programming(ALPS-GP)algorithm is introduced to determine an analytical polynomial model of the three-dimensional attack area,considering real-time requirements.The accuracy of the polynomial model is enhanced through the coefficient correction using an improved gradient descent algorithm.The study reveals a remarkable combination of high accuracy and efficient real-time computation,with a mean error of 91.89 m using the analytical polynomial model of the three-dimensional attack area solved in just 10^(-4)s,thus meeting the requirements of real-time combat scenarios.
基金the National Natural Science Foundation of China(Grant Nos.52371342,52271338,52101378 and 51979277)。
文摘This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabricated fragments are employed to examine the damage under blast shock waves and combined blast and fragments loading on various liquid-filled cylindrical shell structures.The test results are compared to numerical calculations and theoretical analysis for the structure's deformation,the liquid medium's movement,and the pressure waves'propagation characteristics under different liquid-filling methods.The results showed that the filling method influences the blast protection and the struc-ture's energy absorption performance.The external filling method reduces the structural deformation,and the internal filling method increases the damage effect.The gapped internal filling method improves the structure's energy absorption efficiency.The pressure wave loading on the liquid-filled cylindrical shell structure differs depending on filling methods.Explosive shock waves and high-speed fragments show a damage enhancement effect on the liquid-filled cylindrical shell structure,depending on the thickness of the internal liquid container layer.The specific impulse on the inner surface of the cylindrical shell positively correlates to the radial deformation of the cylindrical shell structure,and the external liquid layer limits the radial structural deformation.
基金supported by the National Natural Science Foundation of China(Grant Nos.52475166,52175148)the Regional Collaboration Project of Shanxi Province(Grant No.202204041101044).
文摘Modern warfare demands weapons capable of penetrating substantial structures,which presents sig-nificant challenges to the reliability of the electronic devices that are crucial to the weapon's perfor-mance.Due to miniaturization of electronic components,it is challenging to directly measure or numerically predict the mechanical response of small-sized critical interconnections in board-level packaging structures to ensure the mechanical reliability of electronic devices in projectiles under harsh working conditions.To address this issue,an indirect measurement method using the Bayesian regularization-based load identification was proposed in this study based on finite element(FE)pre-dictions to estimate the load applied on critical interconnections of board-level packaging structures during the process of projectile penetration.For predicting the high-strain-rate penetration process,an FE model was established with elasto-plastic constitutive models of the representative packaging ma-terials(that is,solder material and epoxy molding compound)in which material constitutive parameters were calibrated against the experimental results by using the split-Hopkinson pressure bar.As the impact-induced dynamic bending of the printed circuit board resulted in an alternating tensile-compressive loading on the solder joints during penetration,the corner solder joints in the edge re-gions experience the highest S11 and strain,making them more prone to failure.Based on FE predictions at different structural scales,an improved Bayesian method based on augmented Tikhonov regulariza-tion was theoretically proposed to address the issues of ill-posed matrix inversion and noise sensitivity in the load identification at the critical solder joints.By incorporating a wavelet thresholding technique,the method resolves the problem of poor load identification accuracy at high noise levels.The proposed method achieves satisfactorily small relative errors and high correlation coefficients in identifying the mechanical response of local interconnections in board-level packaging structures,while significantly balancing the smoothness of response curves with the accuracy of peak identification.At medium and low noise levels,the relative error is less than 6%,while it is less than 10%at high noise levels.The proposed method provides an effective indirect approach for the boundary conditions of localized solder joints during the projectile penetration process,and its philosophy can be readily extended to other scenarios of multiscale analysis for highly nonlinear materials and structures under extreme loading conditions.
基金Project(51208178)supported by the National Natural Science Foundation of ChinaProject(2012M520991)supported by China Postdoctoral Science Foundation
文摘A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
基金Project(50539100) supported by the National Natural Science Foundation of ChinaProject(BK2006171) supported by the Jiangsu Natural Science Foundation
文摘The searching method of failure surface which consists of complex geological structures in high and steep rock slopes was studied. Based on computer simulation technology and Monte-Carlo method, three dimensional multi-scale geological structures such as engineering scale and statistical scale structures of the slope were simulated. The searching method of failure route which consists of joints and rock bridges was determined via simulation annealing method by considering the shear strength of joints or rock bridges in one supposed route. When shear strengths of all the supposed routes were computed, the least shear strength route was considered failure route. Then, the inclined slice of joint slices and rock bridge slices were separated according to the position of joints and rock bridges. For the rock bridge slices, by distinguishing the failure model, the force direction to the next slice was defined. Finally, the limit equilibrium equations for every slice were established, and the slope stability factor was obtained. One practical example indicates that the discussed method is more closely to the real condition.
基金the National Natural Science Foundation of China (5963140) Doctor Point Fund of National Education Committee Parent Company Fund of Aviation Industry
文摘In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.
基金Projects(51178263,51378307)supported by the National Natural Science Foundation of China
文摘The static performance of inflatable structures has been well studied and the dynamic deployment simulation has received much attention. However, very few studies focus on its deflation behavior. Although there are several dynamic finite element algorithms that can be applied to the deflation simulation, their computation costs are expensive, especially for large scale structures. In this work, a simple method based on classic thermodynamics and the analytical relationship between air and membrane was proposed to efficiently analyze the air state variables under the condition of ventilation. Combined with failure analysis of static bearing capacity, a fast incremental analytical method was presented to predict both elastic and post wrinkling deflation process of inflatable structures. Comparisons between simplified analysis, dynamic finite element simulation, and a full-scale experimental test are presented and the suitability of this simple method for solving the air state and predicting the deflation behavior of inflatable structures is proved.
基金Supported by the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
基金Project(2011CB013504) supported by the National Basic Research Program(973 Program)of ChinaProject(2013BAB06B01) supported by the National Science&Technology Pillar Program during the Twelfth Five-year Plan Period+2 种基金Projects(11772118,51479049,51709282) supported by the National Natural Science Foundation of ChinaProject(2017M620838) supported by the Postdoctoral Science Foundation of ChinaProject(487237) supported by the Natural Sciences and Engineering Research Council of Canada
文摘Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.
文摘In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.
基金the National Defense Science and Technology Research Projects of China (51421060505DZ0155)the National Science Foundation of Shaanxi Province of China (2005A009)
文摘The dynamic characteristic analysis model of antenna structures is built,in which the structural physical parameters and geometrical dimensions are all considered as unascertained variables.And a structure dynamic characteristic analysis method based on the unascertained factor method is given.The computational expression of structural characteristic is developed by the mathematics expression of unascertained factor and the principles of unascertained rational numbers arithmetic.An example is given,in which the possible values and confidence degrees of the unascertained structure characteristics are obtained.The calculated results show that the method is feasible and effective.
基金Project(51674265) supported by the National Natural Science Foundation of ChinaProjects(2018YFC0603705,2016YFC0600901) supported by the State Key Research Development Program of ChinaProject supported by the Yueqi Outstanding Scholar Award Program of China University of Mining&Technology,Beijing,China。
文摘For the 110 mining method,it is challenging to accurately calculate the support resistance of the roadway due to the lack of understanding of the dynamic movement of the overlying strata in this method.The consequential excessive support results in a significant increase in the cost of roadway support.The authors explored the overlying strata movement and roadway deformation of the gob-entry retaining in the 110 mining method to solve this problem.First,the typical stages of the roof-cutting gob-side entry were defined.Second,the mechanical model and calculation formula of the support resistance on the roof were explored.Then,using numerical simulation software,the starting ranges of the specific supports at different stages were verified and the feasibility of the support scheme was examined.Finally,combined with the field measurement data,the stress and the deformation of the gob roadway at different stages under the influence of two mining processes in the 110 mining method were obtained.The numerical simulation results obtained are consistent with the field test results,providing a theoretical basis for precision support at different stages by the 110 mining method.
基金Project(U1334208)supported by the National Natural Science Foundation of ChinaProject(2013GK2001)supported by the Fund of Hunan Provincial Science and Technology Department,China
文摘In order to optimize the crashworthy characteristic of energy-absorbing structures, the surrogate models of specific energy absorption （SEA） and ratio of SEA to initial peak force （REAF） with respect to the design parameters were respectively constructed based on surrogate model optimization methods （polynomial response surface method （PRSM） and Kriging method （KM））. Firstly, the sample data were prepared through the design of experiment （DOE）. Then, the test data models were set up based on the theory of surrogate model, and the data samples were trained to obtain the response relationship between the SEA ＆amp; REAF and design parameters. At last, the structure optimal parameters were obtained by visual analysis and genetic algorithm （GA）. The results indicate that the KM, where the local interpolation method is used in Gauss correlation function, has the highest fitting accuracy and the structure optimal parameters are obtained as： the SEA of 29.8558 kJ/kg （corresponding toa=70 mm andt= 3.5 mm） and REAF of 0.2896 （corresponding toa=70 mm andt=1.9615 mm）. The basis function of the quartic PRSM with higher order than that of the quadratic PRSM, and the mutual influence of the design variables are considered, so the fitting accuracy of the quartic PRSM is higher than that of the quadratic PRSM.
基金Project(2001G022) supported by the Science and Technology Development Program of Railways Department, China
文摘Finite element modeling methods of steel-concrete composite structure with overlap slab were investigated. A two-step finite element method was presented. It was applied to analyze an extra long span composite bridge. The conversion of structure system and the mechanical behavior of the bridge were analyzed with two different construction methods. The stresses of steel beams, precast slabs and in-situ-place concrete under the total load were compared. The results show that steel-concrete composite structure with overlap slab has many advantages, the construction method that the top in-situ concrete and the concrete in construction joints are cast respectively is rather reasonable than the one that the top in-situ concrete and the concrete in construction joints are cast at the same time, and the two-step finite element method is affective to such large-scale structures.
基金Project(51275542) supported by the National Natural Science Foundation of Chinaproject(CDJXS12110010) supported by the Fundamental Research Funds for the Central Universities of China
文摘Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing mechanism of fluid structure interaction between circulation valve and liquid of hydraulic shock absorber. The fluid mesh distortion was controlled by the CEL language, and the fluid struc^tre interaction mathematical model was established. The finite element model was established by ANSYS CFX software and was analyzed by dynamic mesh technique. The local sensitive computational area was meshed by prismatic grid, which could reduce the negative volume problem during the simulation. The circulation valve and liquid of hydraulic shock absorber were simulated and analyzed under the condition of sinusoidal inlet velocity loads. Flow characteristic and dynamics characteristic were obtained. The pressure distribution and the displacement of circulation value were obtained, and the acceleration curve of circulation valve was simulated and analyzed. The conformity of the final simulation results with the experimental datum indicates that this method is accurate and reliable to analyze the dynamics characteristic between circulation valve and liquid of hydraulic shock absorber, which can provide a theoretical foundation for optimizing hydraulic shock absorber in the future.
基金Projects(51409167,51139001,51179066)supported by the National Natural Science Foundation of ChinaProjects(201401022,201501036)supported by the Ministry of Water Resources Public Welfare Industry Research Special Fund,ChinaProjects(GG201532,GG201546)supported by the Scientific and Technological Research for Water Conservancy,Henan Province,China
文摘The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.
