Disneyland ride final design report

Structural Dynamics Research Corporation (SDRC)

Daniel Belongia
Adam Mayer
Donny Kuettel
Nasser M. Abbasi

EMA 542 Advanced dynamics
University Of Wisconsin, Madison
Fall 2013

Abstract

Dynamic analysis was completed for a new spinning ride as requested by Walt Disney Corporation. Detailed derivation of model was completed for the main structural elements using rigid body dynamics.
Critical section was identified and maximum stress calculated to insure that the member does not fail during operations and passengers acceleration does not exceed 6  g.
Large software simulation program was completed to verify the model used and to allow selection of optimal design parameters.
Prepared by:
Dynamic design team
Structural Dynamics Research Corporation (SDRC)

Contents

1 Introduction
 1.1 Gantt chart and history of design project
2 Safety considerations
 2.1 Locations of possible failure in the structure
3 Mathematical model of system dynamics
 3.1 Review of the model structure used in the design
 3.2 Setting up the mathematical model
  3.2.1 Summary of design input and design output
  3.2.2 System dynamic loads and free body diagram
 3.3 Beam to column analysis
  3.3.1 Finding Mbeam  (beam dynamic moment)
  3.3.2 Finding Mcol  (column dynamic moment)
  3.3.3 Finding Mcabinet  (cabinet dynamic moment)
  3.3.4 Finding F
  cabinet  (cabinet dynamic linear force)
  3.3.5 Finding F
  beam  (beam dynamic translational force)
  3.3.6 Using free body diagram and solving for constraint forces
 3.4 Column dynamic analysis
4 Simulation of the dynamic equations found
 4.1 Review of the simulation
 4.2 Simulation output, time histories and discussion of results
 4.3 Discussion and analysis of results
 4.4 Cost analysis
5 Conclusions of results and future work
 5.1 Future work and possible design improvement
6 Appendix
 6.1 Use of simulator to validate different design parameters
 6.2 References

List of Figures

1 Artist rendering of ride after construction
2 Illustrating typical dynamic movement over four time instances for one revolution
3 Time line used by the design team in the development of the final project report
4 Identification of critical sections in the structure
5 Main parts of the ride structure
6 Main dimensions of ride structure
7 relation between rotating coordinates system, body fixed coordinates system, and body principal axes.
8 Beam dynamics. Balancing dynamic forces to external forces and reactions
9 Column dynamics. Balance with and external loads and beam transferred loads.
10 Configuration used for finding torque and force at beam/column joint
11 Rotating coordinates system xyz  used to find passenger acceleration
12 View of passenger head in the rotating coordinates system xyz
13 Rotating coordinates system xyz  used to find beam center of mass acceleration
14 Finding the bending moment at different locations along the span of the beam
15 Dynamic load balance between column and external loads
16 overview of simulator user interface
17 dynamic loads at the end of ride using optimal design values
18 critical section current and maximum moments and stresses
19 optimal set of parameters obtained from simulation.
20 simulator keeps track of maximum g  felt by passenger to insure it does not exceed 6g
21 acceleration and velocity of passenger time history and angular velocity time history of beam and column
22 time histories using the ramp down option used after reaching 6g goal
23 Changing the structure dimensions to select optimal design using simulation

List of Tables

1 design input parameters
2 design output
3 design output for loading and forces using optimal parameters found
4 cost estimate
5 ride statistics based on optimal design parameters

1 Introduction

A four-member team at Structural Dynamics Research Corporation (SDRC) has completed the final design for a new spinning ride for Disneyland.

The ride features two non-collinear components of angular velocity. The head of each passenger will experience a maximum of 6g  acceleration. Just before this acceleration is reached, the ride will enter steady state. During steady state, passengers will experience a small periodic fluctuation of acceleration that ranges between 4.8g  and 6g  but will not exceed 6g  . The ride can then enter the ramp down phase and starts to decelerate until it stops with smooth landing. All three phases of the ride have been simulated to insure the passengers will not exceed 6g  during any of the phases. The ride is specifically designed to be light, safe, affordable, and fun. The following is an artist rendering showing loading the passengers in the cabinet before starting the ride


PIC
Figure 1: Artist rendering of ride after construction


Once the cabinet has reached the top of the support column, the ride will start. Extensive simulation of the mathematical model of the dynamics of the model was performed to achieve an optimal set of design parameters in order to meet the design goals as specified in the customer requirements of a minimum weight and cost and at the same time insuring the structural members do not fail and that the passengers will safely achieve the 6g  acceleration in reasonable amount of time. The conclusion section outlines the final design parameters found. The following diagram illustrates typical one revolution ride for illustrations that was generated by the simulator developed specifically for this design contract


PIC
Figure 2: Illustrating typical dynamic movement over four time instances for one revolution


1.1 Gantt chart and history of design project

The design team followed the following timeline in the development of the report and the design. This is illustrated below using Gantt chart


PIC
Figure 3: Time line used by the design team in the development of the final project report


2 Safety considerations

The flight simulator will be equipped with multiple safety measures to ensure that passengers will have a fun and exciting ride. In order to ride the flight simulator, each passenger must be at least 5  feet tall. This insures that the riders can be securely fastened into the seat. Assuming an average rider weight of 175  pounds, one single rider cannot weigh more than 350  pounds.

Any more weight will induce a moment on the main arm that might be considered unsafe. A factor of safety was factored into the building of the arm in case two riders combined weight to be more than 350  pounds. This additional weight accounts for the seating weight and the frame of the cabinet as well.

While the ride is in motion, each passenger will be harnessed into his or her seat via a 3-point harness. The harness will let the passengers fly upside-down while still secured in the cockpit. Since the flight simulator will be subject to 6g  acceleration, complementary sick bags will be provided upon starting.

In case of a medical emergency of a passenger or if it has been determined that it is unsafe to ride mid-flight, an emergency stop will be activated which will bring the ride to an end. When activated, the ride will right itself upwards while bringing itself to a stop about the center of the ride. This is so when the ride stops, the passengers are not hanging upside down which would be unsafe.

2.1 Locations of possible failure in the structure

Four critical sections in the structure were identified as possible failure sections. These are shown in the following diagram. They ranked from 1  to 4  in order of possible first to fail. Hence section 1  is the one expected to fail first.


PIC
Figure 4: Identification of critical sections in the structure


From bending moment diagram generated during initial runs of simulation it was clear that the bending moment at section 1 was much larger than section 3. This agrees with typical cantilever beam model which the above have very close similarity when considering the cabinet as additional distributed load on the beam. However, this is a dynamic design and not static, hence time dependent bending moment and shear force diagrams are used to validate this. These diagrams were not included in the final simulation software due to time limitation to fully implement them in an acceptable manner. Due to also time limitations analysis for section 2 and 4 were not completed. The design team felt that protecting against failure in section 1 was the most important part at this design stage as this is the most likely failure section. If awarded the design, the team will include full analysis of all sections using finite element methods for most accurate results.

3 Mathematical model of system dynamics

This section explains and shows the derivation of the mathematical model and dynamic equations. These equations are used in the implementation of the software simulator in order to test and validate the design and select the final optimal design parameters.

3.1 Review of the model structure used in the design

There are two rigid bodies: the beam and the supporting column. The cabinet is part of the beam but was analyzed as a rigid body on its own in order to simplify the design by avoiding the determination of moments of inertia for a composite shaped body. The following architectural drawing shows the ride structure.


PIC
Figure 5: Main parts of the ride structure


The ride consists of the main support vertical column attached to a spinning base. Attached to one side of the column is an aluminum beam connected to the column using a drive shaft coupling that allow the beam to spin while attached to the column. A motor supplies the power needed to spin the shaft.

The cabinet is mounted and welded on the beam. The location of the cabinet on the beam is a configurable parameter in the design, and was adjusted during simulation to find an optimal location for the seating cabinet. In final design the cabinet was located at the far end of the beam to achieve maximum passenger felt acceleration.

The passengers are modeled as one rigid body of an equal side solid cube of a mass that represents the total mass of the passengers (maximum of 2 persons) with additional mass to account for the seating weight and a factor of safety. The factor of safety was also an adjustable parameter in the simulation. The following diagram shows the main dimensions of the structure used in the design.


PIC
Figure 6: Main dimensions of ride structure


3.2 Setting up the mathematical model

Euler rigid body dynamic equations of motion are used to determine the dynamic moments due to the rotational motion of the rigid bodies. Principal Body axes, with its origin at the center of mass of each rigid body was used as the local body fixed coordinates system. Newton method is used to obtain the dynamics forces due to translation motion of the beam center of mass and also the center of mass of the cabinet. The column has rotational motion only and no translation motion.

After finding the dynamic forces, the unknown reaction forces at the joint between the beam and the column are solved for. Since these forces are functions of time, simulation was required to check that they remain below yield strength of Aluminium during the ride duration. Analytical solution is difficult due to the nonlinearity of the equations of motion, but a numerical solution of the equations of motion would have been possible.

From beam bending moment diagrams generated for this design, the cross section at the beam/column joint was determined to be the critical section. This is the section which will have the maximum bending moment as well maximum shear force.

During simulation, the current values of the bending moment and shear force at the joint were tracked for each time step taken. The maximum values of these are used to determine the corresponding maximum stress concentration on the section to insure they do not reach 0.55  of yield strength of Aluminum. 0.55  was used to protect against failure in shear which can occur before failure in tension.

In order to minimize the number of parameters to vary in the design, the width of the cabinet was set to be the same as the beam width. The stresses in the beam are calculated based on simple beam theory and not plate theory. Due to time limitation, finite element analysis would was not performed. Finite element analysis would give more accurate stress calculations which would have allowed the design to be free to use less material by using thin plate for the platform and not thick beam as was used.

The following is a summary of the main steps used in the dynamic analysis process

1.
Break the system into 3 separate rigid bodies
2.
Use Euler and Newton methods to determine dynamic loads on each body. Principal body fixed axes are used with the reference point being the center of mass. (called case one analysis or ω = Ω  ).
3.
Draw free body diagram for each body and balance the dynamic loading found in the above step in order to solve for unknown reaction forces.
4.
Apply these reactions forces to the second rigid body connected to the first body by reversing the sign on all vector. These new vectors now act as external loads on the second rigid body.
5.
Perform Euler and Newton analysis on the second body to find its dynamic loads needed to cause it motion.
6.
Make free body diagram for the second body to balance the external forces with the dynamic loads and remembering to use the loads found in step 3 as external loads to this second body.

This diagram below illustrate the different coordinates axes used. The rotating coordinates system that all forces and resolved for is the xyz  . This has its origin at the joint between the beam and the column. This coordinates system is attached to the column and rotates with the column at an absolute angular velocity ωp  . Each rigid body has its own local body fixed coordinates system x′y′z′ . In this design, x′y′z′ have the origin at the center of mass of each rigid body and are aligned with the body principal axes. Hence x′y′z′ is the same as the e1,e2,e3  axes commonly used to mean the principal axes. Therefore ω = Ω  in all cases.  Once dynamic loads are found using x′y′z′ the results are transformed back to the xyz  coordinates system. This way all the results from different rigid bodies are resolved with respect to a common coordinates system xyz  (which is itself a rotating coordinates system).


PIC
Figure 7: relation between rotating coordinates system, body fixed coordinates system, and body principal axes.


3.2.1 Summary of design input and design output

The following tables summarize the input and the output of the overall design. The tables list all the design parameters and the meaning and usage of each. They show what is known at the start of the design and the output from the design and simulation




Parameter name Meaning and usage


ρ  Density of Aluminum 2700  -kg-
m3   , E = 69  GPa, Max tensile 125 MPa, Max yield strength 55 MPa


q  Mass per unit length of the beam


L  Length of the beam


Ls  Distance to the center of cabinet from the left edge of the beam


h  Thickness of the beam (rectangular cross section beam)


b  Width of the beam and cabinet


ω˙p  Angular acceleration of vertical column (zero at steady state)


ω˙s  Angular acceleration of platform and cabinet (zero at steady state)


m  Total mass of cabinet. 175  lbs per person, total of two persons including additional  200  lbs for seats


M  Mass of main support column. Fixed in design


gLimit  Maximum acceleration felt by rider. Must not exceed 6  g


σyield  Yield tensile stress for Aluminum. 55  MPa



Table 1: design input parameters

The following table shows the output of the design based on the above input. Simulation was used to find an optimal set of input parameters in order to achieve the customer specifications




Parameter name Meaning and usage


am  Acceleration time history experienced by passenger. Not to exceed 6g
F
  weld  Reaction forces at joint connecting the beam with the column
Mweld  Reaction moment at joint connecting the beam with the column


ωp  Column angular velocity time history


ωs  Beam angular velocity time history


σ  Direct stress tensor at critical section (joint between beam and column)


τ  Shear stress tensor at critical section (joint between beam and column)


σmax  Maximum direct stress recorded, must remain below yield stress for Aluminium


τmax  Maximum shear stress recorded, must remain below 0.55  of tensile yield stress


amax  Maximum acceleration reached by riders. Must be as close as possible to 6g


vmax  Maximum velocity reached by riders. Typical value from simulation was 180  m.p.h.



Table 2: design output

3.2.2 System dynamic loads and free body diagram

Before starting the derivation, the following two diagrams are given to show the dynamic loads to be balanced with constraint forces. Two free body diagrams used. One for the beam and one for the column.


PIC
Figure 8: Beam dynamics. Balancing dynamic forces to external forces and reactions


After Mweld  and Fweld  are solved for, they are used (with negative signs) as known constraint forces on the column in order to solve for the column's own constraint forces and any external loads. The free body diagram for the column is given below


PIC
Figure 9: Column dynamics. Balance with and external loads and beam transferred loads.


The analysis below shows all five derivations. The first obtains MBeam  (dynamic moment to rotate the beam) using Euler method. The second finds Mcabinet  (dynamic moment to rotate the cabinet) using Euler method, the third uses Newton method to find linear acceleration of center of mass Fcabinet  (dynamic force to translate the cabinet), the fourth finds the linear acceleration of the center of the beam and FBeam  and the final derivation finds Mcolumn  (dynamic moment to rotate the column).

3.3 Beam to column analysis

3.3.1 Finding Mbeam  (beam dynamic moment)

The platform is modeled as a rectangular beam. Its principal moments of inertia are given below.


PIC
Figure 10: Configuration used for finding torque and force at beam/column joint


Let ω be the absolute angular velocity of the local body rotating coordinates x′y′z′ . Let Ω  be the beam (the body) absolute angular velocity. Hence

ωcs = ωpk+ ωsj

But ωcs = Ωbody  , therefore

pict

In component form

pict

Taking time derivative

pict

In component form

pict

The moments of inertia of the beam using its principal axes at the center or mass are

pict

Since ρc = 0  (center of mass is used as reference point) then

M ρc× ¨rp = 0

Moments of inertia cross products are all zero since principal axes is used. The relative angular momentum of the beam becomes

     (         ) (  )
      I1  0   0   Ω1
h  = || 0  I   0|| ||Ω ||
  p  (     2   ) (  2)
       0  0  I3   Ω3

The rate of change of the relative angular momentum of the beam using Euler equations is

     (  )   (                 )
      h˙1      I1Ω˙1 +Ω2Ω3 (I3− I2)
h˙p = ||h˙2|| = || I2Ω˙2 +Ω1Ω3 (I1− I3)||
     ( ˙)   (   ˙             )
      h3      I3Ω3 +Ω1Ω2 (I2− I1)

Therefore, the moment needed to rotate the beam with the angular velocity specified is

Mp = ˙hp

The above components are expressed using in the beam body fixed coordinates system x′y′z′ (which is the same as e1,e2,e3  in this case). These are converted back to the xyz  coordinates system using the following transformation

pict
3.3.2 Finding M
  col  (column dynamic moment)

The main support column has one degree of freedom as it only spins around its z  axes with angular velocity ωp  . Its center of mass does not translate in space. The column has a square cross section. Its height and sectional area were fixed in the design to allow changing the beam and cabinet parameters freely and see the effect on the joint stresses between the beam and the column as the failure point in the design was considered to be the the joint between the beam and the column


PIC

This is a case of one body rotating around its own axes. Therefore,

Mz = I3˙ωp

Where

pict

Where mcol  is the mass of the column. Hence

M     = 1Mr2 ˙ω
 column  6     p

3.3.3 Finding Mcabinet  (cabinet dynamic moment)

The passengers including the cabinet are modeled as solid cube rigid body. The cabinet and the beam rotate with the same absolute angular velocity and act as one solid body. They were analyzed separately as it is easier to find the moment of inertias of each body separately than if both were combined.

The center of mass of the cabinet is at a distance h
2  above the beam where h  is the width of cube which is the same as the beam width.


PIC

Since the cabinet is attached to the platform and is a rigid body as well, the same exact analysis that was made to the beam above can be used for the cabinet. The only difference is that the moments of inertia I1,I2,I3  are different. In this case they are

            1   (     )
I1 = I2 = I3 = 12m b2+ h2

Therefore, the body dynamic moments are

pict

The above components are expressed using the cabinet own principal axes coordinates system  ′ ′′
x yz (local body coordinate systems) which is its principal axes in this case. These are converted back to the xyz  coordinates using the same transformation used for the beam

pict
3.3.4 Finding Fcabinet  (cabinet dynamic linear force)

To find Fm = ma  for the cabinet, Newton method is used as follows


PIC
Figure 11: Rotating coordinates system xyz  used to find passenger acceleration


The rotating coordinates system xyz  has its origin at the beam column joint. xyz  is attached to the column and rotates with the column with angular velocity ωpk  . The center of mass of the cabinet shown above as the circle p  , is at distance Ls  from the origin O.

All calculations are expressed using unit vectors of the rotating coordinates system and are valid for all time. In the rotating coordinates system, point p,  the center of mass of cabinet, appears as shown in the following diagram. In this diagram 𝜃  is the angle p  makes with the z  axes, where 𝜃 = ωst  and 𝜃˙= ωs


PIC
Figure 12: View of passenger head in the rotating coordinates system xyz


Using the above diagrams, the absolute velocity of p  is found as follows

pict

Hence the absolute velocity of p  is

pict

The absolute acceleration of p  is found from

pict

Therefore the absolute acceleration of the passenger is

pict

Simplifying gives

pict

The above is expressed using the common xyz  rotating coordinate system

3.3.5 Finding Fbeam  (beam dynamic translational force)

The linear acceleration of the center of mass of platform, which is located at distance L
2  from the origin o  of the xyz  rotating coordinates system.


PIC
Figure 13: Rotating coordinates system xyz  used to find beam center of mass acceleration


Therefore

pict

Hence

pict

Therefore

pict

The above is expressed using the xyz  rotating coordinates system.

3.3.6 Using free body diagram and solving for constraint forces

The dynamic forces have been found from above. The are balanced with constraint forces and any external loads using free body diagram. The following diagram shows the balance between dynamic forces and moments and external forces. Mweld  below is used to represent all constraint moments at the joint between the beam the column, including the extra torque needed to rotate the beam


PIC

Taking moments at point o  , the left end of the beam which is the origin of the rotating coordinates system xyz

pict

Hence

       ( L          )                   ( L        )  (     b  )
Mweld =  2Mg  +Lsmg  i+ Mbeam + Mcabinet+  2j× Fbeam  +  Lsj+ -2k  ×Fcabinet

The force vector at the joint is

pict
Bending moment and shear force calculations  Now that the constraint forces are solved for from the above analysis, the bending moment and shear force diagram are formulated. The moments will be a function of distance from the beam/column joint. Let BM (ξ)  be the moments vector at distance ξ  along the beam length. There will be 3 components to this moment. Bending Mx  , torsional My  and twisting Mz  .
PIC
Figure 14: Finding the bending moment at different locations along the span of the beam

Let the weight be per unit length of the beam which isM-g
L  be q  . In the following, the notation ⟨ξ − x⟩ is used to indicate that the term is effective only when ⟨ξ − x⟩ is positive. Let the distance to start of the cabinet be
        b
α= Ls − 2
Where b  is the width of the cabinet.
                           (ξ       )   (ξ− α     mg        )
BM (ξ)= Mweld+ (ξj×Fweld)+  2j× −qξk  +  --2-j× − -b-(ξ − α )k ⟨ξ− α⟩
In component form, the bending moment will be BMx  (ξ)  and The torsion moment will be BMy (ξ)  and the twisting moment will be BMz  (ξ).  Let SF (ξ)  be the shear force vector at distance ξ  . Hence
                   mg-
SF (ξ)= Fweld− qξk −  b (ξ − α )⟨ξ− α⟩K
The above completes the mathematical derivation of the dynamics of the system. The next step is to implement this model and use simulation to validate it and design for an optimal set of parameters. Finding shear and direct stress from bending and shear forces  The result of the above calculations is the moments and forces at the joint between the beam and the column and using BM  (ξ)  and SF (ξ)  at any other section in the beam. The next step is to use these to obtain complete description of stress state at the section. Due to lack of time finite element analysis was not performed. Therefore, basic beam theory equations were used for stress calculation. Care was taken to insure that the beam cross section selected had thickness not less that its width. Having a thin beam would require analysis using plate theory making it much more complicated. The disadvantages of this method is that the beam was much heavier than needed if thin beam was used, but the advantage is that the stress equations used are known to be valid in this case. Given the moments Mx, My,Mz  and the forces Fx,Fy,Fz  all the cross section, the following equations were used. These equations assume a rectangle beam cross section of thickness h  and width b  and that h ≥b
pict
Torsional stress was not fully developed in this design since it is a rectangular cross section and would require finite element analysis. The beam is expected to fail due to bending moment Mx  and this is what the rest of the analysis address. Future analysis of stress concentration will use finite element analysis and will take torsion stress into account.

3.4 Column dynamic analysis

In the above section the constraint forces in the beam/column joints were found. These are now used as external forces on the column with an opposite sign. Free body diagram is used for the column in order to find the constraint forces and external loads acting on the column. The following diagram shows the free body diagram used


PIC
Figure 15: Dynamic load balance between column and external loads


Taking moments at the joint between the column and the ground

                   (  H-        )
T + Mweld2− Mweld+  − 2 k × −Fweld = Mcolumn

Solving for the unknown constraint force N  and the external torque T

                    (H         )
Mweld2 +T = Mcolumn−  -2 k× Fweld + Mweld

The torque T  is unknown at this stage and has to be determined by other means to obtain complete solution. This is the external torque needed to accelerate the column during ramp up and to decelerate it during ramp down phases. Combining all the unknowns into one term called Mweld3  , the above reduces to

                 (          )
                  H-
Mweld3 = Mcolumn−  2 k ×Fweld + Mweld

The balance equation for forces gives

pict

Now that all loads acting on the column are found, bending moment and shear force diagrams can be also be made or finite element analysis used in order to determine the stress state inside the column at every section.

4 Simulation of the dynamic equations found

4.1 Review of the simulation

The simulation accepts as input all the parameters shown in table 1 on page 9. The goal of the simulation is to verify visually the dynamics and to allow the selection of correct sizes for the structure and to insure that the acceleration does not exceed 6g  using the selected parameters. Based on the simulation, one optimal set of values was selected and given in the conclusion section. The simulator displays tables showing all the current values for stress and moments found at the beam/column joint. It keeps track of the maximum stress values reached and uses these to determine the maximum stress using the equations shown above.

This diagram shows an overview of the user interface. This software can be run from the project web site located at http://12000.org/my_notes/mma_demos/EMA542_project/index.htm


PIC
Figure 16: overview of simulator user interface


4.2 Simulation output, time histories and discussion of results

All these tables and results below are generated from the final design using the selected final optimal parameters.


PIC
Figure 17: dynamic loads at the end of ride using optimal design values



PIC
Figure 18: critical section current and maximum moments and stresses



PIC
Figure 19: optimal set of parameters obtained from simulation.



PIC
Figure 20: simulator keeps track of maximum g  felt by passenger to insure it does not exceed 6g



PIC
Figure 21: acceleration and velocity of passenger time history and angular velocity time history of beam and column



PIC
Figure 22: time histories using the ramp down option used after reaching 6g goal


4.3 Discussion and analysis of results

The following table gives the optimal design parameters found by simulation of the derived model in order to achieve the customer requirements.





parameter value description



beam mass 10.4  ton one ton is 2000 lbs



beam width 1 meters



beam thickness 1  meters



beam length 3.5 meters



cabinet mass 560  lbs includes 2 passengers, seating, frame and factor of safety



cabinet height 1  meters



cabinet width 1  meters



column mass 17.1 ton



column cross section 3  by 3  meters



maximum bending moment Mx  175  KNm



maximum torsion moment My  7.6  KNm



maximum twisting moment Mz  28  KNm



maximum shear force Fx  − 15.96  KN



maximum shear force Fy  − 85  KN



maximum shear force Fz  107  KN



maximum direct stress σx  1.055  MPa Below tensile yield. Pure Aluminium has 10  MPa. and



maximum direct stress σy  0.046  MPa Aluminium 6061-O yields at 200  MPa.



maximum direct stress σz  0.172  MPa



maximum shear stress τx  − 23.94  KPa



maximum shear stress τy  − 128.6  KPa



maximum shear stress τz  150.5  KPa




Table 3: design output for loading and forces using optimal parameters found

It was found that in order to be able to achieve the 6g  limit and not exceed it, the acceleration have to put turned off well before the 6g  is detected. This can be seen by examining the passenger acceleration expression from above, which is

pict

We can see that, by letting ˙ωs  and ˙ωp  then the acceleration becomes

    (                  )   (                 )
a = i − bω2 sin𝜃− ω2b sin𝜃 + j 2bω ω cos𝜃− ω2L  − kbω2 cos𝜃
       2  s      p2           2 p s       p s    2  s

Even though from now on the angular velocities ωs  and ωp  are constant, this does not imply that a  will become constant. Since 𝜃  is still changing in time, then a  will still fluctuate in periodic fashion from now on. Hence the passenger acceleration can still exceed 6g  if we were to turn off the ramp up acceleration too close to 6g  . For this reason the value the acceleration was turned off at 5.8g  in order to final value of 5.98g  as felt by the passengers.

4.4 Cost analysis

Based on the above result and using the mass needed, the following table gives a summary of cost for construction of the ride





item cost description



cost of Aluminum alloy 6061-O $0.8 per lb. can depend on market conditions
beam material cost (10.4 ton) $16,000 (10.4)(2000)(0.8)
column material cost (17.1 ton) $27,360 (17.1)(2000)(0.8)



cabinet material cost (500 lb.) $446.5 PIC



Labor cost for construction $12,000 300 labor hrs @ 40 per hr.



Equipment and labor insurance $10,000



Management cost (one manager) $4,000 50 hrs @ $80 per hr.



Electric spindel motors for column and beam $10,000 2 @ $5,000 PIC



Total cost $79,806




Table 4: cost estimate

The major part of the cost is for material. This is due to the use of thick beam and column. This allowed the use of basic beam theory stress analysis. This cost however can be reduced by the use of plate theory or numerical finite elements methods in order to be able to safely used less material and reduce the thickness of the beam and column while insuring accurate stress calculations.

5 Conclusions of results and future work

The final design given above meets the requirement specification that the customer provided. Using simulation, it was possible to validate the equations found and to confirm that the beam/column section is safe for the selected optimal parameters.

The selected parameters allow the passengers to reach almost 6g  in 12  seconds using a ride that consist of two noncollinear angular velocities. There are many different profiles that could have been selected to achieve this goal. The set selected reached the closest to 6g  without crossing over and that is why it was selected. The following is the final design used





parameter value description



maximum g  reached 5.98  g After many simulations this was selected.
time to reach maximum g  5.8  sec.
maximum passenger velocity reached 180  m.p.h. calculated using finite difference from acceleration data



steady state ωp  reached 2.16  rad/sec. This is the column angular velocity in steady state



steady state ωs  reached 9.7  rad/sec. This is the beam angular velocity in steady state



initial ramp up ˙ω
 p  0.2  rad/sec.2  column supplied ramp up angular acceleration



initial ramp up ˙ωs  0.9  rad/sec.2  beam supplied ramp up angular acceleration



ramp down ω˙p  0.2  rad/sec.2  symmetrical shape to ramp-up as seen in above plot.



ramp down ω˙s  − 0.9  rad/sec.2  symmetrical shape to ramp-up as seen in above plot.




Table 5: ride statistics based on optimal design parameters

The cost estimate is $79,800. The material cost was the major part of this cost. This was due to the use of simple beam theory for stress analysis equations which required the use of a thick beam in order for the stress equations to be valid. The maximum stress of σ   = 1.055
 max  MPa reached is well below the yield strength of Aluminum. Therefore, the use of finite element stress analysis or advanced plate theory would have allowed the reduction of the size of the beam while at the same time using accurate stress calculations. This would have resulted in lower cost in material. If awarded this contract, finite element would be used in order to lower the cost of material.

5.1 Future work and possible design improvement

The following are items that can be improved in the current design given additional time to perform

1.
The beam and column weight can be reduced significantly by using plate shell stress analysis. This should reduce the material cost. This design used simple beam theory stress analysis which required the use of thick beam. This caused the beam to become too thick. It will be possible to have thinner beam and still not reach the yield strength. Using finite element method will allow this investigation.
2.
There are additional possible cross sections to consider for failure analysis. This design concentrated on the most likely section based on beam theory. Using finite element software will allow one to more easily analyze the full structure more easily than was done in current design based on simple beam theory.
3.
Torsional and twisting stress analysis were not addressed in this design due to time limitation. It is however expected that the beam will fail in bending.

6 Appendix

6.1 Use of simulator to validate different design parameters

These are selected screen shots showing different configurations tested during simulation in order to find an optimal one. These show the effect of changing the dimensions of the structure and the spin rates.


PIC
Figure 23: Changing the structure dimensions to select optimal design using simulation


6.2 References

1.
Daniel C. Kammer, EMA 542 lecture notes and printed book. University of Wisconsin, Madison. 2013
2.
R. C. Hibbeler. Engineering Mechanics, Dynamics. Prentice Hall, NJ  2011
3.
Aluminium page at Wikipedia http://en.wikipedia.org/wiki/Aluminium
4.
Moments of inertia page at Wikipedia http://en.wikipedia.org/wiki/List_of_moments_of_ inertia
5.
Density of materials page http://physics.info/density/
6.
Beam design formulas with shear and moment diagrams book, AWC council, 2007, Washington, DC.