PentaCollapseIntoTriaVector

PURPOSE ^

PENTACOLLAPSEINTOTRIAVECTOR - collapse a pentaelem into a tria and create a load vector for ths element

SYNOPSIS ^

function Pe=PentaCollapseIntoTriaVector(pentaelem,grids,materials,inputs,analysis_type);

DESCRIPTION ^

PENTACOLLAPSEINTOTRIAVECTOR - collapse a pentaelem into a tria and create a load vector for ths element

   This routine creates a fake triaelem from the pentaelem and then call the routine
   to create a load vector fot this triaelem.
   It is usefull when running MacAyeal's model on a 3d mesh.

   Usage:
      Pe=PentaCollapseIntoTriaVector(pentaelem,grids,materials,inputs,analysis_type)

   See also PENTACOLLAPSEINTOTRIAMATRIX, CREATEPVECTOR

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function Pe=PentaCollapseIntoTriaVector(pentaelem,grids,materials,inputs,analysis_type);
0002 %PENTACOLLAPSEINTOTRIAVECTOR - collapse a pentaelem into a tria and create a load vector for ths element
0003 %
0004 %   This routine creates a fake triaelem from the pentaelem and then call the routine
0005 %   to create a load vector fot this triaelem.
0006 %   It is usefull when running MacAyeal's model on a 3d mesh.
0007 %
0008 %   Usage:
0009 %      Pe=PentaCollapseIntoTriaVector(pentaelem,grids,materials,inputs,analysis_type)
0010 %
0011 %   See also PENTACOLLAPSEINTOTRIAMATRIX, CREATEPVECTOR
0012 
0013     %we collapse the penta into its base tria, and use the tria to build the stiffness matrix.
0014     element=triaelem;
0015     element.type='triaelem';
0016     element.id=NaN; %not needed here, we are going to destroy this triaelem very soon.
0017     element.matid=pentaelem.matid; %same materials, not dependent on element.
0018     element.g=pentaelem.g(1:3); %first three grids correspond to the base of the penta element.
0019     element.h=pentaelem.h(1:3);
0020     element.s=pentaelem.s(1:3);
0021     element.b=pentaelem.b(1:3);
0022     element.friction_type=pentaelem.friction_type;
0023     element.k=pentaelem.k(1:3);
0024     element.p=pentaelem.p;
0025     element.q=pentaelem.q;
0026     element.shelf=pentaelem.shelf;
0027     element.meanvel=pentaelem.meanvel;
0028     element.epsvel=pentaelem.epsvel;
0029     element.acceleration=pentaelem.acceleration;
0030 
0031     %Call CreateKMatrix for this new element.
0032     Pe=CreatePVector(element,grids,materials,inputs,analysis_type);
0033

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003