Metafor

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== HowTo ======
==== General Consideration ====
Frequency Analysis are computed on a the whole model loaded though his ''Domain''.
That means that geometry, mesh, elements, material, fixations has to be defined.
Time integration, Iteration Manager, Solver Manager may not be defined.

Frequency Analysis may be conduced the same way :
  * on a domain not loaded
  * on a loaded domain 
  * during the load of the domain

==== FrequencyAnalysisMethod ====
Actually two Frequency Analysis Method are implemented in Metafor
  * a Power Iterations Frequency Analysis Method
  * a Lanczos Frequency Analysis Method (the most efficient one)
Common commands (according to the fact freqAnaly is a FrequencyAnalysisMethod : 
<code>
    freqAnaly.setNumberOfEigenValues(nbEigenVal)
    freqAnaly.setSpectralShifting(spectralShifting)
    freqAnaly.setSymmetrizeK(True)
    freqAnaly.setSolver(LinearSolver)
    freqAnaly.save2Fafac(filename)    
    freqAnaly.setWriteMatrix2Matlab(True)
    freqAnaly.setComputeErrors(True)
    freqAnaly.setPrec(True)
    freqAnaly.setMaxIte(True)
    freqAnaly.setVerbose(True)
    freqAnaly.execute()
    freqAnaly.setTSCLabel("EigenFrequency")
    freqAnaly.writeTSC()
    freqAnaly.getEigenValue(i)    
</code>

=== PowerIterationsFrequencyAnalysisMethod ===
Defining a PowerIterationsFrequencyAnalysisMethod allows to build the few first eigen values / eigen vectors (eigen vectors are intrinsiquely computed during orthogonalisation process).

<code>
    freqAnaly = PowerIterationsFrequencyAnalysisMethod(domain)    
</code>
=== LanczosFrequencyAnalysisMethod ===
The Lanczos Frequency Analysis Method benefit 
<code>
    freqAnaly = LanczosFrequencyAnalysisMethod(domain)
    freqAnaly.setInitByPowerMethod(False)
    freqAnaly.setGrahamSchmidt(true)
    freqAnaly.setLanczosSizeVSNbModesFactor(10)
    freqAnaly.setSymTriDiagEigenSolver(DSTE_VR)
</code>
==== Direct Use ====
As soon as a domain is defined (geometry, mesh, elements, materials, loadings), 
a frequency analysis can be defined and called using 
<code>
    freqAnaly.execute()
    freqAnaly.setTSCLabel("EigenFrequency")
    freqAnaly.writeTSC()
    freqAnaly.getEigenValue(i)    
    freqAnaly.save2Fafac(filename)    
</code>
==== FrequencyAnalysisValueExtractor ====
In some particular situation, it can be interesting to compute the eigen frequencies according to the load of the Domain. A good example is the measurement of the centrifugal stiffening...

To archieve this operation, a FrequencyAnalysisValueExtractor can be defined and values can be computed at each time step of a simulation : 
<code>
  freqAnalyValExt = FrequencyAnalysisValueExtractor(freqAnaly)
  freqAnalyValExt.setSaveFaFac(False)
</code>
(be care to the potential CPU cost of this operation as FrequencyAnalysis will be computed at each time step)

FrequencyAnalysisValueExtractor is a common ''ValueExtractor'', so values can be visualized during computation, used as TSC, ... refer to  [[doc:user:results:courbes_res]] and  [[doc:user:results:viz_courbes]] for more informations on how to use ''ValueExtractor''.

==== FrequencyAnalysisObjectiveFunction ====
In some other case, user may only want the eigen values at the end of a complex simulation. FrequencyAnalysisObjectiveFunction  is done for this
<code>
  freqAnalyFObj = FrequencyAnalysisObjectiveFunction(no, freqAnaly)
  freqAnalyFObj.saveToFaFac(False)
</code>
refer to [[doc:user:results:fobj|Objective functions]] for more information on how to use Objective Functions

==== Vizualisation of eigen modes ====

<code>
  try:  
        win = VizWin()
        for i in range(domain.getInteractionSet().size()):
            win.add(domain.getInteractionSet().getInteraction(i))
        for i in range(nbEigenVal):
            freqAnaly.showEigenVector(i)
            win.update()
            print( "Eigen Vector ", i , "EigenValue = ", freqAnaly.getEigenValue(i) )
            input("press enter to continue")
    except NameError:
        pass
</code>