The vibrational mode frequencies and mode shapes of ancient stone chimes are analyzed and their dependence on stone shapes are discussed. Mode shapes and frequencies of several chime models are determined by using finite element methods, and these show good agreement with mode shapes and frequencies observed in Korean pyeongyeong chime stones using holographic interferometry and experimental modal testing. The dependence of mode shapes and frequencies on vertex angle and base curvature suggests that the geometries used in late Chinese bianqing and Korean pyeongyeong may have been selected to give the best sound.
I. Introduction
Stone chimes have been cherished musical instruments for centuries in several Asian countries. They have taken a number of different forms, but generally they have two legs which meet at a vertex. The longer leg is often called the drum, because it is where the chime is struck, and the shorter leg is called the femur or thigh.
Several sets of stone chimes have been found in ancient Chinese tombs dating from 400 BC and before. A bianqing of 32 stones was found in the tomb of the Marquis Yi of Zheng (the same tomb that contained the renowned set of 65 bells), dating from about 433 BC. Although many stones were found broken, Chinese scholars were able to determine their dimensions, and from these, Lehr determined scaling laws for their various dimensions.1 Stones are scaled in lengths as well as thicknesses.
Another ancient bianqing from the Warring States period (450–221 BC) of the Zhou dynasty was found more or less intact, so that the fundamental frequencies of most of the stones, as well as the dimensions, could be determined.2 Vertex angles of these stones ranged from .
The Korean pyeongyeong consists of a set of 16 hard stones cut in the shape of an inverted letter L. The stones are hung from a wooden stand and struck, near the lower end of the longer side, by a mallet tipped with horn. The upper ridge of the longer side of a stone is about 1.5 times longer than the shorter side. Unlike most Chinese stone chimes, the pyeongyeong stones in a set have nearly the same shape except for the highest stone in which the long edge may be slightly shortened, but are scaled in thickness.3
Old pyeongyeongs are preserved in museums as national treasures, but replicas are built and played in concerts. Earlier we analyzed a pyeongyeong built by Hyungon Kim and played by the National Traditional Orchestra of Korea. The fundamental frequency was shown to be directly proportional to thickness over most of the range.4 The second mode was found to be tuned to 1.5 times the nominal frequency, the third mode to about 2.3 times the nominal frequency, and the fourth mode about 3 times the nominal frequency up to the 12th stone.4
In this paper, we report on studies designed to show how the vibrational frequencies of L-shaped chime stones, such as are found in the Korean pyeongyeong and the Chinese bianqing, depend upon the details of their geometry.
II. Shapes of Pyeongyeong and Bianqing chime stones
The shape of a pyeongyeong stone is said to be the “shape of the heaven that curves to cover the earth.” The 115° angle between the drum part and the femur part and the concave curved base form a smooth L shape. In Korea, the pyeongyeong has been built in this smooth L shape for some five centuries, whereas Chinese bianqing stone chimes have varied in shape from time to time. Old stones from before the Christian era mostly have a straight base, or an arc-shaped base, for example, whereas stones from the Qing Dynasty (1644–1911) are L shaped.
According to Kuttner, stones with a straight base [Figs. 1(a) and 1(b)] are older than those with an arc-shaped base [Fig. 1(c)], and he concluded that the curved base resulted from experiments to improve the acoustical qualities of the stones. Also, the angle between the drum part and the femur part became narrower in stones from the Qing Dynasty and later.5
(Color online) Bianqing stones from different ages: (a) ancient stones from Shanxi and Henan provinces; (b) ancient stones from Henan; (c) stone chime from Marquis Yi tomb (433 BC); (d) chime from 1717 AD.
(Color online) Bianqing stones from different ages: (a) ancient stones from Shanxi and Henan provinces; (b) ancient stones from Henan; (c) stone chime from Marquis Yi tomb (433 BC); (d) chime from 1717 AD.
In this study, we examine the effects of geometry on the tuning and the sound quality of stone chimes, such as the Korean pyeongyeong and the Chinese bianqing. Modal shapes and modal frequencies are estimated by means of finite element methods (FEM) analysis, and these results are compared to modal shapes and frequencies of existing stones, determined by holographic interferometry and by experimental modal testing. Modal shapes and frequencies of stones of different shapes are calculated by finite element methods. Of particular interests are the effects of varying the vertex angle and the curvature of the base.
III. Modeling of stones and finite element analysis
Since the sound spectra and mode shapes of several pyeongyeong replica stones have previously been measured,4 we used these stones as models for finite element analyses. The stone of lowest frequency is tuned to C5 , and we use that stone for (FEM) modeling and designate it as PGC 115, since the vertex angle is 115°. The modal frequencies calculated for this stone are compared to the modal frequencies determined from the sound spectrum, from holographic interferometry, and from experimental modal testing using the (STAR) system. The agreement is very good, as shown in Table I.
Calculated (FEM) and measured mode frequencies in the first gyeong. Measured values are the average of the frequencies obtained by holographic interferometry, experimental modal testing, and analysis of the sound spectra.
Modeno. . | PGC116Calculated databy ANSYS . | PGCAverage ofexperimentaldata . | Error . | ||
---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | ||
1 | 527 | 1.00 | 526 | 1.00 | 0.30 |
2 | 810 | 1.54 | 810 | 1.54 | 0 |
3 | 1233 | 2.34 | 1242 | 2.36 | |
4 | 1625 | 3.08 | 1634 | 3.10 | |
5 | 2203 | 4.18 | 2207 | 4.20 | |
6 | 2704 | 5.13 | 2704 | 5.14 | 0 |
7 | 3554 | 6.74 | 3553 | 6.75 | 0.04 |
8 | 3993 | 7.57 | 3957 | 7.52 | 1.42 |
9 | 4527 | 8.59 | 4489 | 8.53 | 1.32 |
Modeno. . | PGC116Calculated databy ANSYS . | PGCAverage ofexperimentaldata . | Error . | ||
---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | ||
1 | 527 | 1.00 | 526 | 1.00 | 0.30 |
2 | 810 | 1.54 | 810 | 1.54 | 0 |
3 | 1233 | 2.34 | 1242 | 2.36 | |
4 | 1625 | 3.08 | 1634 | 3.10 | |
5 | 2203 | 4.18 | 2207 | 4.20 | |
6 | 2704 | 5.13 | 2704 | 5.14 | 0 |
7 | 3554 | 6.74 | 3553 | 6.75 | 0.04 |
8 | 3993 | 7.57 | 3957 | 7.52 | 1.42 |
9 | 4527 | 8.59 | 4489 | 8.53 | 1.32 |
The upper seventh stone from the tomb of Marquis Yi is supposed to be tuned to C, but it was broken when it was excavated.6 The size of each side, angle, and thickness are reported by the Museum of Hubei Province.6 Since the vertex angle is 163°, the model for this stone is designated as BQC 163. Keeping the principle lengths A, B, E, and F of stone the same, vertex angle was then changed to 140° and 115°, and these models are designated BQC 140 and BQC 115, respectively. FEM was used to determine the mode frequencies and mode shapes of the three bianqing models and the pyeongyeong model. The four models are shown in the upper row in Fig. 2.
Bianqing and pyeongyeong models used to test geometry dependence on tuning
With the same angle , the base lines were changed from resembling the smooth arc of the pyeongyeong (BQ115 in Fig. 2) to a near right angle shape (BQ 115ra), three oblique lines (BQ 115to) and a straight line (BQ 115sb) as shown in the second row in Fig. 2. FEM was again used to determine the mode frequencies and mode shapes of the three new bianqing models.
The commercial program ANSYS Mechanical 8.1 was used for the FEM analysis. Since a gyeong is a flat plate of constant thickness, we adapted the three-dimensional eight-node Shell 93 element for modeling it. The shell element 93 was found to give the most accurate results. For the PGC 115 model, 13 key points were created to make the keypoints for modeling by ANSYS. We chose quadrilateral element shapes to construct the mesh.
Systematic ways of varying the vertex angle and base line were introduced. The first set of 19 alternative gyeong models was constructed by adding 5° increments to the vertex angle from 90°, the shape of an L, to 180°, the shape of a bar. In this series of shapes the baseline curvature was kept the same. Next we estimated the effect of baseline curvature on mode frequencies and shapes as the baseline was systematically varied.3
IV. Results and discussion
Table II shows the dependence of the modal frequencies on vertex angle as calculated by FEM. The models correspond to those in Fig. 2. Also shown are the frequencies relative to the lowest mode of vibration. Note that decreasing the vertex angle lowers the ratios of the higher modes and also brings them nearer to a harmonic relationship. This may have made them more sonorous. The ratios for BQC 115 are reasonably close to those of the Korean pyeongyeong model PGC 115.
Model frequencies and frequency ratios for models with varying angle .
No. . | BQC163 . | BQC140 . | BQC115 . | PGC115 . | ||||
---|---|---|---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | |
1 | 591 | 1.00 | 756 | 1.00 | 972 | 1.00 | 527 | 1.00 |
2 | 1357 | 2.29 | 1378 | 1.82 | 1471 | 1.51 | 809 | 1.54 |
3 | 1558 | 2.63 | 1781 | 2.35 | 2114 | 2.18 | 1232 | 2.34 |
4 | 1607 | 2.72 | 2963 | 3.92 | 3109 | 3.20 | 1624 | 3.08 |
5 | 3356 | 5.67 | 3657 | 4.84 | 4016 | 4.13 | 2202 | 4.18 |
6 | 5022 | 8.64 | 4893 | 6.47 | 5230 | 5.38 | 2701 | 5.13 |
No. . | BQC163 . | BQC140 . | BQC115 . | PGC115 . | ||||
---|---|---|---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | |
1 | 591 | 1.00 | 756 | 1.00 | 972 | 1.00 | 527 | 1.00 |
2 | 1357 | 2.29 | 1378 | 1.82 | 1471 | 1.51 | 809 | 1.54 |
3 | 1558 | 2.63 | 1781 | 2.35 | 2114 | 2.18 | 1232 | 2.34 |
4 | 1607 | 2.72 | 2963 | 3.92 | 3109 | 3.20 | 1624 | 3.08 |
5 | 3356 | 5.67 | 3657 | 4.84 | 4016 | 4.13 | 2202 | 4.18 |
6 | 5022 | 8.64 | 4893 | 6.47 | 5230 | 5.38 | 2701 | 5.13 |
Table III shows the dependence of the modal frequencies on baseline curvature. Also shown are the frequencies relative to the lowest mode of vibration. The differences between the smooth arc (BQC 115), the near right angle (BQC 115ra), and the three oblique lines (BQC 115to) are small, but in the stone with the straight base (BQC 115sb) the second, third, and fourth mode frequencies are raised substantially compared to the fundamental. The sound of this stone model would be expected to be quite different from the others and also different from the pyeongyeong stone PGC 115.
Modal frequencies and frequency ratios for models with varying baseline.
No. . | BQC115ra . | BQC115 . | BQC115to . | BQC115sb . | ||||
---|---|---|---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | |
1 | 940 | 1.00 | 972 | 1.00 | 984 | 1.00 | 941 | 1.00 |
2 | 1403 | 1.49 | 1471 | 1.51 | 1532 | 1.56 | 1521 | 1.62 |
3 | 2036 | 2.17 | 2114 | 2.18 | 2312 | 2.35 | 2971 | 3.16 |
4 | 3122 | 3.32 | 3109 | 3.20 | 3075 | 3.12 | 3807 | 4.04 |
5 | 4031 | 4.29 | 4016 | 4.13 | 3987 | 4.05 | 4010 | 4.26 |
6 | 5179 | 5.51 | 5230 | 5.38 | 5641 | 5.73 | 5166 | 5.49 |
No. . | BQC115ra . | BQC115 . | BQC115to . | BQC115sb . | ||||
---|---|---|---|---|---|---|---|---|
Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | Freq.(Hz) . | Ratioto fund. . | |
1 | 940 | 1.00 | 972 | 1.00 | 984 | 1.00 | 941 | 1.00 |
2 | 1403 | 1.49 | 1471 | 1.51 | 1532 | 1.56 | 1521 | 1.62 |
3 | 2036 | 2.17 | 2114 | 2.18 | 2312 | 2.35 | 2971 | 3.16 |
4 | 3122 | 3.32 | 3109 | 3.20 | 3075 | 3.12 | 3807 | 4.04 |
5 | 4031 | 4.29 | 4016 | 4.13 | 3987 | 4.05 | 4010 | 4.26 |
6 | 5179 | 5.51 | 5230 | 5.38 | 5641 | 5.73 | 5166 | 5.49 |
More details of the dependence of mode shapes and frequencies on vertex angle were provided by varying the vertex angle in five-degree increments from 90° to 180°. Figure 3 shows how the mode shape and frequency of the lowest mode changes as the vertex angle changes from 90° to 180°. Similar plots for the higher modes are shown online in file . Figure 4(a) shows how the mode frequencies depend on the vertex angle. A polynomial trend line calculated for best fit, also shown, is of the form . Figure 4(a) shows the fundamental frequency, and Fig. 4(b) shows the frequency ratios of the higher modes to the fundamental.
(Color online) First mode shapes and frequencies of chimelike models with varying vertex angles. (This figure is intended for color viewing.)
(Color online) First mode shapes and frequencies of chimelike models with varying vertex angles. (This figure is intended for color viewing.)
(Color online) Frequency dependence on vertex angle in chimelike models: (a) fundamental frequency; (b) ratio of higher mode frequencies to the fundamental. (This figure is intended for color viewing.)
(Color online) Frequency dependence on vertex angle in chimelike models: (a) fundamental frequency; (b) ratio of higher mode frequencies to the fundamental. (This figure is intended for color viewing.)
shows the dependence of mode shapes and frequencies on vertex angles (10.4MB). This is a file of the type avi. [URL: http://dx.doi.org/10.1121/1.2218983.1]
shows the dependence of mode shapes and frequencies on vertex angles (10.4MB). This is a file of the type avi. [URL: http://dx.doi.org/10.1121/1.2218983.1]
The frequency ratios of the various partials might be expected to determine the timbre of the sound. A vertex angle of about 115° was adopted by both the makers of Korean pyeongyoung and the Chinese bianqing of the Qing dynasty and later.
From examining existing Chinese bianqing, it appears that the L shape with a gently curved base gradually evolved over several centuries. The 18th century stone in Fig. 1(d) has a vertex angle of approximately 113°, as measured from the photograph. This is very close to the angle of 115° measured in Korean pyeongyeong stones of about the same age, which is not surprising since the earliest pyeongyeong were imported from China. This vertex angle was probably arrived at by attempting to maximize sonority.
The dependence of modal frequency, and especially the ratios of the higher mode frequencies to the fundamental, on chime shape suggests that the ancient Chinese and Korean chime makers tuned their chimes by adjusting the shape of the chimes to obtain the most pleasing sound. This was probably done through centuries in China, and when the Koreans began making stone chimes, they were guided by the Chinese chimes having the best sound. The most important parameter is the vertex angle, with the curvature of the base having a much smaller influence.
V. Conclusion
By means of finite element methods (FEM) analysis, we have shown how changing the vertex angle and the curvature of the base changes the relative frequencies of the modes of vibration of stone chimes (lithophones) and hence the character of the sound. Pyeongyeong stone chimes in Korea have maintained pretty much the same shape during the past five centuries, and the same appears to be true of Chinese bianqing. Our calculations show that relative modal frequencies are changed considerably by changing the vertex angle, but are less sensitive to changes in curvature of the base. The L shape with vertex angle around 115° appears to produce stones with maximum sonority.
The fabricators of ancient lithophones no doubt discovered the same thing, albeit by years and years of trial and error.
The agreement between mode frequencies and mode shapes calculated by FEM and measured by holographic interferometry and by experimental modal testing under playing conditions was found to be quite satisfactory, and this validated the use of FEM to study the effects of varying the model shapes.