Scribd
Upload
13 views40 pages

Vibrationdata Module

The document presents a two-degree-of-freedom (2-DOF) numerical example for determining modal parameters of a system, including mass and stiffness matrices, natural frequencies, and mode shapes. It also discusses a beam subjected to base acceleration, detailing material properties, responses to be calculated, and methods for achieving accuracy through finite element analysis. The results include comparisons of classical solutions and finite element analysis outcomes, highlighting the importance of element count for accurate stress calculations.

Uploaded by

moyojr60
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
13 views40 pages

Vibrationdata Module

The document presents a two-degree-of-freedom (2-DOF) numerical example for determining modal parameters of a system, including mass and stiffness matrices, natural frequencies, and mode shapes. It also discusses a beam subjected to base acceleration, detailing material properties, responses to be calculated, and methods for achieving accuracy through finite element analysis. The results include comparisons of classical solutions and finite element analysis outcomes, highlighting the importance of element count for accurate stress calculations.

Uploaded by

moyojr60
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 40

https://StudentLabs.uOttawa.

ca
1. Launch Matlab
https://vibrationdata.wordpress.com/2022/09/28/vibrationdata-matlab-stand-
alone-executable-gui-package/
Two-DOF Numerical Example 1

Determine modal parameters for system

m2 = 1.0 lbm

k2 = 10,000 lbf/in

m1 = 2.0 lbm

k1 = 12,000 lbf/in

14
15
16
17
18
19
mass =

0.0052 0
0 0.0026

stiffness =

22000 -10000
-10000 10000

Natural Participation Effective


Mode Frequency Factor Modal Mass
1 182.6 Hz 0.08631 0.007449
2 414.7 Hz -0.01796 0.0003227

modal mass sum = 0.007772

ModeShapes =

9.4707 -10.1640
14.3740 13.3935

Output arrays:
two_dof_mass, two_dof_stiffness, ModeShapes

20
Beam Subjected to Base Acceleration Example 2

• Material:
aluminum with E=1.e+07 lbf/in2

• Length: 24 inch

y(x, t) • Circular cross section with 1.0 inch


diameter
w(t)
• Area: 0.785 in2

• Area Moment of Inertia: 0.0409 in4


• Responses to determine: • Mass/length: 0.0785 lbm/in
2.035e-04 lbf sec2/in2
Acceleration at free end
Bending stress at fixed end • 5% damping all modes (Q=10)
• Solve for “exact” solution using classical methods • Stress concentration factor = 1
• Also solve using finite element analysis

• Determine how many elements are needed for < 3% accuracy for each response
parameter
Step 1

Step 2
Step 1

Step 2

Step 3
• Excite beam at fixed end

• Calculate bending stress


at fixed end

Base excitation
= 2 G at 47.72 Hz

Response at x= 0 inch

Transfer Magnitude
= 1 G/G Accel

Distance from neutral


axis = 0.5 in

Bending Stress =
4105 psi at 0 in
• Excite beam at fixed end

• Calculate acceleration at
free end

Base excitation
= 2 G at 47.72 Hz

Response at x=24 inch

Transfer Magnitude
= 15.69 G/G

Accel = 31.38 G

Rel Vel = 40.32 in/sec

Rel Disp = 0.1345 in


Beam Response to
Base ExcitationOne
Element Model

Step 1

Step 2

Step 3

Step 4
dof = 2
Finite element model
stiff =

1.0e+04 * E1

0.0426 -0.5113 N2
N1
-0.5113 8.1812

mass = Displacement vector

X T =  y1 1 2 
0.0018 -0.0061
-0.0061 0.0268
y2

Total mass = 1.885 lbm


y1 is constrained as a dof at “fixed” end
Mode Freq (Hz) 1

1 47.94 1 is constrained to zero because no


2 472.4 rotation at fixed end
• Excellent agreement for first
mode

• Large frequency error for


second mode

• Higher modal frequencies


missing from FEA curve
Beam Free End Accel
= 31.39 G

Beam Free End Pseudo Velocity


= 37.83 in/sec

Beam Fixed End Stress


= 3766 psi

Response Classical One Element FEA Error


Free End Accel (G) 31.38 31.39 0.03%
Fixed End Bending Stress (psi) 4105 3766 -8.26%
• Continue with the previous example

• Steady-state sine: 2 G input at 47.72 Hz

• Recreate the numbers in the table

Number of Fixed End Bending Stress


Elements (psi)
1 3766
2 4084
4 4104
8 4105
16 4105

• A higher number of elements theoretically yields a more accurate stress


level, but numerical error can occur if too many elements are used

• Typically, a higher number of elements are needed for accurate stress


calculation as compared to acceleration and displacement calculation

You might also like