


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

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