6WUXFWXUDO'\QDPLFV 7RROER[ )RU8VHZLWK0$7/$% 8VHU V*XLGH (WLHQQH%DOPqV 9HUVLRQ
+RZWR&RQWDFW6'722/6 Phone
Fax
6'722/6 Mail UXH9HUJQLDXG 3DULV )UDQFH KWWS:ZZZVGWRROVFRP Web KWWS:ZZZLQWHUQDWLRQDOWHFFRP FRPSVRIWV\VPDWODE Newsgroup VXSSRUW#VGWRROVFRP Technical support VXJJHVW#VGWRROVFRP Product enhancement suggestions LQIR#VGWRROVFRP Sales, pricing, and general information 6WUXFWXUDO'\QDPLFV7RROER[8VHU V*XLGH &RS\ULJKWE\6'722/6 7KHVRIWZDUHGHVFULEHGLQWKLVGRFXPHQWLVIXUQLVKHGXQGHUDOLFHQVHDJUHHPHQW 7KHVRIWZDUHPD\EHXVHGRUFRSLHGRQO\XQGHUWKHWHUPVRIWKHOLFHQVHDJUHHPHQW 1RSDUWRIWKLVPDQXDOLQLWVSDSHU3')DQG+70/YHUVLRQVPD\EHFRSLHGSULQWHGSKRWRFRSLHG RUUHSURGXFHGLQDQ\IRUPZLWKRXWSULRUZULWWHQFRQVHQWIURP6'722/6 Structural Dynamics Toolbox is a registered trademark of Etienne Balmès MATLAB is a registered trademark of The MathWorks, Inc. Other products or brand names are trademarks or registered trademarks of their respective holders. Contents
1Preface 7 1.1 Getting started ...... 8 1.2 Understanding the Toolbox architecture ...... 11 1.2.1 Layers of code ...... 11 1.2.2 Graphical User Interfaces (GUI) ...... 12 1.2.3 User Interfaces (UI) Command functions ...... 15 1.3 Typesetting conventions and scientific notations ...... 16 1.4 Release notes for version 4.1 ...... 18 1.4.1 Key features ...... 18 1.4.2 Detail by function ...... 18 1.5 Release notes for version 4.0 ...... 21 1.5.1 Key features ...... 21 1.5.2 Detail by function ...... 21
2 Fundamental concepts 25 2.1 I/O shape matrices ...... 26 2.2 Normal mode models ...... 28 2.3 Damping ...... 30 2.3.1 Viscous damping in the normal mode model form ...... 31 2.3.2 Damping in finite element models ...... 32 2.4 State space models ...... 32 2.5 Complex mode models ...... 34 2.6 Pole/residue models ...... 36 2.7 Parametric transfer function ...... 37 2.8 Non-parametric transfer function ...... 38 2.9 Nodes ...... 39 2.10 Finite element model description matrices ...... 40 2.11 Material properties ...... 42 2.12 Element property matrices ...... 43
1 CONTENTS
2.13 DOF definition vector ...... 44 2.14 DOF selection ...... 45
3 Experimental modal model identification 47 3.1 Getting started with the iiplot interface ...... 48 3.2 Importing FRF data ...... 50 3.3 The standard identification procedure using id rc ...... 51 3.3.1 The id rc procedure step by step ...... 53 3.3.2 Background theory ...... 58 3.3.3 When id rc fails ...... 59 3.4 Transformations of the basic pole/residue model ...... 60 3.4.1 Multiplicity (minimal state-space model) ...... 61 3.4.2 Reciprocal models of structures ...... 62 3.4.3 Normal mode form ...... 64 3.5 Other identification algorithms ...... 68 3.5.1 Direct system parameter identification algorithm ...... 68 3.5.2 Orthogonal polynomial identification algorithm ...... 69
4 Pre- and post- processing a modal test 71 4.1 Preparing a modal test ...... 73 4.1.1 Wire-frame display ...... 73 4.1.2 Basic sensor/shaker configurations ...... 75 4.1.3 General test setups ...... 76 4.1.4 Sensor/shaker placement ...... 78 4.2 Expansion methods ...... 79 4.2.1 Unified perspective on expansion methods ...... 79 4.2.2 Basic interpolation methods for unmeasured DOFs ...... 80 4.2.3 Subspace based expansion methods ...... 81 4.2.4 Model based expansion methods ...... 83 4.3 Correlation criteria ...... 84 4.3.1 Shape based criteria ...... 84 4.3.2 Energy based criteria ...... 85 4.3.3 Correlation of FRFs ...... 86 4.4 Direct uses of experimentally derived models ...... 87
2 5 Basic finite element analysis 89 5.1 Declaring finite element models ...... 90 5.1.1 Basic model declaration (truss example) ...... 90 5.1.2 Building models with femesh ...... 91 5.1.3 Importing models from other codes ...... 92 5.1.4 Importing element matrices from other codes ...... 93 5.2 Computing the response of a model ...... 94 5.2.1 Properties, assembly, boundary conditions ...... 94 5.2.2 Static responses ...... 96 5.2.3 Normal modes (partial eigenvalue solution) ...... 96 5.2.4 Models of structural dynamics ...... 97 5.2.5 Manipulating large finite element models ...... 99 5.3 Learning feplot features ...... 100
6 Model reduction 103 6.1 Theoretical aspects ...... 104 6.1.1 General framework ...... 104 6.1.2 Static correction to normal mode models ...... 105 6.1.3 Static correction with rigid body modes ...... 107 6.1.4 Other standard reduction bases ...... 108 6.1.5 Substructuring ...... 109 6.1.6 Reduction for parametrized problems ...... 110 6.2 Sample applications ...... 111 6.2.1 Component mode synthesis ...... 111 6.2.2 Substructuring using superelements ...... 113
7 Parametrized finite element problems 115 7.1 The upcom interface ...... 116 7.1.1 General framework ...... 116 7.1.2 upcom parametrization for full order models ...... 117 7.1.3 Getting started with upcom ...... 118 7.1.4 Reduction for variable models ...... 119 7.1.5 Predictions of the response using upcom ...... 120 7.2 Finite element model updating ...... 121 7.2.1 Error localization/parameter selection ...... 122 7.2.2 Update based on frequencies ...... 122 7.2.3 Update based on FRF ...... 123
3 CONTENTS
8 Supported Element functions 125 bar1 ...... 127 beam1 ...... 128 celas ...... 130 elem0 ...... 131 hexa8, hexa20, penta6 ...... 136 mass1,mass2 ...... 137 quad4, quadb ...... 138 rigid ...... 140 tria3, tria6 ...... 142 tetra4 ...... 144
9 Function reference 145 basis ...... 151 commode ...... 154 comstr ...... 156 db, phaseb ...... 158 fe2ss ...... 159 fecom ...... 161 femesh ...... 172 feplot ...... 183 fesuper ...... 188 feutil ...... 191 fe c ...... 193 fe ceig ...... 195 fe coor ...... 196 fe eig ...... 197 fe exp ...... 200 fe load ...... 203 fe mat ...... 204 fe mk ...... 205 fe norm ...... 207 fe reduc ...... 208 fe sens ...... 211 fe stres ...... 215 fe super ...... 217 idcom ...... 221 idiplot ...... 226 idopt ...... 228 id dspi ...... 231
4 id nor ...... 232 id poly ...... 234 id rc, id rcopt ...... 235 id rm ...... 239 iicom ...... 241 iimouse ...... 248 iiplot ...... 252 ii cost ...... 255 ii mac ...... 256 ii mmif ...... 266 ii plp ...... 269 ii poest ...... 270 ii pof ...... 272 nasread ...... 274 naswrite ...... 277 nor2res, nor2xf ...... 278 nor2ss ...... 280 psi2nor ...... 283 qbode ...... 285 res2nor ...... 286 res2ss, ss2res ...... 287 res2tf, res2xf ...... 289 rms ...... 290 setlines ...... 291 sdtdef ...... 292 sdth ...... 293 skyline ...... 294 sp util ...... 296 ufread ...... 298 ufwrite ...... 304 upcom ...... 305 up freq, up ifreq ...... 313 up ixf ...... 314 xfopt ...... 315
Bibliography 319
Index 323
5 CONTENTS
6 1
Preface
1.1 Getting started ...... 8 1.2 Understanding the Toolbox architecture ..... 11 1.2.1 Layers of code ...... 11 1.2.2 Graphical User Interfaces (GUI) ...... 12 1.2.3 User Interfaces (UI) Command functions ...... 15 1.3 Typesetting conventions and scientific notations .16 1.4 Release notes for version 4.1 ...... 18 1.4.1 Key features ...... 18 1.4.2 Detail by function ...... 18 1.5 Release notes for version 4.0 ...... 21 1.5.1 Key features ...... 21 1.5.2 Detail by function ...... 21 1Preface 1.1 Getting started
The Structural Dynamics Toolbox provides solutions in many areas of structural dynamic modeling. The division shown below gives a good summary of the Toolbox areas. This section is intended for people who don’t want to read the manual. It summarizes what you should know before going through the SDT demos to really get started.
The SDT demonstrations are located in the sdtdemos directory which for a proper installation should be in your Matlab path. Executing demosdt at the Matlab prompt will also add the demo directory to your path if needed. Many of these demonstrations are associated to manual pages. You can easily access the proper page with your favorite web browser by typing the doc commands listed in the demos at the Matlab prompt. The series of gart.. demos covers a great part of the typical uses of the SDT.These demos are based on the test article used by the GARTEUR Structures & Materi- als Action Group 19 which organized a Round Robin exercise where 12 European laboratories tested a single structure between 1995 and 1997.
8 gartfe builds the finite element model using the femesh pre-processor gartte shows how to prepare the visualization of test results and perform basic correlation gartid does the identification on a real data set gartco shows how to use fe sens and fe exp to perform modeshape expansion and more advanced correlation gartup shows how the upcom interface can be used to further correlate/update the model
Area 1: Frequency domain identification
Starting with transfer functions, the identification phase of a modal analysis seeks to determine system poles and associated modeshapes. The SDT introduces a unique model tuning procedure to gradually build very accurate models. Key theoreti- cal notions are pole/residue models (page 36), residual terms (page 105), and the relation between residues and modeshapes (page 34). The iiplot/idcom graphical user interface is designed to help you going through the identification procedure. The d iiplot demo and the associated manual section (page 48) walk you through basic capabilities of the interface. The gartid script gives real data and an identification result for the GARTEUR example. The demo id script analyses a simple identification example. Associated details are given in section 3.3 (identification in the pole/residue format) and sec- tion 3.4 (transformations from the pole/residue model to other model formats). d pre describes basic of wire frame display of identification results.
9 1Preface
Area 2: Experimental modal analysis
Experimental modal analysis combines techniques related to system identification (data acquisition and signal processing, done outside the SDT, followed paramet- ric identification as detailed above) with information about the spatial position of multiple sensors/actuators. One of the first steps of a modal test is to build an experimental geometry allowing the visualization of the actual motion of structural response at the various sensors. The steps needed to visualize experimental deformations in the SDT are described in section 4.1 and demonstrated in gartte and d pre. Correlation between test results and finite element predictions is typically the second step. Usual correlation criteria and their implementation are discussed in section 4.3 anddemonstratedinthegartte and gartco. Direct uses of experimentally derived models are finally discussed in section 4.4.
Area 3: Basic finite element analysis
The .mex files introduced in version 3.1 give good performance for models with thousands of elements and tens of thousands of DOFs. While the performance will still increase over time, flexibility and performance are often conflicting objectives and the SDT will always favor flexibility and modularity. The library of elements (see chapter 8) is also steadily improving but will not, for a while longer, match the quality of element libraries found in major commercial finite element codes (NASTRAN, ANSYS, . . . ). See chapter 8 for the SDT policy statement on this subject. To manipulate finite element models, you will probably want to start by getting a working understanding of the femesh user interface by taking a look at section 5.1 and d truss, d plate, d ubeam, gartfe, ... You will then want to learn how to assemble finite element models and make basic predictions (static and dynamic re- sponse, . . . ) by reading section 5.2 and running demo fe. You will finally want to learn more about the graphical capabilities of the SDT for the visualization of deformed structures by reading section 5.3.
10 Area 4: Advanced FE analysis (model reduction, component mode synthesis, families of models)
Advanced model reduction methods are one of the key applications of the SDT.To learn more about model reduction in structural dynamics read section 6.1.Tosee how the SDT gives you simple access to all traditional component mode synthesis methods see section 6.2.1 and d cms. Reduction methods are often applied to components. Easy manipulation of com- ponent models is enabled by superelements. section 6.2.2 and d cms2 show how SDT superelements can be used for the traditional application of coupled prediction with reduction of component models. The superelements supported by fesuper and fe super are however more general in particular because of their ability to handle parametrized models. The SDT supports many tools necessary for finite element model updating. The upcom user interface function allows flexible parametrization of large finite element models (5000 element, 50000 DOF range). The fe sens (see section 4.1.2)and fe exp (see section 4.2) functions provide methods for easily correlating industrial size models. The updating strategy is however not stabilized and optimization rou- tines implemented with the up freq, up ifreq, up ixf functions (see section 7.2) are not very robust.
1.2 Understanding the Toolbox architecture
1.2.1 Layers of code
The SDT has three layers of code.
• Graphical user interfaces (feplot, iiplot, ii mac) provide a layer of pre- defined operations for Frequency Response Function (FRF) visualization and analysis, identification, 3-D deformation animation, and test/analysis correla- tion. Graphically supported operations (interactions between the user and plots/ menus/mouse movements/key pressed) form a subset of commands provided by user interface functions. • User interface (UI) functions provide high level solutions to problems in iden- tification, finite element mesh handling, model reduction, sensor placement, su- perelement handling or parametrized models for FE model update. The first argument to these functions is a string command which is parsed to know what
11 1Preface
operations to perform. See commode for conventions linked to parsed commands. • Scientific functions implement standard and state of the art methods in exper- imental modal analysis, Finite Element analysis, and to some extent in structural design and FE model update. These functions are open and can be easily extended to suit particular needs using the scientific environment provided by Matlab.
1.2.2 Graphical User Interfaces (GUI)
Graphical User Interfaces (GUI) in the Structural Dynamics Toolbox are designed to share a number of features in particular to allow inter-operability. The GUI layer
• generates/refreshes standard plots using plot functions: feplot (for deformed structures), idiplot (for identification results) iiplot (for measured/predicted time/frequency responses), ii mac (for test/analysis correlation). • allows modifications of plots through menus, buttons, mouse and key movements supported by related user interface (UI) functions (fecom, feplot, iicom, idcom, iimouse, ii mac).
The policy of the GUI layer is to let the user free to perform his own operations at any point. Significant efforts are made to ensure that this does not conflict with the continued use of GUI functions, but it is accepted that it may exceptionally do so. In most such cases, clearing the figure (using clf) or in the worst case closing it (use close or delete) and replotting (feplot or iiplot) will solve the problem.
Command and plot functions
Plot functions (iiplot, idiplot, feplot, ii mac) create and maintain plots. Each plot is contained in a Matlab axes object, initialized with the iicom sub commands and called a drawing axis. All plot functions can handle arbitrary numbers of axes within the same figure (window) as well as several figures. Plot functions use global variables and/or graphical object data to update the axes linked to a particular interface. Command functions parse text commands to provide the real access to function- ality of the user interface. The unified format used by command functions (see commode for details) allows command chaining, command mode and command script- ing through the general purpose parser commode. idcom and iicom are the two command functions using the iiplot interface. fecom and feplot can be used interchangably for feplot related commands. ii mac is
12 used to modify its own plots. iiplot serves as a common entry point to update plots of all types and handles plots linked to the visualization of test data, idiplot plots linked to identification, feplot plots of deformed structures. iimouse allows a number of mouse operations in the drawing axes and is designed to work with all Matlab generated axes.
Graphical inputs
Graphical objects (menus and buttons) provide access to functionality implemented in the command func- tions. To allow inter-operability of all the interfaces, these objects are all created and maintained by comgui. The functionality of different objects are described un- der command functions iicom, idcom, fecom.
The menus and buttons (controls) can be placed in the same figure as the drawing axes (iicom GuiOne command) or in a different one (iicom GuiTwo command). You can open a number of feplot and iiplot figures. The IIplot:Command Figure and FEplot:Command Figure menus let you swich between the various configurations. Context menus are associated with various objects and opened by cliking on the object witht the right mouse button. These menus are used to give easier access to functionality associated with each object. Starting with SDT 4 you can no longer mix feplot and iiplot drawing axes but you can operate both interfaces simultaneously and open more than one figure of each type. iimouse allows operations on plots using mouse (zoom and information about par- ticular objects or points in the plot, etc.) and keyboard inputs (zooming, rotation, next and previous channel selection, etc.).
Data handling
User interfaces require a knowledge of the current state of the interface and appro- priate data. The policy of the Toolbox is to let the user free to perform his own operations at any point. Significant efforts are made to ensure that this does not conflict with the continued use of GUI functions, but it is accepted that it may exceptionally do so. This flexibility resulted in the use of both global variables
13 1Preface
(for information that the user is likely to modify) and graphical objects (for other information). The user interface for data visualization and identification (iicom, idcom, iiplot, idiplot) uses a number of standard global variables shown below
Frequency response data XF standard data base wrapper pointing to the global variables with the data and storing the characteristics of FRF data sets (see xfopt for details) IIw vector of frequency points (same as XF(1).w) IIxf MIMO set of measured FRF data (XF(1).xf) IIxe identified data set (XF(2).xf) IIxh, IIxi other sets of FRF data (XF(3).xf and XF(4).xf) IDopt identification options for the current model (see idopt,sameas XF(1).idopt) XFdof global variable storing options describing each column of XF(1).xf
Identified model data IDopt identification options for the current model (same variable than for FRF data) IIpo set of poles of the main identified model (see ii pof,sameasXF(5).po) IIpo1 set of poles for the alternate identified model (XF(6).po) IIres residues of the main identified model (see idcom,sameasXF(5).res) IIres1 residues of alternate identified model (see idcom,sameasXF(5).res)
The femesh user interface for finite element mesh handling uses a number of standard global variables shown below
FEnode main set of nodes (also used by feplot) FEn0 selected set of nodes FEn1 alternate set of nodes FEelt main finite element model description matrix FEel0 selected finite element model description matrix FEel1 alternate finite element model description matrix By default, femesh and iiplot automatically use base worspace definitions of the standard global variables: base workspace variables with the correct name are trans- formed to global variables even if you did not dot it initially. When using the stan- dard global variables within functions, you should always declare them as global at
14 the begining of your function. If you don’t declare them as global modifications that you perform will not be taken into account, unless you call femesh, iiplot, ... from your function which will declare the variables as global there too. The only thing that you should avoid is to use clear and not clear global within a function and then reinitialize the variable to something non-zero. In such cases the global variable is used and a warning is passed.
1.2.3 User Interfaces (UI) Command functions
UI command functions provide a low number of access points to numerous com- plex operations. They use text string commands to specify what should be done. The commode function provides a general purpose command parser for all the UI functions (take a look at the associated help for conventions in giving the help on possible commands). commode allows command chaining, scripting, etc. The SDT currently provides the following UI command functions fecom modifications of properties of deformed structure plots generated by feplot. Assumes initialization using feplot. femesh finite element meshing utilities (FE pre-processor). Calls feplot, fecom to visualize the results. fesuper manipulations of superelements which can be combined with or used instead of element functions. fe reduc computations of vector basis used for model reduction and component mode synthesis. fe sens sensor configuration declaration and sensor/actuator placement tools fe stres energy and stress computations (limited capabilities) fe load creation of distributed load patterns (limited capabilities) idcom standard procedures for frequency domain identification. Lets you per- form the steps of an identification using short text commands and over- lays model/test comparisons when possible (iicom commands can be used to select the comparisons shown). iicom modifications of properties of standard plots generated by iiplot.Also supports general purpose commands use by all plot functions (sub, show, ...). upcom standard computations and manipulations of the Up superelement used for finite element update and optimization problems.
15 1Preface 1.3 Typesetting conventions and scientific notations
The following typesetting conventions are used in this manual
courier commands, function names, variables Italics Matlab Toolbox names, mathematical notations, and new terms when they are defined Bold key names, menu names and items Small print comments Conventions used to specify string commands used by user interface functions are detailed under commode. The following conventions are used to indicate elements of a matrix
(1,2) the element of indices 1, 2 of a matrix (1,:) the first row of a matrix (1,3: ) elements 3 to whatever is consistent of the first row of a matrix Usual abbreviations are
CMS Component Mode Synthesis (see section 6.2.1) COMAC Coordinate Modal Assurance Criterion (see ii mac) DOF,DOFs degree(s) of freedom (see section 2.13) FE finite element MACModal Assurance Criterion (see ii mac) MMIF Multivariate Mode Indicator Function (see ii mmif) POCPseudo-orthogonality check (see ii mac) For mathematical notations, an effort was made to comply with the notations of the International Modal Analysis Conference (IMAC) which can be found in Ref. [1]. In particular one has
16 [],{} matrix, vector ¯ conjugate [b] input shape matrix for model with N DOFs and NA inputs (see T T th section 2.1). φj b , ψj b modal input matrix of the j normal / complex mode [c] sensor output shape matrix, model with N DOFs and NS out- th puts (see section 2.1). {cφj} , {cψj} modal output matrix of the j normal / complex mode
[E]NS×NA correction matrix for high frequency modes (see section 2.6) [F ]NS×NA correction matrix for low frequency modes (see section 2.6) M,C,K mass, damping and stiffness matrices N,NM numbers of degrees of freedom, modes NS,NA numbers of sensors, actuators { } p NM×1 principal coordinate (degree of freedom of a normal mode model) (see section 2.2) { } q N×1 degree of freedom of a finite element model s Laplace variable (s = iω for the Fourier transform) { } T th [Rj]=cψj ψj b residue matrix of the j complex mode (see sec- tion 2.6) { } T th [Tj]=cφj φj b residue matrix of the j normal mode (used for proportionally damped models) (see section 2.6) { } u(s) NA×1 inputs (coefficients describing the time/frequency content of applied forces) { } y(s) NS×1 outputs (measurements, displacements, strains, stresses, etc.) [Z(s)] dynamic stiffness matrix (equal to Ms2 + Cs + K ) [α(s)] dynamic compliance matrix (force to displacement transfer func- tion) p, α design parameters of a FE model (see section 7.1.1) ∆M,∆C, ∆K additive modifications of the mass, damping and stiffness matrices (see section 7.1.1) [Γ] non-diagonal modal damping matrix (see section 2.3) λj complex pole (see section 2.5) ≤ [ φ]N×NM real or normal modes of the undamped system(NM N) \ 2 Ω \ modal stiffness (diagonal matrix of modal frequencies squared) ma- trices (see section 2.2)
[θ]N×NM NM complex modes of a first order symmetric structural model (see section 2.5)
[ψ]N×NM NM complex modes of damped structural model (see section 2.5)
17 1Preface 1.4 Release notes for version 4.1
1.4.1 Key features
Key features of this release are
• The new idcom Graphical User Interface gives much easier access to identification methods. • Many new context menus and enhancements to the SDT Handle objects to allow easier user manipulations of plots generated by iiplot and feplot. • Extensions to the femesh FindElt and FindNode commands allow a uniform treatement FE model part selection for display feplot, test/analysis correlation fe sens, energy computations fe stres, modification parametrization upcom. • Transformations to/from state-space models now support LTI objects of the Con- trol Toolbox for easier incorporation into SIMULINK.
SDT 4.1 is developed for Matlab 6.0 and later. It is compatible with Matlab 5.3.1 but not with earlier versions of Matlab. Some fe sens commands are not compatible with earlier versions of the same func- tion.
1.4.2 Detail by function
basis was extended to support local coordinate system definitions for nodes and degrees of freedom. comstr this string utility function is now provided as compiled function. This allows significant speed enhancement in certain GUI and IO operations. femesh Element selection was significantly extended. In particular, it now supports complex node declarations. For example femesh(’selelt group4 & innode {y>-230 & y<-215}’). Node selection was ex- tended and now supports complex element declarations. For example femesh(’FindNode y>-230 & inElt {egid -1}’). Problems in global variable checking when calling femesh from functions were fixed. Problems with unequal division of volume elements were fixed.
18 feplot See dfeplot for an introduction to the main current features. The notion of element selection has been fully revised to use femesh FindElt syntax. The text commands now accept node selection based on femesh FindNode syntax. TextNode commands automatically search element selection if command results in no displayed node. For wire frame meshes, problems when showing isolated points were fixed. MatID is set to be the line number so that fecom(’colordatamat’) generates colored wire frame displays. ProID is used to store color information if defined. Significant enhancements to menus and context menus are introduced. Buttons are hidden for copying to clipboard. Support for local coordinate system definitions for nodes and degrees of freedom is introduced. The new show test, fem, links, 2def commands allow easier creation of standard plots. Bugs with wire-frame animation of tests with arbitrary directions (5 column format) were fixed. The new AnimMoviei command returns an i step Matlab movie. feutil now supports EltId, MatId, ProId extraction, conversion of deforma- tions to/from mode form, ... fe sens was fully revised with a new strategy to for topological correlation tools. Rotational sensors are now supported. The result is not compatible with the previous version of fe sens. idcom now comes with a new graphical user interface which gives much easier access to all standard commands. The new estlocal commands makes building models by band easier. ii mac supports improved table output with tabs for easier incorporation into Word and tex output for incorporation in LaTeX. Some modifications in the data structure of vector correlation objects. iiplot has many new context menus and supports new standard plots (sum of imaginary parts, ...) mass2 is a new concentrated mass element with offset support and more.
19 1Preface
nasread was extented to support EltID as well as additional elements. More consistent support of output2 reading including importing of element matrix dictionaries into superelements used by upcom (see nas2up). Sig- nificant speed improvement are obtained by calling mex functions. nas2up is a new function, offered in the FEMLINK extension to the SDT, al- lowing direct reading of model definition and element matrices from NASTRAN output2 obtained when using PARAM,POST,-4 to type 3 superelements supported by upcom. nor2ss was extended to support acceleration input and direct ouput to LTI objects of the Control Toolbox for easier incorporation into SIMULINK. Support for direct input from FEM modes is also provided. nor2xf supports direct input from FEM modes and allows direct output to database wrapper format with all options set. quad4 was corrected to solve problems with the drilling DOF stiffness. This removes numerical conditioning problems previously found for tilted flat elements. tria3 was corrected to solve problems with the drilling DOF stiffness. This removes numerical conditioning problems previously found for tilted flat elements. ufread Compilation of file start search speeds reading significantly. Color in- formation for wire frame displays is now properly imported. For UFF 58, delayed data loading and block separation are now handled more consistently. upcom the whole function has been revised to allow the simultaneous use of multiple type 3 superelements. The save command is much faster. For energy display with fecom, you should now use the eners command as demonstrated in gartup xfopt problems with option check and assignation for database wrappers with global variables were fixed.
20 1.5 Release notes for version 4.0
1.5.1 Key features
Key features of this release are
• A major revision of the Graphical User Interfaces, with a complete rewritting of iiplot/idiplot (FRF visualization), feplot (deformed structure) and the introduction of a GUI for vector correlation tools (ii mac extensions). • Introduction of SDT Handle objects which give easier and more mnemonic access to complex functionality. Current SDT Handle types are feplot figures, IDopt options, iiplot figures, vector correlation objects (see ii mac), and database wrappers (see xfopt). • A reprint of the paper manual. A new PDF version of the manual. Enhanced HTML documentation with more links, compatibility with the Matlab doc com- mand, direct access to tutorial sections using doc, listing of related sections when using the Matlab help.
SDT 4.0 is developped for Matlab 5.3 and later. It is not compatible with Matlab 4.2c, 5.0 and 5.1. There are known graphical incompatibilities with Matlab 5.2.1.
1.5.2 Detail by function bar1 Fixed a problem with the mass that was only considered in the axial direction. elem0 Since feplot no longer uses line definitions, element functions no longer need to declare a line connectivity. Negative element group identifiers are used for element groups that are for display only. fe2ss Fully revised to combine modeshape and static correction computations in the same call. Now supports free or fixed interface modes. Backward compatibility will disappear in future releases. femesh A number of minor modifications, extensions and bug fixes. The find commands now accept an arbitrary number of numeric arguments. El- ement selection supports display only element groups for test/analysis correlation.
21 1Preface
feplot Fully revised. Now supports SDT handles, multiple figures, enhanced mouse operation, new scaling options, easier manipulation of plots with multiple objects, ... The Feplot menu is fully revised and gives direct access to many more operations. The figure is associated to a context menu (open with the right mouse button). Information on nodes is now provided in the cursor mode (click on an axis and type the c key to turn this mode on/off). Animation supports multiple axes and color animation. All plots are now based on patches so that element line definitions are no longer used. A special treatement of true segments (beams) and isolated nodes (concentrated masses) is done. Variable number of nodes are now supported for real faces (not just 3 or 4). Text generation commands are fully revised. Text is placed on defor- mation. Directly uses Matlab 3-D camera based graphic properties which pro- vides significant speedup in particular for object rotations and transla- tions. Read the function help for more details. fecom Fully revised (see details under feplot) fegui No longer needed (global variables declared in workspace when calling femesh). fe exp Significant documentation and optimisation efforts. Now fully compati- ble with sensor configurations defined with fe sens. Minimum residual dynamic expansion with and without measurement error is supported. fe mk Now properly handles matrices with more than 46340 (sqrt(2^31)) DOFs. Ignores element groups that are for display only. Optimized speed. fe reduc Direct support for Craig-Bampton and free interface model reduction methods. fe super Fixed problem linked to element DOF assembly for unique superele- ments idopt Fully revised. Options are stored in a SDT Handle which gives more user friendly access to the options. iigui No longer needed (global variables are declared in workspace when call- ing idcom,oriiplot). iiplot Fully revised. Now supports SDT handles, multiple figures, enhanced Iiplot menu and mouse operation (in particular cursor and context menu), contextmenu/mouse editing of title/labels, automatic computa- tion of the CMIF, ... Supports multiple data sets with non matching abcissa values. 22 ii mac Introduction of a interactive visualization tool (view rotation, mouse cursor to display values, ...). Menu based selection of criterion and reference data set. Supports MAC, MAC-M, POC, relative error, MACCO, COMAC, eCO- MAC, ... nasread Major revision. Improved reading of NASTRAN bulk input files (major speed improvements, coordinate transformations, include cards, more elements, ...). Better support for OUTPUT2 and OUTPUT4 binary formats. skyline Minor speed and robustness enhancements. ufread, Completely rewritten to take advantage of database wrappers. This ufwrite allows easier access to various datasets and handles to data set types. upcom Significant optimization and robustness enhancements. Compiled ele- ment energy computations. xfopt Now used to manipulate of database wrapper objects which introduce much more consistence and flexibility in the handling of options describ- ing data sets.
23 1Preface
24 2
Fundamental concepts
2.1 I/O shape matrices ...... 26 2.2 Normal mode models ...... 28 2.3 Damping ...... 30 2.3.1 Viscous damping in the normal mode model form .31 2.3.2 Damping in finite element models ...... 32 2.4 State space models ...... 32 2.5 Complex mode models ...... 34 2.6 Pole/residue models ...... 36 2.7 Parametric transfer function ...... 37 2.8 Non-parametric transfer function ...... 38 2.9 Nodes ...... 39 2.10 Finite element model description matrices .... 40 2.11 Material properties ...... 42 2.12 Element property matrices ...... 43 2.13 DOF definition vector ...... 44 2.14 DOF selection ...... 45 2 Fundamental concepts
This chapter is intended as a reference for the fundamental notions and associated variables used throughout the Toolbox. This information is grouped here and hy- pertext reference is given in the HTML version of the manual. Models of dynamic systems are used for identification phases and links with control applications supported by other Matlab toolboxes and Simulink. Key concepts and variables are
b,c (p. 26) input/output shape matrices (b,c,pb,cp variables) nor (p. 28) normal mode models (freq,damp,cp,pb variables) damp (p. 30) damping for full and reduced models cpx (p. 34) complex mode models (lambda, psi variables) res (p. 36) pole/residue model (res,po variables) ss (p. 32) state space model (a,b,c,d variables) tf (p. 37) parametric transfer function (num,den variables) xf (p. 38) non-parametric transfer function (w,xf variables)
The description of geometrical and material properties of structures are used for finite element modeling and to process experimental modal tests of structures. Key concepts and variables are
node (p. 39) node matrices elt (p. 40) model description matrices containing the elements pl (p. 42) material property matrices il (p. 43) element property matrices mdof (p. 44) DOF definition vector adof (p. 45) active DOF definition vector and wild cards tr (p. 103) model reduction for structures
2.1 I/O shape matrices
Dynamic loads applied to a discretized mechanical model can be decomposed into a { } { } product F q =[b] u(t) where
• the input shape matrix [b] is time invariant and characterizes spatial properties of the applied forces and • the vector of inputs {u} allows the description of the time/frequency properties.
26 Similarly it is assumed that the outputs {y} (displacements but also strains, stresses, etc.) are linearly related to the model coordinates {q} through the sensor output shape matrix ({y} =[c] {q}). Input and output shape matrices are typically generated with fe c or fe load. Understanding what they represent and how they are transformed when model DOFs/states are changed is essential for a good use of the Toolbox.
Mechanical models considered in the SDT take the general forms 2 { } { } Ms + Cs + K N×N q(s) =[b]N×NA u(s) NA×1 { } { } (2.1) y(s) NS×1 =[c]NS×N q(s) N×1 in the frequency domain (note Z(s)=Ms2 + Cs + K), and [M] {q¨} +[C] {q˙} +[K] {q} =[b] {u(t)} (2.2) {y(t)} =[c] {q(t)} in the time domain. In the model form (2.1), the first set of equations the describe the evolution of {q}. The components of q are called Degrees Of Freedom (DOFs) by mechanical engineers and states in control theory. The second observation equation is rarely considered by mechanical engineers (hopefully the SDT may change this). The purpose of this distinction is to lead to the block diagram representation of the structural dynamics
{u(s)} {F (s)} {q(s)} {y(s)} ✲ ✲ 2 −1 ✲ ✲ [b] Ms + Cs + K [c] which is very useful for applications in both control and mechanics.
In the simplest case of a point force input at a DOF ql, the input shape matrix is equal to zero except for DOF l where it takes the value 1 . . 0 [bl]= 1 ← (2.3) l 0 . .
27 2 Fundamental concepts
T Since [bl] {q} = {ql}, the transpose this Boolean matrix is often called a localization matrix. Boolean input/output shape matrices are easily generated by fe c (see the section on DOF selection page 45). Input/output shape matrices become really useful when not Boolean. For applica- tions considered in the SDT they are key to
• test analysis correlation. Since you often have measurements that do not directly correspond to DOFs (accelerations in non global directions at positions that do not correspond to finite element nodes, see section 4.1.2). • model reduction. To allow the changes to the DOFs q while retaining the physical meaning of the I/O relation between {u} and {y} (see page 103).
2.2 Normal mode models
The spectral decomposition is a key notion for the resolution of linear differential equations and the characterization of system dynamics. Predictions of the vibrations of structures are typically done for linear elastic structures or, for non-linear cases, refer to an underlying tangent elastic model. Spectral decomposition applied to elastic structures leads to modal analysis.The main objective is to correctly represent low frequency dynamics by a low order model whose size is typically orders of magnitude smaller than that of the finite element model of an industrial structure. The use of normal modes defined by the spectral decomposition of the elastic model and corrections (to account for the restricted frequency range of the model) is funda- mental in modal analysis. Associated models are used in the normal mode model format