Group-Theoretic Exploitations of Symmetry in Novel Prestressed Structures

: In recent years, group theory has been gradually adopted for computational problems of solid and structural mechanics. This paper reviews the advances made in the application of group theory in areas such as stability, form-ﬁnding, natural vibration and bifurcation of novel prestressed structures. As initial prestress plays an important role in prestressed structures, its contribution to structural stiffness has been considered. General group-theoretic approaches for several problems are presented, where certain stiffness matrices and equilibrium matrices are expressed in symmetry-adapted coordinate system and block-diagonalized neatly. Illustrative examples on structural stability analysis, force-ﬁnding analysis, and generalized eigenvalue analysis on cable domes and cable-strut structures are drawn from recent studies by the authors. It shows how group theory, through symmetry spaces for irreducible representations and matrix decompositions, enables remarkable simpliﬁcations and reductions in the computational effort to be achieved. More importantly, before any numerical computations are performed, group theory allows valuable and effective insights on the behavior or intrinsic properties of a prestressed structure to be gained.


Introduction
Symmetry is one of the most common and important features in nature. Different forms of symmetry can be easily observed from a microscopic view to macroscopic view, from atoms and crystals to large-scale space structures. At the same time, innovative methods that consider the inherent symmetry of a system will be systematic and effective [1][2][3]. As an important branch of mathematics and vector algebra, group theory is a powerful tool for systematic analysis on symmetric systems. It simplifies the complex computation process. Moreover, it has a direct and qualitative understanding on the intrinsic properties of the system. Thus, group theory has been extensively applied in many fields such as computational physics, electromagnetics, crystal chemistry [4,5] and molecular vibration [6,7].
Notably, applying group theory into structural analysis can significantly improve computational efficiency [8,9]. In recent years, group theory has been gradually adopted by some researchers from structural engineering to solve the involved computational mechanics of symmetric structures [10][11][12][13][14]. Based on group theory, Zlokoviâc [15] and Zingoni [16] showed how to evaluate irreducible representations and symmetry subspaces for symmetric structures. They also summarized the main advantage of group theory through its applications in structural stability, vibration and control. Kangwai et al. [3] briefly described how to utilize group theory to establish a symmetry-adapted coordinate system and perform static analysis on symmetric structures. They concluded that the Fourier method was a special case of the symmetry method using group theory. During nonlinear buckling analysis or frequency analysis of a symmetric structure, if symmetry subspaces associated with certain irreducible representations are established, the tangential

Group Theory and Its Matrix Representations
A symmetry group G describes all the inherent symmetry of a structure by a set of elements G = {g i , i ∈ [1, τ]}, where each element g i denotes an independent symmetry operation, and τ denotes the total number of symmetry operations. This group must obey the following four group properties: i.
Identity: ∃E ∈ G, for any element g i ∈ G, it satisfies g i · E = E · g i = g i . ii. Inverses: iv. Associativity: ∀g i , g j , g k ∈ G ⇒ (g i · g j ) · g k = g i · (g j · g k ).
where · denotes the multiplication between different elements. In fact, each element g i in the symmetry group describes an independent symmetry operation for the structure [2]. Under every symmetry operation, the structure can be transformed into a new configuration that is equivalent or indistinguishable from the original configuration. For structures with finite nodes and members, all symmetry operations can be divided into the following five types [14,30]: (a) the identity E; (b) rotation C n about an n-fold symmetry axis; (c) reflection σ n along a symmetry plane; (d) inverse S 2 ; and (e) improper rotation S n .
Mathematically speaking, each symmetry operation can be taken as a linear transformation in 2D/3D space, and thus the configuration can be expressed as a matrix form: where the matrices X and X S are the generalized nodal coordinates for the initial configuration and the transformed configuration by the symmetry operation S, and R S is the corresponding transformation matrix. As the base vectors are variable for different coordinate systems, the transformation matrix R S is not unique. Fortunately, group theory can be utilized to neatly construct a symmetry-adapted coordinate system, and to decompose the reducible transformation matrix R S into a series of irreducible representations [23,29]: where the parameters α 1 and α 2 are the coefficients for the first and second irreducible representations Γ (1) and Γ (2) . Generally, a symmetry group G has µ ≥ 1 independent irreducible representations, . The character χ Γ (i) for each irreducible representation keeps constant in different coordinate systems. Specific values associated with different kinds of symmetry operations can be directly read from group theory books [15,28].

Stability Analysis on Prestressable Pin-Jointed Structures
Note that novel prestressed structures described in this study refer to various types of prestressable pin-jointed structures, including cable domes, cable nets, tensegrity structures and cable-strut tension structures [31][32][33]. These structures have innovative and attractive geometries, and often contain internal mechanism modes. Thus, the stability of a prestressable pin-jointed structure, that is, whether it can keep a stable state by the initial prestresses, is the basic and key problem that structural engineers need to solve [29,34,35].

Positive Definiteness of the Tangent Stiffness Matrix for a Prestressed Pin-Jointed Structure
It is known that the potential energy Π R of a stable structure states at a minimum [36]. In other words, for an arbitrary virtual nodal displacement vector δd, the second-order variation of the potential energy δ 2 Π R should satisfy which is expressed by the quadratic form of the tangent stiffness matrix K T of the structure. Therefore, stability analysis of the structure is transformed into the positive definiteness problem of the tangent stiffness matrix [34,35,37]. Guest [38,39] has established a unifying approach for the tangent stiffness matrix of a general prestressed pin-jointed structure: In Equation (4), K E is the elastic stiffness matrix, and K G is the geometric stiffness matrix contributed by the initial prestresses. Unlike conventional prestressed structures with positive definite matrix K E , novel prestressed structures have internal mechanisms M, which come from the null space of the positive semi-definite matrix K E and satisfy For a symmetric structure, the tangent stiffness matrix K T can be block-diagonalized in the symmetry-adapted coordinate system.'According to the great orthogonal theorem [3], the transformation matrix U for expressing generalized nodal displacements and external loads is established: where U (i−h) denotes the symmetry subspace associated with the irreducible representation Γ (i) , µ is the number of irreducible representations of a symmetry group, and l i is the dimension of the irreducible representation Γ (i) . Then, the symmetry-adapted tangent stiffness matrix K T of the structure is where diag [ ] denotes a diagonal matrix, and K represents a typical small-sized block matrix along the diagonal of the matrix K T . Equation (7) shows that the original tangential stiffness matrix is decomposed into a series of block matrices, and the total number of the block matrices is computed by On condition that the structural configuration and initial prestresses are fully symmetric, the stiffness matrices K E and K G can be similarly decomposed into the block-diagonalized ones: where K E and K G are symmetry-adapted elastic stiffness matrix and geometric stiffness matrix. For a positive integer h ∈ [1, l i ], the block stiffness matrices K Therefore, to evaluate the positive definiteness of the original tangent stiffness matrix can be independently evaluated by numerical approaches, e.g., the eigenvalue decomposition where λ (i−h) and Ψ (i−h) are the eigenvalue and eigenvector extracted from the symmetry subspace associated with the irreducible representation Γ (i) . Because any similarity transformations do not alter the eigenvalues of a matrix [16,40], the eigenvalues obtained from Equation (11) belong to the eigenvalues λ of the original tangent stiffness matrix and the corresponding eigenvector Ψ (i−h) can be expressed as where the symmetry subspace (see Equation (6)) reveals that the eigenvector Ψ (i−h) holds the same symmetry associated with the irreducible representation Γ (i) . Subsequently, the full eigenvector matrix Ψ can be obtained by Hence, by using group theory, the positive definiteness problem of the original tangent stiffness matrix is transformed into many independent sub-problems. Consequently, the computational complexity required for eigenvalue decomposition will be significantly reduced.

Necessary Conditions for Structural Stability of Prestressed Pin-Jointed Structures
In the symmetry-adapted coordinate system, Equation (3) can be rewritten as where δd = Mα is the symmetry-adapted nodal displacement vector, M is the symmetry-adapted mechanism mode, and K E M = 0 (refer to Equation (5)). To satisfy the structural stability, the matrix K G should be positive definite or positive semi-definite, i.e., ∑ λ(K G ) > 0. Then, a necessary condition for stability of the symmetric prestressed structure can be obtained [36]. The structure must have at least a single state of self-stress with full symmetry, given by where Γ s is the symmetry representation of the self-stress states. On the other hand, when the structure has only self-stress states without full symmetry (Γ (1) Γ s ), the symmetry of the self-stress states is associated with the i-th type irreducible representation Γ (i) , where 1 < i ≤ µ. Thereafter, the symmetry-adapted geometric stiffness matrix K G cannot be decomposed into block-diagonalized matrices. However, it has certain regularity, and all the diagonal blocks are zero matrices. At this point, the trace of the matrix K G is 0, ∑ λ(K G ) =0. Therefore, this matrix is neither a positive definite nor a negative definite matrix, and does not satisfy the necessary condition for the stability (see Equations (3) and (15)). In this case, the structure must be an unstable (finite) structure.

Example: C 2v Symmetric Cable-Strut Structures
Stability analysis on two simple two-dimensional symmetric cable-strut structures is presented, to explain the basic group-theoretic process. As shown in Figure 1, both of the structures consist of two vertical struts and six cables. The length of struts is L c = 2000 mm, the distance between the struts is 2000 mm, and the minimum distance from the boundary nodes (nodes 5 and 6) to the struts is 1000 mm. The elastic modulus of the cables and that of the struts is E t = 1.9 × 10 5 MPa and E c = 2 × 10 5 MPa; the cross-sectional area of the cables and that of the struts is A t = 500 mm 2 and A c = 3000 mm 2 . Both structures have a single mechanism mode (m = 1) and a self-stress state (s = 1). Note that the initial prestress of the vertical struts is −100 kN, and prestresses of other members can be uniquely determined by the self-stress state.
is 2000 mm, and the minimum distance from the boundary nodes (nodes 5 and 6) to the struts is 1000 mm. The elastic modulus of the cables and that of the struts is Both structures have a single mechanism mode ( 1 m = ) and a self-stress state ( 1 s = ). Note that the initial prestress of the vertical struts is −100 kN, and prestresses of other members can be uniquely determined by the self-stress state.  The structures are C 2v symmetric as they keep equivalent by the following four symmetry operations: Identity, rotation by π, and the reflections along X and Y axes. These four symmetry (6) and (7), the tangential stiffness matrices can be decomposed into four 2 × 2 block matrices, as listed in Table 1. Therefore, the eigenvalues and eigenvectors can be independently extracted from each of the block matrices. Figure 2 plots the eigenvalue decomposition results of the tangent stiffness matrix of the structure shown in Figure 1a. The endpoints of the polyline and the discrete dots are the solutions obtained by the symmetry method and the conventional method, respectively. The structures are 2v C symmetric as they keep equivalent by the following four symmetry operations: Identity, rotation by π , and the reflections along X and Y axes. These four symmetry operations constitute a 2v  (6) and (7), the tangential stiffness matrices can be decomposed into four 2 2× block matrices, as listed in Table 1.

Irreducible Representation
Therefore, the eigenvalues and eigenvectors can be independently extracted from each of the block matrices. Figure 2 plots the eigenvalue decomposition results of the tangent stiffness matrix of the structure shown in Figure 1a. The endpoints of the polyline and the discrete dots are the solutions obtained by the symmetry method and the conventional method, respectively. As the minimum eigenvalue is positive, min ( ) 0 , this structure is stable. However, the minimum eigenvalue is much smaller than the others, as shown in Figure 2. This is because the structure has an internal and infinitesimal mechanism mode, and the first-order eigenvalue is mainly determined by the initial prestress [40,41]. Besides, the results obtained by the symmetry method are exactly the same as those of the conventional method. Because the tangent stiffness matrix has been diagonalized into four small-sized blocks, computation cost of the symmetry method is only 0.004 s. However, the conventional method concerns the original matrix and takes 0.0156 s. In Figure 2, the first, third, sixth, and eighth modes of the eigenvectors obtained from different symmetry spaces are also plotted, where A1 indicates full symmetry and the others indicate lowerorder symmetry. The short arrows indicate the motion trend of the free nodes. It turns out the As the minimum eigenvalue is positive, λ min (K T ) > 0, this structure is stable. However, the minimum eigenvalue is much smaller than the others, as shown in Figure 2. This is because the structure has an internal and infinitesimal mechanism mode, and the first-order eigenvalue is Symmetry 2018, 10, 229 7 of 16 mainly determined by the initial prestress [40,41]. Besides, the results obtained by the symmetry method are exactly the same as those of the conventional method. Because the tangent stiffness matrix has been diagonalized into four small-sized blocks, computation cost of the symmetry method is only 0.004 s. However, the conventional method concerns the original matrix and takes 0.0156 s.
In Figure 2, the first, third, sixth, and eighth modes of the eigenvectors obtained from different symmetry spaces are also plotted, where A 1 indicates full symmetry and the others indicate lower-order symmetry. The short arrows indicate the motion trend of the free nodes. It turns out the eigenvectors keep consistent with the inherent symmetry properties of specific symmetry subspaces. For instance, the fifth and eighth eigenvectors come from the full symmetry subspace and thus they have full symmetry of C 2v . Similarly, since the first and fourth eigenvectors are from the symmetry subspace A 2 , they exhibit rotational symmetry of C 2v . The third and seventh eigenvectors keeps C v symmetry along the X axis, while the second and the sixth eigenvectors keeps C v symmetry along the Y axis.
On the other hand, stability of the structure shown in Figure 1b is also studied, where the upper and lower horizontal cables are replaced by two cross cables. As shown in Figure 3, this structure has a single state of self-stress with full symmetry (i.e., C 2v symmetry) and an internal mechanism mode with lower-order symmetry (i.e., C v symmetry). However, it has been verified that the internal mechanism mode cannot be rigidified by prestressing (see Figure 3), as the minimum eigenvalue λ min (K T ) = −0.037 < 0 . Consequently, this structure is unstable, although it satisfies the necessary conditions given by Equation (16). It should be pointed out that the results for this structure are consistent with the reported ones [42]. On the other hand, stability of the structure shown in Figure 1b is also studied, where the upper and lower horizontal cables are replaced by two cross cables. As shown in Figure 3, this structure has a single state of self-stress with full symmetry (i.e., C2v symmetry) and an internal mechanism mode with lower-order symmetry (i.e., Cv symmetry). However, it has been verified that the internal mechanism mode cannot be rigidified by prestressing (see Figure 3), as the minimum eigenvalue . Consequently, this structure is unstable, although it satisfies the necessary conditions given by Equation (16). It should be pointed out that the results for this structure are consistent with the reported ones [42].

Form-Finding Analysis on Tensegrity Structures
A significance of novel prestressed structures mentioned in this study is that these structures rely on initial prestresses to obtain or enhance the structural stiffness. Therefore, form-finding (or force-finding) analysis on these structures is important [43,44]. For a structure with multiple selfstress states and complex geometry, Yuan et al. [45] considered the symmetry of the structure and proposed the concepts of integral self-stress state and feasible prestress to seek an effective and proper distribution for initial prestresses. They proposed the double singular value decomposition(DSVD) method for form-finding, which manually classified the members into different types and imposed symmetry constraints on the internal forces in the members. Zhang et al. [46] have utilized the energy method to obtain the integral self-stress state. Recently, some researchers [17,18,32] have combined the force density method with the structural symmetry to simplify the form-finding analysis for novel prestressed structures.
Admittedly, these methods rely on the correct classification for the members. When the configuration becomes complicated or the number of members rises [47], the involved form-finding becomes difficult. Notably, in a symmetry-adapted coordinate system, each block matrix is associated with a certain symmetry subspace, and thus the null space of the block with full symmetry necessarily contains independent self-stress states that satisfy the requirements of full symmetry [36,48]. In this section, the authors will describe how to adopt group theory for computing the integral self-stress states.

Integral Self-Stress State Obtained from the Block with Full Symmetry
According to the force equilibrium condition for each node of a structure, the equilibrium equation of the structure can be established:

Form-Finding Analysis on Tensegrity Structures
A significance of novel prestressed structures mentioned in this study is that these structures rely on initial prestresses to obtain or enhance the structural stiffness. Therefore, form-finding (or force-finding) analysis on these structures is important [43,44]. For a structure with multiple self-stress states and complex geometry, Yuan et al. [45] considered the symmetry of the structure and proposed the concepts of integral self-stress state and feasible prestress to seek an effective and proper distribution for initial prestresses. They proposed the double singular value decomposition(DSVD) method for form-finding, which manually classified the members into different types and imposed symmetry constraints on the internal forces in the members. Zhang et al. [46] have utilized the energy method to obtain the integral self-stress state. Recently, some researchers [17,18,32] have combined the force density method with the structural symmetry to simplify the form-finding analysis for novel prestressed structures.
Admittedly, these methods rely on the correct classification for the members. When the configuration becomes complicated or the number of members rises [47], the involved form-finding becomes difficult. Notably, in a symmetry-adapted coordinate system, each block matrix is associated with a certain symmetry subspace, and thus the null space of the block with full symmetry necessarily contains independent self-stress states that satisfy the requirements of full symmetry [36,48]. In this section, the authors will describe how to adopt group theory for computing the integral self-stress states.

Integral Self-Stress State Obtained from the Block with Full Symmetry
According to the force equilibrium condition for each node of a structure, the equilibrium equation of the structure can be established: where H is the equilibrium matrix, t is the internal force vector, and P is the external load vector. Because most novel prestressed structures are statically indeterminate, the self-stress states S come from the null space of the singular matrix H (i.e., HS = 0). In the symmetry-adapted coordinate system, the symmetry-adapted equilibrium matrix H can be transformed from the original matrix H: where V P is the transformation matrix associated with the load vector P, and V t is the transformation matrix associated with the internal force vector t; the positive integer i ∈ [1, µ], h ∈ [1, l i ], and l i and l µ are the dimensions of the irreducible representations Γ (i) and Γ (µ) . Equation (18) where V The integral self-stress state of the structure is given by [48] Note that the integral self-stress state in Equation (21) has considered the whole symmetry of the structure, which allows the members of the same type to retain equal prestress. Thus, the solution space for the form-finding problem of prestressed structures can be effectively reduced, especially for the structures with multiple self-stress states [48,49]. Moreover, this symmetry method does not need to manually classify the members into different groups in advance, and it avoids repeating calculations of conventional methods. Thus, the involved computation process is simple and convenient, and particularly suitable for structures with complex geometry [50]. Figure 4 shows the geometric configuration and symmetry rotations for a simple prismatic tensegrity structure, which consists of six pin-joints, six tension cables and three compression struts. The twist angle between the bottom and top triangles is π/6; the radius of the circumcircle of the triangles is unit length, and the vertical height of the structure is unit length. out that this symmetry method just computes the null space of the first block matrix and extract integral self-stress states from all the independent self-stress. Because the dimension of the first block matrix is significantly smaller than that of the original matrix, the computational complexity is effectively reduced. In comparison with the conventional DSVD method, this method avoids manual classification on the members and repeated computations on the matrices. As far as a highly symmetric structure with many self-stress states and members is concerned, the computational efficiency of this symmetry method is considerably improved [48].

Symmetry-Adapted Frequency Analysis
Generalized eigenvalue problems can be frequently found among structural analysis [16,40], which generally is a challenge from a computational perspective. For example, frequency analysis This symmetric structure has τ = 6 independent symmetry operations: the identity, two rotations around the vertical C 3 axis, and three two-fold rotations C 2 around the axes connecting the center of the structure and the mid-points of the diagonal cables (indicated by the dotted lines in Figure 4a). Thus, this structure belongs to the D 3 group [34,51]. Thereafter, the full symmetry subspaces corresponding to the internal force vector and the load vector are Because this 3 × 3 matrix H (1−1) is a rank-deficient (its rank r (1−1) = 2), it has a s (1−1) = 1-dimensional null space Therefore, the integral self-stress state of the structure is: As shown in Figure 4b, the obtained integral self-stress state satisfies the feasibility condition for the members, where six horizontal cables, three vertical cables, or three struts respectively belong to the same type. All the cables are in tension while the struts are in compression. It should be pointed out that this symmetry method just computes the null space of the first block matrix and extract integral self-stress states from all the independent self-stress. Because the dimension of the first block matrix is significantly smaller than that of the original matrix, the computational complexity is effectively reduced. In comparison with the conventional DSVD method, this method avoids manual classification on the members and repeated computations on the matrices. As far as a highly symmetric structure with many self-stress states and members is concerned, the computational efficiency of this symmetry method is considerably improved [48].

Symmetry-Adapted Frequency Analysis
Generalized eigenvalue problems can be frequently found among structural analysis [16,40], which generally is a challenge from a computational perspective. For example, frequency analysis can be taken as a well-known generalized eigenvalue problem in terms of the tangent stiffness matrix K T and the mass matrix which describes free vibrations of a general structure. In Equation (27), ω denotes a circular frequency and Ψ is the corresponding mode shape. On condition that this structure is symmetric, group-theoretic process considers smaller block matrices rather than the full matrices. Using the transformation matrix U in Equation (7), the mass matrix M can be decomposed into similar block-diagonalized forms

Illustrative Example: A 12v C Symmetric Cable Dome Structure
The well-known Levy cable dome is a typical prestressed cable-strut structure [40,53]. Figure 5 shows a highly symmetric Levy cable dome with a diameter of 100 m. This structure is composed of 84 pin-joints, 36 compression struts and 168 tension cables, where the boundary nodes 1-12 are constrained in three directions. Note that cross-sectional areas and initial prestresses of different types of members of this cable dome structure have been reported [40]. In the view of symmetry, this cable dome structure keeps unshifted by twelve rotations along the 12-fold rotation axis 12 C , and twelve reflections 1 1 2 σ σ − , as shown in Figure 5a. Thus, this structure is 12v C symmetric. According to the point-group theory table [28], the 12v C symmetry group has four one-dimensional ( [1,4]  In the view of symmetry, this cable dome structure keeps unshifted by twelve rotations along the 12-fold rotation axis C 12 , and twelve reflections σ 1 − σ 12 , as shown in Figure 5a. Thus, this structure is C 12v symmetric. According to the point-group theory table [28], the C 12v symmetry group has four one-dimensional (l i∈ [1,4] = 1) irreducible representations A 1 , A 2 , B 1 and B 2 , and five two-dimensional (l i∈ [5,9] = 2) irreducible representations E 1 − E 5 . To establish the symmetry-adapted coordinate system, the 216 × 216 orthogonal transformation matrix U for the cable dome is obtained from Equation (6). This matrix is composed by independent vectors associated with fourteen symmetry subspaces, where the first vector U (A 1 ) for the full symmetry subspace is a 216 × 12 matrix, U (A 2 ) is a 216 × 6 matrix, U (B 1 ) is a 216 × 10 matrix, U (B 2 ) is a matrix, and all the vectors U (E 11 ) -U (E 52 ) for the lower-order symmetry subspaces are 216 × 18 matrices. These 14 independent symmetry subspaces are utilized to decompose the 216 × 216 stiffness matrix or the mass matrix into 14 block matrices with smaller dimensions. For example, the distribution patterns of nonzero entries of the original tangent stiffness matrix and the block-diagonalized matrix are shown in Figure 6. As shown in Figure 6b, all the block matrices are independent and distribute along the main diagonal. two-dimensional ( [5,9]

U
for the lower-order symmetry subspaces are 216 ×18 matrices. These 14 independent symmetry subspaces are utilized to decompose the 216 216 × stiffness matrix or the mass matrix into 14 block matrices with smaller dimensions. For example, the distribution patterns of nonzero entries of the original tangent stiffness matrix and the block-diagonalized matrix are shown in Figure 6. As shown in Figure 6b, all the block matrices are independent and distribute along the main diagonal. As far as the frequency analysis is concerned, the original problem is neatly simplified into fourteen sub-problems involved with small-sized matrices. Figure 7 shows the first 100 frequencies obtained by the proposed symmetry method, which is compared with the results by the conventional numerical method. In addition, to study the influence of initial prestresses on the natural frequencies, the C12v symmetric cable dome with different prestress levels are analyzed by the proposed method (the symmetry subspaces keep invariant). The prestress levels are respectively 0, 0.25t, 0.5t, 2t and 4t, and the initial prestresses t are determined by the feasible prestress modes [40]. As far as the frequency analysis is concerned, the original problem is neatly simplified into fourteen sub-problems involved with small-sized matrices. Figure 7 shows the first 100 frequencies obtained by the proposed symmetry method, which is compared with the results by the conventional numerical method. In addition, to study the influence of initial prestresses on the natural frequencies, the C 12v symmetric cable dome with different prestress levels are analyzed by the proposed method (the symmetry subspaces keep invariant). The prestress levels are respectively 0, 0.25t, 0.5t, 2t and 4t, and the initial prestresses t are determined by the feasible prestress modes [40]. As far as the frequency analysis is concerned, the original problem is neatly simplified into fourteen sub-problems involved with small-sized matrices. Figure 7 shows the first 100 frequencies obtained by the proposed symmetry method, which is compared with the results by the conventional numerical method. In addition, to study the influence of initial prestresses on the natural frequencies, the C12v symmetric cable dome with different prestress levels are analyzed by the proposed method (the symmetry subspaces keep invariant). The prestress levels are respectively 0, 0.25t, 0.5t, 2t and 4t, and the initial prestresses t are determined by the feasible prestress modes [40].  It can be noticed from Figure 7 that the obtained results are in good agreement with those of the numerical method, where the maximum error is 10 −9 . Lower-order natural frequencies are low and intensively distributed. Notably, the first 14 frequencies are dominated by the prestress level, as they are associated with m = 14 infinitesimal mechanism modes. In fact, the initial prestresses also improve the structural stiffness of the cable dome. After increasing the initial prestress level, the first 50 frequencies of this structure increase significantly.
On the other hand, Figure 8 depicts the first eight vibration shapes of the C 12v symmetric cable dome [40]. It can be noticed that this symmetric dome structure has many repeated eigenvalues and equivalent eigenvectors, such as modes 1-2, modes 3-4, and modes 5-6. In fact, this phenomenon is ubiquitous for most symmetric structures. This is because the roots computed from the symmetry spaces for multi-dimensional irreducible representations are identical, and the generalized eigenvalues for these symmetry subspaces are exactly the same. In addition, each vibration shape obtained from lower-order symmetry subspace U (i−h) does not maintain its full symmetry (i.e., C 12v ). Then, it may be reduced to a low-order symmetry. Because the seventh vibration shape is obtained from the symmetry subspace Γ (2−1) = A 2 , it has rotational symmetry (i.e., C 12 ). In other words, the symmetry of all the vibration modes can be predicted from symmetry subspaces in advance [16,29], without numerical computing. It can be noticed from Figure 7 that the obtained results are in good agreement with those of the numerical method, where the maximum error is 9 10 -. Lower-order natural frequencies are low and intensively distributed. Notably, the first 14 frequencies are dominated by the prestress level, as they are associated with 14 m = infinitesimal mechanism modes. In fact, the initial prestresses also improve the structural stiffness of the cable dome. After increasing the initial prestress level, the first 50 frequencies of this structure increase significantly.
On the other hand, Figure 8 depicts the first eight vibration shapes of the C12v symmetric cable dome [40]. It can be noticed that this symmetric dome structure has many repeated eigenvalues and equivalent eigenvectors, such as modes 1-2, modes 3-4, and modes 5-6. In fact, this phenomenon is ubiquitous for most symmetric structures. This is because the roots computed from the symmetry spaces for multi-dimensional irreducible representations are identical, and the generalized eigenvalues for these symmetry subspaces are exactly the same. In addition, each vibration shape obtained from lower-order symmetry subspace ( ) does not maintain its full symmetry (i.e., 12v C ).
Then, it may be reduced to a low-order symmetry. Because the seventh vibration shape is obtained from the symmetry subspace (2 1) 2 A − Γ = , it has rotational symmetry (i.e., 12 C ). In other words, the symmetry of all the vibration modes can be predicted from symmetry subspaces in advance [16,29], without numerical computing.

Conclusions
This paper has described some new developments of symmetry analysis on novel prestressed structures using group theory. Through theoretical group-theoretic approaches and illustrative examples, basic process and key features of group-theoretic approaches in stability analysis, formfinding analysis, frequency analysis and buckling analysis have been identified. A significant advantage of these group-theoretic approaches is that the involved large-sized matrices can be neatly

Conclusions
This paper has described some new developments of symmetry analysis on novel prestressed structures using group theory. Through theoretical group-theoretic approaches and illustrative examples, basic process and key features of group-theoretic approaches in stability analysis, form-finding analysis, frequency analysis and buckling analysis have been identified. A significant advantage of these group-theoretic approaches is that the involved large-sized matrices can be neatly block-diagonalized into many smaller-sized and independent matrices. Thus, the computational complexity for dealing with the original problem will be significantly reduced. More importantly, these methods can qualitatively reveal certain physical meanings and obtain effective insights into the involved problems. It is worth mentioning that group-theoretic approaches are systematic and general. When a structure has only the lowest-order symmetry, the method still works and degenerates into a conventional method. This research is helpful for enriching developments in the field of symmetric structures.