idmanipulation.cpp

有限元计算程序用于计算平面板结构单元的应力和位移。

#include <iostream>
using namespace std;
#include "Plate.h"

//---------------------------------------------------------------------------

double** TPlate::ID_Matrix(){

        for (int i = 0; i < DOF; i++){
            for (int j = 0; j < Nodes; j++){
                ID[i][j] = Freedom[j][i];
            }
        }
        return ID;
}

//---------------------------------------------------------------------------

double** TPlate::Convert_ID(){

        NEQ = 0;
        for (int i = 0; i < Nodes; i++){
            for (int j = 0; j < DOF; j++){
                if (ID[j][i] > 0) goto esc1;
                   NEQ = NEQ + 1;
                   ID[j][i] = NEQ;
                   goto esc2;
esc1:           ID[j][i] = 0;
esc2:      }
        }

        return ID;
}

//---------------------------------------------------------------------------