论坛汇总第一部分
另外，seansheu 大侠，我按你的方法将非线性的开关打开后，结果就正常了，可是我不明白的是既然我已经将混凝土的塑性曲线数据按点输入了，肯定就不是线性的了，为什么还要另外加上非线性开关它才明白呢。
几何非线性即通常所说的大变形分析，是指系统的微分方程等都建立在结构当前变形基础上的；不打开NLgeom即意味着系统的控制方程等都建立结构的初始构型状态。
嗯，明白，就是说控制方程是采用前一步的结果作为起始量，而不是象线性可以采用初始值。
關於willams_lee的問題，我重新回去檢視這個問題，當初因為初始定義用的是c3d8r單元，當線性錯誤時，直覺反應是認為沒開啟nlgeom來搭配材料非線性。但是當我回去確認應變時，發現它也沒大到那麼離譜...我切開模型，發現了致命的錯誤，鋼筋中心截面應力值不合理，變形分佈也不合理。顯然先前我作的假設有誤，所以計算錯誤...因為它網格密度並非充足，且有些單元品質不好，所以採用縮減積分的c3d8r單元計算會有問題，因此我改用了c3d8I，重新計算。此時，nlgeom開啟與關閉的結果都一樣。結果顯示，鋼筋中心截面應力值結果較為合理，變形分佈也合理，所以用c3d8r可能才是錯誤原因。非關幾何非線性...當然，我還要提高網格密度，去確認結果的收斂性。對willams_lee的誤導，深感抱歉。也感謝willams_lee的提問，讓我有機會發現錯誤。
- 5 -
论坛汇总第一部分
在有限元分析中，应该尽量避免采用一次减缩积分单元（否则应该加密网格，有时也不一定管用）；而应该采用高次积分单元。 C3D8I是非协调单元，和高次单元一样，都是为了防止或者沙漏现象的影响。估计你采用同样的网格密度，而选用C3D20单元，结果也应该正确。类似问题在这里也曾讨论过：另外，个人觉得，既然钢筋和混凝土之间的摩擦不是考察的重点，也就是说不考虑钢筋和混凝土之间的粘结滑移，那么，这个问题完全没有必要如此复杂，定义混凝土一个Part（或者在CAE中组装成Assembly），然后直接利用*Rebar定义钢筋即可，既可以定义为单独的钢筋，也可定义为钢筋层。在6.3版本中，这一功能非常强大（比较予6.2而言），Benchmark Manual和Example Manual中有相应的例子可供参考。
谢谢回应，现附上该文件。
74768-Beam-deform.inp (0.0k)
內存不足，應該修改 abaqus_v6.env檔中的內存設定，加大 pre_memory值或是standard_memory值，端視你是在 前處理 或是 求解 時出現內存不足之警訊，再擇一放大。通常，abaqus中之內定值都會過小，碰到大問題時常會出現內存不足問題BTW,你的input file不完整，無法執行...
输入文件因为包含太多节点信息，足足有六七兆大，没办法只好删去了很多节点的信息，用省略号代替的。内存不足我也想过，不过我觉得还是存在我接触以及约束定义不正确的可能，否则不会一个简单的问题就需要如此多的内存，不知道这样的理解对不对。因为我今天下午修改取消了所有约束，以及接触的摩擦定义，就很快正常运行，不过结果是混凝土梁被压变形，钢筋跑了出来，这显然是不符合实际的，但是运行还是正常，真不知道怎么去定义。
sorry, 我不是很懂兆的單位，英文是指 MB 還是 GB呢?如果可以，你可以壓縮檔案，因為input file是ASCII格式，壓縮比可以很大。其實，你可先照我的方式調整 abaqus_v6.env 中內存設定值，把它調整至abaqus的建議值，如此應該可行...此外，忠於你問題的真實情況假設才是重點，怎可因為內存不足而因噎廢食...如果為了遷就內存，而胡亂假設問題，豈不GIGO!!所以你可能還要先好好想清楚你的分析問題，一點淺見供你參考...
若我没记错，它应该是以实体内存为优先，若不足，再以硬盘来当虚拟内存你用的是学校服务器也无妨，首先你将 abaqus_v6.env 拷贝一份至你的工作目录，然后将拥有权限设为自己，修改它的内存设定，就可以执行了通常在工作站下，abaqus 对环境设定档(即abaqus_v6.env)的读取优先数序为：工作目录下的 env file > User目录下的 env file > ABAQUS site 目录下的 env file
上述前两者你可以修改，最后者只有 Administrator 才能修改
- 6 -
论坛汇总第一部分
地应力问题
一、在桩土接触分析中，地应力步中，老是提示我this step can't continue by the select increment .我只把*geostatic改成*static就能够计算到30多个increment.为什么会出现这种问题？请各位大虾赐教。 你的geostatic使用不恰当,因为在有接触面的情况下,要注意力的方向。
楼上大哥，我还是不太明白，我的inp文件是由hypermesh生成的，既然方向不对，应该用*static也不能计算阿，请说的明白一点，非常感谢，我现在把我的inp文件传上，望赐教。
268204-1.inp (0.0k)
二、不平地面初始地应力的处理方法
245572-geo3.ps 命令流文件
245574-geo3.inp (0.0k) 我看ABAQUS的帮助里，初始地应力好像是通过给定两点的高程，及相应的z方向应力值、水平方向应力系数来计算单元集的初始地应力，这样对于地面不是水平的怎么办？能不能通过外部文件给单元赋初始地应力,如何操作 可以！！通过inp文件
还有，如果按重力计算出来的应力作为初始应力，请问侧向压力系数怎么体现？急！ 请看本人通过文件进行初始应力场平衡的工程事例（也是刚刚才发出的）吧！ 我有这样一个想法：
我们在一些情况下有了地应力，但是，我们对地表的初始变形不满意，那么我们先在第一步计算出地应力合初始变形，然后我们将地应力写入为一个文件，在这个文件里我们得到单元积分点上保持平衡的应力（或者其他，我还在想）然后我们通过USER SUBROUTINE 或者INPUT读入，那样我们不就得到了应力却没有如何初始变形吗！网上也有这样的帖子，但是过程不是很明了。
我看了帮助，INPUT好像不行，因为它输入的是单元号和几个应力，然后ABAQUS将这几个应力平均的施加到单元的面上。 USER SUBROUTINE 好像可以它输入的是单元号，积分点号，然后就是应力分量。而我们可以在第一步就输出这几个量，只需稍微操作一下，加一些逗号就可以。只是还有一个问题，USER SUBROUTINE 一次只能返回一个单元号一个积分点的应力分量，但是我们需要一批单元的应力分量，暂时还不知道如何处理！！！不知有无大侠知道如何处理！！！ 纠正一点，不是INPUT，是指忽略GEOSTATIC,而直接用*initial conditions,type=stress输入应力分量。 很好的想法，我也在想通过子程序施加，但是ABAQUS的例子太简单，看完 不知所云，：）
我仔细看了帮助，现在面对的难题是：SIGINI里输入的是一个ELEMENT的SRESS，注意是一个ELEMENT而不是一个ELEMENT的一个积分点的STRESS，但是我找了一下，好像找不到一个ELEMENT的STRESS输出，不知兄台有无办法输出！！！如果能找到，就应该能实现目标！
努力值得赞赏，但是你的物质的弹模过高，土的弹模一般达不到GPa的量级，而且如果采用弹塑性本构后引入非线性，收敛情况将会恶化，呵呵
假定你计算的是个土质边坡，假设土的粘结力为零，内摩擦角为20，那么此时边坡的安全系数为1.2。和想和原作者一起探讨边坡稳定性分析的问题 请教zhangpeiwen 兄： 你的INP文件里 *GEOSTATIC 1E-6,1E-6
这个1E-6是什么参数啊？我在HELP里看到*GEOSTATIC里是没有参数的。 是初始时间增量，和最终时间，一些例题里采用这样的参数，应该没有问题．
如果土的粘结力c=0的话，首先采用的屈服准则要改变：好象应该用cam-clay准则（可以看一下帮助确认一下）．我做地应力的目的主要是想考虑有水位变化作用下（例如库水位骤升，骤降）边坡的稳定问题，实际上是耦合问题．这些天一直在摸索采用abaqus来模拟水位的骤升、降问题。模拟结果不是很理想。不知道你有什么好的处理方法。请指教。多谢 我有这样一个想法：
- 7 -
论坛汇总第一部分
我们在一些情况下有了地应力，但是，我们对地表的初始变形不满意，那么我们先在第一步计算出地应力合初始变形，然后我们将地应力写入为一个文件，在这个文件里我们得到单元积分点上保持平衡的应力（或者其他，我还在想）然后我们通过USER SUBROUTINE 或者INPUT读入，那样我们不就得到了应力却没有如何初始变形吗！网上也有这样的帖子，但是过程不是很明了。
但是我现在还不明白生成这个初始应力时应该生成哪一些结果，INPUT里的格式和生成的格式是否相同。这些都是问题 三、【讨论】初始地应力的施加问题 我们在分析一个边坡或者其他的岩土问题的时候，一定要考虑初始地应力问题的，比如边坡稳定性分析问题，一个自然边坡在我们开挖之前是出于自然平衡状态的，但是我们在模拟的时候，给它一个初始地应力的时候，它就不出于平衡状态了，和自然中的实际情况时不一样的。岩体中的应力为荷载产生的应力和初始地应力的和，因此初始地应力的大小直接决定我们的计算结果，如何正确的模拟这个问题到现在还没有很好的妥善的处理。由于初始地应力的存在，岩体的开挖将导致部分岩体卸载，但是在某些部位却是加载的，这个问题也导致了我们在模拟时的难度，也产生如何正确模拟这个过程的问题。然而通常我们都是这样处理的：沿开挖面作用着与初始地应力等价的“开挖释放荷载”的方式来模拟，把由此所得的位移即作为由于开挖产生的岩体位移，由此所得的应力场同初始应力场迭加即为开挖后的应力场。这种方法其实质就是等价原则的实施而已。我的问题的就是对上面阐述在我们的ABAQUS中的实现问题，我们在计算的时候总是要先 施加初始地应力的，当我们只是在这个力的作用下进行计算，我们发现我们的计算对象已经发生了严重的变相，这和我们的实际情况是不相同的，这个问题我们不管是在使用什么软件的时候都出现过，由此我本人请教各位高手对这个问题的理解是怎么样的？我们在实现模拟计算的时候如何处理这个问题？还有就是初始地应力的作用是如何体现的？对我指教一二，本人将不胜感激，我想这个问题也许是很多人想知道的吧。谢谢！！！
呵呵，zhuxunguo所说问题确实存在，而且也没有很好的处理方法。看经验和问题的重要程度了，化些功夫还是可以做的比较满意的。通常的初始地应力场的施加（最简单的就是垂直向加以自重，水平向加以地质勘测给出的侧压力系数乘以自重）常常要引起网格的变形，而这个变形在事实上天然状态下是不存在的，因为我们以天然状态下已经发生变形的尺寸数据来做造型，现在让它再变形一次，显然引起我们的造型网格尺寸的变化，比较恼火呀。一次最好就是探讨探讨怎么能即加进初始应力场，又不产生变形，保证后续计算的正确性。 欢迎一起讨论。。
感谢楼上的小姐，我看了你告诉我的这个网址。受益匪浅，但是还是不是完全的理解。所谓初始地应力是时间和空间的函数。因此我们可以用场这个概念来反映。也就是初始地应力场。初始应力一般情况是由下面的几部分组成的：重力，构造力，动水压力，温度应力和结晶力等。在岩体或者土体开挖中破坏了原来的应力平衡，引起周围岩体的应力重新分布，形成一个新的应力场，称为二次应力场。它与工程断面是由关系的，并随开挖工作的进行再时间上，空间上不断变化。
上面是初始地应力和二次地应力的概念，从他们的概念中我们可以发现，初始地应力破坏后就会形成二次地应力，也就是我们为什么在施加了初始地应力后计算的边坡和隧道等支护体为什么破坏的原因。但是事实上再很多的情况下，我们开挖后支护体并没有破坏，这个原因是我需要知道的，在你告诉我的网址上大家并没有得出很好的结论，而且有很多的理解，尤其是对初始地应力和二次地应力的理解都不正确怎么可能知道这个原因呢，希望你能对我的提问给与我帮助，谢谢。也许我的理解有所偏颇，还请指点一二。谢谢
感谢斑竹的鼓励，我的觉得和你们这些高手交流实在是受益匪浅。可能大家把的意思理解不是很对，我的意思是在简化的情况下，当然关键的因素都要考虑到不然我们计算的意义就失去了，斑竹你说呢？比如我们对材料使用DRUCKER--PRAGER 屈服准则，模型使用的摩尔库仑模型的时候，我们也一般都是使用这样的模型和屈服准则，当然由于材料的不同也有使用邓肯张模型的，剑桥模型以及剑桥修改模型的，都是必须要考虑初始地应力的问题，然后再施加我们要给它施加的载荷，以求解材料的受力情况。所以我提问的目的是为什么我们在施加初始地应力计算（先不施加二次载荷）的时候，材料发生了变形，这和我们实际情况中是不同的，因为一个边坡工程在我们开挖之前是达到平衡状态的，为什么我们计算的时候会出现这个现象呢？我不是很明白。所以提问，请教谢谢！！！
我曾经计算过水工地下厂房，也涉及到原始地应力场的问题。当时做的时候好像是两种处理办法。当然首先要做回归计算，根据测点的应力数据，根据问题的性质假定机组力，用最小二乘法拟合得到计算模型的边界条件。一是将此边界条件施加到计算模型上计算一次，得到的应力场就作为初始地应力场，当然会出现变形，一般我们回忽略这个变形引起的造型误差，后续计算由于边界改变引起的应力场变化叠加到初始应力场上就是二次应力场，每次计算位移要减掉第一次计算得到的位移作为由于条件改变（如开挖、荷载变化）引起的位移。第二种方法是在模型上同时施加应力场和相应边界条件，由于内力与外力平衡，不会有位移产生，但是这样加好像麻烦一些。
- 8 -