Issues with a parallel model in OpenSeesMP
Posted: Tue May 17, 2022 7:53 am
I want to implement a parallel model and do a dynamic analysis for a 3d building using OpenSeesMP (for more details, please see the structure of the code at the end of this message). I got the following error message when the gravity load analysis was performed. Please, could you give me a hint about how to solve this problem?
WARNING MumpsParallelSolver::solve(void)- Error WARNING MumpsParallelSolver::solve(void)- WARNING MumpsParallelSolver::solve(void)- Error -10 returned in substitution dmumps()
cause: Matrix is Singular Numerically
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor -10 returned in substitution dmumps()
cause: Matrix is Singular Numerically
My model decomposes the building’s elements in 4 processors using the next criteria
if {$pid==0} {
Transverse beams definition
} elseif {$pid==1} {
Longitudinal beams definition
} elseif {$pid==2} {
Columns definition
} elseif {$pid==3} {
Members representing the joints
}
The gravity load was applied only in $pid==2
integrator LoadControl 1
test NormDispIncr 1.0e-12 10
algorithm Newton
numberer ParallelRCM
constraints Transformation
system Mumps
analysis Static
analyze 1
Thank you
WARNING MumpsParallelSolver::solve(void)- Error WARNING MumpsParallelSolver::solve(void)- WARNING MumpsParallelSolver::solve(void)- Error -10 returned in substitution dmumps()
cause: Matrix is Singular Numerically
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor -10 returned in substitution dmumps()
cause: Matrix is Singular Numerically
My model decomposes the building’s elements in 4 processors using the next criteria
if {$pid==0} {
Transverse beams definition
} elseif {$pid==1} {
Longitudinal beams definition
} elseif {$pid==2} {
Columns definition
} elseif {$pid==3} {
Members representing the joints
}
The gravity load was applied only in $pid==2
integrator LoadControl 1
test NormDispIncr 1.0e-12 10
algorithm Newton
numberer ParallelRCM
constraints Transformation
system Mumps
analysis Static
analyze 1
Thank you