JP3466583B2 - Simulation method for predicting the performance of products made of viscoelastic materials - Google Patents

Description

Detailed Description of the Invention

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a simulation method for predicting the performance of a product made of a viscoelastic material, and more particularly to a simulation method for accurately predicting the performance of a product made of a viscoelastic material by simulation. is there.

[0002]

Viscoelastic materials represented by polymer materials such as rubber and elastomer have been widely applied to various products such as tires, various balls used in sports competitions, and rolls used in printing machines.

With respect to various industrial products such as the above viscoelastic materials and metallic materials, product development using simulation is being carried out in various fields in order to save the cost and time for trial manufacture. For example, in order to predict the resilience performance of a golf ball, a simulation method of an actual impact test has been proposed.

Conventionally, in the case of performing the above-mentioned simulation, the constituent materials of the ball generally measured by a viscoelasticity spectrometer for measuring the rigidity and viscosity of the material, a tensile tester for measuring the longitudinal elastic modulus (Young's modulus), and the like. The physical properties of are used as input values to the simulation.
In particular, the viscoelasticity spectrum meter measures the physical property values of a test piece to which a dynamic strain is applied, and is therefore useful for simulating a product made of a viscoelastic material.

[0005]

However, a large amount of deformation cannot be given to the test piece by the measurement with the above-mentioned viscoelasticity spectrum meter, a tensile tester for measuring the longitudinal elastic modulus, etc. The maximum strain rate applied to the test piece made of is from 0.001 / s to 1.
It is as small as 0 / s and the maximum compressive strain is 0.00.
The value is as small as 01 to 0.02.

On the other hand, a product made of a viscoelastic material may undergo large deformation at high speed due to the influence of external force when it is actually used. For example, in the material forming the golf ball, the maximum strain rate is 500 / s to 5000 / s when actually hit.
s, and the maximum compressive strain is mainly 0.05 to 0.
It is about 50, which is a very large value.

Thus, the viscoelasticity spectrum meter,
In the measurement with a tensile tester that measures the longitudinal elastic modulus, etc., it is not possible to measure the physical property value of the viscoelastic material under the same high speed and large deformation conditions as the above-mentioned actual use conditions. , The maximum strain rate and the maximum compressive strain are very different. Therefore, there is a problem in that it is not possible to perform an accurate simulation in consideration of the physical properties of the viscoelastic material in the conventional simulation in which the physical property values of the material obtained by a viscoelastic spectrum meter, a tensile tester, etc. are input.

That is, it is known that the deformation behavior of a viscoelastic material subjected to an impact load is significantly influenced by the amount of deformation or the deformation speed, unlike the case of receiving a static load. In particular, a polymer material typified by rubber or elastomer causes a very high-speed deformation of 10,000 minutes or a few thousandths of a second due to an impact load, and the amount of deformation also becomes a very large value of several tens of percent. There are many viscoelastic materials with such high-speed and large-deformation behavior, and highly accurate simulations are required for efficient product development. Specifically, in products such as golf balls that are subject to impact during use, the dynamic behavior under high-speed and large-deformation conditions and material properties affect the performance of the product. Simulation is essential for product development.

Further, as the viscoelastic material, there is a material in which physical values such as a loss coefficient change according to the magnitude of strain and strain rate when an external force such as an impact load is applied. That is, since the viscoelastic material has a wide range of deformation speed and size, the physical property values are not linearly linear depending on the deformation speed and size.
There is a property with a strong nonlinear change. Specifically, as the material deforms due to external force and the strain increases, S
-The loop area of the S (strain-stress) curve increases,
The physical properties such as loss factor are the deformation state of the material (deformation speed,
It changes with the size) and shows non-linearity. Many viscoelastic materials have strong non-linearity, and highly accurate simulation is required for products using such viscoelastic materials.

However, there has been no method capable of accurately expressing a physical property showing the non-linearity of a viscoelastic material, for example, a phenomenon in which the loss coefficient strongly changes non-linearly depending on the speed and size of deformation of the material. Conventionally, simulations have been performed on products such as golf balls made of a viscoelastic material assuming that the physical properties of the viscoelastic material hardly change. As a result, there is a problem that the performance of a product made of viscoelastic material in actual use cannot be accurately predicted by simulation, and in reality, the performance of the product must be evaluated by making a prototype. There is a problem that there is no.

The present invention has been made in view of the above problems, and it is an object of the present invention to accurately predict the performance of a product made of a viscoelastic material whose material properties are non-linear under actual use conditions by simulation. I am trying.

[0012]

In order to solve the above-mentioned problems, the present invention provides a strain, a strain rate, and a stress generated in the viscoelastic material under measurement conditions assuming a practical use state of a product made of the viscoelastic material. The time value data of the viscoelastic material is derived from the time history data of the strain, strain rate, and stress, and the viscoelastic model considering the viscosity of the viscoelastic material. , A product made of the viscoelastic material is set as a product model to be analyzed, and the relationship between the strain, strain rate, and viscous resistance is input to the product model, and the change in viscous resistance due to the difference in strain and strain rate is considered. A simulation method for predicting the performance of a product made of a viscoelastic material, characterized by predicting the performance of a product model made of the viscoelastic material To have.

As described above, the viscoelastic resistance is derived using the viscoelastic model considering the viscosity of the viscoelastic material, and the relationship between strain, strain rate, and viscous resistance is input to the product model for simulation, so It is possible to accurately represent a phenomenon in which a physical property value indicating the non-linearity of a material changes non-linearly depending on the speed and size of deformation of the material. In addition, since the strain, strain rate, and stress values used in the simulation are measured under conditions that assume the actual usage state of products made of viscoelastic materials, they can be used for various deformation states of actual viscoelastic materials. The simulation can be performed. Therefore, the relationship between strain and strain rate changes depending on the deformed state of the material, and it is possible to accurately predict the performance by simulation even for products made of viscoelastic materials whose material properties such as loss coefficient show non-linearity. ing.

In this simulation method, the momentary values of strain, strain rate, and stress occurring in the viscoelastic material are measured under measurement conditions that assume the actual use state of the product made of the viscoelastic material. Specifically, the measurement conditions are determined assuming that the viscoelastic material is deformed by the external force applied to the product when the product is actually used. Under the above measurement conditions, the values of strain, strain rate, and stress occurring in the viscoelastic material are measured every moment to obtain time history data. Therefore, it is possible to obtain information on the state of deformation of the viscoelastic material when an external force assuming the actual use state is applied to the viscoelastic material from the time history data of strain, strain rate, and stress. As a result, it is possible to accurately predict the physical properties of a viscoelastic material accompanied by high-speed and large-deformation behavior due to impact load.

It is preferable that the values of the strain, strain rate, and stress are measured under a plurality of different measurement conditions. By changing the magnitude of the external force applied to the product and setting the measurement conditions for multiple patterns, it is possible to obtain data for various patterns regarding strain, strain rate, and stress, and improve the accuracy of the simulation input values. be able to. In addition, as much as possible, in order to obtain data for as many patterns as possible, the strain, strain rate, and stress values are
It is preferable to measure until the strain becomes almost zero after the strain is generated by applying an external force to the viscoelastic material, and the shorter measurement time interval is better from the viewpoint of simulation accuracy.

As the viscoelastic model considering the viscosity of the viscoelastic material, a viscoelastic model including a spring and a dashpot is preferable. Since such a viscoelastic model can simplify the viscosity of the viscoelastic material, it is possible to easily consider the influence of the viscosity on the deformed state of the material.
Specifically, a Maxwell model, a Voigt model, or a combination of a plurality of springs and dashpots may be used. The two-element model is preferable from the viewpoint of model simplification. In such a viscoelastic material model, the viscous resistance η of the dashpot and the rigidity (longitudinal elastic modulus E or shear coefficient) of the spring are made variable. The shear modulus is a physical property value determined by the longitudinal elastic modulus (Young's modulus) E and the Poisson's ratio.

The time history data of the viscoelastic resistance of the viscoelastic material is derived from the time history data of the strain, strain rate, and stress and the viscoelastic model considering the viscosity of the viscoelastic material. Specifically, the above-mentioned viscoelasticity model formulates the respective relationships among the strain generated in the viscoelastic material, the strain rate, the longitudinal elastic modulus of the viscoelastic material, and the viscous resistance. By formulating in this way, the viscous resistance is expressed as a function of strain and strain rate. Then, the longitudinal elastic modulus of the viscoelastic material is derived from the strain and stress generated in the viscoelastic material obtained by the above measurement, and along with the strain and the strain rate, the value of the viscous resistance is derived by substituting into the above function. . Since the time history data regarding strain, strain rate, and stress is obtained in the above measurement, corresponding time history data can be similarly obtained regarding viscous resistance. in this way,
The viscous resistance is determined by the values of strain and strain rate, and the viscous resistance also changes according to the time change.

A product made of the above viscoelastic material is set as a product model to be analyzed, and the product model including the product model
Input data for calculation including constraint conditions (or input
The relationship between the strain, strain rate, and viscous resistance is input to (data), and simulation is performed in consideration of the change in viscous resistance due to the difference in strain and strain rate. The above relationship of strain, strain rate, and viscous resistance is input to the product model during simulation calculation. Specifically, each two-dimensional data of the relationship between strain and strain rate and the relationship between strain and viscous resistance is input. Can be calculated.
In addition, the relationship between strain, strain rate and viscous resistance is made into three-dimensional data, and the data is input as a function of viscous resistance.
It can also be calculated.

For each two-dimensional data of the relationship between strain and strain rate and the relationship between strain and viscous resistance, the strain, strain rate and corresponding viscous resistance are written as input data using the above relationships. . Specifically, strain, strain rate, etc. are measured under a plurality of measurement conditions, and for each pattern under different measurement conditions, the correspondence between strain and strain rate is recorded from the time series data of strain and strain rate. To do. The value of the viscous drag corresponding to each curve is also recorded. By arranging these three correspondences of strain, strain rate, and viscous resistance, the viscous resistance at a certain strain and strain rate is accurately derived and calculated.

Since the physical properties of the viscoelastic material can be more accurately realized as the number of measurements under different strains and strain rates is increased, the strains, strain rates, etc. can be measured under a plurality of different measurement conditions as described above. It is preferable to carry out. However,
Since more measurement data requires more calculation time during simulation, it is preferable to calculate the viscous resistance using interpolation when the strain and strain rate are not the same as the measurement data measured under the specified measurement conditions. . The above-mentioned interpolation is, for example, various methods such as linear interpolation from two values of viscous resistance under the conditions of close strain and strain rate, or interpolation using individual measured values under all measured conditions. Can be done by. By performing such an interpolation operation, the viscous resistance can be calculated corresponding to the changes in strain and strain rate that occur in the material due to the difference in measurement conditions.

According to the present simulation method, even if the viscoelastic material is a material whose physical properties are non-linear,
It is possible to accurately simulate the physical properties of materials and the deformation behavior assuming actual usage. Even in a material in which the physical properties of the viscoelastic material are not linearly linear and exhibit strong nonlinearity accompanied by changes in the deformation rate and the magnitude of deformation, as described above, the viscous resistance is strained and the strain rate is
Non-linearity can be considered by calculating the viscous drag using a viscoelastic model that is determined by two values. In this simulation method, particularly in a material exhibiting a strong nonlinearity with a loss coefficient, it is possible to accurately perform a simulation and predict the performance of a product assuming an actual use state.

The simulation of the present invention is preferably performed by the finite element method as shown below.

When performing the simulation by the finite element method analysis, a large number of nodes and elements are set in the product model. That is, when simulating a product composed of a viscoelastic material by the finite element method and predicting the material properties, the stiffness (longitudinal elastic modulus or shear coefficient) of the spring of the viscoelastic model is determined, and the viscous resistance indicating the viscosity is determined. As described above, by determining the strain and strain rate generated in each element for each element, it is possible to represent the material physical properties under the conditions of appropriate strain and strain rate for each element. It should be noted that instead of inputting the longitudinal elastic modulus, a shearing coefficient may be input from the relationship with the Poisson's ratio. Whether to input the longitudinal elastic coefficient or the shear coefficient can be changed according to the specifications of the finite element method program.

The above strain, strain rate and stress are measured by a split Hopkinson bar tester. Using the Split Hopkinson bar tester, a test piece can be subjected to high-speed and large-deformation strain at a high speed of 10,000 minutes or several thousandths of a second, and the amount of deformation is also tens of percent.
Time series data of strain, strain rate and stress of viscoelastic material under high speed and large deformation condition can be obtained. If the measurement conditions of the split Hopkinson rod tester are strains that occur in the viscoelastic material when an impact load is applied to the product, and conditions equivalent to the strain rate, high-speed large deformation state of the viscoelastic material caused by external force, etc., It is possible to obtain the material properties of the viscoelastic material corresponding to various states. Therefore, the accuracy of the simulation can be improved by using the material properties measured by this split Hopkinson bar tester.

The split Hopkinson rod tester can measure the material properties of the test piece in various strain and strain rate regions simply by changing the impact speed of the striking rod that impacts the test piece. Therefore, it is possible to easily obtain material properties in various strain and strain rate patterns.

The split Hopkinson bar tester was originally used to evaluate the impact behavior of metallic materials.
In the present invention, the split Hopkinson bar tester is modified and used for evaluation of viscous viscoelastic materials. The measuring method and the like of the split Hopkinson bar tester will be described later.

If it is possible to apply a strain of large deformation to the test piece at a high speed and the physical properties of the material can be measured under the conditions of strain and strain rate under the conditions in which the product is actually used, It goes without saying that the physical properties of the material may be determined by a measuring method other than the split Hopkinson bar.

The maximum value of the strain generated in the viscoelastic material at the time of measuring the strain, strain rate and stress under the measurement conditions assuming the actual use state of the product made of the viscoelastic material is 0.05 to 0. And / or the maximum value of the strain rate is 500 / s to 10000 / s,
It is preferably in the range of 500 / s to 5000 / s. The maximum value of the above-mentioned strain, the range of the maximum value of the strain rate, the strain that occurs during high-speed large deformation of the viscoelastic material, since the condition of the strain rate, to predict the performance of the product during high-speed large deformation, It is preferable to use three time-series data of strain, strain rate, and stress under the conditions.

The viscoelastic material is a golf ball material, and the product model is a golf ball. Golf balls are products made of viscoelastic materials that are subject to external forces such as impact loads when subjected to high-speed deformation or large deformation during actual use. The high-speed, large-deformation state has a large effect on the performance of the golf ball itself. Exert. Therefore, the analysis by the simulation method of the present invention is very useful for predicting the performance of the golf ball, and the performance of the golf ball can be predicted with high accuracy without trial production. Can be realized.

Simulation of a collision phenomenon between the golf ball and a hitting object assuming a golf club head,
The behavior of the golf ball at the time of the collision is predicted. In practice, the golf club head can predict the material properties of the viscoelastic material that composes the golf ball under the same strain, strain rate, and stress conditions as when the golf ball is hit. It is possible to predict the behavior of the golf ball such as the coefficient of restitution and the state of deformation of the golf ball when hit.

The above-mentioned golf ball may be a so-called one-piece golf ball consisting of a single body such as a cross-linked rubber layer, or a so-called two-piece golf ball having a cover such as a cross-linked rubber layer covered with a cover. It may be a so-called multi-piece golf ball composed of layers or more, and is applicable to golf balls of any structure made of a viscoelastic material.

The viscoelastic material includes all materials having viscoelasticity. Examples thereof include thermoplastic resins, thermosetting resins, various elastomers, rubbers, and the like, and these can be used alone or as a mixture. Further, these simple substances and mixtures include a colorant, a deterioration inhibitor,
Various additives such as a cross-linking agent may be blended as necessary. It can be applied to any material as long as it can measure strain, strain rate, and stress under actual use conditions of the product.

Specifically, synthetic resins such as ionomer resins used as golf ball constituent materials, polybutadiene, natural rubber, polyisoprene, styrene-butadiene copolymers, ethylene-propylene-diene copolymers (EPDM), etc. Is mentioned.

Examples of the product made of the viscoelastic material include, in addition to golf balls, rubber rollers used in printing machines such as printers, tires, and other sports goods made of the viscoelastic material. In particular, it is suitable for a product that receives an impact load and exhibits a large amount of deformation at high speed during use, and the performance y and dynamic behavior of the product can be accurately predicted by this simulation method.

In the present simulation method, the strain, strain rate, longitudinal elastic modulus or shear coefficient is input to the input data for calculation including the product model, and the longitudinal elastic coefficient or shear due to the difference in strain or strain rate is input. The accuracy of simulation can be further improved by considering the coefficient.

Further, in the present simulation method, the above-mentioned viscous resistance is divided into a "load state" in which the strain applied to the viscoelastic material increases in the compression direction and a "restoration state" in which the amount of compression gradually decreases. Then, the simulation can be performed by changing it, and thereby, even when the strain has the same value, different values are derived in the strain load state and the strain unload state, and the simulation accuracy can be further improved.

[0037]

DETAILED DESCRIPTION OF THE INVENTION Embodiments of the present invention will be described below with reference to the drawings. In the simulation method of the first embodiment of the present invention, a material containing butadiene rubber, which is a constituent material of a golf ball, as a main component is used as a viscoelastic material exhibiting nonlinearity. Using a material containing the butadiene rubber as a main component as a test piece, a split Hopkinson rod tester improved for measuring viscoelastic materials was used under high-speed, large-deformation measurement conditions assuming actual use of a golf ball, The strain, strain rate, and stress generated in the material containing butadiene rubber as the main component are measured every moment. The measuring method using the split Hopkinson bar tester will be described later.

The split Hopkinson rod tester can measure the physical properties of materials in various strain and strain velocity regions by changing the impact velocity of the striking rod. In this embodiment, four patterns of collision speeds (7 m / s, 14 m / s)
s, 20 m / s, 25 m / s) is adopted, and measurement is performed under four measurement conditions, and strain, strain rate, and time history data of stress are obtained for each collision pattern. The time history data of strain ε, strain rate ε ′, and stress σ are shown in FIGS. 1, 2, and 3, respectively.

FIG. 4 shows a strain-stress diagram drawn from the strain time history data and the stress time history data. In the strain-stress diagram of FIG. 4, using the strain ε _max when the strain value is maximum and the value of the stress σ _a at that time, the longitudinal elastic modulus E of the test piece for each collision pattern is calculated by the following mathematical formula 1. Is calculated. (Equation 1) E = [sigma] _a / [epsilon] _max --- (1) In the material containing butadiene rubber as the main component used in the present embodiment, the impact velocity of the striking rod was changed in a split Hopkinson rod testing machine. Also, the longitudinal elastic modulus in each collision pattern shows the same value.

For a product made of a viscoelastic material, a viscoelastic model considering the viscosity of the viscoelastic material is set in order to perform a simulation in which the viscosity is taken into consideration. Specifically, in this embodiment, a basic two-element Voigt model is used as a viscoelastic model including a spring and a dashpot as shown in FIG. That is, in the viscoelastic material model, the dashpot viscous resistance η and the spring rigidity (longitudinal elastic modulus E or shear coefficient) are used in a variable manner.

As shown in FIG. 5, a two-element Voigt (Vo
igt) model, the response generated by the spring
Force σ₁, The stress generated in the dashpot is σ_TwoAnd
The overall stress σ can be expressed by the following Equation 2.
it can. (Formula 2) σ = σ₁+ Σ_Two = Eε + ηε ′ −−− (2)

Therefore, from the above equation (2), the viscous resistance η can be expressed by the following equation (3) in such a viscoelastic model. (Equation 3) η = (σ−Eε) / ε ′ −−− (3)

Therefore, the viscous resistance η corresponding to the strain ε and the strain rate ε ′ is calculated from the time history data of strain ε, strain rate ε ′ and stress σ shown in FIGS. It can be calculated moment by moment, and as shown in FIG. 6, time history data of viscous resistance is obtained at each collision velocity.

As described above, after the longitudinal elastic modulus E of the viscoelastic material is determined from the strain and stress obtained by the measurement under the high-speed large deformation condition, the above mathematical formula derived from the viscoelastic model considering the viscosity is obtained. The viscous resistance η is calculated according to 3, and this viscous resistance η is applied to the simulation.

From the strain time history data and the viscous resistance time history data, the relationship between the strain ε and the viscous resistance η is shown in FIG.
Shown in. Further, the relationship between the strain ε and the strain rate ε ′ is shown in FIG. 8 based on the time history data of strain and the time history data of strain rate. As described above, the viscous resistance η changes depending on the strain ε and the strain rate ε ′, so the viscous resistance can be determined by the strain ε and the strain rate ε ′. It is a function.

The relationship between strain and viscous resistance shown in FIG. 7 and the relationship between strain and strain rate shown in FIG. 8 are written as input data in each collision pattern, as will be described later, and the above is obtained by finite element method analysis. Simulation is performed considering the change in viscous resistance due to the difference in strain and strain rate.

In this embodiment, a golf ball model is set as a product model made of a viscoelastic material, and a simulation is performed assuming a collision of a golf club head (hit object) with a golf ball. As shown in FIG. 9, the golf ball model 10 assumes a one-piece ball containing high-cis polybutadiene rubber as a main component, and has a diameter of 42.8.
It has a spherical shape of mm.

When performing the finite element method analysis, initial conditions are set in the product model. That is, the initial conditions such as the size, shape, and structure of the golf ball model 10 are set, and the golf ball model 10 is divided into a large number of elements 11
And a large number of nodes 12 are obtained. The number of elements in all ranges of the golf ball model is preferably 1000 to 100,000, more preferably 2500 to 20000 in the whole model. The upper limit is
This is because the capacity of the computer at the present stage is taken into consideration. As the capacity of the computer is improved in the future, the upper limit range can be easily imagined because the analysis time is shortened.

Based on the above set conditions, the data of the golf ball model 10 and the relationship between the strain, the strain rate, and the viscous resistance are written in the input data during the simulation calculation. When the calculation is executed, an appropriate viscous resistance is calculated for each element every moment, and the formula is calculated using the viscous resistance. In this simulation, each two-dimensional data of the relation between strain and strain rate and the relation between strain and viscous resistance is input and calculated. From the time series data of strain and strain rate for each of the four patterns with different measurement conditions, record the correspondence between strain and strain rate, and write the strain, strain rate and corresponding viscous resistance as input data. Become. In the case of strains and strain rates that are not the same as the strain and strain rate measurement data measured under the four measurement conditions, linear interpolation is performed using viscous resistance binary values under the conditions of adjacent strain and strain rate. ing.

More specifically, when one element is focused on, information on the strain and strain rate occurring in that element is obtained. Next, from the relationship between the strain and strain rate of each collision case (measurement case) obtained from the measurement, refer to the value of the strain rate when the strain of the same value occurs in each case, and We search for two cases with strain rate values close to the strain rate values. When the strains corresponding to these two cases have the same value, the values of the strain rate and the viscous resistance can be used to calculate the appropriate value of the viscous resistance corresponding to the strain and strain rate generated in the element by interpolation. it can. In the present embodiment, interpolation is performed by such simple primary interpolation. In this way, the viscous resistance is determined for each element by the strain and strain rate generated in that element, thereby expressing the material physical properties under the conditions of appropriate strain and strain rate for each element.

Specifically, as shown in FIGS. 10A, 10B, and 10C, a golf ball model 10 is provided with a cylindrical aluminum hollow of 200 g (same as the weight of a golf club head) as an impacting object. Stick model 20 speed 45m / s
The physical properties of the golf ball model 10 at the time of the collision are analyzed by the finite element method and a simulation is performed. As a result, each element 11 of the golf ball model 10 in a predetermined time from the collision of the impacting object to the predetermined time later.
Calculates the amount of strain that occurs in.

Mass is distributed to each node constituting one element, each node is replaced with a mass point, and the velocity of each node is divided by the velocity of the mass point to obtain the total velocity.
That is, the velocity Vbi of the ball after impact (i = x,
y, z) is calculated by the following equation (4).
The total momentum of the ball is considered to be the sum of the momentums of all masses, and the total momentum divided by the total weight is the ball ejection velocity Vbi.
It is defined as

(Formula 4) Here, N and M are the number of nodes of the ball, the total mass, Vni,
Mn is the n-th translational velocity, the mass of the element containing the node divided by the number of nodes contained in that element.

The circular surface 20a of the aluminum hollow bar model 20 is the collision surface, and the collision surface is flat.
It is a flat frontal collision with the golf ball model 10. In addition, the circular surface 2 of the aluminum hollow bar model 20
The center point 20b of 0a is the first golf ball model 10
The collision is made so that it becomes a collision point with.

By the above method, the speeds of the aluminum hollow rod model 20 and the golf ball 10 before and after the collision are calculated, and the coefficient of restitution of the golf ball 10 is calculated and predicted from the respective speeds and weights.

Thus, the viscous resistance is determined as a function of strain and strain rate, and the relationship among the three factors of strain, strain rate, and viscous resistance is input to the golf ball model, and simulation is performed by the finite element method analysis. Therefore, the viscous resistance corresponding to each strain and strain rate of each element can be calculated, and the characteristics of the viscoelastic material having strong nonlinearity can be accurately predicted even under the high-speed large-deformation condition.

Therefore, the performance such as the coefficient of restitution of a golf ball model made of a viscoelastic material, the dynamic behavior, etc., under the same conditions as when hit by an actual golf club head,
It can be grasped accurately without making a prototype.

Further, as shown in FIG. 11, the phase angle δ is calculated from the stress-strain curve using the following equation (5), and the loss coefficient (tan δ) can be calculated from the phase angle δ. it can. (Equation 5) δ = sin ⁻¹ ((σ _a −σ _b ) / σ _max ) −−− (5)

(Measurement of Material Properties by Split Hopkinson Bar Testing Machine) FIG. 12 is a schematic front view showing a split Hopkinson bar testing machine improved for measuring viscoelastic materials.

The split Hopkinson rod tester shown in FIG. 12 includes a striking rod 51, an input rod 53 and an output rod 55, which are arranged on a straight line. Input rod 53
The first strain gauge 57 and the second strain gauge 59.
Is attached. A third strain gauge 61 and a fourth strain gauge 63 are attached to the output rod 55. The rear end 53a of the input rod 53 and the front end 55a of the output rod 55
A columnar test piece 70 is sandwiched between and.

The test piece 70 may be formed by molding the viscoelastic material to be measured into the shape of the test piece, or may be cut out from the molded product of the viscoelastic material. In this embodiment, the length of the test piece 70 (that is, the distance between the rear end 53a of the input rod 53 and the front end 55a of the output rod 55) is 4 mm, and the cross-sectional diameter of the test piece 70 is 18 mm.

The striking rod 51, the input rod 53 and the output rod 55 are cylinders made of polymethylmethacrylate, and have a cross-sectional diameter of 20 mm, a longitudinal elastic modulus of 5300 MPa and a specific gravity of 1.
It is 19. The striking rod 51 has a length of 100 mm.
The length of the input rod 53 and the output rod 55 (hereinafter, the input rod 53 and the output rod 55 are also referred to as “stress rods”) is 2
It is 000 mm. The first strain gauge 57 is the input rod 53.
It is attached at a position 900 mm from the rear end 53a, and the second strain gauge 59 is attached at a position 600 mm from the rear end 53a of the input rod 53. In addition, the third strain gauge 61 is connected to the front end 55a of the output rod 55 by 300
It is attached at the position of mm and the fourth strain gauge 6
3 is attached at a position 600 mm from the front end 55a of the output rod 55.

As described above, the split Hopkinson rod tester of FIG. 12 has the striking rod 51, the input rod 53, and the output rod 5.
5 is a synthetic resin made of polymethylmethacrylate, the striking rod 51 and the input rod 53 are as large as 2000 mm, and the first strain gauge 57 and the rear end 5 of the input rod 53.
3a and the distance between the second strain gauge 59 and the rear end 53a of the input rod 53 are large, so that strain, strain rate, and stress of viscoelastic materials such as crosslinked rubber used for golf balls can be measured. Suitable for

The above-mentioned first strain gauge 57, second strain gauge 59, third strain gauge 61, and fourth strain gauge 6
A strain gauge for uniaxial plastic is used as 3, and in this embodiment, KFP-5-350 manufactured by Kyowa Denki Co., Ltd.
-C1-65 is used to attach the input rod 53 and the output rod 55 to the above-mentioned positions. The mounting positions of the first strain gauge 57 to the fourth strain gauge 63 on the input rod 53 and the output rod 57 are on the same line in the length direction.

In the measurement of strain, strain rate and stress by this split Hopkinson rod tester, first, the striking rod 51 collides with the front end 53a of the input rod 53 at a predetermined velocity. As a result, an incident distorted wave is generated, and this incident distorted wave travels toward the rear end 53a of the input rod 53. A part of the incident strain wave is reflected at the rear end 53a of the input rod 53 and becomes a reflected strain wave, which proceeds toward the front end 53b of the input rod 53. A part of the incident strain wave passes through the test piece 70 from the rear end 53a of the input rod 53, further propagates to the output rod 55 to become a transmitted strain wave, and the rear end 55b of the output rod 55.
Head towards.

The incident strain wave is actually measured by the first strain gauge 57 and the second strain gauge 59. The measured incident distorted wave is passed through a low-pass filter for 10
High frequencies above kHz are removed. Further, the time history of the incident distorted wave is subjected to zero correction with its baseline value set to zero. First strain gauge 57 thus obtained
And the time axis strains in the second strain gauge 59 are Fourier transformed to obtain the frequency axis strain. The first strain gauge 57 and the second strain gauge 5
From the frequency axis strain at 9, the transfer function is derived. While considering the ratio (X1: X2) of the distance X1 between the first strain gauge 57 and the rear end 53a of the input rod 53 and the distance X2 between the second strain gauge 59 and the rear end 53a of the input rod 53, Based on the transfer function, the rear end 53 of the input rod 53
The frequency axis strain at a is estimated. This frequency axis distortion is subjected to inverse Fourier transform, so that the input rod 53
The time axis strain (time history of strain) εi of the incident strain wave at the rear end 53a is obtained.

Similarly, the reflected strain wave reflected by the rear end 53a of the input rod 53 and traveling toward the front end 53b is measured by the second strain gauge 59 and the first strain gauge 57. From the measured reflected distortion wave, the rear end 53 of the input rod 53
The time base strain (time history of strain) εr of the reflected strain wave at a is obtained.

Further, the third strain gauge 61 of the output rod 55
Also, the transmitted strain wave propagated to the output rod 55 through the test piece 70 is measured by the fourth strain gauge 63, and the time-axis strain of the transmitted strain wave at the front end 55a of the output rod 55 is measured from the measured transmitted strain wave ( Time history of strain) ε
t is obtained.

From the thus obtained εi, εr and εt, the strain rate ε ′ of the test piece 70 is calculated by the following mathematical expression (6). (Equation 6) ε '= (C ₀ / L) · (εi-εr-εt) = ((E / ρ) ^1/2 / L) · (εi-εr-εt) ----- (6) (Equation 6) In (6), C ₀ represents the propagation velocity (m / s) of the strain wave in the input rod and the output rod (stress rod), and L
Represents the length (m) of the test piece, E represents the longitudinal elastic modulus (N / m ² ) of the stress rod, and ρ is the density of the stress rod (kg / m ³ ).
Represents)

Further, the strain ε of the test piece 70 is calculated from εi, εr and εt by the following mathematical expression (7).

(Formula 7) (In Equation (7), C ₀ represents the propagation velocity (m / s) of the strain wave in the input rod and the output rod (stress rod), and L
Represents the length (m) of the test piece, E represents the longitudinal elastic modulus (N / m ² ) of the stress rod, and ρ is the density of the stress rod (kg / m ³ ).
Represents)

Further, the stress σ of the test piece 70 is calculated from εi, εr, and εt by the following equation (8). (Equation 8) σ = (E · A / (2As)) · (εi + εr + εt) = (E · D ² / (2 (Ds) ² )) · (εi + εr + εt) −−− (8) (in Equation (8) , E is expressed on the vertical elasticity coefficient of the stress rods made of the input bar and output bar (N / m ^2), a represents the cross-sectional area of the stress rods (m ^2), As the cross-sectional area of the test piece (m ² ), D represents the diameter (m) of the stress rod, and Ds
Represents the diameter (m) of the test piece)

The strain time history of the test piece 70 thus obtained is shown in the graph of FIG. As shown in FIG. 13, the curve is gentle for a while after the peak P, but thereafter becomes uneven. Select a point S at a smooth stage after the peak P, draw a tangent to the curve at this point S, derive the relaxation time λ from the intersection of this tangent and the time axis, and calculate the curve by the following formula (9). Is a curve after the point S, the entire strain time history can be a smooth curve (shown by a dotted line in FIG. 13). This makes it possible to eliminate the influence of noise on the finally obtained material properties. (Formula 9) ε (t) = ε ₀ · e ^{−t / λ} −−− (9) (In the formula (9), ε ₀ represents the strain at the contact)

Similarly, the following equation (10) makes it possible to make the entire stress time history into a gentle curve, thereby eliminating the influence of noise on the finally obtained material properties. (Formula 10) σ (t) = σ ₀ · e ^{−t / λ} −−− (10) (In Formula (10), σ ₀ represents stress at the contact)

As described above, the strain time history and the stress time history of the test piece 70 are corrected.

By the above method, the time history data of strain, the time history data of strain rate, and the time history data of stress at high speed large deformation are obtained by the split Hopkinson bar tester.

Examples of the present invention and comparative examples will be described in detail below. First, a golf ball was manufactured using a material containing urethane as a main component. A material containing urethane as a main component was compression-molded at 160 ° C. for 30 minutes with a mold having a diameter of 42.8 mm to obtain a one-piece ball.

Strain, strain rate, and stress of the material containing urethane as a main component were measured by the split Hopkinson rod tester described above, and the impact velocity of the striking rod was 7 m / s.
It measured about 4 patterns of m / s, 20 m / s, and 25 m / s (room temperature 23 degreeC, relative humidity 50%). The maximum strain and maximum strain rate during measurement are shown below. Collision velocity 7m / s (maximum strain 0.12, maximum strain velocity 1378 / s) Collision velocity 14m / s (maximum strain 0.24, maximum strain velocity 2703 / s) Collision velocity 20m / s (maximum strain 0.35, Maximum strain velocity 3898 / s) Collision velocity 25m / s (maximum strain 0.43, maximum strain velocity 4716 / s) Since the measurement is performed with four patterns of collision velocity, the value of the phase angle δ at each collision velocity is as follows. It shows in Table 1.

[0079]

[Table 1]

(Example) Time history data of strain, strain rate, and stress measured by a split Hopkinson bar tester using a material containing urethane as a main component and a viscoelastic model similar to that of the first embodiment. Using and, the simulation was performed considering the viscous resistance. In the simulation method according to the first embodiment, the relationship between the strain, the strain rate, and the longitudinal elastic coefficient is further input to the product model, and the change in the longitudinal elastic coefficient due to the difference in the strain and the strain rate is taken into consideration. The viscous resistance is changed when the strain state generated in the viscoelastic material is in the loaded state and in the unloading or restoring state. The phase angle δ predicted by the simulation of the example is shown in Table 1 above.

(Comparative Example) A simulation was performed using a conventional viscoelastic material with a constant value at which the viscous resistance did not change, with reference to the loss coefficient measured by a split Hopkinson bar tester using a material containing urethane as a main component as a test piece. It was
The phase angle δ predicted by the simulation of the comparative example is shown in Table 1 above.

According to the finite element method analysis, 200 g of a golf ball made of a material containing urethane as a main component is used.
Simulating the performance of the golf ball and the deformed state of the material when the aluminum hollow rod model (same as the weight of the head) is collided at three initial velocities of 35, 40 and 45 m / s, The coefficient of restitution in the analysis of the golf ball was calculated by the method described above. Table 2 below shows the coefficient of restitution in the analysis of Examples and Comparative Examples of the golf balls made of the above-mentioned urethane-based material.

[0083]

[Table 2]

Further, according to an experiment using a real golf ball molded from the above-mentioned material containing urethane as a main component,
The value of the coefficient of restitution of the actual golf ball was measured by the following method. Table 2 shows the difference (%) between the coefficient of restitution in the analysis of the above Examples and Comparative Examples and the coefficient of restitution in the experiment using the following actual products.

(Measurement of Restitution Coefficient by Experiment Using Actual Golf Ball) As a method of measuring the restitution coefficient of a golf ball, a golf ball made of the above material was used in place of the head of a golf club at room temperature of 23 ° C. As 20
A hollow rod made of aluminum (0 g (same as the weight of the head)) was made to collide at three speeds of 35, 40 and 45 m / s. The velocities of the hollow rod and the golf ball before and after the collision were measured, and the coefficient of restitution of the golf ball was calculated from the respective velocities and weights. The collision surface of the aluminum hollow rod was flat, and the collision surface was a flat frontal collision with the golf ball. Since the golf ball does not rotate at the time of collision because the collision surface does not have an angle like a golf club head, only the resilience performance of the golf ball can be evaluated.

As shown in Table 1, the phase angle δ in the analysis of the example has almost the same value as the phase angle δ of the experimental result at each collision velocity. On the other hand, the phase angle δ in the analysis of the comparative example
At each collision velocity, a large difference was found with the phase angle δ of the experimental result, and it was confirmed that the example could accurately simulate the actual experimental result.

Further, as shown in Table 2, the restitution coefficient value in the analysis of the examples is the restitution coefficient value in the experiment obtained by actually making a golf ball as a trial, and which hollow bar Even in the case of speed, there was a good agreement. The difference between the example and the actual test result was -4.82% to -5.98%, and it was confirmed that the actual test result was accurately simulated. On the other hand, the value of the coefficient of restitution in the analysis of the comparative example was greatly deviated from the value of the coefficient of restitution in the experiment obtained by actually making a golf ball as a prototype and in the case of any hollow rod speed. . The difference between the comparative example and the actual experimental result was −24.70% to −27.02%, which was significantly different from the experimental result.

From the above, it was confirmed that by performing a simulation in consideration of viscous resistance, the performance of a product made of a viscoelastic material whose material properties show non-linearity under actual use conditions can be accurately predicted by simulation. .

[0089]

As is apparent from the above description, according to the present invention, the viscoelastic resistance is derived by using the viscoelastic model considering the viscosity of the viscoelastic material, and the relationship between strain, strain rate and viscous resistance is calculated. Since the simulation is performed by inputting into the model, it is possible to accurately represent the phenomenon that the physical property value indicating the non-linearity of the viscoelastic material changes non-linearly depending on the speed and size of the deformation of the material.

The values of strain, strain rate, and stress used in the simulation are measured under the condition that the product made of the viscoelastic material is actually used.
It is possible to perform simulations corresponding to various deformation states of the actual viscoelastic material.

Therefore, depending on the deformed state of the material, strain,
Even in the case of a product made of a viscoelastic material in which the relationship between strain rates changes and the material properties such as the loss coefficient show nonlinearity, the performance and the dynamic behavior of the product can be accurately analyzed by simulation.

In an actual golf club head, an actual impact test can be accurately simulated under the conditions of strain and strain rate equivalent to those when the golf ball is impacted, and in a state close to the actual impact. The performance and deformation behavior of the golf ball can be accurately predicted. This makes it easy to understand the physical properties of the material that affect the performance of the golf ball, which not only helps improve the performance of the product,
In the design stage of a golf ball, it is possible to reduce the number of actual ball trials and reduce the cost and time required for trial production.

[Brief description of drawings]

FIG. 1 is a diagram showing time history data of strain ε.

FIG. 2 is a diagram showing time history data of strain rate ε ′.

FIG. 3 is a diagram showing time history data of stress σ.

FIG. 4 is a diagram showing a strain-stress diagram and a method of calculating a longitudinal elastic modulus.

FIG. 5 is a diagram showing a two-element Voigt model used as a viscoelastic model of the present embodiment.

FIG. 6 is a diagram showing time history data of viscous resistance.

FIG. 7 is a diagram showing a relationship between strain and viscous resistance.

FIG. 8 is a diagram showing a relationship between strain and strain rate.

FIG. 9 is a diagram showing how a golf ball model is divided by meshes.

FIG. 10 shows a collision state of an aluminum hollow rod model with a golf ball model, (A) before collision, and (B).
Is a diagram during a collision, and (C) is a diagram after the collision.

FIG. 11 is a diagram showing a method of calculating a loss coefficient.

FIG. 12 is a schematic front view showing a split Hopkinson bar tester.

FIG. 13 is a diagram showing a state of strain time history of a test piece.

[Explanation of symbols]

10 golf club models 11 elements 12 nodes 20 Aluminum hollow bar model

─────────────────────────────────────────────────── ─── Continuation of the front page (56) References JP 2002-55034 (JP, A) JP 2001-13059 (JP, A) (58) Fields investigated (Int.Cl. ⁷ , DB name) G01N 19 / 00 G01N 3/00-3/62 G06F 19/00 JISST file (JOIS)

Claims (7)

(57) [Claims] 1. The strain, strain rate, and momentary values of stress occurring in the viscoelastic material are measured under measurement conditions that assume the actual use state of a product made of a viscoelastic material, and the strain, strain rate, From the time history data of stress and the viscoelastic model considering the viscosity of the viscoelastic material, derive the time history data of the viscous resistance of the viscoelastic material, and use the product made of the viscoelastic material as the product model to be analyzed. Set the strain, strain rate,
Input the relationship of viscous resistance, perform a simulation considering the change of viscous resistance due to the difference of strain and strain rate, and predict the performance of the product model composed of the viscoelastic material. Simulation method for product performance prediction. 2. The simulation method for predicting the performance of a product made of a viscoelastic material according to claim 1, wherein the viscoelastic material is a material whose physical properties are non-linear. 3. The simulation method for predicting the performance of a product made of a viscoelastic material according to claim 1, wherein the simulation is performed by a finite element method analysis. 4. The strain, strain rate, and stress are measured by a split Hopkinson bar tester.
A simulation method for predicting the performance of a product made of the viscoelastic material according to claim 3. 5. The maximum value of the strain generated in the viscoelastic material during the measurement is in the range of 0.05 to 0.50, and / or the maximum value of the strain rate is 500 / s to.
A simulation method for predicting the performance of a product made of the viscoelastic material according to any one of claims 1 to 4, which is in the range of 10000 / s. 6. The viscoelastic material is a golf ball material, and the product model is a golf ball.
A simulation method for predicting the performance of a product made of the viscoelastic material according to claim 5. 7. The viscoelastic material according to claim 6, wherein the behavior of the golf ball at the time of the collision is predicted by simulating a collision phenomenon between the golf ball and an impacting object assuming a golf club head. Simulation method for product performance prediction.