A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path...A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.展开更多
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.展开更多
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.展开更多
Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on ...Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.展开更多
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%.展开更多
特高压直流输电系统阀厅内部结构复杂,采用有限元(FEM)法模拟计算其内部电场时剖分困难。为此,结合某±800 k V特高压换流站阀厅的3维全模型,基于模块化与独立化的思想,提出了实体模型区域化网格剖分方法。该方法克服了复杂模型下...特高压直流输电系统阀厅内部结构复杂,采用有限元(FEM)法模拟计算其内部电场时剖分困难。为此,结合某±800 k V特高压换流站阀厅的3维全模型,基于模块化与独立化的思想,提出了实体模型区域化网格剖分方法。该方法克服了复杂模型下网格剖分困难、对计算机硬件要求高的问题,剖分后可获得质量较高的网格。针对模型单元数较多的情况,选用不完全Cholesky分解预处理共轭梯度法(ICCG)求解复杂模型的有限元方程组。最后结合区域网格剖分与子模型技术,实现了阀厅内局部区域电场强度的准确计算。结果表明:将求解域尺寸设置为内部模型尺寸的1.3倍以上时,所得到的切割边界符合子模型计算要求,此时边界上的最大计算误差≤2%,子模型电场计算精度满足要求。展开更多
基金supported by the National Science Fund for Distinguished Young Scholars(52425211)BIT Research Fund Program for Young Scholars(XSQD-202201005).
文摘A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.
基金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.
基金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(51008188)supported by National Natural Science Foundation of ChinaProject(KLE-TJGE-B1302)supported by Key Laboratory Fund of Geotechnical and Underground Engineering of Ministry of Education,ChinaProject(SKLGDUEK1205)supported by Open Program of State Key Laboratory for Geomechanics and Deep Underground Engineering,China
文摘Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.
基金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%.
文摘特高压直流输电系统阀厅内部结构复杂,采用有限元(FEM)法模拟计算其内部电场时剖分困难。为此,结合某±800 k V特高压换流站阀厅的3维全模型,基于模块化与独立化的思想,提出了实体模型区域化网格剖分方法。该方法克服了复杂模型下网格剖分困难、对计算机硬件要求高的问题,剖分后可获得质量较高的网格。针对模型单元数较多的情况,选用不完全Cholesky分解预处理共轭梯度法(ICCG)求解复杂模型的有限元方程组。最后结合区域网格剖分与子模型技术,实现了阀厅内局部区域电场强度的准确计算。结果表明:将求解域尺寸设置为内部模型尺寸的1.3倍以上时,所得到的切割边界符合子模型计算要求,此时边界上的最大计算误差≤2%,子模型电场计算精度满足要求。
