Cad system and cad program

TECHNICAL FIELD

The present invention relates to a computer aided design system and a computer aided design program which transform the shape of a member into an objective curved surface shape.

BACKGROUND ART

Today, there is a desire to shorten processes from planning, to design and production to respond to consumer demand. In order to improve the efficiency of design and production processes, the use of CG (Computer Graphics) and CAD (Computer Aided Design) systems is popular. In order to depict shapes having complex curved lines or curved surface shapes, such as for motor vehicles, domestic electric appliances or the like, on a computer, the following processing methods have conventionally existed.

The first is solid modeling, where simple shapes called primitives, are held in a computer, and operations to combine the shapes with each other are repeated in order to express complex shapes. A primitive is for example a column, a cube, a hexahedron, a torus, a ball, or the like, and in the solid modeling, shapes are represented by set operations on these primitives. Therefore, in order to produce a complex shape many steps are required and precise calculations are required.

The second is surface modeling, which utilizes an algorithm such as a bezier, b-spline, rational bezier, NURBS (Non-Uniform Rational b-spline), or the like in order to perform operations such as cutting or connecting lines or surfaces, and by repeating these operations complex free curved lines or curved surfaces are represented.

However, even with a model where there are no problems from the point of representation with the solid model or surface model described above, in some cases problems may occur in the case where it is used by a downstream application such as CAM, CAE, or the like. This is caused by differences between the support element to be supported by the produced computer graphics, and the support element to be supported by the other computer graphics, computer aided design, and downstream applications, and differences in shape definition, or the like. The model is corrected via an application such as a translator which modifies these differences (Japanese Unexamined Patent Publication Nos. 2001-250130, Hei 11-65628, Hei 10-69506, Hei 4-134571, Hei 4-117572, and Hei 1-65628).

DISCLOSURE OF INVENTION

However, the above described correcting operations are extremely ineffective for shortening the design and production processes. The reasons for requiring the corrections vary depending on each case, but the point which becomes a problem, particularly in the production stage is that the representations of all curved lines and curved surfaces are approximated by Euclidean geometry in a conventional computer graphics or computer aided design system. For example, in the case where a saddle type tabulated cylinder surfaces shown in FIG. 6 is generated by a sweep operation, a long line in the lower slope part of the saddle and a short line in the central part of the saddle exist. Therefore, this sweep operation is a transformation accompanied with graphical expansion and contraction in order to maintain the continuity of the curved surface generated. However, in the conventional computer graphics or computer aided design system, this expansion and contraction is not considered, and the internal representation is approximately represented as a cylinder type. Therefore, if the computer graphics model or computer aided design model which are actually approximately represented by such Euclidean geometry, are passed to a CAE, the error occurring here becomes a problem in production.

The present invention has been devised to solve such problems, with an object of providing a computer aided design system and a computer aided design program which can greatly utilize a computer graphics model or a computer aided design model, and can improve the efficiency of design and production processes.

A computer aided design system of the present invention comprises: a point sequence information extraction device which extracts a plurality of point sequences on a curved surface; a dividing device which generates a curved surface from the point sequences using another computer aided design system, and divides the curved surface into a predetermined number of meshes; a first fundamental form computing device for computing coefficients of the first fundamental form defined by a tangent vector which forms a tangent plane of the mesh; a second fundamental form computing device for computing coefficients of the second fundamental form defined by the tangent vector and a normal vector of the mesh; and a memory device which stores the point sequence information, the coefficients of the first fundamental form and the coefficients of the second fundamental form.

Furthermore, the computer aided design system of the present invention further comprises: a principal curvature computing device which computes a principal curvature of the mesh based on the coefficients of the first fundamental form and coefficients of the second fundamental form; a line of curvature computing device which computes a line of curvature showing a principal direction of the mesh based on the principal curvature; a feature point/feature line analyzing device which extracts a point or a line which become a reference point or a reference line of transformation defined by changing patterns of one or more feature quantities among five feature quantities showing features of the curved surface comprising a Gaussian curvature and a mean curvature computed based on the principal curvature, the principal direction, the line of curvature, and the coefficients of the first fundamental form and coefficients of the second fundamental form; and a girth length computing device which computes a girth length based on a curvature computed from the coefficients of the first fundamental form and coefficients of the second fundamental form.

Moreover, the computer aided design system of the present invention further comprises: a reproducing device which transforms the line of curvature for the girth length in the line of curvature direction, with the feature point or feature line as a transformation reference, and reproduces the mesh or the curved surface.

Furthermore, the computer aided design system of the present invention further comprises: a converting device which extracts a plurality of point sequences on a curved surface from the reproduced mesh or curved surface, and converts the point sequences according to a graphical representation algorithm in another computer aided design system.

A computer aided design program of the present invention executes on a computer: a point sequence information extraction process for extracting a plurality of point sequences on a curved surface; a dividing process for generating a curved surface from the point sequences using another computer aided design system, and dividing the curved surface into a predetermined number of meshes; a first fundamental form computing process for computing coefficients of the first fundamental form defined by a tangent vector which forms a tangent plane of the mesh; a second fundamental form computing process for computing coefficients of the second fundamental form defined by the tangent vector and a normal vector of the mesh; and a storage process for storing the point sequence information, the coefficients of the first fundamental form and the coefficients of the second fundamental form.

Moreover, the computer aided design program of the present invention is a computer aided design program for further executing on a computer: a principal curvature computing process for computing a principal curvature of the mesh based on the coefficients of the first fundamental form and coefficients of the second fundamental form; a line of curvature computing process for computing a line of curvature showing a principal direction of the mesh based on the principal curvature; a feature point/feature line analyzing process for extracting a point or a line which become a reference point or a reference line of transformation defined by changing patterns of one or more feature quantities among five feature quantities showing features of the curved surface comprising a Gaussian curvature and a mean curvature computed based on the principal curvature, the principal direction, the line of curvature, and the coefficients of the first fundamental form and coefficients of the second fundamental form; and a girth length computing process for computing a girth length based on a curvature computed from the coefficients of the first fundamental form and coefficients of the second fundamental form.

Furthermore, the computer aided design program of the present invention is a computer aided design program for further executing on a computer: a reproducing process for transforming the line of curvature for the girth length in the line of curvature direction, with the feature point or feature line as a transformation reference, and reproducing the mesh or the curved surface.

Moreover, the computer aided design program of the present invention is a computer aided design program for further executing on a computer: a converting process for extracting a plurality of point sequences on a curved surface from the reproduced mesh or curved surface, and converting the point sequences according to a graphical representation algorithm in another computer aided design system.

A computer graphics system of the present invention comprises: a point sequence information extraction device which extracts a plurality of point sequences on a curved surface; a dividing device which generates a curved surface from the point sequences using another computer graphics system, and divides the curved surface into a predetermined number of meshes; a first fundamental form computing device for computing coefficients of the first fundamental form defined by a tangent vector which forms a tangent plane of the mesh; a second fundamental form computing device for computing coefficients of the second fundamental form defined by the tangent vector and a normal vector of the mesh; and a memory device which stores the point sequence information, the coefficients of the first fundamental form and the coefficients of the second fundamental form.

Moreover, the computer graphics program of the present invention is a computer graphics program for executing on a computer: a point sequence information extraction process for extracting a plurality of point sequences on a curved surface; a dividing process for generating a curved surface from the point sequences using another computer graphics system, and dividing the curved surface into a predetermined number of meshes; a first fundamental form computing process for computing coefficients of the first fundamental form defined by a tangent vector which forms a tangent plane of the mesh; a second fundamental form computing process for computing coefficients of the second fundamental form defined by the tangent vector and a normal vector of the mesh; and a storage process for storing the point sequence information, the coefficients of the first fundamental form and the coefficients of the second fundamental form.

The present invention demonstrates the following effects.

Since it comprises: the point sequence information extraction device which extracts a plurality of point sequences on a curved surface; the dividing device which generates a curved surface from the point sequences using another computer graphics system or computer aided design system, and divides the curved surface into a predetermined number of meshes; the first fundamental form computing device for computing coefficients of the first fundamental form defined by a tangent vector which forms a tangent plane of the mesh; the second fundamental form computing device for computing coefficients of the second fundamental form defined by the tangent vector and a normal vector of the mesh; and the memory device which stores the point sequence information, the coefficients of the first fundamental form and the coefficients of the second fundamental form, then by adopting a curved surface theory which ensures the continuity of a free-form line/free-form surface, a computer graphics model or a computer aided design model can be widely utilized, and the efficiency of design and production processes can be improved.

Moreover, since it further comprises: the principal curvature computing device which computes a principal curvature of the mesh based on the coefficients of the first fundamental form and the coefficients of the second fundamental form; the line of curvature computing device which computes a line of curvature showing a principal direction of the mesh based on the principal curvature; the feature point/feature line analyzing device which extracts a point or a line which become a reference point or a reference line of transformation defined by changing patterns of one or more feature quantities among five feature quantities showing features of the curved surface comprising a Gaussian curvature and a mean curvature computed based on the principal curvature, the principal direction, the line of curvature, and the coefficients of the first fundamental form and the coefficients of the second fundamental form; and the girth length computing device which computes a girth length based on a curvature computed from the coefficients of the first fundamental form and the coefficients of the second fundamental form, then a computer graphics or computer aided design model analyzed by the curved surface theory can be reproduced and converted into another computer graphics or computer aided design model.

Furthermore, since it further comprises: the reproducing device which transforms the line of curvature for the girth length in the line of curvature direction, with the feature point or the feature line as the transformation reference, and reproduces the mesh or the curved surface, then a computer graphics or computer aided design model analyzed by the curved surface theory can be reproduced.

Moreover, since it further comprises: the converting device which extracts a plurality of point sequences on a curved surface from the reproduced mesh or curved surface, and converts the point sequences according to a graphic expression algorithm in another computer graphics or computer aided design system, then a computer graphics or computer aided design model analyzed by the curved surface theory can be converted into another computer graphics or computer aided design model.

BRIEF DESCRIPTION OF DRAWINGS

<ul id="ul0001" list-style="none" compact="compact"><li>FIG. 1 is a block diagram showing a configuration of a computer aided design system of the present embodiment.</li><li>FIG. 2 is an explanatory diagram showing a situation for dividing a curved surface into an m x n mesh and defining fundamental vectors Su and Sv.</li><li>FIG. 3 is an explanatory diagram showing planes in which a unit tangent vector t and an unit normal vector n extend.</li><li>FIG. 4 is a flowchart showing a processing flow from free-form surface analysis to data transfer, by an analysis program 1.</li><li>FIG. 5 is an explanatory diagram showing an aspect of curvature change.</li><li>FIG. 6 is an explanatory diagram showing classifications of mean curvature and Gaussian curvature.</li><li>FIG. 7 is an explanatory diagram showing isoclinic orthogonal lines.</li><li>FIG. 8 is an explanatory diagram showing principal curvature extremum lines.</li><li>FIG. 9 is an explanatory diagram showing isoclinic extremum lines and aspects of Gaussian curvature distribution.</li><li>FIG. 10 is an explanatory diagram showing lines of curvature.</li></ul>

BEST MODE FOR CARRYING OUT THE INVENTION

Hereunder is a description of an embodiment of a computer aided design system of the present invention, with reference to the drawings. FIG. 1 is a block diagram showing a configuration of the computer aided design system of the present embodiment. The computer aided design system of the present embodiment comprises; a central processor such as central processing unit (CPU) or the like (not shown), a storage memory such as a ROM, RAM, or the like (not shown), a database 10, a graphic display processing section 11, a display section 12, an output section 13, and a communication section (not shown).

The CPU reads out an analysis program 1, a converting program 2, and a reproducing program 3 stored in a ROM and executes a series of processes related to free-form surface analysis, conversion, and reproduction. The RAM is a semiconductor memory in which the CPU primarily stores data.

The analysis program 1 is a program which executes in the CPU, a process for reading in actual measurement value data 20 of a three-dimensional shaped body by CAT (computer aided testing) or the like, or other computer aided design format data 21 (graphics data represented for example by a surface model such as a solid model, a bezier, a b-spline, a rational bezier, or a NURBS), creating a point sequence information table 30, a table 31 of coefficients of the first fundamental form, and a table 32 of coefficients of the second fundamental form, and then storing this in the database 10.

The point sequence information table 30 comprises point sequence information (u, v) for on a curved surface expressed in a parameter form of;<maths id="math0001" num="(equation 1)"><math display="block"><mrow><mtext>s(u, v) = {x(u, v), y (u, v), z (u, v)} 0≤u, v≤1</mtext></mrow></math><img file="EP1548618A1_D0001.tif" /></maths> as shown in FIG. 2. For example, assuming that u = 0, 1/m, 2/m, to m-1/m (m is a natural number) and v = 0, 1/n, 2/n, to n-1/n (n is a natural number), the curved surface shown in FIG. 2 is divided into an m x n mesh. In this case, the point sequence information (u, v) becomes an mn data sequences from the mesh ID1 to IDmn.

The table 31 of coefficients of the first fundamental form comprises coefficients of the first fundamental form E, F, and G derived from the following equations. In the case where u and v described above have a functional relation, then s (u, v) denotes a curved line on the curved surface, the partial derivative ∂s/∂u = Su denotes a tangent vector of a curved line of u = constant, and the partial derivative ∂s/∂v = Sv denotes a tangent vector of a curved line of v = constant. At this time, the fundamental vectors Su and Sv form a tangent plane of the curved surface. Moreover, a vector ds linking two points on the curved surface, from s (u, v) to s (u+du, v+dv), is represented by:<maths id="math0002" num="(equation 2)"><math display="block"><mrow><msub><mrow><mtext>ds = s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext>du + s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><mtext>dv</mtext></mrow></math><img file="EP1548618A1_D0002.tif" /></maths>

Here, the square of the absolute value of ds is represented by:<maths id="math0003" num="(equation 3)"><math display="block"><mrow><msup><mrow><mtext>(ds)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext> = ds · ds = s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> (du)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext> + 2s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> · s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext>dudv + s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> (dv)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1548618A1_D0003.tif" /></maths>

The coefficients of the first fundamental form described above are defined from the fundamental vector of the curved surface by the following equation:<maths id="math0004" num="(equation 4)"><math display="block"><mrow><msub><mrow><mtext>E = s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>, F = s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> · s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext>, G = s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1548618A1_D0004.tif" /></maths>

The coefficients of the first fundamental form E, F, and G described above are uniquely determined for the respective meshes in this way. The table 31 of coefficients of the first fundamental form stores the values for the respective meshes ID1 to IDmn.

Moreover, combining the abovementioned equation 3 and equation 4 gives:<maths id="math0005" num="(equation 5)"><math display="block"><mrow><msup><mrow><mtext>ds</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> = E (du)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> + 2Fdudv + G(dv)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1548618A1_D0005.tif" /></maths>

The table 32 of coefficients of the second fundamental form comprises coefficients of the second fundamental form L, M, and N derived from the following equations. Assuming that ω is an angle between the fundamental vectors Su and Su, then their inner product F, and the absolute value H of the vector product of the fundamental vectors, are represented as follows using the coefficients of the first fundamental form:<maths id="math0006" num="(equation 6)"><math display="block"><mrow><msub><mrow><mtext>F = |s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext>| · |s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><mtext>|cosω = (√EG)cosω</mtext></mrow></math><img file="EP1548618A1_D0006.tif" /></maths><maths id="math0007" num="(equation 7)"><math display="block"><mrow><msub><mrow><mtext>H = |s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext>×s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext>| = |s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext>| · |s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><mtext>| sinω</mtext><mspace linebreak="newline" /><mtext> = </mtext><msqrt><msup><mrow><mtext>EG(1 - cos</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>ω)</mtext></msqrt><mtext> = </mtext><msqrt><msup><mrow><mtext>EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></msqrt></mrow></math><img file="EP1548618A1_D0007.tif" /></maths>

Then using this calculated value H, the unit normal vector n on the curved surface is represented by:<maths id="math0008" num="(equation 8)"><math display="block"><mrow><mtext>n = </mtext><mfrac><mrow><msub><mrow><mtext>(s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> × s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><mtext>H</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0008.tif" /></maths>

Moreover, as shown in FIG. 3, the pencil of lines of the tangent vectors at a point P on the curved surface exists in this tangent plane, and a unit tangent vector t of these is represented by the following equation:<maths id="math0009" num=""><img file="EP1548618A1_D0009.tif" /></maths>

The plane determined by t and n as shown in FIG. 3 is called a normal plane.

The curvature κ at the point P on this normal section plane is called a normal curvature. Differentiating t for the arc length s on the normal section plane gives:<maths id="math0010" num=""><img file="EP1548618A1_D0010.tif" /></maths>

Multiplying both equations by the normal vector, and introducing the following coefficients of the second fundamental form:<maths id="math0011" num="(equation 11)"><math display="block"><mrow><msub><mrow><mtext>L = n · s</mtext></mrow><mrow><mtext>uu</mtext></mrow></msub><msub><mrow><mtext>, M = n · s</mtext></mrow><mrow><mtext>uv</mtext></mrow></msub><msub><mrow><mtext>, N = n · s</mtext></mrow><mrow><mtext>vv</mtext></mrow></msub></mrow></math><img file="EP1548618A1_D0011.tif" /></maths>    gives:<maths id="math0012" num=""><img file="EP1548618A1_D0012.tif" /></maths>

The coefficients of the second fundamental form L, M, and N described above are uniquely determined for the respective meshes in this way. The table 32 of coefficients of the second fundamental form stores the values for the respective meshes ID1 to IDmn.

If equation 5 is substituted in equation 12, the following equation is obtained.<maths id="math0013" num="(equation 13)"><math display="block"><mrow><mtext>κ = </mtext><mfrac><mrow><msup><mrow><mtext>L (du)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> + 2 Mdudv + N (dv)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><msup><mrow><mtext>E (du)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> + 2 Fdudv + G (dv)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP1548618A1_D0013.tif" /></maths>

From the above, the normal curvature is computed from the coefficients of the first fundamental form and the coefficients of the second fundamental form.

The converting program 2 is a program which executes on a computer a process for; reading out the necessary information for a free-form surface from the point sequence information table 30, the table 31 of coefficients of the first fundamental form, and the table 32 of coefficients of the second fundamental form, then creating free-form surface data, and transforming this into a form which an other computer aided design application can interpret.

The reproducing program 3, similarly to the converting program 2, is a program which executes on a computer a process for; reading out the necessary information for a free-form surface from the point sequence information table 30, the table 31 of coefficients of the first fundamental form, and the table 32 of coefficients of the second fundamental form, then creating free-form surface data, and outputting to the graphic display processing section 11.

The database 10 stores the above described point sequence information table 30, the table 31 of coefficients of the first fundamental form, and the table 32 of coefficients of the second fundamental form, and writes in the output result of the analysis program 1 in association with a mesh ID described later.

The graphic display processing section 11 performs graphic display processing on the output results from the reproducing program, and other computer aided design applications.

The display section 12 displays the output results of the graphic display processing section 11.

The output section 13 outputs the output results of the graphic display processing section 11 to the communication section, other recording media, or the like. The communication section transfers data such as the point sequence information, the coefficients of the first fundamental form, and the coefficients of the second fundamental form stored in the database 1 to other severs or clients via a network such as a LAN, the Internet, or the like.

Next is a description of a series of processing flows related to the free-form surface analysis, conversion, and reproduction by the computer aided design system of the present embodiment, with reference to the drawings. FIG. 4 is a flowchart showing a processing flow from free-form surface analysis to data transfer, by the analysis program 1.

By the user's operation, the CPU receives an analyze command for the actual measurement value data 20 or other computer aided design format data 21, reads out the analysis program 1 from ROM, and executes the free-form surface analyzing process. Firstly, the CPU performs a process for extracting a plurality of point sequences on a curved surface such as a two-dimensional NURBS surface, bicubic surface, or the like, held by the actual measurement value data 20 or the other computer aided design format data 21. Then, a curved surface is generated from this point sequence using the other computer aided design system (step S1 in FIG. 4), and the curved surface is divided into a predetermined mn mesh number as shown in FIG. 2, after which the respective mesh parts are standardized by fundamental vectors Su and Sv. The point sequence information (u, v) generated during the standardization is written in the point sequence information table 30 held by the database 10, in association with the mesh ID.

Next, the CPU executes differential geometric analysis processing. That is, it performs processing for computing the coefficients of the first fundamental form E, F, and G defined by the fundamental vectors Su and Sv which form the tangent plane of the mesh. The computed coefficients of the first fundamental form E, F and G, similarly to the point sequence information, are written in the table 31 of coefficients of the first fundamental form held by the database 10, in association with the mesh ID. Moreover, the CPU performs a process for computing the coefficients of the second fundamental form L, M, and N defined by the fundamental vectors Su and Sv and an unit normal vector n of the mesh. The computed coefficients of the second fundamental form L, M, and N, similarly to the coefficients of the first fundamental form E, F and G, are written in the table 32 of coefficients of the second fundamental form held by the database 10, in association with the mesh ID.

Moreover, the CPU performs processing for computing an integrable condition which is a condition where the differential equation representing the above described mesh is continuous at the respective boundaries of the mesh, in other words, a condition where this differential equation has a unique solution.

Now, the curved surface coordinates (u, v) described above are substituted by (u1, u2) and the point is made p (u1, u2). If a curved line formed by fixing u2 and moving u1 is called a u1 curved line, and a curved line formed by fixing u1 and moving u2 is called a u2 curved line, then assuming that p (u1, u2) on the curved surface is the initial point, then the tangent vector along the u1 curved line and the u2 curved line can be calculated as follows:<maths id="math0014" num="(equation 14)"><math display="block"><mrow><msub><mrow><mtext>e</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>∂p</mtext></mrow><mrow><msup><mrow><mtext>∂u</mtext></mrow><mrow><mtext>1</mtext></mrow></msup></mrow></mfrac><msub><mrow><mtext>, e</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>∂p</mtext></mrow><mrow><msup><mrow><mtext>∂u</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP1548618A1_D0014.tif" /></maths>

Then, the unit normal vector n can be calculated from e1 and e2 as follows:<maths id="math0015" num="(equation 15)"><math display="block"><mrow><mtext>n = </mtext><mfrac><mrow><msub><mrow><mtext>e</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> × e</mtext></mrow><mrow><mtext>2</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext>∥e</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> × e</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>∥</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0015.tif" /></maths>

In this way, three vectors {e1, e2, n} are defined for the respective points on the curved surface.

For the respective points, first fundamental quantities E, F and G are defined as follows:<maths id="math0016" num="(equation 16)"><math display="block"><mrow><msub><mrow><mtext>E = ∥e</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext>∥</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>, F = (e</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>, e</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>), G = ∥e</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext>∥</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1548618A1_D0016.tif" /></maths>

Then, a first fundamental tensor (g_ij, i, j = 1, 2) is defined as follows:<maths id="math0017" num="(equation 17)"><math display="block"><mrow><msub><mrow><mtext>g</mtext></mrow><mrow><mtext>11</mtext></mrow></msub><msub><mrow><mtext> = E, g</mtext></mrow><mrow><mtext>12</mtext></mrow></msub><msub><mrow><mtext> = g</mtext></mrow><mrow><mtext>21</mtext></mrow></msub><msub><mrow><mtext> = F, g</mtext></mrow><mrow><mtext>22</mtext></mrow></msub><mtext> = G</mtext></mrow></math><img file="EP1548618A1_D0017.tif" /></maths>

Moreover, four numerical sets g^ij, i, j = 1, 2 are defined as follows.<maths id="math0018" num="(equation 18)"><math display="block"><mrow><msup><mrow><mtext>g</mtext></mrow><mrow><mtext>11</mtext></mrow></msup><mtext> = </mtext><mfrac><mrow><mtext>G</mtext></mrow><mrow><msup><mrow><mtext>EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><msup><mrow><mtext>, g</mtext></mrow><mrow><mtext>12</mtext></mrow></msup><msup><mrow><mtext> = g</mtext></mrow><mrow><mtext>21</mtext></mrow></msup><mtext> =-</mtext><mfrac><mrow><mtext>F</mtext></mrow><mrow><msup><mrow><mtext>EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><msup><mrow><mtext>, g</mtext></mrow><mrow><mtext>22</mtext></mrow></msup><mtext> = </mtext><mfrac><mrow><mtext>E</mtext></mrow><mrow><msup><mrow><mtext>EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP1548618A1_D0018.tif" /></maths>

Furthermore, for the respective points, second fundamental quantities L, M and N are defined as follows:<maths id="math0019" num=""><img file="EP1548618A1_D0019.tif" /></maths>

Then, a second fundamental tensor (h_ij, i, j = 1, 2) is defined as follows:<maths id="math0020" num="(equation 20)"><math display="block"><mrow><msub><mrow><mtext>h</mtext></mrow><mrow><mtext>11</mtext></mrow></msub><msub><mrow><mtext> = L, h</mtext></mrow><mrow><mtext>12</mtext></mrow></msub><msub><mrow><mtext> = h</mtext></mrow><mrow><mtext>21</mtext></mrow></msub><msub><mrow><mtext> = M, h</mtext></mrow><mrow><mtext>22</mtext></mrow></msub><mtext> = N</mtext></mrow></math><img file="EP1548618A1_D0020.tif" /></maths>

Now, if the dynamic frame {e1, e2, n} is differentiated by the curved surface coordinates (u1, u2), structural equations of a curved surface shown by the following two equations (Gaussian equation of equation 21 and Weingarten's equation of equation 22) are obtained:<maths id="math0021" num=""><img file="EP1548618A1_D0021.tif" /></maths><maths id="math0022" num="(equation 22)"><math display="block"><mrow><mfrac><mrow><mtext>∂n</mtext></mrow><mrow><msup><mrow><mtext>∂u</mtext></mrow><mrow><mtext>i</mtext></mrow></msup></mrow></mfrac><msup><mrow><mtext> = -g</mtext></mrow><mrow><mtext>jk</mtext></mrow></msup><msub><mrow><mtext>h</mtext></mrow><mrow><mtext>ij</mtext></mrow></msub><msub><mrow><mtext>e</mtext></mrow><mrow><mtext>k</mtext></mrow></msub></mrow></math><img file="EP1548618A1_D0022.tif" /></maths><maths id="math0023" num=""><img file="EP1548618A1_D0023.tif" /></maths> where equation 23 exhibits the Christoffel symbol.

The integrable condition of these structural equations 21 and 22 is shown by the following two equations (the Gaussian equation of equation 24 and the Mainardi-Codazzi's equation of equation 25) :<maths id="math0024" num="(equation 24)"><math display="block"><mrow><msubsup><mrow><mtext>R</mtext></mrow><mrow><mtext>jkl</mtext></mrow><mrow><mtext>i</mtext></mrow></msubsup><msup><mrow><mtext> = g</mtext></mrow><mrow><mtext>im</mtext></mrow></msup><msub><mrow><mtext> (h</mtext></mrow><mrow><mtext>ik</mtext></mrow></msub><msub><mrow><mtext>h</mtext></mrow><mrow><mtext>lm</mtext></mrow></msub><msub><mrow><mtext> - h</mtext></mrow><mrow><mtext>jl</mtext></mrow></msub><msub><mrow><mtext>h</mtext></mrow><mrow><mtext>km</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP1548618A1_D0024.tif" /></maths><maths id="math0025" num=""><img file="EP1548618A1_D0025.tif" /></maths><maths id="math0026" num=""><img file="EP1548618A1_D0026.tif" /></maths> where equation 26 exhibits the Riemann-Christoffel curvature tensor.

In the case where the first fundamental tensor (g_ij, i, j = 1, 2) and the second fundamental tensor (h_ij, i, j = 1, 2) are applied as the function of the curved surface coordinate (u1, u2) and they satisfy the Gaussian equation and the Mainardi-Codazzi's equation described above, the shape of the curved surface having such g_ij and h_ij is uniquely determined (refer to Bonnet's fundamental theory). Therefore, the respective meshes become C2 continuous.

The CPU performs these arithmetic processes and calculates the integrable condition described above (step S2).

Next, the CPU executes a line of curvature analyzing process, a feature line analyzing process, and a curvature/girth length converting process (step S3). Firstly, the principal curvatures κ₁ and κ₂ in the mesh are calculated based on the coefficients of the first fundamental form E, F and G, and the coefficients of the second fundamental form L, M and N by the line of curvature analyzing process (step S4). That is to say, firstly the extremum of the curvature κ₁ described above is calculated. The shape of the normal section plane, which is the line of intersection of the normal plane and the curved surface, changes together with the tangential direction, and accompanied by this, the normal curvature also changes. This shape returns to the initial condition when the normal plane is half rotated. Now, assuming that;<maths id="math0027" num="(equation 27)"><math display="block"><mrow><mtext>γ = </mtext><mfrac><mrow><mtext>dv</mtext></mrow><mrow><mtext>du</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0027.tif" /></maths> and rewriting κ as the function κ(γ) of γ, gives:<maths id="math0028" num="(equation 28)"><math display="block"><mrow><msup><mrow><mtext>{L - κ(γ) · E} + 2{M - κ(γ) · F}γ + {N - κ(γ) · G}γ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> = 0</mtext></mrow></math><img file="EP1548618A1_D0028.tif" /></maths> From this quadratic equation of γ, for dκ(γ)/dy = 0, κ(γ) becomes the extremum. Then, if equation 15 is differentiated in this condition for the extremum, and κ and γ are rewritten as (<img file="EP1548618A1_D0029.tif" />) and (<img file="EP1548618A1_D0030.tif" />),<maths id="math0029" num=""><img file="EP1548618A1_D0031.tif" /></maths> is obtained. Then, if it is substituted in equation 16,<maths id="math0030" num=""><img file="EP1548618A1_D0032.tif" /></maths> is obtained. The following relations are obtained from these equations:<maths id="math0031" num=""><img file="EP1548618A1_D0033.tif" /></maths><maths id="math0032" num=""><img file="EP1548618A1_D0034.tif" /></maths>

If equation 18 is modified,<maths id="math0033" num=""><img file="EP1548618A1_D0035.tif" /></maths> is obtained. The coefficient <img file="EP1548618A1_D0036.tif" /> of is positive from equation 7. Assuming that the roots are κ₁ and κ₂, the value becomes the principal curvature as shown in FIG. 5.

Next, the Gaussian curvature or the mean curvature are calculated based on the principal curvature (step S5). That is, from the relation of the roots and the coefficients of the quadratic equation,<maths id="math0034" num="(equation 34)"><math display="block"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>m</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><msub><mrow><mtext>(κ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>+κ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>) = </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mfrac><mrow><mtext>(EN + LG - 2MF)</mtext></mrow><mrow><msup><mrow><mtext>(EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>)</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0037.tif" /></maths><maths id="math0035" num="(equation 35)"><math display="block"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>g</mtext></mrow></msub><msub><mrow><mtext> = κ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>κ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><msup><mrow><mtext>(LN - M</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>)</mtext></mrow><mrow><msup><mrow><mtext>(EG - F</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>)</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0038.tif" /></maths> is expressed. Here, K_m is the mean curvature and K_g is the Gaussian curvature. When K_g = 0, this is the case where the curved surface becomes the developable surface as shown in FIG. 6, and the line of curvature on the curved surface becomes a straight line. In the present embodiment, the point where this Gaussian curvature becomes zero is assumed to be the reference point of transformation described later.

As an appropriate point for the reference point of transformation other than this point, for example, lines of curvature, borderlines (ridgelines), isoclinic orthogonal lines shown in FIG. 7, principal curvature extremum lines shown in FIG. 8, isoclinic extremum lines shown in FIG. 9, or umbilicus points may be selected. These are points or lines which become a reference point or a reference line of transformation defined by changing patterns of one or more feature quantities among the principal curvature, the principal direction, the Gaussian curvature, the mean curvature and the line of curvature, which are feature quantities showing the feature of the curved surface. It is possible to calculate these based on the coefficients of the first fundamental form and the coefficients of the second fundamental form.

Moreover, the line of curvature showing the principal direction of the mesh is calculated based on the principal curvature. That is, eliminating <img file="EP1548618A1_D0039.tif" /> from equation 19 gives:<maths id="math0036" num=""><img file="EP1548618A1_D0040.tif" /></maths>    or<maths id="math0037" num="(equation 37)"><math display="block"><mrow><msup><mrow><mtext>(MG - NF)dv</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msup><mrow><mtext> + (GL - NE)dudv + (FL - ME)du</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> = 0</mtext></mrow></math><img file="EP1548618A1_D0041.tif" /></maths> Both these two equations are equations of the line of curvature and the quadratic equation so that γ₁ and γ₂ have the following relations:<maths id="math0038" num="(equation 38)"><math display="block"><mrow><msub><mrow><mtext>γ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> + γ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>-(GL - NE)</mtext></mrow><mrow><mtext>(MG -NF)</mtext></mrow></mfrac><msub><mrow><mtext>, γ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>γ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>=</mtext><mfrac><mrow><mtext>(FL - ME)</mtext></mrow><mrow><mtext>(MG - NF)</mtext></mrow></mfrac></mrow></math><img file="EP1548618A1_D0042.tif" /></maths>

At a point on the curved surface, the curvature becomes the extremum in the direction determined by γ₁ and γ₂. The tangent vector on the curved surface is (Sudu + Svdv) and the inner product of the two tangent vectors corresponding to γ₁ and γ₂ becomes:<maths id="math0039" num="(equation 39)"><math display="block"><mrow><msub><mrow><mtext>(ds)</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> · (ds)</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext> = {(s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> + s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext>γ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>) · (s</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> + s</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext>γ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>)} (du)</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> (dv)</mtext></mrow><mrow><mtext>2</mtext></mrow></msub></mrow></math><img file="EP1548618A1_D0043.tif" /></maths> If inside the bracket {} is converted then:<maths id="math0040" num="(equation 40)"><math display="block"><mrow><mtext>{E(MG - NF) - F(GL - NE) + G(FL - ME)}/(MG - NF)</mtext></mrow></math><img file="EP1548618A1_D0044.tif" /></maths> becomes zero. That is, it is found that the two tangential directions of the normal section plane of the principal curvature become orthogonal. This direction is called the principal direction. In the case where this direction and the tangent line on the curved surface are matched, this becomes the lines of curvature shown in FIG. 10.

From the above, the calculation process of the line of curvature showing the principal direction of the mesh is performed.

Next, the curvature/girth length converting process is performed (step S6). That is, the CPU calculates the girth length based on the curvature which is calculated based on the coefficients of the first fundamental form E, F, and G and the coefficients of the second fundamental form L, M, and N. Along the line of curvature calculated by the line of curvature calculating process described above, the radius of curvature is calculated from the curvature (1/r), and the girth length of the line of curvature is expanded and contracted in each calculation interval.

From the above, the analyzing process is performed.

Next, after the point sequence information, the coefficients of the first fundamental form, and the coefficients of the second fundamental form, which were generated and extracted at step S1 and step S2, have been collected (Yes in step S7), the CPU performs the curved surface data transferring process (step S9). On the other hand, in the case where such information has not been completed, the database evaluating process is performed (No in step S7). That is, the shape reproduced based on the principal direction, the reference position (point, line, or the like), the transformed amount, that were calculated at steps S4 to S6, and the shape reproduced based on the point sequence information, the coefficients of the first fundamental form, and the coefficients of the second fundamental form are compared, and in the case where they match (Yes in step S8), the curved surface data transferring process is performed (step S9). Moreover, in the case where they do not match (No in step S8), an accuracy improving process by approximation and interpolation is performed. That is, the initial curved surface is approximated and interpolated so that it becomes doubly differentiable, and the processes described above are repeatedly performed again from step SI. Then, at the stage where the comparative evaluation in step S8 is matched, the flow shifts to the curved surface data transferring process.

The curved surface data is transferred to the converting program 2 or the reproducing program 3 shown in FIG. 1. If the CPU receives a convert command, it executes the converting program 2. That is, firstly assuming that a point selected as a feature point or a feature line, where the Gaussian curvature becomes zero, is the transformation reference, the line of curvature is expansion and contraction transformed by the girth length in the line of curvature direction so that the mesh or the curved surface is reproduced. Then, a plurality of point sequences on the curved surface are extracted from the reproduced mesh or curved surface, and the point sequences are converted according to a graphical representation algorithm in an other computer aided design system. The converted graphics data is reproduced by the other computer aided design application 22 and then output to the graphic display processing section 11. The graphic display processing section 11 performs graphic display processing of the data output from the computer aided design application 22 and outputs this to the display section 12. The display section 12 receives the input of the display data and displays this.

Moreover, if the CPU receives a reproduction command, it executes the reproducing program 3. The reproducing program makes the CPU to execute the processes in the converting program except for the converting process. That is, assuming that a point where the Gaussian curvature becomes zero, is the transformation reference, the line of curvature is expansion and contraction transformed by the girth length in the line of curvature direction so that the mesh or the curved surface is reproduced. Then, the reproduced graphics data is output to the graphic display processing section 11, and after display processing, it is displayed in the display section 12.

As described above, according to the computer aided design system of the present embodiment, an effect can be obtained where a free-form surface can be analyzed, converted and reproduced while retaining the continuity of the C2 continuous. Therefore, an effect can be obtained where the utility value of a computer aided design model can be greatly increased, and the efficiency of the design and production processes can be improved.

In the computer aided design system of the present embodiment, the description is for a series of processes related to free-form surface analysis, conversion, and reproduction in a computer aided design model. However, the computer aided design system of the present invention is not limited to this, and is applicable to a computer graphics system, or a system and a program which performs graphical representation using a computer.

Moreover, in the computer aided design system of the present embodiment, as shown in FIG. 2, as a suitable example, a curved surface is divided into meshes, and then standardized by the fundamental vectors Su and Sv, so that the free-form surface analysis, conversion, and reproduction are performed by a u, v parameter form which uses point sequence information (u, v). However, the computer aided design system of the present invention is not limited to this, and coordinate values for (x, y, z) coordinate parameters may be used.

The computer aided design system described above contains a computer system inside. Moreover steps of a series of processes related to the aforementioned free-form surface analysis, conversion, and reproduction are stored in a computer readable recording media in program format. A computer reads out and executes this program, to thereby perform the above processes. Here, the computer readable recording media is for example a magnetic disk, magnetooptical disk, CD-ROM, DVD-ROM, semiconductor memory or the like. Moreover, the arrangement may be such that this computer program is delivered to a computer by a communication line, and the computer which receives this delivery, executes the program.