In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of ...In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of volume visualization is explored, based on the characteristics of explosion and shock. Based on this, a visualization system for 3-D explosion--ViSC3D is designed. Approaches for the visualization of 3-D field of explosion are presented. The algorithm and the functions of ViSC3D are also presented. ViSC3D is thus a useful tool to observe and analyze either the full picture or the details of a 3-D field of explosion, that are difficult to observe and analyze directly. With ViSC3D, the field of explosion between the hill slopes is visualized. The cutaway views and 2-D slices are also given. The full picture and partial details of 3-D field of explosion can be observed clearly. Furthermore, ViSC3D can be used to visualize other similar 3-D data fields.展开更多
Recent studies have demonstrated the importance of LUCC change with climate and ecosystem simulation, but the result could only be determined precisely if a high-resolution underlying land cover map is used. While the...Recent studies have demonstrated the importance of LUCC change with climate and ecosystem simulation, but the result could only be determined precisely if a high-resolution underlying land cover map is used. While the efforts based satellites have provided a good baseline for present land cover, what the next advancement in the research about LUCC change required is the development of reconstruction of historical LUCC change especially spatially-explicit historical dataset. Being different from other similar studies, this study is based on the analysis of historical land use patterns in the traditional cultivated region of China. Taking no account of the less important factors, altitude, slope and population patterns are selected as the major drivers of reclamation in ancient China, and used to design the HCGM (Historical Cropland Gridding Model, at a 60 km×60 km resolution), which is an empirical model for allocating the historical cropland inventory data spatially to grid cells in each political unit. Then we use this model to reconstruct cropland distribution of the study area in 1820, and verify the result by prefectural cropland data of 1820, which is from the historical documents. The statistical analyzing result shows that the model can simulate the patterns of the cropland distribution in the historical period in the traditional cultivated region efficiently.展开更多
Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement...Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp://lsec. cc. ac. cn/phg/), a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.展开更多
Glacier area changes in the Qangtang Plateau are analyzed during 1970-2000 using air photos,relevant photogrammetric maps and satellite images based on the multi-temporal grid method.The results indicate that the melt...Glacier area changes in the Qangtang Plateau are analyzed during 1970-2000 using air photos,relevant photogrammetric maps and satellite images based on the multi-temporal grid method.The results indicate that the melting of glaciers accelerated,only a few of glaciers in an advancing state during 1970-2000 in the whole Qangtang Plateau.However,the glaciers seemed still more stable in the study area than in most areas of western China.We estimate that glacier retreat was likely due to air temperature warming during 1970-2000 in the Qangtang Plateau.Furthermore,the functional model of glacier system is applied to study climate sensitivity of glacier area changes,which indicates that glacier lifespan mainly depends on the heating rate,secondly the precipitation,and precipitation increasing can slow down glacier retreat and make glacier lifespan prolonged.展开更多
The elasticity, viscosity, and the relationships derived from rheology weakness properties are taken into account in mechanics. Comparing with the corresponding relationships derived from damage mechanics, we find the...The elasticity, viscosity, and the relationships derived from rheology weakness properties are taken into account in mechanics. Comparing with the corresponding relationships derived from damage mechanics, we find the weakness factor has the same significance as the damage factor. We simulate the wave field using a staggered-grid pseudospectral method to show the influence of the weakness factor qualitatively. Applying the analytical solution of plane waves, we give the velocity and attenuation coefficient of three body waves, which are affected by the wave frequency and the weakness factor of saturated discrete media. Our results show that velocity decreases with increasing weakness factor, the attenuation coefficient increases with an increase in the weakness factor, and that the influence of weakness depends on the mode of the body waves.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where fo...An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where for the purpose of stability and efficiency a Lagrangian smoothing strain at nodal point is constructed and thereafter the internal energy is evaluated nodally. This formulation ensures the linear exactness, efficiency and spatial stability in a unified manner and it makes the conventional Galerkin meshfree method affordable for three dimensional simulation. The three dimensional implementation of stabilized conforming nodal integration is discussed in details. To model the failure evolution in soil medium a coupled elasto-plastic damage model is used and an objective stress integration algorithm in combination of elasto-damage predictor and plastic corrector method is employed for stress update. Two typical numerical examples are shown to demonstrate the effectiveness of the present method for modeling large deformation soil failure.展开更多
We present a novel algorithm for adaptive triangular mesh coarsening. The algorithm has two stages. First, the input triangular mesh is refined by iteratively applying the adaptive subdivision operator that performs a...We present a novel algorithm for adaptive triangular mesh coarsening. The algorithm has two stages. First, the input triangular mesh is refined by iteratively applying the adaptive subdivision operator that performs a so-called red-green split. Second, the refined mesh is simplified by a clustering algorithm based on centroidal Voronoi tessellations (CVTs). The accuracy and good quality of the output triangular mesh are achieved by combining adaptive subdivision and the CVTs technique. Test results showed the mesh coarsening scheme to be robust and effective. Examples are shown that validate the method.展开更多
Mesh parameterization is one of the fundamental operations in computer graphics(CG) and computeraided design(CAD). In this paper, we propose a novel local/global parameterization approach, ARAP++, for singleand multi-...Mesh parameterization is one of the fundamental operations in computer graphics(CG) and computeraided design(CAD). In this paper, we propose a novel local/global parameterization approach, ARAP++, for singleand multi-boundary triangular meshes. It is an extension of the as-rigid-as-possible(ARAP) approach, which stitches together 1-ring patches instead of individual triangles. To optimize the spring energy, we introduce a linear iterative scheme which employs convex combination weights and a fitting Jacobian matrix corresponding to a prescribed family of transformations. Our algorithm is simple, efficient, and robust. The geometric properties(angle and area)of the original model can also be preserved by appropriately prescribing the singular values of the fitting matrix. To reduce the area and stretch distortions for high-curvature models, a stretch operator is introduced. Numerical results demonstrate that ARAP++ outperforms several state-of-the-art methods in terms of controlling the distortions of angle, area, and stretch. Furthermore, it achieves a better visualization performance for several applications, such as texture mapping and surface remeshing.展开更多
Surface remeshing is widely required in modeling, animation, simulation, and many other computer graphics applications. Improving the elements' quality is a challenging task in surface remeshing. Existing methods ...Surface remeshing is widely required in modeling, animation, simulation, and many other computer graphics applications. Improving the elements' quality is a challenging task in surface remeshing. Existing methods often fail to efficiently remove poor-quality elements especially in regions with sharp features. In this paper, we propose and use a robust segmentation method followed by remeshing the segmented mesh. Mesh segmentation is initiated using an existing Live-wire interaction approach and is further refined using local mesh operations. The refined segmented mesh is finally sent to the remeshing pipeline, in which each mesh segment is remeshed independently. An experimental study compares our mesh segmentation method as well as remeshing results with representative existing methods. We demonstrate that the proposed segmentation method is robust and suitable for remeshing.展开更多
Understanding the mechanisms of hard–soft material interaction under impact loading is important not only in the defense industry but also in daily life.However,traditional mesh-based spatial discretization methods t...Understanding the mechanisms of hard–soft material interaction under impact loading is important not only in the defense industry but also in daily life.However,traditional mesh-based spatial discretization methods that are time consuming owing to the need for frequent re-meshing,such as the finite element method and finite difference method,can hardly handle large deformation involving failure evolution in a multi-phase interaction environment.The objective of this research is to develop a quasi-meshless particle method based on the material point method for the model-based simulation of the hard–soft material interaction response.To demonstrate the proposed procedure,scenarios of a hard–soft material impact test are considered,where a force is applied to layers of materials and a hard bar with an initial velocity impacts a target with layers of different materials.The stress wave propagation and resulting failure evolution are simulated and compared with available data.Future research tasks are then discussed on the basis of the preliminary results.展开更多
In fracture simulation,how to model the pre-existing cracks and simulate their propagation without remeshing is an important topic.The newly developed triangular element partition method(TEPM)provides an efficient app...In fracture simulation,how to model the pre-existing cracks and simulate their propagation without remeshing is an important topic.The newly developed triangular element partition method(TEPM)provides an efficient approach to this problem.It firstly meshes the cracked body regardless of the geometry integrity of the interesting object with triangular elements.After the meshing procedure is completed,some elements are intersected by cracks.For the element intersected by a crack,the TEPM takes the element partition technique to incorporate the discontinuity into the numerical model without any interpolation enrichment.By this approach,the TEPM can simulate fracture without mesh modification.In the TEPM,all the cracked elements are treated as the usual partitioned elements in which the crack runs through.The virtual node pairs(the intersection points of crack faces and elements)at the opposite faces of the crack move independently.Their displacements are respectively determined by their neighbor real nodes(nodes formatted in the original mesh scheme)at the same side of the crack.However,among these cracked elements,the element containing a crack tip,referred to as the crack tip element thereafter,behaves differently from those cut through by the crack.Its influence on the singular field at the vicinity of the fracture tip becomes increasingly significant with the element size increasing.In the crack tip element,the virtual node pair at the crack tip move consistently before fracture occurs while the virtual node pair separate and each virtual node moves independently after the fracture propagates.Accordingly,the crack tip element is automatically transformed into the usual partitioned element.In the present paper,the crack tip element is introduced into the TEPM to account for the effect of the crack tip.Validation examples indicate that the present method is almost free from the element size effect.It can reach the same precision as the conventional finite element method under the same meshing scheme.But the TEPM is much more efficient and convenient than the conventional finite element method because the TEPM avoids the troubles that the conventional finite element method suffers,e.g.,the meshing problem of cracked body,modification of mesh scheme,etc.Though the extended finite element method can also avoid these troubles,it introduces extra degrees of freedom due to node interpolation enrichment.Due to the simplicity of the present TEPM,it is believed that its perspective should be highly inspiring.展开更多
文摘In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of volume visualization is explored, based on the characteristics of explosion and shock. Based on this, a visualization system for 3-D explosion--ViSC3D is designed. Approaches for the visualization of 3-D field of explosion are presented. The algorithm and the functions of ViSC3D are also presented. ViSC3D is thus a useful tool to observe and analyze either the full picture or the details of a 3-D field of explosion, that are difficult to observe and analyze directly. With ViSC3D, the field of explosion between the hill slopes is visualized. The cutaway views and 2-D slices are also given. The full picture and partial details of 3-D field of explosion can be observed clearly. Furthermore, ViSC3D can be used to visualize other similar 3-D data fields.
基金Natiional Natural Science Foundation of China,No.40471007Innovation Knowledge Project of CAS,No.KZCX2-YW-315
文摘Recent studies have demonstrated the importance of LUCC change with climate and ecosystem simulation, but the result could only be determined precisely if a high-resolution underlying land cover map is used. While the efforts based satellites have provided a good baseline for present land cover, what the next advancement in the research about LUCC change required is the development of reconstruction of historical LUCC change especially spatially-explicit historical dataset. Being different from other similar studies, this study is based on the analysis of historical land use patterns in the traditional cultivated region of China. Taking no account of the less important factors, altitude, slope and population patterns are selected as the major drivers of reclamation in ancient China, and used to design the HCGM (Historical Cropland Gridding Model, at a 60 km×60 km resolution), which is an empirical model for allocating the historical cropland inventory data spatially to grid cells in each political unit. Then we use this model to reconstruct cropland distribution of the study area in 1820, and verify the result by prefectural cropland data of 1820, which is from the historical documents. The statistical analyzing result shows that the model can simulate the patterns of the cropland distribution in the historical period in the traditional cultivated region efficiently.
基金supported by the 973 Program of China 2005CB321702China NSF 10531080.
文摘Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp://lsec. cc. ac. cn/phg/), a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.
基金supported by the National Natural Science Foundation of China (Nos.40871043,40801025)the Project of National Scientific Basic Special Fund on the Ministry of Science and Technology of China (No.2006FY110200)the Key Construction Disciplines of Hunan Province (No.40652001)
文摘Glacier area changes in the Qangtang Plateau are analyzed during 1970-2000 using air photos,relevant photogrammetric maps and satellite images based on the multi-temporal grid method.The results indicate that the melting of glaciers accelerated,only a few of glaciers in an advancing state during 1970-2000 in the whole Qangtang Plateau.However,the glaciers seemed still more stable in the study area than in most areas of western China.We estimate that glacier retreat was likely due to air temperature warming during 1970-2000 in the Qangtang Plateau.Furthermore,the functional model of glacier system is applied to study climate sensitivity of glacier area changes,which indicates that glacier lifespan mainly depends on the heating rate,secondly the precipitation,and precipitation increasing can slow down glacier retreat and make glacier lifespan prolonged.
基金0ur work is supported by the 0pen Fund of the CNPC Key Lab of Geophysical Exploration (GPKL0202), the 0pen Fund of the State Key Laboratory of 0il and Gas Reservoir Geology and Exploitation (PLC200304), and the Natural Science Foundation of Hubei Province (2002AB018).
文摘The elasticity, viscosity, and the relationships derived from rheology weakness properties are taken into account in mechanics. Comparing with the corresponding relationships derived from damage mechanics, we find the weakness factor has the same significance as the damage factor. We simulate the wave field using a staggered-grid pseudospectral method to show the influence of the weakness factor qualitatively. Applying the analytical solution of plane waves, we give the velocity and attenuation coefficient of three body waves, which are affected by the wave frequency and the weakness factor of saturated discrete media. Our results show that velocity decreases with increasing weakness factor, the attenuation coefficient increases with an increase in the weakness factor, and that the influence of weakness depends on the mode of the body waves.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10972188, 10602049)the Program for New Century Excellent Talents in University from China Education Ministry (Grant No. NCET-09-0678)the Fundamental Research Funds for the Central Universities of China (Grant No. 2010121073)
文摘An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where for the purpose of stability and efficiency a Lagrangian smoothing strain at nodal point is constructed and thereafter the internal energy is evaluated nodally. This formulation ensures the linear exactness, efficiency and spatial stability in a unified manner and it makes the conventional Galerkin meshfree method affordable for three dimensional simulation. The three dimensional implementation of stabilized conforming nodal integration is discussed in details. To model the failure evolution in soil medium a coupled elasto-plastic damage model is used and an objective stress integration algorithm in combination of elasto-damage predictor and plastic corrector method is employed for stress update. Two typical numerical examples are shown to demonstrate the effectiveness of the present method for modeling large deformation soil failure.
基金supported by the National Natural Science Foundation of China (No. 60773179)the National Basic Research Program (973) of China (No. 2004CB318000)
文摘We present a novel algorithm for adaptive triangular mesh coarsening. The algorithm has two stages. First, the input triangular mesh is refined by iteratively applying the adaptive subdivision operator that performs a so-called red-green split. Second, the refined mesh is simplified by a clustering algorithm based on centroidal Voronoi tessellations (CVTs). The accuracy and good quality of the output triangular mesh are achieved by combining adaptive subdivision and the CVTs technique. Test results showed the mesh coarsening scheme to be robust and effective. Examples are shown that validate the method.
基金supported by the National Natural Science Foundation of China(Nos.61432003,61572105,11171052,and 61328206)
文摘Mesh parameterization is one of the fundamental operations in computer graphics(CG) and computeraided design(CAD). In this paper, we propose a novel local/global parameterization approach, ARAP++, for singleand multi-boundary triangular meshes. It is an extension of the as-rigid-as-possible(ARAP) approach, which stitches together 1-ring patches instead of individual triangles. To optimize the spring energy, we introduce a linear iterative scheme which employs convex combination weights and a fitting Jacobian matrix corresponding to a prescribed family of transformations. Our algorithm is simple, efficient, and robust. The geometric properties(angle and area)of the original model can also be preserved by appropriately prescribing the singular values of the fitting matrix. To reduce the area and stretch distortions for high-curvature models, a stretch operator is introduced. Numerical results demonstrate that ARAP++ outperforms several state-of-the-art methods in terms of controlling the distortions of angle, area, and stretch. Furthermore, it achieves a better visualization performance for several applications, such as texture mapping and surface remeshing.
基金the National Natural Science Foundation of China(Nos.61772523,61372168,61620106003,and 61331018)supported by a Chinese Government Scholarship
文摘Surface remeshing is widely required in modeling, animation, simulation, and many other computer graphics applications. Improving the elements' quality is a challenging task in surface remeshing. Existing methods often fail to efficiently remove poor-quality elements especially in regions with sharp features. In this paper, we propose and use a robust segmentation method followed by remeshing the segmented mesh. Mesh segmentation is initiated using an existing Live-wire interaction approach and is further refined using local mesh operations. The refined segmented mesh is finally sent to the remeshing pipeline, in which each mesh segment is remeshed independently. An experimental study compares our mesh segmentation method as well as remeshing results with representative existing methods. We demonstrate that the proposed segmentation method is robust and suitable for remeshing.
基金The National Natural Science Foundation of China(Grant Nos.51476150,11102185 and 11232003)U.S.Defense Threat Reduction Agency(Grant No.HDTRA1-10-1-0022)+1 种基金International Joint Research Program of Shanxi Province,China(Grant No.2014081028)Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi
文摘Understanding the mechanisms of hard–soft material interaction under impact loading is important not only in the defense industry but also in daily life.However,traditional mesh-based spatial discretization methods that are time consuming owing to the need for frequent re-meshing,such as the finite element method and finite difference method,can hardly handle large deformation involving failure evolution in a multi-phase interaction environment.The objective of this research is to develop a quasi-meshless particle method based on the material point method for the model-based simulation of the hard–soft material interaction response.To demonstrate the proposed procedure,scenarios of a hard–soft material impact test are considered,where a force is applied to layers of materials and a hard bar with an initial velocity impacts a target with layers of different materials.The stress wave propagation and resulting failure evolution are simulated and compared with available data.Future research tasks are then discussed on the basis of the preliminary results.
基金supported by the National Natural Science Foundation of China (Grant No. 11172172)the National Basic Research Program of China ("973" Project) (Grant No. 2011CB013505)
文摘In fracture simulation,how to model the pre-existing cracks and simulate their propagation without remeshing is an important topic.The newly developed triangular element partition method(TEPM)provides an efficient approach to this problem.It firstly meshes the cracked body regardless of the geometry integrity of the interesting object with triangular elements.After the meshing procedure is completed,some elements are intersected by cracks.For the element intersected by a crack,the TEPM takes the element partition technique to incorporate the discontinuity into the numerical model without any interpolation enrichment.By this approach,the TEPM can simulate fracture without mesh modification.In the TEPM,all the cracked elements are treated as the usual partitioned elements in which the crack runs through.The virtual node pairs(the intersection points of crack faces and elements)at the opposite faces of the crack move independently.Their displacements are respectively determined by their neighbor real nodes(nodes formatted in the original mesh scheme)at the same side of the crack.However,among these cracked elements,the element containing a crack tip,referred to as the crack tip element thereafter,behaves differently from those cut through by the crack.Its influence on the singular field at the vicinity of the fracture tip becomes increasingly significant with the element size increasing.In the crack tip element,the virtual node pair at the crack tip move consistently before fracture occurs while the virtual node pair separate and each virtual node moves independently after the fracture propagates.Accordingly,the crack tip element is automatically transformed into the usual partitioned element.In the present paper,the crack tip element is introduced into the TEPM to account for the effect of the crack tip.Validation examples indicate that the present method is almost free from the element size effect.It can reach the same precision as the conventional finite element method under the same meshing scheme.But the TEPM is much more efficient and convenient than the conventional finite element method because the TEPM avoids the troubles that the conventional finite element method suffers,e.g.,the meshing problem of cracked body,modification of mesh scheme,etc.Though the extended finite element method can also avoid these troubles,it introduces extra degrees of freedom due to node interpolation enrichment.Due to the simplicity of the present TEPM,it is believed that its perspective should be highly inspiring.