回复 #4 zhangmeng 的帖子
现在我想做一个碎石桩复合地基的数值模拟分析，但是存在下面的困难，首先，在装配时，Assembly具体的意思是什么？不同的部件在装配时，如果各个部件的网格分割的不相等，在没有设置接触单元时，这几个部件能认为它们已经连在一起了吗？如果不行的话，请问大哥：我应该怎么做？求教：
我现在做一个碎石桩复合地基的数值模拟，现在出了点问题，下面是一些警告：***WARNING: THERE ARE 5 UNCONNECTED REGIONS IN THE MODEL.
EQUATION ARE BEING REORDERED TO MINIMIZE WAVEFRONT
COLLECTING MODEL CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS
COLLECTING STEP CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS
INCREMENT 1 STARTS. ATTEMPT NUMBER1, TIME INCREMENT 1.00
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
IN AN EIGENVALUE EXTRACTION STEP THE NUMBER OF NEGATIVE EIGENVALUES IS
EQUAL TO THE NUMBER OF EIGENVALUES BELOW THE CURRENT SHIFT POINT.
THIS MAY BE USED TO CHECK THAT EIGENVALUES HAVE NOT BEEN MISSED.
NOTE: THE LANCZOS EIGENSOLVER APPLIES AN INTERNAL SHIFT WHICH WILL
RESULT IN NEGATIVE EIGENVALUES.
IN A DIRECT-SOLUTION STEADY-STATE DYNAMIC ANALYSIS, NEGATIVE
EIGENVALUES ARE EXPECTED. A STATIC ANALYSIS CAN BE USED TO VERIFY THAT THE SYSTEM IS STABLE.
IN OTHER CASES, NEGATIVE EIGENVALUES MEAN THAT THE SYSTEM MATRIX IS NOT
POSITIVE DEFINITE:
FOR EXAMPLE, A BIFURCATION (BUCKLING) LOAD MAY HAVE BEEN EXCEEDED.
NEGATIVE EIGENVALUES MAY ALSO OCCUR IF QUADRATIC ELEMENTS ARE
USED TO DEFINE CONTACT SURFACES.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_1_1.
EQUILIBRIUM ITERATION 1
AVERAGE FORCE 0.441 TIME AVG. FORCE 0.441
LARGEST RESIDUAL FORCE -11.1 AT NODE 15 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 1.107E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.107E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_2_1.
EQUILIBRIUM ITERATION 2
AVERAGE FORCE 0.513 TIME AVG. FORCE 0.513
LARGEST RESIDUAL FORCE -14.1 AT NODE 31 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 2.215E+05 AT NODE 2 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.107E+05 AT NODE 33 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_3_1.
EQUILIBRIUM ITERATION 3
AVERAGE FORCE 0.529 TIME AVG. FORCE 0.529
LARGEST RESIDUAL FORCE -11.2 AT NODE 8 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 3.322E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.107E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_4_1.
EQUILIBRIUM ITERATION 4
AVERAGE FORCE 0.629 TIME AVG. FORCE 0.629
LARGEST RESIDUAL FORCE -13.8 AT NODE 32 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 4.430E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.107E+05 AT NODE 3 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_5_1.
EQUILIBRIUM ITERATION 5
AVERAGE FORCE 0.665 TIME AVG. FORCE 0.665
LARGEST RESIDUAL FORCE 22.9 AT NODE 22 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 5.537E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.107E+05 AT NODE 17 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: THE SOLUTION APPEARS TO BE DIVERGING.
***NOTE: THE INCREMENT WILL BE ATTEMPTED AGAIN WITH A TIME INCREMENT OF 0.25000
INCREMENT 1 STARTS. ATTEMPT NUMBER2, TIME INCREMENT0.250
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_1_2.
EQUILIBRIUM ITERATION 1
AVERAGE FORCE 0.110 TIME AVG. FORCE 0.110
LARGEST RESIDUAL FORCE -2.76 AT NODE 15 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 2.769E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 2.769E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_2_2.
EQUILIBRIUM ITERATION 2
AVERAGE FORCE 0.128 TIME AVG. FORCE 0.128
LARGEST RESIDUAL FORCE -3.52 AT NODE 31 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 5.537E+04 AT NODE 2 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 2.769E+04 AT NODE 33 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_3_2.
EQUILIBRIUM ITERATION 3
AVERAGE FORCE 0.132 TIME AVG. FORCE 0.132
LARGEST RESIDUAL FORCE -2.81 AT NODE 8 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 8.306E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 2.769E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_4_2.
EQUILIBRIUM ITERATION 4
AVERAGE FORCE 0.157 TIME AVG. FORCE 0.157
LARGEST RESIDUAL FORCE -3.44 AT NODE 32 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 1.107E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 2.769E+04 AT NODE 3 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_5_2.
EQUILIBRIUM ITERATION 5
AVERAGE FORCE 0.166 TIME AVG. FORCE 0.166
LARGEST RESIDUAL FORCE 5.73 AT NODE 22 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 1.384E+05 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 2.769E+04 AT NODE 17 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: THE SOLUTION APPEARS TO BE DIVERGING.
***NOTE: THE INCREMENT WILL BE ATTEMPTED AGAIN WITH A TIME INCREMENT OF
6.25000E-02
INCREMENT 1 STARTS. ATTEMPT NUMBER3, TIME INCREMENT6.250E-02
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_1_3.
EQUILIBRIUM ITERATION 1
AVERAGE FORCE 2.755E-02 TIME AVG. FORCE 2.755E-02
LARGEST RESIDUAL FORCE -0.691 AT NODE 15 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 6.922E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 6.922E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_2_3.
EQUILIBRIUM ITERATION 2
AVERAGE FORCE 3.204E-02 TIME AVG. FORCE 3.204E-02
LARGEST RESIDUAL FORCE -0.879 AT NODE 31 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 1.384E+04 AT NODE 2 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 6.922E+03 AT NODE 33 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
这是什么意思？请各位大哥帮帮忙？到底问题是什么造成的？本人急需要解决的问题，谢谢！
求教：
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_3_3.
EQUILIBRIUM ITERATION 3
AVERAGE FORCE 3.307E-02 TIME AVG. FORCE 3.307E-02
LARGEST RESIDUAL FORCE -0.702 AT NODE 8 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 2.076E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 6.922E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_4_3.
EQUILIBRIUM ITERATION 4
AVERAGE FORCE 3.932E-02 TIME AVG. FORCE 3.932E-02
LARGEST RESIDUAL FORCE -0.861 AT NODE 32 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 2.769E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 6.922E+03 AT NODE 3 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_5_3.
EQUILIBRIUM ITERATION 5
AVERAGE FORCE 4.157E-02 TIME AVG. FORCE 4.157E-02
LARGEST RESIDUAL FORCE 1.43 AT NODE 22 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 3.461E+04 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 6.922E+03 AT NODE 17 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: THE SOLUTION APPEARS TO BE DIVERGING.
***NOTE: THE INCREMENT WILL BE ATTEMPTED AGAIN WITH A TIME INCREMENT OF
1.56250E-02
INCREMENT 1 STARTS. ATTEMPT NUMBER4, TIME INCREMENT1.563E-02
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_1_4.
EQUILIBRIUM ITERATION 1
AVERAGE FORCE 6.887E-03 TIME AVG. FORCE 6.887E-03
LARGEST RESIDUAL FORCE -0.173 AT NODE 15 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 1.730E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.730E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_2_4.
EQUILIBRIUM ITERATION 2
AVERAGE FORCE 8.009E-03 TIME AVG. FORCE 8.009E-03
LARGEST RESIDUAL FORCE -0.220 AT NODE 31 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 3.461E+03 AT NODE 2 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.730E+03 AT NODE 33 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_3_4.
EQUILIBRIUM ITERATION 3
AVERAGE FORCE 8.266E-03 TIME AVG. FORCE 8.266E-03
LARGEST RESIDUAL FORCE -0.175 AT NODE 8 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 5.191E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.730E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_4_4.
EQUILIBRIUM ITERATION 4
AVERAGE FORCE 9.829E-03 TIME AVG. FORCE 9.829E-03
LARGEST RESIDUAL FORCE -0.215 AT NODE 32 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 6.922E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.730E+03 AT NODE 3 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 22
INSTANCE BASE3-1 D.O.F. 2 RATIO = 1.65422E+15
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 26
INSTANCE CZB1-1 D.O.F. 2 RATIO = 8.00767E+14
***WARNING: SOLVER PROBLEM.NUMERICAL SINGULARITY WHEN PROCESSING NODE 13
INSTANCE DIANCENG1-1 D.O.F. 2 RATIO = 1.71125E+15
***WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: Solver problem. Numerical singularity when processing D.O.F. 2 of
3 nodes. The nodes have been identified in node set
WarnNodeSolvProbNumSing_2_1_1_5_4.
EQUILIBRIUM ITERATION 5
AVERAGE FORCE 1.039E-02 TIME AVG. FORCE 1.039E-02
LARGEST RESIDUAL FORCE 0.358 AT NODE 22 DOF2
INSTANCE: CZB1-1
LARGEST INCREMENT OF DISP. 8.652E+03 AT NODE 4 DOF2
INSTANCE: CZB1-1
LARGEST CORRECTION TO DISP. 1.730E+03 AT NODE 17 DOF2
INSTANCE: CZB1-1
FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
***WARNING: THE SOLUTION APPEARS TO BE DIVERGING. 数值奇异，
估计是边界问题或assemble part后，各个part间没有发生作用，导致
刚度矩阵非正定。 本人做RTM充模过程的模拟也和渗透有关。需要定义一个注入口，从入口按照一定的流量，把多孔介质填充满，由非饱和状态达到饱和状态。请问如何定义注入这一条件？敬请高手指教。 本人也在做岩土分析，
用ADD， REMOVE ELEMENT 出了点问题，
REMOVE SOME ELEMENT 后，
计算结果无变化，感觉奇怪：SCRIPT 如下：
"""
This module constructs and analyzes amodel.
"""
from abaqus import *
from abaqusConstants import *
#形成岩石区域
session.viewports['Viewport: 1'].setValues(displayedObject=None)
session.viewports['Viewport: 1'].setValues(displayedObject=None)
session.viewports['Viewport: 1'].partDisplay.setValues(mesh=OFF)
session.viewports['Viewport: 1'].partDisplay.meshOptions.setValues(
meshTechnique=OFF)
s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__',
sheetSize=200.0)
g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints
s.setPrimaryObject(option=STANDALONE)
s.rectangle(point1=(-30.0, 30.0), point2=(30.0, -30.0))
#s.CircleByCenterPerimeter(center=(0.0, 0.0), point1=(5.0, 0.0))
#s.RadialDimension(curve=g, textPoint=(7.93981409072876, -9.88425731658936),
# radius=5.0)
p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=TWO_D_PLANAR,
type=DEFORMABLE_BODY)
p = mdb.models['Model-1'].parts['Part-1']
p.BaseShell(sketch=s)
s.unsetPrimaryObject()
p = mdb.models['Model-1'].parts['Part-1']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Model-1'].sketches['__profile__']
#形成衬砌
s1 = mdb.models['Model-1'].ConstrainedSketch(name='__profile__',
sheetSize=50.0)
g, v, d, c = s1.geometry, s1.vertices, s1.dimensions, s1.constraints
s1.setPrimaryObject(option=STANDALONE)
s1.CircleByCenterPerimeter(center=(0.0, 0.0), point1=(5.0, 0.0))
s1.RadialDimension(curve=g, textPoint=(6.07638835906982, -4.69907331466675),
radius=5.0)
p = mdb.models['Model-1'].Part(name='Part-2', dimensionality=TWO_D_PLANAR,
type=DEFORMABLE_BODY)
p = mdb.models['Model-1'].parts['Part-2']
p.BaseWire(sketch=s1)
s1.unsetPrimaryObject()
p = mdb.models['Model-1'].parts['Part-2']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Model-1'].sketches['__profile__']
#形成材料
#define material and section
session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=ON,
engineeringFeatures=ON)
mdb.models['Model-1'].Material(name='concrate')
mdb.models['Model-1'].materials['concrate'].Density(table=((2600.0, ), ))
mdb.models['Model-1'].materials['concrate'].Elastic(table=((25000000000.0,
0.3), ))
mdb.models['Model-1'].Material(name='rock')
mdb.models['Model-1'].materials['rock'].Density(table=((2300.0, ), ))
mdb.models['Model-1'].materials['rock'].Elastic(table=((14000000000.0, 0.25),
))
#define section
#beam sectionn
mdb.models['Model-1'].RectangularProfile(name='Profile-1', a=1.0, b=0.6)
mdb.models['Model-1'].BeamSection(name='Section-beam', profile='Profile-1',
integration=BEFORE_ANALYSIS, poissonRatio=0.3, density=3900.0,
referenceTemperature=20.0, thermalExpansion=OFF, temperatureDependency=OFF,
dependencies=0, table=((320000000000.0, 270000000000.0), ),
alphaDamping=0.0, betaDamping=0.0, compositeDamping=0.0, centroid=(0.0,
0.0), shearCenter=(0.0, 0.0), outputPts=((4.0, 4.0), ))
# rock section
mdb.models['Model-1'].HomogeneousSolidSection(name='Section-rock',
material='rock', thickness=1.0)
import regionToolset
p = mdb.models['Model-1'].parts['Part-1']
f = p.faces
faces = f.getSequenceFromMask(mask=('[#1 ]', ), )
region = regionToolset.Region(faces=faces)
p = mdb.models['Model-1'].parts['Part-1']
p.SectionAssignment(region=region, sectionName='Section-rock', offset=0.0)
p1 = mdb.models['Model-1'].parts['Part-2']
session.viewports['Viewport: 1'].setValues(displayedObject=p1)
p = mdb.models['Model-1'].parts['Part-2']
e = p.edges
edges = e.getSequenceFromMask(mask=('[#1 ]', ), )
region = regionToolset.Region(edges=edges)
p = mdb.models['Model-1'].parts['Part-2']
p.SectionAssignment(region=region, sectionName='Section-beam', offset=0.0)
#asambly part
a = mdb.models['Model-1'].rootAssembly
a.DatumCsysByDefault(CARTESIAN)
p = mdb.models['Model-1'].parts['Part-1']
a.Instance(name='Part-1-rock', part=p, dependent=OFF)
#a = mdb.models['Model-1'].rootAssembly
#p = mdb.models['Model-1'].parts['Part-2']
#a.Instance(name='Part-2-concrete', part=p, dependent=OFF)
p = mdb.models['Model-1'].parts['Part-2']
e = p.edges
edges = e.getSequenceFromMask(mask=('[#1 ]', ), )
region=regionToolset.Region(edges=edges)
p = mdb.models['Model-1'].parts['Part-2']
p.assignBeamSectionOrientation(region=region, method=N1_COSINES, n1=(0.0, 0.0,
-1.0))
#-------------------------------------------------------
import step
# Create a step. The time period of the static step is 1.0,
# and the initial incrementation is 0.1; the step is created
# after the initial step.
#a = mdb.models['Model-1'].rootAssembly
mdb.models['Model-1'].StaticStep(name='Step-grave',previous='Initial',
timePeriod=1.0, initialInc=0.0001,
description='Load the top of the beam.')
mdb.models['Model-1'].steps['Step-grave'].setValues(nlgeom=ON)
mdb.models['Model-1'].StaticStep(name='Step-grave2',previous='Step-grave',
timePeriod=1.0, initialInc=0.0001,
description='Load the top of the beam.')
mdb.models['Model-1'].steps['Step-grave2'].setValues(nlgeom=ON)
#output
#mdb.models['Model-1'].fieldOutputRequests['F-Output-cal'].setValues(variables=(
# 'CDISP', 'CF', 'CSTRESS', 'LE', 'PE', 'PEEQ', 'PEMAG', 'RF', 'S', 'U' ))
mdb.models['Model-1'].FieldOutputRequest(name='F-Output-1',
createStepName='Step-grave', variables=('S', 'MISESMAX', 'TSHR', 'CTSHR',
'ALPHA', 'VS', 'PS'))
#create curve setx
#a = mdb.models['Model-1'].rootAssembly
#e1 = a.instances['Part-1-rock'].edges
#edges1 = e1.getSequenceFromMask(mask=('[#1 ]', ), )
#a.Set(edges=edges1, name='Set-rock-curvx')
#: The set 'Set-rock-curv' has been created (1 edge).
#a = mdb.models['Model-1'].rootAssembly
#e1 = a.instances['Part-2-concrete'].edges
#edges1 = e1.getSequenceFromMask(mask=('[#1 ]', ), )
#a.Set(edges=edges1, name='Set-concrete-curvx')
#: The set 'Set-concrete-curv' has been created (1 edge).
#mdb.models['Model-1'].Equation(name='Constraint-r-cx', terms=((1.0,
# 'Set-rock-curvx', 1), (-1.0, 'Set-concrete-curvx', 1)))
#create curve sety
#a = mdb.models['Model-1'].rootAssembly
#e1 = a.instances['Part-1-rock'].edges
#edges1 = e1.getSequenceFromMask(mask=('[#1 ]', ), )
#a.Set(edges=edges1, name='Set-rock-curvy')
#: The set 'Set-rock-curv' has been created (1 edge).
#a = mdb.models['Model-1'].rootAssembly
#e1 = a.instances['Part-2-concrete'].edges
#edges1 = e1.getSequenceFromMask(mask=('[#1 ]', ), )
#a.Set(edges=edges1, name='Set-concrete-curvy')
#: The set 'Set-concrete-curv' has been created (1 edge).
#mdb.models['Model-1'].Equation(name='Constraint-r-cy', terms=((1.0,
# 'Set-rock-curvy', 2), (-1.0, 'Set-concrete-curvy', 2)))
a = mdb.models['Model-1'].rootAssembly
e1 = a.instances['Part-1-rock'].edges
edges1 = e1.getSequenceFromMask(mask=('[#8 ]', ), )
region = regionToolset.Region(edges=edges1)
mdb.models['Model-1'].EncastreBC(name='BC-fixbot', createStepName='Initial',
region=region)
mdb.saveAs(pathName='C:/Temp/yy.cae')
#: The model database has been saved to "C:\Temp\yy.cae".
session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Step-grave')
a = mdb.models['Model-1'].rootAssembly
s1 = a.instances['Part-1-rock'].edges
side1Edges1 = s1.getSequenceFromMask(mask=('[#2 ]', ), )
region = regionToolset.Region(side1Edges=side1Edges1)
mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-grave',
region=region, distributionType=UNIFORM, field='', magnitude=1000.0,
amplitude=UNSET)
#session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Step-seismic')
#mdb.models['Model-1'].Gravity(name='Load-2', createStepName='Step-grave',
# comp2=0.8)
#CREATE PARTITION
#session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=ON)
#session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues(
# meshTechnique=ON)
#a = mdb.models['Model-1'].rootAssembly
#f1 = a.instances['Part-1-rock'].faces
#pickedFaces = f1.getSequenceFromMask(mask=('[#1 ]', ), )
#e1 = a.instances['Part-1-rock'].edges
#a.PartitionFaceByShortestPath(faces=pickedFaces,
# point1=a.instances['Part-1-rock'].InterestingPoint(edge=e1,
#rule=MIDDLE), point2=a.instances['Part-1-rock'].InterestingPoint(
# edge=e1, rule=MIDDLE))
#a = mdb.models['Model-1'].rootAssembly
#f1 = a.instances['Part-1-rock'].faces
#pickedFaces = f1.getSequenceFromMask(mask=('[#3 ]', ), )
#e11 = a.instances['Part-1-rock'].edges
#a.PartitionFaceByShortestPath(faces=pickedFaces,
# point1=a.instances['Part-1-rock'].InterestingPoint(edge=e11,
#rule=MIDDLE), point2=a.instances['Part-1-rock'].InterestingPoint(
#edge=e11, rule=MIDDLE))
#attribute mesh type
import mesh
elemType1 = mesh.ElemType(elemCode=CPS4R, elemLibrary=STANDARD,
secondOrderAccuracy=OFF, hourglassControl=DEFAULT,
distortionControl=DEFAULT)
elemType2 = mesh.ElemType(elemCode=CPS3, elemLibrary=STANDARD)
a = mdb.models['Model-1'].rootAssembly
f1 = a.instances['Part-1-rock'].faces
faces1 = f1.getSequenceFromMask(mask=('[#1 ]', ), )
pickedRegions =(faces1, )
a.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2))
a = mdb.models['Model-1'].rootAssembly
e1 = a.instances['Part-1-rock'].edges
pickedEdges = e1.getSequenceFromMask(mask=('[#f ]', ), )
a.seedEdgeByNumber(edges=pickedEdges, number=8)
a = mdb.models['Model-1'].rootAssembly
partInstances =(a.instances['Part-1-rock'], )
a.generateMesh(regions=partInstances)
a = mdb.models['Model-1'].rootAssembly
a.PartFromInstanceMesh(name='Part-1-mesh-1')
p1 = mdb.models['Model-1'].parts['Part-1-mesh-1']
session.viewports['Viewport: 1'].setValues(displayedObject=p1)
#create kill node
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
n = p.nodes
nodes = n.getSequenceFromMask(mask=('[#c0000000 #381 ]', ), )
p.Set(nodes=nodes, name='Set-kill-node')
#: The set 'Set-kill-node' has been created (6 nodes).
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
e = p.elements
elements = e.getSequenceFromMask(mask=('[#ffffffff:2 ]', ), )
p.Set(elements=elements, name='Set-all-element')
#: The set 'Set-all-element' has been created (64 elements).
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
region = p.sets['Set-all-element']
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
p.SectionAssignment(region=region, sectionName='Section-rock', offset=0.0)
a1 = mdb.models['Model-1'].rootAssembly
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
a1.Instance(name='Part-rock-mesh', part=p, dependent=ON)
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p = mdb.models['Model-1'].parts['Part-1-mesh-1']
e = p.elements
elements = e.getSequenceFromMask(mask=('[#3c3c0000 #3c3c ]', ), )
p.Set(elements=elements, name='Set-kill-e')
#: The set 'Set-kill-e' has been created (16 elements).
mdb.models['Model-1'].keywordBlock.synchVersions(storeNodesAndElements=FALSE)
mdb.models['Model-1'].keywordBlock.insert(59, """
*MODEL CHANGE ,remove
Part-rock-mesh.Set-kill-e""")
#Part-rock-mesh.Set-all-element""")
#Part-rock-mesh.Set-kill-e""")
jobName = 'tun-exp'
a = mdb.models['Model-1']
myJob = mdb.Job(name=jobName, model=a,
description='tunneltutorial') 我以前提出的问题解决了，
有兴趣的来交流。
地震分析应力平衡如何实现？！！
在做地震分析。提一个高难度的：
地震边界：由静力到动力设置粘弹性边界，
如何保证在去除以前静力边界后，保持静力分析时候得到的力不变，以便维持其初始平衡，
否则，即使动力分析没有力，也会出现很大变形？ wxyz1234 你好.我也是做地震作用下结构动力分析的, 终于找到一个做地震方面的了. 如果不介意加我QQ402802221 标注:aba 地震
页:
1
[2]