Search found 5 matches

Go to advanced search

by cmckee
Wed Apr 28, 2021 12:13 pm
Forum: Framework
Topic: save/restore commands with OpenSeesSP
Replies: 1
Views: 9000

save/restore commands with OpenSeesSP

When using OpenSeesSP, I am unable to use the save/restore commands. I am getting:

Code: Select all

ShadowSubdomain::sendSelf()  - NOT YET IMPLEMENTED
Domain::send - element with tag 2 failed in sendSelf
FE_Datastore::commitState - domain failed to sendSelf
WARNING - database failed to commitState
When running the exact input script with sequential OpenSees, there are no issues — as expected. I have searched through as many resources as I could find, and I have not found any documentation indicating that the database commands are incompatible with OpenSeesSP.

If these functions are not possible with the parallel interpreters, what is missing and/or still needs to be implemented and is there a particular reason why it hasn't been implemented yet?
by cmckee
Mon Apr 12, 2021 3:05 pm
Forum: Parallel Processing
Topic: MPI_Channel::recvVector() - incorrect number of entries for Vector received
Replies: 4
Views: 11450

Re: MPI_Channel::recvVector() - incorrect number of entries for Vector received

I believe it is related to the NormUnbalance convergence test. If I change to NormDispIncr, the issue goes away.
by cmckee
Mon Apr 12, 2021 1:57 pm
Forum: Parallel Processing
Topic: MPI_Channel::recvVector() - incorrect number of entries for Vector received
Replies: 4
Views: 11450

Re: MPI_Channel::recvVector() - incorrect number of entries for Vector received

Code: Select all

vectX.size() - 24
vectX[0] = -4.43559e-06 -0.000700246 1.16597e-06 -2.14741e-06 -0.00130121 9.16138e-07 2.14741e-06 -0.00130121 -9.16138e-07 4.43559e-06 -0.000700246 -1.16597e-06 3.62509e-06 -0.00051289 3.78252e-07 1.75734e-06 -0.000940812 3.64984e-07 -1.75734e-06 -0.000940812 -3.64984e-07 -3.62509e-06 -0.00051289 -3.78252e-07
vectX[1] =
vectX[2] =
vectX[3] =
vectX[4] = 0 0 0
vectX[5] =
vectX[6] =
vectX[7] =
vectX[8] = 0.000155822 -3.1354 1.16593e-06 0.000147964 -6.34034 9.16101e-07 -0.000147964 -6.34034 -9.16101e-07 -0.000155822 -3.1354 -1.16593e-06 0.15083 -3.81564 2.17735 -0.1511 -6.99583 1.58268 0.000127963 -7.01149 -3.6497e-07 0.000141536 -3.79998 -3.78237e-07
vectX[9] =
vectX[10] =
vectX[11] =
vectX[12] = -0.00145004 -0.205156 -0.000151301
vectX[13] =
vectX[14] =
vectX[15] =
vectX[16] = 0 0 0
vectX[17] =
vectX[18] =
vectX[19] =
vectX[20] = 4.43559e-06 -0.000700246 -1.16597e-06
vectX[21] =
vectX[22] =
vectX[23] =
by cmckee
Mon Apr 12, 2021 12:59 pm
Forum: Parallel Processing
Topic: MPI_Channel::recvVector() - incorrect number of entries for Vector received
Replies: 4
Views: 11450

MPI_Channel::recvVector() - incorrect number of entries for Vector received

Howdy,

When using the Mumps system, I am receiving the following error:

Code: Select all

MPI_Channel::recvVector() - incorrect number of entries for Vector received: -32766 expected: 24
The call stack leading up to this failure is:

Code: Select all

	OpenSeesSP.exe!MPI_Channel::recvVector(int dbTag, int commitTag, Vector & theVector, ChannelAddress * theAddress) Line 340	C++
 	OpenSeesSP.exe!MumpsParallelSOE::solve() Line 353	C++
 	OpenSeesSP.exe!NewtonRaphson::solveCurrentStep() Line 187	C++
 	OpenSeesSP.exe!DirectIntegrationAnalysis::analyzeStep(double dT) Line 215	C++
 	OpenSeesSP.exe!DirectIntegrationAnalysis::analyze(int numSteps, double dT) Line 173	C++
 	OpenSeesSP.exe!analyzeModel(void * clientData, Tcl_Interp * interp, int argc, const char * * argv) Line 1902	C++
 	[External Code]	
 	OpenSeesSP.exe!OPS_SourceCmd(void * dummy, Tcl_Interp * interp, int objc, Tcl_Obj * const * objv) Line 1179	C++
 	[External Code]	
 	OpenSeesSP.exe!OPS_SourceCmd(void * dummy, Tcl_Interp * interp, int objc, Tcl_Obj * const * objv) Line 1179	C++
 	[External Code]	
 	OpenSeesSP.exe!g3TclMain(int argc, char * * argv, int(*)(Tcl_Interp *) appInitProc, int rank, int np) Line 443	C++
 	OpenSeesSP.exe!main(int argc, char * * argv) Line 275	C++
Error occurs when the MumpsParallelSOE is receiving the *vectX from each channel. A quick try with SuperLU shows the same issue.

Error occurs using both IMPI and OpenMpi.

Is there anything obvious that may be the issue here? I can provide more details if necessary.

Thank you
by cmckee
Sun Apr 04, 2021 9:28 am
Forum: Framework
Topic: Accessing Section and Fiber data from material
Replies: 1
Views: 9129

Accessing Section and Fiber data from material

Is it possible to access the information related to section type or fiber area from within a derived UniaxialMaterial?

For example, I have a macro-mechanical material model that requires the cross-sectional area of the bar - I can input that parameter no problem. Within this model, I also include a time-dependent cross-section and strength reduction due to corrosion.

When used in a straight layer fiber, F = stress *A. Without redefining the area of the fiber, to get the appropriate resulting force in the fiber I need to additionally scale the stress from the material model relative to the fiber area. But this is not necessary if the material is used with a uniaxial section. I would also like to be able to include a check to make sure the area input into the material model matches the area assigned to the fiber.

Go to advanced search