WO2018131304A1 - 曲面生成装置および曲面生成用プログラム - Google Patents

曲面生成装置および曲面生成用プログラム Download PDF

Info

Publication number
WO2018131304A1
WO2018131304A1 PCT/JP2017/042636 JP2017042636W WO2018131304A1 WO 2018131304 A1 WO2018131304 A1 WO 2018131304A1 JP 2017042636 W JP2017042636 W JP 2017042636W WO 2018131304 A1 WO2018131304 A1 WO 2018131304A1
Authority
WO
WIPO (PCT)
Prior art keywords
cross
curved surface
sectional curve
curve
sectional
Prior art date
Application number
PCT/JP2017/042636
Other languages
English (en)
French (fr)
Inventor
庄一 土江
卓 前川
Original Assignee
日本ユニシス株式会社
国立大学法人横浜国立大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 日本ユニシス株式会社, 国立大学法人横浜国立大学 filed Critical 日本ユニシス株式会社
Priority to EP17891292.9A priority Critical patent/EP3550456A4/en
Priority to US16/477,061 priority patent/US11145120B2/en
Publication of WO2018131304A1 publication Critical patent/WO2018131304A1/ja
Priority to US17/476,163 priority patent/US20220005274A1/en

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30—Polynomial surface description
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/12—Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]

Definitions

  • the present invention relates to a curved surface generation apparatus and a curved surface generation program, and is particularly suitable for use in a curved surface generation apparatus for generating a curved surface having a smooth curvature change required in the field of design design.
  • product shape data is generated by CAD (Computer Aided Design).
  • CAD Computer Aided Design
  • a method of reconstructing a curved surface model of an object from measurement data of the object as a model by CAD is widely used.
  • CAD data is created from measurement data of a clay model that is modeled by a clay modeler using a clay spatula (curve ruler).
  • a curved surface model by CAD for a clay model is reconstructed by regarding the curve ruler as a cross-sectional curve and sweeping the cross-sectional curve along a certain flow. If both the cross-sectional curve and the flow of the cross-sectional curve are correctly expressed, a high-quality curved surface can be generated. However, in particular, it is extremely difficult to correctly express a curve indicating the flow of a cross-sectional curve from measurement data, and there is a problem that a high-quality curved surface cannot be easily generated.
  • a curved surface reproduction device that reproduces curved surface data based on a curvature line indicating the main direction of the mesh (see, for example, Patent Document 4).
  • a curved surface in real space is mapped to a parameter space, and a primary basic quantity and a secondary basic quantity are calculated.
  • the real space coordinates of are calculated.
  • sweep-based methods have the merit that control points can be intentionally arranged regularly.
  • a method for calculating an appropriate guide line has not been established. Therefore, the CAD modeler has to repeatedly adjust the guide line by trial and error until the quality of the generated curved surface can be ensured. It took a lot of effort.
  • the present invention has been made to solve such a problem, and an object thereof is to efficiently generate a high-quality curved surface without repeating many trials and errors.
  • a cross-sectional curve fitted to a part of model shape data is set, and the cross-sectional curve is translated and rotated so as to satisfy a predetermined condition.
  • a curved surface is generated by interpolation from the boundary of the curved surface defined by the trajectory of the moved sectional curve.
  • the predetermined condition is a condition that a locus formed by one point (for example, a middle point) on the cross-sectional curve becomes a curved line of a curved surface to be generated.
  • the cross section curve corresponding to the guide line can be obtained by performing the processing of the present invention configured as described above.
  • a surface corresponding to a sweep of the set cross-sectional curve is generated along the locus of the upper point.
  • the trajectory corresponding to the guide line is given as a curvature line representing the flow of the curved surface, and a high-quality curved surface that uses the curvature line representing the flow of the curved surface as a guide line for the sweep method is subjected to many trials and errors. Can be generated efficiently without any problems.
  • FIG. 1 is a block diagram illustrating a functional configuration example of the curved surface generation apparatus according to the first embodiment.
  • the curved surface generation apparatus according to the first embodiment is configured to include a shape data input unit 1, a cross section curve setting unit 2, a cross section curve moving unit 3, and a curved surface generation unit 4 as its functional configuration.
  • the cross-sectional curve moving unit 3 includes a first movement processing unit 31 and a second movement processing unit 32.
  • the above functional blocks 1 to 4 can be configured by any of hardware, DSP (Digital Signal Processor), and software.
  • DSP Digital Signal Processor
  • each of the functional blocks 1 to 4 is actually configured with a computer CPU, RAM, ROM, and the like, and is stored in a recording medium such as RAM, ROM, hard disk, or semiconductor memory. Is realized by operating.
  • the shape data input unit 1 inputs model shape data.
  • the shape data input unit 1 inputs measurement data (point cloud data or polygon data) of a shape such as a clay model.
  • the cross-section curve setting unit 2 sets a cross-section curve fitted to a part of the shape data input by the shape data input unit 1.
  • the cross-sectional curve set here is created by imitating a curve ruler used when a design designer (Clay Modeler) creates a clay model of an automobile.
  • FIG. 2 is a diagram illustrating an outline of processing of the curved surface generation apparatus according to the first embodiment. In FIG. 2, an example of a cross-sectional curve C (s) set by imitating the curve ruler 200 by the cross-sectional curve setting unit 2 is shown.
  • the cross section curve C (s) fitted to a part of the shape data by the cross section curve setting unit 2 is moved as described later, and is defined from the locus of the cross section curve C (s) formed by the movement.
  • a target curved surface (hereinafter referred to as a target curved surface) S (u, v) is generated by interpolation from the boundary line.
  • the target curved surface S (u, v) is represented by a B-spline curved surface.
  • the shape of the target curved surface S (u, v) in the u direction is completely characterized by the cross-sectional curve C (s) set by the cross-sectional curve setting unit 2.
  • two boundary lines Cs and Ce with respect to the v direction of the target curved surface S (u, v) are defined as trajectories when both end points of the cross-sectional curve C (s) are moved.
  • the section curve moving unit 3 performs a process of moving the section curve C (s) initially set by the section curve setting unit 2 so as to satisfy a predetermined condition.
  • the movement of the cross-sectional curve C (s) is performed by the first movement processing unit 31 and the second movement processing unit 32.
  • the first movement processing unit 31 performs only parallel movement
  • the second movement processing unit 32 performs both parallel movement and rotational movement. Details of these processes will be described below.
  • the first movement processing unit 31 translates the sectional curve C (s) initially set by the sectional curve setting unit 2 so that an error from the shape data input by the shape data input unit 1 is minimized. Then, a tangent vector indicating the direction of the locus in which one point on the cross-sectional curve C (s) is translated is obtained at predetermined intervals.
  • One point on the cross-sectional curve C (s) is, for example, the midpoint of the cross-sectional curve C (s), but is not limited thereto.
  • the predetermined interval may be a fixed value set in advance in the curved surface generation apparatus, or may be a variable value that can be specified by the user.
  • FIG. 3 is a diagram for explaining the processing contents of the first movement processing unit 31.
  • the first movement processing unit 31 translates the cross-sectional curve C (s) set by the cross-sectional curve setting unit 2 so as to fit a part of the shape data in the following procedure.
  • the initially set sectional curve C (s) is defined as C 0 (s), and at a predetermined interval with respect to the direction of the normal vector T 0 ⁇ (not shown in FIG. 3) of the plane P 0 including this.
  • the cross-sectional curve C 0 (s) is translated, and the obtained cross-sectional curve C (s) after the movement is defined as C 1 (s) ⁇ (not shown in FIG. 3). Note that if the sectional curve C 0 (s) is simply translated in the direction of the normal vector T 0 ⁇ , the sectional curve C 1 (s) ⁇ after the movement may be separated from the shape data and cause an error. There is.
  • the first movement processing unit 31 translates the section curve C 0 (s) in the direction of the normal vector T 0 ⁇ , and then includes a plane including the moved section curve C 1 (s) ⁇ . by translating the ⁇ sectional curve C 1 (s) in the P 1 by a C 1 (s), so that the error of the cross-sectional curve C 1 (s) and the shape data is minimized.
  • the sectional curve C 1 (s) is given at a predetermined interval with respect to the direction of the normal vector T 1 ⁇ (not shown in FIG. 3) of the plane P 1 including the moved sectional curve C 1 (s).
  • the cross-sectional curve C (s) obtained after the parallel movement is defined as C 2 (s) ⁇ (not shown in FIG. 3).
  • the first moving unit 31 in the plane P 2 containing the sectional curve C 2 (s) ⁇ after movement, the cross-sectional curve C 2 (s) ⁇ and is moved parallel C 2 (s)
  • the error between the cross-sectional curve C 2 (s) and the shape data is minimized.
  • C i + 1 (s) ⁇ in the plane P i + 1 the error between the sectional curve C i + 1 (s) and the shape data is minimized.
  • each sectional curve C i (s) is fitted to the shape data as much as possible.
  • the tangent vector T i obtained for each predetermined interval corresponds to a vector indicating the direction of the locus in which one point on the cross-sectional curve C i (s) is translated.
  • the curve C i (s) is further moved by translation and rotation, respectively.
  • the second moving unit 32 updates the tangent vector T i when moving parallel sectional curve C i (s), performs a rotational movement on the basis of the tangent vector T i after updating.
  • the first condition corresponds to the first condition in the claims, and is a point on the cross-sectional curve C i (s) ′ for each predetermined interval moved by the second movement processing unit 32.
  • pseudo guide line Sc In order to make the trajectory formed by (corresponding to the guide line when the cross-sectional curve C (s) is swept, referred to below as pseudo guide line Sc (see FIG. 2)) close to the curvature line with small twist. This is the condition.
  • the second condition corresponds to the third condition in the claims, and is input by the target curved surface S (u, v) generated by the curved surface generation unit 4 and the shape data input unit 1. This is a condition for minimizing an error from a curved surface represented by shape data.
  • the second movement processing unit 32 uses the point C i (s * ) on the cross-sectional curve C i (s) as the rotation center so that the tangent vector T i-1 matches the next tangent vector T i.
  • the parameter s * is a value that gives one point (for example, a midpoint) on the cross-sectional curve C i (s).
  • the rotation matrix R shown in (Equation 1) can be considered as approximating a rotation minimizing frame (RMF) in which there is no rotation in the direction of the tangent vector T i .
  • RMF rotation minimizing frame
  • FIG. 4 is a schematic diagram illustrating an example of the movement (parallel movement and rotation movement) of the cross-sectional curve C i (s) executed by the second movement processing unit 32.
  • FIG. 4A shows a result of the cross-sectional curve C i (s) being translated by the first movement processing unit 31 (the same state as FIG. 3B).
  • FIG. 4B shows a result obtained by further translating and rotating the sectional curve C i (s) shown in FIG.
  • a cross-sectional curve C i (s) ′ shown in FIG. 4B shows the result of movement so as to satisfy the two conditions described above (the first condition and the third condition in the claims). is there.
  • the pseudo guide line Sc means an intentional flow on the target curved surface S (u, v) to be generated, and the flow is specified as a curvature line.
  • the curvature line indicates the flow of the curved surface in the main direction with respect to the minimum main curvature in the target curved surface S (u, v).
  • P2 The pseudo guide wire Sc is less twisted.
  • the condition P1 is that the pseudo guide line Sc, which is a trajectory formed by connecting one point C i (s * ) ′ on the cross-sectional curve C i (s) ′ for each predetermined interval moved based on the parallel movement and the rotational movement, is This is a condition that the curvature line of the target curved surface S (u, v) generated by the curved surface generation unit 4 is obtained.
  • Condition P2 is that the angle from the normal vector of the target curved surface S (u, v) generated by the curved surface generation unit 4 to the main normal vector of the pseudo guide line Sc is constant on the locus of the pseudo guide line Sc. This is the condition.
  • the curved surface generation unit 4 generates a curved surface by interpolation from the boundary line defined by the trajectory of the sectional curve C (s) ′ moved by the sectional curve moving unit 3. That is, in FIG. 4B, the curved surface generation unit 4 has two cross-sectional curves C 0 (s) ′ and C d (s) ′ with respect to the u direction at the start and end points of movement, and two boundary lines with respect to the v direction.
  • the objective curved surface S (u, v) is generated by, for example, a known Coons patch by interpolation from these four sides.
  • the second movement processing unit 32 performs a process of bringing the pseudo guide line Sc closer to the curvature line in order to achieve the condition P1.
  • the curvature line is given by the following (Equation 2).
  • (EM ⁇ FL) du 2 + (EN ⁇ GL) dudv + (FN ⁇ GM) dv 2 0
  • E, F, and G are coefficients of the first basic form of the curved surface (first basic quantity), respectively
  • L, M, and N are coefficients of the second basic form of the curved face (second basic quantity), respectively. .
  • V N represents a normal vector of the target curved surface S (u, v).
  • the geodetic torsion ⁇ g of a certain curve C on the curved surface is given by the following (formula 4).
  • ⁇ g (d ⁇ / ds) ⁇ (Formula 4)
  • ⁇ is the normal torsion of the curve C
  • ⁇ is the angle from the normal vector V N of the curved surface to the main normal vector of the curve C.
  • the pseudo guide line Sc is a curvature line with respect to the target curved surface S (u, v), and the torsion ⁇ of the pseudo guide line Sc is small. Evaluation is performed by the function ⁇ shown in (Expression 5).
  • w 1 and w 2 are weighting factors.
  • N 2 is the number of sampling points extracted from the pseudo guide line Sc. This sampling point is preferably extracted uniformly from the pseudo guide line Sc.
  • the third condition is that the error between the target curved surface S (u, v) generated by the curved surface generating unit 4 and the curved surface represented by the shape data input by the shape data input unit 1 is minimized. is there.
  • Equation 6 the error between the target curved surface S (u, v) and the sampling point Q i on the shape data is small.
  • N is the number of sampling points Q i extracted from the shape data.
  • the sampling points Q i are preferably extracted from the entire shape data without deviation, and are appropriately extracted so as to cover the entire curved surface.
  • the second movement processing unit 32 considers the first condition and the third condition described above, and the sectional curve C (s) so as to minimize the objective function J defined by the following (Expression 7). Move.
  • the target curved surface S (u
  • a trajectory (pseudo guide line Sc) defined by the trajectory of one point on the cross-sectional curve C (s) is generated based on the model shape. It can be expressed as a curvature line representing the flow of the curved surface in the main direction with respect to the minimum main curvature of the target curved surface S (u, v) to be obtained.
  • the curve C The trajectory of the point C (s * ) due to the movement of (s) becomes the curvature line of the target curved surface S (u, v) to be generated and has the property that the twist is small. ) Is generated. That is, it is possible to efficiently generate a good target curved surface S (u, v) in which the pseudo guide line Sc becomes a curvature line with a small torsion without repeating many trials and errors.
  • the condition P1 in which the pseudo guide line Sc becomes a curvature line of the target curved surface S (u, v) by setting the weighting coefficient w 2 in (Equation 5) to zero is the first condition.
  • the first condition may be used. Even if only the condition P1 is used, the section curve C (s) is moved by the rotation approximating the RMF, so that it is possible to avoid the torsion of the pseudo guide line Sc from being increased explicitly.
  • the target curved surface S (u, v) is generated from the measurement data such as the car clay model.
  • the shape data used as the input data is the measurement data of the clay model. It is not limited to.
  • the target curved surface S (u, v) may be generated from measurement data of a commercial vehicle.
  • the shape data is not limited to the measurement data, and may be CAD data or its discrete approximation.
  • the target model is not limited to a car, and can be any object.
  • the generated target curved surface S (u, v) and shape data such as a clay model are used.
  • the position of each cross-sectional curve C i (s) is adjusted so as to satisfy the third condition for minimizing the error with respect to the first term representing the third condition in Equation (7).
  • a second term representing the first condition may be set so that the weights of both conditions can be adjusted.
  • CAD data representing a concept of a design to be newly created is input as input data from the shape data input unit 1 instead of measurement data for the actual model, and the target curved surface S (u, v) is generated. In some cases, it is possible to reduce the weight of the third condition for processing.
  • FIG. 5 is a block diagram illustrating a functional configuration example of the curved surface generation apparatus according to the second embodiment.
  • the curved surface generation apparatus according to the second embodiment includes a second movement processing unit 32 ′ instead of the second movement processing unit 32.
  • the sectional curve C (s) when the sectional curve C (s) is moved in order to improve the quality of the target curved surface S (u, v) reconstructed from the shape data (measurement data such as a clay model).
  • the first condition is satisfied so that the trajectory of the point C i (s * ) ′ on each possible section curve C i (s) ′ becomes a curvature line with a small twist on the target curved surface S (u, v).
  • An example of adjusting the position of each cross-sectional curve C i (s) ′ has been described.
  • each sectional curve C i (s) is adjusted by imposing the second condition on the boundary lines Cs and Ce in the v direction of the target curved surface S (u, v).
  • the second condition is a condition that the curves of the boundary lines Cs ′ and Ce ′ have curvature monotonicity.
  • the curvature monotony required for design design is not always satisfied. Therefore, the optimum tangent vector T i is calculated by finely adjusting each cross-sectional curve C i (s) in the plane P i so that both the boundary lines Cs and Ce have curvature monotonicity.
  • the second movement processing unit 32 ′ sets the cross section curve C (s) so as to minimize the objective function J defined as in (Expression 10) instead of (Expression 7). Move.
  • ⁇ 1 and ⁇ 2 are weighting factors.
  • the locus of one point C i (s * ) ′ (for example, the middle point) on the cross-sectional curve C i (s) ′ is a desired curve as the pseudo guide line Sc.
  • the boundary lines Cs ′ and Ce ′ formed by the trajectories of the end points of the cross-sectional curve C i (s) ′ are curves having a monotone curvature.
  • the pseudo guide line Sc also has curvature monotonicity. For this reason, the quality of the target curved surface S (u, v) reconstructed from the shape data can be further improved.
  • Equation 5 by setting the weighting coefficient w 2 zero (Equation 5), may be used only the condition P1 as a first condition.
  • the shape data as input data is not limited to measurement data such as a clay model, and may be CAD data or discrete approximation thereof.
  • Equation 10 a weighting factor may be set for the first term representing the third condition.
  • the target curved surface S (u, v) generated by the curved surface generating unit 4 is expressed by a B-spline curved surface.
  • the present invention is not limited to this.
  • the target curved surface S (u, v) may be expressed by a Bezier curved surface.
  • the cross-sectional curve C (s) may be set according to the shape of the target curved surface S (u, v) to be generated.
  • the present invention is not limited to this. That is, the cross-sectional curve C (s) may be moved as a one-step process by rotating the cross-sectional curve C (s) while being translated at predetermined intervals.
  • the rotational movement of the cross-sectional curve C (s) is defined by (Expression 1), but is not limited to (Expression 1). That is, there are several other alternatives that approximate RMF. Further, (Equation 1) can cope with rotation around an arbitrary axis by the setting method of the rotation axis A (i) .
  • the objective curved surface S (u, v) is generated using interpolation from the boundary line by the Coons patch.
  • the present invention is not limited to this.
  • the locus of one point C (s *) (for example, the middle point) on the cross section curve C (s) is approximated by a B-spline curve or the like, and the cross section curve C (s) is swept along the approximate curve.
  • all the above arguments are applicable.
  • the sectional curve C i (s) is translated in the normal vector direction of the plane P i and then the parallel movement is performed.
  • the translation may be performed only in the normal vector direction.
  • the tangent vector T i in which the error between the cross-sectional curve C i + 1 (s) and the shape data is minimized is not generated as the initial vector, the first movement processing unit 32 performs the optimization calculation to perform the first calculation.
  • a target curved surface S (u, v) satisfying the conditions (1) to (3) can be generated.
  • a tangent vector T i that minimizes the error between the cross-sectional curve C i + 1 (s) and the shape data is generated as an initial vector, so that the second There is an advantage that the processing time of the optimization calculation by the movement processing unit 32 can be shortened.
  • each of the first and second embodiments described above is merely an example of a specific example for carrying out the present invention, and the technical scope of the present invention should not be interpreted in a limited manner. It will not be. That is, the present invention can be implemented in various forms without departing from the gist or the main features thereof.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Mathematical Physics (AREA)
  • Computer Graphics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Human Computer Interaction (AREA)
  • Numerical Control (AREA)
  • Image Generation (AREA)
  • Processing Or Creating Images (AREA)

Abstract

モデルの形状データにフィッティングした断面曲線を設定する断面曲線設定部２と、所定の条件を満たすように断面曲線を移動させる断面曲線移動部３と、所定の条件を満たす断面曲線が動く軌跡によって定義された曲面を生成する曲面生成部４とを備え、上記所定の条件として、断面曲線上の一点（例えば、中点）の軌跡が、生成しようとする曲面の曲率線になるという条件を用いることにより、曲面の流れを表す曲率線をスイープ法のガイド線とするような高品質の曲面を、多くの試行錯誤を繰り返すことなく効率的に求めることができるようにする。

Description

曲面生成装置および曲面生成用プログラム
　本発明は、曲面生成装置および曲面生成用プログラムに関し、特に、意匠デザイン分野で求められる曲率変化の滑らかな曲面を生成するための曲面生成装置に用いて好適なものである。
　一般に、製品の設計を行う際には、ＣＡＤ（Computer Aided Design）によって製品の形状データを生成することが行われている。形状データを生成する１つの方法として、モデルとなる物体の測定データからその物体の曲面モデルをＣＡＤによって再構築する方法が広く利用されている。例えば、自動車の曲面形状の場合、クレイモデラが粘土ヘラ（カーブ定規）を使って造形したクレイモデルの測定データからＣＡＤデータを作成することが行われている。
　測定データからＣＡＤで曲面モデルを再構築する場合、(i)生成する曲面は指定されたトレランス以内で測定データを十分に近似できていること、(ii)生成する曲面はその曲率変化が滑らかであること、の２点が要求される。特に、(ii)については、意匠デザインにおける曲面の美的性質にとって重要であるのみならず、実際に消費者が商品の購入を検討する際のインパクトにも繋がる点で極めて重要である。
　自動車の曲面生成の場合、カーブ定規を断面曲線とみなし、当該断面曲線をある流れに沿ってスイープさせることにより、クレイモデルに対するＣＡＤによる曲面モデルを再構築することが行われている。断面曲線と、当該断面曲線の流れとの双方が正しく表現されていれば、高品質の曲面を生成することが可能である。しかしながら、特に、断面曲線の流れを示す曲線を測定データから正しく表現することが極めて困難であり、高品質の曲面を生成することが容易にはできないという問題があった。
　なお、従来、ＣＡＤ等において、ガイドとなる曲線に沿って断面曲線をスイープさせることによって曲面を生成する技術がいくつか知られている（例えば、特許文献１～３参照）。ガイドとなる曲線は文献により異なる名称（制御曲線、軌道曲線、特徴線など）が与えられているが、以下、それらを総称して、ガイド線とよぶ。
　特許文献１に記載の曲面生成装置では、座標入力装置からの指示により、ガイド線に沿って入力曲線をスイープしながら、スイープ位置での指定回転量・倍率で入力曲線を回転・拡大することにより、自由曲面を生成する。特許文献２に記載の３次元曲面形状生成装置では、ガイド線および断面曲線に加え、補間曲線と制御ベクトルとを用いてスイープ処理をすることにより、軌道の曲率変化が激しい場合でも正確なスイープ曲面を生成可能としている。特許文献３に記載の曲面作成装置では、ガイド線に沿って断面曲線が変形移動した曲面で、ガイド線上で支持曲面と滑らかにつながり、かつ支持曲線上を通過する曲面を作成する。
　また、メッシュの主方向を示す曲率線に基づいて曲面データを再生するようになされた曲面再生装置も知られている（例えば、特許文献４参照）。特許文献４に記載の曲面再生装置では、実空間上の曲面をパラメータ空間に写像し、一次基本量と二次基本量とを算出する。次に、曲面の曲率線および曲率線が通過する実空間座標を算出し、パラメータ空間上の逆写像目標点を通過する線とパラメータ空間上に写像された曲率線とのパラメータ空間上での交点の実空間座標を算出する。そして、算出した交点の実空間座標を通過する実空間上の曲線を算出し、これに基づいて曲線上の逆写像目標点の実空間座標値を算出して、実空間上で曲面データを再生する。
特許第２９３８９０９号公報 特開平５－２５０４４４号公報 特開２０００－１１２１０号公報 特開２００５－１４９２４５号公報
　一般にスイープに基づく方法では、意図的に規則正しく制御点を配置できるというメリットがある。しかし、断面曲線のスイープに必要なガイド線を測定データから正しく特定することが難しいという問題がある。スイープ法によって生成される曲面の品質は、特定するガイド線の品質に大きく依存するため、ガイド線を適切に表現することが極めて重要となる。しかし、適切なガイド線を算出する方法が確立されておらず、そのためＣＡＤモデラは、生成した曲面の品質が確保できるまで、試行錯誤によりガイド線の調整を繰り返し行わなければならず、多大な時間と労力がかかっていた。
　本発明は、このような問題を解決するために成されたものであり、高品質の曲面を、多くの試行錯誤を繰り返すことなく効率的に生成できるようにすることを目的とする。
　上記した課題を解決するために、本発明では、モデルの形状データの一部に対してフィッティングされた断面曲線を設定し、所定の条件を満たすように、断面曲線を平行移動および回転移動させ、当該移動された断面曲線の軌跡によって定義される曲面の境界からの内挿により曲面を生成する。ここで、上記所定の条件は、断面曲線上の一点（例えば、中点）により形成される軌跡が、生成しようとする曲面の曲率線になるという条件としている。
　モデルの形状データと断面曲線とが設定されていて、断面曲線のスイープに必要なガイド線が不明の場合、上記のように構成した本発明の処理を行えば、ガイド線に相当する、断面曲線上の一点の軌跡に沿って、設定された断面曲線をスイープしたものに相当する面が生成される。ガイド線に相当する軌跡は曲面の流れを表す曲率線として与えられることになり、曲面の流れを表す曲率線をスイープ法のガイド線とするような高品質の曲面を、多くの試行錯誤を繰り返すことなく効率的に生成することができる。
第１の実施形態による曲面生成装置の機能構成例を示すブロック図である。 第１の実施形態による曲面生成装置の処理の概要を示す図である。 第１の実施形態による第１の移動処理部の処理内容を説明するための図である。 第１の実施形態による第２の移動処理部によって実行される断面曲線の移動の一例を示す模式図である。 第２の実施形態による曲面生成装置の機能構成例を示すブロック図である。
（第１の実施形態）
　以下、本発明による第１の実施形態を図面に基づいて説明する。図１は、第１の実施形態による曲面生成装置の機能構成例を示すブロック図である。図１に示すように、第１の実施形態による曲面生成装置は、その機能構成として、形状データ入力部１、断面曲線設定部２、断面曲線移動部３および曲面生成部４を備えて構成されている。断面曲線移動部３は、第１の移動処理部３１および第２の移動処理部３２を備えている。
　上記各機能ブロック１～４は、ハードウェア、ＤＳＰ（Digital Signal Processor）、ソフトウェアの何れによっても構成することが可能である。例えばソフトウェアによって構成する場合、上記各機能ブロック１～４は、実際にはコンピュータのＣＰＵ、ＲＡＭ、ＲＯＭなどを備えて構成され、ＲＡＭやＲＯＭ、ハードディスクまたは半導体メモリ等の記録媒体に記憶されたプログラムが動作することによって実現される。
　形状データ入力部１は、モデルの形状データを入力する。第１の実施形態では、形状データ入力部１は、クレイモデルなどの形状の測定データ（点群データまたはポリゴンデータ）を入力する。
　断面曲線設定部２は、形状データ入力部１により入力された形状データの一部にフィッティングした断面曲線を設定する。ここで設定する断面曲線は、意匠デザイナ（クレイモデラ）が自動車のクレイモデルを作成する際に使用するカーブ定規を模して作成したものである。図２は、第１の実施形態による曲面生成装置の処理の概要を示す図である。この図２の中に、断面曲線設定部２によりカーブ定規２００を模して設定された断面曲線Ｃ(s)の一例が示されている。
　第１の実施形態では、この断面曲線設定部２により形状データの一部にフィッティングした断面曲線Ｃ(s)を後述するように移動させ、移動によってできる断面曲線Ｃ(s)の軌跡から定義される境界線からの内挿によって、目的の曲面（以下、目的曲面という）Ｓ(u, v)を生成する。
　第１の実施形態では、目的曲面Ｓ(u, v)をＢ－ｓｐｌｉｎｅ曲面で表現する。図２に示すように、目的曲面Ｓ(u, v)のｕ方向の形状は、断面曲線設定部２により設定された断面曲線Ｃ(s)によって完全に特徴付けられる。また、目的曲面Ｓ(u, v)のｖ方向に対する２つの境界線Ｃs，Ｃeは、断面曲線Ｃ(s)の両端点が移動した際の軌跡として定義する。
　断面曲線移動部３は、断面曲線設定部２により初期設定された断面曲線Ｃ(s)を、所定の条件を満たすように移動させる処理を行う。第１の実施形態では、この断面曲線Ｃ(s)の移動を、第１の移動処理部３１および第２の移動処理部３２によって行う。このときの移動は、第１の移動処理部３１では平行移動のみ、第２の移動処理部３２では平行移動および回転移動の双方を行う。これらの処理の詳細を以下に説明する。
　第１の移動処理部３１は、形状データ入力部１により入力された形状データとの誤差が最小となるように、断面曲線設定部２により初期設定された断面曲線Ｃ(s)を平行移動させ、断面曲線Ｃ(s)上の一点が平行移動した軌跡の方向を示す接線ベクトルを所定の間隔毎に求める。断面曲線Ｃ(s)上の一点は、一例としては断面曲線Ｃ(s)の中点であるが、これに限定されるものではない。所定の間隔は、曲面生成装置にあらかじめ設定された固定値としてもよいし、ユーザが指定可能な可変値としてもよい。
　図３は、第１の移動処理部３１の処理内容を説明するための図である。以下、この図３を用いて第１の移動処理部３１の処理例を具体的に説明する。すなわち、第１の移動処理部３１は、断面曲線設定部２により形状データの一部にフィッティングするように設定された断面曲線Ｃ(s)を、次の手順で平行移動させる。
　まず、初期設定された断面曲線Ｃ(s)をＣ0(s)とし、これを含む平面Ｐ0の法線ベクトルＴ0~（図３において図示せず）の方向に対し、所定の間隔で断面曲線Ｃ0(s)を平行移動し、得られた移動後の断面曲線Ｃ(s)をＣ1(s)~（図３において図示せず）とする。なお、法線ベクトルＴ0~の方向に断面曲線Ｃ0(s)を単純に平行移動させるだけだと、移動後の断面曲線Ｃ1(s)~は形状データから離れて誤差が生じる可能性がある。そこで、第１の移動処理部３１は、法線ベクトルＴ0~の方向に断面曲線Ｃ0(s)を平行移動させた後、さらに、移動後の断面曲線Ｃ1(s)~を含む平面Ｐ1内で断面曲線Ｃ1(s)~を平行移動させてＣ1(s)とすることにより、断面曲線Ｃ1(s)と形状データとの誤差が最小となるようにする。
　次に、移動後の断面曲線Ｃ1(s)を含む平面Ｐ1の法線ベクトルＴ1~（図３において図示せず）の方向に対し、所定の間隔で断面曲線Ｃ1(s)を平行移動し、得られた移動後の断面曲線Ｃ(s)をＣ2(s)~（図３において図示せず）とする。さらに、第１の移動処理部３１は、移動後の断面曲線Ｃ2(s)~を含む平面Ｐ2内で、断面曲線Ｃ2(s)~を平行移動させてＣ2(s)とすることにより、断面曲線Ｃ2(s)と形状データとの誤差が最小となるようにする。
　以下、この処理を繰り返し実行する。図３の例では、この処理をｄ回行った例を示している。これにより、第１の移動処理部３１は、所定の間隔毎に存在する複数の断面曲線Ｃi(s)（ｉ＝０～ｄ）と、各断面曲線Ｃi(s)を含む平面Ｐi（ｉ＝０～ｄ）と、断面曲線Ｃi(s)の移動方向に対する複数の接線ベクトルＴi（ｉ＝０～ｄ）とが得られる。
　上述のように、ｉ＝０～ｄの全てについて、同じ法線ベクトルＴ0~の方向に断面曲線Ｃi(s)~を単純に平行移動させるだけだと（Ｔi（ｉ＝０～ｄ）＝Ｔ0~の場合）、平行移動を重ねるにつれて（ｉが大きくなるにつれて）、断面曲線Ｃi(s)~は形状データから離れて誤差が大きくなってしまう可能性がある。
　これに対し、第１の移動処理部３１は、平面Ｐi（ｉ＝０～ｄ）の法線ベクトル方向に断面曲線Ｃi(s)を平行移動させた後、その平行移動後の断面曲線Ｃi+1(s)~を平面Ｐi+1内でさらに平行移動させることにより、断面曲線Ｃi+1(s)と形状データとの誤差が最小となるようにする。
　これにより、図３（ｂ）に示すように、必ずしもＴi（ｉ＝１～ｄ）＝Ｔ0~になるとは限らなくなり、各断面曲線Ｃi(s)が形状データに極力フィッティングするような条件を満たす接線ベクトルＴi（ｉ＝０～ｄ）が求められる。こうして所定の間隔毎に求められる接線ベクトルＴiは、断面曲線Ｃi(s)上の一点が平行移動した軌跡の方向を示すベクトルに相当する。
　第２の移動処理部３２は、第１の移動処理部３１により算出された接線ベクトルＴi（ｉ＝０～ｄ）を初期ベクトルとして、所定の条件を満たすように、所定の間隔毎の断面曲線Ｃi(s)をそれぞれ平行移動および回転移動によってさらに移動させる。ここで、第２の移動処理部３２は、断面曲線Ｃi(s)を平行移動させたときに接線ベクトルＴiを更新し、更新後の接線ベクトルＴiに基づいて回転移動を行う。
　第１の実施形態では、所定の条件として次の２つの条件を設定する。１つ目の条件は、特許請求の範囲における第１の条件に相当するものであり、第２の移動処理部３２によって移動された所定の間隔毎の断面曲線Ｃi(s)’上の一点により形成される軌跡（断面曲線Ｃ(s)を仮にスイープさせるとした場合のガイド線に相当するもの。以下、疑似ガイド線Ｓc（図２参照）という）を、捩れの小さな曲率線に近づけるための条件である。２つ目の条件は、特許請求の範囲における第３の条件に相当するものであり、曲面生成部４により生成される目的曲面Ｓ(u, v)と、形状データ入力部１により入力された形状データで表される曲面との誤差を最小化するための条件である。
　第２の移動処理部３２が行う処理の詳細に関して、まず、断面曲線Ｃi(s)の回転移動について説明する。第２の移動処理部３２は、断面曲線Ｃi(s)上の点Ｃi(s*)を回転中心として、接線ベクトルＴi-1が次の接線ベクトルＴiに一致するように断面曲線Ｃi(s)を回転する。なお、パラメータs*は、断面曲線Ｃi(s)上の一点（例えば、中点）を与える値である。
　ここで、断面曲線Ｃi(s)の回転軸Ａ(i)および回転角θはそれぞれ、Ａ(i)＝Ｔi-1×Ｔiおよびθ＝ｃｏｓ-1（Ｔi-1・Ｔi）とする。Ａ(i)＝Ｔi-1×Ｔiは、接線ベクトルＴi-1と次の接線ベクトルＴiとの外積を回転軸Ａ(i)とするという意味である。したがって、回転軸Ａ(i)は、接線ベクトルＴi-1，Ｔiの両方に直交する方向を示すことになるから、接線ベクトルＴiの方向には回転しないということになる。
　このとき、回転行列Ｒは、次の（式１）で与えられる。なお、ｌ＝Ａx (i)、ｍ＝Ａy (i)、ｎ＝Ａz (i)とする。θ＝０の場合、Ａ(i)はゼロベクトルとなるが、Ｒは恒等変換（単位行列）となるため、当該変換を適用すれば問題はない。
Figure JPOXMLDOC01-appb-M000001
　上記（式１）に示す回転行列Ｒは、接線ベクトルＴiの方向の回転がない回転最小枠（ＲＭＦ: rotation minimizing frame）を近似していると考えることができる。なぜなら、ある曲線の接線ベクトルＴに垂直な方向を回転軸Ａとする任意回転は、Ａ・Ｔ＝０であり、これはＲＭＦの定義だからである。本実施形態の回転Ａ(i)に関して、接線ベクトルＴi-1，Ｔiの間隔を狭くして、Ｔiが想定している疑似ガイド線Ｓcの接線ベクトルを十分近似できる場合は、行列ＲはＲＭＦによる回転を意味する。
　図４は、第２の移動処理部３２によって実行される断面曲線Ｃi(s)の移動（平行移動および回転移動）の一例を示す模式図である。図４（ａ）は、第１の移動処理部３１により断面曲線Ｃi(s)が平行移動された結果（図３（ｂ）と同じ状態）を示している。これに対し、図４（ｂ）は、図４（ａ）に示す断面曲線Ｃi(s)をさらに平行移動および回転移動させた結果を示している。図４（ｂ）に示す断面曲線Ｃi(s)’は、上述した２つの条件（特許請求の範囲における第１の条件および第３の条件）を満たすように移動された結果を示すものである。
　次に、１つ目の条件（第１の条件）について説明する。ここでは、クレイモデラがカーブ定規２００を滑らかに動かすという行為を目的曲面Ｓ(u, v)がもつ特徴に関連付けるために、次の２つ条件を想定する。
　Ｐ１）疑似ガイド線Ｓcは、生成しようとする目的曲面Ｓ(u, v)上の意図的な流れを意味し、その流れは曲率線として特定される。曲率線は、目的曲面Ｓ(u, v)における最小の主曲率に対する主方向による曲面の流れを示すものである。
　Ｐ２）疑似ガイド線Ｓcは捩れが小さい。
　条件Ｐ１は、平行移動および回転移動に基づき移動された所定の間隔毎の断面曲線Ｃi(s)’上の一点Ｃi(s*)’をそれぞれ結んでできる軌跡である疑似ガイド線Ｓcが、曲面生成部４によって生成される目的曲面Ｓ(u, v)の曲率線になるという条件である。
　条件Ｐ２は、曲面生成部４によって生成される目的曲面Ｓ(u, v)の法線ベクトルから疑似ガイド線Ｓcの主法線ベクトルへの角度が、疑似ガイド線Ｓcの軌跡上において一定になるという条件である。
　なお、曲面生成部４は、断面曲線移動部３により移動された断面曲線Ｃ(s)’の軌跡によって定義される境界線からの内挿により曲面を生成する。すなわち、曲面生成部４は、図４（ｂ）において、移動の始点および終点におけるｕ方向に対する２つの断面曲線Ｃ0(s)’，Ｃd(s)’と、ｖ方向に対する２つの境界線Ｃs’，Ｃe’とを用いて、これらの４辺からの内挿により目的曲面Ｓ(u, v)を、例えば公知のＣｏｏｎｓパッチにより生成する。
　第２の移動処理部３２は、条件Ｐ１を達成するために、疑似ガイド線Ｓcを曲率線に近づける処理を行う。ここで、曲率線は、次の（式２）で与えられる。
　（ＥＭ－ＦＬ）du2 ＋（ＥＮ－ＧＬ）dudv＋（ＦＮ－ＧＭ）dv2＝０　・・・（式２）
　ここで、Ｅ，Ｆ，Ｇはそれぞれ曲面の第１基本形式の係数（第１基本量）であり、Ｌ，Ｍ，Ｎはそれぞれ曲面の第２基本形式の係数（第２基本量）である。疑似ガイド線Ｓcは、断面曲線の一点Ｃ(s*)の軌跡に対応する目的曲面Ｓ(u, v) ｖ方向の等パラメータ線Ｓ(u*, v)（Ｃ(s*)＝Ｓ(u*, 0)）であることから、条件Ｐ１はｄｕ＝０により次の（式３）で表すことができる。この（式３）において、ＶNは目的曲面Ｓ(u, v)の法線ベクトルを示す。
Figure JPOXMLDOC01-appb-M000002
　次に、条件Ｐ２について説明する。曲面上のある曲線Ｃの測地捩率τgは、次の（式４）で与えられる。
　τg＝（ｄα／ｄｓ）－τ　・・・（式４）
ここで、τは曲線Ｃの通常の捩率であり、αは曲面の法線ベクトルＶN から曲線Ｃの主法線ベクトルへの角度である。さらに、曲率線はτg＝０である。したがって、条件Ｐ２を実現するためにｄα／ｄｓをできるだけ小さくすることにより、疑似ガイド線Ｓcは捩率τの小さい曲線となる。
　以上のことから、第１の条件に関しては、目的曲面Ｓ(u, v)に対して、疑似ガイド線Ｓcが曲率線になっており、疑似ガイド線Ｓcの捩率τが小さいことを、次の（式５）に示す関数Ψにより評価する。ここで、ｗ1，ｗ2は重み係数である。また、Ｎ2は疑似ガイド線Ｓc上から抽出するサンプリング点の数である。このサンプリング点は、疑似ガイド線Ｓc上から一様に抽出するのが好ましい。
Figure JPOXMLDOC01-appb-M000003
　次に、２つ目の条件（第３の条件）について説明する。第３の条件は、曲面生成部４によって生成される目的曲面Ｓ(u, v)と、形状データ入力部１により入力された形状データで表される曲面との誤差を最小化するという条件である。
　この第３の条件に関しては、目的曲面Ｓ(u, v)と形状データ上のサンプリング点Ｑiとの誤差が小さいことを、次の（式６）により評価する。ここで、Ｎは形状データ上から抽出するサンプリング点Ｑiの数である。なお、サンプリング点Ｑiは、形状データの全体から偏りなく抽出するのが好ましく、曲面の全体にいきわたるように適度に抽出する。
Figure JPOXMLDOC01-appb-M000004
　第２の移動処理部３２は、上述した第１の条件および第３の条件を考慮して、次の（式７）で定義される目的関数Ｊを最小化するように断面曲線Ｃ(s)を移動させる。
Figure JPOXMLDOC01-appb-M000005
　以上詳しく説明したように、クレイモデルなどの形状データと、当該形状データにフィッティングした断面曲線Ｃ(s)とを入力して第１の実施形態の処理を実行すれば、目的曲面Ｓ(u, v)を生成するために断面曲線Ｃ(s)を移動させたときに当該断面曲線Ｃ(s)上の一点の軌跡により定義される軌跡（疑似ガイド線Ｓc）が、モデル形状に基づいて生成しようとする目的曲面Ｓ(u, v)における最小の主曲率に対する主方向による曲面の流れを表す曲率線として表現できる。
　すなわち、第１の実施形態によれば、ガイド線を入力データとして事前に与えることが要求される従来のスイープ法とは異なり、形状データと断面曲線Ｃ(s)とを与えれば、当該曲線Ｃ(s)の移動による点Ｃ(s*)の軌跡が、生成しようとする目的曲面Ｓ(u, v)の曲率線となり、かつ捩れも小さいという性質をもつような目的曲面Ｓ(u, v)が生成される。つまり、疑似ガイド線Ｓcが捩率の小さい曲率線となるような良好な目的曲面Ｓ(u, v)を、多くの試行錯誤を繰り返すことなく効率的に生成することができる。
　なお、上記第１の実施形態において、（式５）の重み係数ｗ2をゼロとすることにより、疑似ガイド線Ｓcが目的曲面Ｓ(u, v)の曲率線になるという条件Ｐ１のみを第１の条件として用いるようにしてもよい。条件Ｐ１のみとしても、ＲＭＦを近似する回転によって断面曲線Ｃ(s)の移動が行われているので、疑似ガイド線Ｓcの捩率が陽に大きくなることを避けることができる。
　また、上記第１の実施形態では、自動車のクレイモデルなどの測定データから目的曲面Ｓ(u, v)を生成する例について説明したが、入力データとして使用する形状データは、クレイモデルの測定データに限定されるものではない。例えば、市販車の測定データから目的曲面Ｓ(u, v)を生成するようにしてもよい。また、形状データは測定データに限定されるものではなく、ＣＡＤデータやその離散近似したものでも構わない。さらに、対象とするモデルは自動車に限らず、あらゆる物を対象とすることが可能である。
　また、上記第１の実施形態では、疑似ガイド線Ｓcが捩率の小さな曲率線となるための第１の条件に加え、生成される目的曲面Ｓ(u, v)とクレイモデルなどの形状データとの誤差を最小化するための第３の条件を満たすように各断面曲線Ｃi(s)の位置を調整する例について説明したが、式（７）において第３の条件を表す第１項と第１の条件を表す第２項とに重み係数を設定し、両条件の重みを調整できるようにしてもよい。例えば、入力データとして、現物のモデルに対する測定データではなく、新たに作成しようとするデザインのコンセプトを表したＣＡＤデータを形状データ入力部１より入力して目的曲面Ｓ(u, v)を生成する場合などは、第３の条件の重みを小さくして処理を行うことも可能である。
（第２の実施形態）
　次に、本発明による第２の実施形態を図面に基づいて説明する。図５は、第２の実施形態による曲面生成装置の機能構成例を示すブロック図である。なお、この図５において、図１に示した符号と同一の符号を付したものは同一の機能を有するものであるので、ここでは重複する説明を省略する。第２の実施形態による曲面生成装置は、第２の移動処理部３２に代えて第２の移動処理部３２’を備えている。
　上記第１の実施形態では、形状データ（クレイモデルなどの測定データ）から再構築する目的曲面Ｓ(u, v)の品質を向上させるために、断面曲線Ｃ(s)を移動させたときにできる各断面曲線Ｃi(s)’上の点Ｃi(s*)’の軌跡が、目的曲面Ｓ(u, v)における捩れの小さい曲率線となるように、第１の条件を満たすように各断面曲線Ｃi(s)’の位置を調整する例について説明した。
　第２の実施形態では、これに加え、目的曲面Ｓ(u, v)のｖ方向の境界線Ｃs，Ｃeに、第２の条件を課すことによって各断面曲線Ｃi(s)の位置を調整するものである。第２の条件は、境界線Ｃs’，Ｃe’の曲線が曲率単調性を持つという条件である。
　すなわち、第１の移動処理部３１によって行われる断面曲線Ｃ(s)の平行移動、および、第２の移動処理部３２によって行われる断面曲線Ｃ(s)の平行移動および回転移動では、それによって生成される曲面Ｓ(u, v)のｖ方向を定義する２つの境界線Ｃs，Ｃeは、複数の断面曲線Ｃi(s)（ｉ＝０～ｄ）の両端点を補間あるいは近似したものに過ぎず、意匠デザインで求められる曲率単調性が満たされているとは限らない。そこで、両境界線Ｃs，Ｃeが曲率単調性を持つように平面Ｐi内で各断面曲線Ｃi(s)を微調整することによって、最適な接線ベクトルＴiを計算する。
　曲面Ｓ(u, v)のｖ方向の２つの境界線Ｃs（＝Ｓ(０, v)），Ｃe（＝Ｓ(１, v)）の曲率が単調増加するとした場合、その評価を行うために、次の（式８）に示す関数Φを導入する。ここで、κiは各境界線Ｃs，Ｃe上で等間隔となるように設定された各サンプリング点における曲率を示し、δ(κi) は（式９）で定義されたペナルティ関数である。（式８）の第１項が一方の境界線Ｃsに関する曲率単調性を定義し、第２項が他方の境界線Ｃeに関する曲率単調性を定義している。Ｎ1はサンプリング点の数を示している。
Figure JPOXMLDOC01-appb-M000006
　第２の実施形態では、第２の移動処理部３２’は、（式７）に代えて（式１０）のように定義される目的関数Ｊを最小化するように断面曲線Ｃ(s)を移動させる。ここで、ε1，ε2は重み係数である。
Figure JPOXMLDOC01-appb-M000007
　以上のように構成した第２の実施形態によれば、断面曲線Ｃi(s)’上の一点Ｃi(s*)’（例えば、中点）の軌跡が、疑似ガイド線Ｓcとして望ましい曲線（捩れの小さい曲率線）となるだけでなく、断面曲線Ｃi(s)’の両端点の軌跡によってできる各境界線Ｃs’，Ｃe’が曲率単調性を有する曲線となる。各境界線Ｃs’，Ｃe’が曲率単調性を有する結果、疑似ガイド線Ｓcも曲率単調性を有するものとなる。このため、形状データから再構築する目的曲面Ｓ(u, v)の品質を更に向上させることができる。
　なお、第２の実施形態においても第１の実施形態と同様、（式５）の重み係数ｗ2をゼロとすることにより、条件Ｐ１のみを第１の条件として用いるようにしてもよい。また、入力データとなる形状データは、クレイモデルなどの測定データに限定されるものではなく、ＣＡＤデータやその離散近似したものでも構わない。また、（式１０）において、第３の条件を表す第１項にも重み係数を設定するようにしてもよい。
　また、上記第１および第２の実施形態では、曲面生成部４により生成する目的曲面Ｓ(u, v)をＢ－ｓｐｌｉｎｅ曲面で表現する例について説明したが、本発明はこれに限定されない。例えば、目的曲面Ｓ(u, v)をベジエ曲面で表現するようにしてもよい。
　また、上記第１および第２の実施形態では、クレイモデルを作成する際に使用するカーブ定規２００を模した断面曲線Ｃ(s)を設定する例について説明したが、本発明はこれに限定されない。すなわち、断面曲線Ｃ(s)は、生成したい目的曲面Ｓ(u, v)の形状に合わせて設定すればよい。
　また、上記第１および第２の実施形態では、断面曲線Ｃ(s)の移動を、第１の移動処理部３１と第２の移動処理部３２とに分けて２段階で行う例について説明したが、本発明はこれに限定されない。すなわち、断面曲線Ｃ(s)を所定の間隔で平行移動させながら回転移動させることによって、１段階の処理として断面曲線Ｃ(s)を移動させるようにしてもよい。
　また、上記第１および第２の実施形態では、断面曲線Ｃ(s)の回転移動を、（式１）で定義したが、この（式１）に限定されない。すなわち、ＲＭＦを近似する代替式は他にも幾つかある。また、（式１）は、回転軸Ａ(i)の設定方法により、任意の軸まわりの回転にも対応可能である。
　また、上記第１および第２の実施形態では、目的曲面Ｓ(u, v)の生成にはＣｏｏｎｓパッチによる境界線からの内挿を利用したが、本発明はこれに限定されない。例えば、断面曲線Ｃ(s)上の一点Ｃ(s*)（例えば、中点）の軌跡をB-spline曲線などで近似し、当該近似曲線に沿って断面曲線Ｃ(s)をスイープしても、上記の議論はすべて適用できる。
　また、上記第１および第２の実施形態では、第１の移動処理部３１の処理に関して、平面Ｐiの法線ベクトル方向に断面曲線Ｃi(s)を平行移動させた後、その平行移動後の断面曲線Ｃi+1(s)~を平面Ｐi+1内でさらに平行移動させる例について説明したが、法線ベクトル方向への平行移動のみとしてもよい。この場合、断面曲線Ｃi+1(s)と形状データとの誤差が最小化された接線ベクトルＴiが初期ベクトルとして生成されないものの、第２の移動処理部３２による最適化計算によって、第１の条件～第３の条件を満たす目的曲面Ｓ(u, v)を生成することができる。なお、第１の移動処理部３１において、断面曲線Ｃi+1(s)と形状データとの誤差が最小となるような接線ベクトルＴiを初期ベクトルとして生成しておくことにより、第２の移動処理部３２による最適化計算の処理時間を短くできる等のメリットを有する。
　その他、上記第１および第２の実施形態は、何れも本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその要旨、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
　１　形状データ入力部
　２　断面曲線設定部
　３　断面曲線移動部
　４　曲面生成部
　３１　第１の移動処理部
　３２，３２’　第２の移動処理部

Claims (12)

  1. 　モデルの形状データを入力する形状データ入力部と、
    　上記形状データの一部にフィッティングした断面曲線を設定する断面曲線設定部と、
    　上記断面曲線設定部により設定された上記断面曲線を、所定の条件を満たすように平行移動および回転移動させる断面曲線移動部と、
    　上記断面曲線移動部により上記断面曲線が動く軌跡によって定義された曲面を生成する曲面生成部とを備え、
    　上記所定の条件は、上記断面曲線移動部により移動された上記断面曲線上の一点により形成される軌跡が、生成しようとする曲面の曲率線になるという第１の条件であることを特徴とする曲面生成装置。
  2. 　上記断面曲線移動部は、
    　上記断面曲線を平行移動させ、上記断面曲線上の一点が平行移動した軌跡の方向を示す接線ベクトルを所定の間隔毎に求める第１の移動処理部と、
    　上記第１の移動処理部により算出された上記接線ベクトルを初期ベクトルとして、上記所定の条件を満たすように、上記所定の間隔毎の上記断面曲線をそれぞれ平行移動および回転移動によってさらに移動させる第２の移動処理部とを備えたことを特徴とする請求項１に記載の曲面生成装置。
  3. 　上記第１の移動処理部は、上記形状データとの誤差が最小となるように上記断面曲線を平行移動させることを特徴とする請求項２に記載の曲面生成装置。
  4. 　上記第１の条件は、上記断面曲線移動部により移動された上記断面曲線上の一点により形成される軌跡が、上記生成しようとする曲面の曲率線になるとともに、上記曲面生成部によって生成される曲面の法線ベクトルから上記軌跡の主法線ベクトルへの角度が、上記軌跡上において一定になるという条件であることを特徴とする請求項１～３の何れか１項に記載の曲面生成装置。
  5. 　上記断面曲線移動部は、上記断面曲線を回転最小枠による回転移動によって移動させることを特徴とする請求項１～４の何れか１項に記載の曲面生成装置。
  6. 　上記所定の条件は、上記断面曲線移動部によって上記断面曲線を移動させたときに上記断面曲線の両端が移動してできる曲線が曲率単調性を持つという第２の条件を更に含むことを特徴とする請求項１～５の何れか１項に記載の曲面生成装置。
  7. 　上記所定の条件は、上記曲面生成部によって生成される曲面と、上記形状データ入力部により入力された上記形状データで表される曲面との誤差を最小化するという第３の条件を更に含むことを特徴とする請求項１～６の何れか１項に記載の曲面生成装置。
  8. 　モデルの形状データを入力する形状データ入力手段、
    　上記形状データの一部にフィッティングした断面曲線を設定する断面曲線設定手段、
    　上記断面曲線設定手段により設定された上記断面曲線を、所定の条件を満たすように平行移動および回転移動させる断面曲線移動手段、および
    　上記断面曲線移動手段により上記断面曲線が動く軌跡によって定義された曲面を生成する曲面生成手段、
    としてコンピュータを機能させ、
    　上記所定の条件が、上記断面曲線移動手段により移動された上記断面曲線上の一点により形成される軌跡が、生成しようとする曲面の曲率線になるという第１の条件であることを特徴とする曲面生成用プログラム。
  9. 　上記第１の条件は、上記断面曲線移動手段により移動された上記断面曲線上の一点により形成される軌跡が、上記生成しようとする曲面の曲率線になるとともに、上記曲面生成手段によって生成される曲面の法線ベクトルから上記軌跡の主法線ベクトルへの角度が、上記軌跡上において一定になるという条件であることを特徴とする請求項８に記載の曲面生成用プログラム。
  10. 　上記断面曲線移動手段は、上記断面曲線を回転最小枠による回転移動によって移動させることを特徴とする請求項８または９に記載の曲面生成用プログラム。
  11. 　上記所定の条件は、上記断面曲線移動手段によって上記断面曲線を移動させたときに上記断面曲線の両端が移動してできる曲線が曲率単調性を持つという第２の条件を更に含むことを特徴とする請求項８～１０の何れか１項に記載の曲面生成用プログラム。
  12. 　上記所定の条件は、上記曲面生成手段によって生成される曲面と、上記形状データ入力手段により入力された上記形状データで表される曲面との誤差を最小化するという第３の条件を更に含むことを特徴とする請求項８～１１の何れか１項に記載の曲面生成用プログラム。
PCT/JP2017/042636 2017-01-11 2017-11-28 曲面生成装置および曲面生成用プログラム WO2018131304A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP17891292.9A EP3550456A4 (en) 2017-01-11 2017-11-28 DEVICE FOR GENERATING A CURVED SURFACE AND PROGRAM FOR GENERATING A CURVED SURFACE
US16/477,061 US11145120B2 (en) 2017-01-11 2017-11-28 Curved surface generation device and curved surface generation program
US17/476,163 US20220005274A1 (en) 2017-01-11 2021-09-15 Curved surface generation device and curved surface generation program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017-002602 2017-01-11
JP2017002602A JP6751029B2 (ja) 2017-01-11 2017-01-11 曲面生成装置および曲面生成用プログラム

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US16/477,061 A-371-Of-International US11145120B2 (en) 2017-01-11 2017-11-28 Curved surface generation device and curved surface generation program
US17/476,163 Continuation US20220005274A1 (en) 2017-01-11 2021-09-15 Curved surface generation device and curved surface generation program

Publications (1)

Publication Number Publication Date
WO2018131304A1 true WO2018131304A1 (ja) 2018-07-19

Family

ID=62840014

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/042636 WO2018131304A1 (ja) 2017-01-11 2017-11-28 曲面生成装置および曲面生成用プログラム

Country Status (4)

Country Link
US (2) US11145120B2 (ja)
EP (1) EP3550456A4 (ja)
JP (1) JP6751029B2 (ja)
WO (1) WO2018131304A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022030104A1 (ja) * 2020-08-07 2022-02-10 日本ユニシス株式会社 曲面生成装置および曲面生成用プログラム

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6998278B2 (ja) 2018-06-13 2022-02-10 三洋化成工業株式会社 樹脂集電体、積層型樹脂集電体、及び、リチウムイオン電池
EP3675063A1 (en) * 2018-12-29 2020-07-01 Dassault Systèmes Forming a dataset for inference of solid cad features
EP3675062A1 (en) 2018-12-29 2020-07-01 Dassault Systèmes Learning a neural network for inference of solid cad features
TWI760675B (zh) * 2020-01-06 2022-04-11 財團法人工業技術研究院 加工路徑缺陷檢測方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS575109A (en) * 1980-06-10 1982-01-11 Fanuc Ltd Curved surface forming method
JPH05250444A (ja) 1992-03-06 1993-09-28 Ricoh Co Ltd ３次元曲面形状生成装置
JP2938909B2 (ja) 1989-12-28 1999-08-25 株式会社日立製作所 曲面生成方法及びその装置
JP2000011210A (ja) 1998-06-19 2000-01-14 Sharp Corp 曲面作成装置および方法ならびに曲面作成プログラムを記録した記録媒体
JP2005149245A (ja) 2003-11-17 2005-06-09 Mitsubishi Heavy Ind Ltd Ｃａｄシステム、曲面解析装置、曲面再生装置、その方法及びそのプログラム
JP2007156604A (ja) * 2005-12-01 2007-06-21 Mitsubishi Heavy Ind Ltd 曲面作成方法及び曲面作成プログラム並びに３次元形状処理装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5982490A (en) * 1997-02-04 1999-11-09 Nikon Corporation Apparatus and method for wavefront absolute calibration and method of synthesizing wavefronts
US7196702B1 (en) * 1998-07-23 2007-03-27 Freedesign, Inc. Geometric design and modeling system using control geometry

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS575109A (en) * 1980-06-10 1982-01-11 Fanuc Ltd Curved surface forming method
JP2938909B2 (ja) 1989-12-28 1999-08-25 株式会社日立製作所 曲面生成方法及びその装置
JPH05250444A (ja) 1992-03-06 1993-09-28 Ricoh Co Ltd ３次元曲面形状生成装置
JP2000011210A (ja) 1998-06-19 2000-01-14 Sharp Corp 曲面作成装置および方法ならびに曲面作成プログラムを記録した記録媒体
JP2005149245A (ja) 2003-11-17 2005-06-09 Mitsubishi Heavy Ind Ltd Ｃａｄシステム、曲面解析装置、曲面再生装置、その方法及びそのプログラム
JP2007156604A (ja) * 2005-12-01 2007-06-21 Mitsubishi Heavy Ind Ltd 曲面作成方法及び曲面作成プログラム並びに３次元形状処理装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3550456A4

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022030104A1 (ja) * 2020-08-07 2022-02-10 日本ユニシス株式会社 曲面生成装置および曲面生成用プログラム
JP7431123B2 (ja) 2020-08-07 2024-02-14 Ｂｉｐｒｏｇｙ株式会社 曲面生成装置および曲面生成用プログラム

Also Published As

Publication number Publication date
US20190333272A1 (en) 2019-10-31
US20220005274A1 (en) 2022-01-06
US11145120B2 (en) 2021-10-12
JP6751029B2 (ja) 2020-09-02
JP2018112869A (ja) 2018-07-19
EP3550456A1 (en) 2019-10-09
EP3550456A4 (en) 2020-08-05

Similar Documents

Publication Publication Date Title
WO2018131304A1 (ja) 曲面生成装置および曲面生成用プログラム
CN113228115A (zh) 网格几何体到水密边界表示的转换
KR100707841B1 (ko) 3차원 목표곡선을 이용한 곡면 변형 장치 및 그 방법
WO2019202813A1 (ja) 曲面生成装置および曲面生成用プログラム
JP5271825B2 (ja) 案内曲線で定義される面の設計をなすコンピュータ実行方法
US20120221297A1 (en) Global Deformation for a Modeled Object
JP4418390B2 (ja) ３次元形状処理装置及び曲面生成プログラム並びに方法
IL167509A (en) Cad system and cad program
CN112396690B (zh) 基于改进型向心参数化法的曲面高精重构方法
JP2020102215A (ja) ユーザインタラクションを介した３ｄモデルオブジェクトのデサイニング
JP7446579B2 (ja) 曲面のフィッティング処理方法、フィッティング処理装置およびフィッティング処理プログラム、並びに、該フィッティング処理プログラムを記憶したコンピュータ読取可能な記憶媒体
Biermann et al. Direct free-form deformation of NC programs for surface reconstruction and form-error compensation
WO2022030104A1 (ja) 曲面生成装置および曲面生成用プログラム
JP4981313B2 (ja) ３次元形状処理装置及び曲面作成プログラム並びに曲面作成方法
JP4175470B2 (ja) Ｃａｄシステム、曲面解析装置、曲面再生装置、その方法及びそのプログラム
JP2016207055A (ja) 形状変形装置および形状変形用プログラム
JP6449703B2 (ja) 形状変形装置および形状変形用プログラム
CN108628255A (zh) 一种指令点纠正平滑处理方法
Cheng et al. Complex CAD surface shape design using semantic features
KR101452374B1 (ko) 변형 축을 이용한 이미지 변형 방법 및 장치
JP3408866B2 (ja) 曲線と曲面間の抜き勾配面生成方法及びその装置
JP7036207B2 (ja) メッシュ生成システム、メッシュ生成プログラムおよびメッシュ生成方法
Hahmann et al. Hierarchical modeling with parametric surfaces
US20110148870A1 (en) Method and apparatus for generating graphic tentacle motions
Rath et al. DESIGN OF FREEFORM CURVE MODELING BY IMPORTING NEUTRAL FILE GEOMETRIC DATA THROUGH THE DEVELOPED IGES READER

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17891292

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2017891292

Country of ref document: EP

Effective date: 20190626