JPH0638269B2 - Method and device for supporting coordinate grid creation - Google Patents

Description

Description: TECHNICAL FIELD The present invention relates to a coordinate grid creation support method and apparatus for an analysis target in a numerical analysis using a computer, and particularly, to a working time required to create a coordinate grid. The present invention relates to a coordinate grid creation support method and device suitable for creating a coordinate grid optimum for analysis in accordance with the purpose of analysis.

[Conventional technology]

In the fields of structural analysis, fluid analysis, electromagnetic field analysis, etc., general-purpose numerical analysis programs using the finite element method, the difference method, etc. have been created to enable highly accurate analysis. However, when a designer wants to analyze a physical phenomenon using these analysis programs, most of the analysis work time is spent on creating coordinate grids, numbering, calculation of nodal coordinates, and data such as characteristic data. There is a problem that. Especially in three-dimensional analysis, data creation requires 80
~ 90% is spent. Therefore, various preprocessors have been developed in order to save the data creation. Preprocessor for mechanical CAE of the finite element method (for example, Journal of Japan Society of Mechanical Engineers, Vol. 88, No. 794, pp.
For example, the shape modeling system section of the CAE system introduced in 29-35) is used to construct the geometric shape of the analysis target in the computer using the graphical display in an online interactive form, and based on this. , The surface of the shape and the internal coordinate grid points are determined by interpolation, and data such as coordinate grid creation, numbering, and physical property data are automatically created. Thus, in the conventional CAE system method,
Generally, the analysis area is displayed on the graph display, and the displayed figure is picked with a light pen to input the data such as the number of coordinate grid points and physical property data. Therefore, the coordinate grid points on the surface and inside of the shape were determined.

[Problems to be solved by the invention]

In numerical analysis, as the number of coordinate grid points increases, the calculation accuracy improves, but at the same time, the required calculation time and memory capacity also increase, so the number of coordinate grid points that can be analyzed is limited by the limited computer resources. There is. Therefore, in order to improve the accuracy and convergence of numerical analysis under a finite number of grid points, a region important for analysis according to the purpose of analysis (for example, near the boundary in fluid analysis, stress concentration part in structural analysis). ,
In electromagnetic field analysis, it is necessary to concentrate the coordinate grid in a place where the potential gradient is large, etc., and especially when using the difference method, it is necessary to orthogonalize the coordinate grid in the vicinity. However, in the above-mentioned prior art, after constructing the geometric shape of the analysis object in the computer, the coordinate grid is created by interpolation based on the input data of the number of divisions, so that the degree of concentration, orthogonality, etc. are controlled. However, there is a problem that it is not always easy to create an optimal coordinate grid according to the purpose of analysis.

An object of the present invention is to provide a coordinate grid creation support method and apparatus that can reduce the working time of creating a coordinate grid and easily create the optimum coordinate grid for analysis according to the analysis purpose.

[Means for solving problems]

In order to achieve the above object, the present invention interactively inputs basic data such as the number of grid points of a coordinate grid and a fixed position (boundary position), and at the same time, the distance between adjacent coordinate grid points, 1
Optimization that controls the layout characteristics such as the smoothness, concentration, and orthogonality of the coordinate grid expressed by the bending angle at the grid points of the two coordinate grid lines and the angle at which the two coordinate grid lines intersect at the coordinate grid points. Input parameters, perform coordinate grid creation calculation, and evaluate the values of various evaluation functions (smoothness, concentration, orthogonality) of the coordinate grid obtained as a result, together with the coordinate grid diagram. The present invention proposes a coordinate grid creation support method and apparatus for displaying it as guidance, and for the analysis operator to view the guidance to further change the value of the optimization parameter and repeat the coordinate grid creation calculation.

The present invention particularly displays, in parallel, a coordinate grid diagram obtained by the creation calculation, a value of a calculated evaluation function, a menu for selecting the optimization parameter, and a control menu of the drawing to be displayed, and a screen is displayed. The present invention proposes a coordinate grid creation support device capable of quickly and easily creating a coordinate grid optimum for analysis while interacting.

[Action]

By repeating the coordinate grid creation calculation while looking at such evaluation guidance, an optimum coordinate grid can be created according to the purpose of the target numerical simulation.

For example, when solving the spatial distribution of a physical quantity, the coordinate grid is concentrated at a place where the gradient of the physical quantity is considered to be large, orthogonality is maintained to improve convergence, and problems with internal boundaries are solved. It is possible to smoothly connect the coordinate grids near the boundary line by using the above-mentioned optimization parameters.

In addition, interactively changing the parameter values and creating a coordinate grid while viewing the results and the evaluation guidance, the know-how of the analysis operator can be incorporated with a small amount of work, and the numerical simulation can be performed in a short time. It is possible to create an optimum coordinate grid.

Then, by using this optimum coordinate grid, the accuracy of numerical simulation can be improved.

〔Example〕

Hereinafter, the present invention will be described in more detail with reference to preferred embodiments. FIG. 1 shows the configuration of a coordinate grid creation support device 1 of the present invention. The support device 1 of the present invention includes a display device (for example, a CRT) 2, an image processing device 3, an arithmetic processing device (for example, an electronic computer) 4, a storage device 5, and an operation panel (for example, a keyboard) 6. The display device 2 is connected to the arithmetic processing device 4 via the image processing device 3, and the storage device 5 and the operation panel 6 are also connected to the arithmetic processing device 4.

FIG. 2 shows the processing procedure 4a stored in the arithmetic processing unit 4.
The contents of 4f are shown. These processing procedures correspond to the flow of processing procedures shown in the flow chart of FIG.
4a is a procedure for fetching parameters for optimizing the coordinate grid, 4b is a procedure for numerically calculating the coordinate grid using the parameters of 4a, and 4c is calculating an evaluation function for issuing evaluation guidance from the result of 4b. Procedure 4d is 4b
And 4c, a procedure for performing graphic processing of the result, 4e, a procedure for outputting the result of 4d to the display device 2 through the image processing device 3,
4f is a procedure for outputting the finally obtained coordinate grid data to the storage device 5.

FIG. 3 shows the contents of the storage device 5, where 5a is coordinate grid data and 5b is a portion for storing evaluation function data for 5a.

FIG. 4 shows a flow chart of a coordinate grid creation support method and a coordinate grid creation method using the apparatus of the present invention. In the figure, the initial coordinate grid creating portion 7 can be worked by using the supporting apparatus of the present invention, but is usually a general geometric modeling CA.
You can also use the E system. Now, assuming that the initial coordinate grid is given by such a system, the procedure thereafter will be described as a coordinate grid creating method 8 of the present invention. First,
At step 4a, optimization parameters are set in order to create an optimum coordinate grid for the target numerical simulation. This parameter is a parameter for controlling the smoothness, the degree of concentration, the orthogonality, and the layout characteristics of the coordinate grid, the details of which will be described later. In step 4b, the coordinate grid creation calculation is performed using this parameter. Next, in step 4c, an evaluation function representing the smoothness, the degree of concentration, the orthogonality, and the arrangement characteristic is calculated for the created coordinate grid. And
In step 4d, the value of the evaluation function is graphically processed (for example, contour line processing) and displayed together with the created coordinate grid. The analysis operator determines whether or not to repeat the optimization based on this guidance, and if it repeats, returns to the setting part of the optimization parameter, and if completed, at step 4f,
Output the final coordinate grid. By performing such an operation interactively, the analysis worker can efficiently create a coordinate grid according to the purpose of analysis by incorporating the know-how of the operator, reducing the work time and coordinating the coordinates. The accuracy of the numerical simulation can be improved by optimizing the lattice.

Here, an example of the coordinate grid creation calculation formula, its evaluation function, and optimization parameters adopted in the present invention will be shown. The coordinate grid creating method of the present invention uses the coordinates x,
Corresponding y, z and coordinates ξ, η, ζ in the mapping space in a one-to-one correspondence, solving a partial differential equation that associates the two, and curving coordinates in the real space corresponding to the Cartesian coordinate grid in the mapping space. This is a method of numerically obtaining the grid. According to this method, using the functions f, g, h, x = f (ξ, η, ζ) (1) y = g (ξ, η, ζ) (2) z = h ( ξ, η, ζ) ... One-to-one correspondence can be made as in (3).

Now, these functions are, for example, JUBrackbill (J.Com
According to p.Phys. 46 , 342-368 (1982), the two-dimensional case is expressed as follows.

As the evaluation function of optimization, the smoothness of the coordinate grid is I
If s, the concentration characteristic is Iv, and the orthogonality is Io, they can be expressed as follows.

_Is = ∫ Ω ^{[(▽ ξ) 2 + (▽} η) 2 ] dΩ ... (4) Iv = ∫ Ω WJdΩ ...... (5) Io = ∫ Ω ^{[(▽ ξ · ▽ η) 2} J 3 dΩ ... ( 6) Here, J is the Jacobian, and Ω is the target area. W is a weighting function, which is a function of J when only simple lattice concentration is performed. Qualitatively, these evaluation functions are: the distance between adjacent coordinate grid points in the curved coordinate system in the real space, the bending angle of one coordinate grid line at the coordinate grid point, and the two coordinate grid lines as coordinates. It can be considered to be represented by the angle of intersection in the grid. Here, if the optimization parameters are λs, λv, and λo, a desired coordinate grid can be created by creating a coordinate grid that optimizes I represented by I = λs Is + λv Iv + λo Io (7) The coordinate grid creation equation at this time is expressed as follows when λs = 1.

Here, the subscripts represent partial differentials, and a _{1 to} a ₃ and b ₁ to
b _{3, c} 1 ~c ₃ are constants. However, these constants are ak = a _sk + λv _avk + λo a _ok (10) bk = b _sk + λv b _vk + λo b _ok (11) (k = 1, 2) using the optimization parameters described above. , 3) ck = c _sk + λ v cv _k + λ oc _ok (12) Also, d ₁ and d ₂ are Here, P and Q are functions that control the coordinate interval, and JFThom
pson et al. (J. Comp. Phys. 15 , 299, (1974)).

Several other control methods for such a coordinate grid have been proposed (CDMobley, J. Comp. Phys. 34, 124).
(1980), DEPapantonis, Int.J.Num.Meth.Fluid
s, 5 , 245 (1985)), all of which relate to the smoothness, the degree of concentration, the orthogonality, etc. of the coordinate grid.

Now, in addition to the above-mentioned optimization parameters, when creating a coordinate grid, it is necessary to specify the number of coordinate grid points (number of mesh divisions), fixed grid points (boundary grid points) and their positions as basic operations. Is. These values are the boundary conditions when solving the equations (8) and (9) to obtain the coordinate grid points of the curved coordinate system. Therefore, the smoothness, concentration degree, and
The number of grid points and the positions of fixed grid points can be considered as optimization parameters of the coordinate grid in a broad sense, as well as the parameters for controlling the orthogonality and the arrangement characteristics.

The outline of the coordinate grid creating method according to the present invention has been described above. Next, the embodiment will be described in more detail with reference to an example of a screen that actually appears on the display device. FIG. 5 is a three-dimensional bird's-eye view showing the initial coordinate grid. Although this coordinate grid is shown by taking the piping of the heat exchanger as an example, it can be created by the supporting apparatus of the present invention or a general geometric modeling system as shown in 7 of FIG. Now, when the support device of the present invention is activated, a screen 9 as shown in FIG. 6 appears with FIG. 5 as the initial coordinate grid. Here, a coordinate grid diagram 9a similar to FIG. 5 is displayed on the left side of both sides, a menu 9b for selecting the optimization parameter and displaying the evaluation guidance is displayed on the right side, and a menu 9b is displayed on the lower right side. A control menu 9c for displaying a three-dimensional bird's-eye view and a two-dimensional cross-section is displayed. By using this drawing control menu, various drawings can be displayed according to the purpose. Two
Taking the three-dimensional sectional view as an example, a two-dimensional sectional view 9a at an arbitrary position can be displayed as shown in FIGS. 7 and 8. FIG. 7 shows a vertical sectional view of FIG. 6, and FIG. 8 shows a horizontal sectional view of a lower disk portion of the heat exchanger of FIG. In addition,
Here, the two-dimensional “cross-sectional view” is considered to include a plan view, an elevation view, a side view, and the like as a special example in which the virtual cutting plane does not intersect with the object.

FIG. 9 shows the flow of processing when the menu of optimization parameters is selected. There are four types of menus: division / fixation, smoothness, concentration, and orthogonality. When a menu is selected, each evaluation function and coordinate grid diagram are displayed as guidance. The analysis operator sequentially changes the optimization parameters while looking at the guidance, and executes the calculation for creating the coordinate grid. Then, the menu is selected again and the respective guidance is displayed. Repeat these operations,
Create an optimal coordinate grid according to the purpose of analysis.

Next, a method of inputting each optimization parameter and a method of displaying the evaluation guidance will be described in detail with reference to screens by taking an example of semiconductor device analysis. Figure 10 shows a two-dimensional structure MOSFET
3 is an overall view of the initial coordinate grid of FIG. FIG. 11 shows a flow chart of processing when the menu for displaying and changing the number of divisions and local grid fixation is selected. In this menu, when the display is made by selecting whether to display the number of divisions, when the display range is input, the number of coordinate grid divisions within the range is displayed on the screen. On the other hand, when dividing, input the division range and increase / decrease of the number of divisions, change the number of division data for coordinate grid creation calculation, and use linear interpolation to set the initial grid for coordinate grid creation calculation. create. When the initial grid is displayed by selecting whether to display it on the screen, it is displayed on the screen. Next, by inputting whether or not to fix the local grid, when fixed, the input of the range to be fixed is received, and the input data for the coordinate grid creation calculation is changed. FIG. 12 shows the details of the coordinate grid when the local portion of the calculated coordinate grid (near the upper right EFGH of the coordinate grid of FIG. 10) is further subdivided (particularly when the subdivision is centered around the curve EH). The display screen of an enlarged view is illustrated. FIG. 13 illustrates a display screen of the coordinate grid calculation result when the curve EH is a fixed coordinate grid line. By adjusting the number of divisions in this way, it is possible to easily increase the number of grid points at a position where the change of the variable is expected to be large. Further, when a locally fixed grid is set, a coordinate grid fitted to contour lines of environmental conditions and variables can be obtained. Therefore, the accuracy of numerical simulation is improved.

FIG. 14 shows the processing when the smoothness menu of the coordinate grid is selected. Depending on the selection input of whether or not to display the smoothness evaluation function, the type of figure,
In response to the input of the display range, the evaluation function of the smoothness of each grid position (for example, Is of the above-mentioned formula (4)) is displayed on the screen as contour lines. This contour line is very easy to understand, for example, if each area is colored. Next, by selecting whether or not to change the optimization parameter for controlling the smoothness, in the case of changing, the optimization parameter for controlling the smoothness is received and changed. Then, the calculation formula for creating the coordinate grid is changed. FIG. 15 exemplifies a detailed enlarged view of the coordinate grid created by emphasizing the smoothness of the coordinate grid and a screen for displaying contour lines of the evaluation function of the smoothness. In this figure, since the coordinate grid points are concentrated on the curve GJ, the curve GJ
It is shown that the contour line becomes denser as it approaches J, and it deviates slightly from the smooth state. From this figure, a guide can be obtained to confirm whether distortion or non-uniformity of the coordinate grid is acceptable.

The flow of processing when the menu of the concentration degree of the coordinate grid is selected is the same as in FIG. First, in the case of displaying, depending on the selection input of whether or not to display the concentration degree, an input of a display range is received, and an evaluation function (eg,
Contour lines of Iv) in equation (5) are displayed on the screen. Next, by selecting whether or not to change the optimization parameter for controlling the degree of concentration, when changing it, the optimization parameter for controlling the degree of concentration, the center of concentration, the range, etc. are received and a coordinate grid is created. Change the formula for the calculation. FIG. 16 exemplifies a detailed enlarged view of the coordinate grid and a screen for displaying the contour lines of the degree of concentration for a calculation example in which the coordinate grid points are concentrated on the fixed curve EH and the curve GJ. By concentrating the coordinate grids, it is possible to efficiently distribute the grid points when the gradient of the physical quantity is greatly different in the analysis region as described above, and it is possible to improve the accuracy of the numerical analysis with a limited number of grid points.

The flow of processing when the orthogonality menu of the coordinate grid is selected is the same as in FIG. By selecting whether or not to display the orthogonality evaluation function, when the display is made, the input of the lattice range to be displayed is received, and the orthogonality evaluation function in the unit lattice at each lattice point (for example, (6) The formula calculates Io),
The evaluation function of orthogonality is displayed in contour lines. Next, by selecting whether or not to change the parameter that controls the orthogonality of the lattice, when changing it, enter the optimization parameter that controls the orthogonality, the range that specifies orthogonality, etc. Change the calculation formula for. FIG. 17 exemplifies a detailed enlarged view of a coordinate grid when the number of coordinates is orthogonal to the upper surface of the system of FIG. 17 and a contour line display of the evaluation function of orthogonality. Maintaining the orthogonality of the coordinate grid at the location where the coordinate grid spacing is large or where the variable changes greatly is effective for improving the convergence method and accuracy of the numerical solution.

A specific effect obtained when a coordinate grid is created using the coordinate grid creation support method of the present invention will be briefly described. FIG. 18 is an example of calculating a simple potential field in a two-dimensional circular region where an exact solution exists. 1 in the figure
0a to 10c are coordinate grids used for calculation, 11a to 11c
Are error distributions for exact solutions of the results calculated using the coordinate grids of 10a to 10c, respectively. 10a, 10b
Is a coordinate grid before being optimized by the supporting method and apparatus of the present invention. As is clear from 11a and 11b, the errors are concentrated and large near the boundary, and the maximum errors are 4.2% and 16.8%, respectively. It has become. When the coordinate grid is optimized by the supporting method and apparatus of the present invention from the viewpoints of smoothness, concentration, and orthogonality, a grid of 10c is obtained in the figure, and the error is distributed evenly over the entire region as shown by 11c. , The maximum error is reduced to 0.7%. In this way, the accuracy of the numerical simulation is improved by optimizing the coordinate grid according to the target problem.

In the embodiment of the present invention, an evaluation function concerning the smoothness, the degree of concentration, and the orthogonality of the created coordinate grid is given to the analysis operator as evaluation guidance, and the optimization parameter is changed based on the experience of the analysis operator. I showed an example
It would be unfamiliar if a so-called expert system could be constructed that uses the knowledge engineering technique to store the know-how of the analyst as knowledge in a computer and more specifically indicate the next operation (parameter change) from this knowledge and the evaluation function. It is possible for an analysis operator to accurately generate an optimum coordinate grid according to the purpose of analysis. In this case, the automatic processing procedure of the determination of the evaluation function and the like and the instruction of the next operation is added to the portion indicated by 4g in the processing flow of FIG. 4, for example.

〔The invention's effect〕

As described above, by using the coordinate grid creation support method and apparatus of the present invention to interactively create a coordinate grid, the analysis target can be acquired while incorporating the know-how of the analysis operator.
The coordinate grid can be optimized according to the purpose of analysis. As a result, the working time required for creating the coordinate grid can be reduced, and the accuracy of the numerical simulation can be improved.

[Brief description of drawings]

FIG. 1 is a diagram showing the configuration of a coordinate grid creation support device of the present invention, FIG. 2 is a diagram showing the contents of an arithmetic processing unit constituting the coordinate grid creation support device of the present invention, and FIG. 3 is a coordinate of the present invention. FIG. 4 is a diagram showing the contents of a storage device constituting the grid creation support device, FIG. 4 is a diagram showing the flow of processing of the coordinate grid creation support method of the present invention, and FIG. 5 is a three-dimensional bird cage of the initial coordinate grid of the object to be analyzed. FIG. 6, FIG. 7, FIG. 7 and FIG. 8 are diagrams showing a menu selection screen of the coordinate grid creation support device of the present invention, and FIG. 9 is a diagram showing a flow of processing when the optimization parameter menu is selected,
FIG. 10 is a general view of the initial coordinate grid of the two-dimensional structure MOSFET,
FIG. 11 is a diagram showing the flow of processing when the split / fixed menu is selected, FIG. 12 is a detailed enlarged view of the coordinate grid when the coordinate grid is subdivided, and FIG. 13 is a case where fixed grid points are included. Fig. 14 shows the result of coordinate grid calculation, Fig. 14 shows the processing flow when smoothness, concentration, and orthogonality menu are selected, and Fig. 15 shows a coordinate grid diagram with smoothness emphasized and smoothness. 16 is a contour map of the evaluation function, FIG. 16 is a coordinate grid diagram when the coordinate grid is concentrated, and FIG. 17 is a contour map of the evaluation function of the concentration degree. FIG. 17 is a coordinate grid diagram when the coordinate lines are orthogonal to the boundary. And a contour display diagram of the orthogonality evaluation function, and FIG. 18 is a coordinate grid diagram and an error distribution diagram showing the effect of the supporting method and apparatus of the present invention. 1 ... Coordinate grid creation support device, 2 ... Display device, 3 ...
Image processing device, 4 ... Arithmetic processing device, 4a ... Optimization parameter input setting unit, 4b ... Coordinate grid generation / calculation unit, 4
c ... Evaluation function calculation unit, 4d ... Graphic processing unit, 4e ...
Result and evaluation guidance display unit, 4f ... coordinate grid output unit, 5 ... storage device, 5a ... coordinate grid data storage unit,
5b ... Evaluation function data storage unit, 6 ... Operation panel, 7 ...
Initial coordinate grid creation part, 8 ... Coordinate grid optimization part, 9
... Display screen, 9a ... Coordinate grid diagram, 9b ... Optimization parameter selection menu, 9c ... Drawing control menu, 9d ... Contour plot of evaluation function, 10a-10c ...
... coordinate grid diagram, 11a to 11c ... error distribution diagram.

 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Ichino Tango 1168 Moriyama-cho, Hitachi, Ibaraki Prefecture 1168 Hiritsu Seisakusho Energy Research Co., Ltd. (72) Yoichi Kose 1168 Moriyama-cho, Hitachi, Ibaraki Hitsuritsu (72) Inventor Hiroki Sano 1168 Moriyama-cho, Hitachi, Hitachi, Ibaraki

Claims (6)

[Claims] 1. A coordinate grid creation support method for creating a plurality of coordinate grid points and coordinate grid lines connecting the coordinate grid points in a target analysis area in order to analyze a physical phenomenon using a computer. As the number of coordinate grid points in the analysis area, the position of the coordinate grid points that are fixed as boundaries, the distance between adjacent coordinate grid points, the bending angle of one coordinate grid point at the coordinate grid points, and the two coordinate grid lines Input the parameters that control the angle of intersection at the coordinate grid point, execute the coordinate grid creation calculation based on the input data, and the coordinates of the coordinate grid point with respect to the created coordinate grid. The step of calculating the value of the evaluation function representing the arrangement characteristic of the coordinate grid in the analysis area calculated using the values, and the coordinate grid diagram created in each step and the value of the calculated evaluation function And displaying on the same screen as the steps of modifying the production conditions of the coordinate grid of the display contents based on the steps of re-executing the coordinate grid creation calculation,
It is characterized in that it comprises a step of recalculating the evaluation function, a step of updating and displaying the display content of the evaluation guidance, and a step of repeating such a correction operation until the value of the evaluation function matches the analysis purpose. A method for supporting coordinate grid creation. 2. The evaluation function relating to the arrangement characteristics of the coordinate grid according to claim 1, is represented by the distance between adjacent coordinate grid points and the bending angle of one coordinate grid line on the coordinate grid points. A coordinate grid creation support method characterized by smoothness of coordinate grid lines. 3. The coordinate grid creation according to claim 1, wherein the evaluation function relating to the arrangement characteristics of the coordinate grid is the concentration degree of the coordinate grid points represented by the distance between adjacent coordinate grid points. How to help. 4. The coordinate grid according to claim 1, wherein the evaluation function relating to the arrangement characteristics of the coordinate grid is represented by an angle at which two coordinate grid lines intersect at a coordinate grid point or an inner product of coordinate grid line vectors. A coordinate grid creation support method characterized by the orthogonality of lines. 5. The optimization function according to claim 1, wherein the evaluation function relating to the arrangement characteristics of the coordinate grid is the optimization parameters for the smoothness of the coordinate grid line, the concentration degree of the coordinate grid point, and the orthogonality of the coordinate grid line. A coordinate grid creation support method characterized by being a sum of values multiplied by. 6. A coordinate grid creation support device for creating a plurality of coordinate grid points and a coordinate grid line connecting these coordinate grid points in a target analysis area for analyzing a physical phenomenon. The number of coordinate grid points, the position of the coordinate grid points that are fixed as boundaries, the distance between adjacent coordinate grid points, the bending angle of one coordinate grid line at the coordinate grid points, and the two coordinate grid points are coordinate grid points. Input the parameters that control the angle of intersection and the input device to correct and input them, and execute the coordinate grid creation calculation based on the input data, and the coordinate grid points for the created coordinate grid. The value of the evaluation function representing the arrangement characteristic of the coordinate grid in the analysis region calculated using the coordinate value of is calculated, and the coordinate grid creation calculation is re-executed in accordance with the correction input and evaluated. From an arithmetic processing device for recalculating the number, and a display device for displaying in parallel the coordinate lattice diagram obtained by the creation calculation, the value of the calculated evaluation function, the menu for selecting the parameter and the selection menu for the drawing to be displayed. A coordinate grid creation support device characterized by: