Modal Analysis of a Beam (SI Units)

Transcription

1 APPENDIX 1a Modal Analysis of a Beam (SI Units) Objectives Perform normal modes analysis of a cantilever beam. Submit the file for analysis in MSC.Nastran. Find the first three natural frequencies and mode shapes of the beam. A-1a-1

2 A-1a-2

3 APPENDIX 1a Modal Analysis of a Beam (SI Units) Model Description: The goal of this example is to find the first 3 modes of a beam pinned at both ends. Figure A-1a.1 below is a finite element representation of the beam. One end is contrained in all translations and the other is free to move in the X. Both ends are held in the X-rotation. Figure A-1a.1-Grid Coordinates and Element Connectivities A-1a-3

4 Table A-1a.1 Length Elastic Modulus 1.0 x 10 3 mm x 10 5 MPa Density x 10-9 N-sec 2 /mm 4 Poisson s Ratio 0.32 Area 5 x 10 3 mm 2 I x 10 6 mm 4 I mm 4 Force 1 x 10 7 N A-1a-4

5 APPENDIX 1a Modal Analysis of a Beam (SI Units) Hand Calculations fn = K n EIg 1 2 2π Wl 4 fn = Kn π ( ) ( 5 10 )( ) fn = Kn( ) From Theory Mode Kn fn Hz Hz Hz A-1a-5

6 Suggested Exercise Steps Explicitly generate a finite element representation of the beam structure. (i.e., the grids (GRID) and element connectivities (CBAR) should be defined manually.) Define material (MAT1) and element (PBAR) properties. the fixed boundary constraints (SPC1). Prepare the model for a normal modes analysis (SOL 103 and PARAMS). PARAM, COUPMASS, 1 EIGRL (To select Lanczos) Generate an input file and submit it to the MSC.Nastran solver for normal modes analysis. Review the results, specifically the eigenvalues. A-1a-6

7 APPENDIX 1a Modal Analysis of a Beam (SI Units) ID SEMINAR,PROB1 CEND BEGIN BULK A-1a-7

8 A-1a-8

9 APPENDIX 1a Modal Analysis of a Beam (SI Units) ENDDATA A-1a-9

10 Exercise Procedure: 1. Users who are not utilizing MSC.Patran for generating an input file should go to Step 10, otherwise, proceed to step Create a new database named probap1.db. File/New Database New Database Name probap1 In the New Model Preference form set the following: Tolerance Analysis Code: Default MSC/NASTRAN 3. Activate the entity labels by selecting the Show Labels icon on the toolbar. Show Labels 4. Create a curve. Geometry Action: Create Object: Curve Method: XYZ Vector Coordinates List <1000, 0, 0> Origin Coordinates List [ 0, 0, 0] A-1a-10

11 APPENDIX 1a Modal Analysis of a Beam (SI Units) 5. Create the finite element model and mesh the surface. Finite Elements Action: Create Object: Mesh Type: Curve Global Edge Length 100 Curve List Curve 1 6. Create nodal displacements. Loads/BCs Action: Create Object: Displacement Type: Nodal New Set Name disp1 Input Data... Translations <T1 T2 T3> <0 0 0> Rotations <R1 R2 R3> <0,, > Select Application Region... Geometry Select Geometry Entities Point 1 Add New Set Name disp2 Input Data... Translations <T1 T2 T3> <, 0 0> A-1a-11

12 Rotations <R1 R2 R3> <0,, > Select Application Region... Geometry Select Geometry Entities Point 2 Add New Set Name permanent_constraint Input Data... Translations <T1 T2 T3> <,, 0> Rotations <R1 R2 R3> <0, 0, > Select Application Region... Select Geometry Entities Curve 1 Add 7. Create a set of material properties for the bar. Materials Action: Create Object: Isotropic Method: Manual Input Material Name mat_1 Input Properties... Elastic Modulus = E5 Poisson Ratio = 0.32 A-1a-12

13 APPENDIX 1a Modal Analysis of a Beam (SI Units) Density = Cancel E-9 8. Define the bar properties. Properties Action: Dimension: 1D Create Type: Beam Property Set Name Input Properties... bar Material Name m:mat_1 (Select from Material Property Sets box) Bar Orientation Coord 0.2 Area 5E3 [Inertia 1,1] E6 [Inertia 1,2] 1.0 Select Members Curve 1 Add 9. Now, you will generate the input file for analysis. Analysis Action: Analyze Object: Entire Model Method Analysis Deck Job Name Solution Type... probap1 A-1a-13

14 Solution Type: Solution Parameters... <deselect Automatic Constraints> Mass Calculation: Data Deck Echo: Subcase Create... Available Subcases Subcase Parameters... Number of Desired Roots = 3 Cancel NORMAL MODES Automatic Constraints Coupled None Default <to close form> An MSC.Nastran input file called probap1.bdf will be generated. The process of translating your model into an input file is called Forward Translation. The Forward Translation is complete when the Heartbeat turns green. MSC.Patran Users should proceed to step 11. A-1a-14

15 APPENDIX 1a Modal Analysis of a Beam (SI Units) Generating an input file for MSC.Nastran Users: MSC.Nastran users can generate an input file using the data from Table A- 1a.1. The result should be similar to the output below. 10. MSC.Nastran Input File: probap1.dat SOL 103 TIME 600 CEND TITLE = Normal Modes Example (SI UNITS) SUBCASE 1 METHOD = 1 SPC = 1 VECTOR(SORT1,REAL)=ALL BEGIN BULK PARAM COUPMASS1 EIGRL PBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR MAT GRID GRID GRID GRID GRID GRID GRID GRID GRID GRID GRID SPC SPC ENDDATA A-1a-15

16 Submit the input file for analysis 11. Submit the input file to MSC.Nastran for analysis. 11a. To submit the MSC.Patran.bdf file for analysis, find an available UNIX shell window. At the command prompt enter: nastran probap1.bdf scr=yes. Monitor the run using the UNIX ps command. 11b. To submit the MSC.Nastran.dat file for analysis, find an available UNIX shell window. At the command prompt enter: nastran probap1 scr=yes. Monitor the run using the UNIX ps command. 12. When the run is completed, edit the probap1.f06 file and search for the word FATAL. If no matches exist, search for the word WARNING. Determine whether existing WARNING messages indicate modeling errors. 13. While still editing probap1.f06, search for the word: E I G E N (spaces are necessary) What are the first three modes? 1st = Hz 2nd = Hz 3rd = Hz A-1a-16

17 A-1a-17 R E A L E I G E N V A L U E S MODE EXTRACTION EIGENVALUE RADIANS CYCLES GENERALIZED GENERALIZED NO. ORDER MASS STIFFNESS E E E E E E E E E E E E E E E Compare the results obtained in the.f06 file with the following results: Comparison of Results APPENDIX 1a Modal Analysis of a Beam (SI Units)

18 15. MSC.Nastran Users have finished this exercise. MSC.Patran Users should proceed to the next step. 16. Proceed with the Reverse Translation process, that is importing the probap1.op2 results file into MSC.Patran. To do this, return to the Analysis form and proceed as follows: Analysis Action: Read Output2 Object: Result Entities Method Translate Select Results File... Select Results File probap1.op2 When the translation is complete bring up the Results form. Results Action: Create Object: Deformation Select Results Cases Select Deformation Result Default, Mode 1:Freq= Eigenvectors, Translational To reset the graphics, click on this icon: Reset Graphics You can go back and select any Results Case, Fringe Results or Deformation Results you are interested in. Quit MSC.Patran when you are finished with this exercise. A-1a-18