【寻找】请教这个程序有什么问题

luckier

一个简单的模型
算了好几变了,还是有问题,实在是找不出问题所在,好打击人呀
希望各位给我看一下,指点以下,多谢
  
*heading
example
*NODE
         1,  0.0            ,  0.0            
         2,  1.0            ,  0.0            
         3,  0.0            ,  1.0            
         4,  1.0            ,  1.0            
         6,  0.1            ,  0.0            
         7,  0.0            ,  0.1            
        10,  0.0            ,  0.1            
        11,  0.0257621358799,  0.0961816697772
        12,  0.0497941299468,  0.0862272919038
        13,  0.0707106781186,  0.0707106781186
        14,  0.0862272919038,  0.0497941299468
        15,  0.0961816697772,  0.0257621358799
        16,  0.1            ,  0.0            
        17,  0.173125       ,  0.0            
        18,  0.2575         ,  0.0            
        19,  0.353125       ,  0.0            
        20,  0.46           ,  0.0            
        21,  0.578125       ,  0.0            
        22,  0.7075         ,  0.0            
        23,  0.848125       ,  0.0            
        24,  1.0            ,  0.0            
        25,  1.0            ,  0.1111111111111
        26,  1.0            ,  0.2444444444444
        27,  1.0            ,  0.4            
        28,  1.0            ,  0.5777777777777
        29,  1.0            ,  0.7777777777777
        30,  1.0            ,  1.0            
        31,  0.7777777807081,  1.0            
        32,  0.5777777833455,  1.0            
        33,  0.4000000079121,  1.0            
        34,  0.2444444544078,  1.0            
        35,  0.1111111228327,  1.0            
        36,  0.0            ,  1.0            
        37,  0.0            ,  0.848125        
        38,  0.0            ,  0.7075         
        39,  0.0            ,  0.578125        
        40,  0.0            ,  0.46            
        41,  0.0            ,  0.353125        
        42,  0.0            ,  0.2575         
        43,  0.0            ,  0.173125        
        44,  0.4142802531718,  0.4142802498529
        45,  0.4845821062629,  0.2628824764996
        46,  0.5384513847209,  0.1242631825913
        47,  0.3283178035596,  0.3283178005499
        48,  0.2520739467954,  0.2520739443222
        49,  0.1832591848984,  0.1832591830869
        50,  0.1227436241303,  0.1227436230457
        51,  0.3874069231229,  0.2094736286703
        52,  0.4304726737686,  0.0999440778324
        53,  0.298372138405 ,  0.1632310573885
        54,  0.3323788546297,  0.0784154050604
        55,  0.2187871874063,  0.1208968835849
        56,  0.2436235068829,  0.0591955803788
        57,  0.1479269486158,  0.0829198873516
        58,  0.1648128001739,  0.0415664973921
        59,  0.5803488838486,  0.4216958396939
        60,  0.7751371056415,  0.4173295909556
        61,  0.6340881587775,  0.2653431534298
        62,  0.8054453906697,  0.2592109064973
        63,  0.6750178626847,  0.1239903687951
        64,  0.8289436717362,  0.1195843272868
        65,  0.4216958441944,  0.580348881354  
        66,  0.4173295968308,  0.7751371042456
        67,  0.5886615726277,  0.5886615693598
        68,  0.5878210434346,  0.780918887149  
        69,  0.7809188889075,  0.5878210392453
        70,  0.7817969257914,  0.7817969235756
        71,  0.2628824812809,  0.4845821042643
        72,  0.1242631889068,  0.5384513836107
        73,  0.2653431596023,  0.6340881572403
        74,  0.1239903766265,  0.6750178617407
        75,  0.2592109142204,  0.8054453897403
        76,  0.119584336736 ,  0.8289436710466
        77,  0.2094736327812,  0.3874069212722
        78,  0.1632310606243,  0.2983721368217
        79,  0.1208968858581,  0.218787186186  
        80,  0.0829198885868,  0.1479269478035
        81,  0.0999440830026,  0.4304726726852
        82,  0.0784154086439,  0.332378853664  
        83,  0.0591955824717,  0.2436235060875
        84,  0.0415664983286,  0.1648127995317
*ELEMENT,TYPE=CPE4I,ELSET=body_2d
         1,        44,        47,        51,        45
         2,        45,        51,        52,        46
         3,        46,        52,        20,        21
         4,        47,        48,        53,        51
         5,        51,        53,        54,        52
         6,        52,        54,        19,        20
         7,        48,        49,        55,        53
         8,        53,        55,        56,        54
         9,        54,        56,        18,        19
        10,        49,        50,        57,        55
        11,        55,        57,        58,        56
        12,        56,        58,        17,        18
        13,        50,        13,        14,        57
        14,        57,        14,        15,        58
        15,        58,        15,        16,        17
        16,        44,        45,        61,        59
        17,        59,        61,        62,        60
        18,        60,        62,        26,        27
        19,        45,        46,        63,        61
        20,        61,        63,        64,        62
        21,        62,        64,        25,        26
        22,        46,        21,        22,        63
        23,        63,        22,        23,        64
        24,        64,        23,        24,        25
        25,        44,        59,        67,        65
        26,        65,        67,        68,        66
        27,        66,        68,        32,        33
        28,        59,        60,        69,        67
        29,        67,        69,        70,        68
        30,        68,        70,        31,        32
        31,        60,        27,        28,        69
        32,        69,        28,        29,        70
        33,        70,        29,        30,        31
        34,        44,        65,        73,        71
        35,        71,        73,        74,        72
        36,        72,        74,        38,        39
        37,        65,        66,        75,        73
        38,        73,        75,        76,        74
        39,        74,        76,        37,        38
        40,        66,        33,        34,        75
        41,        75,        34,        35,        76
        42,        76,        35,        36,        37
        43,        44,        71,        77,        47
        44,        47,        77,        78,        48
        45,        48,        78,        79,        49
        46,        49,        79,        80,        50
        47,        50,        80,        12,        13
        48,        71,        72,        81,        77
        49,        77,        81,        82,        78
        50,        78,        82,        83,        79
        51,        79,        83,        84,        80
        52,        80,        84,        11,        12
        53,        72,        39,        40,        81
        54,        81,        40,        41,        82
        55,        82,        41,        42,        83
        56,        83,        42,        43,        84
        57,        84,        43,        10,        11
*SOLID SECTION, ELSET=body_2d, MATERIAL=stell
*NSET, NSET=fix1
      10,      36,      37,      38,      39,      40,
      41,      42,      43
*NSET, NSET=fix2
      16,      17,      18,      19,      20,      21,  
      22,      23,      24
*NSET, NSET=force
      30,      31,      32,      33,      34,      35,
      36
*MATERIAL, NAME=stell
*ELASTIC, TYPE = ISOTROPIC
30000000.0,0.3       ,0.0        
*STEP
*STATIC
*CLOAD
force,2,1.0
*BOUNDARY
fix1,1,1,0.0
fix2,2,2,0.0
*EL FILE
SINV
*NODE FILE
U
*END STEP

chenhu

也是刚学abaqus,导入cAE看了一下，好像是你的边界条件出了问题，在abaqus默认的initial步中（先于step-1），也你将两条边上的结点（fix1，fix2）约束就应该存在！而你命令中，是将*boundary放入step中，所以可能出现问题！试着将其提前，再算一下试试吧！不知大侠们的意见呢？刚学，有错误请不要笑！对的话，赏个分先！：L)

chenhu

做了出来，看一下是不是你要的结果！：）

wjytj

chenhu wrote:
也是刚学abaqus,导入cAE看了一下，好像是你的边界条件出了问题，在abaqus默认的initial步中（先于step-1），也你将两条边上的结点（fix1，fix2）约束就应该存在！而你命令中，是将*boundary放入step中，所以可能出现问题！试着将其提前，再算一下试试吧！不知大侠们的意见呢？刚学，有错误请不要笑！对的话，赏个分先！：L)

  
不用改也可以运行，结果和你的一样。
在ABAQUS中，不改变的固定边界条件中最好在Initial步中加，但在荷载步中加也是可以的。

luckier

chenhu大侠,你的意思是把*boundary那到*step外面,如下吗?
*step在这里不是用于指明分析类型吗那需要把力也拿出来呢？
  
*BOUNDARY  
fix1,1,1,0.0  
fix2,2,2,0.0  
*STEP  
*STATIC  
*CLOAD  
force,2,1.0  
*EL FILE  
SINV  
*NODE FILE  
U  
*END STEP

chenhu

1.wjytj 老大指教的对，又学了点东西！但我有个疑问：如果不加在step外或initial中，像你说的加在step中，那么是不是多个step中对于此边界条件都要定义一次还是只在第一个step中就可以了呢？另外，如果我在后面的step中要改动边界条件（比如说取消此边条，怎么办泥？）在caE中好做，但命令还请老大指点一下！
2.to luckier，力不是边条，当然不能拿到外面吧！step是指定分析类型，但boundary也是呀，各司其职，我还是老老实实那么写吧！：）

cafe0519

chenhu wrote:
1.wjytj 老大指教的对，又学了点东西！但我有个疑问：如果不加在step外或initial中，像你说的加在step中，那么是不是多个step中对于此边界条件都要定义一次还是只在第一个step中就可以了呢？另外，如果我在后面的step中要改动边界条件（比如说取消此边条，怎么办泥？）在caE中好做，但命令还请老大指点一下！
  

  
*BOUNDARY OP=new or mod:)

luckier

哦＊boundary的
这个问题wjytj 老大以前给我批评过（又反出以前的帖子看了），当时还不是很清楚，现在明白了
不过着段程序运行的结果是chenhu那样吗？看来我的后处理出了一点毛病
*BOUNDARY OP默认是mod把，所以这个＊step中定义了，如果下会
还是这个条件，就不许要修改．如果要删除所有的用new，如果
删除部分（比较少），可能还需要用mod方便，是不是这样

chenhu

我想也应该是这样，对应于cae中：new---created，mod---modified！

cafe0519

Set OP=MOD (default) to modify existing boundary conditions or to add boundary conditions to degrees of freedom that were previously unconstrained.
  
Set OP=NEW if all boundary conditions that are currently in effect should be removed. To remove only selected boundary conditions, use OP=NEW and respecify all boundary conditions that are to be retained.