Changeset 4002 for issm/trunk/src/m/solutions/jpl/CreateFemModel.m

Timestamp:

06/02/10 17:42:51 (15 years ago)

Author:

Eric.Larour

Message:

Added new Verticesdof and NodesDof modules. Deleted Dof module. Created new counters for loads and constraints. Presorting grouped together. Update of counters. CreateDataSets simplified.

File:

• issm/trunk/src/m/solutions/jpl/CreateFemModel.m (modified) (1 diff)

issm/trunk/src/m/solutions/jpl/CreateFemModel.m

@@ -6 +6 @@
 %-----------------------------------------------------------------------
 
-function  m=CreateFEMModel(md)
+function  m=CreateFEMModel(md,analyses)
         
         displaystring(md.verbose,'\n   reading data from model %s...',md.name);
-        [m.elements,m.nodes,m.vertices,m.constraints,m.loads,m.materials,m.parameters]=ModelProcessor(md);
+        [m.elements,m.nodes,m.vertices,m.constraints,m.loads,m.materials,m.parameters]=ModelProcessor(md,analyses);
 
-        displaystring(md.verbose,'%s','   generating degrees of freedom...');
-        [m.nodes,m.vertices,m.part,m.tpart]=Dof(m.elements,m.nodes,m.vertices, m.parameters);
+        for i=1:length(analyses),
 
-        displaystring(md.verbose,'%s','   generating single point constraints...');
-        [m.nodes,m.yg]=SpcNodes(m.nodes,m.constraints);
+                displaystring(md.verbose,'%s','   generating degrees of freedom...');
+                if ~isfield(m,'part') [m.vertices,m.part,m.tpart]=VerticesDof(vertices, m.parameters);
+                [m.nodes]=NodesDof(m.nodes,m.parameters);
 
-        displaystring(md.verbose,'%s','   generating rigid body constraints...');
-        [m.Rmg,m.nodes]=MpcNodes(m.nodes,m.constraints);
+                displaystring(md.verbose,'%s','   generating single point constraints...');
+                [m.nodes,m.yg]=SpcNodes(m.nodes,m.constraints);
 
-        displaystring(md.verbose,'%s','   generating node sets...');
-        m.nodesets=BuildNodeSets(m.nodes);
+                displaystring(md.verbose,'%s','   generating rigid body constraints...');
+                [m.Rmg,m.nodes]=MpcNodes(m.nodes,m.constraints);
 
-        displaystring(md.verbose,'%s','   reducing single point constraints vector...');
-        m.ys=Reducevectorgtos(m.yg.vector,m.nodesets);
-        
-        displaystring(md.verbose,'%s','   normalizing rigid body constraints matrix...');
-        m.Gmn = NormalizeConstraints(m.Rmg,m.nodesets);
+                displaystring(md.verbose,'%s','   generating node sets...');
+                m.nodesets=BuildNodeSets(m.nodes);
 
-        displaystring(md.verbose,'%s','   configuring element and loads...');
-        [m.elements,m.loads,m.nodes,m.parameters] = ConfigureObjects( m.elements, m.loads, m.nodes, m.vertices,m.materials,m.parameters);
+                displaystring(md.verbose,'%s','   reducing single point constraints vector...');
+                m.ys=Reducevectorgtos(m.yg.vector,m.nodesets);
+                
+                displaystring(md.verbose,'%s','   normalizing rigid body constraints matrix...');
+                m.Gmn = NormalizeConstraints(m.Rmg,m.nodesets);
 
-        displaystring(md.verbose,'%s','   processing parameters...');
-        m.parameters= ProcessParams(m.parameters,m.part.vector);
+                displaystring(md.verbose,'%s','   configuring element and loads...');
+                [m.elements,m.loads,m.nodes,m.parameters] = ConfigureObjects( m.elements, m.loads, m.nodes, m.vertices,m.materials,m.parameters);
+
+                displaystring(md.verbose,'%s','   processing parameters...');
+                m.parameters= ProcessParams(m.parameters,m.part.vector);
+        end
 
 end