Method initialize

SRC/element/RockingBC/RockingBC.cpp:452–506

}

int
RockingBC::initialize(Node *nodeIPointer, Node *nodeJPointer)
{
	int error;

	nodeIPtr = nodeIPointer;
	nodeJPtr = nodeJPointer;

	if ((!nodeIPtr) || (!nodeJPtr))
	{
		opserr << "\nRockingBC::initialize";
		opserr << "\ninvalid pointers to the element nodes\n";
		return -1;
	}

	// get element length and orientation
	if ((error = this->compElemtLengthAndOrient()))
		return error;

	fr_calc();

	k1 = 1.0 / fr(2, 2);
	k2 = fr(2, 1) / fr(2, 2);
	frr(0, 0) = fr(0, 0);
	frr(1, 1) = fr(1, 1) - fr(1, 2) * fr(2, 1) / fr(2, 2);

	dw1_due = Vector(6); dw1_due(2) = 1.0;
	dr_due = Vector(6); dr_due(4) = 1.0 / L; dr_due(1) = -1.0 / L;
	dw2_due = Vector(6); dw2_due(5) = 1.0;

	durth_due = Matrix(2, 6);
	durth_dW = Matrix(2, Nw);

	CC = Matrix(Nw, 2);
	for (int i = 0; i != Nw; i++) {
		CC(i, 0) = 1.0;
		CC(i, 1) = Yw(i);
	}

	BB(0, 0) = 1. / b;
	BB(1, 1) = 1.;

	CB = CC*BB;

	dutar_due = Matrix(2, 6);
	dutar_dW = Matrix(2, Nw);

	TF1(0, 0) = 1;
	TF1(1, 1) = b;
	TF1(2, 1) = -b;
	TF1(2, 2) = -L;
	TF = b*w*TF1*E;

	return 0;
}

int
RockingBC::compElemtLengthAndOrient(void)