Index of /reinerwt/tonguedata

Icon  Name                    Last modified      Size  Description
[DIR] Parent Directory                             -   
[   ] FEM3x3Gausspts.m        07-Jul-2005 11:22  338   
[   ] FEM3x3Lobatto.m         07-Jul-2005 11:22  331   
[TXT] FEM4x4Gausspts.m        07-Jul-2005 11:22  964   
[TXT] FEMDbar.m               07-Jul-2005 11:22  856   
[TXT] FEMFields.m             07-Jul-2005 11:22  2.7K  
[TXT] FEMForces.m             07-Jul-2005 11:22  2.1K  
[TXT] FEMGauss3on27.m         07-Jul-2005 11:22  1.8K  
[TXT] FEMGrad.m               07-Jul-2005 11:22  4.7K  
[TXT] FEMKonstructor.m        07-Jul-2005 11:23  120K  
[TXT] FEMMatrices.m           07-Jul-2005 11:23  1.1K  
[TXT] FEMReCalc.m             07-Jul-2005 11:23  728   
[TXT] FEMgradoperator.m       07-Jul-2005 11:23  4.4K  
[TXT] FEMinterpol.m           07-Jul-2005 11:23  1.3K  
[   ] FEMintoperator.m        07-Jul-2005 11:23  312   
[   ] FEMtestK.m              07-Jul-2005 11:23  1.6K  
[   ] InitFEM.m               07-Jul-2005 11:23  1.1K  
[   ] Q27index.m              07-Jul-2005 11:23  452   
[   ] RunFEM.m                07-Jul-2005 11:23  779   
[TXT] Total.FEM               07-Jul-2005 11:24  795K  
[TXT] fourthindicies.m        07-Jul-2005 11:23  1.2K  
[TXT] lagrangeindicies.m      07-Jul-2005 11:23  1.3K  
[TXT] readFEMfile.m           07-Jul-2005 11:23  3.5K  
[   ] rewttongue.zip          07-Jul-2005 11:33  339K  
[   ] rewttonguemodel.gz      07-Jul-2005 11:33  319K  
[TXT] splitaname.m            07-Jul-2005 11:23  819   

Added July 2005. Author: Reiner Wilhelms-Tricarico (ReWT).



These files allow loading the data structures of the tongue model, using
Matlab. I used version 7.0 or R14.

In a matlab window, run

InitFEM 

It loads the file Total.FEM, which contains the data defining the tongue
model. It then generates various structures that used to all make sense to me
;-)

The structure that is generated is in the variable (global ) Fe. It contains
contains node points, element structures, and fiber fields for the 12 muscles
of the model. The fiber indices refere to the Gauss points in each element. In
each element there are 27 Gauss points (Gauss Lobatto points). Each Gauss
point has a global index which is the sum of 27 times element-number plus
local index minus one.  (the local index is a number between 1 and 27, see
Q27index.m).


To see what's in the data structure Fe, type Fe, and the following structures
are listed. For some I don't remember what they are for.


Fe =

             xyz: [2487x3 double]        % node coordinates
          constr: [2487x3 double]        % constraints on nodes (indices)
            elem: [250x27 double]        % element assembly
     musclenames: {1x12 cell}            
       fiberindx: {1x12 cell}
        fiberdir: {1x12 cell}
    fiberdensity: {1x12 cell}
              Nk: [6750x6750 double]     % Interpolation function at Gauss pts.
              Ax: [6750x6750 double]     % gradient x of Interpol. functions
              Ay: [6750x6750 double]     %          y 
              Az: [6750x6750 double]     %          z
             Ass: [2487x6750 double]     % assembly matrix
           nelem: 250                      
          ngauss: 6750
          nnodes: 2487
               W: [6750x6750 double]     % weights
              Wv: [6750x1 double]        % forgot what it is
              J0: [6750x6750 double]     % det(Jacobians)
             Nk1: [6750x2487 double]     % transformed gradients of interpol. 
             Nk2: [6750x2487 double]     % (1,2,3) corresponding to x,y,z
             Nk3: [6750x2487 double]     % 
            Mass: [2487x2487 double]     % mass matrix
               U: [2487x3 double]
           lamda: 164.4295               % const in Neohookean model?
              mu: 3.3557                 % 
           kappa: 10                     % forget about this (was bulk compressibility ).