DESIGN METHOD FOR A PNEUMATIC TIRE

Abstract: A design method for a pneumatic tire according to the present invention comprising the steps of determining a tire basic model representing tire cross-sectional configurations including the internal construction and divided into a plurality of elements, an objective function representing a physical value for evaluating a tire performance, a design variable for determining the cross-sectional configurations of the tire, and limiting conditions for limiting the cross-sectional configurations of the tire, obtaining a corrected model of the tire in which the design variable is changed by a unit amount Δr i , calculating values of the objective function and limiting conditions of the corrected model, calculating sensitivity for every design variable and calculating an estimated value of a changed amount of the design variable that maximizes the value of the objective function by using sensitivity so calculated while allowing for the limiting conditions, repeating the above steps until the values of the objective function converge, and determining the configuration of a tire based on the design variable when the values of the objective function converge upon the estimated value.

Description: The present invention relates to a method of designing a pneumatic tire, and more particularly to a method of designing a pneumatic tire which makes it possible to efficiently and easily design and develop the structure and shape of a tire for attaining single-purpose performance or mutually incompatible performances of the tire, as well as to determine the best structure and shape of the tire and design a tire having improved cost performance. Conventional methods of designing a tire have been based on experiential rules obtained through repeated numerical experiments using experiments and computers. For this reason, the number of trial manufactures and tests required for development has been enormously large, resulting in increased development cost and making it impossible to readily reduce a development period. (1) Background Art Concerning the Design of a Tire Carcass Configuration In general, methods of designing a tire carcass configuration are broadly classified into the following two design methods. Presently most tires are designed, developed, and marketed in accordance with either one of the two methods. One method of designing a tire carcass configuration is called a natural equilibrated configuration design method, in which the cross-sectional shape of a carcass (a cross-sectional shape of a tire cut in a plane including the rotational axis of the tire) is determined in such a manner that tension in tire carcass cords becomes constant particularly when the tire is fitted onto a rim and is inflated with air so as to bring the internal pressure of the tire to a standard level of internal pressure. Namely, the cross-sectional shape of the carcass at a time when the tire is vulcanized and cured in a mold is made substantially identical to the cross-sectional shape of the carcass which is under the standard internal pressure, so as to prevent the occurrence of a failure due to internal strains or the like occurring locally between laminates or in laminated layers of various reinforcing layers of the tire when the tire is inflated with air. The other method is called a non-natural equilibrated configuration design method. As a typical method of designing a non-natural equilibrated configuration, a designing method is known which has already been marketed under a trade mark called "R COT" proposed earlier by the present applicant and has been commercially successful. Unlike the above-mentioned method of designing a natural equilibrated configuration, this method of designing a non-natural equilibrated configuration is a designing method which, conversely, controls the tension of the carcass. By making use of the force which tends to shift to a balanced shape owing to inflation with air, the tension is locally controlled so as to obtain a predetermined performance. (2) Background Art Concerning the Design of the Configuration of a Tire Crown Portion The configuration of the tire crown portion is designed by a number of circular arcs in a cross section including the rotational axis of the tire. Conventional methods of determining values of the circular arcs include, among others, (1) a method of determining the same on the basis of the values of the circular arcs and a data base on tire performance; (2) a method wherein a number of molds are fabricated, and the determination is made from data on the trial manufacture and evaluation of a tire fabricated by the molds; and (3) a method wherein the determination is made by conducting numerical experiments a large number of times. (3) Background Art Concerning the Design of a Belt Structure The belt structure is determined by the belt width, the belt angle, the type of cord, the placement density, and the gauge distribution of the inter-belt rubber. Methods of designing a belt structure include, among others, (1) a method wherein the strain at a belt end is measured through experiments and numerical experiments; (2) a method wherein the rigidity of the overall belt is measured; and (3) a method wherein the belt structure is altered, and indoor evaluations and actual-vehicle evaluations are carried out, so as to extract a draft design. For this reason, there has been the problem that the period of evaluation is long. (4) Background Art Concerning the Design of a Pattern Since the pattern design permits a large degree of freedom, in conventional practice, after tires are grooved in accordance with a drafted basic pattern or after a mold is actually fabricated and the tires are trial-manufactured, an actual-vehicle evaluation is made. Problems encountered with the actual vehicle are overcome by fine-correcting the drafted basic pattern before the drafted basic pattern is finished as a final pattern. For this reason, the pattern design generally requires more manhours as compared to the design of the tire configuration and structure. (5) Background Art Concerning Other Designs Problems similar to those of the above-described design methods have also been experienced with the gauge distribution of a tire rubber member (a gauge distribution of a bead filler, a gauge distribution of a rubber chafer, a gauge distribution of side rubber, a gauge distribution of tread rubber, a gauge distribution of tread base rubber, a gauge distribution of inner-surface reinforcing rubber, a gauge distribution of belt-end rubber, and a gauge distribution of belt end rubber). It has been conventional practice to undertake design and development by trial-manufacturing and testing the tire in accordance with these designing methods and repeating trial manufacture and testing until targeted performance is obtained. With the design and development of tires using conventional designing methods, a targeted value is set with respect to a certain performance. Design and development are considered to be completed for the present if this targeted value is cleared. In other words, the conventional design and development are not based on a viewpoint of obtaining the best performance with the given resources. In addition, conventional design and development are aimed neither at designing mutually incompatible performances, nor at determining the best shape and structure. Moreover, according to either designing method, development is carried out by repeating trial manufacture and testing on a trial-and-error basis, so that there have been problems in that development is very inefficient and cost performance is high. Accordingly, it is an object of the present invention to provide a method of designing a pneumatic tire which, when an attempt is made to obtain a certain single performance or mutually incompatible performances, makes it possible to design a best mode for the tire under given conditions, increase efficiency in the design and development of the tire, and remarkably lower the cost performance of the tire. To attain the above object, the present inventors took note of the fact that an "optimization design technique" which is being used in different fields is possibly applicable to the special field of tires, and conducted various studies, and established that technique specifically as a method of designing a pneumatic tire. Specifically, the present invention provides a computer implemented method of designing a pneumatic tire, which comprises the steps of: (a) defining a basic model of a tire representing a cross-sectional shape of the tire including an internal structure, an objective function representing a physical amount for evaluating tire performance, a design variable for determining at least one of the cross-sectional shape of the tire and a tire structure, and a constraint other than the design variable for constraining at east one of the cross-sectional shape of the tire, the tire structure, a physical amount for evaluating performance, and a tire size; (b) varying a value of the design variable while taking the constraint into consideration and carrying out an optimization computation until an optimum value of the objective function is obtained; and (c) designing the tire on the basis of the design variable which gives the optimum value of the objective function. The aforementioned step (b) may comprise: estimating an amount of change of the design variable which gives the optimum value of the objective function while taking the constraint into consideration on the basis of a sensitivity of the objective function, which is a ratio of an amount of change of the objective function to an amount of unit change of the design variable, and a sensitivity of the constraint, which is a ratio of an amount of change of the constraint to an amount of unit change of the design variable; calculating a value of the objective function when the design variable is changed by an amount corresponding to the estimated amount and a value of the constraint when the design variable is changed by the amount corresponding to the estimated amount; and determining the value of the design variable which gives the optimum value of the objective function on the basis of the estimated value and the calculated values while taking the constraint into consideration. In addition, as a result of conducting various studies, the present inventors took note of the fact that the "technique of genetic algorithms" which is being utilized in different fields can be applied to the special field of tires, and the present inventors attempted various studies, and specifically established it as a method of designing a tire. Specifically, step (a) may comprise: determining a group subject to selection, consisting of a plurality of basic models of the tire each representing a cross-sectional shape of the tire including an internal structure thereof; and determining, with respect to the basic models of the tire of the group subject to selection, the objective function representing the physical amount for evaluating tire performance, the design variable for determining the cross-sectional shape of the tire and/or the tire structure, the constraint for constraining at least one of the cross-sectional shape of the tire, the tire structure, the physical amount for evaluating performance, and the tire size, and a fitness function which can be evaluated from the objective function and the constraint; and step (b) may comprise: selecting two basic models of the tire from the group subject to selection on the basis of the fitness function; effecting at least one of generating a new basic model of the tire by allowing design variables of the basic models of the tire to cross over each other at a predetermined probability and of generating a new basic model of the tire by changing (mutating) a portion of the design variable of at least one of the basic models of the tire; determining the objective function, the constraint, and the fitness function of the basic model of the tire with the design variable changed, preserving the basic model of the tire and the basic model of the tire with the design variable not changed, the aforementioned process of step (b) being repeated until the basic models of the tire preserved reaches a predetermined number; determining whether or not a new group consisting of the predetermined number of the preserved basic models of the tire satisfies a predetermined convergence criterion; if the convergence criterion is not satisfied, repeating the foregoing process of step (b) until the group subject to selection satisfies a predetermined convergence criterion by setting the new group as the group subject to selection; and if the predetermined convergence criterion is satisfied, determining the value of the design variable which gives the optimum value of the objective function among the predetermined number of the preserved basic models of the tire while taking the constraint into consideration. In this step (b), with respect to the basic model of the tire with the design variable changed, the amount of change of the design variable which gives the optimum value of the objective function may be estimated while taking the constraint into consideration on the basis of the sensitivity of the objective function, which is the ratio of the amount of change of the objective function to the amount of unit change of the design variable, and on the basis of the sensitivity of the constraint, which is the ratio of the amount of change of the constraint to the amount of unit change of the design variable, the value of the objective function when the design variable is changed by the amount corresponding to the estimated amount and the value of the constraint when the design variable is changed by the amount corresponding to the estimated amount may be calculated, the fitness function may be determined from the value of the objective function and the value of the constraint, the basic model of the tire and the basic model of the tire with the design variable not changed may be preserved, and the aforementioned process of step (b) may be repeated until the preserved basic models of the tire reaches the predetermined number. The design variable may include at least one of: a function representing the shape of at least one line selected from a carcass line, a turn-up ply line, a line representing an outer configuration of the tire, a line representing a tire crown portion, and a reinforcing-material line; a variable representing at least one gauge distribution of a tire rubber member selected from a gauge distribution of a bead filler, a gauge distribution of a rubber chafer, a gauge distribution of side rubber, a gauge distribution of tread rubber, a gauge distribution of tread base rubber, a gauge distribution of inner-surface reinforcing rubber, a gauge distribution of inter-belt rubber, and a gauge distribution of belt end rubber; a variable representing at least one structure of a belt portion selected from an angle of each belt layer, a width thereof, a cord type thereof, and a placing density thereof; and a variable representing at least one configuration of a pattern selected from a configuration of a block, a position of a sipe, a number of sipes, and a length of the sipe. A plurality of lines may be set as the design variable by representing the shape of at least one line selected from a carcass line, a turn-up ply line, a line representing an outer configuration of the tire, a line representing a tire crown portion, and a reinforcing-material line, and, in order to prevent the plurality of lines from intersecting with each other when the design variable is varied in step (b), a mapping function of a function representing shapes of the lines and a constraint for preventing the lines from intersecting with each other may be determined, and the value of the design variable which gives the optimum value of the objective function while satisfying the constraint may be determined by executing step (b) in a space mapped by the mapping function. In step (a) of the present invention, a basic model of a tire representing a cross-sectional shape of the tire including an internal structure, an objective function representing a physical amount for evaluating tire performance, a design variable for determining at least one of the cross-sectional shape of the tire and a tire structure, and a constraint other than the design variable for constraining at least one of the cross-sectional shape of the tire, the tire structure, a physical amount for evaluating performance, and a tire size are defined. In addition to a line representing the outer configuration of the tire, this basic model of the tire may include a line representing the configuration of a tire crown, a belt line representing a belt in the interior of the tire, a carcass line representing the carcass of the tire, a turn-up ply line representing a turn-up line of the carcass ply in the interior of the tire, and a reinforcing-material line representing a line of each of various reinforcing materials, the gauge distribution of a tire rubber member, the angle, width, cord type, and placement density of each belt layer representing the structure of the belt portion, as well as the configuration of a block, the position of a sipe, the number of sipes, and the length of the sipe which represent the configuration of the pattern. In addition, a technique called the finite element method or an analytical technique may be used for the basic model of the tire for dividing the basic model of the tire into a plurality of elements. As the objective function representing a physical amount for evaluating tire performance, it is possible to use a physical amount for governing the relative merit of the tire performance, including the lateral spring constant, belt tension in the circumferential direction of the tire when inflated to improve steering stability, or the like. As the design variable for determining the cross-sectional shape of the tire, it is possible to use, among others, a function representing at least one line selected from the carcass line, the turn-up ply line, the outer configuration of the tire, the line representing the configuration of the tire crown, and the reinforcing-material line representing the line of each reinforcing material. As the design variable for determining the tire structure, it is possible to use a variable representing at least one gauge distribution of a tire rubber member selected from a gauge distribution of a bead filler, a gauge distribution of a rubber chafer, a gauge distribution of side rubber, a gauge distribution of tread rubber, a gauge distribution of tread base rubber, a gauge distribution of inner-surface reinforcing rubber, a gauge distribution of inter-belt rubber, and a gauge distribution of belt end rubber, as well as a variable representing the structure of the belt portion, the bead and side portions, including an angle of a belt layer, the width of the belt layer, the height of the ply, an amount of the ply turned up, and the angle, width, position, material and the like of a bead-portion reinforcing material. In addition, as the design variable for determining the pattern, it is possible to use a variable representing the configuration of the pattern, including the configuration of a block, the block groove-wall angle, the position of a sipe, the number of sipes, and the length of the sipe. As the constraint for constraining the cross-sectional configuration of the tire and the tire structure, it is possible to cite, for example, a constraint in a periphery value of the carcass line, a constraint in the vertical primary eigenfrequency, a constraint in the angle of the belt layer, a constraint in the width of the belt layer, and constraints such as the tire size, spring constants, the amount of deformation of the tire, tire weight, stress, strain, strain energy, and rolling resistance. It should be noted that the objective function, the design variable, and the constraint are not confined to the aforementioned examples, and various items may be set as such depending on the objective of tire design. In the ensuing step (b), there is varied a value of the design variable while taking the constraint into consideration and there is carried out an optimization computation until an optimum value of the objective function is obtained. In this case, it is effective if an amount of change of the design variable which gives the optimum value of the objective function while taking the constraint into consideration is estimated on the basis of a sensitivity of the objective function, which is a ratio of an amount of change of the objective function to an amount of unit change of the design variable, and a sensitivity of the constraint, which is a ratio of an amount of change of the constraint to an amount of unit change of the design variable, a value of the objective function when the design variable is changed by an amount corresponding to the estimated amount and a value of the constraint when the design variable is changed by the amount corresponding to the estimated amount are calculated, and the value of the design variable which gives the optimum value of the objective function while taking the constraint into consideration is determined on the basis of the estimated value and the calculated values. By so doing, the value of the design variable which gives an optimum value of the objective function while taking the constraint into consideration can be obtained. In addition, in step (a), it is effective to determine a group subject to selection, consisting of a plurality of basic models of the tire each representing a cross-sectional shape of the tire including an internal structure thereof, and determine, with respect to the basic models of the tire of the group subject to selection, the objective function representing the physical amount for evaluating tire performance, the design variable for determining the cross-sectional shape of the tire and/or the tire structure, the constraint for constraining at least one of the cross-sectional shape of the tire, the tire structure, the physical amount for evaluating performance, and the tire size, and a fitness function which can be evaluated from the objective function and the constraint. In step (b), it is effective to carry out the following: selecting two basic models of the tire from the group subject to selection on the basis of the fitness function; effecting at least one of generating a new basic model of the tire by allowing design variables of the basic models of the tire to cross over each other at a predetermined probability and of generating a new basic model of the tire by changing a portion of the design variable of at least one of the basic models of the tire; determining the objective function, the constraint, and the fitness function of the basic model of the tire with the design variable changed; preserving the basic model of the tire and the basic model of the tire with the design variable not changed, the aforementioned process of step (b) being repeated until the basic models of the tire preserved reaches a predetermined number; determining whether or not a new group consisting of the predetermined number of the preserved basic models of the tire satisfies a predetermined convergence criterion; if the convergence criterion is not satisfied, repeating the foregoing process of step (b) until the group subject to selection satisfies a predetermined convergence criterion by setting the new group as the group subject to selection; and if the predetermined convergence criterion is satisfied, determining the value of the design variable which gives the optimum value of the objective function among the predetermined number of the preserved basic models of the tire while taking the constraint into consideration. In this case, in step (b), with respect to the basic model of the tire with the design variable changed, it is further effective to carry out as follows: The amount of change of the design variable which gives the optimum value of the objective function is estimated while taking the constraint into consideration on the basis of the sensitivity of the objective function, which is the ratio of the amount of change of the objective function to the amount of unit change of the design variable, and on the basis of the sensitivity of the constraint, which is the ratio of the amount of change of the constraint to the amount of unit change of the design variable, the value of the objective function when the design variable is changed by the amount corresponding to the estimated amount and the value of the constraint when the design variable is changed by the amount corresponding to the estimated amount are calculated, the fitness function is determined from the value of the objective function and the value of the constraint, the basic model of the tire and the basic model of the tire with the design variable not changed are preserved, and the aforementioned process of step (b) is repeated until the preserved basic models of the tire reaches the predetermined number. Through this measure as well, it is possible to obtain the value of the design variable for optimizing the value of the objective function by taking the constraint into account. It should be noted that, as the fitness function which can be evaluated from the objective function and the constraint, it is possible to use a function for determining adaptability with respect to the tire model from the objective function and the constraint. In addition, the objective function, the design variable, the constraint, and the fitness function are not confined to the above examples, and various ones may be set in accordance with the object of the tire design. Furthermore, in the crossover of the design variables of the basic models of the tire, there is a method wherein a part of the design variable or design variables at or subsequent to a predetermined position are exchanged with respect to the design variables of two tire models selected. Furthermore, in the change of a part of the design variable of the tire model, there is a method wherein the design variable at a position determined at a predetermined probability is changed (mutated). In a case where a plurality of lines are set as the design variable by means of a function representing the shape of at least one line selected from a carcass line, a turn-up ply line, a line representing an outer configuration of the tire, a line representing a tire crown portion, and a reinforcing-material line, if a mapping function of a function representing shapes of the lines and a constraint for preventing the lines from intersecting with each other are determined in order to prevent the plurality of lines from intersecting with each other when the design variable is varied, and if the value of the design variable which gives the optimum value of the objective function while satisfying the constraint is determined by executing step (b) in a space mapped by the mapping function, then the plurality of lines cease to intersect each other. Then in step (c), the tire is designed by changing the basic model of the tire and the like on the basis of the design variable which gives the optimum value of the objective function. In addition, in a case where mapping functions are used, the tire is designed on the basis of the design variable mapped on the original function. As described above, in step (a), in a case where a group subject to selection, consisting of a plurality of basic models of the tire each representing a cross-sectional shape of the tire including an internal structure thereof is determined, and in a case where the objective function representing the physical amount for evaluating tire performance, the design variable for determining the cross-sectional shape of the tire and/or the tire structure, the constraint for constraining at least one of the cross-sectional shape of the tire, the tire structure, the physical amount for evaluating performance, and the tire size, and a fitness function which can be evaluated from the objective function and the constraint are determined with respect to the basic models of the tire of the group subject to selection, step (b) can be comprised of the following steps (1) to (4): (1) the step of selecting two basic models of the tire from the group subject to selection on the basis of the fitness function; (2) the step of effecting at least one of generating a new basic model of the tire by allowing design variables of the tire models to cross over each other at a predetermined probability and of generating a new basic model of the tire by changing a portion of the design variable of at least one of the basic models of the tire; (3) the step of determining the objective function, the constraint, and the fitness function of the basic model of the tire with the design variable changed by crossover or change, preserving the basic model of the tire and the basic model of the tire with the design variable not changed, and repeating the steps (1) to (3) until the basic models of the tire preserved reaches a predetermined number; and (4) the step of determining whether or not a new group consisting of the predetermined number of the preserved basic models of the tire satisfies a predetermined constraint; if the constraint is not satisfied, repeating the steps (1) to (4) until the group subject to selection satisfies a predetermined constraint by setting the new group as the group subject to selection; and, if the predetermined constraint is satisfied, determining the value of the design variable which gives the optimum value of the objective function among the predetermined number of the preserved basic models of the tire while taking the constraint into consideration. In this step (3), it is effective to carry out the following with respect to the basic model of the tire with the design variable changed: The amount of change of the design variable which gives the optimum value of the objective function is estimated while taking the constraint into consideration on the basis of the sensitivity of the objective function, which is the ratio of the amount of change of the objective function to the amount of unit change of the design variable, and on the basis of the sensitivity of the constraint, which is the ratio of the amount of change of the constraint to the amount of unit change of the design variable. Then, the value of the objective function when the design variable is changed by the amount corresponding to the estimated amount and the value of the constraint when the design variable is changed by the amount corresponding to the estimated amount are calculated. Furthermore, the fitness function is calculated from the value of the objective function and the value of the constraint, and the new basic model of the tire with the design variable changed and the basic model of the tire with the design variable not changed are preserved. The steps (1) to (3) are repeatedly executed until the preserved basic models of the tire reaches the predetermined number. By virtue of this measure as well, it is possible to determine the value of the design variable which optimizes the value of the objective function by taking into the constraint into consideration. If design and development are effected in accordance with the designing method of the present invention, unlike the conventional design and development based on trial and error, the process ranging from the design of a best mode to the performance evaluation of the designed tire becomes possible to a certain extent based principally on computer calculation. Hence, it is possible to increase the efficiency remarkably, thereby permitting a reduction in the cost required in development. As described above, in accordance with the present invention, since the design variable which gives an optimum value of the objective function while satisfying the constraint is determined, and the tire is designed on the basis of this design variable, it is possible to obtain an advantage in that a best mode for the tire can be designed. The invention will be further described, by way of example only, with reference to the accompanying drawings, in which: Fig. 1 is a schematic diagram of a personal computer used in the embodiments of the present invention; Fig. 2 is a flowchart illustrating a processing routine in accordance with a first embodiment of the present invention; Fig. 3 is a flowchart illustrating a routine for determining a design variable shown in Fig. 2; Fig. 4 is a diagram illustrating a basic model of a tire; Fig. 5 is a diagram illustrating a state in which phantom lines passing through a reference point P are drawn for each d on the basic model of the tire shown in Fig. 4; Fig. 6 is a diagram illustrating relationships among selected nodes, a distance r i to each of these nodes, and an ascending angle  i ; Fig. 7 is a flowchart illustrating a processing routine in accordance with a second embodiment of the present invention; Fig. 8 is a flowchart illustrating a routine for determining the design variable shown in Fig. 7; Fig. 9 is a diagram in which a starting point S, an ending point E, a point of inflection C, and a point of a maximum-width position M are shown in a basic model of a tire; Fig. 10 is a diagram illustrating radii of curvature of circular arcs assumed between respective points, i.e., the starting point S, the ending point E, the point of inflection C, and the point of a maximum-width position M; Fig. 11 is a diagram illustrating an amount of offset between an original node in a basic model of a tire and an optimal point determined by a method of least squares; Fig. 12 is a diagram illustrating r-z coordinates; Fig. 13 is a diagram illustrating ξ 1 -η 1 coordinates; Fig. 14 is a diagram illustrating ξ 2 -η 2 coordinates; Fig. 15 is a diagram illustrating configurations of a bead filler of a tire and its peripheral components; Fig. 16 is a diagram illustrating lines of belt layers; Fig. 17 is a diagram illustrating the approximation of the lines of the belt layers by means of curves; Fig. 18 is a diagram illustrating the distortion of the shapes of elements due to the movement of nodes; Fig. 19 is a flowchart illustrating a routine for processing division into elements; Figs. 20A to 20D are diagrams illustrating the processing of division into elements, wherein Fig. 20A is a diagram illustrating a model at the time of defining regions, Fig. 20B is a diagram illustrating the nodes generated in the respective regions, Fig. 20C is a diagram illustrating the generation of triangular elements, and Fig. 20D is a diagram illustrating that the smoothing of the nodes has been carried out; Fig. 21 is a flowchart illustrating a Lagrange interpolation processing routine; Fig. 22 is a diagram illustrating a crown portion being represented by curves; Fig. 23 is a diagram illustrating the crown portion being divided by the number of the degree of the Lagrange interpolation; Fig. 24 is a diagram illustrating lines normal to the line of the crown portion; Fig. 25 is a diagram illustrating the configuration of the crown portion being represented by circular arcs; Fig. 26 is a diagram illustrating the configuration of the crown portion being represented by circular arcs due to independent central coordinates; Fig. 27 is a diagram illustrating a pattern surface being divided into the form of a lattice in accordance with the degree of the Lagrange interpolation; Fig. 28 is a diagram illustrating the subdivision of the pattern surface; Fig. 29 is a flowchart illustrating a processing routine in accordance with a tenth embodiment of the present invention; Fig. 30 is a flowchart illustrating an intersection processing routine; Figs. 31A and 31B are diagrams illustrating chevron-shaped mapping functions, wherein Fig. 31A is a diagram illustrating a continuous chevron-shaped mapping function, and Fig. 31B is a diagram illustrating a linear chevron-shaped mapping function; Figs. 32A and 32B are diagrams illustrating trough-shaped mapping functions, wherein Fig. 32A is a diagram illustrating a continuous trough-shaped mapping function, and Fig. 32B is a diagram illustrating a linear trough-shaped mapping function; Fig. 33 is a flowchart illustrating a mutation processing routine; Figs. 34A and 34B are diagrams illustrating a basic model of a tire in accordance with an eleventh embodiment, wherein Fig. 34A is a diagram illustrating a basic model of the tire during fabrication, and Fig. 34B is an image diagram explaining the presence of the belts; Fig. 35 is a flowchart illustrating an intersection processing routine; Fig. 36 is a flowchart illustrating a mutation processing routine; Fig. 37 is a flowchart illustrating a processing routine in accordance with a twelfth embodiment of the present invention; Figs. 38A and 38B are diagrams illustrating the modeling of a block, wherein Fig. 38A is a diagram illustrating the configuration of the block for modeling, and Fig. 38B is a diagram illustrating a modeling example of a finite element model; Fig. 39 is a flowchart illustrating a processing routine in accordance with a thirteenth embodiment of the present invention; Fig. 40 is a flowchart illustrating a processing routine in accordance with a fourteenth embodiment of the present invention; and Fig. 41 is a diagram illustrating heel-and-toe wear. Fig. 1 shows an outline of a personal computer for carrying out the method of designing a pneumatic tire in accordance with the present invention. This personal computer comprises a keyboard 10 for entering data and the like; a main unit 12 for calculating design variables for satisfying a constraint in accordance with a program stored in advance, and for optimizing, e.g., maximizing or minimizing, an objective function; and a CRT 14 for displaying such as the results of calculation by the main unit 12 of the computer. Next, a description will be given of a first embodiment for designing the configuration of a carcass line for setting the belt tension in the circumferential direction of a tire when inflated with air to a maximum value, i.e., an optimum value, so as to improve steering stability. Fig. 2 shows a processing routine of a program in accordance with this first embodiment. In Step 100, the cross-sectional shape of the tire in a state of natural equilibrium is set as a reference shape. This reference shape is subjected to modeling by a technique, such as a finite element method, which makes it possible to numerically or analytically determine the belt tension in the circumferential direction of the tire when inflated with air, so as to determine a basic model of the tire which represents a cross-sectional shape of the tire including its internal structure, and in which the cross-sectional shape is divided into a plurality of elements by meshing. It should be noted that the basic shape is not confined to the cross-sectional shape of the tire in a state of natural equilibrium, but may be an arbitrary shape. The term modeling referred to herein means the numerical representation of the shape, structure, materials, and pattern of the tire into a format of input data for a computer program prepared according to a numerical and analytical technique. Fig. 4 shows a basic model of the tire, in which CL denotes a carcass line; OL, a line representing an outer configuration of the tire; PL, a turn-up ply line; and B1 and B2 denote lines indicating belts. In addition, this basic model of the tire is divided into a plurality of elements by a plurality of lines NL 1 , NL 2 , NL 3 , ... normal to the carcass line CL. It should be noted that although, in the above, a description has been given of an example in which the basic model of the tire is divided into a plurality of elements by a plurality of lines normal to the carcass line, the basic model of the tire may be divided into a plurality of elements by a plurality of lines normal to a line representing an outer configuration of the tire or by a plurality of lines normal to a turn-up ply line, or may be divided into arbitrary shapes such as triangles depending on a design objective. In an ensuing Step 102, an objective function representing a physical amount for evaluating tire performance, a constraint for constraining the cross-sectional shape of the tire, and a design variable determining the cross-sectional shape of the tire are determined. In this embodiment, the objective function OBJ and the constraint G are defined as follows in order to design the configuration of the carcass line for maximizing the belt tension in the circumferential direction of the tire when inflated with air: Objective function OBJ: the sum total, in the widthwise direction of the tire, of components of the belt tension acting in the circumferential direction of the tire Constraint G: A periphery value of the carcass line is within ±5% of the periphery value of the carcass line of the basic model of the tire. It should be noted that the aforementioned periphery value of the carcass line can be calculated as the sum total of distances between nodes (points of intersection between the carcass line and the respective normal line) in the carcass line which are located in a domain that changes the tire shape. In addition, the configuration of the carcass line, which is a design variable, is determined by a Lagrange interpolation routine, shown in Fig. 3, for approximating a curve. In Step 130 for this Lagrange interpolation routine, a reference point P is set in the interior of the tire, as shown in Fig. 4. In an ensuing Step 132, the domain from a node q 1 in the vicinity of a belt end to a node q 2 in a vicinity constrained by the rim is specified as the domain for changing the shape of the tire. In Step 134, by using as a reference line a straight line connecting the node q 1 and the reference point P, an ascending angle , which is an angle formed by this reference line and the straight line connecting the node q 2 and the reference point P, is calculated. In Step 136, an angular increment d is calculated in accordance with the following formula: d =  / the order of Lagrange interpolation It should be noted that the order of Lagrange interpolation is inputted in advance by the user. In Step 138, as shown in Fig. 5, phantom lines L 1 , L 2 , L 3 , ... are assumed for respective angular increments d with the reference line set as a reference, and normal lines nl 1 , nl 2 , nl 3 , ... which are closest to the respective phantom lines are selected. In an ensuing Step 140, as shown in Fig. 6, distances r 1 , r 2 , r 3 , ... (hereafter expressed as r i in a general formula, where i = 1, 2, 3, ...; the order of Lagrange interpolation-1) between the reference point P and innermost nodes Q 1 , Q 2 , Q 3 , ... on the selected normal lines nl 1 , nl 2 , nl 3 , as well as ascending angles  1 ,  2 ,  3 , ... (hereafter expressed as  i in a general formula, where i = 1, 2, 3, ...; the order of Lagrange interpolation-1) at the nodes Q 1 , Q 2 , Q 3 , ... are calculated. Then, in Step 142, the distance r i is set as a design variable. After the objective function OBJ, the constraint G, and the design variable r i are thus determined, in Step 104 in Fig. 2, an initial value OBJo of the objective function OBJ and an initial value Go of the constraint G when the design variable r i is at an initial value ro are calculated. In an ensuing Step 106, the design variable r i is changed by Δr i each so as to change the basic model of the tire. In Step 108, a distance r m between the reference point P and an innermost node other than an innermost node Q i corresponding to the design variable is calculated in accordance with the following formulae: where  m is an ascending angle formed between the reference line and a straight line connecting that node and the reference point P. In addition, in Step 108, a distance between the innermost node Q i and a node on the normal line other than the innermost node (a node on the carcass line), i.e., each node on the line OL representing the outer configuration of the tire, on the turn-up ply line PL, and on the lines B1, B2 representing the belts, is determined. The distance thus determined is added to the coordinates at the innermost node Q i so as to determine the coordinates of the node other than the innermost node on the normal line. Thus, the cross-sectional shape of the tire after the design variable is changed by Δr i , i.e., a corrected model of the tire, is determined. In Step 110, a value OBJ i of the objective function and a value G i of the constraint after the design variable is changed by Δr i are calculated with respect to the corrected model of the tire determined in Step 108. In Step 112, a sensitivity dOBJ/dr i of the objective function, which is a ratio of an amount of change of the objective function to an amount of unit change of the design variable, as well as a sensitivity dG/dr i of the constraint, which is a ratio of an amount of change of the constraint to an amount of unit change of the design variable, are calculated for each design variable in accordance with the following formulae: dOBJ dr i = OBJ i - OBJ o Δr i dG dr i = G i - G o Δr i By means of these sensitivities, it is possible to estimate to what extent the value of the objective function and the value of the constraint change when the design variable is changed by Δr i . It should be noted that these sensitivities are sometimes determined analytically depending on the technique used in modeling the tire and the nature of the design variable, in which case the calculation of Step 110 becomes unnecessary. In an ensuing Step 114, an amount of change of the design variable which maximizes the objective function while satisfying the constraint is estimated by means of mathematical programming by using the initial value OBJo of the objective function, the initial value Go of the constraint, the initial value ro of the design variable, and the sensitivities. In Step 115, by using this estimated value of the design variable, a corrected model of the tire is determined in a method similar to that of Step 108, and a value of the objective function is calculated. In Step 116, a determination is made as to whether or not the value of the objective function has converged by comparing a threshold value inputted in advance with the difference between the value OBJ of the objective function calculated in Step 115 and the initial value OBJo of the objective function calculated in Step 104. If the value of the objective function has not converged, Steps 104 through 116 are executed repeatedly by setting as the initial value the value of the design variable determined in Step 114. If it is determined that the value of the objective function has converged, the value of the design variable at this juncture is set as the value of the design variable which maximizes the objective function while satisfying the constraint, and the shape of the tire is determined in Step 120 by using this value of the design variable. Next, a description will be given of a second embodiment of the present invention. In this embodiment, design variables different from the design variable of the first embodiment are used. Fig. 7 shows a processing routine of a program in accordance with the second embodiment, in which steps corresponding to those of Fig. 2 will be denoted by the same numerals. In Step 160, the objective function and the constraint are determined in the same manner as described above, and the design variables are determined in Step 162. Fig. 8 shows the details of Step 162, wherein in Step 150, a starting point S, an ending point E, a point of inflection C, and a point of a maximum-width position M are selected from nodes on the carcass line in the domain for changing the shape of the basic model of the tire, and data on these points are inputted, as shown in Fig. 9. In an ensuing Step 152, as shown in Fig. 10, one or a plurality of circular arcs are assumed between the respective points, i.e., the starting point S, the ending point E, the point of inflection C, and the point of a maximum-width position M, and radii of curvature R 1 , R 2 , R 3 , and R 4 of these circular arcs and ascending angles  1 ,  2 ,  3 ,  4 are calculated. It should be noted that if it is desired to increase the number of the design variables, it suffices if the number of the circular arcs is increased. In an ensuing Step 154, the coordinates of a point of intersection between the circular arc and the normal line are calculated, and, as shown in Fig. 11, an amount of offset ε (expressed as ε i in a general formula) between the node on the carcass line and the point of intersection between the circular arc and the normal line is calculated. In Step 156, a radii of curvature R i(opt) of the circular arc and an ascending angle  i(opt) in which Σε i 2 becomes minimum are calculated by mathematical programming, and these values R i(opt) and  i(opt) are set as design variables R i and  i . It should be noted that, although not shown in the drawing, after the objective function, the constraint, and the design variables are determined, their initial values are calculated in the same manner as described above. In Step 164 in Fig. 7, the node on the carcass line is moved to an optimal point which is defined by R i(opt) and  i(opt) determined in Step 156. In addition, in an ensuing Step 166, a distance between the node on the carcass line and a node other than the node on the carcass line, i.e., each node on the line OL representing the outer configuration of the tire, on the turn-up ply line PL, and on the lines B1, B2 representing the belts, is determined. The distance thus determined is added to the coordinates at the node on the carcass line so as to determine the coordinates of the node other than the node on the carcass line. Thus, the cross-sectional shape of the tire when each node on the carcass line is moved to the optimal point is set as a basic model of the tire. Thereafter, in the same manner as in Steps 104 through 106, the following processing is carried out: The initial value OBJo of the objective function and the initial value Go of the constraint when the design variables are at initial values Ro and o are calculated. The sensitivity is calculated for each design variable by using the corrected model of the tire after the design variables are changed by ΔR i and Δ i . An amount of change of the design variable which maximizes the objective function while satisfying the constraint is estimated by means of mathematical programming. A corrected model of the tire is determined at the estimated value of each of the design variables, and the objective function OBJ is calculated. A determination is made as to whether or not the value of the objective function has converged by comparing a threshold value with the difference with the initial value OBJo. If the value of the objective function has not converged, Steps 104 through 106 are executed repeatedly by setting as initial values the values of the design variables determined in Step 114. If it is determined that the value of the objective function has converged, the shape of the tire is determined in the same manner as described above. It should be noted that although in the above-described embodiment the shape of the carcass line is adopted as the design variable, the shape of the turn-up ply line, the line representing the outer configuration of the tire, the lines of reinforcing members, and the like may be determined as the design variables. Next, a description will be given of a third embodiment of the present invention. In this embodiment, the present invention is used in the determination of the shapes of a plurality of lines, and is designed to optimize mutually incompatible performances. In this embodiment, the shape of the carcass line, the shape of the turn-up ply line, and the outer configuration of the tire are determined for improving steering stability without impairing the riding comfort of the occupant in the vehicle. If the riding comfort is made satisfactory, the steering stability deteriorates, so that the riding comfort of the occupant in the vehicle and the steering stability are incompatible. The procedure of this embodiment is substantially similar to the above-described embodiments, but the selection of the objective function, the constraint, and the design variables differs. In this embodiment, lateral stiffness, i.e., a physical amount for improving the steering stability, is used as the objective function. As the constraint, a condition that a vertical 1st order natural frequency, i.e., a physical amount for governing the riding comfort, is fixed is used so as to determine the shape of the carcass line, the shape of the turn-up ply line, and the outer configuration of the tire which maximize the lateral stiffness under the condition of the vertical 1st order natural frequency being fixed. As for the design variables, if the methods of the above-described embodiments are applied as they are to three kinds of lines including the carcass line, the turn-up ply line, and the line representing the outer configuration of the tire, there arises the problem that the lines intersect with each other in the process of optimization. Therefore, the design variables are determined as follows: With respect to the carcass line, the design variable r i is determined by using the Lagrange interpolation method explained in the first embodiment. With respect to the turn-up ply line and the line representing the outer configuration of the tire, the design variables are determined as follows: First, as shown in Fig. 12, a region RE 1 having the carcass line as one side is assumed in the original r-z coordinates. Then, a mapping function in which this region RE 1 is formed into a square having the carcass line as one side is prepared in ξ 1 -η 1 coordinates which are different from the r-z coordinates. Further, lines representing the turn-up ply line and the outer configuration of the tire are mapped. Also, a region RE 2 having the turn-up ply line as one side is similarly assumed in the ξ 1 -η 1 coordinates. A mapping function in which a square having the turn-up ply line as one side is then prepared in ξ 2 -η 2 coordinates which are different from the ξ 1 -η 1 coordinates. Further, a line representing the outer configuration of the tire is mapped. Here, the mapping functions are expressed as follows: r = i N i n (ξ 1 , η 1 )r i z = i N i n (ξ 1 , η 1 )z i ξ 1 = i N i n (ξ 2 , η 2 )ξ 1i η 1 = i N i n (ξ 2 , η 2 )n 1i Then, the turn-up ply line in the ξ 1 -η 1 coordinates is determined as a design variable, and the line representing the outer configuration of the tire in the ξ 2 -η 2 coordinates is determined as another design variable. In addition, a constraint that ξ 1 , ξ 2 > -1 is added. This prevents the lines from intersecting with each other. The results of tests conducted by actually trial-manufacturing the tires obtained in the first and third embodiments are shown in Table 1 below. Although in the above example the three lines are adopted as design variables, two lines or four or more lines may be adopted as the design variables. Conventional Tire 1st Embodiment 3rd Embodiment Time required for design and development 100 50 60 Performance Steering stability 100 106 104 Riding comfort 100 98 100 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 Testing method: feeling evaluation by an actual-vehicle running test Next, a description will be given of a fourth embodiment of the present invention. In this embodiment, the determination of the shapes of a plurality of lines illustrated in the third embodiment is applied to the determination of the shapes of a tire bead filler and its peripheral members. In this embodiment, the shape of the bead filler and the shape of a rubber chafer are determined for alleviating the rolling resistance without impairing the durability of a bead portion. Although the procedure of this embodiment is substantially similar to that of the third embodiment, the selection of the objective function, the constraint, and the design variables differs. First, the objective function and the constraint are determined as follows: Objective function: rolling resistance value Constraint: A principal strain occurring at a turn-up ply end under a load is within +3% of the initial structure. In addition, as shown in Fig. 15, design variables are lines BLa, BLb defining an outer configuration (gauge distribution) of a bead filler 20 and a boundary line KL between a rubber chafer 22 and side rubber 24. It should be noted that, if the aforementioned lines are applied as they are with respect to the design variables in the same way as in the third embodiment, there arises the problem that the lines intersect with each other in the process of optimization. Therefore, the design variables may be determined by using such as the Lagrange interpolation method and the circular interpolation method described in the foregoing embodiments. In addition, a mapping function may be prepared in which a region assumed in predetermined coordinates is formed into a square or the like in different coordinates, so as to map the line. Then, a description will be given of a fifth embodiment of the present invention. In this embodiment, the determination of the shapes of a plurality of lines illustrated in the third embodiment is applied to the determination of the thickness of belt layers in a belt portion. In this embodiment, the gauge distribution of the belt layers are determined for alleviating the rolling resistance without impairing the durability of the belt portion. Although the procedure of this embodiment is substantially similar to that of the third embodiment, the selection of the objective function, the constraint, and the design variables differs. First, the objective function and the constraint are determined as follows: Objective function: rolling resistance value Constraint: Main strains occurring at a ply end and a belt end under a load are within +3% of the initial structure. In addition, as shown in Fig. 16, the design variables are the lines B1, B2, B3, and B4 representing the respective belt layers. As these lines B1 to B4 are changed, the gauge distribution of the belt layers for optimizing the objective function is determined. In the foregoing embodiments, the Lagrange interpolation method and the circular interpolation method are used as methods of representing the shapes of lines. In addition to these methods, interpolation may be effected by using spline curves, B-spline curves, Bezier curves, NURBS (non-uniform rational B-splines), or the like which are shown in Mathematical Elements for Computer Graphics (by David F. Rogers and J. Alan Adams). Tires were actually trial-manufactured by combining the above-described fourth and fifth embodiments and were tested, and the results of the test are shown in Table 2 below. Conventional Tire 4th/5th Embodiment Time required for design and development 100 55 Rolling resistance value 100 93 Main strain at a ply end 100 102 Main strain at a belt end 100 97 Tire size: 165SR13 Internal pressure: 1.9 kg/cm 2 Load: 425 kgf Method of testing the rolling resistance value: measurement of a resistance value at a speed of 50 km/h during a drum running test under an inertia Next, a description will be given of a sixth embodiment of the present invention. In this embodiment, the present invention is applied to a belt structure, and the durability is improved by minimizing a strain concentration occurring in the belt portion without increasing the weight of the belt portion. Although the procedure of this embodiment is substantially similar to that of the first embodiment shown in Fig. 2, only the portion of the selection of the objective function, the constraint, and the design variables in Step 102 differs. First, the objective function and the constraint are determined as follows: Objective function: a maximum value of a main strain occurring between the respective belt layers under a load Constraint: The total weight of the belt is within +1% of the total weight of the initial structure. In addition, the design variables are the angle, placement (e.g., the number of placements and placing strength), and width of each belt layer. Of these design variables, as for the angle and placement, it suffices if their values are directly set as the design variables, and, as for the width, the design variable is determined as follows. First, when a basic model of the tire such as the one shown in Fig. 4 is prepared, the line of each belt is approximated by means of, for instance, polynomials by using coordinates of nodes constituting each belt layer. For instance, the lines B1 and B2 are approximated by polynomials, as shown in Fig. 17. In this case, the line B2 is approximated by r = a 1 z 3 + b 1 z 2 + c 1 z 1 + d 1 while the line B1 is approximated by r = a 2 z 3 + b 2 z 2 + c 2 z 1 + d 2 In the step of determining the design variables that give an optimum value of the objective function, the nodes constituting the belt layers are moved to points having as z-coordinates values in which z-coordinates in the basic model at the nodes constituting the layer are multiplied by the value of a scale factor set for each belt layer, and having as r-coordinates values obtained by substituting these values into a function approximating that layer. The value of the scale factor for each belt layer is set as the design variable. Here, in the case of this embodiment, in a case where modeling by means of the finite element method is carried out, because the nodes constituting each belt layer move due to a change in the width of each belt layer, there is the problem that the shapes of elements are too distorted, as shown by a circle 32 in Fig. 18, making it impossible to continue the calculation. To overcome this problem, at the time of determining a corrected model of the tire in Step 115 in Fig. 2, it suffices to add the routine shown in Fig. 19 for automatically performing the division of elements from the given positions of the nodes. In Step 180 in Fig. 19, data is inputted to define the region including each belt layer. This data includes coordinates of nodes and a parameter of node connectivity mesh coarseness for defining each region. As a result of this data input, in an example in which three regions are defined, it is possible to define regions AR1, AR2, and AR3, as shown in Fig. 20A. In an ensuing Step 182, nodes are generated in each region, as shown in Fig. 20B. In an ensuing Step 184, triangular elements are generated by connecting all the nodes including the generated nodes, as shown in Fig. 20C. In an ensuing Step 186, the smoothing of the nodes in the regions is effected by moving the positions of the nodes in the regions, as shown in Fig. 20D. In an ensuing Step 188, the reduction processing of the band width concerning internal processing, such as renumbering, in the finite element method is effected, and this routine ends. It should be noted that although, in this embodiment, polynomials are used as the method of representing the configuration of the belt line, it is possible to use the Lagrange interpolation and the circular interpolation used in the first and second embodiments, or spline curves, B-spline curves, Bezier curves, NURBS, or the like. Tires were actually trial-manufactured in accordance with the above-described sixth embodiment and were tested, and the results of the test are shown in Table 3 below. Conventional Tire Tire of the present invention Time required for design and development 100 42 Main strain between belt layers 100 91 Running distance until a failure in drum test 100 120 Total weight of belts 100 101 Tire size: TBR10.00R20 Internal pressure: 7.25 kg/cm 2 Load: 4050 kgf Speed: 60 km/h Drum diameter: 1.7 m It can be appreciated from Table 3 that the durability improves owing to a reduction of the principal strain between belt layers under a load. Next, a description will be given of a seventh embodiment of the present invention. In this embodiment, the present invention is applied to the determination of the configuration of the tire crown portion. The pressure distribution of a contact area is made uniform and the wear resistance performance is improved without changing the configuration of the area of contact between the tire and the ground. Although the procedure of this embodiment is substantially similar to that of the first embodiment shown in Fig. 2, only the portion of the selection of the objective function, the constraint, and the design variables in Step 102 differs. First, the objective function and the constraint are determined as follows: Objective function: a standard deviation of a pressure distribution in the contact area Constraint: The contact length in the circumferential direction of the tire at the center of the crown portion and the belt end is within ±5% of the contact length in the initial configuration. Next, the configuration of the crown portion which is a design variable is determined by the Lagrange interpolation routine shown in Fig. 21. In Step 170 of this Lagrange interpolation routine, ranges in which the configuration of the crown portion is changed as well as the number of curves for approximating the ranges are specified. Fig. 22 shows an example in which the configuration of the crown portion is represented by two curves, one ranging from a node q3 at the center of the tire to a node q4 at a groove end and the other ranging from a node q5 at the other groove end to a node q6 in the vicinity of a belt end. However, it goes without saying that, for instance, only one curve ranging from the node q3 to the node q6 may be used, or the region between the node q6 and the node q5 and/or between the node q3 and the node q4 may be further divided to represent the configuration of the crown portion by means of three or more curves. In an ensuing Step 172, as shown in Fig. 23, the respective ranges are divided equally by the number of the order of Lagrange interpolation inputted in advance, and nodes which are closest to the equally divided points are selected. In the example shown in Fig. 23, the tire is divided equally by dividing lines Ls1, Ls2, Ls3, Ls4, and Ls5, and nodes qs1, qs2, qs3, qs4, and qs5 which are closest to the dividing lines Ls1 to Ls5 are selected. In an ensuing Step 174, lines normal to the crown portion line at the positions of the selected nodes are calculated, and these normal lines are assumed to be lanes in which the respective nodes move. In this case, as shown in Fig. 24, normal lines concerning the nodes q3, q4, q5, and q6 representing the ranges of the curves and concerning the nodes qs1, qs2, qs3, qs4, and qs5 selected in Step 172 are calculated to determine lines Lp1 to Lp8. Then, in an ensuing Step 176, amounts of movement of the nodes on these lines from their initial positions are selected as design variables, and this routine ends. It should be noted that the amounts of movement of the nodes in the crown portion other than the selected nodes are determined by the Lagrange interpolation in the same way as shown in the first embodiment. Next, a description will be given of an eighth embodiment of the present invention. In this embodiment, the present invention is applied to the determination of the configuration of the tire crown portion in the same way as in the seventh embodiment, but only the design variables differ from those of the seventh embodiment. In addition, although the procedure of this embodiment is substantially similar to that of the second embodiment shown in Fig. 7, only the portion of the selection of the objective function, the constraint, and the design variables in Step 162 differs. First, the objective function and the constraint are determined as follows: Objective function: a standard deviation of a pressure distribution in the contact area Constraint: The contact length in the circumferential direction of the tire at the center of the crown portion and the belt end is within ±5% of the contact length in the initial configuration. Next, as for the configuration of the crown portion which is the design variable, the ranges of the crown portion specified in advance are approximated by a plurality of circular arcs. For instance, as shown in Fig. 25, the configuration of the crown portion is represented by three circular arcs with respect to the ranges of an approximatable crown portion from a node q7 at the tire center to a node q8 in the vicinity of the belt end, i.e., in this case, three regions including ranges CR1, CR2, and CR3. The range CR1 can be approximated at a radius of R1 and an angle of 1; the range CR2 can be approximated at a radius of R2 and an angle of 2; and the range CR3 can be approximated at a radius of R3 and an angle of 3. The above-described example in Fig. 25 shows that the configuration of the crown portion is represented by three circular arcs in the range from the node q7 at the tire center to the node q8 in the vicinity of the belt end. However, it goes without saying that the number of the circular arcs may be 1 or 2, or 4 or more. In addition, the circular arcs need not necessarily be smoothly continuous, and coordinates at the center of each circular arc may be handled as independent variables, as shown in Fig. 26. That is, the range CR1 is approximated at a radius of R1' from a coordinate point O1 (r1, z1) and an angle of 1'; the range CR2 is approximated at a radius of R2' from a coordinate point O2 (r2, z2) and an angle of 2'; and the range CR3 is approximated at a radius of R3' from a coordinate point O3 (r3, z3) and an angle of 3'. Furthermore, a radius Ri and an ascending angle 1 of each circular arc and coordinates ri, zi at the center of each circular arc in the case of the example shown in Fig. 26 are additionally included as design variables. It should be noted that although, in the seventh and eighth embodiments, the Lagrange interpolation and the circular interpolation are used as methods of representing the configuration of the crown portion, it is possible to use spline curves, B-spline curves, Bezier curves, NURBS, or the like. Next, a description will be given of a ninth embodiment of the present invention. In this embodiment, the present invention is applied to the determination of the configuration of the pattern surface of a tire. The pressure distribution at a time when the pattern is in contact with the ground is made uniform, and the wear resistance performance is improved. Although the procedure of this embodiment is substantially similar to that of the first embodiment shown in Fig. 2, only the portion of the selection of the objective function, the constraint, and the design variables in Step 102 differs. First, the objective function and the constraint are determined as follows: Objective function: a standard deviation of a pressure distribution in the contact area Constraint: The total volume of the pattern is within ±5% of the initial volume. With respect to the configuration of the pattern surface which is the design variable, the pattern surface is divided into the form of a lattice in accordance with the order of Lagrange interpolation inputted in advance, as shown in Fig. 27. The coordinate r i in the direction of the pattern thickness at each point (16 points in the example of Fig. 27) thus obtained is set as the design variable. In a case where the finite element method is used in the calculation of the objective function, there is a need to further subdivide the pattern surface. For example, in the case of Fig. 28 in which the pattern surface shown in Fig. 27 is subdivided, the coordinates of nodes other than the nodes corresponding to the design variable on the pattern surface may be calculated in accordance with the following formulae: x = i = 1 n M i · x i y = i = 1 n M i · yi r = i = 1 n M i · r i Here, M i ≡ M em = N e (ξ) · N m (η) where, n = k 1 · k 2 k 1 : order of Lagrange interpolation in the x-direction + 1 k 2 : order of Lagrange interpolation in the y-direction + 1 As a result, the determination and correction of the configuration of the pattern surface corresponding to Steps 108 and 115 in Fig. 2 become possible. It should be noted that although the Lagrange interpolation is used as a method of representing the configuration of the pattern surface. it is possible to alternatively use such as the Bezier curved surface, B-spline curved surface, or the Coons curved surface which are shown in Mathematical Elements for Computer Graphics mentioned above. Parameters for controlling them may be directly selected as design variables. The foregoing embodiments need not be carried out independently. For instance, various combinations are possible such as by carrying out the determination of the ply lines in the first embodiment and the determination of the belt structure in the sixth embodiment, or by carrying out the determination of the configuration of the crown portion in the seventh embodiment by using the ply lines determined in the first embodiment. Tires obtained by combining the seventh and ninth embodiments were actually trial-manufactured, and were mounted on a vehicle as right-hand front wheels. The vehicle was run 10,000 km on a general road, and the results of investigating the ratio of amounts of wear (shoulder wear) between the tread center and a tread side end are shown in Table 4 below. Conventional Tire Tire of the present invention Time required for design and development 100 36 Ratio of amounts of wear (wear at side end / wear at center portion) 2.65 1.50 Standard deviation of pressure distribution 100 41 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 It can be appreciated from Table 4 that the shoulder wear resistance improves owing to the fact that the pressure distribution is made uniform. Next, a description will be given of a tenth embodiment in which the configuration of a carcass line is designed genetically by means of an algorithm so as to render the belt tension in the circumferential direction of the tire when inflated with air a maximum value, i.e., an optimum value, in order to improve steering stability. It should be noted that since the arrangement of this embodiment is substantially similar to that of the above-described embodiments, identical portions will be denoted by the same reference numerals, and a detailed description thereof will be omitted. Fig. 29 shows a program processing routine in accordance with the tenth embodiment. In Step 200, cross-sectional configurations of N tires are subjected to modeling by a technique which is capable of numerically and analytically determining the belt tension in the circumferential direction of the tire when inflated with air such as by the finite element method, and a basic model of the tire including its internal structure is determined. The term modeling referred to herein means the numerical representation of the shape, structure, materials, and pattern of the tire into a format of input data for a computer program prepared according to a numerical and analytical technique. It should be noted that N is inputted in advance by the user. The basic model of the tire used in this embodiment is identical to the one shown in Fig. 4 illustrating the first embodiment. It should be noted that, in the same way as in the first embodiment, the basic model of the tire may be divided into a plurality of elements by a plurality of lines normal to a line representing an outer configuration of the tire or by a plurality of lines normal to a turn-up ply line, or may be divided into arbitrary shapes such as triangles depending on a design objective. In an ensuing Step 202, an objective function representing a physical amount for evaluating tire performance, a constraint for constraining the cross-sectional shape of the tire, and design variables determining the cross-sectional shapes of N tire models are determined. In this embodiment, the objective function OBJ and the constraint G are defined as follows in order to design the configuration of the carcass line for maximizing the belt tension in the circumferential direction of the tire when inflated with air: Objective function OBJ: the sum total, in the widthwise direction of the tire, of components of the belt tension acting in the circumferential direction of the tire Constraint G: A periphery value of the carcass line is within ±5% of a set value. It should be noted that the aforementioned periphery value of the carcass line can be calculated as the sum total of distances between nodes (points of intersection between the carcass line and the respective normal line) in the carcass line which are located in a domain that changes the tire shape. In addition, the configuration of the carcass line, which is a design variable, is determined by the Lagrange interpolation routine, shown in Fig. 3, for approximating a curve with respect to the N tire models, respectively. Since this Lagrange interpolation routine is identical to that of the first embodiment, a description thereof will be omitted. By repeating the Lagrange interpolation routine N times, the objective function OBJ, the constraint G, and respective design variables r iJ (J = 1, 2, ..., N) of the N tire models are determined. Subsequently, in Step 204 in Fig. 29, objective functions OBJ J and constraints G J are calculated for the respective design variables r iJ of the N tire models. In an ensuing Step 206, by using the objective functions OBJ J and constraints G J of the N tire models determined in Step 204, adaptive functions F J of the N tire models are calculated in accordance with the following Formulae (4). In this embodiment, the value of the fitness function (degree of adaptation) becomes larger as the belt tension becomes larger so as to maximize the belt tension. Φ j = -OBJ J + γ · max(G J , O) F J = -Φ J or F J = 1/Φ J or F J = -a · Φ J + b where, a = Φ avg (c - 1) (Φ avg - Φ min ) b = Φ avg (c - Φ min ) (Φ avg - Φ min ) Φ avg = j = 1 N Φ j N c: constant γ: penalty coefficient Φ min : min(Φ 1 , Φ 2, ..., Φ N ) Φ j : penalty function of a J-th tire model among the N tire models (J = 1, 2, 3, ..., N) It should be noted that c and γ are inputted in advance by the user. In an ensuing Step 208, two models to be crossed over each other are selected from among the N models. As the method of selection, an adaptation-degree proportional strategy, which is generally known, is used. A probability P 1 with which certain individuals ℓ among the N tire models are each selected in the selection is expressed by the following formula: P 1 = F 1 j = 1 N F j where, F 1 : fitness function of a certain individual ℓ among the N tire models F J : J-th fitness function among the N tire models (J = 1, 2, 3, ..., N) Although, in the above-described embodiment, the adaptation-degree proportional strategy is used as the method of selection, it is possible to alternatively use an expected-value strategy, a rank strategy, an elite preservation strategy, a tournament selection strategy, a GENITOR algorithm, or the like, as shown in "Genetic Algorithms" (edited by Hiroaki Kitano). In an ensuing Step 210, a determination is made as to whether or not the two selected tire models are to be crossed over each other by a probability T inputted in advance by the user. The term crossover used herein means the exchange of certain elements of two tire models, as will be described later. If NO is the answer in the determination and crossover is not be carried out, in Step 212, the present two tire models are kept intact, and the operation proceeds to Step 216. Meanwhile, if YES in the answer in the determination and crossover is to be carried out, in Step 214, the two tire models are crossed over each other, as will be described later. The crossover of the two tire models is carried out in accordance with an intersecting routine shown in Fig. 30. First, the two tire models selected in Step 208 are set as a tire model a and a tire model b , and the design variables of the tire models a and b are expressed by design variable vectors including lists. The design variable vectors of the tire model a are set as Vr a = (r 1 a , r 2 a , ..., r i a , ..., r n-1 a ) , while the design variable vectors of the tire model b are set as Vr b = (r 1 b , r 2 b , ..., r i b , ..., r n-1 b ) . In Step 250 in Fig. 30, predetermined random numbers are generated, and a point of crossover i concerning the design variable vectors of the tire models a , b is determined in accordance with the random numbers. In an ensuing Step 252, a distance d is determined in accordance with the following formula with respect to the design variables r 1 a , r 1 b of the tire models a, b which are determined to cross over each other: d = |r i a - r i b | In an ensuing Step 254, a normalized distance d' is determined in accordance with the following formula by using a minimum value B L and a maximum value B U which fall in ranges that can be assumed by r 1 a , r 1 b . d' = d B U - B L In Step 256, to disperse the values of the normalized distance d' appropriately, a function Z ab is determined in accordance with the following formula by using a chevron-shaped mapping function Z(x) (0 ≤ x ≤ 1, 0 ≤ Z(x) ≤ 0.5) such as the one shown in Fig. 31A or Fig. 31B: Z ab = Z(d') After the function Zab is thus determined, new design variables r i ' a , r i ' b are determined in accordance with the following formulae in Step 258. r i ' a = r i a - min(|r i a - B L |, |r i a - B U |) 0.5 · Z ab r i ' b = r i b + min(|r i b - B L |, |r i a - B U |) 0.5 · Z ab or r i ' a = r i a + min(|r i a - B L |, |r i a - B U |) 0.5 · Z ab r i ' b = r i b - min(|r i b - B L |, |r i a - B U |) 0.5 · Z ab After r i ' a and r i ' b are thus determined, in Step 260, design variable vectors Vr' a and Vr' b which are lists of new design variables are determined as follows: Vr' a = (r 1 a , r 2 a , ..., r i ' a , r i+1 b , ..., r n-1 b ) Vr' b = (r 1 b , r 2 b , ..., r i ' b , r i+1 a , ..., r n-1 a ) It should be noted that the minimum value B L and the maximum value B U which fall in ranges that can be assumed by r 1 a , r 1 b are inputted in advance by the user. In addition, the mapping functions Z(x) may be trough-shaped functions such as those shown in Figs. 32A and 32B. Although there is only one point of crossover i in the above-described example, it is possible to alternatively use a multipoint crossover or uniform crossover such as those shown in "Genetic Algorithms" (edited by Hiroaki Kitano). After two new tire models are generated by such a crossover, in Step 216 in Fig. 29, a determination is made as to whether or not a mutation is to take place at a probability S inputted in advance by the user. The term mutation referred to herein means changing a portion of the design variable by an infinitesimal degree, and is aimed at enhancing the probability of including a population capable of assuming optimum design variables. If NO is the answer in the determination in Step 216 and a mutation is not to take place, in Step 218, the present two tire models are kept intact, and the operation proceeds to an ensuing Step 222. If YES is the answer in the determination and a mutation is to take place, mutation processing is carried out as follows in an ensuing Step 220. This mutation is carried out by a mutation routine shown in Fig. 33. First, in Step 262, random numbers are generated, and the place of mutation i is determined by the random numbers. In an ensuing Step 264, the distance d' is determined by the random numbers in the following range: 0 ≤ d' ≤ 1 In an ensuing Step 266, the function Zd (0 ≤ x ≤ 1, and 0 ≤ Z(x) ≤ 0.5) is determined in accordance with the following formula by using the chevron-shaped mapping function such as the one shown in Fig. 31A or 31B or the trough-shaped function such as the one shown in Fig. 32A or 32B: Zd = Z(d') After the function Zd is thus determined, in Step 268, a new design variable r i ' is determined in accordance with the following formula: r i ' =r i - min(|r i - B L |, |r i - B U |) 0.5 · Zd r i ' =r i + min(|r i - B L |, |r i - B U |) 0.5 · Zd After the design variable r i ' is thus determined, a design variable vector Vr', i.e., a list of new design variables which is determined in Step 270, becomes as follows: Vr' = (r 1 , r 2 , ..., r i ', r i+1 , ..., r n-1 ) With respect to the two tire models newly generated in the above-described manner, the value of the objective function and the value of the constraint are calculated in Step 222 in Fig. 29. In an ensuing Step 224, a fitness function is calculated from the resultant values of the objective function and the constraint by using Formula (4) in the same way as in the preceding embodiments. In an ensuing Step 226, the aforementioned two tire models are preserved. In an ensuing Step 228, a determination is made as to whether or not the number of tire models preserved in Step 226 has reached N, and if it has not reached N, Steps 208 through 228 are executed repeatedly until it reaches N. Meanwhile, if the number of tire models has reached N, a determination is made with respect to convergence in Step 230. If a convergence has not been attained, the N tire models are updated to the tire models preserved in Step 226, and Steps 208 through 230 are executed repeatedly. Meanwhile, if it is determined in Step 230 that a convergence has been attained, the values of the design variables of a tire model which give a maximum value of the objective function while substantially satisfying the constraint among the N tire models are set as values of the design variables which maximize the objective function while substantially satisfying the constraint. In Step 232, the configuration of the tire is determined by using the values of these design variables. It should be noted that, as for the determination with respect to convergence in Step 230, it is assumed that a convergence has been attained if any one of the following conditions is met: 1) The number of generations has reached M. 2) The number of line rows in which the value of the objective function is the largest has accounted for q% of the total. 3) The maximum value of the objective function is not updated in subsequent p generations. It should be noted that M, q, and P are inputted in advance by the user. In addition, the above-described embodiment may be applied to the design variables in the second, third, fourth, fifth, sixth, seventh, eighth, and ninth embodiments. Tires were actually trial-manufactured in accordance with the above-described tenth embodiment and were tested, and the results of the test are shown in Table 5 below. Conventional Tire 10th Embodiment Time required for design and development 100 53 Performance Steering stability 100 108 Riding comfort 100 99 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 Testing method: feeling evaluation by an actual-vehicle running test Thus, in the tenth embodiment, since the amount of calculation increases as compared to the first embodiment shown in Table 1, the time required for design and development increases slightly. However, there is an advantage in that tires of a better performance can be designed. Next, a description will be given of an eleventh embodiment of the present invention. In this embodiment, the present invention is applied to the belt structure, and is aimed at optimizing mutually incompatible performances. In this embodiment, a belt structure for improving the steering stability without impairing the riding comfort of the occupant in the vehicle is determined. Although the procedure of this embodiment is substantially similar to that of the program processing routine of the tenth embodiment, the selection of the objective function, the constraint, and the design variables, as well as the method of crossover and the method of mutation differ. In this embodiment, lateral stiffness, i.e., a physical amount for improving the steering stability, is used as the objective function. A condition that a vertical spring constant, i.e., a physical amount for governing the riding comfort, is fixed is used as the constraint. A belt structure is determined such that the lateral stiffness becomes maximum under the condition of the vertical spring constant being fixed. In addition, when a basic model of the tire such as the one shown in Fig. 34A is fabricated, places b 1 , b 2 , ..., b i , ..., b m where the belt is capable of being present are set as design variables. Further, it is assumed that if b i = 1, the belt is present at that place, and if b i = 0, the belt is not present at that place. For instance, if the belt is present as indicated by bold lines in Fig. 34B, the design variables become 1011001100 in the order of b 1 , b 2 , ..., b n , and become 1111111001 in the order of b n+1 , b n+2 , ..., b m . Next, the crossover is carried out on the basis of the processing routine shown in Fig. 35. First, it is assumed that a-th and b-th tire models among the N tire models are to cross over each other. The lists of the design variables of the respective tire models are set as design variable vectors: VB a = (b 1 a , b 2 a , ..., b i a , b i+1 a , ..., b m a ) VB b = (b 1 b , b 2 b , ..., b i b , b i+1 b , ..., b m b ) In Step 272, random numbers are generated, and the point of crossover i is determined in accordance with the random numbers generated. In an ensuing Step 274, the lists of design variables which are design variable vectors are altered as follows, to determine design variable vectors VB a ' and VB b ' as lists of new design variables. VB a ' = (b i a , b 2 a , ..., b i-1 a , b i b , b i+1 b , ..., b m b ) VB b ' = (b i b , b 2 b , ..., b i-1 b , b i a , b i+1 a , ..., b m a ) In an ensuing Step 276, two new tire models are generated in accordance with the design variable vectors VB a ' and VB b ' obtained. Although there is only one point of crossover i in this example, it is possible to alternatively use a multipoint crossover or uniform crossover such as those shown in "Genetic Algorithms" (edited by Hiroaki Kitano). Next, the mutation is carried out on the basis of the processing routine shown in Fig. 36. It is assumed that the list of design variables of the basic model of the tire are the design variable vectors: VB = (b 1 , b 2 , ..., b i-1 , b 1 , b i+1 , ..., b m ) In Step 278, random numbers are generated, and the place of mutation i is determined in accordance with the random numbers. In an ensuing Step 280, a design variable b i at the place of mutation i is changed as shown below, so as to generate a design variable vector VB' as a list of new design variables. If b i = 0, then b i ' = 1, and If b i = 1, then b i ' = 0. VB' = (b 1 , b 2 , ..., b i-1 , b i ', b i+1 , ..., b m ) In an ensuing Step 282, a new tire model is generated from the design variable vector VB' thus determined. It should be noted that although, in the above-described embodiment, the belt structure is adopted as the design variable, other reinforcing material may be determined as the design variable. In addition, this embodiment is applicable to a case where a decision is made as to whether or not a reinforcing material is to be inserted in the bead portion. Tires were actually trial-manufactured in accordance with the above-described eleventh embodiment and were tested, and the results of the test are shown in Table 6 below. Conventional Tire 11th Embodiment Time required for design and development 100 65 Performance Steering stability 100 105 Riding comfort 100 100 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 Testing method: feeling evaluation by an actual-vehicle running test Thus, in the eleventh embodiment, since a design variable which gives an optimum value of the objective function satisfying the constraint is determined, and the tire is designed from this design variable, there is an advantage in that tires in a best mode can be designed. Next, a description will be given of a twelfth embodiment of the present invention. In this embodiment, the present invention is applied to the determination of the configuration of blocks of a tire having blocks. This embodiment is aimed at reducing the difference in rigidity among blocks having different sizes and arranged in the circumferential direction of the tire, and at improving the irregular wear by a pattern which makes uniform the block rigidity in various directions in one block. Fig. 37 shows a program processing routine in accordance with this embodiment. In Step 300, as shown in Fig. 38A, the configuration of a block 42 is defined by the coordinates of nodes of the block 42, the coordinates of nodes of a sipe 44, the angle of a block groove, the depth of the groove, the number of sipes, the sipe width, the angle at which the sipe is embedded, and the sipe depth, and is subjected to modeling so as to allow analysis by the finite element method and the analytical method. Fig. 38B shows a modeling example of a finite element model. In an ensuing Step 302, the following are determined: an objective function representing a physical amount for evaluating the performance of the block; a constraint for constraining the sipe length, the sipe length × the sipe depth, the number of sipes and the like; and design variables concerning the coordinates of the nodes of the block, the coordinates of the nodes of the sipe, the groove angle at each side of the block, the groove depth at each side of the block, the sipe width, the angle at which the sipe is embedded, and the sipe depth. In this embodiment, the objective function, the constraint and the design variables are determined as follows so as to design a pattern in which the block rigidity in various directions in one block is made uniform: Objective function OBJ: The block rigidity in various directions is made uniform. Constraint G: The amount of change in the sipe length is within ±5% of the sipe length in an initial model. The coordinates of the node of the sipe are within a figure enclosed by the coordinates of the nodes of the block. The distance between the coordinates of a node of the sipe and a straight line representing an outer portion of the block is 2 mm or more. Design variable: Coordinates of the nodes of the sipe (coordinates of the nodes of the block are fixed) In an ensuing Step 304, an initial value OBJ 0 of the objective function OBJ at an initial value X 0 of a design variable X i and an initial value G 0 of the constraint G are calculated. In an ensuing Step 306, the design variable X i is changed by ΔX i each so as to change the basic model of the block. In Step 308, a block which includes changes in the design variable carried out in Step 306 is subjected to modeling. In Step 310, with respect to the corrected model of the block determined in Step 308, the value OBJ i of the objective function and the value Gi of the constraint after the design variable is changed by ΔX i are calculated. In Step 312, a sensitivity dOBJ/dX i of the objective function, which is a ratio of an amount of change of the objective function to an amount of unit change of the design variable, as well as a sensitivity dG/dX i of the constraint, which is a ratio of an amount of change of the constraint to an amount of unit change of the design variable, are calculated for each design variable in accordance with the following formulae: dOBJ dr i = OBJ i - OBJ 0 ΔX i dG dr i = G i - G 0 dX i By means of these sensitivities, it is possible to estimate to what extent the value of the objective function and the value of the constraint change when the design variable is changed by ΔX i . It should be noted that these sensitivities are sometimes determined analytically depending on the technique used in modeling the block and the nature of the design variable, in which case the calculation of Step 312 becomes unnecessary. In an ensuing Step 314, an amount of change of the design variable which maximizes the objective function while satisfying the constraint is estimated by means of mathematical programming by using the initial value OBJ 0 of the objective function, the initial value G 0 of the constraint, the initial value X 0 of the design variable, and the sensitivities. In Step 315, by using this estimated value of the design variable, a corrected model of the block is determined in a method similar to that of Step 308, and a value of the objective function is calculated. In Step 316, a determination is made as to whether or not the value of the objective function has converged by comparing a threshold value inputted in advance with the difference between the value OBJ of the objective function calculated in Step 315 and the initial value OBJ 0 of the objective function calculated in Step 304. If the value of the objective function has not converged, Steps 304 through 316 are executed repeatedly by setting as the initial value the value of the design variable determined in Step 314. If it is determined that the value of the objective function has converged, the value of the design variable at this juncture is set as the value of the design variable which maximizes the objective function while satisfying the constraint, and the shape of the block is determined in Step 320 by using this value of the design variable. Tires were actually trial-manufactured in accordance with the above-described twelfth embodiment and were tested, and the results of the test are shown in Table 7 below. Conventional Tire 12th Embodiment Time required for design and development 100 20 Uniformity in block rigidity 100 40 H & T height 100 80 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 Testing method: 5000 km running on a test course It should be noted that the uniformity (BK) in block rigidity is obtained by actually measuring block rigidity Ki (i = 1 to 6) for each 30° of the block and by using the following formula: BK = i=1 6 (Ki - K ) 2 K K = 1 6 i=1 6 Ki The H & T height (heel-and-toe height) is a phenomenon in which the blocks connected in the circumferential direction of the tire are worn in the form of saw teeth in the circumferential direction. This H & T height is shown in Fig. 41. Next, a description will be given of a thirteenth embodiment. This embodiment is a combination of the first and tenth embodiments. It should be noted that since the arrangement of this embodiment is substantially similar to those of the aforementioned embodiments, identical portions will be denoted by the same reference numerals, and a detailed description will be omitted. In the tenth embodiment, the objective function and the constraint are calculated on the basis of a draft design obtained through a crossover and/or a mutation. In this case, this may not lapse into a local optimal solution as described in Goldberg, D. E., "Genetic Algorithms in Search, Optimization and Machine Learning," Addison-Wesley (1989), and yet there is a problem in that it is difficult to obtain a true optimal solution. Accordingly, as the arithmetic processing in Step 222 in the tenth embodiment, if the methods of the first embodiment and the tenth embodiment are combined by using the processing in Steps 104 to 116 in the first embodiment, the aforementioned problem can be overcome. Fig. 39 show a program processing routine in accordance with this embodiment. Since Steps 200 to 220 are identical to those of the tenth embodiment, the description thereof will be omitted. In Step 340, the objective function and the constraint are calculated by using as an initial draft design the two draft designs obtained as described above. The subsequent Steps 106 to 116 are executed repeatedly until the value of the objective function converges in the same way as in the processing of the first embodiment. If it is determined that the value of the objective function has converged, in an ensuing Step 224, a fitness function is calculated from the resultant values of the objective function and the constraint. Then in an ensuing Step 226, the aforementioned two tire models are preserved. Steps 208 through 228 are executed repeatedly until the number of tire models preserved in Step 226 reaches N. If it has reached N, a determination is made with respect to convergence in Step 230 in the same way as described above. If a convergence has been attained, the values of the design variables of a tire model which give a maximum value of the objective function while substantially satisfying the constraint among the N tire models are set as values of the design variables which maximize the objective function while substantially satisfying the constraint. In Step 232, the configuration of the tire is determined by using the values of these design variables. By the method for obtaining a true optimal solution without lapsing into a local optimal solution, apart from the technique described herein, it is possible to combine the method of the tenth embodiment with a method called simulated annealing which is described in the aforementioned reference. Next, a description will be given of a fourteenth embodiment. In this embodiment, the process leading to an optimal tire configuration is processed speedily. It should be noted that since the arrangement of this embodiment is substantially similar to those of the aforementioned embodiments, identical portions will be denoted by the same reference numerals, and the detailed description thereof will be omitted. Fig. 40 shows a program processing routine in accordance with this embodiment. In Step 400, configuration data in a predetermined number M specified in advance are read from the data base concerning tire configurations. In the tire-configuration data base, the carcass line of the tire is represented by the Lagrange interpolation in Formula (2), the method using circular arcs shown in Fig. 10, spline curves, Bezier curves, NURBS (non-uniform rational B-splines), or the like. In an ensuing Step 402, a carcass line F(x) of the tire is represented by the following Formula (5): F(x) = i = 1 M a i f i (x) where a i (i = 1 to M) is a constant, and f i (x) is read from the tire-configuration data base. Next, as shown in Fig. 4, a plurality of normal lines NL 1 , NL 2 , NL 3 , ... are drawn with respect to the carcass line, and nodes are distributed on these lines by taking into consideration the thickness corresponding to the tire member. Then, the cross-sectional configuration of the tire is divided into a plurality of elements by making use of the nodes. In an ensuing Step 404, an objective function representing a physical amount for evaluating tire performance, a constraint for constraining the cross-sectional shape of the tire, and a design variable determining the cross-sectional shape of the tire are determined. In addition, the configuration of the carcass line which is the design variable is determined by the constant a i of Formula (5). In an ensuing Step 406, the design variable a i is changed by Δa i each to change the basic model of the tire. In an ensuing Step 408, the carcass line of the tire after the design variable a i is changed by Δa i is obtained from Formula (5), and a corrected model of the tire is determined in accordance with the procedure of Step 402. Step 110 and and subsequent steps are similar to those described above, and therefore the description thereof will be omitted. If the tires are designed in accordance with the above-described fourteenth embodiment, it is possible to optimize the configuration by combining the configuration data base of actual tires and the configuration inputted newly by the designer. Hence, it is possible to obtain an optimal solution with less calculating time as compared with the first embodiment. However, since optimal solutions are dependent on the configuration data base which is read initially, there is a tendency that the advantage is slightly lessened in comparison with the first embodiment. Table 8 below shows the results of trial manufacture and test of tires designed with the procedure of the fourteenth embodiment on the basis of the same objective function and constraint as those of the first embodiment. Conventional Tire 1st Embodiment 14th Embodiment Time required for design and development 100 50 30 Performance Steering stability 100 106 105 Riding comfort 100 98 98 Tire size: 205/60HR15 Internal pressure: 2.0 kg/cm 2 Testing method: feeling evaluation by an actual-vehicle running test

Claims: A computer implemented method of designing a pneumatic tire, comprising the steps of: (a) defining a basic model of a tire representing a cross-sectional shape of the tire including an internal structure, an objective function representing a physical amount for evaluating tire performance, a design variable for determining at least one of the cross-sectional shape of the tire and a tire structure, and a constraint other than the design variable for constraining at least one of the cross-sectional shape of the tire, the tire structure, a physical amount for evaluating performance, and a tire size; (b) varying a value of the design variable while taking the constraint into consideration and carrying out an optimization computation until an optimum value of the objective function is obtained; and (c) designing the tire on the basis of the design variable which gives the optimum value of the objective function. A method as claimed in claim 1, characterized in that said step (b) comprises: estimating an amount of change of the design variable which gives the optimum value of the objective function while taking the constraint into consideration on the basis of a sensitivity of the objective function, which is a ratio of an amount of change of the objective function to an amount of unit change of the design variable, and a sensitivity of the constraint, which is a ratio of an amount of change of the constraint to an amount of unit change of the design variable; calculating a value of the objective function when the design variable is changed by an amount corresponding to the estimated amount and a value of the constraint when the design variable is changed by the amount corresponding to the estimated amount; and determining the value of the design variable which gives the optimum value of the objective function on the basis of the estimated value and the calculated values while taking the constraint into consideration. A method as claimed in claim 1, characterized in that said step (a) comprises: determining a group subject to selection, consisting of a plurality of basic models of the tire each representing a cross-sectional shape of the tire including an internal structure thereof; and determining, with respect to the basic models of the tire of the group subject to selection, the objective function representing the physical amount for evaluating tire performance, the design variable for determining the cross-sectional shape of the tire and/or the tire structure, the constraint for constraining at least one of the cross-sectional shape of the tire, the tire structure, the physical amount for evaluating performance, and the tire size, and a fitness function which can be evaluated from the objective function and the constraint; and in that said step (b) comprises: selecting two basic models of the tire from the group subject to selection on the basis of the fitness function; effecting at least one of generating a new basic model of the tire by allowing design variables of the basic models of the tire to cross over each other at a predetermined probability and of generating a new basic model of the tire by changing a portion of the design variable of at least one of the basic models of the tire; determining the objective function, the constraint, and the fitness function of the basic model of the tire with the design variable changed; preserving said basic model of the tire and the basic model of the tire with the design variable not changed, an aforementioned process of step (b) being repeated until the basic models of the tire preserved reaches a predetermined number; determining whether or not a new group consisting of the predetermined number of the preserved basic models of the tire satisfies a predetermined convergence criterion; if the convergence criterion is not satisfied, repeating a foregoing process of step (b) until the group subject to selection satisfies a predetermined convergence criterion by setting the new group as the group subject to selection; and if the predetermined convergence criterion is satisfied, determining the value of the design variable which gives the optimum value of the objective function among the predetermined number of the preserved basic models of the tire while taking the constraint into consideration. A method as claimed in claim 3, characterized in that, in said step (b), with respect to the basic model of the tire with the design variable changed, the amount of change of the design variable which gives the optimum value of the objective function is estimated while taking the constraint into consideration on the basis of the sensitivity of the objective function, which is the ratio of the amount of change of the objective function to the amount of unit change of the design variable, and on the basis of the sensitivity of the constraint, which is the ratio of the amount of change of the constraint to the amount of unit change of the design variable, the value of the objective function when the design variable is changed by the amount corresponding to the estimated amount and the value of the constraint when the design variable is changed by the amount corresponding to the estimated amount are calculated, the fitness function is determined from the value of the objective function and the value of the constraint, said basic model of the tire and the basic model of the tire with the design variable not changed are preserved, and an aforementioned process of step (b) is repeated until the preserved basic models of the tire reaches the predetermined number. A method as claimed in claim 1, characterized in that the design variable includes at least one of: a function representing the shape of at least one line selected from a carcass line, a turn-up line, a line representing an outer configuration of the tire, a line representing a tire crown portion, and a reinforcing-material line; a variable representing at least one gauge distribution of a tire rubber member selected from a gauge distribution of a bead filler, a gauge distribution of a rubber chafer, a gauge distribution of side rubber, a gauge distribution of tread rubber, a gauge distribution of tread base rubber, a gauge distribution of inner-surface reinforcing rubber, a gauge distribution of inter-belt rubber, and a gauge distribution of belt end rubber; a variable representing at least one structure of a belt portion selected from an angle of each belt layer, a width thereof, a cord type thereof; and a placing density thereof; and a variable representing at least one configuration of a pattern selected from a configuration of a block, a position of a sipe, a number of sipes, and a length of the sipe. A method as claimed in claim 1, characterized in that a plurality of lines are set as the design variable by representing the shape of at least one line selected from a carcass line, a turn-up ply line, a line representing an outer configuration of the tire, a line representing a tire crown portion, and a reinforcing-material line, and in that, in order to prevent the plurality of lines from intersecting with each other when the design variable is varied in step (b), a mapping function of a function representing shapes of the lines and a constraint for preventing the lines from intersecting with each other are determined, and the value of the design variable which gives the optimum value of the objective function while satisfying the constraint is determined by executing step (b) in a space mapped by the mapping function. A method as claimed in claim 1, characterized in that the objective function is a sum total, in a widthwise direction of the tire, of components of belt tension acting in a circumferential direction of the tire. A method as claimed in claim 1, characterized in that the constraint is that a periphery value of a carcass line is within a predetermined range which is set with the periphery value of the carcass line of a basic model of the tire as a reference. A method as claimed in claim 1, characterized in that the objective function is a lateral stiffness which is a physical amount for improving steering stability. A method as claimed in claim 1, characterized in that the constraint is a condition that a vertical 1st order natural frequency which is a physical amount for governing riding comfort is fixed. A method as claimed in claim 1, characterized in that, if a plurality of lines are set as design variables, in order to prevent the plurality of lines from intersecting with each other when the design variables are varied in step (b), mapping functions of functions representing the shapes of the lines and a constraint for preventing the lines from intersecting with each other are determined, and values of the design variables which give the optimum value of the objective function while satisfying the constraint are determined by executing step (b) in a space mapped by each of the mapping functions. A computer implemented method of designing a pneumatic tire, comprising the steps of: (a) setting a basic model of the tire representing a cross-sectional configuration of the tire including an internal structure; (b) selecting an objective function among at least one physical amount for evaluating tire performance; (c) defining a range of a design variable from at least one of a portion of the cross-sectional configuration of the tire and a tire structure, the range of the design variable providing for changes in an optimization computation; (d) selecting a constraint other than the design variable for constraining at least one of the cross-sectional configuration of the tire, the tire structure, a physical amount for evaluating performance, and a tire size; (e) varying the value of the design variable while taking the constraint into consideration and carrying out the optimization computation until an optimum value of the objective function is obtained; and (f) designing the tire based on the design variable that provides the optimum value of the objective function.

Inventors: ABE AKIHIKO, KAMEGAWA TATSUHIKO, NAKAJIMA YUKIO, ABE, AKIHIKO, KAMEGAWA, TATSUHIKO, NAKAJIMA,YUKIO
Applicants: BRIDGESTONE CORP, BRIDGESTONE CORPORATION
Application year: 1994
Main classification: B29D30
IPC codes: B60C3, G06F17