From the viewpoint of calculation cost, truss element and membrane element, which do not consider the bending of an element, are often used for analyses of large-scale dynamic FEM (Shirasawa et al., 2011;Ono et al., 2014).However, there are few studies that efficiently detect buckling from the entire structure during response analysis using the element model which does not consider bending. Because solution of buckling analysis is generally acquired from equation of bending moment. Therefore, the authors have proposed and developed a method to detect wrinkle and slack using the membrane elements (Arita et al., 2014) and a method to detect snap through buckling during the dynamic analysis using the truss element (Arita and Miyazaki, 2018). This paper proposes a comprehensive and efficient visualization method of buckling occurrence region and buckling magnitude during dynamic response analysis by using these two theories to progress convenience in actual design of gossamer spacecraft. By visualizing the buckling region, it is expected at actual designing that it will be easier to specify the part where buckling easily occurs, and to control or deal with the buckling part. A general method of visualization of buckling is to display the buckling mode shape in order to examining the buckling load and the shape after buckling (Noguchi and Fujii, 2000;Ikeda dt al., 2003). On the other hand, the purpose of the method proposed in this paper is not to obtain the mode shape but to grasp comprehensively when and what area the buckling displacement occurs during transient response of the deployment. Visualization based on such a concept is a unique point of this research. In a previous study, a method using tension field theory has been proposed as a method of calculating the amount of deformation after buckling (Conci, 2007).However, it is pointed that the assumption of tension field theory does not establish for the buckling of flexible structures (Iwasa et al., 2004).Therefore, the method proposed in this paper is to calculate the buckling displacement amount without using the tension field theory. The theory of the visualization method is shown in Section 2. Section 3 shows the results of verification of the proposed method by a simple truss model, and Section 4 shows an analysis example assuming development of flexible membrane. # II. Theories of Buckling Analysis and Visualization In the proposed method, the dynamic response of the entire structure is analyzed by the truss element, and the amount of the snap through buckling displacement of the node is visualized by color contour. However, since the truss element cannot express the deformation other than in the axial direction, it is necessary to calculate the long column buckling of the element itself in order to grasp the comprehensive buckling of the entire structure. Since a flexible member such as a membrane has a characteristic that compression rigidity is very small compared with tension and it is easy to buckle, it is necessary to devise modeling of the compressive stiffness in post buckling analysis. Hence, we adopt the theory called Mod-SRM, which models appropriately such features, to detect and visualize the buckling of each truss element. In other words, this paper proposes the method to visualize the comprehensive buckling of the entire structure by displaying the color contour of the snap through buckling of nodes and the long column buckling of elements with respect to the dynamic transient response diagram of the truss element. In this chapter, the detection and visualization method of the snap through buckling of nodes is explained in Section 2.1, and the detection and visualization method of the long column buckling of elements by Mod-SRM is explained in Section 2.2. Only the snap through buckling of the node is detectable in the truss element. Since detection and quantification of the snap through buckling are described in detail in the previous work (Arita and Miyazaki, 2018) by the authors, in this section, a schematic view of the method is shown in Fig. 1 Q from each eigen mode of the stiffness matrix K , the eigen mode is judged if buckling or not according to the work of deformation. Multiple buckling modes often appears at the same time. In such a case, quantification is carried out in order to determine the most likely buckling mode considering the motion state of the structure. Assuming that buckling displacement x is caused by external disturbance force F , F is obtained by solving the equation of motion. Then, the norm of the obtained F is defined as a validation parameter "DF value". It is determined that the mode with the smallest DF value is the most likely buckling mode. Moreover, the norm of the buckling displacement x is also defined as another validation parameter "BD value". # ii. Visualization of Buckling of Nodes Based on Buckling Displacement In this section, the authors propose a method of visualization using the detected buckling mode and the BD value explained above. The visualization method is to display the amount of the snap through buckling displacement of each node is visualized by color contour for the most likely buckling mode. The element used in this theory is a two-node truss element having three degrees of freedom of x, y, z per node. We define n as the number of nodes of the entire structure, and define h as the most likely buckling mode. h is a normalized vector, in which the proportion of x , y and z displacements of each node is arranged. Then, the component corresponding to ( , , ) x y z of the j th node in h is defined as ( , , ) j j j x y z , and j h is defined as follows. [ , , ] T j j j j x y z º h (1) buckling by calculating the work for each node. Nodes of which deformation work becomes 0 or less are considered as buckling nodes because the displacement proceed by negative deformation work. To calculate the deformation work, it is necessary to extract only the deformation component from the mode that includes rigid-body motion. The orthogonal component * h to the rigid-body modal space , which means the deformation component, is extracted by Eq. ( 2). The deformation work of each node * j w in the buckling mode j h is calculated by using the stiffness matrix K as Eq. ( 3). ( )2 ( Buckling displacement amount (BD value) is calculated for the node of which * j w is 0 or less by Eq. ( 4). ? is a scalar giving the magnitude of the displacement in the buckling mode direction, and the details are explained in the previous study (Arita and Miyazaki, 2016). Therefore, only the results are shown here. (4) (5) where (6) b is the coefficient of the Newmark-? method, t D is the time step width in the time integration of the Newmark-? method, M is the mass matrix, and are the velocity and acceleration of the previous time step respectively. Note that although Newmark-? method was used in the previous paper (Arita and Miyazaki, 2016), the concept of this method is also effective for other numerical integration methods, and similar equations can be derived by other methods. The magnitude of the displacement by the buckling is visualized by color contour display of nodes according to the BD value. # b) Long Column Buckling of Elements in Truss Elements i. Introduction of Mod-SRM A flexible member is easy to buckle because the compressive stiffness of the flexible member is very small compared to tensile stiffness. Therefore, methods that determine the compressive stiffness by multiplying the tensile stiffness by a coefficient smaller than 1 has been proposed for the element model that does not consider bending (e.g. Miyazaki, 2006). Mod-SRM, which is one of them, proposed to be able to determine the compressive stiffness ratio uniquely according to the amount of the out-of-plane buckling by introducing the stretchable elastic theory into the element model (Arita et al., 2014). Because Mod-SRM is based on the bending deformation of a beam, the visualization method proposed this paper is to apply Mod-SRM to the truss element. The outline of Mod-SRM and the derivation of physical quantities for the visualization are introduced in this section. First of all, supposing one truss element of which total length is defined / 2 l deforms like a slack of a cable, the equilibrium equation is expressed as Eq.( 7) in an infinitesimal line element dx subjected to the load P as shown in Fig. 2.Eventually, Eq.( 8) and ( 9)are obtained from the equilibrium equation. Details of this derivation are written in the previous work (Arita et al., 2014). Fig. 2: Mathematical modeling of Mod-SRM. Supposing one truss element deforms like a slack of a cable, the equilibrium of forces and moments are obtained in an infinitesimal line element dx subjected to the load P . P P 0 q q dx 0 q q 0 M = ds z x y Before deformation After deformation / 2 l d + M M M e P P s x Q * 1 () () n k k kk k hq h h q qq ( ) * * * jj j w h Kh =? j BD h = 2 a b a c a 1 2 1 0 1 0 1 ( [ ] ) 1 ( [ ] ) 1 ( [ ] ) 2 TT TT TT t t a I Q Q Q Q Mh b I Q Q Q Q Mx c I Q Q Q Q Mx 0 x 0 x © 2020 is the equivalent in-plane compressive strain, defined as shown in Fig. 3. * E means Young's modulus after buckling, and denotes the ratio of stiffness in the compressive direction to the tensile direction. l is a material constant and is defined by Eq. (10). is a non dimensional compressive load and is defined by Eq. ( 11). E means Young's modulus, A means crosssectional area and I means moment of inertia of area. C q and ( ) f t are functions defined by Eq. ( 12) and ( 13) respectively. The compressive stiffness ratio is obtained by Eq. ( 14 E is Young's modulus after buckling, and is the compressive stiffness ratio. That is to say, given the equivalent plane compressive strain and the material constant l , the relation between the non-dimensional compressive load and strain can be obtained when the load and angle 0 are determined under Eq.( 8) and ( 9).In the actual calculation, 0 is given and h is obtained by Newton's method so as to satisfy Eq. ( 8), and is obtained by Eq. ( 15): (15) It is convenient to calculate the approximation of as the polynomial of by using the relation obtained by Eq. ( 15) in advance, and in the transient response analysis, is calculated from the polynomial according to the value of at the time step. Hence, the compressive stiffness ratio is decided every time step by Eq. ( 14). # ii. Visualization of Buckling of Elements Based on Buckling Displacement In visualization, the buckling of each element is judged and the buckling elements are displayed in color contour according to the magnitude of . The calculated value of compressive buckling load is used for the judgement of the buckling. In general, when an axial compressive load P is applied to a simply supported beam at both ends of bending stiffness EI and length l shown in Fig. 4, the Euler buckling load 0 cr P is obtained as Eq. ( 17) from Eq. ( 16), which is the equilibrium equation assuming no expansion, contraction and shear deformation. On the other hand, considering expansion and contraction and shear deformation, if the equilibrium equation can be written as Eq. ( 18), the buckling load cr P is written as Eq.(19). The integral term in the Eq. ( 21) can be written as below; (22) Here, the following variable transformation is defined; The following equation is obtained from Eq. (23); 0 1 ( ) d d ( ) cos C u u C u ? ? + = ?(24) Using Eq. ( 21), ( 22) and (23), the following relation is obtained: (25) Furthermore, using Eq. ( 23) and (24), the ( ) h u is written as below; (26) Because 0 0 ? ? in the case of buckling, the values of ( ) C u , ? and ( ) h u are as follows; (cos cos ) ( ) (cos cos ) '( ) cos f d f f f d f f f d f f f d f f d ?? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? = ? ? ? ? = ? ? ? + ? ? + = ? ? ? ?? ? ? ? 00 00 1 cos 1 cos cos ( ) cos 22 : : C u u u ?? ? ? ? ? ?? +? = = + ?? ? ? ?? ? ? ? 0 0 0 d d d d d d 21 (1 ( )) ( ) ( ( ) cos ) '( ) ( ) cos 2 (1 ( )) ( ) ( )uu ss P C u f C u f d EI C u P C u f h u EI ? ? ? ? ? ? ? ? ? ? ? = ?? = + ? ? ?? ? ?? ? + ? ? 0 0 0 0 0 0 0 3 2 0 1 ( ) ( ( ) cos ) '( ) ( ) cos 2 (1 cos )(1 cos ) '( ) (1 cos )(1 cos ) 2 ( ) cos 1 (1 cos ) '( ( )) 1 cos 1 ( ) 1 cos (1 cos ) '( ( )) 2 1 cos 1 ( ) u u h u C u f d Cu u fd u C fd u C f d u C ? ? ? ? ? ? ? ? ? ? ? ?? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? =? ? ?+ = ?+ ? =+ + + ?+ = + + ? ? ? ? © 2020 ( ) 1 , 0 , ( ) 0 C u h u ? ? ? ? (27) Thus, Eq. ( 19) is derived by Eq. (25) as follows; (28) Q.E.D. In the case of Mod-SRM, the equilibrium equation is Eq. ( 7). Therefore, cr P is below; (29 ) where cr P is smaller of the two solutions. From Eq. ( 11), is below; (30) The same solution is obtained by setting in Eq.( 8) and ( 13).When 0 0 = , Thus, applying it to Eq.( 8), Solving this for , Two solutions are obtained. Since the buckling load is the one with a smaller value, cr is obtained as Eq. (30). Incidentally, the buckling strain is obtained from Eq. (15) as follows. (35) # III. Verification of Visualization of Buckling by Numerical Experiment The authors verified the visualization method with a simple truss model consisting of 3 elements and 4 nodes. The truss model is shown in Fig. 5, and its specifications are shown in Table 1. Before transient response and visualization, the relationship and polynomial approximation between and were obtained in advance by Mod-SRM as shown in Fig. 6.The tensile strain is the positive value, and the compressive strain is the negative value. The buckling load of an element was obtained as 1.000129 = that is, 0.005623[N] cr P = and 0.0001285 cr = . Note that the compressive stiffness ratio a is calculated in each time step of the transient response analysis based on the polynomial approximation and Eq. ( 14). The relationship between and , incidentally, is shown in Fig. 7. The result of the buckling analysis is shown in Fig. 8. The horizontal axis is the time step. The first vertical axis at the left side is the strain of the element. The second vertical axis at the right side represents the detection result of the snap through buckling of the nodes. Since the nodes ? and ? are fixed, they are excluded. Figure 8 shows the sharp fluctuation of the In the previous study, Mod-SRM only defined the compressive stiffness ratio for the post buckling analysis of a membrane. The study of the dynamic buckling of a truss element only quantified the ease of buckling and degree of buckling deformation as representative values for the entire structure. On the other hand, the method proposed in this paper enables visually recognizing where and how large the buckling occurs in a transient analysis by calculating the magnitude of buckling displacement for each node and element. # IV. Application Example of Membrane Structure The authors conducted a calculation as an application example with a model assuming that the triangular membrane deploys with booms. The triangular membrane folded into a bellows deploys as shown in Fig. 12 and Table 2. The transient response is shown in Fig. 13. Buckling occurring at the nodes and the elements are visualized, indicating that we can confirm when, where and how large the buckling occurs in the structure. [s] Young's ratio E 3.5×10 # V. Conclusions The authors proposed a comprehensive and efficient visualization method of buckling region during dynamic response analysis based on two previous studies to progress convenience in actual design of gossamer spacecraft. The proposed method is verified by a simple truss model, and the calculation of an application example is conducted. The result indicates that using this method, we can confirm when, where and how large the buckling occurs in the structure. 1![Fig.1: Buckling detection method. After extracting the orthogonal component to the rigid-body modal space Q from each eigen mode of the stiffness matrix K , the eigen mode is judged if buckling or not according to the work of deformation.](image-2.png "Fig. 1 :") 3![Definition of symbols. is the equivalent in-plane compressive strain, *](image-3.png "3 :") 4![Fig. 4: Model of beam buckling2](image-4.png "Fig. 4 :") ![of Buckling Region in Dynamic Transient Analysis of Cable Network Structures strain after Step 750. It can be considered that the drastic decrease of element stiffness due to compressive buckling causes the sharp fluctuation of the strain. Figure8also shows that snap through buckling of the nodes occur from Step 794 onwards. Figure9is the enlarged figure of Fig.8from Step 794 toStep 805 regarding the strain of the elements. Figure10is also the enlarged figure of Fig.8regarding the nodes, however, the vertical axis shows not the detection result but the BD value of the buckling nodes. Transient response is shown in Fig.11as a result of visualization of the buckling. Comparing Fig.9, 10 and 11, it can be confirmed that the color contour is properly visualized in the response diagram corresponding to the strain of the elements and the BD value of the nodes where the buckling occurred.](image-5.png "") 5![Fig. 5: Truss model for verification.](image-6.png "Fig. 5 :") 67![Fig. 6: Relation and polynomial approximation between and . at each time step in transient response is determined by the polynomial according to .](image-7.png "Fig. 6 :Fig. 7 :") 8![Fig.8: Result of the buckling analysis. The first vertical axis at the left side is the strain of the element. The second vertical axis at the right side represents the detection result of the snap through buckling of the nodes. The drastic decrease of elements tiff ness due to compressive buckling causes the sharp fluctuation of the strain after Step 750. The snap through buckling of the nodes occur from Step 794 onwards.](image-8.png "Fig. 8 :") 9![Fig. 9: The strain of the elements from Step 794 to Step 805. Color contour of each element in Fig. 10 is based on the magnitude of the strain.](image-9.png "Fig. 9 :") 10![Fig. 10: BD value of the buckling nodes from Step 794 to Step 805.Color contour of each node in Fig. 10 is based on the magnitude of the BD value.](image-10.png "Fig. 10 :") 12![Fig. 12: Membrane structure model. The red node is fixed and the blue nodes are forcibly displaced every time step. At first, the membrane is folded into a bellows.](image-11.png "Fig. 12 :") ![node and forcibly displace the two blue nodes by 0.75×10 -7 per time step in the +z direction. © 2020 Global Journals lobal Journal of Researches in Engineering ( ) Volume Xx X Issue I Ve rsion I 19 Year 2020 J Gl Efficient Visualization Method of Buckling Region in Dynamic Transient Analysis of Cable Network Structures](image-12.png "") 11![Fig. 11: Transient response and the result of the visualization. Comparing with Fig.8 and 9, it can be confirmed that the color contour is properly visualized corresponding to the strain of the elements and the BD value of the nodes where the buckling occurred.](image-13.png "Fig. 11 :") 13![Fig. 13: Transient response and the result of the visualization. Snap through buckling occurs everywhere on the membrane, and the buckling of the element occurs intensively at the corners at the end of the deployment.](image-14.png "Fig. 13 :") 1 2ParameterSymbol ValueUnittime step sizedt5×10-6 * A study of dynamic evaluation of structural buckling SArita YMiyazaki Mechanical Engineering Letters 2 2016 * Modified method of detecting dynamic buckling SArita YMiyazaki Mechanical Engineering Letters 4 2018 * Numerical model for prediction of wrinkling behavior on a thinmembrane structure SArita TOkumiya YMiyazaki Mechanical Engineering Journal 1 4 2014 * after buckling behavior of thin membranes and textiles based on tension field theory and experiments, Proceedings of Mechanics of Solids in Brazil AConci 2007. 2007 * Simulation and interpretation of diffuse mode bifurcation of elastoplastic solids KIkeda YYamakawa STsutsumi Journal of the Mechanics and Physics of Solids 51 2003 * Evaluation of Tension Field Theory for Wrinkling Analysis With Respect to the Post-Buckling Study TIwasa MCNatori KHiguchi Journal of Applied Mechanics 71 2004 * Wrinkle/slack model and finite element dynamics of membrane YMiyazaki International Journal for Numerical Methods in Engineering 66 7 2006 * Extracting the buckling mode from LDL T -decomposed stiffness matrix HNoguchi FFujii Advances in Computational Engineering and Sciences Satya N. Atluri and Frederick W. Brust Eds. 1 2000 Tech Science Press * Study of Asynchronous Solar Sail Deployment Using Finite Element Method, Transaction of the GOno KShintaku YShirasawa OMori YMiyazaki SMatsunaga Aerospace Technology Japan 12 2014 Japan Society for Aeronautical and Space Sciences * Analysis of Membrane Dynamics using Multi-Particle Method for Solar Sail Demonstrator "IKAROS YShirasawa OMori YMiyazaki HSakamoto MHasome NOkuizumi HSawada HFuruya SMatunaga MNatori Proceedings of 52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference AIAA 2011-1890 (2011