+

WO1994016877A1 - Design method for a pneumatic tire - Google Patents

Design method for a pneumatic tire Download PDF

Info

Publication number
WO1994016877A1
WO1994016877A1 PCT/JP1994/000110 JP9400110W WO9416877A1 WO 1994016877 A1 WO1994016877 A1 WO 1994016877A1 JP 9400110 W JP9400110 W JP 9400110W WO 9416877 A1 WO9416877 A1 WO 9416877A1
Authority
WO
WIPO (PCT)
Prior art keywords
tire
design
value
objective function
line
Prior art date
Application number
PCT/JP1994/000110
Other languages
English (en)
French (fr)
Inventor
Tatsuhiko Kamegawa
Yukio Nakajima
Akihiko Abe
Original Assignee
Bridgestone Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Bridgestone Corporation filed Critical Bridgestone Corporation
Priority to DE69419678T priority Critical patent/DE69419678T2/de
Priority to JP51687394A priority patent/JP4142094B2/ja
Priority to EP94905215A priority patent/EP0642916B1/en
Priority to US08/302,841 priority patent/US5710718A/en
Publication of WO1994016877A1 publication Critical patent/WO1994016877A1/ja

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60CVEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C3/00Tyres characterised by the transverse section
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60CVEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C3/00Tyres characterised by the transverse section
    • B60C3/04Tyres characterised by the transverse section characterised by the relative dimensions of the section, e.g. low profile
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60CVEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C9/00Reinforcements or ply arrangement of pneumatic tyres
    • B60C9/02Carcasses
    • B60C9/0292Carcass ply curvature
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60CVEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C99/00Subject matter not provided for in other groups of this subclass
    • B60C99/006Computer aided tyre design or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10TTECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
    • Y10T152/00Resilient tires and wheels
    • Y10T152/10Tires, resilient
    • Y10T152/10495Pneumatic tire or inner tube

Definitions

  • the present invention relates to a method for designing a pneumatic tire, and more particularly, to efficiently and easily develop and design a structure and shape of a tire that achieves a single purpose performance, a trade-off performance, etc.
  • the present invention relates to a pneumatic tire design method that requires a simple structure and shape and can design a high-performance tire. Background art
  • the conventional tire design method is based on empirical rules obtained by repeating experiments and numerical experiments using experiments and computers. As a result, the number of prototypes and tests required for development became enormous, increasing development costs and making it difficult to shorten the development period.
  • One of the design methods for the tire force-cass shape is a design method called the natural equilibrium shape, which acts on the tire carcass cord, especially when the tire is assembled on the rim and filled with air to make the tire internal pressure normal. It also determines the cross-sectional shape of the carcass (cross-sectional shape when cut along a plane that includes the rotation axis of the shaft) so that the tension is constant. That is, by making the cross-sectional shape of the force “scrap” when the tire is vulcanized and cured in a mold substantially equal to the cross-sectional shape of the carcass at normal internal pressure, the various reinforcing layers of the tire during air filling are formed.
  • non-natural equilibrium shape This is to prevent the occurrence of failure due to internal strain or the like that occurs partially between or within the laminates.
  • the other one is a design method called non-natural equilibrium shape.
  • a typical non-natural flat mouth shape design method is already a trademark called “RC ⁇ T” proposed by the applicant. ⁇ ⁇ ⁇
  • design methods that are sold and are commercially successful.
  • this non-natural equilibrium shape design method is a design method that controls the tension of the carcass. By controlling the tension, it is intended to obtain the required performance.
  • the crown shape of the tire is designed by several arcs in the cross section including the rotation axis of the tire.
  • the method of determining the value of the arc is 1) the method of determining it from the database of the value of the arc and tire performance, 2) the method of creating several molds, and the prototype of tires using this mold. It was a method of conducting many experiments and making decisions. For this reason, development efficiency was not good.
  • Belt structure is determined by belt width, belt angle, cord type, driving density, and rubber gauge distribution between belts.
  • the method of designing the belt structure is as follows: (1) a method of measuring the strain at the end of the belt by experiments and numerical experiments; (2) a method of measuring the rigidity of the entire belt; and (3) an indoor evaluation and an actual vehicle evaluation by changing the belt structure.
  • the method was to extract design proposals. Therefore, there were problems such as a long evaluation period.
  • pattern design Since the pattern design has a high degree of freedom, the basic pattern draft is grouped on the tires or after actually making a mold, a prototype is prototyped, and the actual vehicle is evaluated.If a problem occurs in the actual vehicle, the basic pattern draft is finely corrected. And solved it to make the final pattern. For this reason, pattern design was an area that required the most man-hours compared to tire shape and structural design.
  • the present invention is capable of designing the best mode of a tire under given conditions when trying to obtain a single performance or a plurality of performances that conflict with each other, and at the same time, makes the tire design and development highly efficient.
  • the present inventors conducted various studies to achieve the above-mentioned object, and to provide a method of designing a pneumatic tire capable of significantly reducing the cost of tires. Focusing on the application of the “optimization design method” used in different fields to the special field of tires, every attempt was made to establish a concrete tire design method.
  • the pneumatic tire designing method includes: (a) a tire basic model representing a tire sectional shape including an internal structure, an objective function representing a tire performance evaluation physical quantity, a tire sectional shape or a tire structure. (B) optimizing the objective function while considering the design variables and at least one of the tire cross-sectional shape, tire structure, physical quantity for performance evaluation, and tire dimensions. And (c) giving an optimal value of the objective function. Designing the tire based on the design variables.
  • the sensitivity of the objective function which is the ratio of the change amount of the objective function to the unit change amount of the design variable
  • the sensitivity of the constraint condition which is the ratio of the change amount of the constraint condition to the unit change amount of the design variable
  • a selection target group including a plurality of tire basic models representing a tire cross-sectional shape including an internal structure is determined, and a tire performance evaluation for each tire basic model in the selection target group is determined.
  • An objective function that represents a physical quantity, a design variable that determines a tire cross-sectional shape or tire structure, a constraint that restricts at least one of a tire cross-sectional shape, a tire structure, a physical quantity for performance evaluation and a tire dimension, and an objective function and An adaptive function that can be evaluated from constraint conditions is determined.
  • two tire basic models are selected from the selection target group based on the adaptive function, and design variables of each tire basic model are determined at a predetermined probability.
  • a new tire basic model by crossover and change some design variables of at least one tire basic model (mutation) At least one of generating a new tire basic model and calculating the objective function, constraints and adaptation function of the tire basic model with the design variables changed, and changing the tire basic model and design variables. Iteratively saves the tyre basic model and repeats until the number of tyre basic models stored reaches a predetermined number.It is determined whether the new group consisting of the predetermined number of stored tyre basic models satisfies the predetermined convergence condition. When the convergence condition is not satisfied, the new group is used as the selection target group until the selection target group satisfies a predetermined convergence condition, and the new convergence condition is satisfied.
  • the value of the design variable that gives the optimal value of the objective function is determined in consideration of the constraints among the predetermined number of the tyr basic models stored at that time.
  • the sensitivity of the objective function which is the ratio of the change of the objective function to the unit change of the design variable
  • the unit Based on the sensitivity of the constraint, which is the ratio of the amount of change of the constraint, the amount of change in the design variable that gives the optimal value of the objective function is predicted while considering the constraint, and the design variable is changed by an amount corresponding to the predicted amount.
  • the value of the objective function and the value of the constraint when the design variable is changed by an amount corresponding to the predicted amount are calculated, and the adaptive function is obtained from the value of the objective function and the value of the constraint to obtain the tire basic model and It is possible to store the tier basic model whose design variables are not changed, and to repeat until the stored tier basic model reaches a predetermined number.
  • the design variables include a force distribution line, a folded ply line, a line representing a tire outer shape, a line representing a tire crown shape, and a function representing the shape of at least one of a reinforcing material line, a bead filler gauge distribution, and a rubber chip. At least one of the following: tire gage distribution, side rubber gage distribution, tread rubber gage distribution, tread base rubber gage distribution, inner surface reinforcement rubber gage distribution, belt-to-belt rubber gage distribution, and belt-end rubber gage distribution.
  • a variable that represents the gauge distribution of the members, a variable that represents the structure of at least one belt section of the angle, width, cord type, and driving density of each belt layer, and the shape, position, number, and length of the block and sipe A variable representing the shape of at least one pattern of and At least one can be included.
  • a plurality of lines are determined from the functions representing the shape of at least one of the carcass line, the folded ply line, the line representing the tire outer shape, the line representing the tire crown shape, and the reinforcing material line as the design variables,
  • a mapping function of a function that represents the shape of the line and a constraint condition for preventing multiple lines from intersecting are obtained, and the mapping function is used by the mapping function.
  • a tire basic model representing a tire sectional shape including an internal structure, an objective function representing a physical quantity for evaluating tire performance, a design variable for determining the tire sectional shape or the tire structure, and a tire sectional shape Or, set restrictions that restrict the tire structure.
  • This basic tire model includes a line representing the outer shape of the tire, a line representing the shape of the tire crown, a belt line representing the belt inside the tire, a force line representing the force of the tire, and a folding line of the carcass ply inside the tire.
  • the basic model of the tire may use a method called a finite element method that divides the element into a plurality of elements, or may use an analytical method.
  • the objective function representing the physical quantity for evaluating the tire performance, it is necessary to use the physical quantity that governs the tire performance such as the tire circumferential belt tension and the lateral spring constant during air filling to improve steering stability. it can.
  • the design variables that determine the cross-sectional shape of the tire are functions that represent at least one of a carcass line, a folded ply line, a line that represents the outer shape of the tire, a line that represents the shape of the tire crown, and a reinforcement line that represents the line of various reinforcements.
  • the design variables that determine the tire structure include the bead filer gauge distribution, rubber chafer gauge distribution, side rubber gauge distribution, tread rubber gauge distribution, and tread base rubber gauge distribution.
  • a variable representing the gauge distribution of at least one of the tire rubber members, the gauge distribution of the inner surface reinforcing rubber, the gauge distribution of the interbelt rubber, and the gauge distribution of the beltend rubber, the belt layer angle, the belt layer width, the play height , Ply turnover, bead angle, width, position, Belts portion and a bead of quality, etc., may be used a variable representing the structure of the rhino de section.
  • variables representing the shape of the pattern such as the shape of the block, the angle of the block groove wall, the position of the sipe, the number, and the length can be used.
  • Constraints that limit the cross-sectional shape of the tire and the tire structure include the constraints on the carcass line peripheral value, the constraints on the upper and lower primary natural frequencies, the constraints on the belt layer angle, the belt layer width, the tire dimensions, and the panel size. There are constants, tire deformation, tire weight, stress, strain, strain energy, and restrictions on rolling resistance.
  • the objective function, design variables, and constraints are not limited to the above examples, and various ones can be determined according to the tire design purpose.
  • the value of the design variable that gives the optimal value of the objective function is determined while considering the constraints.
  • the sensitivity of the objective function which is the ratio of the change amount of the objective function to the unit change amount of the design variable
  • the sensitivity of the constraint condition which is the ratio of the change amount of the constraint condition to the unit change amount of the design variable.
  • a selection target group consisting of a plurality of tire basic models representing a tire cross-sectional shape including an internal structure is determined, and a tire performance evaluation physical quantity is calculated for each tire basic model of the selection target group.
  • the design variables that determine the tire cross-sectional shape or tire structure, the constraints that restrict at least one of the tire cross-sectional shape, tire structure, physical quantity for performance evaluation, and tire dimensions, and the objective function and constraints An adaptive function that can be evaluated is determined.
  • two tire basic models are selected from the selection target group based on the adaptive function, and the design variables of each tire basic model are crossed at a predetermined probability to form a new tire basic model.
  • Creating a tire basic model and changing at least one of the design variables of the tire basic model At least one of the generation of the tire basic model, and obtaining the objective function, constraints and adaptation function of the tire basic model in which the design variables are changed, and the tire basic model and the design of the tire base in which the design variables are not changed.
  • step (b) for the tire basic model with the design variables changed, the sensitivity of the objective function, which is the ratio of the change of the objective function to the unit change of the design variables, and the constraint on the unit change of the design variables
  • the change in the design variable that gives the optimal value of the objective function is predicted while considering the constraint, and the design variable is changed by an amount corresponding to the predicted amount.
  • the value of the objective function and the value of the constraint when the design variable is changed by an amount corresponding to the predicted amount are calculated, and the adaptive function is obtained from the value of the objective function and the value of the constraint to obtain the tire basic model. Further, it is more effective to store the tire basic model without changing the design variables and repeat until the stored tire basic model reaches a predetermined number. In this case, too, the value of the design variable when the value of the objective function is optimized in consideration of the constraints can be obtained.
  • a function that calculates the fitness for the tire model from the objective function and the constraint condition can be used as the adaptation function that can be evaluated from the objective function and the constraint condition.
  • the objective function, design variables, constraints and adaptation functions are not limited to the above examples, and various ones can be determined according to the purpose of tire design.
  • the crossover of the design variables of the tire basic model there is a method of exchanging the design variables of a part or a predetermined part or later of the design variables of the two selected tire models.
  • the design variables of the tire model there is a method of changing (mutating) the design variables at positions determined by a predetermined probability or the like.
  • step (c) the tire is designed by changing the tire basic model and the like based on the design variables giving the optimal value of the objective function.
  • the tire is designed based on the design variables mapped to the original function.
  • a selection target group consisting of a plurality of tire basic models representing a tire cross-sectional shape including an internal structure is determined, and a tire performance evaluation for each tire basic model in the selection target group is performed.
  • Objective functions representing physical quantities, design variables that determine the tire cross-sectional shape or tire structure, constraints that restrict at least one of the tire cross-sectional shape, tire structure, physical quantity for performance evaluation, and tire dimensions, and the objective function Step (b) can be composed of the following steps (1) to (4) when an adaptive function that can be evaluated from the constraints and the constraints is determined.
  • step (3) for the basic tire model with the design variables changed, the sensitivity of the objective function, which is the ratio of the change in the objective function to the unit change of the design variables, and the change in the constraint condition for the unit change of the design variable.
  • the amount of change in the design variable that gives the optimal value of the objective function is predicted based on the sensitivity of the constraint, which is the ratio of Calculate the value of the constraint condition when the value of the objective function and the design variable are changed by an amount corresponding to the predicted amount, calculate the adaptive function from the value of the objective function and the value of the constraint condition, and change the design variable.
  • Steps (1) to (3) are repeated until the new tire basic model and the tire basic model that do not change the design variables are saved and the stored tire basic models reach a predetermined number. It is effective to perform Te. In this case, too, the value of the design variable when the value of the objective function is optimized in consideration of the constraints is obtained.
  • design based on conventional trial and error • Unlike development, it is possible to some extent from the design of the best mode to the performance evaluation of the designed tire mainly by computer calculation As a result, significant efficiency gains can be achieved, and development costs can be reduced.
  • FIG. 1 is a schematic diagram of a personal computer used in an embodiment of the present invention.
  • FIG. 2 is a flowchart showing a processing routine of the first embodiment of the present invention.
  • FIG. 3 is a flowchart showing a routine for determining design variables in FIG. 1
  • Fig. 4 is a diagram showing the basic tire model
  • FIG. 5 is a diagram showing a state in which an imaginary line passing through the reference point P is drawn for each d on the tire basic model of FIG. 4,
  • FIG. 6 is a diagram showing the relationship between the selected node, the distance r i and the estimated angle, etc. for this node,
  • FIG. 7 is a flowchart showing a processing routine of the first embodiment of the present invention.
  • FIG. 8 is a flowchart showing a routine for determining design variables in FIG. 7,
  • FIG. 9 is a diagram showing a starting point S, an ending point E, an inflection point (:, and a maximum width position point M in the basic tire model,
  • FIG. 10 is a diagram showing the radius of curvature of an arc assumed between the start point S, end point E, inflection point C, and maximum width position point M,
  • Fig. 11 is a diagram showing the amount of deviation between the original node in the basic tire model and the optimal point determined by the least squares method.
  • FIG. 12 is a diagram showing rz coordinates
  • FIG. 13 is a diagram showing ⁇ V coordinates
  • FIG. 15 is a diagram showing the shape of the tire bead filler and its peripheral members
  • FIG. 16 is a diagram showing the belt layer line
  • FIG. 17 is a diagram showing that the line of the belt layer is approximated by a curve
  • Fig. 18 is a diagram showing the deterioration of element shape due to the movement of nodes.
  • FIG. 19 is a flowchart showing an element division processing routine
  • FIG. 20 is a diagram showing the process of element division, (a) is a diagram showing a model at the time of region definition, (b) is a diagram showing nodes generated in each region, (c) ) Indicates that a triangular element has been generated, (d) indicates that node smoothing has been performed,
  • FIG. 21 is a flowchart showing a Lagrangian interpolation processing routine.
  • FIG. 22 is a diagram showing that the crown portion is represented by a curve
  • FIG. 23 is a diagram showing that the crown is divided by the order of Lagrange interpolation
  • FIG. 24 is a diagram showing a normal line to the line of the crown portion
  • FIG. 25 is a diagram showing that the crown shape is represented by an arc
  • FIG. 26 is a diagram showing that the crown shape is represented by an arc with independent center coordinate values
  • FIG. 27 is a diagram showing that the pattern surface is divided into a lattice shape according to the Lagrange order
  • FIG. 28 is a diagram showing that the pattern surface is subdivided
  • FIG. 29 is a flowchart showing a processing routine of a tenth embodiment of the present invention.
  • FIG. 30 is a flowchart showing a crossover processing routine
  • FIG. 31 is a diagram showing a mountain-shaped mapping function
  • (a) is a diagram showing a continuous mountain-shaped mapping function
  • (b) is a diagram showing a linear mountain-shaped mapping function
  • FIG. 32 is a diagram showing a valley mapping function
  • (a) is a diagram showing a continuous valley mapping function
  • (b) is a diagram showing a linear valley mapping function
  • FIG. 33 is a flowchart showing a mutation processing routine.
  • FIG. 34 is a diagram showing a tire basic model of the eleventh embodiment, (a) is a diagram showing a tire basic model at the time of creation, and (b) is a diagram for explaining the existence of a belt.
  • Image diagram is a diagram showing a tire basic model of the eleventh embodiment, (a) is a diagram showing a tire basic model at the time of creation, and (b) is a diagram for explaining the existence of a belt.
  • FIG. 35 is a flowchart showing a crossover processing routine
  • FIG. 36 is a flowchart showing a mutation processing routine
  • FIG. 37 is a flowchart showing a processing routine according to the 12th embodiment of the present invention.
  • FIG. 38 is a diagram showing the modeling of a block
  • (a) is a diagram showing a block shape for modeling
  • (b) is a diagram showing an example of modeling a finite element model.
  • FIG. 39 is a flowchart showing a processing routine of the thirteenth embodiment of the present invention.
  • FIG. 40 is a flowchart showing a processing routine of a 14th embodiment of the present invention
  • FIG. 41 is an explanatory diagram for explaining heel and toe wear. BEST MODE FOR CARRYING OUT THE INVENTION
  • FIG. 1 shows an outline of a personal computer for implementing the pneumatic tire designing method of the present invention.
  • the personal computer includes a keyboard 10 for inputting data and the like, a computer main body 12 that satisfies constraints in accordance with a program stored in advance, and calculates design variables that optimize the objective function, for example, maximize or minimize the objective function, and It consists of a CRT 14 that displays the calculation results and other information of the computer body 12.
  • FIG. 2 shows a processing routine of the program of the first embodiment.
  • step 100 the tire cross-sectional shape in a natural equilibrium state is used as a reference shape, and the tire circumferential belt tension at the time of air filling is numerically and analytically determined using this reference shape as in the finite element method.
  • the model is modeled using a method that can be used to represent the cross-sectional shape of the tire, including the internal structure, and a basic tire model that is divided into multiple elements by mesh division.
  • the reference shape is not limited to the tire cross-sectional shape in a natural equilibrium state, and may be any shape.
  • modeling refers to converting the tire shape, structure, material, and pattern into numerical data in the form of input data to a computer program created based on numerical and analytical methods. Fig.
  • the basic tire model is divided into a plurality of elements by a plurality of normals NL of a force line CL and a plurality of elements NL 2 , NL 3 .
  • CL is the carcass line
  • OL is the line representing the outer shape of the tire
  • PL is the folded ply line
  • Bl and B2 are the lines representing the belt.
  • the basic tire model is divided into a plurality of elements by a plurality of normals NL of a force line CL and a plurality of elements NL 2 , NL 3 .
  • Normal Normal Therefore, it may be divided into a plurality of elements, or may be divided into an arbitrary shape such as a triangle depending on the design purpose.
  • an objective function representing a physical quantity for tire performance evaluation, constraints that restrict the tire cross-sectional shape, and design variables for determining the tire cross-sectional shape are determined.
  • the objective function OBJ and the constraint condition G are set as follows. It has established.
  • the perimeter value of the force-cass line can be calculated as the sum of the distances between the nodes of the force-cass line (intersection points between the force-cass line and the normal) within the range in which the tire shape is changed.
  • the shape of the carcass line which is a design variable, is determined by the Lagrange interpolation routine shown in FIG. 3 that approximates the curve.
  • a reference point P is set inside the tire as shown in FIG.
  • a next step 1 3 2 specifies the range of from the node in the vicinity of the belt end to the node q 2 in the vicinity of which is constrained to the rim as a range of changing the tire shape.
  • a straight line connecting the node q and the reference point P is used as a reference line, and a prospect angle 6 which is an angle between the reference line and a straight line connecting the node q 2 and the reference point P is calculated.
  • an angle increment d0 is calculated according to the following equation.
  • the user inputs the order of Lagrange interpolation in advance.
  • step 13 as shown in FIG. 5, a virtual line L,. L 2. L 3 ,... Is assumed for each angle increment d ⁇ ⁇ ⁇ with respect to the reference line, and the normal line closest to the virtual line is assumed.
  • Q 3 , ⁇ ⁇ ⁇ And criteria The distance r from the point P is r 2 , r 3.
  • Q 3. ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 2. ⁇ 3 ⁇ ⁇ --(Hereinafter, expressed as 0, where i 1, 2, ⁇ ⁇ ⁇ 1) is calculated.
  • the distance is set as a design variable.
  • step 104 of FIG. 2 the initial value 0 BJ of the objective function 0 BJ at the initial value ro of the design variable r, Calculate the initial value G 0 of 0 and the constraint condition G.
  • Step 106 the design variables r, are each changed by ⁇ , to change the tire basic model.
  • Step 1 0 8 the innermost node Q corresponding to the design variable, the distance r m between ⁇ innermost nodes and reference points other than that calculated by the following equation.
  • e m is the straight line connecting the said node and the reference point P, Mel in viewing angle from a reference line.
  • step 108 nodes other than the innermost node on the normal (nodes on the carcass line), namely, the line OL representing the tire outer shape, the folded back line PL, and the line B1, representing the belt,
  • the distance between the node on B2 and the innermost node Qi is calculated, and the calculated distance is added to the coordinates of the innermost node Q, to obtain the node coordinates other than the innermost node on the normal, and Determine the tire cross section shape after changing the variable ⁇ ⁇ ⁇ , that is, the tire correction model.
  • step 110 the objective function value OB J, after changing the design variables ⁇ ⁇ ⁇ , and the constraint value G, for the tire correction model obtained in step 108, are calculated.
  • the objective function for the unit change of the design variable is calculated according to the following formula in 2
  • the sensitivity d OBJ / dr of the objective function which is the ratio of the change in
  • the sensitivity d G / dri of the constraint which is the ratio of the change of the constraint to the unit change of the design variable
  • step 114 using the initial value of the objective function OBJ ⁇ , the initial value of the constraint condition Go, the initial value of the design variable r0, and the sensitivity, the objective function is satisfied while satisfying the constraint condition by mathematical programming. Predict the amount of change in the design variables to be maximized.
  • a tire correction model is determined in step 115 in the same manner as in step 108, and an objective function value is calculated.
  • step 116 the difference between the objective function value OBJ calculated in step 115 and the initial value OBJ0 of the objective function calculated in step 104 is compared with a threshold value previously input.
  • FIG. 7 shows a processing routine of the program according to the present embodiment. The parts corresponding to those in FIG.
  • step 162 After determining the objective function and the constraint conditions in the same manner as described above in step 160, the design variables are determined in step 162.
  • Fig. 8 shows the details of step 162.
  • step 150 as shown in Fig. 9, the starting point S and the ending point start from the nodes on the carcass line within the range in which the shape of the basic tire model is changed. Select E, inflection point C, and maximum width position point M, and input the data of these points.
  • step 152 one or more arcs are assumed between the start point S, end point E, inflection point C, and maximum width position point M as shown in FIG. Calculate the radius of curvature, R 2 , R 3 , R 4, and the expected angle, ⁇ 2 ⁇ , ⁇ 4 of the arc.
  • the design variables When the design variables are to be increased, the number of arcs may be increased.
  • step 154 the coordinates of the intersection of the arc and the normal are calculated, and as shown in Fig. 11, the deviation ⁇ between the intersection of the arc and the normal and the node on the carcass line (Represented by the general formula as £ i.)
  • step 156 the radius of curvature R i ( .P ., The expected angle ⁇ .P u of the arc where the ⁇ ⁇ , 2 is minimized by the mathematical programming is calculated, and the R, ( ..P t >, S, (. P t> is assumed to be design variables R,, ⁇ ,.
  • step 16 of FIG. 7 the node on the force-cass line is moved to the optimal point determined by R i, ⁇ i (. P t> obtained in step 15 6.
  • nodes other than the nodes on the force-cass line namely, the line ⁇ L representing the tire outer surface shape, the folded ply line PL, and the nodes on the lines B1 and B2 representing the belt and the nodes on the carcass line Is obtained by adding the calculated distance to the coordinates of the nodes on the carcass line to obtain the coordinates of the nodes other than the nodes on the force-cass line.
  • the cross-sectional shape is used as the sunset basic model, and the same as steps 104 to 106 described above.
  • the initial value OBJ 0 of the objective function and the initial value G 0 of the constraint condition at the initial value of the design variable Ro, 0 0 are calculated, and the design variable is changed by ⁇ ⁇ , and the tire correction model is changed.
  • Calculates the sensitivity for each design variable predicts the amount of change in the design variable that maximizes the objective function while satisfying the constraints using mathematical programming, and determines a Tier correction model for the predicted value of the design variable.
  • the objective function 0 BJ is calculated, and it is determined whether or not the value of the objective function has converged by comparing the difference with the initial value 0 BJ 0 with a threshold value. If the convergence has not occurred, steps 104 to 116 are repeatedly executed using the design variable values obtained in step 114 as initial values. When 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.
  • the shape of the force-cass line is adopted as a design variable.
  • the shape of a folded ply line, a line representing a tire outer surface shape, a line of a reinforcing material, or the like may be determined as a design variable.
  • the present invention is applied to the determination of the shape of a plurality of lines, and an attempt is made to optimize the conflicting performance.
  • the shape of the force line, the shape of the folded ply line, and the outer shape of the tire for improving the steering stability without impairing the riding comfort of the vehicle are determined. If the ride comfort of the vehicle is improved, the steering stability is deteriorated, so that the ride comfort and the steering stability of the vehicle are contrary to each other.
  • the procedure of this embodiment is substantially the same as that of the above embodiment, except that the objective function, the constraint conditions, and the design variable selection are different.
  • a lateral spring constant which is a physical quantity for improving steering stability, is used as an objective function
  • a condition that the upper and lower first-order eigenvalues, which are physical quantities that govern ride comfort, are used as constraint conditions is adopted.
  • the shape of the carcass line, the shape of the folded ply line, and the outer shape of the tire that maximize the lateral spring constant are determined.
  • the design variable r i is determined using the Lagrangian interpolation method described in the first embodiment.
  • design variables are determined as follows. First, as shown in Fig. 12, assuming a region RE, which has a carcass line as one side in the original rz coordinates, this region RE, is different from the rz coordinates ⁇ , 77, and the force-cass line is defined by the coordinates. Create a mapping function so that a square with one side is formed, and map the folded ply line and the line representing the outer shape of the tire. Similarly, assuming a region RE 2 that has the folded ply line as one side at the V 1 coordinate, a square having the folded blind line as one side is formed at the f 2 ⁇ ] 2 coordinate different from the ⁇ 1 m coordinate. Such a mapping function is created, and a line representing the outer shape of the tire is mapped.
  • mapping function is expressed as follows.
  • Table 1 below shows the results of actual trial production and tests of the sunseters obtained in the first and third examples.
  • three lines are used as design variables, but two lines or four or more lines may be used as design variables.
  • the shape determination of the plurality of lines shown in the third embodiment is applied to the shape determination of a tire bead filler and its peripheral members.
  • the shape of the bead filter and the shape of the rubber zipper X-filer which reduce the rolling resistance without impairing the durability of the bead portion are determined.
  • the procedure of this embodiment is almost the same as that of the third embodiment, except that the objective function, the constraint conditions, and the design variables are different.
  • the objective function and the constraint conditions are determined as follows.
  • the principal strain under the load at the end of the folded ply is within + 3% of the initial structure.
  • the design variable is a line that defines the outer shape (gauge distribution) of the bead filter 20 as shown in Fig. 15 This is the boundary line KL between BLa, BLb and rubber chain 22 and side rubber 24.
  • the determination of the shape of the plurality of lines shown in the third embodiment is applied to the determination of the thickness between the belt layers in the belt portion.
  • the gauge distribution between the belt eyebrows that reduces the rolling resistance without impairing the durability of the belt portion is determined.
  • the procedure of this embodiment is almost the same as that of the third embodiment, except that the objective function, the constraint condition, and the design function are different.
  • the objective function and the constraint condition are determined as follows.
  • Constraints The principal strain under load at the end of the belt end ply is within + 3% of the initial structure.
  • the design variables are the lines B1, B2, B3, B representing each belt layer as shown in Fig. 16. By changing these lines B1 to B4, the gauge distribution between the belt layers that optimizes the objective function is determined.
  • the Lagrange interpolation method and the circular interpolation method were used as the method of representing the line shape. Interpolation may be performed using a spline curve, B-spline curve, Bezier curve, NURBS (weighted B-spline), or the like. Table 2 below shows the results of actual trial production and testing of a tire obtained by combining the above-described fourth and fifth embodiments. [Table 2]
  • the present invention is applied to a belt structure, and the durability is improved by minimizing the concentration of strain generated in the belt portion without increasing the weight of the belt portion.
  • the procedure of the present embodiment is substantially the same as the procedure of the first embodiment shown in FIG. 2. Only the objective function, the constraint conditions, and the determination of the design function in step 102 are different. And the constraints are defined as follows.
  • Total belt weight is within + 1% of the total weight of the initial structure
  • the design variables are the angle, the driving (for example, the number of drivings and the driving strength) and the width of each belt layer.
  • the values of the angle and the driving may be directly used as design variables, and the design variables of the width are determined as follows.
  • the lines of each belt are approximated by, for example, a cubic function using the coordinate values of the nodes constituting each belt layer.
  • lines B 1 and B 2 are approximated by a cubic function.
  • each node constituting each belt layer is set to a value of the magnification determined for each belt layer, and the z-coordinate in the basic model of each node constituting the layer is obtained.
  • the value multiplied by the value is assumed to be the z coordinate, and this value is substituted for a function approximating the layer.
  • step 115 of FIG. 2 when modeling was performed by the finite element method, the nodes constituting each belt layer moved as the width of each belt layer changed. As shown in the circle 32, there is a problem that the element shape is significantly deteriorated and the calculation cannot be continued.
  • the routine shown in FIG. 19, which automatically divides the element from the given node position is used. Just add it.
  • step 180 of FIG. 19 data is input to define an area including each belt layer.
  • the areas AR1, AR2 and AR3 can be defined as shown in FIG. 20 (a).
  • nodes are generated in each area, and in the next step 184, the nodes are generated as shown in Fig. 20 (c).
  • a triangular element is generated by connecting all nodes including the generated nodes.
  • the nodes in the area are smoothed by moving the positions of the nodes, and in the next step 188, renumbering is performed. Performs bandwidth reduction processing for internal processing in the finite element method such as, and ends this routine.
  • a cubic function is used as a method of expressing the shape of the belt line.
  • Lagrangian interpolation, circular interpolation, or the like used in the first and second embodiments are used.
  • a spline curve, a B-spline curve, a Veggie curve, NURB S, etc. may be used. The results of actual trial production and testing of the tire obtained in the sixth embodiment described above are shown in Table 3 below. .
  • the present invention is applied to the determination of the shape of the tire crown portion.
  • the pressure distribution in the contact region is made uniform without changing the shape of the contact region between the tire and the ground, and the wear performance is improved. is there.
  • the procedure of this embodiment is substantially the same as that of the first embodiment shown in FIG. 2, except that the objective function, the constraint conditions, and the design variable of step 102 are different.
  • the objective function and the constraints are defined as follows.
  • the crown shape which is a design variable, is determined by the Lagrangian interpolation routine shown in FIG.
  • step 170 of this Lagrange interpolation routine a range in which the shape of the crown portion is changed and the number of curves for approximating the range are specified.
  • the crown shape is two from the node Q3 at the center of the tire to the node Q4 at the groove end, and from the other node Q5 at the groove end to the node q6 near the belt end.
  • the curve may be, for example, one from node q3 to node Q6, or nodes Q5 and q6 between nodes Q3 and q4.
  • the space may be further divided and represented by three or more curves.
  • each range is equally divided into the degree of Lagrange interpolation input beforehand, and a node closest to the equal division point is selected.
  • the tire is equally divided by dividing lines Ls1, Ls2, Ls3, Ls4, and Ls5, and each of these dividing lines Ls1 to Ls5 The nearest nodes Qsl, qs2, qs3, qs4, and Qs5 are selected.
  • step 4 the normal to the crown line at the selected node position is calculated, and this is defined as the line on which each node moves.
  • the nodes q3, q4, q5, and Q6 representing the range of the curved line and the nodes Qsl, qs2, and qs3 selected in step 172 are obtained.
  • qs 4 and qs 5 are calculated to determine the lines L pl to L p 8.
  • the amount of movement of each node on each line from the initial position is selected as a design variable, and this routine ends.
  • the amount of movement of the nodes on the crown other than the selected nodes is determined by the Lagrangian trap, as described in the first embodiment.
  • an eighth embodiment of the present invention will be described.
  • the present invention is applied to the determination of the shape of the tire crown as in the seventh embodiment, but differs from the seventh embodiment only in the design variables.
  • the procedure of this embodiment is the same as that of the second embodiment shown in FIG. However, the only difference is that the objective functions, constraints and design variables in steps 160 and 162 are determined.
  • the objective functions and constraints are defined as follows.
  • the contact length in the tire circumferential direction at the center of the crown and the belt end is within 5% of the contact length in the initial shape.
  • the range of the crown portion specified in advance is approximated by a plurality of arcs.
  • the range of the crown portion that can be approximated from the node Q7 at the center of the tire to the node q8 near the belt end in this case, the range CR1, CR2, CR
  • the crown shape is represented by three arcs for three regions.
  • the radius can be approximated by an angle 1 with R 1
  • the radius can be approximated by an angle 02 with a radius of R 2
  • a radius of R 3 with an angle 03 Can be approximated by
  • Fig. 25 shows that the crown shape is represented by three arcs in the range from the node q7 at the center of the tire to the node q8 near the end of the belt.
  • the number of arcs may be 1 to 2 or 4 or more.
  • each arc does not necessarily need to be smoothly continuous, and the center coordinate value of each arc may be treated as an independent variable as shown in FIG.
  • the radius is approximated by R 1 ′ and the angle ⁇ ⁇ ⁇ from the coordinate point 01 (r 1, z 1), and for the range CR 2, the coordinate point 02 (r 2, z 2)
  • the radius is approximated by R2 and the angle 62 '
  • the range CR3 is approximated by the radius R3 and the angle 03' from the coordinate point 03 (r3, z3).
  • the radius R i and the expected angle 0 i of each arc and in the case of the example shown in FIG. 26, the coordinate values r i and z i of the center of each arc are additionally set as design variables.
  • the method of expressing the shape of the crown portion is represented by the Lang ⁇ interpolation and the circular arc interpolation.
  • a spline curve, a B-spline curve, a Bezier curve, NURBS, or the like may be used.
  • the present invention is applied to the determination of the shape of the tire pattern surface, in which the pressure distribution when the pattern comes into contact with the ground is made uniform and the wear performance is improved.
  • the procedure of this embodiment is substantially the same as that of the first embodiment shown in FIG. 2, except that the objective function, the constraint conditions and the design variable determination in step 102 are different.
  • the functions and constraints are defined as follows.
  • the pattern surface which is a design variable, as shown in FIG. 27, the pattern surface is divided into a lattice shape in accordance with the Lagrangian interpolation order input in advance.
  • the coordinate values r i in the pattern thickness direction of each point (16 points in the example of FIG. 27) obtained in this manner are used as design variables.
  • the finite element method it is necessary to further subdivide the pattern surface.
  • the coordinates of the nodes other than the nodes corresponding to the design variables on the pattern surface may be calculated by the following equation.
  • N m (7j) ⁇ (77-7? I) ⁇ (V m — n))
  • Lagrangian interpolation is used as a method of representing the shape of the pattern surface.
  • parameters for controlling these may be directly selected as design variables.
  • Each of the above embodiments need not be performed independently.
  • the determination of the ply line in the first embodiment and the determination of the belt structure in the sixth embodiment are simultaneously performed.
  • Each of them may be combined, for example, the shape of the crown part in the seventh embodiment is determined using the ply line determined in the above.
  • FIG. 29 shows the processing of the program of the tenth embodiment. It shows a routine.
  • step 200 the N tire cross-sectional shapes are modeled by a method that can numerically and analytically determine the circumferential belt tension at the time of air filling, such as the finite element method.
  • the tire basic model including the internal structure.
  • modeling refers to digitizing the shape, structure, and material patterns into an input data format for a combi- nation program created based on a numerical analysis method.
  • N is input by the user in advance.
  • the basic tire model used in this embodiment is the same as that shown in FIG. 4 of the first embodiment.
  • the tire basic model may be divided into a plurality of elements by a plurality of normals of a line representing the outer shape of the tire and a plurality of normals of the folded ply line. It may be divided into any shape such as a triangle depending on the purpose.
  • the objective function representing the physical quantity for tire performance evaluation, the constraint conditions for restricting the tire cross-sectional shape, and the design variables for determining the tire cross-sectional shapes of the N tire models are determined.
  • the objective function ⁇ BJ and the constraint condition G are set as follows. It is determined as follows.
  • the shape of the carcass line which is a design variable, is determined for each of the N tire models by the Lagrangian interpolation routine of FIG. 3 that approximates a curve described in the first embodiment. Since the Lagrangian interpolation routine is the same as that of the first embodiment, the description is omitted.
  • the objective function OBJ and the constraint condition of each of the design variables ri of the ⁇ tire models are calculated.
  • the objective function ⁇ BJ j and the constraint G of each of the N tire models obtained in step 204 are used to calculate the adaptive function F of each of the N tire models. Is calculated according to the following equation (4).
  • the value (fitness) based on the adaptation function increases as the belt tension increases.
  • the user inputs c and a in advance.
  • the probability P at which an individual with N tire models is selected by selection, is expressed by the following equation:
  • F is the adaptive function of an individual ⁇ in N tire models
  • the fitness proportional strategy was used as the selection method.
  • next step 210 it is determined whether or not the two selected tire models are crossed according to the probability T input by the user in advance.
  • the term “crossover” refers to exchanging some of the elements of the two tire models. If a negative judgment is made that no crossover is to be performed, the process proceeds to step 2 16 in step 2 12 with the current two tire models as they are. On the other hand, when crossing is determined in the affirmative judgment, two tire models are crossed as described later in step 214.
  • the intersection of the two tire models is performed by the intersection routine shown in FIG.
  • the two tire models selected in step 208 are referred to as tire models a and b, and the design variables of the respective tire models a and b are represented by design variable vectors including a list.
  • the design variable base-vector of the tire model a V r a two (3, 1 "2 3, , 3, ⁇ ⁇ ⁇ , * n '.' -!, a)
  • the crossover location i of the design variable vector of the tire models a and b is determined.
  • step 2 5 crossing the determined tire model a, design variables r of b, a, with respect to ri b, determine the distance d according to the following equation.
  • r, a with a possible range minimum value B and and maximum value B u of ri b, determine the normalized distance d 'according to the following equation.
  • step 2 5 6 in order to distribute the value of the normalized distance d 'appropriately, Using the mountain-shaped mapping function Z (x) (0 ⁇ x ⁇ 1, 0 ⁇ Z (x) ⁇ 0.5) as shown in Figs. (A) and (b), the function value Z Ask for ab .
  • Step 2 6 0 is a sequence of new design variables design variable vector V r 'a, V r' b are as follows Desired.
  • V r ' 3 (r, 3 , r 2 a ⁇ ⁇ ,' a , r 1 + 1 b , ⁇ ⁇ ⁇ , r,-, b )
  • V r 'b (ri b sr 2 b, ⁇ ⁇ ⁇ ri' b, r 1 + 1 a, ⁇ ⁇ ⁇ , r n -.
  • the minimum value BL and the maximum value of r, possible range B u is input by the user in advance.
  • the mapping function Z (X) may be a valley-type function as shown in FIGS. 32 (a) and (b). In the above example, there is only one crossover point i. However, as shown in the genetic algorithm (edited by Hiroaki Kitano), there are multiple crossover points i. Uniform crossover or the like may be used.
  • step 216 of FIG. 29 it is determined whether or not to mutate with the probability S input by the user in advance.
  • This mutation means that a part of the design variables is minutely changed, in order to increase the accuracy of including a population that can be an optimal design variable.
  • Step 226 if the mutation is not made in the negative judgment, in Step 226, the current two tire models are used and the process proceeds to the next Step 222. If the mutation is suddenly made in the affirmative judgment, mutation processing is performed in the following step 220 as follows.
  • step 2 62 a random number is generated, and the location i of the mutation is determined by the random number.
  • step 2 6 4 the distance d '
  • the IS meter variable vector V r ' is as follows.
  • V r ( ⁇ ⁇ , r ⁇ ⁇ r, r i + ⁇ ⁇ ⁇ , r n -In this way, for two newly generated diamond models, the value of the objective function and the value of the constraint Is calculated in step 222 of Fig. 29.
  • step 222 adaptively using the obtained objective function value and constraint value using equation (4) in the same manner as in the above embodiment. Operate a function.
  • next step 226, save the above two tire models.
  • next step 2208 it is determined whether or not the number of tire models stored in step 226 has reached N, and if not, step 20 is performed until N is reached. Repeat steps 2 2 8 from 8.
  • the convergence is determined in step 230, and if not, the N tire models are stored in step 22 And repeat steps 2 to 8 from step 2.
  • the constraint condition is determined using the values of the design variables of the tire model that maximize the value of the objective function while substantially satisfying the constraint conditions among the N tire models.
  • the value of the design variable that maximizes the objective function while substantially satisfying the condition is set, and the shape of the diamond is determined using the value of the design variable in step 232.
  • the convergence judgment in step 230 is regarded as convergence if any of the following conditions is satisfied.
  • the user inputs M, q, and p in advance.
  • Table 5 shows the results of actual trial production and testing of the tire obtained in the above-described 10th Example. Five ]
  • Test method Filling evaluation by actual vehicle driving test
  • the present invention is applied to a belt structure, and the performance that conflicts with each other is optimized.
  • a belt structure that improves steering stability without impairing the riding comfort of the vehicle is determined.
  • the procedure of this embodiment is almost the same as the processing routine of the program of the tenth embodiment shown in FIG. 29, except that the objective function, the constraint conditions and the design variable selection, the crossover method, and the sudden variation method are different. different.
  • the lateral spring constant which is a physical quantity for improving steering stability
  • the vertical spring constant which is the physical quantity that governs ride comfort
  • the belt structure that maximizes the lateral spring constant is determined.
  • crossover is performed based on the processing routine shown in FIG. First, it is assumed that the a-th and b-th tire models of the N tire models are crossed.
  • a design variable vector is used to list the design variables for each tire model.
  • VB a (b 1 3, b 2 a ⁇ ⁇ ⁇ b, b,. ⁇ a. ⁇ ⁇ ⁇ B m a)
  • VB b (b, b, b 2 b, ⁇ ⁇ ⁇ bi b, b i + 1 b, ⁇ ⁇ ⁇ b m b) to.
  • step 272 a random number is generated, and a crossover location i is determined according to the generated random number.
  • the sequence of design variables which are design variable vectors, is changed as follows, and the design variable vectors VB a ′ and VB b ′ are obtained as new design variable sequences.
  • VB a ' (b, a , b 2 a, ... ⁇ b, -! A, bi b , b i +1 b, ... ⁇ b m b)
  • VB b' (b, b, b 2 b, ... ⁇ b b , b, a , b 1+ , a , ⁇ ra b ra a )
  • two new tires are selected according to the determined design variable vectors VB a ′ and VB b ′. Generate a model.
  • the number of crossover points i is one.
  • a multipoint crossover, a uniform crossover, or the like as shown in a genetic algorithm (edited by Hiroaki Kitano) may be used. Mutation is performed based on the processing routine shown in FIG.
  • the design variable vector shows the sequence of design variables for the basic tire model.
  • VB (b,, b 2 , ... b,, b, bt10,, ... b m )
  • step 278 a random number is generated, and the location i of the mutation is determined according to the random number. .
  • a design variable vector VB ' is generated as a new sequence of design variables by changing the design variables b, at the mutation location i as follows.
  • the new design variable vector VB ' Generate a tire model.
  • the belt structure is adopted as a design variable, but other reinforcing materials and the like may be determined as design variables. Also, the present invention can be applied to whether or not a reinforcing material is put into a bead portion.
  • Table 6 below shows the results of actual trial production and tests of the tire obtained in the eleventh example.
  • Test method Feeling evaluation by actual vehicle driving test
  • the design variable that gives the optimal value of the objective function that satisfies the constraint condition is obtained, and the tire is designed from this design variable, so that the best mode tire design can be performed.
  • This has the effect.
  • a 12th embodiment will be described.
  • the present invention is applied to the determination of the block shape of a block-formed tire.
  • the uneven wear performance is improved by reducing the difference in block stiffness between the blocks and making the block stiffness in each direction uniform in one block.
  • FIG. 37 shows a program processing routine of the present embodiment.
  • step 300 the shape of the block 42 is changed as shown in FIG.
  • Figure 38 (a) by the node coordinates of the block 42, the node coordinates of the sipe 44, the angle of the groove of the block, the depth of the groove, and the sipe. It is defined by the number of slabs, the width of the sipe, the sipe embedding angle, and the sipe depth, and is modeled so that it can be analyzed by the finite element method and analytical methods.
  • Figure 38 (b) shows an example of modeling a finite element model.
  • the constraints that restrict the objective function representing the physical quantity for performance evaluation of the block the area of the block, the length of the sipe, the length of the sipe, the depth of the sipe, the number of the sipe, etc.
  • design variables related to the node coordinates of the block, the node coordinates of the sipe, the groove angle of each side of the block, the groove depth of each side of the block, the sipe width, the sipe embedding angle, and the sipe depth are related to the node coordinates of the block, the node coordinates of the sipe, the groove angle of each side of the block, the groove depth of each side of the block, the sipe width, the sipe embedding angle, and the sipe depth.
  • an objective function, constraints, and design variables are defined as follows in order to design a pattern in which the block rigidity in each direction is uniform in one block. ing.
  • Constraint G 'The change in sipe length is the sipe length in the initial model.
  • sipe node coordinates are within the shape enclosed by the block node coordinates
  • the distance between the sipe node coordinates and a line representing the outside of the block is 2 mm or more
  • the initial value X of the design variable X is set.
  • initial value G of the constraint G. Is calculated.
  • design variables X are changed by ⁇ ⁇ , respectively, to change the basic model of the block. Change.
  • the block including the change in the design variables implemented in step 306 is modeled.
  • step 310 the value of the objective function ⁇ BJ, after the design variable is changed by ⁇ ⁇ and the value of the constraint condition G i are calculated for the block correction model obtained in step 308, and the step 3 12
  • the sensitivity of the objective function, d 0 BJ d X which is the ratio of the change of the objective function to the unit change of the design variable, and the ratio of the change of the constraint condition to the unit change of the design variable,
  • the sensitivity d GZ d X, of a certain constraint condition is calculated for each design variable.
  • this sensitivity it is possible to predict how much the value of the objective function and the value of the constraint condition change when the design variable is changed by ⁇ ⁇ ⁇ . Note that this sensitivity may be obtained analytically depending on the method used for modeling the pro- cess or the nature of the design variables, and in that case, the calculation in step 312 is not required.
  • the initial value of the objective function is OBJ.
  • the initial G of the constraint The initial value of the design variable X.
  • the sensitivity we predict the amount of change in the design variables that maximizes the objective function while satisfying the constraints by mathematical programming.
  • a block correction model is determined in step 315 in the same manner as in step 308, and the objective function value is calculated.
  • step 3 14 Using the design variable values obtained in the above as initial values, steps 304 to 316 are repeatedly executed. When it is determined that the value of the objective function has converged, the objective function is maximized while satisfying the constraints with the value of the design variable at this time. In step 320, the shape of the block is determined using the values of the design variables.
  • Table 7 shows the results of actual trial production and abrasion tests of the tire obtained in the above-mentioned 12th embodiment.
  • Test method 500 km on test course
  • the H & T height (the height of the heel 'and') is a phenomenon in which the tire is worn in a saw blade shape in the circumferential direction. This H & T height is shown in FIG. Next, a thirteenth embodiment will be described. This embodiment is a combination of the first embodiment and the tenth embodiment. Note that this embodiment has substantially the same configuration as the above-described embodiment. The same parts are denoted by the same reference numerals, and detailed description is omitted.
  • the objective function and the constraint condition were calculated based on the design plan obtained by the crossover and mutation in step 118.
  • the method of the first embodiment and the method of the tenth embodiment are performed by using the processing of steps 104 to 116 of the first embodiment as the arithmetic processing of step 222 of the tenth embodiment. The combination can solve the above problems.
  • FIG. 39 shows a program processing routine of the present embodiment. Steps 200 to 220 are the same as those in the tenth embodiment, and a description thereof will be omitted.
  • step 340 the objective function and the constraint condition are calculated using the two design plans obtained as described above as initial design plans.
  • the next steps 106 to 116 are repeatedly executed until the value of the objective function converges, in the same manner as the processing of the first embodiment.
  • an adaptive function is calculated from the value of the objective function and the value of the constraint obtained in the next step 224, and the above two values are calculated in the next step 226. Save the tire model.
  • Repeat steps 2208 to 2228 until the number of tire models saved in step 226 reaches N, and if it reaches N, repeat step 2330 with the above.
  • the convergence judgment is performed, and when convergence is achieved, the constraint condition is substantially satisfied with the design variable value of the tire model in which the value of the objective function is maximized while substantially satisfying the constraint condition among the N tire models.
  • the value of the design variable that maximizes the objective function is determined, and the shape of the diamond is determined in step 232.
  • FIG. 40 shows a program processing routine of the present embodiment.
  • a predetermined number M of shape data specified in advance is read from the tire shape data 'base.
  • the tire shape database expresses the force-cass line of the tire using the Lagrangian intercept of equation (2), the method using the arc shown in Fig. 10, the spline curve, the page II curve, NURBS (weighted B spline), etc. It was done.
  • the carcass line F (X) of the tire is expressed by the following equation (5).
  • a plurality of normals, NL 2 , NL 3,... are drawn on the carcass line, and nodes are arranged on the lines in consideration of the thickness according to the tire member. Then, a plurality of elements are divided using the nodes.
  • the objective function representing the physical quantity for evaluating the tire performance, the constraint conditions for restricting the tire cross-sectional shape, and the design variables for determining the tire cross-sectional shape are determined.
  • the shape of the carcass line, which is a design variable, is determined by the constant a i in Eq. (5).
  • the design variables a are changed by ⁇ a, respectively, to change the tire basic model.
  • the carcass line of the tire after changing only the design variable mu a is obtained from equation (5), and a tire correction model is determined according to the procedure of step 402. Since step 110 and subsequent steps are the same as above, the description is omitted.
  • Test method Filling evaluation by actual vehicle driving test

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Tires In General (AREA)

Description

明細書 空気入りタイヤの設計方法 技術分野
本発明は、 空気入りタイヤの設計方法に関するもので、 特に、 タイヤの単一目 的性能、 二律背反性能等を達成するタイヤの構造、 形状等の設計開発を効率的に かつ容易にし、 しかもタイヤのベストな構造、 形状を求めかつコスト 'パーフォ —マンスの高いタイヤを設計することができる空気入りタイヤの設計方法に関す る。 背景技術
従来のタィャ設計方法は、 実験及び計算機を用 、た数値実験の繰り返しによる 経験則から成り立つていた。 このため、 開発に必要な試作 '試験の件数が膨大な ものとなり、 開発コストがアップし、 開発期間もなかなか短縮できなかった。
( 1 ) タイやの力一カス形状設計に関する背景技術
一般にタイヤのカーカス形状設計方法としては大別して下記の 2つの設計方法 があり、 しかもそのいずれか一方の設計方法で殆どのタイヤが設計開発され、 巿 販されているのが現状である。
タイヤの力一カス形状設計方法の 1つは自然平衡形状と称される設計方法で、 特にタイヤをリム組し、 空気を充塡してタイヤ内圧を正規内圧にしたとき、 タイ ャカーカスコードに作用する張力が一定となるように、 カーカスの断面形状 (夕 ィャの回転軸を含む平面で切断したときの断面形状) を決定するものでもある。 すなわち、 タイヤを金型で加硫硬化するときの力"カスの断面形状を正規内圧時 のカーカスの断面形状と実質上等しい形状とすることにより、 空気充塡時にタイ ャの種々捕強層の積層物間または積層内に部分的に発生する内部歪等による故障 発生を防止するものである。 他の 1つは非自然平衡形状と称される設計方法で、 代表的な非自然平銜形状設 計方法としては、 出願人が先に提案した 「R C〇T」 と称される商標で既に巿 販され、 コマーシャル的に成功をおさめている設計方法がある。 この非自然平衡 形状設計方法は前記自然平衡形状設計方法と異なり、 逆にカーカスの張力をコン トロールする設計方法であり、 空気充塡による釣合い形状へ移行しょうとする力 を利用し、 部分的に張力をコントロールすることにより、 所定性能を得ようとす るものである。
( 2 ) タイヤのクラウン部形状設計に関する背景技術
タイヤのクラウン部形状は、 タイヤの回転軸を含む断面において、 数個の円弧 によって設計されている。 円弧の値を決定する方法は、 ①円弧の値とタイヤ性能 のデータベースから決定する方法、 ②数個のモールドを作成し、 このモールドに よるタイヤを試作 ·評価したデータから決定する方法、 ③数値実験を数多く行い 決定する方法等であった。 このため、 開発効率が良くなかった。
( 3 ) ベルト構造設計に関する背景技術
ベルト構造は、 ベルト幅、 ベルト角度、 コード種類、 打ち込み密度、 ベルト間 ゴムゲージ分布で決定される。 ベルト構造設計の方法は、 ①実験、 数値実験によ つてベルト端部の歪を計測する方法、 ②ベルト全体の剛性を計測する方法、 ③べ ルト構造を変更して室内評価、 実車評価を実施する方法等によって、 設計案を抽 出することであった。 このため、 評価期間が長い等の問題があった。
( 4 ) パターン設計に関する背景技術
パターン設計は、 自由度が大きいので、 基本パターン案をタイヤにグルービン グしたり実際にモールドを作成した後にタィャを試作して実車評価し、 実車で生 じた問題は基本パターン案を微修正することによって解決して最終的パターンに 仕上げていた。 このため、 パターン設計は、 タイヤ形状、 構造設計に比較して最 も工数を要する分野となっていた。
( 5 ) その他の設計に関する背景技術
これまで述べてきた設計方法と同様の問題を、 タイヤゴム部材のゲージ分布 ( ビードフイラ一のゲージ分布、 ゴムチェ一ファーのゲージ分布、 サイ ドゴムのゲ —ジ分布、 トレッ ドゴムのゲージ分布、 トレッ ドベースゴムのゲージ分布、 内面 補強ゴムのゲージ分布、 ベルトエンドゴムのゲージ分布) の設計についても有し ている。 発明の開示
しかしながら従来の設計方法を用いたタイヤの設計開発はある性能について目 標値を定め、 この目標値をクリアすれば一応終了とされ、 与えられた資源でべス ト、 言いかえるならばその最良の性能を得ると言う考え方のものではなかった。 また、 二律背反する性能を設計するものでなく、 そのべストな形状、 構造を決定 するものでもなかった。 その上いずれの設計法も開発が試作 ·試験の試行錯誤の 繰返しで行われるため、 非常に非効率でコスト 'パフォーマンスが高い等の問題 があった。
従つて本発明は、 ある単一の性能または二律背反する複数の性能を得ようとす るとき、 与えられた条件でタイヤのベストモードを設計することができると共に 、 タイヤの設計 '開発を高効率化し、 タイヤのコスト 'パーフォーマンスを著し く低下することができる空気入りタイャの設計方法を提供することを目的とする 上記目的を達成するために、 本発明者等は種々検討を加えた結果、 異分野に利 用されている 「最適化設計手法」 をタイヤと言う特殊分野に応用することに着目 し、 あらゆる検討を試み、 具体的にそれをタイヤ設計方法として確立したもので ある。 具体的には、 第 1の発明の空気入りタイヤの設計方法は、 (a ) 内部構造 を含むタイャ断面形状を表すタイャ基本モデル、 タィャ性能評価用物理量を表す 目的関数、 タイヤ断面形状またはタイヤ構造を決定する設計変数、 及びタイヤ断 面形状、 タイヤ構造、 性能評価用物理量及びタイヤ寸度の少なくとも 1つを制約 する制約条件を定めるステップ、 (b ) 制約条件を考慮しながら目的関数の最適 値を与える設計変数の値を求めるステップ、 及び (c ) 目的関数の最適値を与え る設計変数に基づいてタイヤを設計するステップ、 を含んでいる。
このステップ (b ) では、 設計変数の単位変化量に対する目的関数の変化量の 割合である目的関数の感度及び設計変数の単位変化量に対する制約条件の変化量 の割合である制約条件の感度に基づいて制約条件を考慮しながら目的関数の最適 値を与える設計変数の変化量を予測すると共に、 設計変数を予測量に相当する量 変化させたときの目的関数の値及び設計変数を予測量に相当する量変化させたと きの制約条件の値を演算し、 予測値と演算値とに基づいて、 制約条件を考慮しな がら目的関数の最適値を与える設計変数の値を求めることができる。
また、 本発明者等は種々検討を加えた結果、 異分野に利用されている 「遺伝的 アルゴリズム手法」 をタイヤと言う特殊分野に応用することに着目し、 あらゆる 検討を試み、 具体的にそれをタイヤ設計方法として確立したものである。 具体的 には、 前記ステップ (a ) では、 内部構造を含むタイヤ断面形状を表す複数個の タイヤ基本モデルからなる選択対象集団を定め、 該選択対象集団の各タイヤ基本 モデルについて、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断面形状また はタイヤ構造を決定する設計変数、 タイヤ断面形状、 タイヤ構造、 性能評価用物 理量及びタイヤ寸度の少なくとも 1つを制約する制約条件、 及び目的関数及び制 約条件から評価できる適応関数を定め、 前記ステップ (b ) では、 適応関数に基 づいて前記選択対象集団から 2つのタイヤ基本モデルを選択し、 所定の確率で各 タィャ基本モデルの設計変数を交叉させて新規のタイャ基本モデルを生成するこ と及び少なくとも一方のタイヤ基本モデルの設計変数の一部を変更 (突然変異) させて新規のタイヤ基本モデルを生成することの少なくとも一方を行い、 設計変 数を変化させたタィャ基本モデルの目的関数、 制約条件及び適応関数を求めて該 夕ィャ基本モデル及び設計変数を変化させなかつたタイャ基本モデルを保存しか つ保存したタイャ基本モデルが所定数になるまで繰り返し、 保存した所定数の夕 ィャ基本モデルからなる新規集団が所定の収束条件を満たすか否かを判断し、 収 束条件を満たさないときには該新規集団を前記選択対象集団として該選択対象集 団が所定の収束条件を満たすまで繰り返すと共に、 該所定の収束条件を満たした ときに保存した所定数のタィャ基本モデルのなかで制約条件を考慮しながら目的 関数の最適値を与える設計変数の値を求める。
このステップ (b ) では、 設計変数を変化させたタイヤ基本モデルについて、 設計変数の単位変化量に対す,る目的関数の変化量の割合である目的関数の感度及 び設計変数の単位変化量に対する制約条件の変化量の割合である制約条件の感度 に基づいて制約条件を考慮しながら目的関数の最適値を与える設計変数の変化量 を予測すると共に、 設計変数を予測量に相当する量変化させたときの目的関数の 値及び設計変数を予測量に相当する量変化させたときの制約条件の値を演算し、 目的関数の値及び制約条件の値から適応関数を求めて該タイヤ基本モデル及び設 計変数を変化させなかったタイャ基本モデルを保存しかつ保存したタイャ基本モ デルが所定数になるまで繰り返すことができる。
前記設計変数は、 力一カスライン、 折り返しプライライン、 タイヤ外面形状を 表すライン、 タイヤクラウン形状を表すライン、 及び補強材ラインの少なくとも 1つのラインの形状を表す関数と、 ビードフイラ一のゲージ分布、 ゴムチヱーフ ァ一のゲージ分布、 サイ ドゴムのゲージ分布、 トレッ ドゴムのゲージ分布、 トレ ッ ドベースゴムのゲージ分布、 内面補強ゴムのゲージ分布、 ベルト間ゴムのゲー ジ分布、 及びベルトェンドゴムのゲージ分布の少なくとも 1つのタイヤゴム部材 のゲージ分布を表す変数と、 各ベルト層の角度、 幅、 コード種類、 及び打ち込み 密度の少なくとも 1つのベルト部の構造を表す変数と、 プロックの形状及びサイ プの位置、 本数、 及び長さの少なくとも 1つのパターンの形状を表す変数と、 の 少なくとも 1つを含むことができる。
前記設計変数としてカーカスライン、 折り返しプライライン、 タイヤ外面形状 を表すライン、 タイヤクラウン形状を表すライン、 及び補強材ラインの少なくと も 1つのラインの形状を表す関数から複数のラインを定め、 ステップ (b ) で設 計変数を変化させたときに複数のラインが交差しないように、 ラインの形状を表 す関数の写像関数及び複数のラインが交差しないための制約条件を求め、 該写像 関数により写像された空間においてステップ (b ) を実行することにより制約条 件を満たしながら目的関数の最適値を与える設計変数の値を求めることができる
本発明のステップ (a ) では、 内部構造を含むタイヤ断面形状を表すタイヤ基 本モデル、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断面形状またはタイ ャ構造を決定する設計変数、 及びタイヤ断面形状またはタイヤ構造を制約する制 約条件を定める。 このタイヤ基本モデルには、 タイヤ外面形状を表すラインの他 、 タイヤクラウン形状を表すライン、 タイヤ内部のベルトを表すベルトライン、 タイヤの力一カスを表す力一カスライン、 タイヤ内部のカーカスプライの折り返 しラインを表す折り返しプライライン、 各種補強材のラインを表す補強材ライン 、 タイヤゴム部材のゲージ分布及びベルト部の構造を表す各ベルト層の角度、 幅 、 コード種類、 打ち込み密度、 並びにパターンの形状を表す、 ブロック形状、 ブ ロック溝壁角度、 サイプの位置、 本数、 長さを含ませることができる。 また、 タ ィャ基本モデルは、 複数の要素に分割する有限要素法と呼ばれる手法を用いても 良く解析的手法を用いても良い。 タィャ性能評価用物理量を表す目的関数として は、 操縦安定性を向上させるための空気充塡時のタイヤ周方向ベルト張力や横ば ね定数等のタイャ性能の優劣を支配する物理量を使用することができる。 タイヤ 断面形状を決定する設計変数としては、 カーカスライン、 折り返しプライライン 、 タイヤ外面形状を表すライン、 タイヤクラウン形状を表すライン、 各種補強材 のラインを表す補強材ラインの少なくとも 1つのラインを表す関数等を用いるこ とができ、 タイヤ構造を決定する設計変数としては、 ビー ドフイラ一のゲージ分 布、 ゴムチェ一ファーのゲージ分布、 サイ ドゴムのゲージ分布、 トレッ ドゴムの ゲージ分布、 トレッ ドベースゴムのゲージ分布、 内面補強ゴムのゲージ分布、 ベ ルト間ゴムのゲージ分布、 及びベルトェンドゴムのゲージ分布の少なくとも 1つ のタイヤゴム部材のゲージ分布を表す変数、 ベルト層の角度、 ベルト層の幅、 プ ライの高さ、 プライの折返し量、 ビード部補強材の角度、 幅、 位置、 材質等のベ ルト部及びビード、 サイ ド部の構造を表す変数を用いることができる。 また、 パ ターンを決定する設計変数としては、 ブロックの形状、 ブロック溝壁角度、 サイ プの位置、 本数、 長さ等のパターンの形状を表す変数を用いることができる。 タ ィャ断面形状やタイヤ構造を制約する制約条件としては、 カーカスラインのペリ フェリ値の制約、 上下一次固有振動数の制約、 ベルト層の角度の制約、 ベルト層 の幅、 タイヤ寸度、 パネ定数、 タイヤ変形量、 タイヤ重量、 応力、 歪、 歪ェネル ギー、 転がり抵抗の制約等がある。 なお、 目的関数、 設計変数及び制約条件は、 上記の例に限られるものではなく、 タイヤ設計目的に応じて種々のものを定める ことができる。
次のステップ (b ) では、 制約条件を考慮しながら目的関数の最適値を与える 設計変数の値を求める。 この場合には、 設計変数の単位変化量に対する目的関数 の変化量の割合である目的関数の感度及び設計変数の単位変化量に対する制約条 件の変化量の割合である制約条件の感度に基づいて制約条件を考慮しながら目的 関数の最適値を与える設計変数の変化量を予測すると共に、 設計変数を予測量に 相当する量変化させたときの目的関数の値及び設計変数を予測量に相当する量変 化させたときの制約条件の値を演算し、 予測値と演算値とに基づいて制約条件を 考慮しながら目的関数の最適値を与える設計変数の値を求めると効果的である。 これによつて、 制約条件を考慮し目的関数の値が最適になるときの設計変数の値 が求められる。
また、 前記ステップ (a ) では、 内部構造を含むタイヤ断面形状を表す複数個 のタイャ基本モデルからなる選択対象集団を定め、 該選択対象集団の各タイャ基 本モデルについて、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断面形状ま たはタイヤ構造を決定する設計変数、 タイヤ断面形状、 タイヤ構造、 性能評価用 物理量及びタイヤ寸度の少なくとも 1つを制約する制約条件、 及び目的関数及び 制約条件から評価できる適応関数を定め、 前記ステップ (b ) では、 適応関数に 基づいて前記選択対象集団から 2つのタイヤ基本モデルを選択し、 所定の確率で 各タイヤ基本モデルの設計変数を交叉させて新規のタイヤ基本モデルを生成する こと及び少なくとも一方のタイヤ基本モデルの設計変数の一部を変更させて新規 のタイャ基本モデルを生成することの少なくとも一方を行い、 設計変数を変化さ せたタイャ基本モデルの目的関数、 制約条件及び適応関数を求めて該タイャ基本 モデル及び設計変数を変化させなかったタイヤ基本モデルを保存しかつ保存した 夕ィャ基本モデルが所定数になるまで繰り返し、 保存した所定数のタィャ基本モ デルからな'る新規集団が所定の収束条件を満たすか否かを判断し、 収束条件を満 たさないときには該新規集団を前記選択対象集団として該選択対象集団が所定の 収束条件を満たすまで繰り返すと共に、 該所定の収束条件を満たしたときに保存 した所定数のタイヤ基本モデルのなかで制約条件を考慮しながら目的関数の最適 値を与える設計変数の値を求めることも効果的である。 この場合、 ステップ (b ) において、 設計変数を変化させたタイヤ基本モデルについて、 設計変数の単位 変化量に対する目的関数の変化量の割合である目的関数の感度及び設計変数の単 位変化量に対する制約条件の変化量の割合である制約条件の感度に基づいて制約 条件を考慮しながら目的関数の最適値を与える設計変数の変化量を予測すると共 に、 設計変数を予測量に相当する量変化させたときの目的関数の値及び設計変数 を予測量に相当する量変化させたときの制約条件の値を演算し、 目的関数の値及 び制約条件の値から適応関数を求めて該タイヤ基本モデル及び設計変数を変化さ せなかつたタイャ基本モデルを保存しかつ保存したタィャ基本モデルが所定数に なるまで繰り返すことが更に効果的である。 これによつても、 制約条件を考慮し 目的関数の値が最適になるときの設計変数の値が求められる。 なお、 目的関数及 び制約条件から評価できる適応関数は、 目的関数及び制約条件からタイヤモデル に対する適応度を求める関数を使用することができる。 また、 目的関数、 設計変 数、 制約条件及び適応関数は、 上記の例に限られるものではなく、 タイヤ設計目 的に応じて種々のものを定めることができる。 さらに、 前記のタイヤ基本モデル の設計変数の交叉には、 選択した 2つのタイヤモデルの設計変数についてその一 部または所定部位以降の設計変数を交換する方法がある。 さらにまた、 タイヤモ デルの設計変数の一部の変更には、 予め定めた確率等で定まる位置の設計変数を 変更 (突然変異) する方法がある。 設計変数として力一カスライン、 折り返しプライライン、 タイヤ外面形状を表 すライン、 タイヤクラウン形状を表すライン、 及び補強材ラインの少なくとも 1 つのラインの形状を表す関数から複数のラインを定めたときには、 設計変数を変 化させたときに複数のラインが交差しないように、 ラインの形状を表す関数の写 像関数及び複数のラインが交差しないための制約条件を求め、 該写像関数により 写像された空間においてステップ (b ) を実行することにより制約条件を考慮し ながら目的関数の最適値を与える設計変数の値を求めれば、 複数のラィンが交差 しなくなる。
そしてステップ (c ) では、 目的関数の最適値を与える設計変数に基づいてタ ィャ基本モデル等を変更することによりタイヤを設計する。 また、 写像関数を用 いたときには、 元の関数に写像した設計変数に基づいてタイヤを設計する。
前記のようにステップ (a ) において、 内部構造を含むタイヤ断面形状を表す 複数個のタイヤ基本モデルからなる選択対象集団を定めると共に、 該選択対象集 団の各タイヤ基本モデルについて、 タイヤ性能評価用物理量を表す目的関数、 タ ィャ断面形状またはタイヤ構造を決定する設計変数、 タイヤ断面形状、 タイヤ構 造、 性能評価用物理量及びタイヤ寸度の少なくとも 1つを制約する制約条件、 及 び目的関数及び制約条件から評価できる適応関数を定めたときには、 ステップ ( b ) は次のステップ①乃至ステップ④から構成できる。 ①適応関数に基づいて選 択対象集団から 2つのタイヤ基本モデルを選択するステップ、 ②所定の確率で各 タィャモデルの設計変数を交叉させて新規のタイャ基本モデルを生成すること及 び少なくとも一方のタイヤ基本モデルの設計変数の一部を変更させて新規のタイ ャ基本モデルを生成することの少なくとも一方を行うステップ、 ③交叉や変更に より設計変数を変化させたタイャ基本モデルの目的関数、 制約条件及び適応関数 を求めて該タイヤ基本モデル及び設計変数を変化させなかったタイヤ基本モデル を保存しかつ保存したタイヤ基本モデルが所定数になるまでステップ①乃至ステ ップ③を繰り返すステップ、 ④保存した所定数のタィャ基本モデルからなる新規 集団が所定の収束条件を満たすか否かを判断し、 収束条件を満たさないときには 該新規集団を選択対象集団として該選択対象集団が所定の収束条件を満たすまで 前記ステップ①乃至ステップ④を繰り返すと共に、 該所定の収束条件を満たした ときに保存した所定数の夕ィャ基本モデルのなかで制約条件を考慮しながら目的 関数の最適値を与える設計変数の値を求めるステップ。
このステップ③では、 設計変数を変化させたタイヤ基本モデルについて、 設計 変数の単位変化量に対する目的関数の変化量の割合である目的関数の感度及び設 計変数の単位変化量に対する制約条件の変化量の割合である制約条件の感度に基 づいて制約条件を考慮しながら目的関数の最適値を与える設計変数の変化量を予 測すると共に、 設計変数を予測量に相当する量変化させたときの目的関数の値及 び設計変数を予測量に相当する量変化させたときの制約条件の値を演算し、 目的 関数の値及び制約条件の値から適応関数を演算し、 該設計変数が変化された新規 な夕ィャ基本モデル及び設計変数を変化させなかつたタイャ基本モデルを保存し かつ保存したタイャ基本モデルが所定数になるまで前記ステツプ①乃至ステツプ ③を繰り返して実行することが効果的である。 これによつても、 制約条件を考慮 し目的関数の値が最適になるときの設計変数の値が求められる。
本発明の設計法に基づき設計 ·開発した場合従来の試行錯誤を基本とした設計 •開発と異なり、 コンピューター計算を主体にしてべストモ一ドの設計から設計 されたタイャの性能評価までがある程度可能となり、 著しい効率化を達成でき、 開発にかかる費用が削減可能となるものである。
以上説明したように本発明によれば、 制約条件を満たす目的関数の最適値を与 える設計変数を求め、 この設計変数からタイヤを設計しているので、 ベストモー ドのタイヤ設計を行うことができる、 という効果が得られる。 図面の簡単な説明
第 1図は、 本発明の実施例に使用されるパーソナルコンピュータの概略図、 第 2図は、 本発明の第 1実施例の処理ルーチンを示す流れ図、
第 3図は、 第 1図の設計変数を決定するルーチンを示す流れ図、 第 4図は、 タイヤ基本モデルを示す線図、
第 5図は、 第 4図のタイヤ基本モデルに d 毎に基準点 Pを通る仮想線を引い た状態を示す線図、
第 6図は、 選択された節点、 この節点についての距離 r i 及び見込み角 , 等 の関係を示す線図、
第 7図は、 本発明の第 1実施例の処理ルーチンを示す流れ図、
第 8図は、 第 7図の設計変数を決定するルーチンを示す流れ図、
第 9図は、 タイヤ基本モデルに始点 S、 終点 E、 変曲点 (:、 及び最大幅位置点 Mを示した状態を示す線図、
第 1 0図は、 始点 S、 終点 E、 変曲点 C、 及び最大幅位置点 Mの各点間に想定 した円弧の曲率半径を示す線図、
第 1 1図は、 タイヤ基本モデルにおける元の節点と最小 2乗法により求めた最 適点とのずれ量を示す線図、
第 1 2図は、 r z座標を示す線図、
第 1 3図は、 ξ V 座標を示す線図、
第 1 4図は、 2 7? 2 座標を示す線図、
第 1 5図は、 タイヤのビードフイラ一及びその周辺部材の形状を示す線図、 第 1 6図は、 ベルト層のラインを示す線図、
第 1 7図は、 ベルト層のラインを曲線で近似したことを示す線図、
第 1 8図は、 節点の移動による要素形状悪化を示す線図、
第 1 9図は、 要素分割処理ルーチンを示す流れ図、
第 2 0図は、 要素分割の過程を示す線図であり、 (a ) は領域定義時のモデル を示す線図、 (b ) は各領域内で発生させた節点を示す線図、 (c ) は三角要素 が発生されたことを示す線図、 (d ) は節点のスムージングを行ったことを示す 線図、
第 2 1図は、 ラグランジ 補間処理ルーチンを示す流れ図、
第 2 2図は、 クラウン部を曲線で表すことを示す線図、 第 2 3図は、 クラウン部をラグランジェ補間の次数個で分割したことを示す線 図、
第 2 4図は、 クラウン部のラインに対する法線を示す線図、
第 2 5図は、 クラウン部形状を円弧で表したことを示す線図、
第 2 6図は、 クラウン部形状を独立した中心座標値による円弧で表したことを 示す線図、
第 2 7図は、 パターン表面をラグランジェ次数に従い格子状に分割したことを 示す線図、
第 2 8図は、 パターン表面を細分化したことを示す線図、
第 2 9図は、 本発明の第 1 0実施例の処理ルーチンを示す流れ図、
第 3 0図は、 交叉処理ルーチンを示す流れ図、
第 3 1図は、 山型写像関数を示す線図であり、 (a ) は連続的な山型写像関数 を示す線図、 (b ) は線型的な山型写像関数を示す線図、
第 3 2図は、 谷型写像関数を示す線図であり、 (a ) は連続的な谷型写像関数 を示す線図、 (b ) は線型的な谷型写像関数を示す線図、
第 3 3図は、 突然変異処理ルーチンを示す流れ図、
第 3 4図は、 第 1 1実施例のタイヤ基本モデルを示す線図であり、 (a ) は作 成時のタイヤ基本モデルを示す線図、 (b ) はベルトの存在を説明するためのィ メージ図、
第 3 5図は、 交叉処理ルーチンを示す流れ図、
第 3 6図は、 突然変異処理ルーチンを示す流れ図、
第 3 7図は、 本発明の第 1 2実施例の処理ルーチンを示す流れ図、
第 3 8図は、 ブロックのモデル化を示す線図であり、 (a ) はモデル化するた めのブロック形状を示す線図、 (b ) は有限要素モデルのモデル化例を示す線図 第 3 9図は、 本発明の第 1 3実施例の処理ルーチンを示す流れ図、
第 4 0図は、 本発明の第 1 4実施例の処理ルーチンを示す流れ図、 第 4 1図は、 ヒールアンドトゥ磨耗を説明するための説明図である。 発明を実施するための最良の形態
第 1図には本発明の空気入りタイヤの設計方法を実施するためのパーソナルコ ンピュー夕の概略が示されている。
このパーソナルコンピュータは、 データ等を入力するためのキーボード 1 0、 予め記憶されたプログラムに従って制約条件を満たしかつ目的関数を最適、 例え ば最大または最小にする設計変数を演算するコンピュータ本体 1 2、 及びコンビ ユータ本体 1 2の演算結果等を表示する C R T 1 4から構成されている。
次に、 操縦安定性を向上させるために、 空気充填時のタイヤ周方向ベルト張力 を最適値である最大値にする力一カスラインの形状を設計する第 1実施例につい て説明する。
第 2図は、 第 1実施例のプログラムの処理ルーチンを示すものである。 ステツ プ 1 0 0では、 自然平衡状態のタイヤ断面形状を基準形状とし、 この基準形状を 有限要素法等のように空気充塡時のタイヤ周方向ベルト張力を数値的 ·解析的に 求めることができる手法によりモデル化し、 内部構造を含むタイャ断面形状を表 すと共にメ ッシュ分割によつて複数の要素に分割されたタイャ基本モデルを求め る。 なお、 基準形状は、 自然平衡状態のタイヤ断面形状に限らず任意の形状でよ い。 ここで、 モデル化とは、 タイヤ形状、 構造、 材料、 パターンを、 数値的 ·解 析的手法に基づいて作成されたコンピュータプログラムへのィンプッ トデータ形 式に数値化することをいう。 第 4図はこのタイヤ基本モデルを示すもので、 C L はカーカスライン、 O Lはタイヤ外面形状を表すライン、 P Lは折り返しプライ ライン、 B l, B 2はベルトを表すラインを各々示している。 また、 このタイヤ 基本モデルは、 力一カスライン C Lの複数の法線 N Lし N L 2, N L 3. · · ·によ つて複数の要素に分割されている。 なお、 上記では、 タイヤ基本モデルをカー力 スラインの複数の法線によって複数の要素に分割した例について説明したが、 夕 ィャ外面形状を表すラインの複数の法線や折り返しプライラインの複数の法線に よつて複数の要素に分割してもよく、 また設計目的によつて 3角形等の任意の形 状に分割してもよい。
次のステップ 1 0 2では、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断 面形状を制約する制約条件及びタィャ断面形状を決定する設計変数を決定する。 本実施例では、 操縦安定性を向上させるために、 空気充塡時のタイヤ周方向ベル ト張力を最大にするカーカスラインの形状を設計するため、 目的関数 O B J及び 制約条件 Gを次のように定めている。
目的関数 O B J :ベルト張力のタイヤ周方向成分のタイヤ幅方向の総和 制約条件 G : カーカスラインのペリフヱリ値がタイヤ基本モデルの力一 カスラインのペリフェリ値の土 5 %以内
なお、 上記力一カスラインのペリフヱリ値は、 タイヤ形状を変化させる範囲内 の力一カスラインの節点 (力一カスラインと法線との交点) 間の距離の総和とし て計算することができる。
また、 設計変数であるカーカスラインの形状は、 曲線を近似する第 3図のラグ ランジェ補間ルーチンによって決定される。 このラグランジヱ補間ルーチンのス テツプ 1 3 0では、 第 4図に示すようにタイヤ内部に基準点 Pを設定する。 次の ステップ 1 3 2では、 ベルト端付近の節点 からリムに拘束される付近の節点 q 2 までの範囲をタイヤ形状を変化させる範囲として指定する。 ステップ 1 3 4 では、 節点 q , と基準点 Pとを結ぶ直線を基準線として、 この基準線と、 節点 q 2 と基準点 Pとを結ぶ直線との成す角である見込み角 6を演算し、 ステップ 1 3 6において以下の式に従って角度増分 d 0を演算する。
d 0 = 0 /ラグランジヱ補間の次数 · · · ( 1 )
なお、 ラグランジェ補間の次数は使用者が予めインプッ トしておく。
ステップ 1 3 8では、 第 5図に示すように、 基準線を基準として角度増分 d Θ 毎に仮想線 L , . L 2. L 3, · · ·を想定し、 仮想線に最も近い法線 n 1 , . n 1 2 . n 1 3 , · · ·を選択する。 次のステップ 1 4 0では、 第 6図に示すように、 選択さ れた法線 n 1 , . n 1 2 , n 1 3 . · · '上の最内側の節点 Q し Q 2. Q 3, · · · と基準 点 Pとの距離 rし r 2, r 3. · · · (以下一般式で r , と表す。 ただし、 i = l, 2, · · · ラグランジヱ補間の次数一 1 ) 、 節点 Q,. Q2. Q3. · · ·の見込み角 Θ Θ 2. Θ 3. · · - (以下一般式で 0 , と表す。 ただし、 i = 1, 2, · · · ラ グランジエ補間の次数一 1 ) とを演算する。 そして、 ステップ 1 4 2において距 離 を設計変数として設定する。
このようにして目的関数 OB J、 制約条件 G及び設計変数 r , を決定した後、 第 2図のステップ 1 0 4において、 設計変数 r , の初期値 r oにおける目的関数 0 B Jの初期値 0 B J 0及び制約条件 Gの初期値 G 0を演算する。
次のステップ 1 0 6では、 タイヤ基本モデルを変化させるために設計変数 r , を各々 Δ , ずつ変化させる。 ステップ 1 0 8では、 設計変数に対応する最内側 の節点 Q, 以外の最内側の節点と基準点 Ρとの距離 rm を次式によって演算する 。 ただし、 em は当該節点と基準点 Pとを結ぶ直線の、 基準線からの見込み角で める。
Figure imgf000017_0001
N, = Π Θ m — Θ Π θ , -θ , ) (2) j =l j =l
また、 ごのステップ 1 0 8では、 法線上の最内側の節点 (カーカスライン上の 節点) 以外の節点、 すなわち、 タイヤ外面形状を表すライン OL、 折り返しブラ イライン P L及びベルトを表すライン B 1, B 2上の節点と最内側の節点 Q i と の距離を求め、 最内側の節点 Q, の座標に求めた距離を加算することで法線上の 最内側の節点以外の節点座標を求め、 設計変数を Δ ι· , 変化させた後のタイヤ断 面形状、 すなわちタイヤ修正モデルを決定する。
ステップ 1 1 0では、 ステップ 1 0 8で求めたタイヤ修正モデルについて設計 変数を Δ ι· , 変化させた後の目的関数の値 OB J , 、 制約条件の値 G, を演算し 、 ステップ 1 1 2で以下の式に従って、 設計変数の単位変化量に対する目的関数 の変化量の割合である目的関数の感度 d O B J/d r , 及び設計変数の単位変化 量に対する制約条件の変化量の割合である制約条件の感度 d G/d r i を各設計 変数毎に演算する。
d O B J O B J , -O B J o
d r i Δ r
d G G , — G o
( 3 ) d r i d r
この感度によって、 設計変数を Δ ι· , 変化させたときに目的関数の値及び制約 条件の値がどの程度変化するか予測することができる。 なお、 この感度は、 タイ ャのモデル化に用いた手法や設計変数の性質によっては、 解析的に求められる場 合があるので、 その際にはステップ 1 1 0の演算は不要になる。
次のステップ 1 1 4では、 目的関数の初期値 OB J ο、 制約条件の初期値 G o 、 設計変数の初期値 r 0及び感度を用いて、 数理計画法により制約条件を満たし ながら目的関数を最大にする設計変数の変化量を予測する。 この設計変数の予測 値を用いて、 ステップ 1 1 5でステップ 1 0 8と同様の方法によりタイヤ修正モ デルを決定すると共に、 目的関数値を演算する。 ステップ 1 1 6において、 ステ ップ 1 1 5で演算した目的関数値 OB Jとステップ 1 0 4で演算した目的関数の 初期値 O B J 0との差と、 予めインプッ 卜されたしきい値とを比較することで目 的関数の値が収束したか否かを判断し、 目的関数の値が収束していない場合には ステップ 1 1 4で求められた設計変数値を初期値として、 ステップ 1 0 4からス テツプ 1 1 6を繰り返し実行する。 目的関数の値が収束したと判断されたときに は、 このときの設計変数の値をもつて制約条件を満たしながら目的関数を最大に する設計変数の値とし、 ステップ 1 2 0においてこの設計変数の値を用いてタイ ャの形状を決定する。 次に本発明の第 2実施例を説明する。 本実施例は第 1実施例の設計変数とは異 る設計変数を用いたものである。 第 7図に本実施例のプログラムの処理ルーチン を示すが、 第 2図と対応する部分には同一符号を付して説明を省略する。
ステップ 1 6 0で上記と同様にして目的関数及び制約条件を決定した後、 ステ ップ 1 6 2において設計変数の決定を行う。 第 8図は、 ステップ 1 6 2の詳細を 示すもので、 ステップ 1 5 0において、 第 9図に示すように、 タイヤ基本モデル の形状を変化させる範囲のカーカスライン上の節点から始点 S、 終点 E、 変曲点 C , 及び最大幅位置点 Mを選択し、 これらの点のデータをインプッ 卜する。 次の ステップ 1 5 2では、 第 1 0図に示すように、 始点 S、 終点 E、 変曲点 C、 及び 最大幅位置点 Mの各点間に 1つまたは複数の円弧を想定し、 これらの円弧の曲率 半径 、 R 2 、 R 3 、 R 4 及び見込み角 、 Θ 2 θ , Θ 4 を演算する。 なお、 設計変数を増加させたいときには円弧の数を増加させればよい。
次のステップ 1 5 4では、 円弧と法線との交点の座標を演算し、 第 1 1図に示 すように、 円弧と法線との交点と、 カーカスライン上の節点とのずれ量 ε (—般 式で £ i と表す。 ) を演算する。 ステップ 1 5 6では、 数理計画法により∑ ε , 2 が最小になる円弧の曲率半径 R i (。P。、 見込み角 <。P uを演算し、 この R , ( 。P t >、 S , (。P t >を設計変数 R , 、 θ , とする。 なお、 図示は省略したが目的関数 、 制約条件及び設計変数を決定した後上記と同様にこれらの初期値を演算してお く。
第 7図のステップ 1 6 4では、 力一カスライン上の節点をステップ 1 5 6で求 めた R i , Θ i (。P t >で定まる最適点に移動する。 また、 次のステップ 1 6 6 では、 力一カスライン上の節点以外の節点、 すなわち、 タイヤ外面形状を表すラ イン〇L、 折り返しプライライン P L及びベルトを表すライン B 1, B 2上の節 点とカーカスライン上の節点との距離を求め、 カーカスライン上の節点の座標に 求めた距離を加算することで力一カスライン上の節点以外の節点座標を求め、 力 —カスライン上の節点を最適点に移動させたときのタイヤ断面形状をもって、 夕 ィャ基本モデルとする。 そして、 上記のステップ 1 0 4からステップ 1 0 6と同 様に、 設計変数の初期値 R o、 0 0における目的関数の初期値 O B J 0及び制約 条件の初期値 G 0を演算し、 設計変数を 、 Α θ , 変化させた後のタイヤ修 正モデルにより設計変数毎に感度を演算し、 数理計画法を用いて制約条件を満た しながら目的関数を最大にする設計変数の変化量を予測し、 当該設計変数の予測 値におけるタイャ修正モデルを決定して目的関数 0 B Jを演算し、 初期値 0 B J 0との差としき値とを比較することで目的関数の値が収束したか否かを判断する 。 収束していない場合はステップ 1 1 4で求められた設計変数値を初期値として ステップ 1 0 4からステップ 1 1 6を繰り返し実行する。 目的関数の値が収束し たと判断されたときには上記と同様にしてタイヤの形状を決定する。
なお、 上記実施例では、 設計変数として力一カスラインの形状を採用したが、 折り返しプライラインの形状、 タイヤ外面形状を表すライン、 補強材のライン等 を設計変数として決定してもよい。 次に本発明の第 3実施例を説明する。 本実施例は、 複数のラインの形状決定に 本発明を適用したものであり、 二律背反する性能を最適にしょうとするものであ る。 本実施例では、 車両の乗り心地を損ねることなく操縦安定性を向上させる力 一カスラインの形状、 折り返しプライラインの形状及びタイヤ外面形状を決定す るものである。 この車両の乗り心地を良好にすれば操縦安定性が悪化するため、 車両の乗り心地と操縦安定性とは二律背反する。
本実施例の手順は、 上記実施例と略同様であるが、 目的関数、 制約条件及び設 計変数の選定が異なっている。 本実施例では、 目的関数として操縦安定性を向上 させるための物理量である横ばね定数を用い、 制約条件として乗り心地を支配す る物理量である上下一次固有値が一定という条件を採用しており、 上下一次固有 値が一定という条件の下で横ばね定数が最大になるカーカスラインの形状、 折り 返しプライラインの形状及びタイャ外面形状を決定するようにしている。
設計変数については、 上記実施例の方法をカーカスライン、 折り返しブライラ ィン及びタイヤ外面形状を表すラインの 3種のラインにそのまま適用すると、 各 ラインが最適化の過程で交差するという問題が生じるため、 次のようにして決定 する。
カーカスラインについては第 1実施例で説明したラグランジヱの補間方法を用 いて設計変数 r i を定める。
折り返しプライライン及びタイヤ外面形状を表すラインについては以下のよう にして設計変数を決定する。 まず、 第 1 2図に示すように、 元の r z座標におい てカーカスラインを一辺とする領域 R E , を想定し、 この領域 R E , が r z座標 とは異なる ^ , 77 , 座標で力一カスラインを一辺とする正方形が形成されるよう な写像関数を作成し、 折り返しプライライン及びタイヤ外面形状を表すラインを 写像する。 また、 同様に , V 1 座標において、 折り返しプライラインを一辺と する領域 R E 2 を想定し、 ξ 1 m 座標とは異なる f 2 Ί] 2 座標で折り返しブラ イラインを一辺とする正方形が形成されるような写像関数を作成し、 タイヤ外面 形状を表すラインを写像する。
ここで各写像関数は以下のように表される。
r =∑ N 7? l ) r
∑ N V \ z
77 , =∑ N 77 2 ) 77 そして、 ξ , 7? , 座標における折り返しプライラインを設計変数として決定し 、 2 7? 2 座標におけるタイヤ外面形状を表すラインを設計変数として決定する 。 また、 , 、 ξ 2 〉一 1なる制約条件を加えておく。 これにより、 ライン同士 が交差することがなくなる。
上記の第 1の実施例及び第 3の実施例で得られた夕ィャを実際に試作し試験を 行った結果は以下の表 1のとおりであった。 なお、 上記では 3つのラインを設計変数として採用してが 2つのラインまたは 4つ以上のラインを設計変数として採用してもよい。
1 ]
Figure imgf000022_0001
タイヤサイズ: 2 0 5 6 0 H R 1 5
内圧: 2 . 0 k g / c m 2
試験法:実車走行試験によるフィーリング評価 次に本発明の第 4実施例を説明する。 本実施例は第 3実施例で示した複数のラ ィンの形状決定を、 タイヤのビードフイラ一およびその周辺の部材の形状決定に 適用したものである。 本実施例では、 ビード部分の耐久性を損ねることなく転が り抵抗を軽減するビードフイラ一形状およびゴムチ X—ファ一形状を決定するも のである。 本実施例の手順は第 3実施例と略同様であるが、 目的関数、 制約条件 および設計変数の選定が異なっており、 先ず目的関数および制約条件を次のよう に定めている。
目的関数:転がり抵抗値
制約条件:折り返しプライ端部の荷重時主歪が初期構造の + 3 %以内 また、 設計変数は、 第 1 5図に示すように、 ビードフイラ一 2 0の外形 (ゲー ジ分布) を定義するライン B L a, B L bおよびゴムチェ一ファー 2 2とサイ ド ゴム 2 4との境界ライン K Lである。
なお、 第 3実施例と同様に、 設計変数について、 上記の各ラインをそのまま適 用すると、 各ラインが最適化の過程で交差するという問題が生じるため、 上記実 施例で説明したラグランジェの補間方法や円弧捕間法等を用いて設計変数を定め てもよい。 また、 第 3実施例のように、 所定の座標において想定される領域が異 なる座標系に正方形等で形成されるような写像関数を作成し、 ラインを写像する ようにしてもよい。 次に、 本発明の第 5実施例を説明する。 本実施例は第 3実施例で示した複数の ラインの形状決定を、 ベルト部における各ベルト層間の厚みの決定に適用したも のである。 本実施例では、 ベルト部分の耐久性を損ねることなぐ転がり抵抗を軽 減する各ベルト眉間のゲージ分布を決定するものである。 本実施例の手順は第 3 実施例と略同様であるが、 目的関数、 制約条件および設計関数の選定が異なって おり、 先ず目的関数および制約条件を次のように定めている。
目的関数:転がり抵抗値
制約条件:ベルト端プライ端部の荷重時主歪が初期構造の + 3 %以内 また設計変数は第 1 6図に示すように、 各ベルト層を表わすライン B 1, B 2 , B 3 , B 4であり、 これらのライン B 1〜B 4が変化することにより、 目的関 数を最適にする各ベルト層間のゲージ分布が決定されるものである。 以上の実施例ではラインの形状を表わす方法としてラグランジェ補間法および 円弧補間法を用いたが、 この他 MATHEMATI CAL ELEMENTS FOR COMPUTER GRAPH I CS C Dav i d F. Rogersおよび J. Alan Adams著) に示されている様なスプライン曲線、 B スプライン曲線、 ベジエ曲線あるいは N U R B S (重み付 Bスプライン) 等を用 いて補間するようにしてもよい。 上記の第 4実施例及び第 5実施例を組み合わせて得られたタイャを実際に試作 し試験を行った結果は、 以下の表 2のとおりであった。 [表 2 ]
Figure imgf000024_0001
タイヤサイズ: 1 6 5 S R 1 3
内圧: 1 . 9 k g / c m 2
荷重: 4 2 5 k g f
転がり抵抗値の試験法: ドラム惰性走行による 5 0 k mZ h時の抵抗値測定 次に、 本発明の第 6実施例を説明する。 本実施例はベルト構造に本発明を適用 したものであり、 ベルト部の重量を増すことなく、 ベルト部に生じる歪集中を最 小化することにより耐久性を向上させるものである。 本実施例の手順は第 2図に 示した第 1実施例の手順と略々同様であるカ^ ステップ 1 0 2の目的関数、 制約 条件および設計関数の決定の部分のみ異なっており、 目的関数および制約条件を 次のように定めている。
目的関数:各ベルト層間に発生する荷重時主歪の最大値
制約条件:ベルト総重量が初期構造での総重量の + 1 %以内
また、 設計変数は各ベルト層の角度、 打込み (例えば、 打込み本数や打込み強 度) 及び幅である。 この内、 角度および打込みについては直接その値を設計変数 とすればよく、 幅については以下のようにして設計変数を決定する。
先ず、 第 4図に示すようなタイヤ基本モデルが作成された時、 各ベルト層を構 成する節点の座標値を用いて、 各ベルトのラインを例えば 3次関数等により、 近 似する。 例えば、 第 1 7図に示すようにライン B 1, B 2を 3次関数で近似する この場合、 ライン B 2は、 r = a I z 3 + b i z 2 + C i z 1 + d i
で近似され、 ライ ン B 1は、
r = a 2 z 3 + b 2 z 2 + c 2 z 1 + d 2
で近似される。 目的関数の最適値を与える設計変数を求めるステップにおいて、 各ベルト層を構成する各節点は、 各ベルト層毎に定められた倍率の値を、 その層 を構成する各節点の基本モデルにおける z座標値に乗じた値を z座標とし、 この 値をその層を近似する関数に代入することで求められる値を r座標値とする点に 移動するものとし、 各ベルト層毎の倍率の値を設計変数とする。
ここで、 本実施例の場合、 有限要素法によるモデル化を行った場合には、 各べ ルト層の幅の変化にともなって各ベルト層を構成する節点が移動した結果、 第 1 8図の丸印 3 2の中に示したように要素形状が著しく悪化し計算を続行できなく なる問題がある。 この問題を解決するためには、 第 2図のステップ 1 1 5におい てタイヤ修正モデルを決定ずる際に、 与えられた節点位置から自動的に要素分割 を行う第 1 9図に示したルーチンを追加すればよい。 第 1 9図のステップ 1 8 0 では、 各ベルト層を含む領域を定義するために、 データが入力される。 このデ一 夕には、 節点の座標値及び各領域を定義する節点接続間隔粗さ (節点コネクティヒ 'ティメツ シュ粗さ) のパラメータがある。 このデータ入力により、 3領域を定義する例では 、 第 2 0図 (a ) に示したように、 領域 A R 1, A R 2, A R 3を定義できる。 次のステップ 1 8 2では、 第 2 0図 (b ) に示したように、 各領域内において節 点を発生させて、 次のステップ 1 8 4において第 2 0図 (c ) に示したように、 発生された節点を含む全ての節点を接続することにより三角要素を発生させる。 次のステップ 1 8 6では、 第 2 0図 (d ) に示したように、 領域内の節点位置を 移動させることにより領域内節点のスムージングを行い、 次のステップ 1 8 8に おいてリナンバリング等の有限要素法における内部処理に関するバンド幅の縮小 処理を行い、 本ルーチンを終了する。
なお、 本実施例ではベルトのラインの形状を表わす方法として 3次関数を用い たが、 第 1〜第 2実施例で用いたようなラグランジヱ補間や円弧補間、 あるいは スプライン曲線、 Bスプライン曲線、 ベジヱ曲線、 NURB S等を用いてもよい 上記の第 6実施例で得られたタイヤを実際に試作し試験を行った結果は、 以下 の表 3のとおりであった。
3]
Figure imgf000026_0001
タイヤサイズ: TBR 1 0. 0 0 R 2 0
内圧 7. 2 5 k g/c m2
4 0 5 0 k g f
6 0 km/h
ドラム径 : 1. 7 m
この表 3から、 ベルト層間の荷重時主歪が低減したことにより、 耐久性が向上 したことが理解される。 次に本発明の第 7実施例を説明する。 本実施例は本発明をタイヤクラウン部の 形状決定に適用したものであり、 タイヤと地面との接触領域の形状を変えること なく、 接触領域の圧力分布を均一化し、 磨耗性能を向上させるものである。 本実 施例の手順は第 2図に示した第 1実施例の手順と略々同様であるが、 ステップ 1 0 2の目的関数、 制約条件および設計変数の決定の部分のみ異なっており、 先ず 目的関数および制約条件を次のように定めている。
目的関数:接触領域における圧力分布の標準偏差
制約条件: クラウン部中心およびベルト端部におけるタイヤ周方向接触長が 初期形状における接触長の土 5 %以内
次に、 設計変数であるクラウン部形状については第 2 1図のラグランジヱ補間 ルーチンによって決定される。 このラグランジェ補間ルーチンのステップ 1 7 0 では、 クラウン部の形状を変化させる範囲と、 これを近似するための曲線の数を 指定する。 第 2 2図の例では、 クラウン部形状をタイヤ中心の節点 Q 3から溝端 部の節点 Q 4までと、 溝端部の他方の節点 Q 5からベルト端部付近の節点 q 6ま での 2個の曲線で表わすことを示しているが、 もちろん曲線は例えば節点 q 3か ら節点 Q 6までの 1個でもよいし、 または節点 Q 3と節点 q 4の間の節点 Q 5と 節点 q 6の間を更に分割して 3個以上の曲線で表わしてもよい。
次のステップ 1 7 2では、 第 2 3図に示すように、 各々の範囲を予めインプッ 卜されたラグランジュ補間の次数個に等分割し、 その等分割点に最も近い節点を 選択する。 第 2 3図の例では、 タイヤは分割線 L s 1, L s 2 , L s 3 , L s 4 , L s 5によって等分割され、 これらの分割線 L s 1〜L s 5の各々に最も近い 節点 Q s l, q s 2 , q s 3 , q s 4 , Q s 5が選択される。 次のステップ 1 7
4では、 選択した節点位置におけるクラウン部ラインに対する法線を計算し、 こ れを各々の節点が移動するラインとする。 この場合、 第 2 4図に示すように、 曲 線の範囲を表す節点 q 3, q 4 , q 5 , Q 6と、 ステップ 1 7 2で選択された節 点 Q s l, q s 2 , q s 3 , q s 4 , q s 5に関する法線を演算し、 ライン L p l〜L p 8を決定する。 そして、 次のステップ 1 7 6では、 この各々のライン上 での各節点の、 初期位置からの移動量を設計変数として選択し、 本ルーチンを終 了する。
なお、 選択された節点以外の、 クラウン部上の節点の移動量については、 第 1 実施例で示したと同様に、 ラグランジ 捕間により決定する。 次に、 本発明の第 8実施例を説明する。 本実施例は第 7実施例と同様に本発明 をタイヤクラウン部の形状決定に適用したものであるが、 設計変数のみ第 7実施 例とは異なっている。 また、 本実施例の手順は第 7図に示した第 2実施例の手順 と略々同様であるが、 ステップ 1 6 0および 1 6 2の目的関数、 制約条件および 設計変数の決定の部分のみ異なつており、 先ず目的関数および制約条件を次のよ うに定めている。
目的関数:接触領域における圧力分布の標準偏差
制約条件: クラウン部中心およびベルト端部におけるタイヤ周方向接触長が 初期形状における接触長の土 5 %以内
次に、 設計変数であるクラウン部形状については、 予め指定したクラウン部の 範囲を複数個の円弧で近似する。 例えば、 第 2 5図に示すように、 タイヤ中心の 節点 Q 7からベルト端部付近の節点 q 8までの範囲で、 近似可能なクラウン部の 範囲、 この場合、 範囲 CR 1, CR 2, CR 3の 3領域についてクラウン部形状 を 3個の円弧で表わす。 範囲 CR 1については、 半径が R 1で角度 1で近似で き、 範囲 CR 2については、 半径が R 2で角度 02で近似でき、 範囲 CR 3につ いては、 半径が R 3で角度 03で近似できる。
上記の第 2 5図の例では、 タイヤ中心の節点 q 7からベルト端部付近の節点 q 8までの範囲で、 クラウン部形状を 3個の円弧で表わすことを示しているが、 も ちろん円弧の数は 1 ~2個あるいは 4個以上でも構わない。 また、 各円弧は必ず しも滑らかに連続している必要はなく、 第 2 6図に示すように各円弧の中心の座 標値を独立した変数として扱ってもよい。 すなわち、 範囲 CR 1については、 座 標点 01 (r 1, z 1 ) から半径が R 1 ' かつ角度 θ Γ で近似し、 範囲 CR 2 については、 座標点 02 ( r 2, z 2) から半径が R 2, かつ角度 62' で近似 し、 範囲 CR 3については、 座標点 03 (r 3, z 3) から半径が R 3, かつ角 度 03 ' で近似する。 そして各々の円弧の半径 R i と見込角 0 i、 さらに第 2 6 図の例のような場合は各円弧の中心の座標値 r i, z iをこれに加えて設計変数 とする。
なお、 第 7および第 8実施例ではクラウン部の形状を表わす方法としてラグラ ンジヱ補間および円弧補間を用いたが、 スプライン曲線、 Bスプライン曲線、 ベ ジェ曲線、 NURB S等を用いてもよい。 次の本発明の第 9実施例を説明する。 本実施例は本発明をタイヤのパターン表 面の形状決定に適用したものであり、 パターンと地面とが接触した際の圧力分布 を均一化し、 磨耗性能を向上させるものである。 本実施例の手順は第 2図に示し た第 1実施例の手順と略同様であるが、 ステップ 1 0 2の目的関数、 制約条件お よび設計変数の決定の部分のみ異なつており、 先ず目的関数および制約条件を次 のように定めている。
目的関数:接触領域における圧力分布の標準偏差
制約条件:パターンンの総体積が初期の体積の土 5 %以内
次に設計変数であるパターン表面の形状については、 第 2 7図に示すようにパ ターン表面を、 予めインプッ 卜されたラグランジヱ補間次数に従って格子状に分 割する。 こうして得られた各点 (第 2 7図の例では 1 6点) のパターン厚み方向 の座標値 r iを設計変数とする。 また、 目的関数の計算に有限要素法を用いる場 合、 パターン表面をさらに細分割する必要があるが、 例えば第 2 7図に示すバタ ーン表面を細分割した第 2 8図の例の場合、 パターン表面上の設計変数に対応す る節点以外の節点の座標値については次式により演算すればよい。
n
X ∑ M x
i =l
y = ∑ M y
i =l n
r = ∑ M
i =l で
M , ≡M e ra = N。 (Π · N m ( 7? ) k , k ,
N Π S , ) Π ) ,
j =l j =i
j ≠e j ≠e k2 k2
Nm ( 7j ) = Π ( 77 - 7? i ) Π ( V m — n ) )
j=l j=l
j ≠m j ≠ ただし、 n = k i · k 2
k , : x方向のラグランジヱ補間次数 + 1
k 2 : y方向のラグランジヱ補間次数 + 1 これにより第 2図のステップ 1 0 8及びステップ 1 1 5に対応するパターン表 面形状の決定および修正が可能となる。 なお、 本実施例ではパターン表面の形状 を表わす方法としてラグランジヱ補間を用いたが、 この他に前述の MATHEMATI CAL ELEMENTS FOR COMPUTER GRAPH ICS に示されている様なベジヱ曲面、 Bスプライ ン曲面あるいはクーンズ曲面などを用い、 これらを制御するパラメータを直接設 計変数として選択してもよい。 以上の実施例は各々が単独で実施される必要はなく、 例えば第 1実施例におけ るプライラインの決定と第 6実施例におけるベルト構造の決定とを同時に実施す る、 あるいは第 1実施例で決定したプライラインを用いて第 7実施例におけるク ラウン部形状の決定を実施する等、 各々を組み合わせてもよい。 上記の第 7実施例及び第 9実施例を組み合わせて得られたタイャを実際に試作 し、 乗用車の右前輪に装着し、 一般道路で 1 0 0 0 0 k m走行させた後、 タイヤ のトレツ ド中心と トレツ ド側端部の磨耗量比 (肩落ち磨耗) について調査した結 果は、 以下の表 4のとおりであった。 4 ]
Figure imgf000031_0001
タイヤサイズ: 2 0 5 Z 6 0 H R 1 5
内圧: 2 . 0 k g / c m 2
この表 4から、 圧力分布を均一化したことにより、 肩落ち磨耗性が向上したこ とが理解される。 次に、 操縦安定性を向上させるために、 空気充塡時のタイヤ周方向ベルト張力 を最適値である最大値にするカーカスラインの形状を遺伝的にァルゴリズムによ つて設計する第 1 0実施例について説明する。 なお、 本実施例は、 上記実施例と 略同様の構成であるため、 同一部分には同一符号を付して詳細な説明を省略する 第 2 9図は、 第 1 0実施例のプログラムの処理ルーチンを示すものである。 ス テツプ 2 0 0では、 N個のタイヤ断面形状を有限要素法等のように空気充塡時の 夕ィャ周方向ベルト張力を数値的 ·解析的に求めることができる手法によりモデ ル化し、 内部構造を含むタイヤ基本モデルを求める。 ここで、 モデル化とは、 タ ィャ形状、 構造、 材料パターン数値的 '解析的手法に基づいて作成されたコンビ ユー夕プログラムへのインプッ トデータ形式に数値化することを言う。 なお、 N は予め使用者がインプッ 卜する。 本実施例で用いるタイヤ基本モデルは、 第 1実 施例の第 4図に示したものと同一である。 なお、 タイヤ基本モデルの分割は、 第 1実施例と同様に、 タイヤ外面形状を表すラインの複数の法線や折り返しプライ ラインの複数の法線によって複数の要素に分割してもよく、 また設計目的によつ て 3角形等の任意の形状に分割してもよい。 次のステップ 2 0 2では、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断 面形状を制約する制約条件及び N個のタイヤモデルのタイヤ断面形状を決定する 設計変数を決定する。 本実施例では、 操縦安定性を向上させるために、 空気充塡 時のタイヤ周方向ベルト張力を最大にする力一カスラインの形状を設計するため 、 目的関数〇B J及び制約条件 Gを次のように定めている。
目的関数 O B J :ベルト張力のタイヤ周方向成分のタイヤ幅方向の総和 制約条件 G : 力一カスラインのペリフヱリ値が設定した値の土 5 %以内 なお、 上記カーカスラインのペリフヱリ値は、 タイヤ形状を変化させる範囲内 の力一カスラインの節点 (カーカスラインと法線との交点) 間の距離の総和とし て計算することができる。
また、 設計変数であるカーカスラインの形状は、 第 1実施例で述べた、 曲線を 近似する第 3図のラグランジヱ補間ルーチンによって、 N個のタイヤモデル各々 について決定される。 このラグランジヱ補間ルーチンは、 第 1実施例と同一のた め、 説明を省略する。
ラグランジヱ補間ルーチンを N回繰り返すことにより、 目的関数〇B J、 制約 条件 G及び N個のタイヤモデルの各々の設計変数 r ( J = l , 2 , · · ·, Ν ) を決定した後、 第 2 9図のステップ 2 0 4において、 Ν個のタイヤモデルの各 々の設計変数 r i の各々の目的関数 O B J 及び制約条件 を演算する。 次のステップ 2 0 6では、 ステップ 2 0 4で求めた N個のタイヤモデルの各々 の目的関数〇B J j 及び制約条件 G』 を用いて、 N個のタイヤモデルの各々の適 応関数 F を以下の式 (4 ) に従って演算する。 本実施例では、 例えばベルト張 力を最大にするため、 適応関数による値 (適応度) は、 ベルト張力が大きくなる と大きくなる。
Figure imgf000032_0001
または、
F j = 1 / Φ ; または、
F j =- a Φ: + b
Figure imgf000033_0001
但し、 a =
Figure imgf000033_0002
N ·
∑ Φ;
J=l
N c :定数
r:ペナルティ係数
Omi n =m i η (Φ i . Φ2 、 · · · ΦΝ )
Φ: : Ν個のタイヤモデルの J番目のタイヤモデルのペナルティ関数
( J = 1、 2、 3、 · · · N)
なお、 c及びァは使用者が予め入力する。
次のステップ 2 08では、 N個のモデルの中から交叉させるモデルを 2個選択 する。 選択方法としては、 一般に知られている適応度比例戦略を用い、 N個のタ ィャモデルのある個体 が各々選択で選ばれる確率 P , は以下の式で表わされる
F
P , -
N
∑ F
J =l
但し、 F, : N個のタイヤモデルの中のある個体^の適応関数
FJ : N個のタイヤモデルの J番目の適応関数
J= l、 2、 3、 · · · N
上記実施例では、 選択方法として適応度比例戦略を用いたが、 この他、 遺伝的 アルゴリズム (北野宏明 編) に示されている様な、 期待値戦略、 ランク戦略、 エリート保存戦略、 トーナメント選択戦略、 あるいは G E N I TORァルゴリズ ム等を用いてもよい。
次のステップ 2 1 0では、 選択された 2個のタイヤモデルを、 使用者が予め入 力した確率 Tによって交叉させるか否かを決定する。 ここでいう、 交叉とは、 後 述するように、 2個のタイヤモデルの要素の一部を交換することをいう。 否定判 定で交叉させない場合は、 ステップ 2 1 2において現在の 2個のタイヤモデルを そのままの状態でステップ 2 1 6へ進む。 一方、 肯定判定で交叉させる場合には 、 ステップ 2 1 4において後述するように 2個のタイヤモデルを交叉させる。
2個のタイヤモデルの交叉は、 第 3 0図に示す交叉ルーチンによって行われる 。 先ず、 ステップ 2 0 8において選択された 2個のタイヤモデルをタイヤモデル a及びタイヤモデル bとすると共に、 各々のタイヤモデル a, bの設計変数につ いて並びを含む設計変数べク トルで表し、 タイヤモデル aの設計変数べク トルを V r a 二 ( 3 、 1" 2 3 、 ' . '、 3 、 · · ·、 !* n— , a ) 、 タイヤモデ ル bの設計変数べク トルを V r b = ( r , b , r 2 b 、 · · · r , b 、 · · · r n-. b ) とする。 第 3 0図のステップ 2 5 0では、 予め定めた乱数を生成し、 こ の乱数に応じてタイヤモデル a, bの設計変数べク トルに関する交叉場所 iを決 疋一 9 る。
次のステップ 2 5 2では、 交叉すると決定されたタイヤモデル a, bの設計変 数 r , a , r i b に対して、 以下の式に従って距離 dを求める。
d = I r ー r , b I
次のステップ 2 5 4では、 r , a 、 r i b の取り得る範囲の最小値 Bし 及び最 大値 Bu を用いて、 以下の式に従って正規化距離 d' を求める。
d
d, =
Figure imgf000034_0001
ステップ 2 5 6では、 正規化距離 d' の値を適度に分散させるために、 第 3 1 図 (a) , (b) に示すような山型の写像関数 Z (x) (0≤x≤ 1 , 0≤Z ( x) ≤ 0. 5) を用いて、 以下の式に従って関数値 Zabを求める。
Zab=Z (d' ) このようにして、 関数値 Zabを求めた後、 ステップ 2 5 8において新しい設計 変数 ' a 、 r , ' b を次の式に従って求める。 m i n ( - B r , Bu i )
= V Za
0. 5
m i n (I b - BL -Bu i )
= r + Z
0. 5
または、
m i n (I r B Bu i )
+ Z
0. 5
m i n (I r , b — B L Bu i )
Z a b
0. 5
このようにして、 r , , a 、 r , , b を求めた後、 ステップ 2 6 0で新しい設 計変数の並びである設計変数ベク トル V r' a 、 V r' b は以下のように求めら れる。
V r ' 3 = ( r , 3 , r 2 a ヽ · · · !·, ' a 、 r 1 + 1 b 、 · · ·、 r,-, b )
V r ' b = ( r i b s r 2 b 、 · · · r i ' b 、 r 1 + 1 a 、 · · ·、 r n-. a ) なお、 r , の取り得る範囲の最小値 BL 及び最大値 B uは、 使用者が予め入力 しておく。 また、 写像関数 Z (X) は第 3 2図 (a) , (b) に示すような、 谷 型の関数でもよい。 また、 上記例では交叉場所 iは 1 ケ所であるが、 この他に遺 伝的アルゴリズム (北野 宏明 編) に示されているような、 複数点交叉または 一様交叉等を用いてもよい。
このような交叉によって新規な 2個のタイヤモデルを生成した後、 第 2 9図の ステップ 2 1 6では、 使用者が予め入力した確率 Sで、 突然変異させるか否かを 決定する。 この突然変異は、 後述するように、 設計変数の一部を微小に変更する ことをいい、 最適な設計変数となりうる母集団を含む確度を高くするためである 。 ステップ 2 1 6で、 否定判定で突然変異させない場合には、 ステップ 2 2 6で は現在の 2個のタイヤモデルのまま、 次のステップ 2 2 2へ進む。 肯定判定で突 然変異させる場合には、 次のステップ 2 2 0で以下のようにして突然変異処理を 行う。
この突然変異は、 第 3 3図に示す突然変異ルーチンによって行われる。 先ず、 ステップ 2 6 2では乱数を生成し、 乱数によって突然変異の場所 iを決定する。 次のステップ 2 6 4では、 距離 d' を
0≤ά' ≤ 1
の範囲で乱数により決定する。
次のステップ 2 6 6では、 第 3 1図 (a) , (b) に示すような山型の写像関 数 Z (X) (0≤χ≤ 1で、 O Z (X) ≤ 0. 5 ) あるいは第 3 2図 (a) ,
(b) に示すような谷型の写像関数 Z (x) を用いて、 以下の式に従って、 関数 値 Z dを求める。
Z d = Z (d' )
このようにして、 関数値 Z dを求めた後、 ステップ 2 6 8において新しい設計 変数 ' を以下の式に従って求める。
m i n U r . —BL I、 i r , — B u | )
r , = r Z d
0. 5
または、
m n ( B B u i )
r , + Z d
0. 5
のようにして、 設計変数 r , ' を求めた後、 ステップ 2 7 0で求められる、 新しい設計変数の並びである IS計変数べク トル V r ' は以下のようになる。 V r = ( Γ ι 、 r · · · r , r i + \ · · ·、 r n - このようにして、 新たに生成された 2個のダイヤモデルについて、 目的関数の 値と制約条件の値を第 2 9図のステップ 2 2 2で演算する。 次のステップ 2 2 4 では、 得られた目的関数の値と制約条件の値から前記実施例例と同様に式 (4 ) を用いて適応関数を演算する。
次のステップ 2 2 6では、 上記 2個のタイヤモデルを保存する。 次のステップ 2 2 8では、 ステップ 2 2 6で保存したタイヤモデルの数が、 N個に達したか否 かを判断し、 N個に達していない場合は、 N個になるまでステップ 2 0 8からス テツプ 2 2 8を繰り返し実行する。 一方、 タイヤモデルの数が N個に達した場合 には、 ステップ 2 3 0で収束判定をし、 収束していない場合には、 N個のタイヤ モデルをステップ 2 2 6で保存されたタイヤモデルに更新し、 ステップ 2 0 8か らステップ 2 3 0を繰り返し実行する。 一方、 ステップ 2 3 0で収束したと判断 された場合には、 N個のタイヤモデルの中で制約条件を略満たしながら目的関数 の値が最大となるタイヤモデルの設計変数の値をもって制約条件を略満たしなが ら目的関数を最大にする設計変数の値とし、 ステップ 2 3 2においてこの設計変 数の値を用いてダイヤの形状を決定する。
なお、 ステップ 2 3 0の収束判定は以下の条件のいずれかを満足したら収束と みなす。
1 ) 世代数が M個に達した
2 ) 一番目的関数の値が大きい線列の数が全体の q %以上になつた
3 ) 最大の目的関数の値が、 続く p回の世代で更新されない。
なお、 M、 q、 pは使用者が予め入力しておく。
また、 上記実施例を、 第 2、 第 3、 第 4、 第 5、 第 6、 第 7、 第 8及び第 9実 施例の設計変数に適用してもよい。
上記の第 1 0実施例で得られたタイヤを実際に試作し試験を行った結果は、 以 下の表 5のとおりであった。 5 ]
Figure imgf000038_0001
タイヤサイズ: 2 0 5 / 6 0 H R 1 5
内圧: 2 . 0 k g / c m 2
試験法:実車走行試験によるフィ一リング評価
このように、 第 1 0実施例では、 表 1に示した第 1実施例対比計算量が増加す るため、 設計開発に要した時間は若干増加するが、 より良い性能のタイヤ設計を 行うことができる、 という効果がある。 次に、 本発明の第 1 1実施例を説明する。 本実施例はベルト構造に本発明を適 用したものであり、 二律背反する性能を最適にしょうとするものである。 本実施 例では、 車両の乗り心地を損ねることなく操縦安定性を向上させるベルト構造を 決定するものである。
本実施例の手順は、 第 2 9図に示す第 1 0実施例のプログラムの処理ルーチン と略同様であるが、 目的関数、 制約条件及び設計変数の選定、 交叉方法、 突然変 異の方法が異なる。 本実施例では、 目的関数として操縦安定性を向上させるため の物理量である横ばね定数を用い、 制約条件として乗り心地を支配する物理量で ある縦ばね定数が一定という条件を採用しており、 縦ばね定数が一定という条件 の下で横ばね定数が最大になるベル卜構造を決定するようにしている。
また、 第 3 4図 (a ) に示すようなタイヤ基本モデルが作成された時、 ベルト が存在する可能性のある場所 b , 、 b 2 、 · · ·、 b , 、 · · · , b m を設計変 数とし、
b , = 1ならばその場所にベル卜が存在する b , = 0ならばその場所にベルトが存在しない
とする。
例えば、 第 34図 (b) に太線で示したように、 ベルトが存在するときに、 場 所 b , 、 b 2 、 · · ·、 bn の順序で 1 0 1 1 00 1 1 00となり、 場所 bn + 1 、 bn + 2 、 · · ·、 bm の順序で 1 1 1 1 1 1 1 00 1となる。
次に交叉は第 35図に示す処理ルーチンに基づいて実施される。 まず、 N個の タイヤモデルの第 a番目と第 b番目のタイヤモデルを交叉させるとする。 各々の タイヤモデルの設計変数の並びを設計変数べク トル
VBa = (b 1 3 , b2 a · · · b , b , . ι a . · · · bm a )
VBb = (b, b 、 b2 b 、 · · · bi b 、 bi + 1 b 、 · · · bm b ) とする。
ステップ 272では乱数を生成し、 生成された乱数に応じて交叉場所 iを決定 する。 次のステップ 274では以下のように設計変数べク トルである設計変数の 並びを変更し、 新しい設計変数の並びとして設計変数ベク トル VBa ' 、 VBb ' を求める。
VBa ' = (b, a 、 b2 a 、 …ヽ b ,-! a 、 bi b 、 b i +1 b 、 …ヽ bm b ) VBb ' = (b, b 、 b2 b 、 …ヽ b b 、 b, a 、 b 1+ , a 、 …ヽ b ra a ) 次のステップ 27 6では、 求めた設計変数べク トル VBa ' 、 VBb ' に応じ て、 2個の新しいタイヤモデルを生成する。
なお、 本実施例では交叉場所 iは 1ケ所であるが、 この他遺伝的アルゴリズム (北野宏明 編) に示されている様な、 複数点交叉、 一様交叉等を用いてもよい 次に、 突然変異は第 3 6図に示す処理ルーチンに基づいて実施される。 タイヤ 基本モデルの設計変数の並びを設計変数べク トル
VB= (b , 、 b2 、 · · · b , 、 b , 、 b t十, 、 · · · bm )
とする。
ステップ 27 8では乱数を生成し、 乱数に応じて突然変異の場所 iを決定する 。 次のステップ 2 8 0では突然変異の場所 iの設計変数 b , を以下のように変化 させて新しい設計変数の並びとして設計変数べク トル VB' を生成する。
b , = 0ならば b , ' = 1
b , = 1ならば b , ' = 0
VB' = (b , 、 b2 、 · · · b , 、 b i ' 、 b i +1 、 · · · bm ) 次のステップ 2 8 2では、 求めた設計変数べク トル VB' から、 新しいタイヤ モデルを生成する。 なお、 上記実施例では設計変数としてベルト構造を採用したが、 その他の補強 材等を設計変数として決定してもよい。 また、 ビ一ド部へ捕強材等を入れるか否 かに適用可能である。
上記の第 1 1実施例で得られたタイヤを実際に試作し試験を行った結果は、 以 下の表 6のとおりであった。
6]
Figure imgf000040_0001
タイヤサイズ: 2 0 5Z6 0 HR 1 5
内圧: 2. 0 k g/cm2
試験法:実車走行試験によるフィーリング評価
このように、 第 1 1実施例では、 制約条件を満たす目的関数の最適値を与える 設計変数を求め、 この設計変数からタイヤを設計しているので、 ベストモードの タイヤ設計を行うことができる、 という効果がある。 次に、 第 1 2実施例を説明する。 本実施例はブロック化したタイヤのブロック 形状決定に本発明を適用したものであり、 タイャ周方向に並んだ大きさの異なる プロック間のプロック剛性差を少なくすると共に一つのブロックにおいて各方向 のプロック剛性を均一にするパターンにより偏磨耗性能を向上させるものである 。 第 3 7図には、 本実施例のプログラム処理ルーチンを示した。 ステップ 3 0 0 では、 ブロック 4 2の形状を第 3 8図 (a ) に示すようにブロック 4 2の節点座 標、 サイプ 4 4の節点座標、 ブロックの溝の角度、 溝の深さ、 サイプの個数、 サ イブの幅、 サイプの埋込み角度、 サイプの深さによって定義し、 有限要素法、 解 析的手法で解析できるようにモデル化する。 第 3 8図 (b ) に有限要素モデルの モデル化例を示す。
次のステップ 3 0 2では、 ブロックの性能評価用物理量を表す目的関数、 プロ ックの面積、 サイプの長さ、 サイプの長さ Xサイプの深さ、 サイプの個数等を制 約する制約条件及びブロックの節点座標、 サイプの節点座標、 ブロックの各辺の 溝角度、 ブロックの各辺の溝深さ、 サイプの幅、 サイプの埋込み角度、 サイプの 深さに関した設計変数を決定する。
本実施例では、 偏磨耗性能を向上させるために、 一つのブロックに於いて各方 向のプロック剛性が均一にするパターンを設計するために目的関数、 制約条件、 設計変数を次のように定めている。
目的関数 O B J : ·各方向のブロック剛性を均一
制約条件 G : 'サイプの長さの変化量が初期モデルでのサイプ長さの
± 5 %以内
•サイプの節点座標がプロックの節点座標で囲まれた図形内 に在る
•サイプ節点座標とブロックの外部を表わす直線間の距離が 2删以上
設計変数 :サイブの節点座標 (ブロックの節点座標は固定)
次のステップ 3 0 4では、 設計変数 X , の初期値 X。 における目的関数 O B J の初期値 O B J。 及び制約条件 Gの初期値 G。 を演算する。 次のステップ 3 0 6 では、 ブロックの基本モデルを変化させるために設計変数 X , を各々 Δ Χ , ずつ 変化させる。 次のステップ 3 0 8では、 前記ステップ 3 0 6で実施した設計変数 の変化を含んだプロックをモデル化する。
ステップ 3 1 0では、 ステップ 3 0 8で求めたブロック修正モデルについて設 計変数を Δ Χ 変化させた後の目的関数の値〇B J , 、 制約条件の値 G i を演算 し、 ステップ 3 1 2で以下の式に従って、 設計変数の単位変化量に対する目的関 数の変化量の割合である目的関数の感度 d 0 B J d X , 及び設計変数の単位変 化量に対する制約条件の変化量の割合である制約条件の感度 d GZ d X , を各設 計変数毎に演算する。
d 0 B J O B J i - O B J o
d r Δ Χ
d G G G
d r! Δ Χ ι
この感度によって、 設計変数を Δ Χ ί 変化させたときに目的関数の値及び制約 条件の値がどの程度変化するか予測することができる。 なお、 この感度は、 プロ ックのモデル化に用いた手法や設計変数の性質によっては、 解析的に求められる 場合があるので、 その際にはステップ 3 1 2の演算は不要になる。
次にステップ 3 1 4では、 目的関数の初期値 O B J。 、 制約条件の初期 G。 、 設計変数の初期値 X。 及び感度を用いて、 数理計画法により制約条件を満たしな がら目的関数を最大にする設計変数の変化量を予測する。 この設計変数の予測値 を用いて、 ステップ 3 1 5でステップ 3 0 8と同様の方法によりプロック修正モ デルを決定すると共に、 目的関数値を演算する。 ステップ 3 1 6において、 ステ ップ 3 1 5で演算した目的関数値 O B Jとステップ 3 0 4で演算した目的関数の 初期値 O B J。 との差と、 予めインプッ 卜されたしきい値とを比較することで目 的関数の値は収束したか否かを判断し、 目的関数の値が収束していない場合には ステップ 3 1 4で求められた設計変数値を初期値として、 ステップ 3 0 4からス テツプ 3 1 6を繰り返し実行する。 目的関数の値が収束したと判断されたときに は、 このときの設計変数の値をもつて制約条件を満たしながら目的関数を最大に する設計変数の値とし、 ステップ 3 2 0においてこの設計変数の値を用いてプロ ックの形状を決定する。
上記の第 1 2実施例で得られたタイャを実際に試作し磨耗試験を行った結果は 、 以下の表 7のとおりであった。
7]
Figure imgf000043_0001
タイヤサイズ: 2 0 5 / 6 0 HR 1 5
内圧: 2. O k g/ cm2
試験法: テストコースで 5 0 0 0 km走行
なお、 プロック剛性の均一化度 (B K) は、 ブロックを 3 0度毎にプロック剛 性 K i ( i = 1〜6) を実測し、 以下の式により得られる。
6
∑ (K i一 TO
i=l
B K =
K
6
∑ K (平均ブロック剛性)
i-1
H&T高さ (ヒール 'アンド ' ト一高さ) とは、 タイヤの周方向に繋がるプロ ック力 周方向に鋸刃状に磨耗する現象である。 この H&T高さは、 第 4 1図に 示した。 次に、 第 1 3実施例を説明する。 本実施例は、 第 1実施例と第 1 0実施例を組 合わせたものである。 なお、 本実施例は、 上記実施例と略同様の構成であるため 、 同一部分には同一符号を付して詳細な説明を省略する。
第 1 0実施例では、 ステップ 1 1 8において交叉、 突然変異によって得られた 設計案をもとに目的関数、 制約条件を演算していた。 この場合には、 Gol dberg, D . E. , "Gene t i c Algor i thms i n Search, Opt imi zat i on and Machi ne ィ Learning", Add i son-Wesl ey(1989)に記載されているように局所的な最適解に落ち込まないも のの、 真の最適解を求めることが難しいという問題点がある。 そこで、 第 1 0実 施例のステップ 2 2 2の演算処理として、 第 1実施例のステップ 1 0 4〜 1 1 6 の処理を用いて、 第 1実施例と第 1 0実施例の方法を組み合わせれば、 上記問題 点を解決できる。
第 3 9図には、 本実施例のプログラム処理ルーチンを示した。 ステップ 2 0 0 〜2 2 0は第 1 0実施例と同一であるため説明を省略する。
ステップ 3 4 0では、 上記のようにして得られた 2個の設計案を初期設計案と して目的関数及び制約条件の演算をする。 次のステップ 1 0 6〜 1 1 6は上記第 1実施例の処理と同様にして、 目的関数の値が収束するまで繰り返し実行する。 目的関数の値が収束したと判断されたときには、 次のステップ 2 2 4において得 られた目的関数の値と制約条件の値から適応関数を演算し、 次のステップ 2 2 6 で上記 2個のタイヤモデルを保存する。 このステップ 2 2 6で保存したタイヤモ デルの数が、 N個に達するまでステップ 2 0 8からステップ 2 2 8を繰り返し実 行し、 N個に達した場合には、 ステップ 2 3 0で上記と同様にして収束判定をし 、 収束した場合に、 N個のタイヤモデルの中で制約条件を略満たしながら目的関 数の値が最大となるタイャモデルの設計変数の値をもつて制約条件を略満たしな がら目的関数を最大にする設計変数の値とし、 ステップ 2 3 2においてこの設計 変数の値を用いてダイヤの形状を決定する。
このような、 局所的な最適解に落ち込まず、 真の最適解を得る方法は、 ここで 述べた手法以外に、 第 1 0実施例の方法に前記参考文献に記載されてある焼きな まし法 (S imulated Anneal ing)と呼ばれる方法を組み合わせることもできる。 次に、 第 1 4実施例を説明する。 本実施例は、 最適なタイヤ形状に至る過程を 迅速に処理するものである。 なお、 本実施例は、 上記実施例と略同様の構成であ るため、 同一部分には同一符号を付して詳細な説明を省略する。
第 40図には、 本実施例のプログラム処理ルーチンを示した。 ステップ 40 0 では、 タイヤ形状のデータ 'ベースから予め指定した所定数 Mだけの形状データ を読み込む。 タイヤ形状データベースは、 タイヤの力一カスラインを式 (2) の ラグランジヱ捕間、 第 1 0図に示す円弧を用いる方法、 スプライン曲線、 ベ一ジ ヱ曲線、 NURBS (重き付 Bスプライン) 等によって表現したものである。 次のステップ 4 02では、 タイヤのカーカスライン F (X) を次の式 (5) で 表現する。
M
F (X) =∑ a i f , (X) · · · (5)
i=l
ここで、 a, (i = l〜M) は定数、 f , (x) はタイヤ形状データベースか ら読んだものである。
次に、 第 4図に示すようにカーカスラインに複数の法線 、 NL2 、 NL 3 、 · · ·を引きその線上にタイヤ部材に応じた厚みを考慮して節点を配する。 そして節点を利用して複数の要素の分割する。
次のステップ 4 04では、 タイヤ性能評価物理量を表す目的関数、 タイヤ断面 形状を制約する制約条件及びタイヤ断面形状を決定する設計変数を決定する。 また、 設計変数であるカーカスラインの形状は、 (5) 式の定数 a i によって 決定される。 次のステップ 406ではタイヤ基本モデルを変化させるために設計 変数 a , を各々 Δ a , ずつ変化させる。 次のステップ 40 8では、 設計変数厶 a , だけ変化させた後のタイヤのカーカスラインを (5) 式より求め、 ステップ 4 0 2の手順に従ってタイヤ修正モデルを決定する。 ステップ 1 1 0以降は、 上記 と同様のため、 説明を省略する。
上記の第 1 4実施例でタイヤを設計すると、 実際のタイヤの形状データベース 及び設計者が新た入力した形状を組み合わせて、 形状の最適化を図ることができ るため、 第 1実施例に比較して少ない計算時間で最適解を得ることができる。 敷 かし最初に読み込む形状データベースに最適解が依存するため、 第 1実施例に比 較して効果が若干少なくなる傾向がある。 以下の表 8には、 第 1実施例と同じ目 的関数、 制約条件の基に第 1 4実施例の手順で設計したタイヤの試作、 試験結果 を示した。
8]
Figure imgf000046_0001
タイヤサイズ: 2 0 5 / 6 0 HR 1 5
内圧: 2. 0 k g/c m2
試験法:実車走行試験によるフィ一リング評価

Claims

請 求 の 範 囲 1 . 次の各ステップを含む空気入りタイヤの設計方法。
( a ) 内部構造を含むタイヤ断面形状を表すタイヤ基本モデル、 タイヤ性能評価 用物理量を表す目的関数、 タイヤ断面形状またはタイヤ構造を決定する設計変数 、 及びタイヤ断面形状、 タイヤ構造、 性能評価用物理量及びタイヤ寸度の少なく とも 1つを制約する制約条件を定めるステップ、 (b ) 制約条件を考慮しながら 目的関数の最適値を与える設計変数の値を求めるステップ。
( c ) 目的関数の最適値を与える設計変数に基づいてタイャを設計するステツプ
2 . 前記ステップ (b ) では、 設計変数の単位変化量に対する目的関数の変化量 の割合である目的関数の感度及び設計変数の単位変化量に対する制約条件の変化 量の割合である制約条件の感度に基づいて制約条件を考慮しながら目的関数の最 適値を与える設計変数の変化量を予測すると共に、 設計変数を予測量に相当する 量変化させたときの目的関数の値及び設計変数を予測量に相当する量変化させた ときの制約条件の値を演算し、 予測値と演算値とに基づいて、 制約条件を考慮し ながら目的関数の最適値を与える設計変数の値を求める請求項 1の空気入りタイ ャの設計方法。
3 . 前記ステップ (a ) では、 内部構造を含むタイヤ断面形状を表す複数個のタ ィャ基本モデルからなる選択対象集団を定め、 該選択対象集団の各タイヤ基本モ デルについて、 タイヤ性能評価用物理量を表す目的関数、 タイヤ断面形状または タイヤ構造を決定する設計変数、 タイヤ断面形状、 タイヤ構造、 性能評価用物理 量及びタイヤ寸度の少なくとも 1つを制約する制約条件、 及び目的関数及び制約 条件から評価できる適応関数を定め、
前記ステップ (b ) では、 適応関数に基づいて前記選択対象集団から 2つのタ ィャ基本モデルを選択し、 所定の確率で各タイャ基本モデルの設計変数を交叉さ せて新規のタィャ基本モデルを生成すること及び少なくとも一方のタイャ基本モ デルの設計変数の一部を変更させて新規のタイャ基本モデルを生成することの少 なくとも一方を行い、 設計変数を変化させたタイヤ基本モデルの目的関数、 制約 条件及び適応関数を求めて該タイヤ基本モデル及び設計変数を変化させなかった タィャ基本モデルを保存しかつ保存したタイャ基本モデルが所定数になるまで繰 り返し、 保存した所定数のタィャ基本モデルからなる新規集団が所定の収束条件 を満たすか否かを判断し、 収束条件を満たさないときには該新規集団を前記選択 対象集団として該選択対象集団が所定の収束条件を満たすまで繰り返すと共に、 該所定の収束条件を満たしたときに保存した所定数のタィャ基本モデルのなかで 制約条件を考慮しながら目的関数の最適値を与える設計変数の値を求める請求項 1の空気入りタイヤの設計方法。
4 . 前記ステップ (b ) において、 設計変数を変化させたタイヤ基本モデルにつ いて、 設計変数の単位変化量に対する目的関数の変化量の割合である目的関数の 感度及び設計変数の単位変化量に対する制約条件の変化量の割合である制約条件 の感度に基づいて制約条件を考慮しながら目的関数の最適値を与える設計変数の 変化量を予測すると共に、 設計変数を予測量に相当する量変化させたときの目的 関数の値及び設計変数を予測量に相当する量変化させたときの制約条件の値を演 算し、 目的関数の値及び制約条件の値から適応関数を求めて該タイヤ基本モデル 及び設計変数を変化させなかったタイャ基本モデルを保存しかつ保存したタイヤ 基本モデルが所定数になるまで繰り返す請求項 3に記載の空気入りタイャの設計 方法。
5 . 前記設計変数は、
カーカスライン、 折り返しプライライン、 タイヤ外面形状を表すライン、 タイ ャクラウン形状を表すライン、 及び補強材ラインの少なくとも 1つのラインの形 状を表す関数と、
ビードフイラ一のゲージ分布、 ゴムチヱ一ファーのゲージ分布、 サイ ドゴムの ゲージ分布、 トレッ ドゴムのゲージ分布、 トレッ ドベースゴムのゲージ分布、 内 面補強ゴムのゲージ分布、 ベルト間ゴムのゲージ分布、 及びベルトエンドゴムの ゲージ分布の少なくとも 1つのタイヤゴム部材のゲージ分布を表す変数と、 各ベルト層の角度、 幅、 コード種類、 及び打ち込み密度の少なくとも 1つのべ ルト部の構造を表す変数と、
ブロックの形状及びサイプの位置、 本数、 及び長さの少なくとも 1つのパター ンの形状を表す変数と、
の少なくとも 1つを含む請求項 1の空気入りタイヤの設計方法。
6 . 設計変数として力一カスライン、 折り返しプライライン、 タイヤ外面形状を 表すライン、 タイヤクラウン形状を表すライン、 及び補強材ラインの少なくとも 1つのラインの形状を表す関数から複数のラインを定め、 ステップ (b)で設計変数 を変化させたときに複数のラインが交差しないように、 ラインの形状を表す関数 の写像関数及び複数のラインが交差しないための制約条件を求め、 該写像関数に より写像された空間においてステップ (b)を実行することにより制約条件を満たし ながら目的関数の最適値を与える設計変数の値を求める請求項 1の空気入りタイ ャの設計方法。
7 . 前記目的関数は、 タイヤのベルト張力の周方向成分のタイヤ幅方向の総和量 を表す請求項 1の空気入りタイャの設計方法。
8 . 前記制約条件は、 カーカスラインのペリフヱリ値がタイヤ基本モデルの力一 カスラインのペリフェリ値を基準とした所定範囲内にあることからなる請求項 1 の空気入りタイヤの設計方法。
9 . 前記目的関数は、 操縦安定性を向上させるための物理量である横ばね定数か らなる請求項 1の空気入りタイヤの設計方法。
1 0 . 前記制約条件は、 乗り心地を支配する物理量である上下一次固有値が一定 という条件からなる請求項 1の空気入りタイャの設計方法。
1 1 . 設計変数として複数のラインを定めたときには、 前記 (b ) で設計変数を 変化させたときに複数のラインが交差しないように、 ラインの形状を表す関数の 写像関数及び複数のラインが交差しないための制約条件を求め、 該写像関数によ り写像された空間において前記 (b ) を実行することにより制約条件を満たしな がら目的関数の最適値を与える設計変数の値を求める請求項 1の空気入りタイヤ の設計方法。
PCT/JP1994/000110 1993-01-27 1994-01-27 Design method for a pneumatic tire WO1994016877A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE69419678T DE69419678T2 (de) 1993-01-27 1994-01-27 Entwurfsverfahren für luftreifen
JP51687394A JP4142094B2 (ja) 1993-01-27 1994-01-27 空気入りタイヤの設計方法
EP94905215A EP0642916B1 (en) 1993-01-27 1994-01-27 Design method for a pneumatic tire
US08/302,841 US5710718A (en) 1993-01-27 1994-01-27 Method of designing a pneumatic tire to achieve a best mode under given conditions

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP1166993 1993-01-27
JP5/11669 1993-01-27

Publications (1)

Publication Number Publication Date
WO1994016877A1 true WO1994016877A1 (en) 1994-08-04

Family

ID=11784398

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP1994/000110 WO1994016877A1 (en) 1993-01-27 1994-01-27 Design method for a pneumatic tire

Country Status (6)

Country Link
US (1) US5710718A (ja)
EP (1) EP0642916B1 (ja)
JP (1) JP4142094B2 (ja)
DE (1) DE69419678T2 (ja)
ES (1) ES2134341T3 (ja)
WO (1) WO1994016877A1 (ja)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0647911A3 (en) * 1993-10-06 1997-01-22 Bridgestone Corp Pneumatic tire design process.
US5710718A (en) * 1993-01-27 1998-01-20 Bridgestone Corporation Method of designing a pneumatic tire to achieve a best mode under given conditions
WO1998006550A1 (fr) * 1996-08-08 1998-02-19 Bridgestone Corporation Procede de conception de materiau multicomposant, analyseur d'optimisation et support d'information sur lequel le programme d'analyse d'optimisation de materiau multicomposant est enregistre
WO1998029270A1 (fr) * 1996-12-27 1998-07-09 Bridgestone Corporation Procedes pour conception de pneumatiques
WO1998029269A1 (fr) * 1996-12-27 1998-07-09 Bridgestone Corporation Procedes de conception de pneumatiques
WO1999007543A1 (fr) * 1997-08-08 1999-02-18 Bridgestone Corporation Methode de conception de pneumatique, module analyseur d'optimisation et support de memorisation sur lequel est enregistre le programme d'analyse d'optimisation
WO2000014504A1 (fr) * 1998-09-07 2000-03-16 Bridgestone Corporation Prediction de la performance d'un pneu
JP2001009838A (ja) * 1998-09-21 2001-01-16 Bridgestone Corp 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP3133738B2 (ja) 1998-09-07 2001-02-13 株式会社ブリヂストン タイヤ性能予測方法、流体シミュレーション方法、タイヤ設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、タイヤ性能予測プログラムを記録した記録媒体
JP2001121540A (ja) * 1999-10-25 2001-05-08 Bridgestone Corp タイヤ設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、タイヤ設計プログラムを記録した記録媒体
JP2001124667A (ja) * 1999-10-25 2001-05-11 Bridgestone Corp 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP2002144829A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP2002144817A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP2002144816A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP2004502589A (ja) * 2000-07-11 2004-01-29 ピレリ・プネウマティチ・ソチエタ・ペル・アツィオーニ タイヤの設計方法及びタイヤの製造方法
JP2009269557A (ja) * 2008-05-09 2009-11-19 Toyo Tire & Rubber Co Ltd タイヤの設計方法及びそのプログラム
JP2011225178A (ja) * 2010-04-22 2011-11-10 Bridgestone Corp シミュレーション方法及びシミュレーション装置
JP2017094504A (ja) * 2015-11-18 2017-06-01 住友ゴム工業株式会社 生タイヤモデルの作成方法
JP2019119399A (ja) * 2018-01-10 2019-07-22 住友ゴム工業株式会社 タイヤの設計方法
WO2023167050A1 (ja) * 2022-03-01 2023-09-07 日本製鉄株式会社 コア設計装置、コア設計方法、およびプログラム

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6199026B1 (en) * 1997-11-25 2001-03-06 Sumitomo Rubber Industries, Ltd. Method of and apparatus for simulating rolling tire
US6083268A (en) * 1998-04-27 2000-07-04 Bridgestone/Firestone, Inc. Method for designing pneumatic tires for rolling conditions
JP3238376B2 (ja) * 1998-12-28 2001-12-10 住友ゴム工業株式会社 空気入りタイヤ
US7133736B2 (en) * 2000-07-11 2006-11-07 Pirelli Pneumatici S.P.A. Method for designing a tire and method for producing at least one tire
US7050952B2 (en) * 2000-10-16 2006-05-23 The Yokohama Rubber Co., Ltd. Method of designing rubber composite
US7228732B2 (en) 2001-01-26 2007-06-12 Bridgestone Firestone North American Tire, Llc Tire wear analysis method
JP3498064B2 (ja) * 2001-03-29 2004-02-16 住友ゴム工業株式会社 タイヤ有限要素モデルの作成方法
US7373287B2 (en) * 2004-08-09 2008-05-13 Bridgestone Firestone North American Tire, Llc Finite element analysis tire footprint smoothing algorithm using multiple load cases
US20060118220A1 (en) * 2004-12-06 2006-06-08 The Goodyear Tire & Rubber Company Pneumatic tire with elliptical shoulder
JP2006172344A (ja) * 2004-12-20 2006-06-29 Sumitomo Rubber Ind Ltd 製品形状の設計方法
US20070175565A1 (en) * 2006-01-27 2007-08-02 Brown Jack E Jr Continuous variable pitching methods
JP4057041B1 (ja) * 2006-08-22 2008-03-05 横浜ゴム株式会社 タイヤの設計方法及びタイヤの設計用コンピュータプログラム、並びにタイヤの製造方法
DE102008016625A1 (de) * 2007-04-13 2008-12-04 Toyo Tire & Rubber Co., Ltd. Verfahren zum Entwerfen von Reifen
JP2008276469A (ja) 2007-04-27 2008-11-13 Toyo Tire & Rubber Co Ltd タイヤの設計方法
JP5160146B2 (ja) 2007-06-08 2013-03-13 東洋ゴム工業株式会社 タイヤの設計方法
JP5160147B2 (ja) * 2007-06-08 2013-03-13 東洋ゴム工業株式会社 タイヤの設計方法
WO2009137056A2 (en) * 2008-05-07 2009-11-12 Bridgestone Americas Tire Operations, Llc Method of designing a tire having a target residual aligning torque
CN102803927B (zh) * 2009-06-29 2015-08-19 米其林研究和技术股份有限公司 用于评价轮胎翻新的表面修整的系统和方法
US8900388B2 (en) * 2010-09-30 2014-12-02 Michelin Recherche Et Technique S.A. Tire with tread profile having continuous curvature and method for forming tire
WO2015037440A1 (ja) * 2013-09-10 2015-03-19 住友ゴム工業株式会社 空気入りタイヤ
RU2627890C2 (ru) * 2015-12-30 2017-08-14 Денис Валентинович Никишин Способ конструирования съемного протектора для пневматической шины и съемный протектор для пневматической шины, полученный в соответствии с этим способом
CN109849393B (zh) * 2019-02-28 2021-02-05 安徽佳通乘用子午线轮胎有限公司 一种标准化轮胎的制造方法
CN118596749B (zh) * 2024-08-08 2024-10-22 中策橡胶集团股份有限公司 轮胎性能提升方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0187003U (ja) * 1987-11-30 1989-06-08
JPH03189205A (ja) * 1989-12-19 1991-08-19 Yokohama Rubber Co Ltd:The 乗用車用空気入りラジアルタイヤ

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS558903A (en) * 1978-06-29 1980-01-22 Yokohama Rubber Co Ltd:The Pneumatic tire
US4341249A (en) * 1979-12-03 1982-07-27 The Goodyear Tire & Rubber Company Compression sidewall tire reinforced for deflated operation
US4564737A (en) * 1983-07-11 1986-01-14 F.M.P. Enterprise, Inc. Automatic layout machine for tire tread patterns
US4614562A (en) * 1984-08-13 1986-09-30 The B. F. Goodrich Company Tires building machine with rigid radially expansible drum
JPS649002A (en) * 1985-04-24 1989-01-12 Sumitomo Rubber Ind High performance (low profile radial) tire
US4667722A (en) * 1985-11-07 1987-05-26 The Goodyear Tire & Rubber Company Pneumatic tire
US5151142A (en) * 1986-01-13 1992-09-29 Bridgestone Corporation Heavy duty pneumatic radial tires using rubber reinforcing fiber cords with improved adhesion
JPS6341204A (ja) * 1986-08-06 1988-02-22 Bridgestone Corp 重荷重用空気入りラジアルタイヤ
CA1309934C (en) * 1986-11-05 1992-11-10 Jolan Fredda Lobb Radial ply aircraft tire and rim
US4840213A (en) * 1987-04-03 1989-06-20 Bridgestone Corporation Heavy duty low profile pneumatic radial tires
JPS6487003A (en) * 1987-09-29 1989-03-31 Nippon Steel Corp Method for preventing upward camber in steel rolling
KR940005866B1 (ko) * 1987-10-24 1994-06-24 가부시끼가이샤 브리지스톤 중하중용 공기 래디얼 타이어
US5042546A (en) * 1987-11-16 1991-08-27 The Goodyear Tire & Rubber Company Radial ply pneumatic tire with reverse curvature carcass ply
JPH0820228B2 (ja) * 1990-01-31 1996-03-04 住友ゴム工業株式会社 タイヤ外周面のプロファイル測定方法及びその装置
JPH0627680B2 (ja) * 1990-10-11 1994-04-13 株式会社ブリヂストン タイヤ騒音シミュレート方法及びシミュレータ
US5318085A (en) * 1990-10-22 1994-06-07 General Tire Radial tire having reduced treadwear
US5327952A (en) * 1991-03-08 1994-07-12 The Goodyear Tire & Rubber Company Pneumatic tire having improved wet traction
US5263526A (en) * 1992-09-30 1993-11-23 The Goodyear Tire & Rubber Company Pneumatic tire having specified bead structure
JP2643065B2 (ja) * 1992-10-19 1997-08-20 住友ゴム工業株式会社 空気入りタイヤ
JP3645277B2 (ja) * 1992-12-24 2005-05-11 株式会社ブリヂストン 空気入りタイヤ
US5710718A (en) * 1993-01-27 1998-01-20 Bridgestone Corporation Method of designing a pneumatic tire to achieve a best mode under given conditions
US5392830A (en) * 1993-03-10 1995-02-28 General Tire, Inc. Protective barrier for tire sidewall
JP3686107B2 (ja) * 1993-10-06 2005-08-24 株式会社ブリヂストン 空気入りタイヤの設計方法
US5464051A (en) * 1993-11-05 1995-11-07 The Goodyear Tire & Rubber Company Radial ply tire with specified bead portion design

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0187003U (ja) * 1987-11-30 1989-06-08
JPH03189205A (ja) * 1989-12-19 1991-08-19 Yokohama Rubber Co Ltd:The 乗用車用空気入りラジアルタイヤ

Non-Patent Citations (1)

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

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5710718A (en) * 1993-01-27 1998-01-20 Bridgestone Corporation Method of designing a pneumatic tire to achieve a best mode under given conditions
EP0647911A3 (en) * 1993-10-06 1997-01-22 Bridgestone Corp Pneumatic tire design process.
US6411945B1 (en) 1996-08-08 2002-06-25 Bridgestone Corporation Method and apparatus for designing multi-component material, optimization analyzer and storage medium using learning process
WO1998006550A1 (fr) * 1996-08-08 1998-02-19 Bridgestone Corporation Procede de conception de materiau multicomposant, analyseur d'optimisation et support d'information sur lequel le programme d'analyse d'optimisation de materiau multicomposant est enregistre
WO1998029270A1 (fr) * 1996-12-27 1998-07-09 Bridgestone Corporation Procedes pour conception de pneumatiques
WO1998029269A1 (fr) * 1996-12-27 1998-07-09 Bridgestone Corporation Procedes de conception de pneumatiques
US6531012B2 (en) 1996-12-27 2003-03-11 Bridgestone Corporation Pneumatic tire designing method
WO1999007543A1 (fr) * 1997-08-08 1999-02-18 Bridgestone Corporation Methode de conception de pneumatique, module analyseur d'optimisation et support de memorisation sur lequel est enregistre le programme d'analyse d'optimisation
US7369976B1 (en) 1997-08-08 2008-05-06 Bridgestone Corporation Method of designing tire, optimization analyzer and storage medium on which optimization analysis program is recorded
WO2000014504A1 (fr) * 1998-09-07 2000-03-16 Bridgestone Corporation Prediction de la performance d'un pneu
JP3133738B2 (ja) 1998-09-07 2001-02-13 株式会社ブリヂストン タイヤ性能予測方法、流体シミュレーション方法、タイヤ設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、タイヤ性能予測プログラムを記録した記録媒体
JP2001009838A (ja) * 1998-09-21 2001-01-16 Bridgestone Corp 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP2001121540A (ja) * 1999-10-25 2001-05-08 Bridgestone Corp タイヤ設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、タイヤ設計プログラムを記録した記録媒体
JP2001124667A (ja) * 1999-10-25 2001-05-11 Bridgestone Corp 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP2004502589A (ja) * 2000-07-11 2004-01-29 ピレリ・プネウマティチ・ソチエタ・ペル・アツィオーニ タイヤの設計方法及びタイヤの製造方法
JP4656811B2 (ja) * 2000-07-11 2011-03-23 ピレリ・タイヤ・ソチエタ・ペル・アツィオーニ タイヤの設計方法及びタイヤの製造方法
JP4580087B2 (ja) * 2000-11-06 2010-11-10 住友ゴム工業株式会社 タイヤトレッドプロファイル展開方法
JP2002144829A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP4567172B2 (ja) * 2000-11-06 2010-10-20 住友ゴム工業株式会社 タイヤトレッドプロファイル展開方法
JP2002144817A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP4580086B2 (ja) * 2000-11-06 2010-11-10 住友ゴム工業株式会社 タイヤトレッドプロファイル展開方法
JP2002144816A (ja) * 2000-11-06 2002-05-22 Sumitomo Rubber Ind Ltd タイヤトレッドプロファイル展開方法、及びそれによって決定された空気入りタイヤ
JP2009269557A (ja) * 2008-05-09 2009-11-19 Toyo Tire & Rubber Co Ltd タイヤの設計方法及びそのプログラム
JP2011225178A (ja) * 2010-04-22 2011-11-10 Bridgestone Corp シミュレーション方法及びシミュレーション装置
JP2017094504A (ja) * 2015-11-18 2017-06-01 住友ゴム工業株式会社 生タイヤモデルの作成方法
JP2019119399A (ja) * 2018-01-10 2019-07-22 住友ゴム工業株式会社 タイヤの設計方法
WO2023167050A1 (ja) * 2022-03-01 2023-09-07 日本製鉄株式会社 コア設計装置、コア設計方法、およびプログラム
KR20240124430A (ko) * 2022-03-01 2024-08-16 닛폰세이테츠 가부시키가이샤 코어 설계 장치, 코어 설계 방법 및 프로그램
JP7564497B2 (ja) 2022-03-01 2024-10-09 日本製鉄株式会社 コア設計装置、コア設計方法、およびプログラム
KR102796532B1 (ko) 2022-03-01 2025-04-18 닛폰세이테츠 가부시키가이샤 코어 설계 장치, 코어 설계 방법 및 프로그램
US12299355B2 (en) 2022-03-01 2025-05-13 Nippon Steel Corporation Core design apparatus, core design method, and program

Also Published As

Publication number Publication date
EP0642916A1 (en) 1995-03-15
JP4142094B2 (ja) 2008-08-27
EP0642916A4 (en) 1995-07-05
US5710718A (en) 1998-01-20
EP0642916B1 (en) 1999-07-28
DE69419678T2 (de) 1999-12-02
DE69419678D1 (de) 1999-09-02
ES2134341T3 (es) 1999-10-01

Similar Documents

Publication Publication Date Title
WO1994016877A1 (en) Design method for a pneumatic tire
JP3686107B2 (ja) 空気入りタイヤの設計方法
EP1030170B1 (en) Tire performance predicting method
JP4393589B2 (ja) 空気入りタイヤの設計方法
JP4393595B2 (ja) タイヤの設計方法、最適化解析装置及び最適化解析プログラムを記憶した記憶媒体
JP2001050848A (ja) 空気入りタイヤの設計方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP7205160B2 (ja) タイヤの金型形状設計方法、タイヤの金型形状設計装置、およびプログラム
JP4393590B2 (ja) 空気入りタイヤの設計方法
US8086428B2 (en) Tire design method
JP5630147B2 (ja) タイヤモデル作成方法、及び、それを用いたタイヤ設計方法
JP4234857B2 (ja) 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP2001287516A (ja) タイヤの設計方法、タイヤ用加硫金型の設計方法、タイヤ用加硫金型の製造方法、タイヤの製造方法、タイヤの最適化解析装置及びタイヤの最適化解析プログラムを記録した記憶媒体
JP4401698B2 (ja) タイヤ性能のシミュレーション方法及びタイヤ設計方法
JP3673323B2 (ja) 空気入りタイヤの設計方法
JP2001124667A (ja) 空気入りタイヤの設計方法、タイヤ用加硫金型設計方法、タイヤ用加硫金型製造方法、空気入りタイヤの製造方法、最適化解析装置及びタイヤの最適化解析プログラムを記憶した記憶媒体
JP4318971B2 (ja) タイヤ性能のシミュレーション方法及びタイヤ設計方法
JP4275991B2 (ja) タイヤ性能のシミュレーション方法及びタイヤ設計方法
JP6123159B2 (ja) 構造体断面形状作成方法、構造体断面形状決定方法、構造体製造方法、構造体断面形状決定装置、及びプログラム
JP7328527B2 (ja) タイヤモデル作成方法、タイヤ形状最適化方法、タイヤモデル作成装置、タイヤ形状最適化装置、およびプログラム
JP4755015B2 (ja) タイヤの設計方法
JP2003320824A (ja) タイヤモデルの作成方法
JP7576447B2 (ja) タイヤ設計方法、プログラム及びタイヤ設計装置
JP2020201736A (ja) タイヤの初期形状設計方法、タイヤの初期形状設計装置、およびプログラム
Stabile et al. Multi-objective structural optimization of vehicle wheels
JP2005075296A (ja) タイヤの性能予測方法及びタイヤの設計方法

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): JP US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FR GB GR IE IT LU MC NL PT SE

WWE Wipo information: entry into national phase

Ref document number: 08302841

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 1994905215

Country of ref document: EP

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWP Wipo information: published in national office

Ref document number: 1994905215

Country of ref document: EP

WWG Wipo information: grant in national office

Ref document number: 1994905215

Country of ref document: EP

点击 这是indexloc提供的php浏览器服务,不要输入任何密码和下载