Merge pull request #343 from ONSAS/jorge-vtk

fix bug truss vtk
ONSAS · Nov 24, 2021 · 06b7247 · 06b7247
2 parents 0806093 + a7ff8e3
commit 06b7247
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 11 deletions.
diff --git a/examples/staticVonMisesTruss/onsasExample_staticVonMisesTruss.m b/examples/staticVonMisesTruss/onsasExample_staticVonMisesTruss.m
@@ -140,7 +140,7 @@
 analysisSettings.methodName   = 'arcLength'                      ;
 analysisSettings.finalTime     = 4.5                             ;
 analysisSettings.incremArcLen = 0.1                             ;
-analysisSettings.iniDeltaLamb = boundaryConds(2).loadsTimeFact(.2)/100
+analysisSettings.iniDeltaLamb = boundaryConds(2).loadsTimeFact(.2)/100 ;
 analysisSettings.posVariableLoadBC = 2 ;
 %md
 [matUs, loadFactorsMat] = ONSAS( materials, elements, boundaryConds, initialConds, mesh, analysisSettings, otherParams ) ;
@@ -153,7 +153,7 @@
 
 verifBoolean =  ( ( norm( difLoadEngRot    ) / norm( loadFactorsNREngRot  ) ) <  1e-4 ) ...
              && ( ( norm( difLoadGreen     ) / norm( loadFactorsNRGreen   ) ) <  1e-4 ) ...
-             && ( ( norm( difLoadGreenNRAL ) / norm( loadFactorsNRALGreen ) ) <  1e-4 )
+             && ( ( norm( difLoadGreenNRAL ) / norm( loadFactorsNRALGreen ) ) <  1e-4 ) ;
 
 %md### Plots
 %md and solutions are plotted.

diff --git a/src/vtk/trussVtkData.m b/src/vtk/trussVtkData.m
@@ -25,21 +25,21 @@
     dispsElem          = U( dofsElem ) ;
     thetaLocIniSubElem = zeros(3,1)    ;
     thetaLocEndSubElem = zeros(3,1)    ;
-    dispLocIniSubElem  = zeros(3,1)    ;
-    dispLocEndSubElem  = zeros(3,1)    ;
 
     [ R0, Rr, locDisp ] = elementBeamRotData( coordsElemNodes, dispsElem ) ;
 
+    dispLocIniSubElem  = zeros(3,1)    ;
+    dispLocEndSubElem  = [ locDisp(1); 0; 0 ]    ;
+
     coordLocSubElem = [ 0; elemLength ] ;
 
     [ Nodesvtk, Conecvtk, Dispsvtk ] = vtkBeam2SolidConverter( coordsElemNodes, ...
        dispsElem, coordLocSubElem, dispLocIniSubElem, dispLocEndSubElem, thetaLocIniSubElem, thetaLocEndSubElem, sectPar, Rr, R0 ) ;
 
-
        Conecvtk( :, 2:end ) = Conecvtk(:,2:end)+counterNodes ;
-       vtkNodes             = [ vtkNodes ;     Nodesvtk ] ;
-       vtkConec             = [ vtkConec ;     Conecvtk ] ;
-       vtkNodalDisps        = [ vtkNodalDisps; Dispsvtk ] ;
+       vtkNodes             = [ vtkNodes ;     Nodesvtk ]    ;
+       vtkConec             = [ vtkConec ;     Conecvtk ]    ;
+       vtkNodalDisps        = [ vtkNodalDisps; Dispsvtk ]    ;
 
        counterNodes = counterNodes + 8 ;
 

diff --git a/src/vtk/vtkBeam2SolidConverter.m b/src/vtk/vtkBeam2SolidConverter.m
@@ -23,8 +23,7 @@
 
   typeSolid = sectPar(1) ;
 
-  dispIniSection = dispsElem(1:2:5)  ;
-  dispEndSection = dispsElem(7:2:11) ;
+  dispIniSection = dispsElem(1:2:5) ;    dispEndSection = dispsElem(7:2:11) ;
 
   if typeSolid == 12 % vtkHexa
 

diff --git a/src/vtk/vtkDataConversion.m b/src/vtk/vtkDataConversion.m
@@ -98,7 +98,7 @@
     vtkConec      = [ vtkConec ;  [currVtkConec(:,1) currVtkConec(:,2:end)+totalNodes]]           ;
     vtkNodalDisps = [ vtkNodalDisps ;  currVtkNodalDisps ] ;
 
-		totalNodes = totalNodes + size(currVtkNodes,1) 
+		totalNodes = totalNodes + size(currVtkNodes,1) ; 
 
   end % for: elemTypeInds