声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: hero2005

(SOS)有用ABAQUS做岩土工程的吗,大家交流一下?

[复制链接]
发表于 2007-1-25 11:34 | 显示全部楼层

回复 #4 zhangmeng 的帖子

现在我想做一个碎石桩复合地基的数值模拟分析,但是存在下面的困难,首先,在装配时,Assembly具体的意思是什么?不同的部件在装配时,如果各个部件的网格分割的不相等,在没有设置接触单元时,这几个部件能认为它们已经连在一起了吗?如果不行的话,请问大哥:我应该怎么做?
回复 支持 反对
分享到:

使用道具 举报

发表于 2007-1-25 16:14 | 显示全部楼层

求教:

我现在做一个碎石桩复合地基的数值模拟,现在出了点问题,下面是一些警告:
***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 NUMBER  1, 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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         1.107E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.107E+05   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         2.215E+05   AT NODE          2   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.107E+05   AT NODE         33   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         3.322E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.107E+05   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         4.430E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.107E+05   AT NODE          3   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         5.537E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.107E+05   AT NODE         17   DOF  2
   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 NUMBER  2, TIME INCREMENT  0.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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         2.769E+04   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        2.769E+04   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         5.537E+04   AT NODE          2   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        2.769E+04   AT NODE         33   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         8.306E+04   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        2.769E+04   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         1.107E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        2.769E+04   AT NODE          3   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         1.384E+05   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        2.769E+04   AT NODE         17   DOF  2
   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 NUMBER  3, TIME INCREMENT  6.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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         6.922E+03   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        6.922E+03   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         1.384E+04   AT NODE          2   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        6.922E+03   AT NODE         33   DOF  2
   INSTANCE: CZB1-1                                                                          
          FORCE     EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.

***WARNING: SOLVER PROBLEM.  NUMERICAL SINGULARITY WHEN PROCESSING NODE 22

这是什么意思?请各位大哥帮帮忙?到底问题是什么造成的?本人急需要解决的问题,谢谢!
发表于 2007-1-25 16:15 | 显示全部楼层

求教:

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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         2.076E+04   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        6.922E+03   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         2.769E+04   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        6.922E+03   AT NODE          3   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         3.461E+04   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        6.922E+03   AT NODE         17   DOF  2
   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 NUMBER  4, TIME INCREMENT  1.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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         1.730E+03   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.730E+03   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         3.461E+03   AT NODE          2   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.730E+03   AT NODE         33   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         5.191E+03   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.730E+03   AT NODE          4   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         6.922E+03   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.730E+03   AT NODE          3   DOF  2
   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   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST INCREMENT OF DISP.         8.652E+03   AT NODE          4   DOF  2
   INSTANCE: CZB1-1                                                                          
LARGEST CORRECTION TO DISP.        1.730E+03   AT NODE         17   DOF  2
   INSTANCE: CZB1-1                                                                          
          FORCE     EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.

***WARNING: THE SOLUTION APPEARS TO BE DIVERGING.
发表于 2007-1-25 16:40 | 显示全部楼层
数值奇异,
估计是边界问题或assemble part后,各个part间没有发生作用,导致
刚度矩阵非正定。
发表于 2007-2-1 15:48 | 显示全部楼层
本人做RTM充模过程的模拟也和渗透有关。需要定义一个注入口,从入口按照一定的流量,把多孔介质填充满,由非饱和状态达到饱和状态。请问如何定义注入这一条件?敬请高手指教。
发表于 2007-2-17 10:50 | 显示全部楼层
本人也在做岩土分析,
用ADD, REMOVE ELEMENT 出了点问题,
REMOVE SOME ELEMENT 后,
计算结果无变化,感觉奇怪:SCRIPT 如下:

"""
This module constructs and analyzes a  model.
"""


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[6], 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[2], 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[1],
  #  rule=MIDDLE), point2=a.instances['Part-1-rock'].InterestingPoint(
   # edge=e1[3], 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[9],
  #  rule=MIDDLE), point2=a.instances['Part-1-rock'].InterestingPoint(
  #  edge=e11[2], 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='tunnel  tutorial')
发表于 2007-3-5 21:54 | 显示全部楼层
我以前提出的问题解决了,
有兴趣的来交流。
发表于 2007-3-5 22:02 | 显示全部楼层

地震分析应力平衡如何实现?!!

在做地震分析。
提一个高难度的:
地震边界:由静力到动力设置粘弹性边界,
如何保证在去除以前静力边界后,保持静力分析时候得到的力不变,以便维持其初始平衡,
否则,即使动力分析没有力,也会出现很大变形?
发表于 2007-3-8 20:56 | 显示全部楼层
wxyz1234 你好.我也是做地震作用下结构动力分析的, 终于找到一个做地震方面的了. 如果不介意加我QQ402802221 标注:aba 地震
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-15 01:38 , Processed in 0.075259 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表