We propose a combined shape and topology optimization approach in this research for 3D acoustics by using the isogeometric boundary element method with subdivision surfaces.The existing structural optimization methods...We propose a combined shape and topology optimization approach in this research for 3D acoustics by using the isogeometric boundary element method with subdivision surfaces.The existing structural optimization methods mainly contain shape and topology schemes,with the former changing the surface geometric profile of the structure and the latter changing thematerial distribution topology or hole topology of the structure.In the present acoustic performance optimization,the coordinates of the control points in the subdivision surfaces fine mesh are selected as the shape design parameters of the structure,the artificial density of the sound absorbing material covered on the structure surface is set as the topology design parameter,and the combined topology and shape optimization approach is established through the sound field analysis of the subdivision surfaces boundary element method as a bridge.The topology and shape sensitivities of the approach are calculated using the adjoint variable method,which ensures the efficiency of the optimization.The geometric jaggedness and material distribution discontinuities that appear in the optimization process are overcome to a certain degree by the multiresolution method and solid isotropic material with penalization.Numerical examples are given to validate the effectiveness of the presented optimization approach.展开更多
In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Mill...In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.展开更多
A novel method is produced to evaluate the energy of the Catmull-Clark subdivision surface including extraordinary points in the control mesh. A closed-form analytic formula for thin plate energy of the Catmull-Clark ...A novel method is produced to evaluate the energy of the Catmull-Clark subdivision surface including extraordinary points in the control mesh. A closed-form analytic formula for thin plate energy of the Catmull-Clark subdivision surface of arbitrary topology is derived through translating the Catmull-Clark subdivision surface into bi-cubic B-spline surface pieces. Using this method, both the membrane energy and the thin plate energy can be evaluated without requiring recursive subdivision. Therefore, it is more efficient and more accurate than the existing methods for calculating the energy of the Catmull-Clark subdivision surface with arbitrary topology. The example of surface fairing demonstrates that this method is efficient and successful for evaluating the energy of subdivision surfaces.展开更多
PHT-splines are defined as polynomial splines over hierarchical T-meshes with very efficient local refinement properties.The original PHT-spline basis functions constructed by the truncation mechanism have a decay phe...PHT-splines are defined as polynomial splines over hierarchical T-meshes with very efficient local refinement properties.The original PHT-spline basis functions constructed by the truncation mechanism have a decay phenomenon,resulting in numerical instability.The non-decay basis functions are constructed as the B-splines that are defined on the 2×2 tensor product meshes associated with basis vertices in Kang et al.,but at the cost of losing the partition of unity.In the field of finite element analysis and topology optimization,forming the partition of unity is the default ingredient for constructing basis functions of approximate spaces.In this paper,we will show that the non-decay PHT-spline basis functions proposed by Kang et al.can be appropriately modified to form a partition of unity.Each non-decay basis function is multiplied by a positive weight to form the weighted basis.The weights are solved such that the sum of weighted bases is equal to 1 on the domain.We provide two methods for calculatingweights,based on geometric information of basis functions and the subdivision of PHT-splines.Weights are given in the form of explicit formulas and can be efficiently calculated.We also prove that the weights on the admissible hierarchical T-meshes are positive.展开更多
Although significant progress has been made in precision machining of free-form surfaces recently, inspection of such surfaces remains a difficult problem. In order to solve the problem that no specific standards for ...Although significant progress has been made in precision machining of free-form surfaces recently, inspection of such surfaces remains a difficult problem. In order to solve the problem that no specific standards for the verification of free-form surface profile are available, the profile parameters of free-form surface are proposed by referring to ISO standards regarding form tolerances and considering its complexity and non-rotational symmetry. Non-uniform rational basis spline(NURBS) for describing free-form surface is formulated. Crucial issues in surface inspection and profile error verification are localization between the design coordinate system(DCS) and measurement coordinate system(MCS) for searching the closest points on the design model corresponding to measured points. A quasi particle swarm optimization(QPSO) is proposed to search the transformation parameters to implement localization between DCS and MCS. Surface subdivide method which does the searching in a recursively reduced range of the parameters u and v of the NURBS design model is developed to find the closest points. In order to verify the effectiveness of the proposed methods, the design model is generated by NURBS and the measurement data of simulation example are generated by transforming the design model to arbitrary position and orientation, and the parts are machined based on the design model and are measured on CMM. The profile errors of simulation example and actual parts are calculated by the proposed method. The results verify that the evaluation precision of freeform surface profile error by the proposed method is higher 10%-22% than that by CMM software. The proposed method deals with the hard problem that it has a lower precision in profile error evaluation of free-form surface.展开更多
On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the...On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software "Auto mapping system of submarine topography and geomorphology MBChart". This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China' s next Five-year Plan.展开更多
By using some elementary inequalities, authors in this paper makes further improvement for estimating the heights of Bézier curve and rational Bézier curve. And the termination criterion for subdivision of t...By using some elementary inequalities, authors in this paper makes further improvement for estimating the heights of Bézier curve and rational Bézier curve. And the termination criterion for subdivision of the rational Bézier curve is also improved. The conclusion of the extreme value problem is thus further confirmed.展开更多
In this paper we introduce a new reverse Loop subdivision method. In contrast to current wavelets based Loop subdivision scheme, our method applies the same rules to both regular and extraordinary vertices and reconst...In this paper we introduce a new reverse Loop subdivision method. In contrast to current wavelets based Loop subdivision scheme, our method applies the same rules to both regular and extraordinary vertices and reconstructs the sharp features easily. Furthermore, our method runs faster because it does not need analysis and synthesis procedural. Our main goal is the design of a reverse subdivision method that can reconstruct the coarser mesh from a finer subdivision surface with sharp features for multiresolution representation, The proposed method only needs a little memory storage and brings little error, and it is easy to implement.展开更多
The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built b...The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.展开更多
Arbitrary topological curve network has no restriction in topology structure,so it has more powerful representing ability in defining complex surfaces.A complex surface modeling system is presented based on arbitrary ...Arbitrary topological curve network has no restriction in topology structure,so it has more powerful representing ability in defining complex surfaces.A complex surface modeling system is presented based on arbitrary topological curve network and the improved combined subdivision method,its functions including creating and editing curve network,and generating and modifying curve network's interpolated surface.This modeling system can be used to the process of products'concept design,and its applications is also significant to the development of subdivision method.展开更多
In this article, the objective is to introduce an algorithm to produce the quaternary m-point (for any integer m>1) approximating subdivision schemes, which have smaller support and higher smoothness, comparing to ...In this article, the objective is to introduce an algorithm to produce the quaternary m-point (for any integer m>1) approximating subdivision schemes, which have smaller support and higher smoothness, comparing to binary and ternary schemes. The proposed algorithm has been derived from uniform B-spline basis function using the Cox-de Boor recursion formula. In order to determine the convergence and smoothness of the proposed schemes, the Laurent polynomial method has been used.展开更多
Design of watertight subdivision inherently involves its optimization with the objective to increase the index "A" above its minimum required value. In view of a big popularity of probabilistic search methods such a...Design of watertight subdivision inherently involves its optimization with the objective to increase the index "A" above its minimum required value. In view of a big popularity of probabilistic search methods such as genetic algorithms, this task is intrinsically time consuming. Thus, even when an optimal subdivision layout (i.e. topology) is determined, it can be found that the optimal bulkhead positions can be a great challenge time-wise, often forcing designers to satisfy with suboptimal solutions. The fundamental reason why this happens is that the nature of the optimized function (e.g., index "A" as a function of bulkhead positions) is unknown and hence it has no effect upon the choice of optimization strategy, which therefore reflects subjective but not factual preferences. In this paper we study the nature of functional dependency between the subdivision index and bulkhead positions, as a simplest case, and indicate pertinent optimization strategies that consequently reduce the optimization time. In our study we use a cruise ship model to demonstrate the application results of our findings.展开更多
In this paper, we propose and analyze a subdivision scheme which unifies 3-point approximating subdivision schemes of any arity in its compact form and has less support, computational cost and error bounds.? The usefu...In this paper, we propose and analyze a subdivision scheme which unifies 3-point approximating subdivision schemes of any arity in its compact form and has less support, computational cost and error bounds.? The usefulness of the scheme is illustrated by considering different examples along with its comparison with the established subdivision schemes. Moreover, B-splines of degree 4and well known 3-point schemes [1, 2, 3, 4, 6, 11, 12, 14, 15] are special cases of our proposed scheme.展开更多
This paper proposes a novel strategy based on fragment meshes of Shiue et al. for GPU rendering of compos- ite subdivision surfaces. Two enumeration systems are established to label the primitives of each fragment mes...This paper proposes a novel strategy based on fragment meshes of Shiue et al. for GPU rendering of compos- ite subdivision surfaces. Two enumeration systems are established to label the primitives of each fragment mesh. A sector-layer-index enumeration system is responsible for retrieving proximities for subdivision masks while a sec- tor-index enumeration system designates a 2D texture buffer in GPU. Recurring to the free conversion between them, our approach may get rid of lookup tables that are designed to record subdivision masks. In addition, relatively small composite subdivision masks make it easy to develop automatically retrieving techniques. Finally, as center vertices are often irregular, their computation is related to an average with alterable number of items. Considering that variable loop is not efficient in GPU, we evaluate the center vertex of each fragment mesh using the linear combination of its level 0, level 1 and limit positions instead of averaging schemes. Experiments demonstrate that our approach generally outper- forms that of Shiue et al. in FPS by a long way.展开更多
In this article,we adopt the C-type spline of degree 2 to model and blend basic shapes including conics and circle arcs.The C-type spline belongs to theωB-spline category of splines that are capable of blending polyn...In this article,we adopt the C-type spline of degree 2 to model and blend basic shapes including conics and circle arcs.The C-type spline belongs to theωB-spline category of splines that are capable of blending polynomial,trigonometric and hyperbolic functions.Commonly used basic shapes can be exactly represented by these types of splines.We derive explicit formulas for the convenience of modeling the basic curves.The entire blending curve is C^1-continuous.In comparison with the existing best blending method by rational G^2 splines,which are rational splines of degree 3,the proposed method allows simpler representation and blending of the basic curves,and it can represent numerous basic shapes including the hyperbolic types.We also design a subdivision method to generate blending curves;this method is precise for the basic curves and approximate for the blending sections.The subdivision process is efficient for modeling and rendering.It has also proven to be C^1-continuous by the asymptotically equivalent theory and the continuity of stationary subdivision method.In addition,we extend the proposed methods to cases involving the modeling and blending of basic surfaces.We provide many examples that illustrate the merits of our methods.展开更多
Based on the butterfly subdivision scheme and the modified butterfly subdivision scheme, an improved butterfly subdivision scheme is proposed. The scheme uses a small stencil of six points to calculate new inserting v...Based on the butterfly subdivision scheme and the modified butterfly subdivision scheme, an improved butterfly subdivision scheme is proposed. The scheme uses a small stencil of six points to calculate new inserting vertex, 2n new vertices are inserted in the 2n triangle faces in each recursion, and the n old vertices are kept, special treatment is given to the boundary, achieving higher smoothness while using small stencils is realized. With the proposed scheme, the number of triangle faces increases only by a factor of 3 in each refinement step. Compared with the butterfly subdivision scheme and the modified butterfly subdivision scheme, the size of triangle faces changes more gradually, which allows one to have greater control over the resolution of a refined mesh.展开更多
Quadrics are of basic importance in Computer Graphics and Computer Aided Design. In this paper,we design a subdivision scheme based on the method suggested by G. Morin and J. Warren to generate conics and quadrics con...Quadrics are of basic importance in Computer Graphics and Computer Aided Design. In this paper,we design a subdivision scheme based on the method suggested by G. Morin and J. Warren to generate conics and quadrics conveniently. Given the control polygon(poly-hedron),the corresponding ellipse (ellipsoid)can be generated. The hyperbolas and hyperboloids are generated based on the generation of ellipses and ellipsoids by a simple transformation. The method in this paper is much simpler and easier to apply than those given by Eugenia Montiel et al.展开更多
In order to study the constitutive behavior of concrete in mesoscopic level, a new method is proposed in this paper. This method uses random polygon particles to simulate full grading broken aggregates of concrete. Ba...In order to study the constitutive behavior of concrete in mesoscopic level, a new method is proposed in this paper. This method uses random polygon particles to simulate full grading broken aggregates of concrete. Based on computational geometry, we carry out the automatic generation of the triangle finite element mesh for the model of random polygon particles of concrete. The finite element mesh generated in this paper is also applicable to many other numerical methods.展开更多
基金supported by the National Natural Science Foundation of China (NSFC)under Grant Nos.12172350,11772322 and 11702238。
文摘We propose a combined shape and topology optimization approach in this research for 3D acoustics by using the isogeometric boundary element method with subdivision surfaces.The existing structural optimization methods mainly contain shape and topology schemes,with the former changing the surface geometric profile of the structure and the latter changing thematerial distribution topology or hole topology of the structure.In the present acoustic performance optimization,the coordinates of the control points in the subdivision surfaces fine mesh are selected as the shape design parameters of the structure,the artificial density of the sound absorbing material covered on the structure surface is set as the topology design parameter,and the combined topology and shape optimization approach is established through the sound field analysis of the subdivision surfaces boundary element method as a bridge.The topology and shape sensitivities of the approach are calculated using the adjoint variable method,which ensures the efficiency of the optimization.The geometric jaggedness and material distribution discontinuities that appear in the optimization process are overcome to a certain degree by the multiresolution method and solid isotropic material with penalization.Numerical examples are given to validate the effectiveness of the presented optimization approach.
基金sponsored by the Graduate Student Research and Innovation Fund of Xinyang Normal University under No.2024KYJJ012.
文摘In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.
文摘A novel method is produced to evaluate the energy of the Catmull-Clark subdivision surface including extraordinary points in the control mesh. A closed-form analytic formula for thin plate energy of the Catmull-Clark subdivision surface of arbitrary topology is derived through translating the Catmull-Clark subdivision surface into bi-cubic B-spline surface pieces. Using this method, both the membrane energy and the thin plate energy can be evaluated without requiring recursive subdivision. Therefore, it is more efficient and more accurate than the existing methods for calculating the energy of the Catmull-Clark subdivision surface with arbitrary topology. The example of surface fairing demonstrates that this method is efficient and successful for evaluating the energy of subdivision surfaces.
基金The work was supported by the NSF of China(No.11801393)the Natural Science Foundation of Jiangsu Province,China(No.BK20180831).
文摘PHT-splines are defined as polynomial splines over hierarchical T-meshes with very efficient local refinement properties.The original PHT-spline basis functions constructed by the truncation mechanism have a decay phenomenon,resulting in numerical instability.The non-decay basis functions are constructed as the B-splines that are defined on the 2×2 tensor product meshes associated with basis vertices in Kang et al.,but at the cost of losing the partition of unity.In the field of finite element analysis and topology optimization,forming the partition of unity is the default ingredient for constructing basis functions of approximate spaces.In this paper,we will show that the non-decay PHT-spline basis functions proposed by Kang et al.can be appropriately modified to form a partition of unity.Each non-decay basis function is multiplied by a positive weight to form the weighted basis.The weights are solved such that the sum of weighted bases is equal to 1 on the domain.We provide two methods for calculatingweights,based on geometric information of basis functions and the subdivision of PHT-splines.Weights are given in the form of explicit formulas and can be efficiently calculated.We also prove that the weights on the admissible hierarchical T-meshes are positive.
基金supported by National Natural Science Foundation of China(Grant No. 51075198)Jiangsu Provincial Natural Science Foundation of China(Grant No. BK2010479)+1 种基金Jiangsu Provincial Project of 333 Talents Engineering of ChinaJiangsu Provincial Project of Six Talented Peak of China
文摘Although significant progress has been made in precision machining of free-form surfaces recently, inspection of such surfaces remains a difficult problem. In order to solve the problem that no specific standards for the verification of free-form surface profile are available, the profile parameters of free-form surface are proposed by referring to ISO standards regarding form tolerances and considering its complexity and non-rotational symmetry. Non-uniform rational basis spline(NURBS) for describing free-form surface is formulated. Crucial issues in surface inspection and profile error verification are localization between the design coordinate system(DCS) and measurement coordinate system(MCS) for searching the closest points on the design model corresponding to measured points. A quasi particle swarm optimization(QPSO) is proposed to search the transformation parameters to implement localization between DCS and MCS. Surface subdivide method which does the searching in a recursively reduced range of the parameters u and v of the NURBS design model is developed to find the closest points. In order to verify the effectiveness of the proposed methods, the design model is generated by NURBS and the measurement data of simulation example are generated by transforming the design model to arbitrary position and orientation, and the parts are machined based on the design model and are measured on CMM. The profile errors of simulation example and actual parts are calculated by the proposed method. The results verify that the evaluation precision of freeform surface profile error by the proposed method is higher 10%-22% than that by CMM software. The proposed method deals with the hard problem that it has a lower precision in profile error evaluation of free-form surface.
基金This study was supported by Projects under contract Nos DY105 China's 0cean-03-01-01 and DY105-03-01-07the National Natural Science Foundation of China under contract No.40506017the Youth Foundation of Marine High-tech Project of China under contract No.2002AA616010.
文摘On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software "Auto mapping system of submarine topography and geomorphology MBChart". This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China' s next Five-year Plan.
文摘By using some elementary inequalities, authors in this paper makes further improvement for estimating the heights of Bézier curve and rational Bézier curve. And the termination criterion for subdivision of the rational Bézier curve is also improved. The conclusion of the extreme value problem is thus further confirmed.
基金Supported by the High Technology Research and Development Progrmnn~ of China (No. 2003AA411310), the National Natural Science Foundation of China (No. 60373070) and Microsoft Research Project 2005-1.
文摘In this paper we introduce a new reverse Loop subdivision method. In contrast to current wavelets based Loop subdivision scheme, our method applies the same rules to both regular and extraordinary vertices and reconstructs the sharp features easily. Furthermore, our method runs faster because it does not need analysis and synthesis procedural. Our main goal is the design of a reverse subdivision method that can reconstruct the coarser mesh from a finer subdivision surface with sharp features for multiresolution representation, The proposed method only needs a little memory storage and brings little error, and it is easy to implement.
文摘The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.
基金Project supported by the Fundamental Research Foundations for the Central Universities (Grant No.2009B30514)
文摘Arbitrary topological curve network has no restriction in topology structure,so it has more powerful representing ability in defining complex surfaces.A complex surface modeling system is presented based on arbitrary topological curve network and the improved combined subdivision method,its functions including creating and editing curve network,and generating and modifying curve network's interpolated surface.This modeling system can be used to the process of products'concept design,and its applications is also significant to the development of subdivision method.
文摘In this article, the objective is to introduce an algorithm to produce the quaternary m-point (for any integer m>1) approximating subdivision schemes, which have smaller support and higher smoothness, comparing to binary and ternary schemes. The proposed algorithm has been derived from uniform B-spline basis function using the Cox-de Boor recursion formula. In order to determine the convergence and smoothness of the proposed schemes, the Laurent polynomial method has been used.
文摘Design of watertight subdivision inherently involves its optimization with the objective to increase the index "A" above its minimum required value. In view of a big popularity of probabilistic search methods such as genetic algorithms, this task is intrinsically time consuming. Thus, even when an optimal subdivision layout (i.e. topology) is determined, it can be found that the optimal bulkhead positions can be a great challenge time-wise, often forcing designers to satisfy with suboptimal solutions. The fundamental reason why this happens is that the nature of the optimized function (e.g., index "A" as a function of bulkhead positions) is unknown and hence it has no effect upon the choice of optimization strategy, which therefore reflects subjective but not factual preferences. In this paper we study the nature of functional dependency between the subdivision index and bulkhead positions, as a simplest case, and indicate pertinent optimization strategies that consequently reduce the optimization time. In our study we use a cruise ship model to demonstrate the application results of our findings.
文摘In this paper, we propose and analyze a subdivision scheme which unifies 3-point approximating subdivision schemes of any arity in its compact form and has less support, computational cost and error bounds.? The usefulness of the scheme is illustrated by considering different examples along with its comparison with the established subdivision schemes. Moreover, B-splines of degree 4and well known 3-point schemes [1, 2, 3, 4, 6, 11, 12, 14, 15] are special cases of our proposed scheme.
基金supported by Natural Science Foundation of China (60973084)Fundamental Research Funds for the Central Universities (2012zz0063)
文摘This paper proposes a novel strategy based on fragment meshes of Shiue et al. for GPU rendering of compos- ite subdivision surfaces. Two enumeration systems are established to label the primitives of each fragment mesh. A sector-layer-index enumeration system is responsible for retrieving proximities for subdivision masks while a sec- tor-index enumeration system designates a 2D texture buffer in GPU. Recurring to the free conversion between them, our approach may get rid of lookup tables that are designed to record subdivision masks. In addition, relatively small composite subdivision masks make it easy to develop automatically retrieving techniques. Finally, as center vertices are often irregular, their computation is related to an average with alterable number of items. Considering that variable loop is not efficient in GPU, we evaluate the center vertex of each fragment mesh using the linear combination of its level 0, level 1 and limit positions instead of averaging schemes. Experiments demonstrate that our approach generally outper- forms that of Shiue et al. in FPS by a long way.
基金This work described in this article was supported by the National Science Foundation of China(61772164,61272032)Provincial Key Platforms and Major Scientific Research Projects in Universities and Colleges of Guangdong(2017KTSCX143)the Natural Science Foundation of Zhejiang Province(LY17F020025).
文摘In this article,we adopt the C-type spline of degree 2 to model and blend basic shapes including conics and circle arcs.The C-type spline belongs to theωB-spline category of splines that are capable of blending polynomial,trigonometric and hyperbolic functions.Commonly used basic shapes can be exactly represented by these types of splines.We derive explicit formulas for the convenience of modeling the basic curves.The entire blending curve is C^1-continuous.In comparison with the existing best blending method by rational G^2 splines,which are rational splines of degree 3,the proposed method allows simpler representation and blending of the basic curves,and it can represent numerous basic shapes including the hyperbolic types.We also design a subdivision method to generate blending curves;this method is precise for the basic curves and approximate for the blending sections.The subdivision process is efficient for modeling and rendering.It has also proven to be C^1-continuous by the asymptotically equivalent theory and the continuity of stationary subdivision method.In addition,we extend the proposed methods to cases involving the modeling and blending of basic surfaces.We provide many examples that illustrate the merits of our methods.
文摘Based on the butterfly subdivision scheme and the modified butterfly subdivision scheme, an improved butterfly subdivision scheme is proposed. The scheme uses a small stencil of six points to calculate new inserting vertex, 2n new vertices are inserted in the 2n triangle faces in each recursion, and the n old vertices are kept, special treatment is given to the boundary, achieving higher smoothness while using small stencils is realized. With the proposed scheme, the number of triangle faces increases only by a factor of 3 in each refinement step. Compared with the butterfly subdivision scheme and the modified butterfly subdivision scheme, the size of triangle faces changes more gradually, which allows one to have greater control over the resolution of a refined mesh.
基金This work is supported by NKBRSF on Mathematical Mechanics(G1998030600),the National Natural Science Foundation of China(19971087,69603009)and the Doctoral Program(20010358003)and TRAPOYT of MOE,China.
文摘Quadrics are of basic importance in Computer Graphics and Computer Aided Design. In this paper,we design a subdivision scheme based on the method suggested by G. Morin and J. Warren to generate conics and quadrics conveniently. Given the control polygon(poly-hedron),the corresponding ellipse (ellipsoid)can be generated. The hyperbolas and hyperboloids are generated based on the generation of ellipses and ellipsoids by a simple transformation. The method in this paper is much simpler and easier to apply than those given by Eugenia Montiel et al.
文摘In order to study the constitutive behavior of concrete in mesoscopic level, a new method is proposed in this paper. This method uses random polygon particles to simulate full grading broken aggregates of concrete. Based on computational geometry, we carry out the automatic generation of the triangle finite element mesh for the model of random polygon particles of concrete. The finite element mesh generated in this paper is also applicable to many other numerical methods.