To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondar...To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.展开更多
Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade compone...Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.展开更多
Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to en...Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.展开更多
To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets ...To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.展开更多
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 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.展开更多
In mammals,the ribbon synapses of cochlear inner hair cells are a synaptic structure of the first sensory nerve in the pathway of acoustical signal transmission to the acoustic center,and it is directly involved in vo...In mammals,the ribbon synapses of cochlear inner hair cells are a synaptic structure of the first sensory nerve in the pathway of acoustical signal transmission to the acoustic center,and it is directly involved in voice coding and neurotransmitter release.It is difficult to quantitatively analyze the ribbon synaptic number only using an electron microscope,because the ribbon synaptic number is relatively limited and their location is deep.In this study,the specific presynaptic structure-RIBEYE,and non-specific postsynaptic structure-GluR 2&3 in C57BL/6J mouse basilar membrane samples were treated by immunofluorescent labeling.Serial section was performed on the samples using a laser scanning confocal microscope,and then the serial sections were used to build three-dimensional models using 3DS MAX software.Each fluorescein color pair indicates one synapse,so the number of ribbon synapses of inner hair cells is obtained.The spatial distribution and the number of ribbon synapses of cochlear inner hair cells were clearly shown in this experiment,and the mean number of ribbon synapses per inner hair cell was 16.10±1.03.Our results have demonstrated the number of ribbon synapses is accurately calculated by double immunofluorescent labeling to presynaptic and postsynaptic structures,serial sections obtained using a laser scanning confocal microscope,and three-dimensional modeling obtained using 3DS MAX software.The method above is feasible and has important significance for further exploring the mechanism of sensorineural deafness.展开更多
In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteri...In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.展开更多
The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this p...The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.展开更多
A three-dimensional （3-D） approach based on the state space method is proposed to study size-dependent mechanical properties of ultra-thin plate-like elastic structures considering surface effects. The structure is ...A three-dimensional （3-D） approach based on the state space method is proposed to study size-dependent mechanical properties of ultra-thin plate-like elastic structures considering surface effects. The structure is modeled as a laminate composed of a bulk bounded with upper and bottom surface layers, which are allowed to have different material properties from the bulk layer. State equations, including the surface properties of the structure, can be established on the basis of 3-D fundamental elasticity to analyze the size-dependent static characteristics of the thin plate-like structure. Compared with two-dimensional plate theories based size-dependent models for thin film structures in literature, the present 3-D approach is exact, which can provide benchmark results to assess the accuracy of 2-D plate theories and various numerical approaches. To show the feasibility of the proposed approach, a 3-D analytical solution for a simply supported plate-like thin structure including surface layers is derived. An algorithm is proposed for the calculation of the state equations obtained to ensure that the numerical results can reveal the surface effects clearly even for extremely thin surface layers. Numerical examples are carried out to exhibit the surface effects and some discussions are provided based on the results obtained.展开更多
With the development of geophysical exploration technology,the anisotropy of underground media has got more and more attention.At present,there are few studies on the anisotropy of the induced polarization method.This...With the development of geophysical exploration technology,the anisotropy of underground media has got more and more attention.At present,there are few studies on the anisotropy of the induced polarization method.This article explores the effect of anisotropy on the underground media of the induced polarization method under three-dimensional complex terrain.The research work transforms the underground electric field control equation into a variational problem,and use the unstructured finite element method to construct a large linear equation system for solving electric potentials.By the sparse matrix compression technique and symmetric successive over-relaxation preconditioned conjugate gradient algorithm(SSOR-PCG)to solve the equation system.Finally,the article uses the classic central gradient array method to obtain the forward apparent polarizability value.The calculation results of the model find that different anisotropic conditions will significantly affect the forward results which show a strong directional correlation,revealing the importance of considering anisotropy in practical work.展开更多
According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing ...According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.展开更多
Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages ha...Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages have difficulties in predicting comprehensively Reynolds number effects on airfoils,matching and characteristics curves.This study proposes Re-correction models for loss,deviation angle and endwall blockage based on classical theories and cascade tests,and loss and deviation models show good agreement in test data of NACA65 and C4 cascades.Throughflow method considering Reynolds number effects is developed by integrating the correction models into a verified Streamline Curvature(SLC)tool.A three-stage axial compressor is investigated through SLC and CFD methods from design Reynolds number(Red=2106)to low Re=4104,and the numerical methods are validated with test data of characteristic curves and spanwise distributions at Red.With Re reduction,SLC method with correction models well predicts variation in overall performances compared with CFD calculations and Wassell's model.Streamwise and spanwise matching such as total pressure and loss distributions in SLC predictions are basically consistent with those in CFD results at near-stall points under design and low Reynolds numbers.SLC and CFD methods share similar detections of stall risks in the third stage(Stg3),and their analyses of diffusion processes deviate to some extent due to different predictions in separated endwall flow.The correction models can be adopted to consider Reynolds number effects in through-flow design and analysis of axial compressors.展开更多
The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.Howeve...The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.展开更多
This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃...This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.展开更多
The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying ge...The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.展开更多
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.展开更多
Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to b...Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to be considered a gold standard.This study aims to introduce novel metrics to differentiate between KCN and healthy corneas using three-dimensional(3D)measurements of surface area and volume.Methods:This retrospective observational study examined KCN patients along with healthy control patients between the ages of 20 and 79 years old at the University of Maryland,Baltimore.The selected patients underwent a nine-line raster scan anterior segment optical coherence tomography(AS-OCT).ImageJ was used to determine the central 6 mm of each image and each corneal image was then divided into six 1 mm segments.Free-D software was then used to render the nine different images into a 3D model to calculate corneal surface area and volume.A two-tailed Mann-Whitney test was used to assess statistical significance when comparing these subsets.Results:Thirty-three eyes with KCN,along with 33 healthy control,were enrolled.There were statistically significant differences between the healthy and KCN groups in the metric of anterior corneal surface area(13.927 vs.13.991 mm^(2),P=0.046),posterior corneal surface area(14.045 vs.14.173 mm^(2),P<0.001),and volume(8.430 vs.7.773 mm3,P<0.001)within the central 6 mm.Conclusions:3D corneal models derived from AS-OCT can be used to measure anterior corneal surface area,posterior corneal surface area,and corneal volume.All three parameters are statistically different between corneas with KCN and healthy corneas.Further study and application of these parameters may yield new methodologies for the detection of KCN.展开更多
Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implemen...Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.展开更多
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.
基金Project supported by the National Natural Science Foundation of China(Nos.12372071 and 12372070)the Aeronautical Science Fund of China(No.2022Z055052001)the Foundation of China Scholarship Council(No.202306830079)。
文摘Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.
基金supported by the National Natural Science Foundation of China(No.92371206)the Postgraduate Scientific Research Innovation Project of Hunan Province,China(No.CX2023063).
文摘Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.
基金financially supported by the Ministry of Science and Technology of China(Nos.2022YFF0801201,2021YFC2900300)the National Natural Science Foundation of China(Nos.41872245,U1911202)the Guangdong Basic and Applied Basic Research Foundation(No.2020A1515010666)。
文摘To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.
基金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.
基金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.
文摘In mammals,the ribbon synapses of cochlear inner hair cells are a synaptic structure of the first sensory nerve in the pathway of acoustical signal transmission to the acoustic center,and it is directly involved in voice coding and neurotransmitter release.It is difficult to quantitatively analyze the ribbon synaptic number only using an electron microscope,because the ribbon synaptic number is relatively limited and their location is deep.In this study,the specific presynaptic structure-RIBEYE,and non-specific postsynaptic structure-GluR 2&3 in C57BL/6J mouse basilar membrane samples were treated by immunofluorescent labeling.Serial section was performed on the samples using a laser scanning confocal microscope,and then the serial sections were used to build three-dimensional models using 3DS MAX software.Each fluorescein color pair indicates one synapse,so the number of ribbon synapses of inner hair cells is obtained.The spatial distribution and the number of ribbon synapses of cochlear inner hair cells were clearly shown in this experiment,and the mean number of ribbon synapses per inner hair cell was 16.10±1.03.Our results have demonstrated the number of ribbon synapses is accurately calculated by double immunofluorescent labeling to presynaptic and postsynaptic structures,serial sections obtained using a laser scanning confocal microscope,and three-dimensional modeling obtained using 3DS MAX software.The method above is feasible and has important significance for further exploring the mechanism of sensorineural deafness.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos..42030806,41774125,41804098,41904104)the Key National Research Project of China(Grant No.2018YFC0603300).
文摘In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.
文摘The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.
基金supported by the Natural Science Foundation of Anhui Province(No.070414190).
文摘A three-dimensional （3-D） approach based on the state space method is proposed to study size-dependent mechanical properties of ultra-thin plate-like elastic structures considering surface effects. The structure is modeled as a laminate composed of a bulk bounded with upper and bottom surface layers, which are allowed to have different material properties from the bulk layer. State equations, including the surface properties of the structure, can be established on the basis of 3-D fundamental elasticity to analyze the size-dependent static characteristics of the thin plate-like structure. Compared with two-dimensional plate theories based size-dependent models for thin film structures in literature, the present 3-D approach is exact, which can provide benchmark results to assess the accuracy of 2-D plate theories and various numerical approaches. To show the feasibility of the proposed approach, a 3-D analytical solution for a simply supported plate-like thin structure including surface layers is derived. An algorithm is proposed for the calculation of the state equations obtained to ensure that the numerical results can reveal the surface effects clearly even for extremely thin surface layers. Numerical examples are carried out to exhibit the surface effects and some discussions are provided based on the results obtained.
文摘With the development of geophysical exploration technology,the anisotropy of underground media has got more and more attention.At present,there are few studies on the anisotropy of the induced polarization method.This article explores the effect of anisotropy on the underground media of the induced polarization method under three-dimensional complex terrain.The research work transforms the underground electric field control equation into a variational problem,and use the unstructured finite element method to construct a large linear equation system for solving electric potentials.By the sparse matrix compression technique and symmetric successive over-relaxation preconditioned conjugate gradient algorithm(SSOR-PCG)to solve the equation system.Finally,the article uses the classic central gradient array method to obtain the forward apparent polarizability value.The calculation results of the model find that different anisotropic conditions will significantly affect the forward results which show a strong directional correlation,revealing the importance of considering anisotropy in practical work.
基金support of the National Natural Science Foundation of China(No.12172023)。
文摘According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.
基金supported by the National Science and Tech-nology Major Project of China(Nos.2017-II-0007-0021 and J2019-II-0017-0038)。
文摘Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages have difficulties in predicting comprehensively Reynolds number effects on airfoils,matching and characteristics curves.This study proposes Re-correction models for loss,deviation angle and endwall blockage based on classical theories and cascade tests,and loss and deviation models show good agreement in test data of NACA65 and C4 cascades.Throughflow method considering Reynolds number effects is developed by integrating the correction models into a verified Streamline Curvature(SLC)tool.A three-stage axial compressor is investigated through SLC and CFD methods from design Reynolds number(Red=2106)to low Re=4104,and the numerical methods are validated with test data of characteristic curves and spanwise distributions at Red.With Re reduction,SLC method with correction models well predicts variation in overall performances compared with CFD calculations and Wassell's model.Streamwise and spanwise matching such as total pressure and loss distributions in SLC predictions are basically consistent with those in CFD results at near-stall points under design and low Reynolds numbers.SLC and CFD methods share similar detections of stall risks in the third stage(Stg3),and their analyses of diffusion processes deviate to some extent due to different predictions in separated endwall flow.The correction models can be adopted to consider Reynolds number effects in through-flow design and analysis of axial compressors.
基金supported by the National Natural Science Foundation of China(Grant Nos.42250103 and 42174090)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(Grant No.GLAB2023ZR02)the MOST Special Fund from the State Key Laboratory of Geological Processes and Mineral Resources(Grant No.MSFGPMR2022-4).
文摘The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.
文摘This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.
基金supported by a grant from the Research Grant Council of Hong Kong Special Administrative Region(Project No.11207724).
文摘The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.
基金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.
文摘Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to be considered a gold standard.This study aims to introduce novel metrics to differentiate between KCN and healthy corneas using three-dimensional(3D)measurements of surface area and volume.Methods:This retrospective observational study examined KCN patients along with healthy control patients between the ages of 20 and 79 years old at the University of Maryland,Baltimore.The selected patients underwent a nine-line raster scan anterior segment optical coherence tomography(AS-OCT).ImageJ was used to determine the central 6 mm of each image and each corneal image was then divided into six 1 mm segments.Free-D software was then used to render the nine different images into a 3D model to calculate corneal surface area and volume.A two-tailed Mann-Whitney test was used to assess statistical significance when comparing these subsets.Results:Thirty-three eyes with KCN,along with 33 healthy control,were enrolled.There were statistically significant differences between the healthy and KCN groups in the metric of anterior corneal surface area(13.927 vs.13.991 mm^(2),P=0.046),posterior corneal surface area(14.045 vs.14.173 mm^(2),P<0.001),and volume(8.430 vs.7.773 mm3,P<0.001)within the central 6 mm.Conclusions:3D corneal models derived from AS-OCT can be used to measure anterior corneal surface area,posterior corneal surface area,and corneal volume.All three parameters are statistically different between corneas with KCN and healthy corneas.Further study and application of these parameters may yield new methodologies for the detection of KCN.
基金supported by grants received by the first author and third author from the Institute of Eminence,Delhi University,Delhi,India,as part of the Faculty Research Program via Ref.No./IoE/2024-25/12/FRP.
文摘Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.
