convert2d23d

This file converts a two-d coord system bilinear mesh ipnode file into a 3-d coordinate system ipnode file with a known uniform thickness.

#load perl libraries for CMISS

use CmUtils::File::Ipnode
use CmUtils::File::Ipdata
use CmUtils::Objects::NodeGroup
use CmUtils::Objects::Field
use CmUtils::Objects::FieldSet
use CmUtils::Objects::Node

use CmUtils::Vector


$thickness=0.3;
$convertFile="displacements_line_23";
$convertedFile="displacements3d_line_23";

fem def para;r;multiaxes

$group2d = readIpnode($convertFile);
$group3d = CmUtils::Objects::NodeGroup->new();
$group3d->name($convertedFile);
$numnodes2d = $group2d->numberOfNodes();
$numnodes3d = $numnodes2d*2;

for($i=1;$i<=$numnodes3d;$i++)
{

	$node3d = CmUtils::Objects::Node->new();
        if($i<=$numnodes2d)
        {
        	$node2d = $group2d->getNode($i);
        }
        else
        {
        	$node2d = $group2d->getNode($i-$numnodes2d);
        }
        $node2dx = $node2d->value(0);
        $node2dy = $node2d->value(1);
	


        #$fieldSet2d = $node2d->fieldSet();
        #$fieldSet2d->list();
        
        
        
        $node3dx = CmUtils::Objects::Component->new();
        $node3dx->name("x");
        $node3dx->valIndex(0);
        $node3dy = CmUtils::Objects::Component->new();
        $node3dy->name("y");
        $node3dy->valIndex(1);
        $node3dz = CmUtils::Objects::Component->new();
        $node3dz->name("z");
        $node3dz->valIndex(2);
        
        
        
        $field3d = CmUtils::Objects::Field->new();
        $field3d->name("coordinates");
        $field3d->addComponent($node3dx);
        $field3d->addComponent($node3dy);
        $field3d->addComponent($node3dz);
       
        $fieldSet3d = CmUtils::Objects::FieldSet->new();
        $fieldSet3d->addField($field3d);
        $fieldSet3d->numValues(3)
       	$fieldSet3d->list();
        $node3d->name($i);
	$node3d->fieldSet($fieldSet3d);        
        $node3d->value(0,$node2dx);
        $node3d->value(1,$node2dy);
        if($i<=$numnodes2d)
        {
        	 $node3d->value(2,0.3);
        }
        else
        {
        	 $node3d->value(2,0);
        }
       
        
        $group3d->addNode($node3d);
}
$group3d->listFields();
writeIpnode($convertedFile,$group3d);