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.展开更多
Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999～2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions...Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999～2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions has been determined. The results show that the pattern of velocity anomalies in the shallower upper crust is somewhat associated with the surface geological tectonics in the region. A relative low-velocity anomaly appears north of the Wangwu-Wenjiao fault zone and a relative high-velocity anomaly appears south of the Wangwu-Wenjiao fault zone,corresponding to the depressed areas in north Hainan Island,where many volcanoes are frequently active and geothermal values are relatively higher,and the uplifted and stable regions in central and south of the Hainan Is-land. In the middle and lower crust velocities are relatively lower in east Hainan than those in west Hainan,possi-bly suggesting the existence of the upwelling of hot materials from the mantle in east Hainan. The pattern of veloc-ity anomalies also indicates that NW faults,i.e.,the Puqian-Qinglan fault,may be shallower,while the E-W Wangwu-Wenjiao fault may be deeper,which perhaps extends down to Moho depth or deeper.展开更多
Passive bionic feet,known for their human-like compliance,have garnered attention for their potential to achieve notable environmental adaptability.In this paper,a method was proposed to unifying passive bionic feet s...Passive bionic feet,known for their human-like compliance,have garnered attention for their potential to achieve notable environmental adaptability.In this paper,a method was proposed to unifying passive bionic feet static supporting stability and dynamic terrain adaptability through the utilization of the Rigid-Elastic Hybrid(REH)dynamics model.First,a bionic foot model,named the Hinge Tension Elastic Complex(HTEC)model,was developed by extracting key features from human feet.Furthermore,the kinematics and REH dynamics of the HTEC model were established.Based on the foot dynamics,a nonlinear optimization method for stiffness matching(NOSM)was designed.Finally,the HTEC-based foot was constructed and applied onto BHR-B2 humanoid robot.The foot static stability is achieved.The enhanced adaptability is observed as the robot traverses square steel,lawn,and cobblestone terrains.Through proposed design method and structure,the mobility of the humanoid robot is improved.展开更多
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.展开更多
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.展开更多
Recent advances in additive manufacturing have enabled the construction of metallic lattice structures with tailored mechanical and functional properties.One potential application of metallic lattice struc-tures is in...Recent advances in additive manufacturing have enabled the construction of metallic lattice structures with tailored mechanical and functional properties.One potential application of metallic lattice struc-tures is in the impact load mitigation where an external kinetic energy is absorbed by the deformation/crushing of lattice cells.This has motivated a growing number of experimental and numerical studies,recently,on the crushing behavior of additively produced lattice structures.The present study overviews the dynamic and quasi-static crushing behavior of additively produced Ti64,316L,and AlSiMg alloy lattice structures.The first part of the study summarizes the main features of two most commonly used additive processing techniques for lattice structures,namely selective-laser-melt(SLM)and electro-beam-melt(EBM),along with a description of commonly observed process induced defects.In the second part,the deformation and strain rate sensitivities of the selected alloy lattices are outlined together with the most widely used dynamic test methods,followed by a part on the observed micro-structures of the SLM and EBM-processed Ti64,316L and AlSiMg alloys.Finally,the experimental and numerical studies on the quasi-static and dynamic compression behavior of the additively processed Ti64,316L,and AlSiMg alloy lattices are reviewed.The results of the experimental and numerical studies of the dynamic properties of various types of lattices,including graded,non-uniform strut size,hollow,non-uniform cell size,and bio-inspired,were tabulated together with the used dynamic testing methods.The dynamic tests have been noted to be mostly conducted in compression Split Hopkinson Pressure Bar(SHPB)or Taylor-and direct-impact tests using the SHPB set-up,in all of which relatively small-size test specimens were tested.The test specimen size effect on the compression behavior of the lattices was further emphasized.It has also been shown that the lattices of Ti64 and AlSiMg alloys are relatively brittle as compared with the lattices of 316L alloy.Finally,the challenges associated with modelling lattice structures were explained and the micro tension tests and multi-scale modeling techniques combining microstructural characteristics with macroscopic lattice dynamics were recommended to improve the accuracy of the numerical simulations of the dynamic compression deformations of metallic lattice structures.展开更多
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.展开更多
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.展开更多
Artificial bone with porous structure is crucial for tissue scaffold and clinic implants.Scaffold provides structure support for cells and guides tissues regeneration for final tissue structure.A computational aided p...Artificial bone with porous structure is crucial for tissue scaffold and clinic implants.Scaffold provides structure support for cells and guides tissues regeneration for final tissue structure.A computational aided process of porous bone modeling was developed which described the design and fabrication of tissue scaffolds by considering intricate architecture,porosity and pore size.To simulate intricate bone structure,different constructive units were presented.In modeling process,bone contour was gotten from computed tomography(CT)images and was divided into two levels.Each level was represented by relatively reconstructive process.Pore size distribution was controlled by using mesh generation.The whole hexahedral mesh was reduced by unit structure,when a 3D mesh with various hexahedral elements was provided.The simulation results show that constructive structure of porous scaffold can meet the needs of clinic implants in accurate and controlled way.展开更多
There have been a great demand for a suitable and convenient method in the field of buckling analysis of stiffened ship structures, which is essential to structural safety assessment and is significantly time-consumin...There have been a great demand for a suitable and convenient method in the field of buckling analysis of stiffened ship structures, which is essential to structural safety assessment and is significantly time-consuming. Modeling, buckling behaviors and ultimate strength prediction of stiffened panels were investigated. The modeling specification including nonlinear finite element model and imperfections generation, and post-buckling analysis procedure of stiffened plates were demonstrated. And a software tool using set-based finite element method was developed and executed in the MSC. Marc environment. Different types of stiffen panels of marine structures have been employed to investigate the buckling behavior and assess the validity in the estimation of ultimate strength. A comparison between results of the generally accepted methods, experiments and the software tool developed was demonstrated. It is shown that the software tool can predict the ultimate capacity of stiffened panels with imperfections with a good accuracy.展开更多
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.展开更多
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.展开更多
In order to establish the baseline finite element model for structural health monitoring,a new method of model updating was proposed after analyzing the uncertainties of measured data and the error of finite element m...In order to establish the baseline finite element model for structural health monitoring,a new method of model updating was proposed after analyzing the uncertainties of measured data and the error of finite element model.In the new method,the finite element model was replaced by the multi-output support vector regression machine(MSVR).The interval variables of the measured frequency were sampled by Latin hypercube sampling method.The samples of frequency were regarded as the inputs of the trained MSVR.The outputs of MSVR were the target values of design parameters.The steel structure of National Aquatic Center for Beijing Olympic Games was introduced as a case for finite element model updating.The results show that the proposed method can avoid solving the problem of complicated calculation.Both the estimated values and associated uncertainties of the structure parameters can be obtained by the method.The static and dynamic characteristics of the updated finite element model are in good agreement with the measured data.展开更多
Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coa...Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of展开更多
The problem considered here is to assess the Bayesian influence on the unknown param-eter matrices in a grwoth curve model with the general covariance structure.Under the non-information prior distribution assumption,...The problem considered here is to assess the Bayesian influence on the unknown param-eter matrices in a grwoth curve model with the general covariance structure.Under the non-information prior distribution assumption,the Kullback-Leibler divergence is employed to eval-uate the effect of a designated response matrix on the posterior distribution of the parameter matrix.展开更多
To examine the similarities and differences in the evolution of cavity,wetting and dynamics of a highspeed,oblique water-entry projectile with different positive angles of attack,a comparative analysis has been conduc...To examine the similarities and differences in the evolution of cavity,wetting and dynamics of a highspeed,oblique water-entry projectile with different positive angles of attack,a comparative analysis has been conducted based on the numerical results of two mathematical models,the rigid-body model and fluid-structure interaction model.In addition,the applicable scope of the above two methods,and the structural response characteristics of the projectile have also been investigated.Our results demonstrate that:(1) The impact loads and angular motion of the projectile of the rigid-body method are more likely to exhibit periodic variations due to the periodic tail slap,its range of positive angles of attack is about α<2°.(2) When the projectile undergone significant wetting,a strong coupling effect is observed among wetting,structural deformation,and projectile motion.With the applied projectile shape,it is observed that,when the projectile bends,the final wetting position is that of Part B(cylinder of body).With the occu rrence of this phenomenon,the projectile ballistics beco me completely unstable.(3) The force exerted on the lower surface of the projectile induced by wetting is the primary reason of the destabilization of the projectile traj ectory and structu ral deformation failure.Bending deformation is most likely to appear at the junction of Part C(cone of body) and Part D(tail).The safe angles of attack of the projectile stability are found to be about α≤2°.展开更多
Marine sediment velocity structural models have strong regional characteristics.Hamilton made two shallow continental shelf sediment velocity structure models,Lu Bo gave a model in accordance with the characteristics ...Marine sediment velocity structural models have strong regional characteristics.Hamilton made two shallow continental shelf sediment velocity structure models,Lu Bo gave a model in accordance with the characteristics of the continental shelf of China.However,no model can contain all geological situations.We got the in-situ velocity data at Zhapu and Jintang near the Hangzhou Bay by using the MFI GeoA(Multi-Frequency In-situ Geoacoustic Measurement),and used these data to make the velocity structure models.Finally,we got two different models.One is Zhapu velocity structural model that we can describe as Lower velocity-Higher velocity -Lower velocity-Higher velocity model simply。展开更多
This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume e...This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume element(RVE)is defined and the bridging model is adopted to establish a theoretical predictive model for its three-dimensional equivalent elastic constants.The results obtained through this method for the previous experimental model are compared with the ones gained respectively by experiments and classical laminate theory to verify the reliability of this model.In addition,the effects of some winding parameters,such as winding angle,on the equivalent elastic behavior of the filament-wound composites are analyzed.The rules gained can provide a theoretical reference for the optimum design of filament-wound composites.展开更多
The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial ...The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.展开更多
This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi...This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.展开更多
基金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.
基金The special project of Detection of Haikou City Earthquake Active Faults from the Tenth Five-year Plan of China Earthquake Administration (0106512)Joint Seismological Science Foundation of China (105086)CAS Key Laboratory of Marginal Sea Geology (MSGL0503).
文摘Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999～2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions has been determined. The results show that the pattern of velocity anomalies in the shallower upper crust is somewhat associated with the surface geological tectonics in the region. A relative low-velocity anomaly appears north of the Wangwu-Wenjiao fault zone and a relative high-velocity anomaly appears south of the Wangwu-Wenjiao fault zone,corresponding to the depressed areas in north Hainan Island,where many volcanoes are frequently active and geothermal values are relatively higher,and the uplifted and stable regions in central and south of the Hainan Is-land. In the middle and lower crust velocities are relatively lower in east Hainan than those in west Hainan,possi-bly suggesting the existence of the upwelling of hot materials from the mantle in east Hainan. The pattern of veloc-ity anomalies also indicates that NW faults,i.e.,the Puqian-Qinglan fault,may be shallower,while the E-W Wangwu-Wenjiao fault may be deeper,which perhaps extends down to Moho depth or deeper.
基金supported by the National Natural Science Foundation of China(Grant No.62073041)the Open Fund of Laboratory of Aerospace Servo Actuation and Transmission(Grant No.LASAT-2023A04)the Fundamental Research Funds for the Central Universities(Grant Nos.2024CX06011,2024CX06079)。
文摘Passive bionic feet,known for their human-like compliance,have garnered attention for their potential to achieve notable environmental adaptability.In this paper,a method was proposed to unifying passive bionic feet static supporting stability and dynamic terrain adaptability through the utilization of the Rigid-Elastic Hybrid(REH)dynamics model.First,a bionic foot model,named the Hinge Tension Elastic Complex(HTEC)model,was developed by extracting key features from human feet.Furthermore,the kinematics and REH dynamics of the HTEC model were established.Based on the foot dynamics,a nonlinear optimization method for stiffness matching(NOSM)was designed.Finally,the HTEC-based foot was constructed and applied onto BHR-B2 humanoid robot.The foot static stability is achieved.The enhanced adaptability is observed as the robot traverses square steel,lawn,and cobblestone terrains.Through proposed design method and structure,the mobility of the humanoid robot is improved.
基金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.
基金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.
基金the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No 101034425 for the project titled A2M2TECHThe Scientific and Technological Research Council of Türkiye (TUBITAK) with grant No 120C158 for the same A2M2TECH project under the TUBITAK's 2236/B program
文摘Recent advances in additive manufacturing have enabled the construction of metallic lattice structures with tailored mechanical and functional properties.One potential application of metallic lattice struc-tures is in the impact load mitigation where an external kinetic energy is absorbed by the deformation/crushing of lattice cells.This has motivated a growing number of experimental and numerical studies,recently,on the crushing behavior of additively produced lattice structures.The present study overviews the dynamic and quasi-static crushing behavior of additively produced Ti64,316L,and AlSiMg alloy lattice structures.The first part of the study summarizes the main features of two most commonly used additive processing techniques for lattice structures,namely selective-laser-melt(SLM)and electro-beam-melt(EBM),along with a description of commonly observed process induced defects.In the second part,the deformation and strain rate sensitivities of the selected alloy lattices are outlined together with the most widely used dynamic test methods,followed by a part on the observed micro-structures of the SLM and EBM-processed Ti64,316L and AlSiMg alloys.Finally,the experimental and numerical studies on the quasi-static and dynamic compression behavior of the additively processed Ti64,316L,and AlSiMg alloy lattices are reviewed.The results of the experimental and numerical studies of the dynamic properties of various types of lattices,including graded,non-uniform strut size,hollow,non-uniform cell size,and bio-inspired,were tabulated together with the used dynamic testing methods.The dynamic tests have been noted to be mostly conducted in compression Split Hopkinson Pressure Bar(SHPB)or Taylor-and direct-impact tests using the SHPB set-up,in all of which relatively small-size test specimens were tested.The test specimen size effect on the compression behavior of the lattices was further emphasized.It has also been shown that the lattices of Ti64 and AlSiMg alloys are relatively brittle as compared with the lattices of 316L alloy.Finally,the challenges associated with modelling lattice structures were explained and the micro tension tests and multi-scale modeling techniques combining microstructural characteristics with macroscopic lattice dynamics were recommended to improve the accuracy of the numerical simulations of the dynamic compression deformations of metallic lattice structures.
基金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(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.
基金Project(2011DFB70230)supported by State International Cooperation Program of ChinaProject(N110403003)supported by Basic Research Foundation of Education Ministry of China
文摘Artificial bone with porous structure is crucial for tissue scaffold and clinic implants.Scaffold provides structure support for cells and guides tissues regeneration for final tissue structure.A computational aided process of porous bone modeling was developed which described the design and fabrication of tissue scaffolds by considering intricate architecture,porosity and pore size.To simulate intricate bone structure,different constructive units were presented.In modeling process,bone contour was gotten from computed tomography(CT)images and was divided into two levels.Each level was represented by relatively reconstructive process.Pore size distribution was controlled by using mesh generation.The whole hexahedral mesh was reduced by unit structure,when a 3D mesh with various hexahedral elements was provided.The simulation results show that constructive structure of porous scaffold can meet the needs of clinic implants in accurate and controlled way.
基金Projects(51575535,51805551)supported by the National Natural Science Foundation of ChinaProject(ZZYJKT2018-15)supported by the of State Key Laboratory of High Performance Complex Manufacturing,China+1 种基金Project(2015CX002)supported by the Innovation-driven Plan in Central South University,ChinaProject(2018BB30501)supported by the Key R&D Program of Liuzhou City,China
文摘There have been a great demand for a suitable and convenient method in the field of buckling analysis of stiffened ship structures, which is essential to structural safety assessment and is significantly time-consuming. Modeling, buckling behaviors and ultimate strength prediction of stiffened panels were investigated. The modeling specification including nonlinear finite element model and imperfections generation, and post-buckling analysis procedure of stiffened plates were demonstrated. And a software tool using set-based finite element method was developed and executed in the MSC. Marc environment. Different types of stiffen panels of marine structures have been employed to investigate the buckling behavior and assess the validity in the estimation of ultimate strength. A comparison between results of the generally accepted methods, experiments and the software tool developed was demonstrated. It is shown that the software tool can predict the ultimate capacity of stiffened panels with imperfections with a good accuracy.
基金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.
基金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(50678052) supported by the National Natural Science Foundation of China
文摘In order to establish the baseline finite element model for structural health monitoring,a new method of model updating was proposed after analyzing the uncertainties of measured data and the error of finite element model.In the new method,the finite element model was replaced by the multi-output support vector regression machine(MSVR).The interval variables of the measured frequency were sampled by Latin hypercube sampling method.The samples of frequency were regarded as the inputs of the trained MSVR.The outputs of MSVR were the target values of design parameters.The steel structure of National Aquatic Center for Beijing Olympic Games was introduced as a case for finite element model updating.The results show that the proposed method can avoid solving the problem of complicated calculation.Both the estimated values and associated uncertainties of the structure parameters can be obtained by the method.The static and dynamic characteristics of the updated finite element model are in good agreement with the measured data.
文摘Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of
文摘The problem considered here is to assess the Bayesian influence on the unknown param-eter matrices in a grwoth curve model with the general covariance structure.Under the non-information prior distribution assumption,the Kullback-Leibler divergence is employed to eval-uate the effect of a designated response matrix on the posterior distribution of the parameter matrix.
基金supported by the Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX24_0714).
文摘To examine the similarities and differences in the evolution of cavity,wetting and dynamics of a highspeed,oblique water-entry projectile with different positive angles of attack,a comparative analysis has been conducted based on the numerical results of two mathematical models,the rigid-body model and fluid-structure interaction model.In addition,the applicable scope of the above two methods,and the structural response characteristics of the projectile have also been investigated.Our results demonstrate that:(1) The impact loads and angular motion of the projectile of the rigid-body method are more likely to exhibit periodic variations due to the periodic tail slap,its range of positive angles of attack is about α<2°.(2) When the projectile undergone significant wetting,a strong coupling effect is observed among wetting,structural deformation,and projectile motion.With the applied projectile shape,it is observed that,when the projectile bends,the final wetting position is that of Part B(cylinder of body).With the occu rrence of this phenomenon,the projectile ballistics beco me completely unstable.(3) The force exerted on the lower surface of the projectile induced by wetting is the primary reason of the destabilization of the projectile traj ectory and structu ral deformation failure.Bending deformation is most likely to appear at the junction of Part C(cone of body) and Part D(tail).The safe angles of attack of the projectile stability are found to be about α≤2°.
文摘Marine sediment velocity structural models have strong regional characteristics.Hamilton made two shallow continental shelf sediment velocity structure models,Lu Bo gave a model in accordance with the characteristics of the continental shelf of China.However,no model can contain all geological situations.We got the in-situ velocity data at Zhapu and Jintang near the Hangzhou Bay by using the MFI GeoA(Multi-Frequency In-situ Geoacoustic Measurement),and used these data to make the velocity structure models.Finally,we got two different models.One is Zhapu velocity structural model that we can describe as Lower velocity-Higher velocity -Lower velocity-Higher velocity model simply。
文摘This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume element(RVE)is defined and the bridging model is adopted to establish a theoretical predictive model for its three-dimensional equivalent elastic constants.The results obtained through this method for the previous experimental model are compared with the ones gained respectively by experiments and classical laminate theory to verify the reliability of this model.In addition,the effects of some winding parameters,such as winding angle,on the equivalent elastic behavior of the filament-wound composites are analyzed.The rules gained can provide a theoretical reference for the optimum design of filament-wound composites.
基金Project(51274192)supported by National Natural Science Foundation of China
文摘The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.
基金Projects(42477162,52108347,52178371,52168046,52178321,52308383)supported by the National Natural Science Foundation of ChinaProjects(2023C03143,2022C01099,2024C01219,2022C03151)supported by the Zhejiang Key Research and Development Plan,China+6 种基金Project(LQ22E080010)supported by the Exploring Youth Project of Zhejiang Natural Science Foundation,ChinaProject(LR21E080005)supported by the Outstanding Youth Project of Natural Science Foundation of Zhejiang Province,ChinaProject(2022M712964)supported by the Postdoctoral Science Foundation of ChinaProject(2023AFB008)supported by the Natural Science Foundation of Hubei Province for Youth,ChinaProject(202203)supported by Engineering Research Centre of Rock-Soil Drilling&Excavation and Protection,Ministry of Education,ChinaProject(202305-2)supported by the Science and Technology Project of Zhejiang Provincial Communication Department,ChinaProject(2021K256)supported by the Construction Research Founds of Department of Housing and Urban-Rural Development of Zhejiang Province,China。
文摘This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.
