e_comac_xml


Purpose

Automated report of differences between test and analysis shapes using e-comac equation

Syntax

EC = e_comac_xml(st,sg,[TAM or CT], [MINVALUE or SHAPELIST])

Description

E_COMAC_XML Automated reporting shape differences between test and analysis shapes using the e_comac equation

This routine can will generate enhanced comac results and place them into an xml file that can be read by Microsoft Excel. The user also has the option of plotting the shapes that are created by this routine.

E_COMAC_XML compares test shapes, ST, and FEM shapes, SG. If TAM is available, then the comparison is based on the test-analysis cross-orthogonality matrix. If no TAM is supplied, then the test-FEM cross-MAC is used based on the supplied coordinate trace, CT. If CT is not supplied, then all non-zero entries in both the FEM and test shapes are used.

Modes that are compared are either all combinations of test and FEM modes that produce a cross-ortho (or cross-MAC) above MINVALUE. If SHAPELIST is supplied, then those shapes are compared directly.

Inputs:

st

imat_shp

sg

imat_shp

TAM

TAM matrices. Must contain TAM.aset, TAM.mass, and TAM.stiffness
If TAM.aset is not supplied, then it is assumed that TAM.aset = alldof(st), and TAM.mass = I.
If TAM.mass is not supplied, then it is assumed that TAM.mass = I
If TAM.stiffness is not supplied, then SE_ASET, and SE_GROUP will not be generated.

ct

imat_ctrace. All coordinates in ct should be in both st and sg.

minvalue

Minimum value for cross-ortho or cross-MAC to compare two sets of shapes.

shapelist

A N by 2 matrix containing a list of test and FEM shapes to compare. The first column corresponds to the test shapes, ST, and the second column corresponds to the FEM shapes, SG.

Outputs:

EC.aset

DOF used in e_comac

EC.shp

E_COMAC results in shape function format

EC.summary

NASET by 1 vector containing summary results for E_COMAC
EC.summary = 1 - [sum(j=1:l)(abs(st(j,:)-sg(j,:)))]/[[sum(j=1:l)max(st(j)][sum(j=1:l)max(sg(j))] ]j corresponds to the mode comparison, and l = number of modes in comparison (NMODES)

EC.matrix (i,j)

NASET by NMODES
i corresponds to EC.aset(i)
j corresponds to jth mode comparison.
EC.matrix (i,j)= (sign(st).*sign(sg)).*abs((st(i,j)-sg(i,j)))/(max(|st(j)|) + max(|sg(j)|));
EC.matrix calculates the e comac between two modes with the following exceptions:
1. A sign flip calculation is put into the equation so that the results will be negative if the two mode shape coefficients have opposite sign.
2. The results are not subtracted from 1. Higher values signify mode shape errors

EC.stfreq (j)

jth test mode frequency

EC.sgfreq (j)

jth FEM mode frequency

Examples

>> EC=e_comac(st,sg,TAM,97)

See Also

ortho, comparison_table, imat_shp, report_test_xml