高等结构动力学大作业

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业Advanced Structural Dynamics ProjectThe dynamic response and stability analysis of the beam under vertical excitationInstructor：Dr. Li WeiName：Student ID：Problem description and the purpose of the project1.1 calculation modelAn Eular bea

2、m subjected to an axial force. Please build the differential equation of motion and use a proper difference method to solve this differential equation. Study the dynamic stability of the beam related to the frequency and amplitude of the force. As shown in the Fig 1.1.Fig1.11.2 purpose and process a

3、rrangementlearning how to create mathematical model of the continuous system and select proper calculation method to solve it.learning how to build beam vibration equation and solve Mathieu equation.using Floquet theory to judge vibration systems stability and analyze the relationship among the freq

4、uency and amplitude of the force and dynamic response.This project will introduce the establishment of the mathematical model of the continuous system in section 2, the movement equation and the numerical solution of using MATLAB in section 3, Applying Floquent theory to study the dynamic stability

5、of the beam related to the frequency and amplitude of the force in section 4. In the last of the project, we get some conclusions in section 5.the mathematical model of the systemThe geometric model of the beam and force-simplified diagram is shown in Fig.2.1.We assume that its stiffness(EI) is cons

6、tant and the deflection of the beam is small, and the boundary conditions is simply support. Now the beam subjected to an axial force. We assume the force is equal to. Fig.2.1We select the length of in any position of the beam, the free-body diagram is shown in Fig.2.2.Fig.2.2Using equations of move

7、ment equilibrium, that is to say: + （1） （2）From equation (1), we will get: （3）Divide equation (3) by and take the limit: （4）Then synthesize equation (2),we can get: （5）Divide equation (5) by and take the limit: （6）Combine equation (4) with (6): （7）And （8）Combine equation (7) with (8): （9）We know EI

8、is a constant, so （10）In equation (10), m is the mass of unit length.Now we will use assumed-modes method. Named ,so: （11） n=1,2,. （12）In the equation (12)And ,so n=1,2,. （13） （14）In the equation (12) Equation (14) is the Mathieu equation. it is difficult to solve the analytical solution directly, t

9、hus, we use the approximate derivative namely an average acceleration method to get the numerical solution from the reference.3. Numerical solution3.1 using MATLAB to solve equationWe will use the Newmark- method 1 to solve equation (14). We can use the initial condition to integrate the move equati

10、on: （15）Fig.3.1As shown in Fig.3.1 （16） （17） （18）From equation (16), (17) and (18), we will get: （19） （20） （21）When applying the MATLAB, we need discrete the processing time t, get time step.When solving the vibration stability interval, there are three variables to participate in the discussion, na

11、mely. So take a particular w first and discuss the remaining two parameters.From Floquent theory2，we can use parameter A to judge stability.Equation （22）Take two sets of special solution: （20）Parameter2 （21）If abs (A) is less than 1, the system is stability. And if abs (A) is greater than 1, the sys

12、tem is instability. When abs(A) is equal to 1, the system is critical state.We use MATLAB Codes to solve equations. We use =2 Mathieu Equation to judge the validity of the codes. From Fig.3.2 and Fig.3.3, we can consider the codes are correct.In these follow figures, =2, the horizontal axis is , ver

13、tical axis is .Fig.3.2. stable domain in reference3and2Fig.3.3. stable domain in MATLAB solutionCompared Fig.3.2 with Fig.3.3, we can see that the stability domain of numerical solutions applying average acceleration method are consistent with the standard solutions. it can concluded that when the s

14、ystem have solution whose cycle is equal to or 2, This chapter discusses the accuracy of the vibration stability determination with Floquent theory. The next chapter will discuss the numerical solution and stable domain and two parameters influences on the stability for this question.4. Parameters i

15、nfluenceIn this part, we only consider two parameters, namely the frequency and amplitude of the force.4.1the influence of the forces frequency4.1.1the stability of the systemWhen we discuss the stability of the system related to the frequency of the force, we should select some different frequencie

16、s, so we choose =1,2,4,6,8 and10. Using MATLAB codes, we can obtain the figs of the stability. We can know the stable region is bigger with the increase of the frequency in Fig.4.1.=1 =2=4 =6=8 =10Fig.4.1 stable domain with different 4.1.2the response of the systemWhen we discuss the response of the

17、 system, the system should be stable. So we choose ，=2,4 and 6. In Fig.4.2, the cycle of the response increase and the range of the reactive amplitude is smaller with the increase of the frequency.Fig.4.2 responses of the system with different 4.2the influence of the forces amplitudeThe is related t

18、o the forces amplitude P. The cycle of the response a little increase and the range of the reactive amplitude is bigger with the increase of the forces amplitude, in Fig.4.3.and Fig4.4.Fig.4.3 vibration response curve with different The red curve is w=2, =12, =1; the blue curve is w=2, =14, =1.This

19、figure state that the vibration cycle is smaller and the amplitude have a little change with the increase of the .Fig.4.3 vibration response curve with different The red curve is w=2, =12, =1; the blue curve is w=2, =12, =5. This figure state that the amplitude is smaller and the vibration cycle hav

20、e a little change with the increase of the .5. Conclusion(1) With the increase of the frequency, the stable region and the cycle of the response are bigger, but the range of the reactive amplitude is smaller.(2)With the increase of the forces amplitude, the cycle of the response a little increase and the range of the reactive amplitude is much bigger.(3)Vibration in the stable reg