This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens ...This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creep/relaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.展开更多
Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low tim...Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the ＂Shape Distribution （D2）＂ algorithm, and a new algorithm named ＂Quick D2＂ is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, ＂3D Searcher＂, a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.展开更多
3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situa...3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.展开更多
This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of com...This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.展开更多
In order to improve the accuracy and efficiency of 3D model retrieval,the method based on affinity propagation clustering algorithm is proposed. Firstly,projection ray-based method is proposed to improve the feature e...In order to improve the accuracy and efficiency of 3D model retrieval,the method based on affinity propagation clustering algorithm is proposed. Firstly,projection ray-based method is proposed to improve the feature extraction efficiency of 3D models. Based on the relationship between model and its projection,the intersection in 3D space is transformed into intersection in 2D space,which reduces the number of intersection and improves the efficiency of the extraction algorithm. In feature extraction,multi-layer spheres method is analyzed. The two-layer spheres method makes the feature vector more accurate and improves retrieval precision. Secondly,Semi-supervised Affinity Propagation ( S-AP) clustering is utilized because it can be applied to different cluster structures. The S-AP algorithm is adopted to find the center models and then the center model collection is built. During retrieval process,the collection is utilized to classify the query model into corresponding model base and then the most similar model is retrieved in the model base. Finally,75 sample models from Princeton library are selected to do the experiment and then 36 models are used for retrieval test. The results validate that the proposed method outperforms the original method and the retrieval precision and recall ratios are improved effectively.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
In this study, the petrophysical parameters such as density, sonic, neutron, and porosity were investigated and presented in the 3D models. The 3D models were built using geostatistical method that is used to estimate...In this study, the petrophysical parameters such as density, sonic, neutron, and porosity were investigated and presented in the 3D models. The 3D models were built using geostatistical method that is used to estimate studied parameters in the entire reservoir. For this purpose, the variogram of each parameter was determined to specify spatial correlation of data. Resulted variograms were non-monotonic. That shows anisotropy of structure. The lithology and porosity parameters are the main causes of this anisotropy. The 3D models also show that petrophysical data has higher variation in north part of reservoir than south part. In addition to, the west limb of reservoir shows higher porosity than east limb. The variation of sonic and neutron data are similar whereas the density data has opposed variation.展开更多
A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary...A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.展开更多
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.展开更多
The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the dat...The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.展开更多
A new horn failure mechanism was constructed for tunnel faces in the soft rock mass by means of the logarithmic spiral curve. The seismic action was incorporated into the horn failure mechanism using the pseudo-static...A new horn failure mechanism was constructed for tunnel faces in the soft rock mass by means of the logarithmic spiral curve. The seismic action was incorporated into the horn failure mechanism using the pseudo-static method. Considering the randomness of rock mass parameters and loads, a three-dimensional (3D) stochastic collapse model was established. Reliability analysis of seismic stability of tunnel faces was presented via the kinematical approach and the response surface method. The results show that, the reliability of tunnel faces is significantly affected by the supporting pressure, geological strength index, uniaxial compressive strength, rock bulk density and seismic forces. It is worth noting that, if the effect of seismic force was not considered, the stability of tunnel faces would be obviously overestimated. However, the correlation between horizontal and vertical seismic forces can be ignored under the condition of low calculation accuracy.展开更多
A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in pol...A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Experimental X-ray crystallography, NMR (Nuclear Magnetic Resonance) spectroscopy, dual polarization interferometry, etc. are indeed very powerful tools to determine the 3-Dimensional structure of a protein (including...Experimental X-ray crystallography, NMR (Nuclear Magnetic Resonance) spectroscopy, dual polarization interferometry, etc. are indeed very powerful tools to determine the 3-Dimensional structure of a protein (including the membrane protein);theoretical mathematical and physical computational approaches can also allow us to obtain a description of the protein 3D structure at a submicroscopic level for some unstable, noncrystalline and insoluble proteins. X-ray crystallography finds the X-ray final structure of a protein, which usually need refinements using theoretical protocols in order to produce a better structure. This means theoretical methods are also important in determinations of protein structures. Optimization is always needed in the computer-aided drug design, structure-based drug design, molecular dynamics, and quantum and molecular mechanics. This paper introduces some optimization algorithms used in these research fields and presents a new theoretical computational method—an improved LBFGS Quasi-Newtonian mathematical optimization method—to produce 3D structures of prion AGAAAAGA amyloid fibrils (which are unstable, noncrystalline and insoluble), from the potential energy minimization point of view. Because the NMR or X-ray structure of the hydrophobic region AGAAAAGA of prion proteins has not yet been determined, the model constructed by this paper can be used as a reference for experimental studies on this region, and may be useful in furthering the goals of medicinal chemistry in this field.展开更多
The design of the cooling system of injection molds directly affects both productivity and the quality of the final part. Using the cooling process CAE system to instruct the mold design, the efficiency and quality ...The design of the cooling system of injection molds directly affects both productivity and the quality of the final part. Using the cooling process CAE system to instruct the mold design, the efficiency and quality of design can be improved greatly. At the same time, it is helpful to confirm the cooling system structure and optimize the process conditions. In this paper, the 3D surface model of mold cavity is used to replace the middle-plane model in the simulation by Boundary Element Method, which break the bottleneck of the application of the injection molding simulation softwares base on the middle-plane model. With the improvements of this paper, a practical and commercial simulation software of injection molding cooling process named as HsCAE3D6.0 is developed.展开更多
This paper presents an experimentally validated weld joint shape and dimensions predictive 3D modeling for low carbon galvanized steel in butt-joint configurations. The proposed modelling approach is based on metallur...This paper presents an experimentally validated weld joint shape and dimensions predictive 3D modeling for low carbon galvanized steel in butt-joint configurations. The proposed modelling approach is based on metallurgical transformations using temperature dependent material properties and the enthalpy method. Conduction and keyhole modes welding are investigated using surface and volumetric heat sources, respectively. Transition between the heat sources is carried out according to the power density and interaction time. Simulations are carried out using 3D finite element model on commercial software. The simulation results of the weld shape and dimensions are validated using a structured experimental investigation based on Taguchi method. Experimental validation conducted on a 3 kW Nd: YAG laser source reveals that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects. The results show great concordance between predicted and measured values for the weld joint shape and dimensions.展开更多
This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuato...This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuator configurations.The 3D governing equations for the magneto-electro-elastic static behavior of plates are explicitly show that are made by the three 3D equilibrium equations,the 3D divergence equation for magnetic induction,and the 3D divergence equation for the electric displacement.The proposed solution involves the exponential matrix in the thickness direction and primary variables’harmonic forms in the in-plane ones.A closed-form solution is performed considering simply-supported boundary conditions.Interlaminar continuity conditions are imposed for displacements,magnetic potential,electric potential,transverse shear/normal stresses,transverse normal magnetic induction and transverse normal electric displacement.Therefore,a layerwise approach is adopted.The results section is composed of an assessment part,where the present model is compared to past 3D electro-elastic or magneto-elastic formulations and a new benchmark part.Benchmarks consider sensor and actuator plate configurations for the fully coupled magneto-electro-elastic cases for different thickness ratios.Tabular and graphical results are presented for displacements,stresses,magnetic potential,electric potential,transverse normal magnetic induction and transverse normal electric displacement.For each presented benchmark,magneto-electro-elastic coupling and thickness and material layer effects are discussed in depth.展开更多
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.展开更多
Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled me...Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.展开更多
The 3D visualization model of slop with structural plane can displayed the characters of structural plane in slop directly, and illustrated the spatial combination. It is a modem and critical question in the field of ...The 3D visualization model of slop with structural plane can displayed the characters of structural plane in slop directly, and illustrated the spatial combination. It is a modem and critical question in the field of geotechnical engineering. Based on the peculiarity of the reconnaissance and the research of the visualization by formers, systemized the method fit for building 3D visualization model of slop with structural plane. Write the special program with Visual C^-＋ computer language and illustrated it by OpenGL, the program can displayed and captured the random section plane. The program has a satisfied result by proving with the real projects.展开更多
文摘This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creep/relaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.
基金supported by National Natural Science Foundation of China(Grant No. 51175287)National Science and Technology Major Project(Grant No. 2011ZX02403)
文摘Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the ＂Shape Distribution （D2）＂ algorithm, and a new algorithm named ＂Quick D2＂ is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, ＂3D Searcher＂, a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.
文摘3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.
基金Sponsored by the National Natural Science Foundation of China (Grant No. 51075083)
文摘In order to improve the accuracy and efficiency of 3D model retrieval,the method based on affinity propagation clustering algorithm is proposed. Firstly,projection ray-based method is proposed to improve the feature extraction efficiency of 3D models. Based on the relationship between model and its projection,the intersection in 3D space is transformed into intersection in 2D space,which reduces the number of intersection and improves the efficiency of the extraction algorithm. In feature extraction,multi-layer spheres method is analyzed. The two-layer spheres method makes the feature vector more accurate and improves retrieval precision. Secondly,Semi-supervised Affinity Propagation ( S-AP) clustering is utilized because it can be applied to different cluster structures. The S-AP algorithm is adopted to find the center models and then the center model collection is built. During retrieval process,the collection is utilized to classify the query model into corresponding model base and then the most similar model is retrieved in the model base. Finally,75 sample models from Princeton library are selected to do the experiment and then 36 models are used for retrieval test. The results validate that the proposed method outperforms the original method and the retrieval precision and recall ratios are improved effectively.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
文摘In this study, the petrophysical parameters such as density, sonic, neutron, and porosity were investigated and presented in the 3D models. The 3D models were built using geostatistical method that is used to estimate studied parameters in the entire reservoir. For this purpose, the variogram of each parameter was determined to specify spatial correlation of data. Resulted variograms were non-monotonic. That shows anisotropy of structure. The lithology and porosity parameters are the main causes of this anisotropy. The 3D models also show that petrophysical data has higher variation in north part of reservoir than south part. In addition to, the west limb of reservoir shows higher porosity than east limb. The variation of sonic and neutron data are similar whereas the density data has opposed variation.
基金This paper is supported by the National Natural Science Foundation ofChina (No .40344002) .
文摘A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.
基金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.
基金This research is sponsored by the National Natural Science Foundation of China （No. 40374024）.
文摘The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.
基金Projects(51804113,51434006,51874130)supported by the National Natural Science Foundation of ChinaProject(E51768)supported by the Doctoral Initiation Foundation of Hunan University of Science and Technology,China+1 种基金Project(E61610)supported by the Postdoctoral Research Foundation of Hunan University of Science and Technology,ChinaProject(E21734)supported by the Open Foundation of Work Safety Key Lab on Prevention and Control of Gas and Roof Disasters for Southern Coal Mines,China
文摘A new horn failure mechanism was constructed for tunnel faces in the soft rock mass by means of the logarithmic spiral curve. The seismic action was incorporated into the horn failure mechanism using the pseudo-static method. Considering the randomness of rock mass parameters and loads, a three-dimensional (3D) stochastic collapse model was established. Reliability analysis of seismic stability of tunnel faces was presented via the kinematical approach and the response surface method. The results show that, the reliability of tunnel faces is significantly affected by the supporting pressure, geological strength index, uniaxial compressive strength, rock bulk density and seismic forces. It is worth noting that, if the effect of seismic force was not considered, the stability of tunnel faces would be obviously overestimated. However, the correlation between horizontal and vertical seismic forces can be ignored under the condition of low calculation accuracy.
文摘A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘Experimental X-ray crystallography, NMR (Nuclear Magnetic Resonance) spectroscopy, dual polarization interferometry, etc. are indeed very powerful tools to determine the 3-Dimensional structure of a protein (including the membrane protein);theoretical mathematical and physical computational approaches can also allow us to obtain a description of the protein 3D structure at a submicroscopic level for some unstable, noncrystalline and insoluble proteins. X-ray crystallography finds the X-ray final structure of a protein, which usually need refinements using theoretical protocols in order to produce a better structure. This means theoretical methods are also important in determinations of protein structures. Optimization is always needed in the computer-aided drug design, structure-based drug design, molecular dynamics, and quantum and molecular mechanics. This paper introduces some optimization algorithms used in these research fields and presents a new theoretical computational method—an improved LBFGS Quasi-Newtonian mathematical optimization method—to produce 3D structures of prion AGAAAAGA amyloid fibrils (which are unstable, noncrystalline and insoluble), from the potential energy minimization point of view. Because the NMR or X-ray structure of the hydrophobic region AGAAAAGA of prion proteins has not yet been determined, the model constructed by this paper can be used as a reference for experimental studies on this region, and may be useful in furthering the goals of medicinal chemistry in this field.
基金Supported by NationalNatural Science FoundationCouncil of the People’s Republic of China (20490224)
文摘The design of the cooling system of injection molds directly affects both productivity and the quality of the final part. Using the cooling process CAE system to instruct the mold design, the efficiency and quality of design can be improved greatly. At the same time, it is helpful to confirm the cooling system structure and optimize the process conditions. In this paper, the 3D surface model of mold cavity is used to replace the middle-plane model in the simulation by Boundary Element Method, which break the bottleneck of the application of the injection molding simulation softwares base on the middle-plane model. With the improvements of this paper, a practical and commercial simulation software of injection molding cooling process named as HsCAE3D6.0 is developed.
文摘This paper presents an experimentally validated weld joint shape and dimensions predictive 3D modeling for low carbon galvanized steel in butt-joint configurations. The proposed modelling approach is based on metallurgical transformations using temperature dependent material properties and the enthalpy method. Conduction and keyhole modes welding are investigated using surface and volumetric heat sources, respectively. Transition between the heat sources is carried out according to the power density and interaction time. Simulations are carried out using 3D finite element model on commercial software. The simulation results of the weld shape and dimensions are validated using a structured experimental investigation based on Taguchi method. Experimental validation conducted on a 3 kW Nd: YAG laser source reveals that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects. The results show great concordance between predicted and measured values for the weld joint shape and dimensions.
文摘This study proposes a three-dimensional(3D)coupled magneto-electro-elastic problem for the static analysis of multilayered plates embedding piezomagnetic and piezoelectric layers by considering both sensor and actuator configurations.The 3D governing equations for the magneto-electro-elastic static behavior of plates are explicitly show that are made by the three 3D equilibrium equations,the 3D divergence equation for magnetic induction,and the 3D divergence equation for the electric displacement.The proposed solution involves the exponential matrix in the thickness direction and primary variables’harmonic forms in the in-plane ones.A closed-form solution is performed considering simply-supported boundary conditions.Interlaminar continuity conditions are imposed for displacements,magnetic potential,electric potential,transverse shear/normal stresses,transverse normal magnetic induction and transverse normal electric displacement.Therefore,a layerwise approach is adopted.The results section is composed of an assessment part,where the present model is compared to past 3D electro-elastic or magneto-elastic formulations and a new benchmark part.Benchmarks consider sensor and actuator plate configurations for the fully coupled magneto-electro-elastic cases for different thickness ratios.Tabular and graphical results are presented for displacements,stresses,magnetic potential,electric potential,transverse normal magnetic induction and transverse normal electric displacement.For each presented benchmark,magneto-electro-elastic coupling and thickness and material layer effects are discussed in depth.
基金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 Major Project of National Science and Technology of China(No.2016ZX05026-002-003)National Natural Science Foundation of China(No.41374108)
文摘Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.
文摘The 3D visualization model of slop with structural plane can displayed the characters of structural plane in slop directly, and illustrated the spatial combination. It is a modem and critical question in the field of geotechnical engineering. Based on the peculiarity of the reconnaissance and the research of the visualization by formers, systemized the method fit for building 3D visualization model of slop with structural plane. Write the special program with Visual C^-＋ computer language and illustrated it by OpenGL, the program can displayed and captured the random section plane. The program has a satisfied result by proving with the real projects.
