1. bookAHEAD OF PRINT
Journal Details
License
Format
Journal
First Published
01 Jan 2016
Publication timeframe
2 times per year
Languages
English
access type Open Access

Constructing Artistic Surface Modeling Design Based on Nonlinear Over-limit Interpolation Equation

Published Online: 27 Jul 2021
Page range: -
Received: 12 Apr 2021
Accepted: 07 Jul 2021
Journal Details
License
Format
Journal
First Published
01 Jan 2016
Publication timeframe
2 times per year
Languages
English
Abstract

The digital and physical methods of establishing minimal curved surfaces are the basis for realizing the design of the minimal curved surface modeling structure. Based on this research background, the paper showed an artistic surface modeling method based on nonlinear over-limit difference equations. The article combines parameter optimization and 3D modeling methods to model the constructed surface modeling. The research found that the nonlinear out-of-limit difference equation proposed in the paper is more accurate than the standard fractional differential equation algorithm. For this reason, the method can be extended and applied to the design of artistic surface modeling.

Keywords

MSC 2010

Introduction

So far, the main breakthroughs of Computer Aided Geometric Design (CAGD) are the Bessel surface patch, Cons surface patch theory, and B-spline method. These methods have become the core algorithms of much well-known computer-aided design (CAD) systems. These algorithms are used to construct surfaces with various shapes ranging from airplanes, cars, and hulls to toys, jewelry, etc., and have achieved significant results. However, the above methods are far from meeting the needs of various applications, from conceptual design to the final, the finished product process is still too long, expensive, and complicated [1]. Ignoring specific algorithms and considering abstractly, both the Non-Uniform Rational B-Splines (NURBS) method and the Bessel method extended Schoenberg's technology to two dimensions 50 years ago. These methods require that the type value points of the constructed surface are relatively evenly distributed and the topology is quadrilateral. However, in actual production, there will be various situations different from this.

The type value points are obtained by calculation or measurement, and the data is not correct, which makes it challenging to construct a smooth surface. Therefore, we must perform data preprocessing or provide a method to control the shape of the character.

In reverse engineering, the model value points come from physical measurement. It isn’t easy to obtain a grid of model value points with uniform distribution and quadrilateral topology [2]. At this time, the hash point interpolation method can be used to construct the surface.

The construction of a non-four-sided surface. To construct an irregular n-sided surface, an appropriate point can be selected on it. Then we divide it into n four-sided surface patches and then use the traditional method to construct. Of course, it can also use the technique proposed by Gregory.

The structure of closed surface entities. At present, to construct a surface entity, we first need to build a single surface and then stitch it together according to the continuous condition. This process is complicated [3]. If we use free-form deformation technology to construct and modify surface entities easily, it will be helpful to make sculpture models.

The construction and control of the entire surface. If we use traditional methods to construct the character, users cannot control the surface flexibly. However, the way of solving partial di

This article intends to explain the basic principles, potential advantages, and problems of hash point interpolation surface, Gregory surface patch, FFD (Free Form Deformation), and PDE (Partial Differential Equation) methods.

Nonlinear overrun hash point interpolation surface

Construct a smooth binary function S for the given hash data (xi, yi, zi), i = 1,2, …, n and make S(xi, yi) = zi, i = 1,2,…, n. There are two methods for hash data interpolation: distance-weighted interpolation and triangular interpolation. The former is suitable for multivariate interpolation. Triangular interpolation is mainly used in geometric modeling. The construction process of hash data triangular interpolation surface patch is as follows

Triangulation and optimization of hash data. It includes data preprocessing, primary triangulation, main triangulation, boundary clipping, and triangulation refinement.

The structure of cubic curve network and triangular surface. First, estimate the average vector at each vertex, and then construct the tangent plane of each vertex according to the intermediate vector. Next, calculate the internal control points of each side close to the vertex on the tangent plane. Thus, a control polygon can be constructed relative to each side [4]. Therefore, a total of 9 control vertices with three boundaries can be obtained for each triangle, and the control vertices inside the triangle can be calculated according to the existing formula.

Splice two adjacent triangular curved surfaces according to the given persistent condition. The triangular surface patch obtained by (2) can only achieve C° continuity when splicing. To achieve C1 or G1 continuity between two adjacent surface patches, it is necessary to ensure that the three triangle pairs shown in Figure 1 are all coplanar. The key is to calculate the internal control points (the circle in Figure 1) based on the cross-boundary guide vector.

The division of triangles. When the three boundaries of the triangular surface patch need to be continuously spliced with the other three surface patches simultaneously, the calculation of the internal control points will be contradictory. To increase the degree of freedom of splicing, we need to subdivide the triangular surface patch into three sub-surface patches. We make each triangular surface patch contact only one sub-triangular surface patch.

Adjustment of the control vertices inside the sub-triangular surface patch [5]. The triangular surface patch is divided into three sub-surface patches. The control vertices in each sub-surface patch can be calculated according to the requirements of cross-border continuity. To ensure the original triangular surface patch is a continuous surface, we should modify the corresponding control vertices of each sub-triangular surface patch.

Generate a straight hash point interpolation surface.

Fig. 1

Joining two adjacent triangular surface patches.

If the triangular surface is required to achieve C2 or G2 continuity, it should be upgraded and divided according to the differential method first, and then all the control vertices that need to be adjusted are calculated according to the continuity requirements. This method is more complex and computationally expensive.

The basic theory of hash point interpolation surface is relatively mature, but the research work is mainly carried out by mathematicians [6]. Therefore, if this method is practical in CAD, all aspects need to work together.

Gregory surface patch

Gregory uses the following expression to define a triangular surface patch (shown in Figure 2). r(V)=i=13γi(V){ei(si)+λiti(si)}\begin{equation} r(V) = \sum_{i = 1}^{3}{\gamma_{i}(V)\{ e_{i}}(s_{i}) + \lambda_{i}t_{i}(s_{i})\}\end{equation}

Fig. 2

Triangular surface patch.

In the formula: represents the coordinate of the center of gravity. represents the parameter of the curve. s1=λ3/(λ2+λ3)s2=λ1/(λ3+λ1)s3=λ2/(λ1+λ2)\begin{equation} \begin{array}{r} s_{1}=\lambda_{3} /\left(\lambda_{2}+\lambda_{3}\right) \\ s_{2}=\lambda_{1} /\left(\lambda_{3}+\lambda_{1}\right) \\ s_{3}=\lambda_{2} /\left(\lambda_{1}+\lambda_{2}\right) \end{array}\end{equation}γi(V)={(1/λi)2/j(1/λj)2λi01λj=0\begin{equation} \gamma_{i}(V)=\left\{\begin{array}{cl} \left(1 / \lambda_{i}\right)^{2} / \sum_{j}\left(1 / \lambda_{j}\right)^{2} & \lambda_{i} \neq 0 \\ 1 & \lambda_{j}=0 \end{array}\right.\end{equation}

ei(si) represents the boundary curve i, sorted in a counterclockwise direction. ti(si) represents the cross-boundary derivative of the boundary curve i, sorted in counterclockwise order. According to the same principle, the n-sided surface patch can be expressed as r(V)=i=13γi(V){ei(si)+λidi(si)},si[0,1]\begin{equation} r(V)=\sum_{i=1}^{3} \gamma_{i}(V)\left\{e_{i}\left(s_{i}\right)+\lambda_{i} d_{i}\left(s_{i}\right)\right\}, s_{i} \in[0,1]\end{equation}

In the formula: V = {λi, i = 1,2,···, n}, si = λi−1/(λi−1 +λi−1). γi(V)=(j1mλi)2/{k=1m(jkmλj)2}\gamma_{i}(V) = (\prod\nolimits_{j \neq 1}^{m}{\lambda_{i})^{2}/\{\sum\nolimits_{k = 1}^{m}(\prod\nolimits_{j \neq k}^{m}{\lambda_{j})^{2}}}\} Represents the convex combination factor. vi represents the expansion factor. ei(si) means the boundary curve i, which is sorted counterclockwise. di(si) represents the cross-border slope of the boundary curve i, which is sorted in a counterclockwise direction [7]. The specific information is shown in Figure 3.

Fig. 3

Gregory patch method.

For irregular surfaces, the Gregory method is more intuitive to describe the character, and countless surface patches can be spliced into a single character. The advantages are apparent, but more practice is needed

Free-form deformation (FFD) of three-dimensional objects

Mathematically, FFD can be defined by the Bernstein polynomial of three variables of tensor product [8]. First, define a local coordinate system X0 − STU on the cuboid, as shown in Figure 4. Then the coordinate X(s.t.u) of any point in the coordinate system is X=X0+sS+tT+uU\begin{equation} X = X_{0} + sS + tT + uU\end{equation}

Fig. 4

Free-form deformation of a three-dimensional object.

The parameter in the formula: s.t.ucan be expressed as a vector s=T×U(XX0)T×USt=U×S(XX0)U×STu=S×T(XX0)S×TU\begin{equation} \begin{gathered} s=\frac{T \times U\left(X-X_{0}\right)}{T \times U S} \\ t=\frac{U \times S\left(X-X_{0}\right)}{U \times S T} \\ u=S \times T\left(X-X_{0}\right) S \times T U \end{gathered}\end{equation}

For any point in the cuboid, 0 < s < 1,0 < t < 1,0 < u < 1 Set (l + 1)(m + 1)(n + 1) points in the cuboid body and on the surface, and generate l + 1, m + 1 and n + 1 planes in the S, T, and U directions through these points. l + 1, m + 2n = 3. The intersections constitute the control vertex grid, and the coordinates of the grid points are Pi,jk=X0+ilS+imT+knU,(i=1,2,,l;j=1,2,,m;k=1,2,,n)\begin{equation} P_{i, j k}=X_{0}+\frac{i}{l} S+\frac{i}{m} T+\frac{k}{n} U,(i=1,2, \cdots, l ; j=1,2, \cdots, m ; k=1,2, \cdots, n)\end{equation}

Move Pi,jk to deform the object. This indicator is to calculate the position Xffd of any point after deformation. (s, t, u) should be calculated first, and then the vector value three-variable Bersnetin polynomial should be calculated. Xffd=i=0l(lil)(1s)liSi[j=0m(mi)(1t)mjtj[k=0m(mk)(1u)mkukPijk]]\begin{equation} X_{f f d}=\sum_{i=0}^{l}\left(l_{i}^{l}\right)(1-s)^{l-i} S^{i}\left[\sum_{j=0}^{m}\left(\begin{array}{l} m \\ i \end{array}\right)(1-t)^{m-j} t^{j}\left[\sum_{k=0}^{m}\left(\begin{array}{c} m \\ k \end{array}\right)(1-u)^{m-k} u^{k} P_{i j k}\right]\right]\end{equation}

Xffd is the Cartesian coordinate vector of the moving point. Pijk is the coordinate vector of the control point. If required to splice two deformed bodies, it should meet a particular continuity requirement [9]. When the control points at the splicing place of the two objects coincide, the C0 continuity requirement can be met. The condition that two things meet the continuous first-order derivative at the joint is X1(1,t,u)s=X2(0,t,u)sX1(1,t,u)t=X2(0,t,u)tX1(1,t,u)u=X2(0,t,u)u\begin{equation} \begin{aligned} \frac{X_{1}(1, t, u)}{\partial s} &=\frac{\partial X_{2}(0, t, u)}{\partial s} \\ \frac{X_{1}(1, t, u)}{\partial t} &=\frac{\partial X_{2}(0, t, u)}{\partial t} \\ \frac{X_{1}(1, t, u)}{\partial u} &=\frac{\partial X_{2}(0, t, u)}{\partial u} \end{aligned}\end{equation}

These conditions are a direct extension of the continuous requirements of curves and surfaces. The FFD method can also be used for local deformation. To achieve the constant local deformation Gk, we only need to ensure that the control points on the K sections adjacent to the joint surface remain unchanged. As shown in Figure 5, the thick line in the figure is the joint surface [10]. Thus, the deformation process is similar to making a sculpture model. In addition, under certain conditions, the object's volume can be guaranteed to remain unchanged before and after deformation.

Fig. 5

Local deformation of FFD method.

From an engineering point of view, the characteristic of FFD is that it can complete global or local deformation under continuous conditions. This method can construct surfaces that require aesthetic appearance, certain transition surfaces, and specific functional surfaces. This method is easy to use. In expression (9), a character is generated when a particular parameter is zero. The curve is generated when the two parameters are zero. Therefore, it can be used to express curves, surfaces and curved bodies uniformly. However, FFD also has its limitations. For example, it is challenging to generate a generally rounded transition surface; simultaneously, it requires a lot of calculation.

Using partial differential equations (PDE) to construct surfaces

From 1992 to 1993, the CAD laboratory of MIT completed the research on the construction of free-form surfaces by solving partial differential equations. The resulting surface, like an elastic film boundary, can be fixed according to any space curve. Various points on the curved surface can also be improved. The curved surface deforms like soap bubbles when subjected to positive or negative pressure. Similar to the physical film surface equation, the weighted sum of area and curvature is minimized. Therefore, the surface generated by it appears smooth and attractive [11]. People can control the shape of the film in different ways: (1) Fix points on the surface or curve. For example, when constructing an engine intake pipe, one can force a curved surface to pass through various sections given along the ridgeline. (2) Apply positive or negative pressure to the curved surface to deform the curved surface. The pressure can be distributed uniformly or non-uniformly. (3) Change the weight of the surface area and curvature. Changing the weight function causes the membrane's stiffness to change, and the curved surface is deformed under a given pressure. To change the shape of the surface, the operator does not need to enter a specific pressure value or weighting factor. They can use the cursor on the graphical interface to adjust the pressure on the surface until a satisfactory surface is obtained.

Physics, mechanics, engineering technology, and other natural sciences have proposed many partial differential equations. An equation containing the partial derivative of an unknown function is called a partial differential equation. If there is more than one equation, these equations are called a partial differential equation system [12].

We suppose that the function u has the continuous partial derivatives of each order appearing in the partial differential equation in the region D. If u is substituted into this partial differential equation to make it an identity in the region D, then u is called the partial differential equation in Solution in area D. u = u(x1, x2,, xn) is a curved surface in the n + 1 dimensional space (u, x1, x2,, xn). We call it the integral surface of the partial differential equation [13]. In geometric modeling, the boundary value problem is concerned. The partial differential equation is a definite solution problem with only boundary value conditions but no initial conditions at the same time. It is also required to generate a smooth surface. Therefore, we should use elliptic partial differential equations. We need that a transition surface z = f (x, y), z is constructed between two semicircles and one semicircle in the region {Ω : 0 ≤ x ≤ 2;0 ≤ y ≤ 2} to meet the following boundary conditions: {z(0,0y2)=+(1(y1)2)1/2z(2,0y1)=+(1/4(y1/2)2)1/2z(2,1y2)=(1/4(y3/2)2)1/2z(0x2,0)=0z(0y2,2)=0\begin{equation} \left\{\begin{array}{l} z(0,0 \leq y \leq 2)=+\left(1-(y-1)^{2}\right)^{1 / 2} \\ z(2,0 \leq y \leq 1)=+\left(1 / 4-(y-1 / 2)^{2}\right)^{1 / 2} \\ z(2,1 \leq y \leq 2)=-\left(1 / 4-(y-3 / 2)^{2}\right)^{1 / 2} \\ z(0 \leq x \leq 2,0)=0 \\ z(0 \leq y \leq 2,2)=0 \end{array}\right.\end{equation}

The surface that meets this condition can be generated by the following formula: x,y2f=(zx2+zy2)f(x,y)=0\begin{equation} \nabla_{x, y}^{2} f=\left(\frac{\partial z}{\partial x^{2}}+\frac{\partial z}{\partial y^{2}}\right) f(x, y)=0\end{equation}

Due to the smoothness of the operator ∇2, the generated surface is also smooth. It has no local maximum and minimum [14]. To ensure the tangential continuity between the generated surface and the adjacent surface, we need to use the fourth-order elliptic partial differential equation D4F=(zu2+zv2)(zu2+a2zv2)\begin{equation} D^{4} F=\left(\frac{\partial z}{\partial u^{2}}+\frac{\partial z}{\partial v^{2}}\right)\left(\frac{\partial z}{\partial u^{2}}+a^{2} \frac{\partial z}{\partial v^{2}}\right)\end{equation}

Where a is the shape control parameter, when the curvature is required to be continuous, a sixth-order elliptic partial differential equation is required. The advantage of using partial differential equations to construct the surface is that the entire surface can be built, the control is flexible, and the operation is convenient. The problem is that the amount of calculation is large, and the current computer can not generate the required surface at a reasonable cost. On the other hand, the modeling method is relatively new, and the application experience is insufficient.

Conclusion

This article briefly describes four surface modeling methods. The PDE method is the latest, and the original literature has been about six years. On the other hand, the hash point interpolation surface has the most extended history, which has reached more than 20 years. On the other hand, these methods have not been widely promoted in CAD, caused by various reasons. On the one hand, the mature application of the new process itself takes time. On the other hand, because it is difficult to modify the mathematical model of the existing CAD art surface modeling, the developers are not interested in the improvement suggestions. This is one reason why the internationally renowned software still uses the traditional method to construct the surface. In this regard, the situation in China is different. Chinese software industry is still immature, so it is possible to apply new modeling methods to CAD art surface modeling based on research and practice. The user interface is essential for the existing CAD art surface modeling's potential functions and should be paid close attention.

Fig. 1

Joining two adjacent triangular surface patches.
Joining two adjacent triangular surface patches.

Fig. 2

Triangular surface patch.
Triangular surface patch.

Fig. 3

Gregory patch method.
Gregory patch method.

Fig. 4

Free-form deformation of a three-dimensional object.
Free-form deformation of a three-dimensional object.

Fig. 5

Local deformation of FFD method.
Local deformation of FFD method.

Yuksel, C. A Class of C 2 Interpolating Splines. ACM Transactions on Graphics (TOG)., 2020. 39(5): 1–14YukselC.A Class of C 2 Interpolating SplinesACM Transactions on Graphics (TOG)2020395114Search in Google Scholar

Nogueira Martins, R., Ferreira Lima Dos Santos, F., De Moura Araújo, G., De Arruda Viana, L., & Fim Rosas, J. T. Accuracy assessments of stochastic and deterministic interpolation methods in estimating soil attributes spatial variability. Communications in Soil Science and Plant Analysis., 2019. 50(20): 2570–2578Nogueira MartinsR.Ferreira Lima Dos SantosF.De Moura AraújoG.De Arruda VianaL.Fim RosasJ. T.Accuracy assessments of stochastic and deterministic interpolation methods in estimating soil attributes spatial variabilityCommunications in Soil Science and Plant Analysis2019502025702578Search in Google Scholar

Habib, M., Alzubi, Y., Malkawi, A., & Awwad, M. Impact of interpolation techniques on the accuracy of large-scale digital elevation model. Open Geosciences., 2020. 12(1): 190–202HabibM.AlzubiY.MalkawiA.AwwadM.Impact of interpolation techniques on the accuracy of large-scale digital elevation modelOpen Geosciences2020121190202Search in Google Scholar

Gahbiche, M. A., Boudhaouia, S., Giraud, E., Ayed, Y., Santo, P. D., & Salem, W. B. A finite element simulation of the incremental sheet forming process: a new method for G-code implementation. International Journal of Materials and Product Technology., 2020. 61(1): 68–86GahbicheM. A.BoudhaouiaS.GiraudE.AyedY.SantoP. D.SalemW. B.A finite element simulation of the incremental sheet forming process: a new method for G-code implementationInternational Journal of Materials and Product Technology20206116886Search in Google Scholar

Ning, Y., Liu, Y., Zhang, Y., & Zhang, C. Adaptive image rational upscaling with local structure as constraints. Multimedia Tools and Applications., 2019. 78(6): 6889–6911NingY.LiuY.ZhangY.ZhangC.Adaptive image rational upscaling with local structure as constraintsMultimedia Tools and Applications201978668896911Search in Google Scholar

Zhao, F., Li, J., Xiao, X., & Feng, X. The characteristic RBF-FD method for the convection-diffusion-reaction equation on implicit surfaces. Numerical Heat Transfer, Part A: Applications., 2019. 75(8): 548–559ZhaoF.LiJ.XiaoX.FengX.The characteristic RBF-FD method for the convection-diffusion-reaction equation on implicit surfacesNumerical Heat Transfer, Part A: Applications2019758548559Search in Google Scholar

Liang, H., Zhang, Q., Fu, C., Liang, F., & Sun, Y. Surface Modelling of Jun Ware Based on Ordinary Differential Equations. Traitement du Signal., 2019. 36(1): 53–58LiangH.ZhangQ.FuC.LiangF.SunY.Surface Modelling of Jun Ware Based on Ordinary Differential EquationsTraitement du Signal20193615358Search in Google Scholar

Xia, S., Chen, D., & Wang, R. A breakline-preserving ground interpolation method for MLS data. Remote Sensing Letters., 2019. 10(12): 1201–1210XiaS.ChenD.WangR.A breakline-preserving ground interpolation method for MLS dataRemote Sensing Letters2019101212011210Search in Google Scholar

Jornet, M. Modeling of Allee effect in biofilm formation via the stochastic bistable Allen–Cahn partial differential equation. Stochastic Analysis and Applications., 2021. 39(1): 22–32JornetM.Modeling of Allee effect in biofilm formation via the stochastic bistable Allen–Cahn partial differential equationStochastic Analysis and Applications20213912232Search in Google Scholar

Talbi, N., Dhahbi, A. B., Boulaaras, S., Baltache, H., & Alnegga, M. A Two Dimensional Mathematical Model of Heat Propagation Equation and its Applications. Computational Mathematics and Modeling., 2020. 31(3): 338–354TalbiN.DhahbiA. B.BoulaarasS.BaltacheH.AlneggaM.A Two Dimensional Mathematical Model of Heat Propagation Equation and its ApplicationsComputational Mathematics and Modeling2020313338354Search in Google Scholar

Acar, E., Zgi, A., & Serenbay, S. K. Note on jakimovski-leviatan operators preserving e–x. Applied Mathematics and Nonlinear Sciences., 2019. 4(2): 543–550AcarE.ZgiA.SerenbayS. K.Note on jakimovski-leviatan operators preserving e–xApplied Mathematics and Nonlinear Sciences201942543550Search in Google Scholar

Rahaman, H., Kamrul Hasan, M., Ali, A. & Shamsul Alam, M. Implicit Methods for Numerical Solution of Singular Initial Value Problems. Applied Mathematics and Nonlinear Sciences., 2020. 6(1): 1–8RahamanH.Kamrul HasanM.AliA.Shamsul AlamM.Implicit Methods for Numerical Solution of Singular Initial Value ProblemsApplied Mathematics and Nonlinear Sciences20206118Search in Google Scholar

Aghili, A. Complete Solution For The Time Fractional Diffusion Problem With Mixed Boundary Conditions by Operational Method. Applied Mathematics and Nonlinear Sciences., 2020. 6(1): 9–20AghiliA.Complete Solution For The Time Fractional Diffusion Problem With Mixed Boundary Conditions by Operational MethodApplied Mathematics and Nonlinear Sciences202061920Search in Google Scholar

Eshbaeva, U., Djalilov, A., Turaev, F., & Dilafruz, S. Mathematical Modeling of Ink Transition of Print Product. European Journal of Molecular & Clinical Medicine., 2021. 8(1): 709–717EshbaevaU.DjalilovA.TuraevF.DilafruzS.Mathematical Modeling of Ink Transition of Print ProductEuropean Journal of Molecular & Clinical Medicine202181709717Search in Google Scholar

Recommended articles from Trend MD

Plan your remote conference with Sciendo