Computational Mechanics Research and Support for Aerodynamics ...

Computational Mechanics Research and Support for Aerodynamics ... Computational Mechanics Research and Support for Aerodynamics ...

07.01.2015 Views

0 0 0 1 0 0 B 0 = E -1 0 0 1 0 1 0 T 0 0 1 1 0 0 B 1 = E -1 [ 0 0 0] T 3.2.1.2. Conversion of wind velocity to wind pressure The conversion of wind velocity to wind pressure is necessary because the wind is described in terms of velocity and in the finite element analysis the velocities are attached to the nodes and pressure is applied to elements. The loading on elements would provide more accurate results compared to the node loading. The pressure due to wind velocity can be computed using the following conversion: Where, P = Pressure (N/m 2 ) C = Drag coefficient D = Density of air (1.25 Kg/m 3 ) V = Speed of air (m/s) The drag coefficient is a dimensionless quantity and it depends upon the surface area under the influence of attacking wind. The drag coefficients for flat and cylindrical body are approximated as follows: Drag coefficient for flat surface = 1 Drag coefficient for cylindrical surface = 0.67. 3.2.1.3. Development with the FEM After understanding the spring-damper system in LS-DYNA, it was necessary to study the motion of a node on a surface. To perform the analysis, a flat surface was created with two beam element and two massless nodes attached to it (Figure 3.13). These nodes were given a velocity and constraint so as to follow the beam element from start to end. TRACC/TFHRC Y1Q3 Page 58

̈ With the success to the above analysis, a spring element was attached to the nodes and a block was mounted on it (Figure 3.14). To verify the proper behavior of the spring, an initial displacement was given and it produced the appropriate result. After verifying the behavior with the above model, the task was to give an appropriate loading to the vertical face of block in order to understand the effects of wind pressure (Figure 3.15). Two types of wind pressure were given to the block, a constant load and a sine load. The properties of the model are: Mass = 5.88 Kg, Stiffness = 60 N/mm, Length of spring = 50 mm, Amplitude = 0.01, Frequency = 100 Hz, Forcing function = Amplitude*SIN(2*3.14*Frequency*TIME) Figure 3.16 shows the graph of horizontal displacement of the block with the application of sinusoidal wave as a forcing function. Theoretically, the equation of motion of the block with vertical springs and horizontal sine forcing function (Figure 3.15) is given by: √ (3.6) Where, m = Mass of the block, k = Stiffness of spring, l = original length, x = Displacement from static equilibrium position F 0 = Forcing function (sine wave) 3.2.1.4. Ford F-800 truck development A Ford F-800 truck model (Figure 3.17) was taken from the National Crash Analysis Center (NCAC) website [2]. The raw file contained a FEM model of the truck with six rigid walls attached to each tire and forward velocity. This model was studied and modified according to the requirements of the TRACC/TFHRC Y1Q3 Page 59

0 0 0 1 0 0<br />

B 0 = E -1 0 0 1 0 1 0 T<br />

0 0 1 1 0 0<br />

B 1 = E -1 [ 0 0 0] T<br />

3.2.1.2. Conversion of wind velocity to wind pressure<br />

The conversion of wind velocity to wind pressure is necessary because the wind is described in terms of<br />

velocity <strong>and</strong> in the finite element analysis the velocities are attached to the nodes <strong>and</strong> pressure is<br />

applied to elements. The loading on elements would provide more accurate results compared to the<br />

node loading.<br />

The pressure due to wind velocity can be computed using the following conversion:<br />

Where,<br />

P = Pressure (N/m 2 )<br />

C = Drag coefficient<br />

D = Density of air (1.25 Kg/m 3 )<br />

V = Speed of air (m/s)<br />

The drag coefficient is a dimensionless quantity <strong>and</strong> it depends upon the surface area under the<br />

influence of attacking wind. The drag coefficients <strong>for</strong> flat <strong>and</strong> cylindrical body are approximated as<br />

follows:<br />

Drag coefficient <strong>for</strong> flat surface = 1<br />

Drag coefficient <strong>for</strong> cylindrical surface = 0.67.<br />

3.2.1.3. Development with the FEM<br />

After underst<strong>and</strong>ing the spring-damper system in LS-DYNA, it was necessary to study the motion of a<br />

node on a surface. To per<strong>for</strong>m the analysis, a flat surface was created with two beam element <strong>and</strong> two<br />

massless nodes attached to it (Figure 3.13). These nodes were given a velocity <strong>and</strong> constraint so as to<br />

follow the beam element from start to end.<br />

TRACC/TFHRC Y1Q3 Page 58

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!