void QTrPlaneStrain :: drawRawGeometry(oofegGraphicContext &gc) { WCRec p [ 3 ]; GraphicObj *go; if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(1); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(2); p [ 0 ].z = 0.; p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveCoordinate(1); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveCoordinate(2); p [ 1 ].z = 0.; p [ 2 ].x = ( FPNum ) this->giveNode(3)->giveCoordinate(1); p [ 2 ].y = ( FPNum ) this->giveNode(3)->giveCoordinate(2); p [ 2 ].z = 0.; go = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void Lattice2d :: drawRawCrossSections(oofegGraphicContext &gc, TimeStep *tStep) { GraphicObj *go; // if (!go) { // create new one WCRec p [ 2 ]; /* poin */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getCrossSectionColor() ); EASValsSetLayer(OOFEG_RAW_CROSSSECTION_LAYER); FloatArray coords; this->giveCrossSectionCoordinates(coords); p [ 0 ].x = ( FPNum ) coords.at(1); p [ 0 ].y = ( FPNum ) coords.at(2); p [ 0 ].z = ( FPNum ) coords.at(3); p [ 1 ].x = ( FPNum ) coords.at(4); p [ 1 ].y = ( FPNum ) coords.at(5); p [ 1 ].z = ( FPNum ) coords.at(6); go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void QTrPlaneStrain :: drawDeformedGeometry(oofegGraphicContext &gc, UnknownType type) { WCRec p [ 3 ]; GraphicObj *go; TimeStep *tStep = domain->giveEngngModel()->giveCurrentStep(); double defScale = gc.getDefScale(); if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ 0 ].z = 0.; p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ 1 ].z = 0.; p [ 2 ].x = ( FPNum ) this->giveNode(3)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ 2 ].y = ( FPNum ) this->giveNode(3)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ 2 ].z = 0.; go = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); }
void LSpace :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) { int i; WCRec p [ 8 ]; GraphicObj *go; double defScale = gc.getDefScale(); if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); EASValsSetFillStyle(FILL_SOLID); for ( i = 0; i < 8; i++ ) { 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); } go = CreateHexahedron(p); EGWithMaskChangeAttributes(WIDTH_MASK | FILL_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); }
void Brick1_ht :: drawRawGeometry(oofegGraphicContext &gc) { WCRec p [ 8 ]; GraphicObj *go; if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); EASValsSetFillStyle(FILL_SOLID); for ( int i = 0; i < 8; i++ ) { 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); } go = CreateHexahedron(p); EGWithMaskChangeAttributes(WIDTH_MASK | FILL_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void InterfaceElem2dQuad :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep) { GraphicObj *go; // if (!go) { // create new one WCRec p [ 2 ]; /* poin */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(1); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(2); p [ 0 ].z = 0.0; p [ 1 ].x = ( FPNum ) this->giveNode(3)->giveCoordinate(1); p [ 1 ].y = ( FPNum ) this->giveNode(3)->giveCoordinate(2); p [ 1 ].z = 0.0; go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); p [ 0 ].x = ( FPNum ) this->giveNode(3)->giveCoordinate(1); p [ 0 ].y = ( FPNum ) this->giveNode(3)->giveCoordinate(2); p [ 0 ].z = 0.0; p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveCoordinate(1); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveCoordinate(2); p [ 1 ].z = 0.0; go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void InterfaceElem2dQuad :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) { GraphicObj *go; // if (!go) { // create new one WCRec p [ 2 ]; /* poin */ if ( !gc.testElementGraphicActivity(this) ) { return; } double defScale = gc.getDefScale(); EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER + 1); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(1, tStep, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(2, tStep, defScale); p [ 0 ].z = 0.0; p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(1, tStep, defScale); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(2, tStep, defScale); p [ 1 ].z = 0.0; go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); p [ 0 ].x = ( FPNum ) this->giveNode(4)->giveUpdatedCoordinate(1, tStep, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(4)->giveUpdatedCoordinate(2, tStep, defScale); p [ 0 ].z = 0.0; p [ 1 ].x = ( FPNum ) this->giveNode(5)->giveUpdatedCoordinate(1, tStep, defScale); p [ 1 ].y = ( FPNum ) this->giveNode(5)->giveUpdatedCoordinate(2, tStep, defScale); p [ 1 ].z = 0.0; go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); }
void TR_SHELL02 :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep) { WCRec p [ 3 ]; GraphicObj *go; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( this->giveMaterial()->isActivated(tStep) ) { EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetFillStyle(FILL_SOLID); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(1); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(2); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveCoordinate(3); p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveCoordinate(1); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveCoordinate(2); p [ 1 ].z = ( FPNum ) this->giveNode(2)->giveCoordinate(3); p [ 2 ].x = ( FPNum ) this->giveNode(3)->giveCoordinate(1); p [ 2 ].y = ( FPNum ) this->giveNode(3)->giveCoordinate(2); p [ 2 ].z = ( FPNum ) this->giveNode(3)->giveCoordinate(3); go = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | FILL_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); } }
void LIBeam3dNL :: drawDeformedGeometry(oofegGraphicContext &gc, UnknownType type) { GraphicObj *go; if ( !gc.testElementGraphicActivity(this) ) { return; } TimeStep *tStep = domain->giveEngngModel()->giveCurrentStep(); double defScale = gc.getDefScale(); // if (!go) { // create new one WCRec p [ 2 ]; /* poin */ EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(3, tStep, EID_MomentumBalance, defScale); p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale); p [ 1 ].z = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(3, tStep, EID_MomentumBalance, defScale); go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); }
void TR_SHELL02 :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) { WCRec p [ 3 ]; GraphicObj *go; double defScale = gc.getDefScale(); if ( !gc.testElementGraphicActivity(this) ) { return; } if ( this->giveMaterial()->isActivated(tStep) ) { EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetFillStyle(FILL_SOLID); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(1, tStep, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(2, tStep, defScale); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(3, tStep, defScale); p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(1, tStep, defScale); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(2, tStep, defScale); p [ 1 ].z = ( FPNum ) this->giveNode(2)->giveUpdatedCoordinate(3, tStep, defScale); p [ 2 ].x = ( FPNum ) this->giveNode(3)->giveUpdatedCoordinate(1, tStep, defScale); p [ 2 ].y = ( FPNum ) this->giveNode(3)->giveUpdatedCoordinate(2, tStep, defScale); p [ 2 ].z = ( FPNum ) this->giveNode(3)->giveUpdatedCoordinate(3, tStep, defScale); go = CreateTriangle3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | FILL_MASK | COLOR_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); } }
void LSpace :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 8 ]; GraphicObj *tr; FloatArray v [ 8 ]; double s [ 8 ], defScale = 0.0; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { for ( i = 1; i <= 8; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], gc.giveIntVarType(), gc.giveIntVarMode(), i, tStep); } if ( result != 8 ) { return; } } else if ( gc.giveIntVarMode() == ISM_local ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 8; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 8; 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); } } gc.updateFringeTableMinMax(s, 8); tr = CreateHexahedronWD(p, s); EGWithMaskChangeAttributes(LAYER_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }
void CohesiveSurface3d :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) { GraphicObj *go1, *go2; if ( !gc.testElementGraphicActivity(this) ) { return; } double defScale = gc.getDefScale(); WCRec p [ 2 ]; EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); // get the displaced particle coordinates Particle *nodeA = ( Particle * ) giveNode(1); Particle *nodeB = ( Particle * ) giveNode(2); p [ 0 ].x = nodeA->giveUpdatedCoordinate(1, tStep, defScale); p [ 0 ].y = nodeA->giveUpdatedCoordinate(2, tStep, defScale); p [ 0 ].z = nodeA->giveUpdatedCoordinate(3, tStep, defScale); p [ 1 ].x = nodeB->giveUpdatedCoordinate(1, tStep, defScale); p [ 1 ].y = nodeB->giveUpdatedCoordinate(2, tStep, defScale); p [ 1 ].z = nodeB->giveUpdatedCoordinate(3, tStep, defScale); // plot the displaced particles EASValsSetMType(FILLED_CIRCLE_MARKER); EASValsSetColor( gc.getNodeColor() ); EASValsSetMSize(6); // plot the first particle go1 = CreateMarker3D(p); EGWithMaskChangeAttributes(COLOR_MASK | LAYER_MASK | MTYPE_MASK | MSIZE_MASK, go1); EMAddGraphicsToModel(ESIModel(), go1); // take into account periodic conditions if ( giveNumberOfNodes() == 3 ) { Node *nodeC = ( Particle * ) giveNode(3); p [ 1 ].x += kxa + kxa * defScale * ( nodeC->giveDofWithID(D_u)->giveUnknown(VM_Total, tStep) ) + kyb * defScale * ( nodeC->giveDofWithID(R_u)->giveUnknown(VM_Total, tStep) ); p [ 1 ].y += kyb + kyb * defScale * ( nodeC->giveDofWithID(D_v)->giveUnknown(VM_Total, tStep) ) + kzc * defScale * ( nodeC->giveDofWithID(R_v)->giveUnknown(VM_Total, tStep) ); p [ 1 ].z += kzc + kzc * defScale * ( nodeC->giveDofWithID(D_w)->giveUnknown(VM_Total, tStep) ) + kxa * defScale * ( nodeC->giveDofWithID(R_w)->giveUnknown(VM_Total, tStep) ); EASValsSetMType(CIRCLE_MARKER); } // plot the second particle go2 = CreateMarker3D(p + 1); EGWithMaskChangeAttributes(COLOR_MASK | LAYER_MASK | MTYPE_MASK | MSIZE_MASK, go2); EMAddGraphicsToModel(ESIModel(), go2); }
// TODO: FIX OOFEG implementation void TrPlaneStress2dXFEM :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep) { if ( !gc.testElementGraphicActivity(this) ) { return; } XfemManager *xf = this->giveDomain()->giveXfemManager(); if ( !xf->isElementEnriched(this) ) { TrPlaneStress2d :: drawRawGeometry(gc, tStep); } else { if ( integrationRulesArray.size() > 1 ) { #if 0 for ( auto &ir: integrationRulesArray ) { // TODO: Implement visualization. PatchIntegrationRule *iRule = dynamic_cast< PatchIntegrationRule * >( ir ); if ( iRule ) { iRule->givePatch()->draw(gc); } } #endif } else { TrPlaneStress2d :: drawRawGeometry(gc, tStep); } } }
// TODO: FIX OOFEG implementation void TrPlaneStress2dXFEM :: drawRawGeometry(oofegGraphicContext &context) { if ( !context.testElementGraphicActivity(this) ) { return; } XfemManager *xf = this->giveDomain()->giveXfemManager(); if ( !xf->isElementEnriched(this) ) { TrPlaneStress2d :: drawRawGeometry(context); } else { if ( numberOfIntegrationRules > 1 ) { int i; // PatchIntegrationRule *iRule; for ( i = 0; i < numberOfIntegrationRules; i++ ) { // TODO: Implement visualization. /* iRule = dynamic_cast< PatchIntegrationRule * >( integrationRulesArray [ i ] ); if ( iRule ) { iRule->givePatch()->draw(context); } */ } } else { TrPlaneStress2d :: drawRawGeometry(context); } } }
void InterfaceElem1d :: drawScalar(oofegGraphicContext &context) { int i, indx, result = 0; GaussPoint *gp; IntegrationRule *iRule = integrationRulesArray [ giveDefaultIntegrationRule() ]; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray gcoord(3), v1; WCRec p [ 1 ]; IntArray map; GraphicObj *go; double val [ 1 ]; if ( !context.testElementGraphicActivity(this) ) { return; } if ( context.getInternalVarsDefGeoFlag() ) { double defScale = context.getDefScale(); p [ 0 ].x = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale) + this->giveNode(2)->giveUpdatedCoordinate(1, tStep, EID_MomentumBalance, defScale) ); p [ 0 ].y = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale) + this->giveNode(2)->giveUpdatedCoordinate(2, tStep, EID_MomentumBalance, defScale) ); p [ 0 ].z = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(3, tStep, EID_MomentumBalance, defScale) + this->giveNode(2)->giveUpdatedCoordinate(3, tStep, EID_MomentumBalance, defScale) ); } else { p [ 0 ].x = ( FPNum )( this->giveNode(1)->giveCoordinate(1) ); p [ 0 ].y = ( FPNum )( this->giveNode(1)->giveCoordinate(2) ); p [ 0 ].z = ( FPNum )( this->giveNode(1)->giveCoordinate(3) ); } result += giveIPValue(v1, iRule->getIntegrationPoint(0), context.giveIntVarType(), tStep); for ( i = 0; i < iRule->getNumberOfIntegrationPoints(); i++ ) { result = 0; gp = iRule->getIntegrationPoint(i); result += giveIPValue(v1, gp, context.giveIntVarType(), tStep); result += this->giveIntVarCompFullIndx( map, context.giveIntVarType() ); if ( result != 2 ) { continue; } if ( ( indx = map.at( context.giveIntVarIndx() ) ) == 0 ) { return; } val [ 0 ] = v1.at(indx); context.updateFringeTableMinMax(val, 1); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); EASValsSetMType(FILLED_CIRCLE_MARKER); go = CreateMarkerWD3D(p, val [ 0 ]); EGWithMaskChangeAttributes(LAYER_MASK | FILL_MASK | MTYPE_MASK, go); EMAddGraphicsToModel(ESIModel(), go); //} } }
void IntElLine2 :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int indx, result = 0; IntegrationRule *iRule = this->giveDefaultIntegrationRulePtr(); FloatArray gcoord(3), v1; WCRec p [ 1 ]; GraphicObj *go; double val [ 1 ]; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { return; } for ( GaussPoint *gp: *iRule ) { result = 0; result += giveIPValue(v1, gp, gc.giveIntVarType(), tStep); if ( result != 1 ) { continue; } indx = gc.giveIntVarIndx(); result += this->computeGlobalCoordinates( gcoord, gp->giveNaturalCoordinates() ); p [ 0 ].x = ( FPNum ) gcoord.at(1); p [ 0 ].y = ( FPNum ) gcoord.at(2); p [ 0 ].z = 0.; val [ 0 ] = v1.at(indx); gc.updateFringeTableMinMax(val, 1); //if (val[0] > 0.) { EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); EASValsSetMType(FILLED_CIRCLE_MARKER); go = CreateMarkerWD3D(p, val [ 0 ]); EGWithMaskChangeAttributes(LAYER_MASK | FILL_MASK | MTYPE_MASK, go); EMAddGraphicsToModel(ESIModel(), go); //} } }
void TrPlaneStress2dXFEM :: drawScalar(oofegGraphicContext &context) { if ( !context.testElementGraphicActivity(this) ) { return; } XfemManager *xf = this->giveDomain()->giveXfemManager(); if ( !xf->isElementEnriched(this) ) { TrPlaneStress2d :: drawScalar(context); } else { if ( context.giveIntVarMode() == ISM_local ) { int indx; double val; FloatArray s(3), v; indx = context.giveIntVarIndx(); TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); PatchIntegrationRule *iRule; for ( int i = 0; i < numberOfIntegrationRules; i++ ) { iRule = dynamic_cast< PatchIntegrationRule * >( integrationRulesArray [ i ] ); #if 0 val = iRule->giveMaterial(); #else val = 0.0; for ( int j = 0; j < iRule->giveNumberOfIntegrationPoints(); j++ ) { GaussPoint *gp = iRule->getIntegrationPoint(0); giveIPValue(v, gp, context.giveIntVarType(), tStep); val += v.at(indx); } val /= iRule->giveNumberOfIntegrationPoints(); #endif s.at(1) = s.at(2) = s.at(3) = val; // TODO: Implement visualization. // iRule->givePatch()->drawWD(context, s); } } else { TrPlaneStress2d :: drawScalar(context); } } }
void Truss2d :: drawRawGeometry(oofegGraphicContext &gc) { int c1, c2; resolveCoordIndices(c1, c2); GraphicObj *go; // if (!go) { // create new one WCRec p [ 2 ]; /* point */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); if ( cs_mode == 0 ) { p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(c1); p [ 0 ].y = 0.; p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveCoordinate(c2); p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveCoordinate(c1); p [ 1 ].y = 0.; p [ 1 ].z = ( FPNum ) this->giveNode(2)->giveCoordinate(c2); } else if ( cs_mode == 1 ) { p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(c1); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(c2); p [ 0 ].z = 0.; p [ 1 ].x = ( FPNum ) this->giveNode(2)->giveCoordinate(c1); p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveCoordinate(c2); p [ 1 ].z = 0.; } else if ( cs_mode == 2 ) { p [ 0 ].x = 0.; p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(c1); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveCoordinate(c2); p [ 1 ].x = 0.; p [ 1 ].y = ( FPNum ) this->giveNode(2)->giveCoordinate(c1); p [ 1 ].z = ( FPNum ) this->giveNode(2)->giveCoordinate(c2); } go = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void CohesiveSurface3d :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep) { if ( !gc.testElementGraphicActivity(this) ) { return; } //WCRec p[4]; GraphicObj *go; Particle *nodeA = ( Particle * ) this->giveNode(1); Particle *nodeB = ( Particle * ) this->giveNode(2); //double rA = nodeA -> giveRadius(); //double rB = nodeB -> giveRadius(); //double r = (rA+rB)/4.; EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH); EASValsSetColor( gc.getElementColor() ); EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); WCRec pl [ 2 ]; // determine coordinates of the particles connected by this element pl [ 0 ].x = ( FPNum ) nodeA->giveCoordinate(1); pl [ 0 ].y = ( FPNum ) nodeA->giveCoordinate(2); pl [ 0 ].z = ( FPNum ) nodeA->giveCoordinate(3); pl [ 1 ].x = ( FPNum ) nodeB->giveCoordinate(1); pl [ 1 ].y = ( FPNum ) nodeB->giveCoordinate(2); pl [ 1 ].z = ( FPNum ) nodeB->giveCoordinate(3); if ( giveNumberOfNodes() == 3 ) { // the second particle should be shifted (periodic arrangement) Particle *nodeC = ( Particle * ) this->giveNode(3); pl [ 1 ].x += kx * ( nodeC->giveCoordinate(1) ); pl [ 1 ].y += ky * ( nodeC->giveCoordinate(2) ); pl [ 1 ].z += kz * ( nodeC->giveCoordinate(3) ); } // plot a line segment connecting the particles go = CreateLine3D(pl); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void LumpedMassElement :: drawRawGeometry(oofegGraphicContext &gc) { GraphicObj *go; WCRec p [ 1 ]; /* point */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetColor( gc.getElementColor() ); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); EASValsSetMType(SQUARE_MARKER); EASValsSetMSize(8); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(1); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveCoordinate(2); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveCoordinate(3); go = CreateMarker3D(p); EGWithMaskChangeAttributes(COLOR_MASK | LAYER_MASK | MTYPE_MASK | MSIZE_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void Lattice2d_mt :: drawYourself(oofegGraphicContext &gc) { OGC_PlotModeType mode = gc.giveIntVarPlotMode(); if ( mode == OGC_rawGeometry ) { this->drawRawGeometry(gc); this->drawRawCrossSections(gc); } else { _error("drawYourself : unsupported mode"); } }
void Lattice2d_mt :: drawYourself(oofegGraphicContext &gc, TimeStep *tStep) { OGC_PlotModeType mode = gc.giveIntVarPlotMode(); if ( mode == OGC_rawGeometry ) { this->drawRawGeometry(gc, tStep); this->drawRawCrossSections(gc, tStep); } else { OOFEM_ERROR("unsupported mode"); } }
void LumpedMassElement :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type) { GraphicObj *go; double defScale = gc.getDefScale(); WCRec p [ 1 ]; /* point */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetColor( gc.getDeformedElementColor() ); EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER); EASValsSetMType(SQUARE_MARKER); EASValsSetMSize(8); p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(1, tStep, defScale); p [ 0 ].y = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(2, tStep, defScale); p [ 0 ].z = ( FPNum ) this->giveNode(1)->giveUpdatedCoordinate(3, tStep, defScale); go = CreateMarker3D(p); EGWithMaskChangeAttributes(COLOR_MASK | LAYER_MASK | MTYPE_MASK | MSIZE_MASK, go); EGAttachObject(go, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), go); }
void TrPlaneStress2dXFEM :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { if ( !gc.testElementGraphicActivity(this) ) { return; } XfemManager *xf = this->giveDomain()->giveXfemManager(); if ( !xf->isElementEnriched(this) ) { TrPlaneStress2d :: drawScalar(gc, tStep); } else { if ( gc.giveIntVarMode() == ISM_local ) { int indx; double val; FloatArray s(3), v; indx = gc.giveIntVarIndx(); for ( auto &ir: integrationRulesArray ) { PatchIntegrationRule *iRule = dynamic_cast< PatchIntegrationRule * >(ir); #if 0 val = iRule->giveMaterial(); #else val = 0.0; for ( GaussPoint *gp: *iRule ) { giveIPValue(v, gp, gc.giveIntVarType(), tStep); val += v.at(indx); } val /= iRule->giveNumberOfIntegrationPoints(); #endif s.at(1) = s.at(2) = s.at(3) = val; // TODO: Implement visualization. // iRule->givePatch()->drawWD(gc, s); } } else { TrPlaneStress2d :: drawScalar(gc, tStep); } } }
void IntElPoint :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep) { GraphicObj *go; // if (!go) { // create new one WCRec p [ 1 ]; /* poin */ if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetColor( gc.getElementColor() ); EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER); EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); EASValsSetColor( gc.getDeformedElementColor() ); p [ 0 ].x = ( FPNum ) ( this->giveNode(1)->giveCoordinate(1) ); p [ 0 ].y = ( FPNum ) ( this->giveNode(1)->giveCoordinate(2) ); p [ 0 ].z = ( FPNum ) ( this->giveNode(1)->giveCoordinate(3) ); EASValsSetMType(CIRCLE_MARKER); go = CreateMarker3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, go); EMAddGraphicsToModel(ESIModel(), go); }
void Brick1_ht :: drawScalar(oofegGraphicContext &context) { int indx, result = 0; WCRec p [ 8 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray v [ 8 ]; double s [ 8 ]; if ( !context.testElementGraphicActivity(this) ) { return; } if ( context.giveIntVarMode() == ISM_recovered ) { for ( int i = 1; i <= 8; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], context.giveIntVarType(), context.giveIntVarMode(), i, tStep); } if ( result != 8 ) { return; } } else if ( context.giveIntVarMode() == ISM_local ) { return; } indx = context.giveIntVarIndx(); for ( int i = 1; i <= 8; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } EASValsSetEdgeColor( context.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( context.getScalarAlgo() == SA_ISO_SURF ) { for ( int i = 0; i < 8; i++ ) { 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); } context.updateFringeTableMinMax(s, 8); tr = CreateHexahedronWD(p, s); EGWithMaskChangeAttributes(LAYER_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }
void Tetrah1_ht :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 4 ]; GraphicObj *tr; FloatArray v [ 4 ]; double s [ 4 ]; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarMode() == ISM_recovered ) { for ( i = 1; i <= 4; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], gc.giveIntVarType(), gc.giveIntVarMode(), i, tStep); } if ( result != 4 ) { return; } } else if ( gc.giveIntVarMode() == ISM_local ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } EASValsSetEdgeColor( gc.getElementEdgeColor() ); EASValsSetEdgeFlag(true); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 4; i++ ) { 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); } gc.updateFringeTableMinMax(s, 4); tr = CreateTetraWD(p, s); EGWithMaskChangeAttributes(LAYER_MASK | EDGE_COLOR_MASK | EDGE_FLAG_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }
void IntElPoint :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int indx, result = 0; IntegrationRule *iRule = this->giveDefaultIntegrationRulePtr(); FloatArray gcoord(3), v1; WCRec p [ 1 ]; GraphicObj *go; double val [ 1 ]; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.getInternalVarsDefGeoFlag() ) { double defScale = gc.getDefScale(); p [ 0 ].x = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(1, tStep, defScale) + this->giveNode(2)->giveUpdatedCoordinate(1, tStep, defScale) ); p [ 0 ].y = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(2, tStep, defScale) + this->giveNode(2)->giveUpdatedCoordinate(2, tStep, defScale) ); p [ 0 ].z = ( FPNum ) 0.5 * ( this->giveNode(1)->giveUpdatedCoordinate(3, tStep, defScale) + this->giveNode(2)->giveUpdatedCoordinate(3, tStep, defScale) ); } else { p [ 0 ].x = ( FPNum ) ( this->giveNode(1)->giveCoordinate(1) ); p [ 0 ].y = ( FPNum ) ( this->giveNode(1)->giveCoordinate(2) ); p [ 0 ].z = ( FPNum ) ( this->giveNode(1)->giveCoordinate(3) ); } result += giveIPValue(v1, iRule->getIntegrationPoint(0), gc.giveIntVarType(), tStep); for ( GaussPoint *gp: *iRule ) { result = 0; result += giveIPValue(v1, gp, gc.giveIntVarType(), tStep); if ( result != 1 ) { continue; } indx = gc.giveIntVarIndx(); val [ 0 ] = v1.at(indx); gc.updateFringeTableMinMax(val, 1); EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); EASValsSetMType(FILLED_CIRCLE_MARKER); go = CreateMarkerWD3D(p, val [ 0 ]); EGWithMaskChangeAttributes(LAYER_MASK | FILL_MASK | MTYPE_MASK, go); EMAddGraphicsToModel(ESIModel(), go); //} } }
void Lattice2d :: drawYourself(oofegGraphicContext &gc, TimeStep *tStep) { OGC_PlotModeType mode = gc.giveIntVarPlotMode(); if ( mode == OGC_rawGeometry ) { this->drawRawGeometry(gc, tStep); this->drawRawCrossSections(gc, tStep); } else if ( mode == OGC_deformedGeometry ) { this->drawDeformedGeometry(gc, tStep, DisplacementVector); } else if ( mode == OGC_eigenVectorGeometry ) { this->drawDeformedGeometry(gc, tStep, EigenVector); } else if ( mode == OGC_scalarPlot ) { this->drawScalar(gc, tStep); } else if ( mode == OGC_elemSpecial ) { this->drawSpecial(gc, tStep); } else { OOFEM_ERROR("unsupported mode"); } }
void Lattice2d :: drawSpecial(oofegGraphicContext &gc, TimeStep *tStep) { WCRec p [ 2 ]; GraphicObj *tr; GaussPoint *gp; FloatArray crackStatuses, cf; if ( !gc.testElementGraphicActivity(this) ) { return; } if ( gc.giveIntVarType() == IST_CrackState ) { gp = integrationRulesArray [ 0 ]->getIntegrationPoint(0); this->giveIPValue(crackStatuses, gp, IST_CrackStatuses, tStep); if ( crackStatuses(0) == 1. || crackStatuses(0) == 2. || crackStatuses(0) == 3 || crackStatuses(0) == 4 ) { FloatArray coords; this->giveCrossSectionCoordinates(coords); p [ 0 ].x = ( FPNum ) coords.at(1); p [ 0 ].y = ( FPNum ) coords.at(2); p [ 0 ].z = ( FPNum ) coords.at(3); p [ 1 ].x = ( FPNum ) coords.at(4); p [ 1 ].y = ( FPNum ) coords.at(5); p [ 1 ].z = ( FPNum ) coords.at(6); EASValsSetLayer(OOFEG_CRACK_PATTERN_LAYER); EASValsSetLineWidth(OOFEG_CRACK_PATTERN_WIDTH); if ( ( crackStatuses(0) == 1. ) ) { EASValsSetColor( gc.getActiveCrackColor() ); } else if ( crackStatuses(0) == 2. ) { EASValsSetColor( gc.getCrackPatternColor() ); } else if ( crackStatuses(0) == 3. ) { EASValsSetColor( gc.getActiveCrackColor() ); } else if ( crackStatuses(0) == 4. ) { EASValsSetColor( gc.getActiveCrackColor() ); } tr = CreateLine3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, tr); EGAttachObject(tr, ( EObjectP ) this); EMAddGraphicsToModel(ESIModel(), tr); } } }