CONSTANT STRESS SOLID DISK ROTOR OF FLYWHEEL FOR FLYWHEEL ENERGY STORAGE SYSTEM AND DESIGN METHOD THEREOF

Abstract

A constant stress solid disk rotor of a flywheel has an outer shape having a plane-symmetric upper surface and lower surface, an outer circumferential radius b, and a rotation center thickness h.sub.0, and includes a thickness decreasing region which decreases monotonously in thickness from a rotation center to a connection radius a and a constant thickness region located on an outer edge of the thickness decreasing region and having a constant thickness h.sub.a from the connection radius a to the outer circumferential radius b. Shape parameters including the outer circumferential radius b, the rotation center thickness h.sub.0, the connection radius a, and the outer edge thickness h.sub.a satisfy an equation below. Here, ? is a Poisson's ratio of a rotor material.

[00001]$\frac{a}{b}=\sqrt{\frac{1}{2}\left(-\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)+\sqrt{{\left(\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)\right)}^2+\frac{4\left(3+v\right)}{1-v}}\right)}$

Claims

1. A constant stress solid disk rotor of a flywheel, having an outer shape having an upper surface and a lower surface which are plane-symmetric with respect to a single-center rotation plane perpendicular to a rotation axis, an outer circumferential radius b, and a rotation center thickness h.sub.0, the constant stress solid disk rotor of the flywheel having a shape including a thickness decreasing region which decreases monotonously in thickness from a rotation center to a connection radius a and a constant thickness region located on an outer edge of the thickness decreasing region and having a constant thickness h.sub.a from the connection radius a to the outer circumferential radius b, wherein shape parameters including the outer circumferential radius b, the rotation center thickness h.sub.0, the connection radius a, and the outer edge thickness h.sub.a satisfy an equation below: $\begin{array}{cc}\frac{a}{b}=\sqrt{\frac{1}{2}\left(-\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)+\sqrt{{\left(\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)\right)}^2+\frac{4\left(3+v\right)}{1-v}}\right)}& \left[\mathrm{Expression}1\right]\end{array}$ where ? is a Poisson's ratio of a rotor material.

2. The constant stress solid disk rotor of the flywheel according to claim 1, wherein the connection radius a and the outer edge thickness h.sub.a are set without depending on a rotation angular velocity.

3. The constant stress solid disk rotor of the flywheel according to claim 1, wherein the thickness decreasing region is formed in a shape in which an in-plane stress of the thickness decreasing region is always invariant entirely in the thickness decreasing region.

4. The constant stress solid disk rotor of the flywheel according to claim 3, wherein the constant thickness region is formed in a shape in which an in-plane stress of the constant thickness region decreases monotonously from a stress value which is invariant in a plane of the thickness decreasing region toward the outer circumferential radius b from the rotation center.

5. The constant stress solid disk rotor of the flywheel according to claim 1, wherein a thickness h of the thickness decreasing region is expressed by an expression below: $\begin{array}{cc}h\left(r\right)=h_0{e}^{\frac{\ln \left(\frac{h_a}{h_0}\right)}{a^2}{r}^2}& \left[\mathrm{Expression}2\right]\end{array}$

6. The constant stress solid disk rotor of the flywheel according to claim 1, wherein when rotating with the thickness decreasing region producing an in-plane stress ?.sub.a, a rotation angular velocity ? is expressed by an expression below: $\begin{array}{cc}?=\sqrt{-\frac{{?}_a\ln \left(\frac{h_a}{h_0}\right)}{?a^2}}& \left[\mathrm{Expression}3\right]\end{array}$ where ? is a density of the rotor material.

7. The constant stress solid disk rotor of the flywheel according to claim 1, wherein assuming that a yield strength of the rotor material is ?.sub.y, a limit energy density D.sub.FR is expressed by an equation below: $\begin{array}{cc}{D}_{\mathrm{FR}}=\frac{-\left({\left(\frac{a}{b}\right)}^4\left(\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)-\left(\frac{h_a}{h_0}-1\right)\right)+\frac{1}{2}\left(\frac{h_a}{h_0}\right){\left(\ln \left(\frac{h_a}{h_0}\right)\right)}^2\left(1-{\left(\frac{a}{b}\right)}^4\right)\right)}{{\left(\frac{a}{b}\right)}^2\left({\left(\frac{a}{b}\right)}^2\left(\frac{h_a}{h_0}-1\right)+\left(\frac{h_a}{h_0}\right)\ln \left(\frac{h_a}{h_0}\right)\left(1-{\left(\frac{a}{b}\right)}^2\right)\right)}?\frac{{?}_y}{?}& \left[\mathrm{Expression}4\right]\end{array}$

8. A method for designing the constant stress solid disk rotor of the flywheel according to claim 1, wherein any three parameters among four parameters of the outer circumferential radius b, the rotation center thickness h.sub.0, the connection radius a, and the outer edge thickness h.sub.a are given, and a remaining one parameter is determined.

9. A method for designing the constant stress solid disk rotor of the flywheel according to claim 7, wherein any three parameters among six parameters of the limit energy density D.sub.FR, a mass, the outer circumferential radius b, the rotation center thickness h.sub.0, the connection radius a, and the outer edge thickness h.sub.a of the constant stress solid disk rotor of the flywheel are given, and remaining three parameters are determined.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0023] FIG. 1 is an explanatory cross-sectional view showing a configuration of a constant stress solid disk rotor of a flywheel of the present invention.

[0024] FIG. 2 is a graph showing an outer shape profile of a constant stress solid disk rotor of a flywheel according to a first embodiment of the present invention.

[0025] FIG. 3 is a graph showing a radial direction limit (fracture) stress distribution of the constant stress solid disk rotor of the flywheel according to the first embodiment of the present invention.

[0026] FIG. 4 is a graph in which a shape basic equation and an energy density (D.sub.FR) equation of a constant stress solid disk rotor of a flywheel according to a second embodiment of the present invention are plotted as a function of a/b and h.sub.a/h.sub.0.

[0027] FIG. 5 is a graph in which a shape basic equation and a rotor mass (m.sub.ab) equation of a constant stress solid disk rotor of a flywheel according to a third embodiment of the present invention are plotted as the function of a/b and h.sub.a/h.sub.0.

DESCRIPTION OF EMBODIMENTS

[0028] Embodiments of a constant stress solid disk rotor of a flywheel of the present invention will be described with reference to the drawings. However, in these drawings, the relationship between thicknesses and planar dimensions, the thickness ratio between respective layers, and the like are drawn in an exaggerated manner for ease of understanding. Identical members are denoted by an identical reference character, and repeated description will be omitted.

[0029] First, a structure of a constant stress solid disk rotor of a flywheel of the present invention, relational expressions to be established by the structure, effects to be obtained, and the like will be described using the cross-sectional view of FIG. 1.

[0030] FIG. 1 is a cross-sectional view of a constant stress solid disk rotor 1 (hereinafter simply referred to as the rotor 1) of a flywheel when cut along a rotation axis 11 at any rotation angle. The same shape as the shape of FIG. 1 is obtained when the rotor 1 is cut at any rotation angle. A line 12 represents a horizontal center plane of the rotor 1, and is perpendicular to the rotation axis 11.

[0031] As shown in FIG. 1, the rotor 1 is a disk having an outer circumferential radius b (a diameter 2b) with the rotation axis 11 placed at the center. The rotor 1 includes three surfaces, that is, an upper surface 13A, a lower surface 13B, and a side surface 14. The upper surface 13A and the lower surface 13B are plane-symmetric (line-symmetric in the cross-sectional view) with respect to the above-described horizontal center plane 12.

[0032] The material of the rotor 1 is a solid material that exhibits isotropy or weak anisotropy. Examples thereof include, but are not limited to, metal, ceramics, polymer (solid resin), and the like.

[0033] Hereinafter, it is assumed that the material to be used for the rotor 1 has a density ?, a yield strength ?.sub.y, and a Poisson's ratio ?, which will be collectively referred to as material parameters.

[0034] Assuming that a thickness between the upper surface 13A and the lower surface 13B at a distance of the radius r from the rotation axis 11 is h(r), the rotor 1 is segmented into a thickness decreasing region 15A in which a thickness h decreases toward the outer circumference and a constant thickness region 15B located on the outer edge of the thickness decreasing region 15A and having the constant thickness h.sub.a. A radius that gives the boundary between the thickness decreasing region 15A and the constant thickness region 15B is a connection radius a (a diameter 2a).

[0035] A thickness at the center of the rotor 1 (rotation center thickness) is represented by h.sub.0, and a thickness of the constant thickness region 15B (outer edge thickness) is represented by h.sub.a. These thicknesses h.sub.0, h.sub.a, the above-described outer circumferential radius b, and the connection radius a will be hereinafter collectively referred to as shape parameters.

[0036] In the rotor 1, the thickness h of the thickness decreasing region 15A (0?r?a) is expressed by Function Expression (1) below:

[00006]$\begin{array}{cc}\left[\mathrm{Expression}5\right]& ?\\ h\left(r\right)=h_o{e}^{\frac{\ln \left(\frac{h_a}{h_0}\right)}{a^2}{r}^2}& \left(1\right)\end{array}$

It is apparent from FIG. 1 that h.sub.a<h.sub.0 holds, and therefore, a coefficient of r.sup.2 in Expression (1) above:

[00007]$\begin{array}{cc}\frac{\ln \left(\frac{h_a}{h_0}\right)}{a^2}& \left[\mathrm{Expression}6\right]\end{array}$

is not a positive but negative constant. Therefore, h is a decreasing function that approaches asymptotically to zero as r increases.

[0037] It is understood that the rotation angular velocity ? is not included in Function Expression (1) above, and therefore, the shape of the thickness decreasing region 15A of the rotor 1 does not depend on the rotation angular velocity.

[0038] When substituting r=a into Expression (1) above, h(a)=h.sub.a holds, and this value agrees with the thickness of the constant thickness region 15B. In other words, it is understood that surfaces of the upper and lower surfaces (13A, 13B) at the connection radius a spot between the thickness decreasing region 15A and the constant thickness region 15B have no steps and are continuous as in FIG. 1.

[0039] In the rotor 1, the shape parameters h.sub.0, h.sub.a, b, and a are associated with one another by Shape Basic Equation (2) below. This point is also one of distinctive characteristics of the rotor of the present invention.

[00008]$\begin{array}{cc}\left[\mathrm{Expression}7\right]& ?\\ \frac{a}{b}=\sqrt{\frac{1}{2}\left(-\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)+\sqrt{{\left(\frac{2}{1-v}\left(1+v-\frac{2}{\ln \left(\frac{h_a}{h_0}\right)}\right)\right)}^2+\frac{4\left(3+v\right)}{1-v}}\right)}& \left(2\right)\end{array}$

However, 0<a/b<1 and 0<h/h.sub.0<1 hold, and ? is the Poisson's ratio of the rotor material.

[0040] Shape Basic Equation (2) above can also be expressed as Equation (2).

[00009]$\begin{array}{cc}\left[\mathrm{Expression}8\right]& ?\\ \frac{h_a}{h_0}={\exp \left(\frac{4}{\left(1-v\right){\left(\frac{a}{b}\right)}^2+2\left(1+v\right)-\frac{3+v}{{\left(\frac{a}{b}\right)}^2}}\right)}^2& {\left(2\right)}^{}\end{array}$

[0041] The rotor 1 satisfies Shape Basic Equation (2), as a result of which an effect is obtained in which a rotation plane stress ?.sub.a (=a circumferential stress ?.sub.a?=a radial stress ?.sub.ar) produced in the thickness decreasing region 15A (0?r?a) does not depend on an in-plane position (the radius r), but ?.sub.a is constant everywhere, that is, ideally averaged.

[0042] The rotation plane stress ?.sub.a in the thickness decreasing region 15A and the rotation angular velocity ? of the rotor are associated with each other in Expression (3) below.

[00010]$\begin{array}{cc}\left[\mathrm{Expression}9\right]& ?\\ {?}_a\left(?\right)=\frac{-{\mathrm{??}}^2{a}^2}{2\ln \left(\frac{h_a}{h_0}\right)}& \left(3\right)\end{array}$

This relationship is established until a moment at which the rotor 1 fractures, in other words, until the rotation plane stress ?.sub.a of the thickness decreasing region 15A reaches the yield strength ?.sub.y of the rotor material. That is, assuming that a fracture rotation angular velocity when the rotor yields is ?.sub.y, Expression (3) below is established.

[00011]$\begin{array}{cc}\left[\mathrm{Expression}10\right]& ?\\ {?}_y=\frac{-{\mathrm{??}}_y^2{a}^2}{2\ln \left(\frac{h_a}{h_0}\right)}& {\left(3\right)}^{}\end{array}$

[0043] On the other hand, as to a rotation plane stress ?.sub.b produced in the constant thickness region 15B (a?r?b) of the rotor 1, a circumferential stress ?.sub.br is expressed by Expression (4) below.

[00012]$\begin{array}{cc}\left[\mathrm{Expression}11\right]& ?\\ {?}_{\mathrm{br}}=\frac{\left(3+v\right)}{8}\frac{-{?}_a}{2{\left(\frac{a}{b}\right)}^2\ln \left(\frac{h_a}{h_0}\right)}\left(1-\left(\frac{1-v}{3+v}\right){\left(\frac{a}{b}\right)}^4-{\left(\frac{r}{b}\right)}^2-\frac{-\left(\frac{1-v}{3+v}\right){\left(\frac{a}{b}\right)}^4}{{\left(\frac{r}{b}\right)}^2}\right).Math.& \left(4\right)\end{array}$

A radial stress ?.sub.b? is expressed by Expression (5) below.

[00013]$\begin{array}{cc}\left[\mathrm{Expression}12\right]& ?\\ {?}_{b?}=\frac{\left(3+v\right)}{8}\frac{-{?}_a}{2{\left(\frac{a}{b}\right)}^2\ln \left(\frac{h_a}{h_0}\right)}\left(1-\left(\frac{1-v}{3+v}\right){\left(\frac{a}{b}\right)}^4-\frac{1+3v}{3+v}{\left(\frac{r}{b}\right)}^2-\frac{-\left(\frac{1-v}{3+v}\right){\left(\frac{a}{b}\right)}^4}{{\left(\frac{r}{b}\right)}^2}\right).Math.& \left(5\right)\end{array}$

Expressions (4) and (5) above include the rotation plane stress of the thickness decreasing region 15A.

[0044] The above-described expressions of ?.sub.br and ?.sub.b? are functions in which Ga is given at r=a and values decrease monotonously toward the outer circumference r=b as shown in FIG. 3 which will be mentioned later, and therefore, the maximum stress indicated by the rotor 1 is the value ?.sub.a of the rotation plane stress of the thickness decreasing region 15A at any rotation angular velocity. Consequently, the rotor 1 fractures when the rotation plane stress of the thickness decreasing region 15A becomes the yield strength of the rotor material, that is, when ?.sub.a=?.sub.y holds.

[0045] Calculating based on this knowledge, the limit energy density D.sub.FR of the rotor 1 is finally described by Expression (6) below only using the shape parameters and the material parameters (?.sub.y and ?).

[00014]$\begin{array}{cc}\left[\mathrm{Expression}13\right]& ?\\ D_{\mathrm{FR}}=K_{\mathrm{SF}}\frac{{?}_y}{?}.Math.& \left(6\right)\end{array}$

However, K.sub.SF in Expression (6) is expressed by an expression below.

[00015]$\begin{array}{cc}{K}_{\mathrm{SF}}=\frac{-\left({\left(\frac{a}{b}\right)}^4\left(\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)-\left(\frac{h_a}{h_0}-1\right)\right)+\frac{1}{2}\left(\frac{h_a}{h_0}\right){\left(\ln \left(\frac{h_a}{h_0}\right)\right)}^2\left(1-{\left(\frac{a}{b}\right)}^4\right)\right)}{{\left(\frac{a}{b}\right)}^2\left({\left(\frac{a}{b}\right)}^2\left(\frac{h_a}{h_0}-1\right)+\left(\frac{h_a}{h_0}\right)\ln \left(\frac{h_a}{h_0}\right)\left(1-{\left(\frac{a}{b}\right)}^2\right)\right)}& \left[\mathrm{Expression}14\right]\end{array}$

Hereinafter, Expression (6) above will be referred to as a limit energy density equation.

[0046] At last, a mass m.sub.AB of the rotor 1 is described as Expression (7) below.

[00016]$\begin{array}{cc}\left[\mathrm{Expression}15\right]& ?\\ m_{\mathrm{AB}}=\frac{??h_0{a}^2}{\ln \left(\frac{h_a}{h_0}\right)}\left(\frac{h_a}{h_0}-1\right)+??h_a\left(b^2-a^2\right).Math.& \left(7\right)\end{array}$

The first term on the right side of Expression (7) corresponds to the mass of the thickness decreasing region 15A, and the second term corresponds to the mass of the constant thickness region 15B. Hereinafter, Expression (7) will be referred to as a rotor mass equation.

[0047] Effects of the rotor 1 having the above configuration will be described. As is apparent with reference to Expression (1), the thickness h(r) of the thickness decreasing region of the rotor of the embodiment of the present invention (corresponding to t(r) in Patent Literature 1) is invariant with respect to the rotation angular velocity ? including the fracture rotation angular velocity ?.sub.y, and is determined only by the shape parameters. Therefore, it can be said that the rotor 1 according to the present embodiment solves the problem in that the structure of the rotor is indeterminate which is the first problem in the conventional technology (Patent Literature 1; the same applies below).

[0048] As is also apparent with reference to Shape Basic Equation (2), in the rotor 1 of the present embodiment, the outer circumferential radius b and the connection radius a are definitely associated with each other in Expression (2), and are not indeterminate like R in Patent Literature 1 (corresponding to b in the present invention). Therefore, it can be said that the problem in that it is difficult to analytically predict the stress distribution and the limit energy density D.sub.FR which is the second problem in the conventional technology is solved.

[0049] Thus, the rotor 1 of the present embodiment achieves a state in which the rotation stress distribution in the thickness decreasing region is made completely flat in addition to the above-described two effects, so that the peak of the rotation stress is ideally reduced. As a result, optimization of the structural parameters and maximization of the energy density are achieved. In other words, it can be said that the rotor 1 of the present embodiment solves the problem in that it is not easy to obtain an optimum structure which is the third problem in the conventional technology.

[0050] Hereinafter, specific examples in accordance with an actual design will be described using the rotor 1 and its relational expressions. Any material that exhibits isotropy or weak anisotropy can be used for the rotor 1 according to the present invention. Herein, a specific material will be described citing a case of using the QCM8 steel (SANYO SPECIAL STEEL Co., Ltd. (Himeji, Hyogo)) evolved from the high-strength SK105 steel as an example. This is merely an example, and anybody can decide the rotor structure according to the present invention by following exactly the same procedure as a procedure which will be described herein even in a case of using another material.

[0051] Typical values of the material parameters of the QCM8 steel are as follows: the density ?=7734 kg/m.sup.3, the Poisson's ratio ?=0.34, and the yield strength ?.sub.y=1.36?10.sup.9 Pa.

Example 1

[0052] Example 1 is an example of giving some of the shape parameters for determining the physical frame of the rotor 1 as the required specifications to determine an optimum structure of the rotor of the present invention. If any three of the four shape parameters (h.sub.0, h.sub.a, b, and a) are designated (as the required specifications), the remaining one is automatically determined because the rotor 1 necessarily needs to satisfy Shape Basic Equation (2) as described above.

[0053] From the perspective of producing a rotor of a flywheel energy storage system, an actually possible combination of the three parameters is (h.sub.0, h.sub.a, and b) because they substantially determine the physical frame of the rotor 1. In view of this point, a description will be provided herein assuming that these (h.sub.0, h.sub.a, and b) are given as the required specifications.

[0054] With the values of the three parameters (h.sub.0, h.sub.a, and b) given, the value of the connection radius a is determined when the three parameters are substituted into Shape Basic Equation (2).

[0055] When the value of the connection radius a is determined, the function h(r) for determining the surface shape of the thickness decreasing region 15A is determined from Expression (1). The structure of the rotor 1 is totally settled accordingly.

[0056] When the structure of the rotor 1 is totally settled, the limit energy density D.sub.FR is calculated by substituting the values of the shape parameter and the material parameters into Relational Expression (6).

[0057] The fracture rotation angular velocity ?.sub.y is obtained from Expression (3).

[0058] The mass m.sub.AB of the rotor 1 can be calculated by Equation (7).

[0059] Suppose that the shape parameters (h.sub.0, h.sub.a, and b) of the rotor 1 have numerical values of: h.sub.0=0.05 m (5 cm); h.sub.a=0.01 m (1 cm); and b=0.15 m (15 cm), the following values and expressions are obtained.

[00017]$a=0.116m\left(11.6\mathrm{cm}\right)$$h\left(r\right)=0.01\exp \left(-119.1r^2\right)\left(m\right)$$D_{\mathrm{FR}}=45\mathrm{Wh}/\mathrm{kg}$${?}_y=6472\mathrm{rad}/s$$m_{\mathrm{AB}}=10.34\mathrm{kg}$

[0060] FIG. 2 is a radial direction surface shape of the cross-section of the rotor 1 of the first embodiment of the present invention determined by the above-described process. It is understood that the respective shape parameters are correctly reflected.

[0061] FIG. 3 shows a radial direction distribution of the rotation stress ? when the rotor 1 of the first embodiment of the present invention reaches the limit energy density. It is clearly recognized that in the thickness decreasing region 15A, the rotation stress is ?.sub.a (=?.sub.br=?.sub.b?)=?.sub.y=constant, and in the constant thickness region 15B, a manner in which the values of the rotation stresses ?.sub.br and ?.sub.b? decrease monotonously from ?.sub.a.

Example 2

[0062] One of property parameters which is most likely to be included in the required specifications of a rotor of an actual flywheel energy storage system is considered to be the limit energy density D.sub.FR. Thus, in Example 2, an example in which the property parameter D.sub.FR is included in the required specifications as a target value is cited.

[0063] When Limit Energy Density Equation (6) above is deformed under conditions that 0<a/b<1 and 0<h/h.sub.0<1, Expression (6) below is obtained.

[00018]$\begin{array}{cc}\left[\mathrm{Expression}16\right]& ?\\ \frac{a}{b}=\sqrt{\frac{1}{2A}\left(-B-\sqrt{B^2-4AC}\right)}.Math.& {\left(6\right)}^{}\end{array}$

Herein, A, B, and C in Equation (6) are expressed by expressions below.

[00019]$\begin{array}{cc}\left[\mathrm{Expression}17\right]& ?\\ A=\left(1-\frac{h_a}{h_0}+\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)\right)\left(1-K_{\mathrm{SF}}\right)-\frac{1}{2}\frac{h_a}{h_0}{\left(\ln \left(\frac{h_a}{h_0}\right)\right)}^2.Math.& {\left(6a\right)}^{}\end{array}$$\begin{array}{cc}B=K_{\mathrm{SF}}\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right).Math.& {\left(6b\right)}^{}\end{array}$$\begin{array}{cc}C=\frac{1}{2}\frac{h_a}{h_0}{\left(\ln \left(\frac{h_a}{h_0}\right)\right)}^2.Math.& {\left(6c\right)}^{}\end{array}$

[0064] In the case of Example 2, the shape parameters (h.sub.0, h.sub.a, b, and a) need to satisfy Shape Basic Equation (2) and Limit Energy Density Equation (6) at the same time. Thus, it is appreciated that if the limit energy density D.sub.FR and any two of the shape parameters are designated as the required specifications, an optimum solution is obtained.

[0065] Consideration of a combination of two parameters which are likely to be selected from among the shape parameters as the required specifications from a viewpoint of producing the rotor of the flywheel energy storage system results in (h.sub.0 and b). This is because they define the physical frame of the rotor 1 most strongly. In view of this point, a description will be provided assuming that D.sub.FR, h.sub.0, and b are given as the required specifications, but the combination of the shape parameters is not limited to (h.sub.0 and b), and any two can be selected.

[0066] The following is a procedure of determining the unknown shape parameters h.sub.a and a. When D.sub.FR is given as a target value, Limit Energy Density Equation (6) is settled. Next, Shape Basic Equation (2) and Limit Energy Density Equation (6) are set up simultaneously to obtain numerical solutions of a/b and h.sub.a/h.sub.0. When substituting the requested values h.sub.0 and b into the solutions, h.sub.a and a are obtained.

[0067] As in Example 1, the fracture rotation angular velocity ?.sub.y of the rotor 1 is calculated from Expression (3), and the mass m.sub.AB is calculated from Expression (7).

[0068] Herein, an optimum design shall be endeavored assuming that the respective requested values of D.sub.FR, h.sub.0, and b are such that D.sub.FR=40 Wh/kg (=40?3600 J/kg), h.sub.0=0.03 m (3 cm), and b=0.15 m (15 cm).

[0069] FIG. 4 is a graph in which Shape Basic Equation (2) and Limit Energy Density Equation (6) of the above-described D.sub.FR value are plotted as curves of a/b vs h.sub.a/h.sub.0. It is obvious that the solutions a/b and h.sub.a/h.sub.0 are obtained from the fact that the two curves cross each other. When actually calculating solutions using a bisection method, a/b=0.641 and h.sub.a/h.sub.0=0.463 hold. When substituting the requested values (the above-described numerical values) of h.sub.0 and b into these relational expressions, the values of h.sub.a and a are calculated.

[0070] Finally, the following values and expressions are obtained.

[00020]$a=0.096m\left(9.6\mathrm{cm}\right)$$h_a=0.014m\left(1.4\mathrm{cm}\right)$$h\left(r\right)=0.03\exp \left(-83.5r^2\right)\left(m\right)$${?}_y=5419\mathrm{rad}/s$$m_{\mathrm{AB}}=9.16\mathrm{kg}$

Example 3

[0071] Another property parameter which is likely to be included in the required specifications of the rotor 1 is considered to be the mass m.sub.AB of the rotor. Thus, an example in which m.sub.AB is included in the required specifications as a target value will be cited as Example 3.

[0072] When deformed under the conditions that 0<a/b<1 and 0<h.sub.a/h.sub.0<1, Rotor Mass Equation (7) above is rewritten as Equation (7) or (7) below, which can be used as a substitute for Expression (7) as appropriate.

[00021]$\begin{array}{cc}\left[\mathrm{Expression}18\right]& ?\\ \frac{a}{b}=\sqrt{\frac{\left(\frac{m_{\mathrm{ab}}}{?h_0?b^2}-\frac{h_a}{h_0}\right)\ln \left(\frac{h_a}{h_0}\right)}{\frac{h_a}{h_0}-1-\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)}}.Math.& {\left(7\right)}^{}\end{array}$$\begin{array}{cc}\frac{a}{b}=\sqrt{\left(\frac{m_{\mathrm{ab}}}{?h_a?b^2}-1\right)\frac{\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)}{\frac{h_a}{h_0}-1-\frac{h_a}{h_0}\ln \left(\frac{h_a}{h_0}\right)}}.Math.& {\left(7\right)}^{}\end{array}$

[0073] In the case of Example 3, the shape parameters (h.sub.0, h.sub.a, b, and a) need to satisfy Shape Basic Equation (2) and Rotor Mass Equation (7) (or (7) or (7)) at the same time. Thus, an optimum structure will be obtained if the rotor mass m.sub.AB and two of the shape parameters are designated as the required specifications.

[0074] From the viewpoint of producing the rotor of the flywheel energy storage system, two shape parameters which are likely to be selected as the required specifications are considered to be h.sub.0 and b as in Example 2. A description will also be provided in Example 3 in view of this point assuming that m.sub.AB, h.sub.0, and b are given as the required specifications. The combination of the shape parameters is not limited to (h.sub.0 and b), and a combination of any two can be selected.

[0075] The following is a procedure of determining the unknown shape parameters h.sub.a and a. When m.sub.AB, h.sub.0, and b are given as target values, Rotor Mass Equation (7) is settled. When setting up Rotor Mass Equation (7) and Limit Energy Density Equation (6) simultaneously to solve the expressions paying attention to the fact that h.sub.0 and b are known numbers, the numerical solutions a/b and h.sub.a/h.sub.0 are obtained. Since h.sub.0 and b are known numbers, h.sub.a and a are obtained.

[0076] As in Example 1 described earlier, the function h(r) of the thickness decreasing region 15A of the rotor 1 can be calculated from Expression (1), the limit energy density can be calculated from Equation (6), the fracture rotation angular velocity ?y can be calculated from Expression (3), and the mass m.sub.AB can be calculated from Equation (7). Herein, an optimum design shall be endeavored assuming that the respective requested values of m.sub.AB, h.sub.0, and b are such that m.sub.AB=12 kg, h.sub.0=0.04 m (4 cm), and b=0.15 m (15 cm).

[0077] FIG. 5 is a graph in which Shape Basic Equation (2) and Rotor Mass Equation (7) of the above-described D.sub.FR value are plotted as curves of a/b vs h.sub.a/h.sub.0 (however, h.sub.0 and b are known numbers). From the fact that the two curves cross each other, it is obvious that the solutions a/b and h.sub.a/h.sub.0 are obtained. When actually calculating solutions using a bisection method, a/b=0.649 and h.sub.a/h.sub.0=0.448 are obtained. When the values of the known numbers h.sub.0 and b are substituted into these relational expressions, the values of h.sub.a and a are calculated.

[0078] Finally, the following rotor specification elements are obtained.

[00022]$a=0.097m\left(9.7\mathrm{cm}\right)$$h_a=0.018m\left(1.8\mathrm{cm}\right)$$h\left(r\right)=0.04\exp \left(-84.8r^2\right)\left(m\right)$${?}_y=5462\mathrm{rad}/s$$D_{\mathrm{FR}}=40\mathrm{Wh}/\mathrm{kg}\left(=40?3600J/\mathrm{kg}\right)$

Other Examples

[0079] For the constant stress solid disk rotor of a flywheel according to the present invention, the limit energy density D.sub.FR and the mass m.sub.AB which are the property parameters as well as one of the shape parameters (e.g., the outer circumferential radius b) can also be designated as the required specifications to determine optimum values of the remaining three shape parameters. In this case, three of the shape basic expression, the limit energy density expression, and the rotor mass expression are set up simultaneously to calculate the remaining three shape parameters.

[0080] This procedure will be simply described citing an example of designating the outer circumferential radius b as the shape parameter. First, Shape Basic Equation (2) and Limit Energy Density Equation (6) are set up simultaneously to calculate numerical solutions of the unknown number a/b and the unknown number h.sub.a/h.sub.0, and the value of b is substituted into the former to settle the value of a.

[0081] Next, the calculated values of h.sub.a/h.sub.0, a, and b are substituted into Rotor Mass Equation (7) to settle the value of h.sub.a, and the value of h.sub.0 is calculated from the value of h.sub.a/h.sub.0 and the value of h.sub.a. In this manner, optimum values of the structural parameters a, h.sub.a, and h.sub.0 which are unknown are all settled.

REFERENCE SIGNS LIST

[0082] 1 constant stress solid disk rotor of flywheel of the present invention [0083] 11 rotation axis [0084] 12 horizontal center plane [0085] 13A upper surface [0086] 13B lower surface [0087] 14 side surface [0088] 15A thickness decreasing region [0089] 15B constant thickness region [0090] b rotor radius (outer circumferential radius) [0091] a connection radius (radius at boundary between thickness decreasing region and constant thickness region) [0092] h.sub.a thickness of constant thickness region [0093] h.sub.0 center thickness [0094] r radius at any spot [0095] h(r) function representing thickness at radius r spot