I am writing an FEM code for a curved box girder bridge and getting an error..Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 9.359430e-32. > In Dec14 (line 579)
1 次查看(过去 30 天)
显示 更早的评论
X_matrix=inv(K_modi)*F_matrix;
1 个评论
madhan ravi
2018-12-15
编辑:madhan ravi
2018-12-15
upload datas of K_modi and F_matrix if there are few values if more attach it as a .mat file
回答(3 个)
Walter Roberson
2018-12-16
with rcond that small you can be sure that your matrix is singular and that no inverse exists . If you check with rank() it is likely to show you that it does not have full rank.
You need to recheck your mesh construction. It could be the case that your solution technique is not applicable to your situation .
You could silence the error by using pinv instead if inv but the answer is not likely to be meaningful for the situation .
0 个评论
John D'Errico
2018-12-16
The trick with FEM is to recognize there are several potential reasons for a singular matrix. The obvious one is you might have constructed it improperly.
Perhaps as likely a reason is you have not properly constrained the problem. This is an easy thing to mistake. For example, consider a very simple truss, but one where you have not anchored at least one node of the truss to a specfific location in space. The result will be you can translate the entire truss anywhere, by adding a simple offset to the entire truss, without introducing any energy penalty at all in the form of deformation of that truss. In that case, the system of equations will be singular, a singular matrix. Is that a problem? In a sense, yes, because the solution to the system is not unique.
Actually, pinv will solve the problem in that case. At least, it will effectively resolve the issue. Better however is to resolve the problem with the correct set of constraints to prevent the truss from penalty free translations, or perhaps from energy free rotations.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Structural Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!