Topology correlation and test preparation

Contents Functions

3.1 Topology correlation and test preparation

Topology correlation is the phase where one correlates test and model geometrical and sensor/shaker configurations. Most of this effort is handled by fe_sens with some use of femesh.

Starting with SDT 6.0, FEM sensors (see section 4.6) can be associated with wire frame model, the strategy where the two models where merged is thus obsolete.

As described in the following sections the three important phases of topology correlation are

combining test and FEM model including coordinate system definition for the test nodes if there is a coordinate system mismatch,
building of an observation matrix allowing the prediction of measurements based on FEM deformations,
sensor and shaker placement.

3.1.1 Defining sensors in the FEM model : data handling

Two types of data are needed to properly associate a test wire frame model to a FEM :

a FEM model (see section 4.5). For this simple example, the FEM model (stored in cf.mdl in the demo) must describe nodes, elements and DOFs
a test wire-frame model (stored in TEST in the demo) with sensors in the .tdof field, as detailed in section 2.7.1 for the geometry and section 2.7.2 for sensors

One then declares the wire frame (with sensors) as SensDof case entry as done below (see also the gartte demo). The objective of this declaration is to allow observation of the FEM response at sensors (see sensor Sens).

 model=demosdt('DemoGartFE') % load FEM
 TEST=demosdt('demo garttewire'); % see sdtweb('pre#presen')
 % Store Test as SensDof (linked test wireframe) in the FEM
 model=fe_case(model,'sensdof','sensors',TEST); 
 cf=feplot(2); cf.mdl=model; % Display the model in feplot
 % Display the superposition of the test wireframe over the FEM
 fecom(cf,'ShowFiCoTopo'); 
 % Open the CoShape Dock from cf, already containing needed data
 fecom(cf,'dockCoShape');

Section 4.6 gives many more details the sensor GUI : the available sensors (sensor trans, sensor triax, laser, ...). Section 4.6.4 discusses topology correlation variants in more details.

If the data come from files, it can be more convenient to load them directly from the GUI.

Here is a tutorial for interactive data loading in DockCoTopo with the TestBas tab.

You will need the garteur example files, which can be found in SDTPath/sdtdemos/gart*.m. If these files are not present, click on the first step on the following tutorial in the HTML version of the documentation or download the patch at the address https://www.sdtools.com/contrib/garteur.zip and unzip the content in the folder SDTPath/sdtdemos.

Execute the command fecom('dockCoTopo') to open an empty dock. You can also click on the button CoTopo on the tree in SDT Root.
Click on Select associated to MasterMesh. This will open the import model window. Select the file to load : for this tutorial, the file is located at SDTPath/sdtdemos/gart_mdl.inp. Data is loaded and displayed in the feplot figure.
Do the same for the SlaveMesh. The test mesh file is located at SDTPath/sdtdemos/gartid.unv. Data is loaded and displayed in the feplot figure. Once selected, the Unv tab is displayed in the feplot('mdl') figure : it shows the content of what is inside the Unv file.
Check the box corresponding to model and click on Import.

The test wireframe is loaded and displayed in the feplot figure in red.

Depending on the loaded data for the the SlaveMesh, it contains already or not the sensor definitions : they are shown as red arrows. It is not the case here.
To retrieve sensors definition from a Unv file, the mesured data need to be loaded.
Click on Select associated to DefineTDof. Select again the Unv file and in the Unv tab, check this time the box corresponding to response and click on Import.
The arrows are then built depending on the measured channels (+X,+Y,+Z,-X,-Y,-Z directions associated to each nodes in the geometry), and displayed.

The system coordinate is not the same between the test wireframe and the FEM : the test geometry needs to be moved and superposed to the FEM (this tutorial continues in the following subsection).

3.1.2 Test and FEM coordinate systems

In many practical applications, the coordinate systems for test and FEM differ. fe_sens supports the use of a local coordinate system for test nodes with the basis command.

Interactive test mesh placement is available in the SDT GUI, using command fe_sensGuiTestBas.

% Loading the interactive test mesh placement GUI
 cf=demosdt('demo garttebasis'); % Load the demo data
 cf.CStack{'sensors'} % contains a SensDof entry with sensors and wireframe
 fecom(cf,'setTestBas'); % Open interactive tab in feplot properties

Operations permitted through the GUI implementation are available in script commands. The modus operandi considers a three steps process.

Phase 1 is used get the two meshes oriented and coarsely aligned. The guess is more precise if a list of paired nodes on the FEM and TEST meshes can be provided.
In phase 2, the values displayed by fe_sens, in phase 1 are fine tuned to obtain the accurate alignment.
In phase 3, the local basis definition is eliminated thus giving a cf.CStack{'sensors'} entry with both .Node and .tdof fields in FEM coordinates which makes checks easier.

In peculiar cases, the FEM and TEST mesh axes differ, and a correction in rotation in the Phase 2 may be easier to use. An additional rotation to apply in the TEST mesh basis can be obtained by fulfilling the field rotation in Phase 2. The rotations are applied after other modifications so that the user can directly interpret the current feplot display. The rotation field corresponds to a basis rotate call. The command string corresponding to a rotation of 10 degrees along axis y is then 'ry=10;'. Several rotations can be combined: 'ry=10; rx=-5;' will thus first perform a rotation along y of 10 degrees and a rotation along x of -5 degrees. These combinations are left to the user's choice since rotation operations are not symmetric (e.g. 'rz=5;rx=10;' is a different call from 'rx=10;rz=5;').

The following example demonstrates the 3 phases in a script.

 cf=demosdt('demo garttebasis'); % Load the demo data
 cf.CStack{'sensors'} % contains a SensDof entry with sensors and wireframe
 
 % Phase 1: initial adjustments done once
  % if the sensors are well distributed over the whole structure
  fe_sens('basis estimate',cf,'sensors');

 % Phase 1: initial adjustments done once, when node pairs are given
  % if a list of paired nodes on the TEST and FEM can be provided
  % For help on 3DLinePick see sdtweb('3DLinePick')
  cf.sel='reset';    % Use 3DLinePick to select FEM ref nodes
  cf.sel='-sensors'; % Use 3DLinePick to select TEST ref
  i1=[62 47 33 39;          % Reference FEM NodeId
     2112 2012 2301 2303]';% Reference TEST NodeId
  cf.sel='reset'; % show the FEM part you seek
  fe_sens('basis estimate',cf,'sensors',i1);
 
 
 %Phase 2 save the commands in an executable form
 % The 'BasisEstimate' command displays these lines, you can
 % perform slight adjustments to improve the estimate
  fecom(cf,'initTestBas') % When you change a value script below displayed
  fe_sens('basis',cf,'sensors', ...
   'x',     [0 1 0], ... % x_test in FEM coordinates
   'y',     [0 0 1], ... % y_test in FEM coordinates
   'origin',[-1 0 -0.005],... % test origin in FEM coordinates
   'scale', [0.01]); % test/FEM length unit change 

 %Phase 3 : Force change of TEST.Node and TEST.tdof to FEM coordinates
 fecom('SetTestBas',struct('BasisToFEM','do'));
 fe_case(cf.mdl,'sensmatch')
 sens=fe_case(cf.mdl,'sens')

Note that FEM that use local coordinates for displacement are discussed in sensor trans.

Here is the continuation of the tutorial for interactive way to superpose and match sensors over the FEM.

If you have not performed previous tutorial (or if you closed everything at the end), click on this link in the HTML version of the documentation to get ready for the following.

To begin with, it is often useful, if the test geometry globally describes well the model geometry, to perform an automatic initial guess for the superposition. To so so, click on the button run associated to basEst.
From this better relative position, one needs to iterate manually with small translations tx, ty, tz and rotations rx, ry, rz until the optimum is reached.
Finally, click on the button Accept associated to BasisToFEM to apply the coordinate transformation to the test wireframe and perform the compute the observation matrix of the FEM at sensors.
Clicking on Finalize will save the result in the corresponding project.

Another strategy using Iterative Closest Point algorithm is also implemented (in the NodePairs subtable). This will be documented in further release.

3.1.3 Sensor/shaker placement

In cases where an analytical model of a structure is available before the modal test, it is good practice to use the model to design the sensor/shaker configuration.

Typical objectives for sensor placement are

Wire frame representations resulting from the placement should allow a good visualization of test results without expansion. Achieving this objective, enhances the ability of people doing the test to diagnose problems with the test, which is obviously very desirable.
seen at sensors, it is desirable that modes look different. This is measured by the condition number of [cφ]^T[cφ] (modeshape independence, see [14]) or by the magnitude of off-diagonal terms in the auto-MAC matrix (this measures orthogonality). Both independence and orthogonality are strongly related.
sensitivity of measured modeshape to a particular physical parameter (parameter visibility)

Sensor placement capabilities are accessed using the fe_sens function as illustrated in the d_cor('TutoSensPlace') demo. This function supports the effective independence [14] and maximum sequence algorithms which seek to provide good placement in terms of modeshape independence.

It is always good practice to verify the orthogonality of FEM modes at sensors using the auto-MAC (whose off-diagonal terms should typically be below 0.1)

 cphi = fe_c(mdof,sdof)*mode; ii_mac('cpa',cphi,'mac auto plot')

For shaker placement, you typically want to make sure that

you excite a set of target modes,
or will have a combination of simultaneous loads that excites a particular mode and not other nearby modes.

The placement based on the first objective is easily achieved looking at the minimum controllability, the second uses the Multivariate Mode Indicator function (see ii_mmif). Appropriate calls are illustrated in the d_cor('TutoSensPlace') demo.