Abaqus-使用中的一些经验.

论坛汇总第一部分Abaqus 使用点滴Part IIFrom Aoxue、SimWe、OkOk etc2005年3月27日- 55 -一个Abaqus CAE分析后问题讨论我目前在做一个很简单的钢筋混凝土梁模型，定义的part分别是混凝土梁，以及三根钢筋，assembly以后开始定义intersection，不知道怎么回事，我定义了混凝土与钢筋之间的摩擦关系，照理说应该不用约束条件了吧。mesh以后运行工作，老是提示内存仅规定256兆，最少需要1029兆，晕，我减少了mesh的element数量，提示还是如此，后来考虑是三根钢筋与混凝土的摩擦关系太复杂，将其中两根变为无摩擦，剩一根有摩擦，这时候运行还是提示出错，不过内存要求减少为417兆，我最后将那根钢筋也变成无摩擦，就没有任何作用，依旧提示内存不足。实在受不了了，不知道哪位高手遇到过这样的问题，或者做过混凝土与钢筋这样的程序，帮忙指教一下，多谢多谢了多谢 seansheu 大侠，摩擦力确实并非我观察的重点，因为我是初学，就先自己想了这么一个自认为比较简单的例子来练练手，结果被彻底搞惨了，我之所以加上摩擦是因为我想不到别的办法可以带动两者一起运动。摩擦系数我需要调整，另外，荷载以及混凝土的一些参数同样需要调整，我会重新检 查一次，按实际的实验数据输入模型看看，然后在发上来和大家讨论。哈.如果摩擦力非您的觀察重點，那您的問題就簡單了.你的模型可視為四個體積，網格連續，材料性質相異如此問題就變為無摩擦的靜力問題.我今天成功修改了内存的设置，多谢，不过，经过长时间的运行，最后还是出错，问题可能是因为步长什么的，不是很明白。所以我只好把我的输入inp文件分成四部份上传到这里，文件名是beam-deform.inp，四部份合在一起就是原文件，帮看一下，谢谢了无摩擦不行啊，因为钢筋在里面会乱跑，起不到承受拉力的作用了，假如不要摩擦，必须加其他的约束，以便可以让混凝土和钢筋固结在一起直到开裂。可以啊，你要把混凝土和鋼筋視為同一個part，四個體積，兩種材料，網格連續即可.abaqus 里面不是可以设钢筋层吗？我记得是：*rebar,里面有钢筋混凝土材料的，在standard mannual里面搜一搜试试：reinforcedconcrete ,或者rebar,不用像ansys里面似的，定钢筋定混凝土，再定interaction,吧。确实，我也看到这部份内容，不过不是很会用，还希望这么用过的大侠指点指点。我现在重新做过一个模型，比较简单，还是定义钢筋(弹性以及塑性)，混凝土(弹性，以及塑性里的混凝土，输入假设实验转化的数据，以及拉伸硬化等)，然后定义界面，interaction，界面定义是 tie，但是运行以后混凝土变形大致符合，不过钢筋还是和混凝土分离开了，很是苦恼，不知道什么原因。79028-Simple-beam.zip (0.0k)你的问题我试了一下，确实是计算无法收敛，初始时间增量无论改多细，第一个时间增量就是无法切入决定发散! 过程中警告讯息都是针对三对接触面的警告。 我看了一下模型，有些个人疑问，1. 钢筋与混凝土间的摩擦力真是您想观察之重点吗？若不是，是否可将两者视为结合一体，网格连续，材料相异 2.若真需计算摩擦，那么网格还需调整，钢筋网格密度不够，这样算出来的结果可能有误。至于，接近摩擦面之混凝土网格密度可密一点，但其余处可渐疏，如此应可减少模型的网格数量 3.这里还有些边界条件的调整会有trick，摩擦系数大小如何得知？此外，钢筋在初始计算时可能需要有用一些特殊拘束让它能与混凝土接触，否则，可能在接触计算的初始，在还没摩擦力产生时，钢筋会因无拘束而造成刚体运动，导致发散(因为您现在只于混凝土上加上拘束)这只是小弟的一些愚见，供您参考。因为最近有点杂务缠身，并无太多时间，不过，您这个问题颇为有趣，大家可再多多讨论，我杂务告一段落后，會再尝试解题How about this.? One part which has TWO cells, one is concrete, and the other is steel. Therefore, here has NO Tie. Follow ur geometry, material properties, load, and BC. But use mapped mesh and Nonlinear analysis. Detail. 能否看看你的inp文件呢，因为我的运行结果出来后，应力只是集中在受力部份附近，没有出现想你的结果附图里面那样，从中央向下端两个支撑处应力的由大到小传递，我从自己的inp文件也看不出原因。 80556-Simple-beam.zip (0.0k)李兄, 在此附上小弟前頁示範帖子之cae檔，共 beam.cae 與 beam.jnl兩個檔案，將其置於同一目錄下，直接執行 beam.cae 即可， cae file好處為保留 model geometry與 model feature，可以於CAE中任意修改.不會像看input file般霧裡看花，您可以從建模到後處理逐步瞧一下過程細節。不過，老實說，我這例子只是小兒科，沒啥大學問，您看看就好，不必太認真李兄, 你原始的 input file計算結果，元素明顯畸形變形，應該是hourglass發生 您應該將非線性分析選項打開才是，先前小弟示範亦是如此 修正如下: *Step, name=Press, nlgeom CAE可於step菜單中turn on nlgeom. 隨帖附上修正後所得計算結果 另外，seansheu 大侠，我按你的方法将非线性的开关打开后，结果就正常了，可是我不明白的是既然我已经将混凝土的塑性曲线数据按点输入了，肯定就不是线性的了，为什么还要另外加上非线性开关它才明白呢。几何非线性即通常所说的大变形分析，是指系统的微分方程等都建立在结构当前变形基础上的；不打开NLgeom即意味着系统的控制方程等都建立结构的初始构型状态。嗯，明白，就是说控制方程是采用前一步的结果作为起始量，而不是象线性可以采用初始值。關於willams_lee的問題，我重新回去檢視這個問題，當初因為初始定義用的是c3d8r單元，當線性錯誤時，直覺反應是認為沒開啟nlgeom來搭配材料非線性。但是當我回去確認應變時，發現它也沒大到那麼離譜.我切開模型，發現了致命的錯誤，鋼筋中心截面應力值不合理，變形分佈也不合理。顯然先前我作的假設有誤，所以計算錯誤.因為它網格密度並非充足，且有些單元品質不好，所以採用縮減積分的c3d8r單元計算會有問題，因此我改用了c3d8I，重新計算。此時，nlgeom開啟與關閉的結果都一樣。結果顯示，鋼筋中心截面應力值結果較為合理，變形分佈也合理，所以用c3d8r可能才是錯誤原因。非關幾何非線性.當然，我還要提高網格密度，去確認結果的收斂性。對willams_lee的誤導，深感抱歉。也感謝willams_lee的提問，讓我有機會發現錯誤。 在有限元分析中，应该尽量避免采用一次减缩积分单元（否则应该加密网格，有时也不一定管用）；而应该采用高次积分单元。 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 才能修改 地应力问题一、在桩土接触分析中，地应力步中，老是提示我this step cant 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来模拟水位的骤升、降问题。模拟结果不是很理想。不知道你有什么好的处理方法。请指教。多谢我有这样一个想法： 我们在一些情况下有了地应力，但是，我们对地表的初始变形不满意，那么我们先在第一步计算出地应力合初始变形，然后我们将地应力写入为一个文件，在这个文件里我们得到单元积分点上保持平衡的应力（或者其他，我还在想）然后我们通过USER SUBROUTINE 或者INPUT读入，那样我们不就得到了应力却没有如何初始变形吗！网上也有这样的帖子，但是过程不是很明了。但是我现在还不明白生成这个初始应力时应该生成哪一些结果，INPUT里的格式和生成的格式是否相同。这些都是问题三、【讨论】初始地应力的施加问题我们在分析一个边坡或者其他的岩土问题的时候，一定要考虑初始地应力问题的，比如边坡稳定性分析问题，一个自然边坡在我们开挖之前是出于自然平衡状态的，但是我们在模拟的时候，给它一个初始地应力的时候，它就不出于平衡状态了，和自然中的实际情况时不一样的。岩体中的应力为荷载产生的应力和初始地应力的和，因此初始地应力的大小直接决定我们的计算结果，如何正确的模拟这个问题到现在还没有很好的妥善的处理。由于初始地应力的存在，岩体的开挖将导致部分岩体卸载，但是在某些部位却是加载的，这个问题也导致了我们在模拟时的难度，也产生如何正确模拟这个过程的问题。然而通常我们都是这样处理的：沿开挖面作用着与初始地应力等价的“开挖释放荷载”的方式来模拟，把由此所得的位移即作为由于开挖产生的岩体位移，由此所得的应力场同初始应力场迭加即为开挖后的应力场。这种方法其实质就是等价原则的实施而已。我的问题的就是对上面阐述在我们的ABAQUS中的实现问题，我们在计算的时候总是要先 施加初始地应力的，当我们只是在这个力的作用下进行计算，我们发现我们的计算对象已经发生了严重的变相，这和我们的实际情况是不相同的，这个问题我们不管是在使用什么软件的时候都出现过，由此我本人请教各位高手对这个问题的理解是怎么样的？我们在实现模拟计算的时候如何处理这个问题？还有就是初始地应力的作用是如何体现的？对我指教一二，本人将不胜感激，我想这个问题也许是很多人想知道的吧。谢谢！呵呵，zhuxunguo所说问题确实存在，而且也没有很好的处理方法。看经验和问题的重要程度了，化些功夫还是可以做的比较满意的。通常的初始地应力场的施加（最简单的就是垂直向加以自重，水平向加以地质勘测给出的侧压力系数乘以自重）常常要引起网格的变形，而这个变形在事实上天然状态下是不存在的，因为我们以天然状态下已经发生变形的尺寸数据来做造型，现在让它再变形一次，显然引起我们的造型网格尺寸的变化，比较恼火呀。一次最好就是探讨探讨怎么能即加进初始应力场，又不产生变形，保证后续计算的正确性。 欢迎一起讨论。感谢楼上的小姐，我看了你告诉我的这个网址。受益匪浅，但是还是不是完全的理解。所谓初始地应力是时间和空间的函数。因此我们可以用场这个概念来反映。也就是初始地应力场。初始应力一般情况是由下面的几部分组成的：重力，构造力，动水压力，温度应力和结晶力等。在岩体或者土体开挖中破坏了原来的应力平衡，引起周围岩体的应力重新分布，形成一个新的应力场，称为二次应力场。它与工程断面是由关系的，并随开挖工作的进行再时间上，空间上不断变化。 上面是初始地应力和二次地应力的概念，从他们的概念中我们可以发现，初始地应力破坏后就会形成二次地应力，也就是我们为什么在施加了初始地应力后计算的边坡和隧道等支护体为什么破坏的原因。但是事实上再很多的情况下，我们开挖后支护体并没有破坏，这个原因是我需要知道的，在你告诉我的网址上大家并没有得出很好的结论，而且有很多的理解，尤其是对初始地应力和二次地应力的理解都不正确怎么可能知道这个原因呢，希望你能对我的提问给与我帮助，谢谢。也许我的理解有所偏颇，还请指点一二。谢谢感谢斑竹的鼓励，我的觉得和你们这些高手交流实在是受益匪浅。可能大家把的意思理解不是很对，我的意思是在简化的情况下，当然关键的因素都要考虑到不然我们计算的意义就失去了，斑竹你说呢？比如我们对材料使用DRUCKER-PRAGER 屈服准则，模型使用的摩尔库仑模型的时候，我们也一般都是使用这样的模型和屈服准则，当然由于材料的不同也有使用邓肯张模型的，剑桥模型以及剑桥修改模型的，都是必须要考虑初始地应力的问题，然后再施加我们要给它施加的载荷，以求解材料的受力情况。所以我提问的目的是为什么我们在施加初始地应力计算（先不施加二次载荷）的时候，材料发生了变形，这和我们实际情况中是不同的，因为一个边坡工程在我们开挖之前是达到平衡状态的，为什么我们计算的时候会出现这个现象呢？我不是很明白。所以提问，请教谢谢！我曾经计算过水工地下厂房，也涉及到原始地应力场的问题。当时做的时候好像是两种处理办法。当然首先要做回归计算，根据测点的应力数据，根据问题的性质假定机组力，用最小二乘法拟合得到计算模型的边界条件。一是将此边界条件施加到计算模型上计算一次，得到的应力场就作为初始地应力场，当然会出现变形，一般我们回忽略这个变形引起的造型误差，后续计算由于边界改变引起的应力场变化叠加到初始应力场上就是二次应力场，每次计算位移要减掉第一次计算得到的位移作为由于条件改变（如开挖、荷载变化）引起的位移。第二种方法是在模型上同时施加应力场和相应边界条件，由于内力与外力平衡，不会有位移产生，但是这样加好像麻烦一些。 好像就是这么做的，当然水工计算相对说要粗一些，不知说的对不对。不吝赐教！四、请教一个初始应力问题，多谢！我是新手，加了初始应力后，可是没有反应， 不知道什么地方出了问题 请各位大哥分析，谢谢了 * INTERACTION PROPERTIES * *Surface Interaction, name=IntProp-1 1., *Friction 0., * *INITIAL CONDITIONS,TYPE=STRESS bulk, 7800., 7800. soil, 1600., 1600. * BOUNDARY CONDITIONS *五、请问如何解决初始地应力问题请教诸位大侠，在基坑分析当中，如何解决初始地应力的问题？1，能否先计算自然土层（开挖前）的每个单元节点的竖向位移（重力作用下），然后再在以后的开挖回填荷载步中，将重力作用下的位移场减去2，或者能否先计算重力作用下每个单元节点的应力状态，再把这个应力场作为初始应力场加在模型上，然后再在这个被加上初始应力场的模型下进行开挖计算这些方法都是在ansys或是adina里面可以用的，可是我不知道如何在abaqus里面如何实现。请教各位，恳请指点一二，谢谢了abaqus的documentation里面的standard/manual中的analysis目录下面的coupled porous fluid flow中有关于地应力场如何处理的说明，abaqus中线建立一个初始平衡场，有应力没有位移，用*geostatic*initial conditions, type=stress, geostatic六、我们在处理土木工程的时候，一般都要遇到初始地应力的处理问题，基本都要使用如下的关键词： *INITIAL CONDITIONS TYPE=STRESS, GEOSTATIC 在手册中队其数据行的解释是如下： Data lines for TYPE=STRESS, GEOSTATIC: First line: 1. Element number or element set label. 2. First value of vertical component of (effective) stress. 3. Vertical coordinate corresponding to the above value. 4. Second value of vertical component of (effective) stress. 5. Vertical coordinate corresponding to the above value. 6. First coefficient of lateral stress. This coefficient defines the X-direction stress components. 7. Second coefficient of lateral stress. This coefficient defines theY -direction stress component in three-dimensional cases and the thickness or hoop direction component in plane or axisymmetric cases. If this value is omitted, it is assumed to be the same as the first lateral stress coefficient given in the previous field. Repeat this data line as often as necessary to define an initial geostatic stress state in various elements or element sets. 对其中我有几点认识不清的地方，首先我对(effective) stress的理解就是初始地应力，在数据行中2为Y方向的垂直地应力，3为其对应的坐标，但是对4种的解释地应力的垂直分量的第二个值如何理解？7侧应力的第二系数，如何理解，在ABAQUS中Y方向不就是垂直方向吗？如何把垂直方向按照是Z方向来理解，则7中的解释可以理解，也就是两个水平应力的计算中的系数，但是按照手册中是如何理解的呢？我现在没有弄明白，希望有朋友可以帮助我一下，谢谢。我看ABAQUS的帮助里，初始地应力好像是通过给定两点的高程，及相应的z方向应力值、水平方向应力系数来计算单元集的初始地应力，可以试着*initial conditions,type=stress（即不用geostatic参数）每个单元所有的应力分量值。Data lines for TYPE=PRESSURE STRESS: First line: 1. Node set or node number. 2. Equivalent pressure stress, .P 楼上的朋友你说的是不是这个啊？如果是的话，也许你的解释有点问题，因为这个参数的数据行只有两个，一个是节点集或节点号，另一个是等效压应力。您可以进一步解释一下你的理解吗？谢谢。七、【讨论】由初始条件地应力引起的接触问题如何处理?考虑重力场，需要在初始条件中加地应力,然后Step中施加Gravity与此相平衡，虽然地表产生的沉降不大，但不可避免总有一些沉降。我在地表放一基础，有Assembly将地基和基础方在一起。由于Grvaity产生的沉降，将两者分离，引起计算不能进行。如何解决这个问题？请高手帮助。谢谢“由于Grvaity产生的沉降，将两者分离，引起计算不能进行。 ” 1如果地基和基础在一个part里，不会产生分离。 2如果是两个part，应该设置接触条件。 3可以把基础单元先remove掉，geosattic后再激活，这样模拟施工顺序要好一点。 4把inp文件传上来看看嘛八、【原创】哈哈，终于建立地应力平衡了，用CAM-CLAY的平衡摸索了好几天，都快崩溃了。最后终于解决。之前也参考了我们这网上的一些文章。但总是解决不了。最后我发现原来主要问题在这里:只要你的模型的密度加速度值你的模型的重力作用方向的尺度应该与你的*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC里设置的应力相对应，想想：如果不对应模型就不平衡，然后就得去寻求平衡的街！我试了一下，改变侧压力系数并不会对初应变有太大的影响。这是我的模型，CONC是用来以后用的 239202-jk1.inp (0.0k)这是变形图 九、本人用ABAQUS软件对某地铁工程的施工过程进行了三维弹塑性数值模拟分析，现将有关内容显示在附件中。 340189-ABAQUS.doc (0.0k)看来大家还是很关注这个工程事例，现把其他几个文件给出，想必对大家关于ABAQUS的实际应用会有很多启发和帮助吧！希望斑竹能够给予鼓励，如果希望还会有更精彩的工程应用出现，我等待。 342333-other.rar (0.0k)你好，我刚运行了你的文件，老是出错，提示：ERROR: SPECIFIED STANDARD_MEMORY VALUE OF 256 MB IS TOO SMALL TO RUN THE ANALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUE IS 436 MB. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHER INFORMATION. 我的电脑配置是512内存，p4处理器，应该不会有问题吧？控制面板管理工具计算机管理，进入后，在计算机管理上点右键，编辑属性高级性能性能选项，更改虚拟内存，把最大和最小都改的大一些侧压系数是通过泊松比来体现的！请教各位：我运行了n次gremms大侠的inp文件，错误提示：ERROR: SPECIFIED STANDARD_MEMORY VALUE OF 256 MB IS TOO SMALL TO RUN THE ANALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUEIS 269 MB. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHERINFORMATION. 我已经改了虚拟内存，而且我的内存是512mb,不知道是什么原因 急盼各位大侠解答，急！你需要修改ABAQUS_V6.ENV中的内存设置,在WIN2000中ABAQUS支持2GB的物理内存，而在WINXP中支持3GB的物理内存，你的ABAQUS_V6.ENV中的默认只有256MB，你改为大于或等于512MB就可以了。goog luck！你在安装的.SITE下面找到该文件后打开,在文件的开始就有,一看便知现在把inp文件传上。 341818-BB8.inp (0.0k)十、【寻找】弹塑性模型中*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC*ELASTIC,TYPE=ISOTROPIC 1.0E7, 0.30 *DRUCKER PRAGER 36, 1.0, 36 *DRUCKER PRAGER HARDENING 5E4 *DENSITY 2.26E3 参数取值如上,为何初始地应力不平衡,而把 *DRUCKER PRAGER HARDENING 5E4 的参数调大,8E5,地应力变平衡,为什么?地应力平衡与DRUCKER PRAGER HARDENING有关吗? 弹性模型也演算过, *ELASTIC,TYPE=ISOTROPIC 1.0E7, 0.30平衡,位移为1E-6以上,为什么?GeoStatic是将模型恢复到初始的真实应力状态的一个平衡过程。你的模型在弹性阶段和塑性阶段对应的应该是同一个初始状态。你是想把达到塑性破坏的状态作为initial condition吗? 斑竹大人,你好,我的计算模型分别采用线弹性材料模型和弹塑性材料(DRUCKER),所加的初始应力相同 *INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC rock,0,0, -1.909E6,-86.2,0.5用弹性材料,初始位移很小(1E-6),地应力平衡,结果合理用弹塑性材料,初始位移较大(1E-1),地应力不平衡,结果不合理接着我就遇到了上面的问题,不知道我说明白了没有?我的意思是要在GEOstatic内保证初始平衡,这个阶段是不用考虑材料的塑性啊,刚性啊,流动性啊的。因为你使用的是TYPE=STRESS对呀,GEOstatic内保证初始平衡,这个阶段是不用考虑材料的塑性啊,刚性啊,流动性啊的。因为使用的是TYPE=STRESS但是后来的计算需要先前定义的模型材料参数材料属性分弹性和弹塑性同一模型,可以用不同材料参数,而其他部分一样当用弹性材料时,地应力很容易平衡,但是我改成弹塑性材料,地应力就不平衡了,是这个意思,不是同一过程的不同阶段而是两个分别的过程在GEOstatic内保证初始平衡我也并没有考虑材料的属性,而是前面定义的材料属性我当然知道你已经是不同材料的不同计算case。但是平衡是使用的材料参数不就是一个density或指定应力吗?谢谢斑竹！问题已解决，在用GEOSTATIC平衡时，我认真阅读了的例子，发现在地应力平衡与DRUCKER PRAGER HARDENING参数有一定的观系，当DRUCKER PRAGER HARDENING较大时，与弹性材料平衡一致，当DRUCKER PRAGER HARDENING较小时 ，需要添加新的边界条件，强制使其初始位移为十一、求助:initial conditions的一个问题不好意思，我看了以前的关于这方面的讨论，但是，还是没有解决自己的问题，实在没有办法，拿出来占大家的时间，谢谢！我做的类似桩基，是土体中加入桶形的桩基础。现在对土进行了iinitial conditions，现在土的初始位移基本为0了，然后我对桶同样用initial conditions，想使桶的初始位移也接近于0，但是，发现这个initial conditions对于桶不起作用，换了很多的参数，题目也不报错，好多天了，没有解决，希望得到大家的帮助，谢谢，谢谢！ 315031-Job-1.inp (0.0k)*Material, name=SOIL *Density 1600., *Drucker Prager 54.1,1.,0. *Drucker Prager Hardening 50000.,0. *Elastic 500000., 0.2 这里的*Drucker Prager 54.1,1.,0. 中的0好象不太对,一般设为你的土的内摩擦角是54.1差不多大fangxiang，不好意思，由于工作太忙，现在在才回答。你给我的那个文件是在我上次的那个基础上多加了一个step吧。你加的是在桶的最左边一点加了一个1000的水平荷载。计算结果没有问题啊，以下是在放大2000倍的位移图。 还有水平位移等值线图。 现在这个问题还没有解决，希望各位大牛给与关注阿我又做了个简单的模型桶，但是通过initial conditions,type=stress,geostatic及initial conditions,type=stress,input=xxx.inp，同时计算，但是发现单桶的竖向初始位移还是恢复不到零，发inp文件，请各位看看，先谢了是不是initial conditions 的位移恢复，和模量有很大的关系，不是以前简单的考虑，希望各位帮助，指点，谢谢，多谢 Job1.inp (52.1k)而且你没有定义初始空隙比和初始静水压力这是一个简化的模型，有些内容不是计算所主要关心的，所以就省略了，现在的主要症结是initial conditions为什么对土有效，而对桶就失效了，多谢，再次感谢帮忙！没有啊，你的桶的initial也是加进去了的。有一句你写错了 把_G14,-764400,0, 40,50,0.5改成 _G14,-764400,40, 0,50,0.5 就可以了，你0,40两个数字的顺序做反了。 下面是平衡后的竖向位移。 下图是在加重力前的初始竖向应力场，由于桶的质量比土大了很多倍，所以初始竖向应力桶的也比较集中。 以下是第一步加入重力后，应力场就变均匀了。 你原来的inp文件的结果的刚好相反，所以会导致在桶的自重作用下产生相当大的位移而不能平衡。 对了，还有一点就是平衡地应力要用*geostatic来做，而不能用*static来做。用static不能让初始位移为0。十二、有关初始应力的问题一根悬臂梁，在自重作用下弯曲，求出在自重作用下的单元应力，然后单独把求得的单元应力加于梁上，此时梁的扰曲线应与自重作用下的扰曲线关于水平面对称，如果把自重应力和自重都加上，作为梁的初始状态，则此时梁应处于自重平衡状态，没有变形。 ansys就是利用这种原理施加初始应力场，作为分析的起始步，但在abaqus中，初应力不能从前面的结果文件中读入，在土工问题中，如果用*initial conditions,type=stress,geostatic来定义初应力，按前面所述，如果只加初应力（不加自重），在二维情况下这种与自重作用下的反向变形很明显，此时的地压平衡也很容易实现；在三维情形，用这种方法加上的初始应力就没有这种与自重作用下的反向变形，因此三维的地压平衡就很难实现（重力场作用下产生变形，而初应力作用下土体没有变形）。因课题涉及三维问题，地压平衡这一步总是实现不了，没法往下做，甚为郁闷，各位有无做过这方面的问题，请指教一二，或是有无这方面的例题。BOW！能否解释解释, 只加初应力（不加自重），产生的反向变形是什么意思?难道说, 只加初应力（不加自重），地表要向上移动?是的，你可以试试，这种情况下地表只产生向上的位移。三维地压平衡问题已顺利解决（在此感谢版主cafe的关心！），问题出在对abaqus 坐标轴的理解上，在二维问题中，坐标轴的竖向是2方向（y轴），初始有效应力就是加在2方向上，当然土体的自重也就加在2方向上了: *dload elset_name,grav,10,0,-1,0而在三维问题中，坐标的竖轴就当为3方向（z轴），而不是我们在cae界面中所见的2方向，所以初始有效应力是加在3方向上（而在二维情形，有效应力是加在2方向），因而，在加土体自重时应当加在3方向上： *dload Elset_name,grav,10,0,0,-1十三、【讨论】再谈初始应力的施加问题岩土问题的研究都要涉及到一个初始应力的问题，其施加问题一般都要遇到。下面我就此问题谈谈个人的观点：初始应力包括垂直和水平两个方向的力，所以我们在施加的时候要注意两个方向都要有作用的力。一般水平和垂直方向的比是K=1-sin .但是我现在有个问题，在abaqus 中施加初始应力的时候我认为是在第一步中施加，可是在create load 对话框中我并没有看到所需要的力，请高手指教。附件中。 176988-a1.rar (0.0k)如果考虑构造应力，水平应力和垂直应力的比一般为2：1对于考虑构造应力的话的，水平应力和垂直应力之间的比值并不是2：1。不知你是哪里的结论，请告之，谢谢。因为据我所知，地应力的主要来源是地质构造运动和岩体自重。他们之间的比值是变化的，不是固定的，当然详细的解释有兴趣的朋友可以参阅有关的书籍，比如高等岩石力学，矿山岩体力学，那里有比较详细的解释。希望和楼上的小姐共勉，和我联系，我的邮箱是我们一起探讨这个问题。我的说法可能让你误会了，一般情况下，若有构造应力的情况（一般为水平方向），根据各国地应力的测试资料，通常来说，在一定的深度内水平应力大于自重应力，他们比值不是固定的，我们在实际计算时，若考虑构造应力，而且没有实测的地应力资料，可以取2：1十四、【讨论】紧急求救：geostatic 分析中如何加初始应力和侧压系数？用AQAQUS CAE时,在step模块选用geostatic 分析时， 如何体现侧压系数？ 初始应力如何体现？是不是必须写keywords, 请给一简单例子。 小妹先谢了cae里不能实现初始应力的施加，可在inp文件里添加，祥见document里abaqus keywords reference manual 里*initial conditions,type=stress,geostatic里的说明。十五、讨论】再谈ABAQUS中开挖工程的初始应力场问题！对于涉及开挖问题来说处理好初始应力问题是至关重要的，目前，常常两种方法，其一就是直接把重力场（也可以有其它力存在）始终作用在模型上，进行开挖模拟，这样求得的应力场为实际应力场，而求得的位移场需要减去初始位移场，其二就是想办法读进初始应力场到模型中（没有初始位移），然后进行开挖，这样所求得的位移场和应力场均为实际的位移和应力场。本人认为第一种方法最好，因为我们的初始应力场既可能为弹性也可能局部进行塑性，而第二中方法似乎只能处理弹性初始场，如果真是这样，不知在ABAQUS中能够实现第一种方法，也就是能否方便实现各开挖步中位移场的相减运算，请高手指点！可以得，在inp里修改或者cae里也行十六、在inp文件中加地基初始应力（not geostatic），怎么读不进去呢？采用C3D10MH单元，命令如下：*Initial conditions, type=geostatic1,1,s11,s22,s33,s12,s13,s231,2,s11,s22,s33,s12,s13,s231,3,s11,s22,s33,s12,s13,s231,4,s11,s22,s33,s12,s13,s232,1.2,2.2,3,.2,4,.该单元四个积分点，不知道是不是应该这么加。abaqus/pre.exe错误提示信息如下：*ERROR: in keyword *INITIALCONDITIONS, file tet10coashMH.inp, line 9248: Unknown assembly id 1.*NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO INTERPRET SOME DATA. SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION*ERROR: AN INITIAL CONDITION HAS BEEN SPECIFIED ON ELEMENT 0 BUT THIS ELEMENT HAS NOT BEEN DEFINEDCARD IMAGE: , 1, -5.4809E+04, 27.20, -2.6842E+04, 5.717, 0.9051, 1.552我试过只按一个积分点加，错误信息同上。那位帮看看可能是怎么回事？多谢了！还有一个问题，每个单元是四个积分点，但是abaqus只能取一个，作为初始应力平均加到整个单元上。请问能不能向输出单元应力那样把不同的初始应力分别加到每个积分点上，或是四个节点上也行啊（最起码保证应力在单元内线性）！你是不是先用abaqus的cae建的模型然后自动剖网格的吧，然后再利用文本编辑器编辑inp文件，是不是，这种情况他都会提示unknown id，他把节点号当成了几何信息里的set，我也不知道更好解决的方法是什么，这个建议你试试吧，你可以把关于inp文件的instance,assembly的信息都删掉，把inp文件里的模型，节点，单元组成形式搞成orphan单元的形势，就像那些example problem manual里的文件一样就可以了这条经验很宝贵，ABAQUS DOCUMENTATION里面恐怕没有。虽然搞不清原因，但确实我改了ABAQUS自动生成的INSTANCE,ASSEMBLY之后就好用了。这是条值得加分的帖子！真是很感谢inchesnail！还有一个问题，每个单元是四个