Changeset 1025 for issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

Timestamp:

06/19/09 09:49:23 (16 years ago)

Author:

Mathieu Morlighem

Message:

Use sigma to update node position

File:

• issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp (modified) (5 diffs)

issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

@@ -99 +99 @@
         int node_partitionborder=0;
         double node_x[3];
+        double node_sigma;
         int node_onbed;
         int node_onsurface;
@@ -392 +393 @@
         ModelFetchData((void**)&model->y,NULL,NULL,model_handle,"y","Matrix","Mat");
         ModelFetchData((void**)&model->z,NULL,NULL,model_handle,"z","Matrix","Mat");
+        ModelFetchData((void**)&model->thickness,NULL,NULL,model_handle,"thickness","Matrix","Mat");
+        ModelFetchData((void**)&model->bed,NULL,NULL,model_handle,"bed","Matrix","Mat");
         ModelFetchData((void**)&model->gridonbed,NULL,NULL,model_handle,"gridonbed","Matrix","Mat");
         ModelFetchData((void**)&model->gridonsurface,NULL,NULL,model_handle,"gridonsurface","Matrix","Mat");
@@ -422 +425 @@
                 #endif
 
-
                 node_x[0]=model->x[i];
                 node_x[1]=model->y[i];
                 node_x[2]=model->z[i];
-
+                node_sigma=(model->z[i]-model->bed[i])/(model->thickness[i]);
                 
                 node_onbed=(int)model->gridonbed[i];
@@ -447 +449 @@
 
                 /*Create node using its constructor: */
-                node=new Node(node_id,node_partitionborder,node_numdofs,node_x,node_onbed,node_onsurface,node_upper_node_id,node_onshelf,node_onsheet);
+                node=new Node(node_id,node_partitionborder,node_numdofs,node_x,node_sigma,node_onbed,node_onsurface,node_upper_node_id,node_onshelf,node_onsheet);
 
                 /*set single point constraints.: */
@@ -474 +476 @@
         xfree((void**)&model->y);
         xfree((void**)&model->z);
+        xfree((void**)&model->thickness);
+        xfree((void**)&model->bed);
         xfree((void**)&model->gridonbed);
         xfree((void**)&model->gridonsurface);