void TR_SHELL02 :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; FloatArray v1, v2, v3; double s [ 3 ], defScale; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( !this->giveMaterial()->isActivated(tStep) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { result += this->giveInternalStateAtNode(v1, gc.giveIntVarType(), gc.giveIntVarMode(), 1, tStep); result += this->giveInternalStateAtNode(v2, gc.giveIntVarType(), gc.giveIntVarMode(), 2, tStep); result += this->giveInternalStateAtNode(v3, gc.giveIntVarType(), gc.giveIntVarMode(), 3, tStep); } else if ( gc.giveIntVarMode() == ISM_local ) { double tot_w = 0.; FloatArray a, v; for ( GaussPoint *gp: *plate->giveDefaultIntegrationRulePtr() ) { this->giveIPValue(a, gp, IST_ShellMomentumTensor, tStep); v.add(gp->giveWeight(), a); tot_w += gp->giveWeight(); } v.times(1. / tot_w); v1 = v; v2 = v; v3 = v; } indx = gc.giveIntVarIndx(); s [ 0 ] = v1.at(indx); s [ 1 ] = v2.at(indx); s [ 2 ] = v3.at(indx); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 3; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(3, tStep, defScale); } else { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(3); } } // //EASValsSetColor(gc.getYieldPlotColor(ratio)); gc.updateFringeTableMinMax(s, 3); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }
void QTrPlaneStrain :: drawScalar(oofegGraphicContext &context) { int t, n [ 3 ], i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray v [ 6 ]; double s [ 6 ], ss [ 3 ], defScale; IntArray map; if ( !context.testElementGraphicActivity(this) ) { return; } if ( context.giveIntVarMode() == ISM_recovered ) { for ( i = 1; i <= 6; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], context.giveIntVarType(), context.giveIntVarMode(), i, tStep); } } else if ( context.giveIntVarMode() == ISM_local ) { return; } if ( result != 6 ) { return; } result = this->giveIntVarCompFullIndx( map, context.giveIntVarType() ); if ( ( !result ) || ( indx = map.at( context.giveIntVarIndx() ) ) == 0 ) { return; } for ( i = 1; i <= 6; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( context.getScalarAlgo() == SA_ISO_SURF ) { for ( t = 1; t <= 4; t++ ) { if ( t == 1 ) { n [ 0 ] = 1; n [ 1 ] = 4; n [ 2 ] = 6; } else if ( t == 2 ) { n [ 0 ] = 2; n [ 1 ] = 5; n [ 2 ] = 4; } else if ( t == 3 ) { n [ 0 ] = 3; n [ 1 ] = 6; n [ 2 ] = 5; } else { n [ 0 ] = 4; n [ 1 ] = 5; n [ 2 ] = 6; } for ( i = 0; i < 3; i++ ) { if ( context.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = context.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(n [ i ])->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ i ].y = ( FPNum ) this->giveNode(n [ i ])->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ i ].z = 0.; } else { p [ i ].x = ( FPNum ) this->giveNode(n [ i ])->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(n [ i ])->giveCoordinate(2); p [ i ].z = 0.; } } //EASValsSetColor(gc.getYieldPlotColor(ratio)); ss [ 0 ] = s [ n [ 0 ] - 1 ]; ss [ 1 ] = s [ n [ 1 ] - 1 ]; ss [ 2 ] = s [ n [ 2 ] - 1 ]; context.updateFringeTableMinMax(ss, 3); tr = CreateTriangleWD3D(p, ss [ 0 ], ss [ 1 ], ss [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } /* } else if (context.getScalarAlgo() == SA_ISO_LINE) { * * EASValsSetColor(context.getActiveCrackColor()); * EASValsSetLineWidth(OOFEG_ISO_LINE_WIDTH); * * for (t=1; t<=4; t++) { * if (t==1) {n[0] = 1; n[1]=4; n[2]=6;} * else if (t==2) {n[0]=2; n[1]=5; n[2]=4;} * else if (t==3) {n[0]=3; n[1]=6; n[2]=5;} * else {n[0]=4; n[1]=5; n[2]=6;} * * * for (i=0; i< 3; i++) { * if (context.getInternalVarsDefGeoFlag()) { * // use deformed geometry * defScale = context.getDefScale(); * p[i].x = (FPNum) this->giveNode(n[i])->giveUpdatedCoordinate(1,tStep,EID_MomentumBalance,defScale); * p[i].y = (FPNum) this->giveNode(n[i])->giveUpdatedCoordinate(2,tStep,EID_MomentumBalance,defScale); * p[i].z = 0.; * * } else { * p[i].x = (FPNum) this->giveNode(n[i])->giveCoordinate(1); * p[i].y = (FPNum) this->giveNode(n[i])->giveCoordinate(2); * p[i].z = 0.; * } * } * sv[0]=s[n[0]-1]; * sv[1]=s[n[1]-1]; * sv[2]=s[n[2]-1]; * * // isoline implementation * oofeg_drawIsoLinesOnTriangle (p, sv); * } */ } }
void Quad10_2D_SUPG :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; FloatArray v1, v2, v3; double s [ 3 ]; if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); // if ((gc.giveIntVarMode() == ISM_local) && (gc.giveIntVarType() == IST_VOFFraction)) { if ( ( gc.giveIntVarType() == IST_VOFFraction ) && ( gc.giveIntVarMode() == ISM_local ) ) { Polygon matvolpoly; //this->formMaterialVolumePoly(matvolpoly, NULL, temp_normal, temp_p, false); EASValsSetColor( gc.getStandardSparseProfileColor() ); //GraphicObj *go = matvolpoly.draw(gc,true,OOFEG_VARPLOT_PATTERN_LAYER); matvolpoly.draw(gc, true, OOFEG_VARPLOT_PATTERN_LAYER); return; } if ( gc.giveIntVarMode() == ISM_recovered ) { result += this->giveInternalStateAtNode(v1, gc.giveIntVarType(), gc.giveIntVarMode(), 1, tStep); result += this->giveInternalStateAtNode(v2, gc.giveIntVarType(), gc.giveIntVarMode(), 2, tStep); result += this->giveInternalStateAtNode(v3, gc.giveIntVarType(), gc.giveIntVarMode(), 3, tStep); } else if ( gc.giveIntVarMode() == ISM_local ) { GaussPoint *gp = integrationRulesArray [ 0 ]->getIntegrationPoint(0); result += giveIPValue(v1, gp, gc.giveIntVarType(), tStep); v2 = v1; v3 = v1; result *= 3; } if ( result != 3 ) { return; } indx = gc.giveIntVarIndx(); s [ 0 ] = v1.at(indx); s [ 1 ] = v2.at(indx); s [ 2 ] = v3.at(indx); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 3; i++ ) { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = 0.; } //EASValsSetColor(gc.getYieldPlotColor(ratio)); gc.updateFringeTableMinMax(s, 3); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } else if ( ( gc.getScalarAlgo() == SA_ZPROFILE ) || ( gc.getScalarAlgo() == SA_COLORZPROFILE ) ) { double landScale = gc.getLandScale(); for ( i = 0; i < 3; i++ ) { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = s [ i ] * landScale; } if ( gc.getScalarAlgo() == SA_ZPROFILE ) { EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetFillStyle(FILL_SOLID); tr = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | FILL_MASK | LAYER_MASK, tr); } else { gc.updateFringeTableMinMax(s, 3); EASValsSetFillStyle(FILL_SOLID); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(FILL_MASK | LAYER_MASK, tr); } EMAddGraphicsToModel(ESIModel(), tr); } }
void Axisymm3d :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; FloatArray v1, v2, v3; double s [ 3 ], defScale; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { result += this->giveInternalStateAtNode(v1, gc.giveIntVarType(), gc.giveIntVarMode(), 1, tStep); result += this->giveInternalStateAtNode(v2, gc.giveIntVarType(), gc.giveIntVarMode(), 2, tStep); result += this->giveInternalStateAtNode(v3, gc.giveIntVarType(), gc.giveIntVarMode(), 3, tStep); } else if ( gc.giveIntVarMode() == ISM_local ) { GaussPoint *gp = integrationRulesArray [ 0 ]->getIntegrationPoint(0); result += giveIPValue(v1, gp, gc.giveIntVarType(), tStep); v2 = v1; v3 = v1; result *= 3; } if ( result != 3 ) { return; } indx = gc.giveIntVarIndx(); s [ 0 ] = v1.at(indx); s [ 1 ] = v2.at(indx); s [ 2 ] = v3.at(indx); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 3; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = 0.; } else { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = 0.; } } //EASValsSetColor(gc.getYieldPlotColor(ratio)); gc.updateFringeTableMinMax(s, 3); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } else if ( ( gc.getScalarAlgo() == SA_ZPROFILE ) || ( gc.getScalarAlgo() == SA_COLORZPROFILE ) ) { double landScale = gc.getLandScale(); for ( i = 0; i < 3; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = s [ i ] * landScale; } else { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = s [ i ] * landScale; } } if ( gc.getScalarAlgo() == SA_ZPROFILE ) { EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetFillStyle(FILL_SOLID); tr = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | FILL_MASK | LAYER_MASK, tr); } else { gc.updateFringeTableMinMax(s, 3); EASValsSetFillStyle(FILL_SOLID); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(FILL_MASK | LAYER_MASK, tr); } EMAddGraphicsToModel(ESIModel(), tr); } }
void QTrPlaneStress2d :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int t, n [ 3 ], i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; FloatArray v [ 6 ]; double s [ 6 ], ss [ 3 ], defScale; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { // ========= plot recovered values ========= for ( i = 1; i <= 6; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], gc.giveIntVarType(), gc.giveIntVarMode(), i, tStep); } if ( result != 6 ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 6; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( t = 1; t <= 4; t++ ) { if ( t == 1 ) { n [ 0 ] = 1; n [ 1 ] = 4; n [ 2 ] = 6; } else if ( t == 2 ) { n [ 0 ] = 2; n [ 1 ] = 5; n [ 2 ] = 4; } else if ( t == 3 ) { n [ 0 ] = 3; n [ 1 ] = 6; n [ 2 ] = 5; } else { n [ 0 ] = 4; n [ 1 ] = 5; n [ 2 ] = 6; } for ( i = 0; i < 3; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(n [ i ])->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(n [ i ])->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = 0.; } else { // use initial geometry p [ i ].x = ( FPNum ) this->giveNode(n [ i ])->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(n [ i ])->giveCoordinate(2); p [ i ].z = 0.; } } //EASValsSetColor(gc.getYieldPlotColor(ratio)); ss [ 0 ] = s [ n [ 0 ] - 1 ]; ss [ 1 ] = s [ n [ 1 ] - 1 ]; ss [ 2 ] = s [ n [ 2 ] - 1 ]; gc.updateFringeTableMinMax(ss, 3); tr = CreateTriangleWD3D(p, ss [ 0 ], ss [ 1 ], ss [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } /* } else if (gc.getScalarAlgo() == SA_ISO_LINE) { * * EASValsSetColor(context.getActiveCrackColor()); * EASValsSetLineWidth(OOFEG_ISO_LINE_WIDTH); * * for (t=1; t<=4; t++) { * if (t==1) {n[0] = 1; n[1]=4; n[2]=6;} * else if (t==2) {n[0]=2; n[1]=5; n[2]=4;} * else if (t==3) {n[0]=3; n[1]=6; n[2]=5;} * else {n[0]=4; n[1]=5; n[2]=6;} * * * for (i=0; i< 3; i++) { * if (gc.getInternalVarsDefGeoFlag()) { * // use deformed geometry * defScale = gc.getDefScale(); * p[i].x = (FPNum) this->giveNode(n[i])->giveUpdatedCoordinate(1,tStep,defScale); * p[i].y = (FPNum) this->giveNode(n[i])->giveUpdatedCoordinate(2,tStep,defScale); * p[i].z = 0.; * * } else { * p[i].x = (FPNum) this->giveNode(n[i])->giveCoordinate(1); * p[i].y = (FPNum) this->giveNode(n[i])->giveCoordinate(2); * p[i].z = 0.; * } * } * sv[0]=s[n[0]-1]; * sv[1]=s[n[1]-1]; * sv[2]=s[n[2]-1]; * * // isoline implementation * oofeg_drawIsoLinesOnTriangle (p, sv); * } */ } } else if ( gc.giveIntVarMode() == ISM_local ) { // ========= plot local values ========= // (so far implemented for 4 Gauss points only) if ( numberOfGaussPoints != 4 ) { return; } IntArray ind(3); WCRec pp [ 6 ]; for ( i = 0; i < 6; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); pp [ i ].z = 0.; } else { // use initial geometry pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); pp [ i ].z = 0.; } } for ( GaussPoint *gp: *integrationRulesArray [ 0 ] ) { //gpCoords = gp->giveNaturalCoordinates(); switch ( gp->giveNumber() ) { case 3: ind.at(1) = 0; ind.at(2) = 3; ind.at(3) = 5; break; case 4: ind.at(1) = 1; ind.at(2) = 4; ind.at(3) = 3; break; case 2: ind.at(1) = 2; ind.at(2) = 5; ind.at(3) = 4; break; case 5: default: ind.at(1) = 3; ind.at(2) = 4; ind.at(3) = 5; } if ( giveIPValue(v [ 0 ], gp, gc.giveIntVarType(), tStep) == 0 ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 3; i++ ) { s [ i - 1 ] = v [ 0 ].at(indx); } for ( i = 0; i < 3; i++ ) { p [ i ].x = pp [ ind.at(i + 1) ].x; p [ i ].y = pp [ ind.at(i + 1) ].y; p [ i ].z = pp [ ind.at(i + 1) ].z; } gc.updateFringeTableMinMax(s, 3); EASValsSetFillStyle(FILL_SOLID); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(FILL_MASK | LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } } }
void TR_SHELL01 :: drawScalar(oofegGraphicContext &context) { int i, indx, result = 0; WCRec p [ 3 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray v1, v2, v3; double s [ 3 ], defScale; if ( !context.testElementGraphicActivity(this) ) { return; } if ( !this->giveMaterial()->isActivated(tStep) ) { return; } if ( context.giveIntVarMode() == ISM_recovered ) { result += this->giveInternalStateAtNode(v1, context.giveIntVarType(), context.giveIntVarMode(), 1, tStep); result += this->giveInternalStateAtNode(v2, context.giveIntVarType(), context.giveIntVarMode(), 2, tStep); result += this->giveInternalStateAtNode(v3, context.giveIntVarType(), context.giveIntVarMode(), 3, tStep); } else if ( context.giveIntVarMode() == ISM_local ) { int nip = plate->giveDefaultIntegrationRulePtr()->giveNumberOfIntegrationPoints(); FloatArray a, v(12); v.zero(); for (int _i=1; _i<= nip; _i++) { this->giveIPValue(a, plate->giveDefaultIntegrationRulePtr()->getIntegrationPoint(_i-1), IST_ShellForceMomentumTensor, tStep); v += a; } v.times(1./nip); v1 = v; v2 =v; v3 =v; } indx = context.giveIntVarIndx(); s [ 0 ] = v1.at(indx); s [ 1 ] = v2.at(indx); s [ 2 ] = v3.at(indx); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( context.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 3; i++ ) { if ( context.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = context.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(3, tStep, defScale); } else { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(3); } } // //EASValsSetColor(gc.getYieldPlotColor(ratio)); context.updateFringeTableMinMax(s, 3); tr = CreateTriangleWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }