平稳，是tid，还是did？
点击联系发帖人
时间：2018-07-02 15:09
数据库实验用表_百度文库
您的浏览器Javascript被禁用，需开启后体验完整功能，
享专业文档下载特权
&赠共享文档下载特权
&10W篇文档免费专享
&每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
数据库实验用表
&&数据库数据
阅读已结束，下载本文需要
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
加入VIP
还剩1页未读，
定制HR最喜欢的简历
你可能喜欢(window.slotbydup=window.slotbydup || []).push({
id: '2014386',
container: s,
size: '234,60',
display: 'inlay-fix'
&&|&&1次下载&&|&&总176页&&|
您的计算机尚未安装Flash，点击安装&
阅读已结束，如需下载到电脑，请使用积分（）
下载：60积分
相关分类推荐
0人评价3页
0人评价3页
0人评价3页
0人评价9页
0人评价9页
所需积分：（友情提示：大部分文档均可免费预览！下载之前请务必先预览阅读，以免误下载造成积分浪费！）
（多个标签用逗号分隔）
文不对题，内容与标题介绍不符
广告内容或内容过于简单
文档乱码或无法正常显示
文档内容侵权
已存在相同文档
不属于经济管理类文档
源文档损坏或加密
若此文档涉嫌侵害了您的权利，请参照说明。
我要评价：
价格：60积分VIP价：Tiddo de Ruiter
SATURNUSSTRAAT 91
2516 AG DEN HAAG
W www.tiddoderuiter.nl[转载]ABAQUS常见问题汇总&-&2.0版（7-9）
good．弹塑性分析
塑性问题不收敛时，msg文件中的常见现象是
1）&& 出现很多EQUILIBRIUM
ITERATION，且TIME INCREMENT 不断减小；
始终出现&& ***WARNING: THE STRAIN
INCREMENT HAS EXCEEDED FIFTY TIMES THE STRAIN TO
CAUSE& FIRST YIELD AT 1 POINTS
&***WARNING: THE STRAIN INCREMENT IS SO LARGE
THAT THE PROGRAM WILL NOT ATTEMPT& THE PLASTICITY
CALCULATION AT 1 POINTS
3）&& 在msg文件的结尾显示
***NOTE: THE SOLUTION APPEARS TO BE DIVERGING. CONVERGENCE IS
JUDGED UNLIKELY.
***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT
解决方法见下面回帖实例。
问：aba_aba老大，你说仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent+j4d&Cd-
用C3D15单元，是15-node quadratic triangular prism,
这还是一种二阶单元，不能用于塑性材料和接触问题。|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent(X&@(p/G$?:J0w
这个是什么原因？www.simwe.com)Te)E#|*@M
-- Getting Started with ABAQUS 12.4.4 Element selection:
When selecting elements for contact analyses in ABAQUS/Standard,
it is better, in general, to use first-order elements for those
parts of a model that will form a slave surface. Second-order
elements can sometimes cause problems in contact simulations
because of the way these elements calculate consistent nodal loads
for a constant pressure.
-- Getting Started with
ABAQUS&&&&
10.3 Selecting elements for elastic-plastic problems:
The fully integrated, second-order, solid elements available in
ABAQUS/Standard are very susceptible to volumetric locking when
modeling incompressible material behavior and, therefore, should
not be used in elastic-plastic simulations.
The second-order reduced-integration elements in ABAQUS/Standard
should be used with caution if the strains exceed 20&40% because at
this magnitude they can suffer from volumetric locking.
问：在我的模拟一个桶土接触问题中，我想加一个和坐标轴成一定角度的集中力，但是直接加集中力会导致不收敛，因此我改用加位移荷载的方法，求出反力，再反推合力，以及合力与坐标轴的夹角，请问这个合力方向是不是与所加的合位移方向相同，本人认为是相同，不知道大家有什其他的看法，请教一个？
有些情况下一样，有些情况下不一样，你试试就知道一样不一样了。施加位移荷载时的反力就等于力载荷。
“直接集中力会导致不收敛”：
一是加载之前要有一个只施加很小载荷的分析步，使各个接触关系都建立起来。
二是如果受集中力的材料是塑性，会因为局部变形过大而不收敛，应该将力施加在一个小区域上，让多个节点共同受力。
当设定stress 和 strain时总不能计算出结果
只有elastic的时候就没问题
当设定时注意什么
要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ. 见附件。
载荷不能太大，改为-1500。
问：二维模型为鋼板接頭實驗，因为是对称的，所以只取其中一半。上、下圓為螺栓（設定為rigid body），Step
1：先讓上；下螺栓接觸鋼板，並施加一微小壓力避免接觸時發生表面震顫。Step
2：下方螺栓固定不動，上方螺栓向上位移直到位移結束。
~7Z-S1_V dwww.simwe.com
qv‑V(y*k8E8oF仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM我的問題是：|Simwe.com|仿真|设计|有限元|虚拟仪器3L%Tfn.L[1]E:T/S
(1)我下方圓孔接觸的地方附近再割密一點，就會發生”too many attemps made for this
increment”，是不是指元素扭曲過大，收斂困難???（因為我把網格切大一點又就可以跑，我有試過aba-aba在別人的回文中的方法，將菜单other
/ Gerneral solution controls / Edit,把Is由默认的12改为30，可結果還是一樣。
t7p仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
y9Ae)bzPQ8B仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM(2)另外，因為我感興趣的部分是在位移過程中（step
2）我想知道鋼板所受的外力，因此我在上方螺栓作reaction force的輸出（REF-1），但是我發現分析出reaction
force來的結果其前面那一部份還蠻符合，但到了後面塑性的部分，力就開始發生不規則的震盪，不知道是哪裡有錯，仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1?
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMe&h5E
2裡，為什麼在一開始位移時卻沒有受力，是因為mesh不夠密嗎（平滑）？，導致一開始並沒有完全接觸到，所以沒受力。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM_I?
|q|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent針對鋼板下方螺栓圓孔的上部頂點位移與上方螺栓的反力，畫出load-displacement關係圖，與真實實驗結果作比較，供大家參考。希望知道問題的人，能不吝嗇指教我。SimWe仿真论坛l`ze&f‑U
www.simwe.comFJ3spY
謝謝了！！
附件: steel-plate.rar ( 20:35, 34.91 K)
msg文件中有
&***WARNING: SOLVER PROBLEM. NUMERICAL
SINGULARITY WHEN PROCESSING NODE STEEL PLATE-1.52 D.O.F. 2 RATIO =
意思是STEEL PLATE-1在2方向上缺乏约束。比如应固定STEEL PLATE-1的底边。
螺栓位移为17.7125，相当于圆孔的半径，圆孔处变形相当大，大变形处的网格需要相当细化，以避免大变形造成畸形单元。细化网格后可能load-displacement關係圖就不震荡了。
螺栓单元也太少，导致主面上有尖角，不容易收敛。每90度圆弧上应至少有10个单元，最好用analytical rigid.
要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ. 见
http://www.simwe.com/forum/viewthread.php?tid=717593
是塑性材料曲线的问题：
- 曲线开始的部分应力上下变化，使ABAQUS很难找到收敛的位置, 应该为一条光滑的曲线；
- 给出的塑性数据点太多，好像也会造成收敛困难，建议改为
&343.964,&&&&&&&
&391.794, 0.0178509
&409.007, 0.0247034
&423.079, 0.0308124
&437.593, 0.0373512
&450.328, 0.0411674
&470.377, 0.0512111
&490.794, 0.0612664
&507.881, 0.0714607
&522.588, 0.0812225
&538.149, 0.0909967
&550.917,& 0.100796
&595.463,& 0.138654
&628.621,& 0.174038
& 654.65,& 0.206945
&675.526,& 0.237475
&695.278,& 0.270728
&706.411,& 0.304198
&700.732,& 0.320962
&643.526,& 0.355134
& 571.81,& 0.374608
& 492.35,&&
470.,&&&&&&&
这样就收敛了（没有最后一行470, 1 ，也可以收敛 ）。在结果中可以看到：
- 大变形处单元扭曲非常厉害，应该在此区域 partition, 使用更加细化的网格。如果有ABAQUS
6.6，可以做自适应网格重划。其他应变小的区域，网格可以粗一些。
- PEEQ最大值为1.719，高于*Plastic最后一行中的塑性应变。对于这样大的变形，难以得到精确的塑性材料曲线：
拉伸试验到一定程度就断了，而且颈缩后的实验结果就开始不精确；压缩试验可以压到很大的应变，但这时摩擦力的影响变大，试样会出现鼓形，结果也不精确。好像没什么太好的办法。
-&&一般的大变形可以用ABAQUS/Standard,
如果想模拟材料破坏就要用ABAQUS/Explicit，见ABAQUS Analysis User’s Manual
第11.6节Progressive damage and failure。
- 弹簧刚度1000太大了。钢板底边的最大y向位移是12，所以弹簧分担的载荷是12*,
而支反力最大值为52070，与弹簧分担的载荷在同一个数量级上。弹簧的作用只是帮助收敛，不应该影响各个部件的受力和变形。可以试试减小弹簧刚度，例如刚度为1，不收敛的话再增大一些。
问：问题 自适应网格重划
之前我用二维模拟钢板接头实验，但发现模拟结果不好，所以改用3D模擬。上、下圆为螺栓（设定为rigid body），Step
1：先让上；下螺栓接触钢板，并施加一微小压力避免接触时发生表面震颤。Step
2：下方螺栓固定不动，上方螺栓向上位移直到位移结束。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMNR
B i&p5oza"H#A
5W~eF2i5a,h-k|Simwe.com|仿真|设计|有限元|虚拟仪器问题是：|Simwe.com|仿真|设计|有限元|虚拟仪器4T.h!o.?Z)]+H"VM:_
(1)还是老问题，我在上方螺栓作reaction force的输出（REF-1），但是我发现分析出reaction
force来的结果其前面那一部份还蛮符合，但到了后面塑性的部分，力发生不规则的震荡SimWe仿真论坛S
Q‑K/C9M%r&M0bT
V4yz&I:b/P"`r(2)我的是ABAQUS6.5，因为下方变形很大，为了避免元素扭曲过大，所以我用自适应网格重化下方大变形螺栓处，但却一开始就发生Abored，不知为什么？
t&J]仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
2D模型的结果怎么不好？你的应力应变在厚度方向几乎没有变化，所以2D和3D的结果应该是一样的，而3D的计算时间长得多。
ABAQUS 6.6才能真正实现自适应网格重划，ABAQUS
6.5的自适应网格主要用于Explicit，以及Standard分析中的表面磨损过程模拟，在ABAQUS
6.5一般的Standard分析中尽管也可以设定自适应网格，但不会起到明显的作用，只是大大增加计算时间。
下列警告都是非线性问题迭代过程中的正常现象，是ABAQUS正在尝试找到正确的解：
&***WARNING: THE SYSTEM MATRIX HAS 8 NEGATIVE
EIGENVALUES.
&***WARNING: THE STRAIN INCREMENT HAS EXCEEDED
FIFTY TIMES THE STRAIN TO CAUSE
&&&&&&&&&&&&
FIRST YIELD AT 34 POINTS
&***WARNING: EXCESSIVE DISTORTION AT A TOTAL OF
2 INTEGRATION POINTS IN SOLID
&&&&&&&&&&&&
(CONTINUUM) ELEMENTS
&***NOTE: ELEMENTS ARE DISTORTING EXCESSIVELY.
CONVERGENCE IS JUDGED UNLIKELY
（以当前的increment不能收敛，自动减小increment，重新迭代）.
在后处理时可以看到大变形而严重扭曲的单元，应在这些地方进行网格细化。
在你的实验中，钢板能够发生这么大的应变而不拉断吗？如果出现了材料破坏，用ABAQUS/standard得到的这种大变形结果就会和实际相差比较大。
螺栓可以用analytical rigid, 以缩短计算时间。
你说的“力发生不规则的震荡”，是说RF出现了下降/上升的变化吗？这好像是正常的，所谓
“局部卸载”。你的塑性材料应力应变曲线有下降段（从数据点
706.445,&& 0.30186
开始），就是说应变达到一定程度时，应力会下降，自然RF也就下降了。
你的试验和FEM已经吻合得不错了，不可能完全一致的。
ABAQUS/Standard不能模拟材料破坏，要用explicit, 见 ABAQUS Analysis User’s
Manual 第11.6节Progressive damage and failure。
ABAQUS 6.6 的自适应网格很不错，可以用于standard分析.
RF震荡，可以用排除法来找原因：
- 怀疑是摩擦造成的： 去掉摩擦，看是否还震荡；
- 怀疑是网格造成的： 螺栓用analytical
rigid，钢材大变形部位用非常非常细的网格，看是否还震荡（你现在大变形部位的单元发生严重扭曲，肯定是太粗了，得到的结果肯定很差）；
-&&&&&&&&&
怀疑是材料塑性曲线的下降段造成的：将材料曲线的下降段改为单调上升，看是否还震荡。
在msg文件中看到反复出现
&&&&&&&&&&&&&&
SEVERE DISCONTINUITY ITERATION&&
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
5 OPENINGS.
&&&&&&&&&&&&&&
SEVERE DISCONTINUITY ITERATION&&
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
0 OPENINGS.
即ABAQUS无法找到正确的接触状态。这是因为x,y方向没有边界条件，只有几个软弹簧，板还是会出现难以确定的位移。
应该尽量让每个部件在每个自由度上都有足够的边界条件：
- x方向： 你的模型是对称的，可以只取1/2，约束x方向位移。
- y方向： 板的顶部(大螺栓的一端）变形很小，可以约束此处的y方向位移，让大螺栓不动，小螺栓向下移动。
increment可以设为
0.001, 1., 1e-8, 1.
请教各位大侠这是怎么回事：
Solver problem. Zero pivot when processing D.O.F. 1 of 1 nodes. The
nodes have been identified in node set
WarnNodeSolvProbZeroPiv_1_1_83_1_2.
OVERCONSTRAINT CHECKS: The model is overconstrained in this
increment due to the interactions among nodes in node set
WarnNodeOverconZeroPivotStep1Inc83Iter1Cut1.
我没有看到OVERCONSTRAINT 信息。如果只是dat文件中有OVERCONSTRAINT
信息，msg文件中没有，就没问题。
刚塑性材料不容易收敛，见Getting Started with ABAQUS 10.4.4 Adding hardening to
the material model。另外要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ.
*Material, name=Material-1
10.0, 100.0
*Material, name=Material-2
10.0, 100.0
把缩放系数设为1，可以看到单元扭曲太大，所以无法收敛。要在大变形部位细化网格，变形小的部位的网格可以很粗。
&&&&&&&&&&&&&&&
问：桩土接触（群桩） 塑性为什么计算不出结果？（弹性可以）
本人一直在研究群桩的计算&&今天计算了一个七根桩的模型（三维）。www.simwe.com#Al%o
模型具体为：
&i8q5z3f[1]L&d3{ @n2}
q!{仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&&&桩：由圆桩（直径2.8m）等效为八面体桩&&桩长100m
全部入土& &梅花形排列www.simwe.com$_
a@**x3Y0r'C
&&&土：共七层土体&&有五层和桩接触&&土体40*40*200m
SimWe仿真论坛
&&&接触：每根桩和五层土体接触（共7*5＝35个接触对
node based surface），桩底采用tie约束&
&共27000左右个单元（c3d8r）
BA@ q-QN)d?F5b^&
&&&我分别用D-P模型和M-C模型都无法算出结果，后来我改为弹性可以算了而且结果还可以
并且平衡了地应力 I|
G|(x7& & &
&& & 讲了这么多
我就是想知道
为什么材料为塑性就算不出结果（经常是单元变形过大&&材料点发生屈服）
是我的塑性参数取的不合理吗？还是别的什么原因？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/_K7p
&& &一般土体yield
stress取多少？第一个step 平衡地应力可以用塑性材料吗？如果用弹性材料 后面计算怎么该为塑性材料？请各位给点建议
我是把yield stress
该大了&&一直不清楚这个值到底该怎么取&&m－c取为c的话
，土体也就几十Kpa&&是不是太小了？ 有没有人知道具体应该怎么取
？ 我200m的土&&如果yield
stress为30Kpa那么在自重作用下不是已经屈服了吗？
问：只在四节点施加一竖向位移，而水平方向怎么会出现应变，且同竖向的应变在同一量级？把弹性单元编为abaqus的uel子程序计算结果也不一致
边界条件对吗？三个节点的位移都固定死了，无法保持体积不变。水平方向应变是因为有泊松比。参见
http://www.simwe.com/forum/post/view?bid=31&id=717759&tpg=1&ppg=1&sty=1&age=0#717759
不同软件要比较单元积分点上的应力，节点应力是插值得到的，不准确。
先做普通弹性材料，正常了再做塑性。
另外，dat文件中有warning:
&***WARNING: A DILATION ANGLE OF 20.000 MAY
RESULT IN DECREASING PLASTIC
WORK&&&&&&&&&&&&&
OR UNSTABLE MATERIAL BEHAVIOR AT HIGH CONFINING STRESS STATES.
SET&&&&&&&&&&&&&
THE DILATION ANGLE LESS THAN 17.831 TO ENSURE STABLE
MATERIAL&&&&&&&&&&&&&
BEHAVIOR UNDER ALL LOADING CONDITIONS
问：在塑性加工中，不同的材料用如下变形抗力模型表示。请问：这个公式指的是整个弹塑性过程还是忽略弹性的？如果包括弹性，那么弹性模量又与该公式有何关系？屈服极限应力在该式中有体现吗？该公式反应的是真实应力应变关系吗？它的适用范围是什么？
不胜感谢！！！！！！
wy_anli答：塑性加工中的具体材料模型有成千上万个，大类可以分为理想刚塑性（忽略弹性、忽略硬化，主要应用于弹性变形相对于塑性变形较小，可以忽略的情况下，如：锻造等）、理想弹塑性（忽略硬化）、弹塑性（在材料成形时其弹性变形不可忽略，例如板料成形）、等等。
你的这个模型属于弹塑性模型，没有忽略弹性。
屈服极限是材料发生塑性变形时的临界应力，是材料的本身性质，在本公式中不能体现。
纵向应力与纵向应变的比例常数就是材料的弹性模量，也叫杨氏模量。也不能在此公式中反映出来。
任何材料模型都不可能完全反映真实的应力应变关系，当然次模型也不例外，只是最大近似。
一般只有根据材料提出材料模型，其适应范围也只能适用于特定材料以及特定场合。
wpengw123 wrote:
非常感谢。还是有些不解，就直接问了，不吝赐教！
那么该公式是否就是整个弹塑性应力应变关系的简化式呢？还是只描述塑性变形阶段？
是塑性应力应变关系的简化式，但还缺少抗拉强度或极限应变（当材料达到这个临界状态时就会发生断裂，继续算下去，就没有了意义）。弹性阶段在塑性阶段前考虑。
wpengw123 wrote:
如果包含弹性，那么杨氏模量是不是和公式有冲突呢？
杨式模量不会和公式有冲突，杨式模量是材料本身的参数，是材料在弹性阶段时应力应变线性关系的系数。Abqus中定义材料时，弹性和塑性是分开来的。弹性时只有杨式模量和波松比两个参数。
wpengw123 wrote:
而且它似乎常用于金属板料成形，优点是什么呢？
板料成形要考虑弹性变形，因为板料成形是大变形小应变情况，弹性应变不能忽略不计，并且在塑性成型结束后弹性变形要完全恢复（也就是回弹）。相比较来说，例如锻造等成形，弹性应变相对于塑性应变非常小，可以忽略不计，也就是要采用刚塑性材料模型。
wpengw123 wrote:
式中常数和变量的具体意义是什么呢？
参见塑性成形原理。
wpengw123 wrote:
怎样和单向拉深的实验数据联系，又怎样在Abaqus中使用呢？
单向拉深的实验数据通常是名义应力和名义应变。要通过名义值和真实值之间的关系式将其转化为真实值也就是等效应力和等效应变。Abaqus中的材料采用真实应力和真实应变。
wpengw123 wrote:
看到很多文章在模拟板料塑性成形时用到此模型，但是在给出杨氏模量及强度系数和硬化指数的取值时却都没有给出屈服极限，屈服极限不需要吗？
屈服极限也就是材料在屈服点上的应力应变状况。材料模型中应变为0时的应力就是屈服极限，大多数金属的屈服应力为材料弹性模量的0.05%-0.1%。
qyuan_13: &关于对球墨铸铁的强度评价标准
通常情况下，铸铁的脆性比较强，延伸率很低，一般以延伸率为0.2%作为屈服极限。对于脆性材料作强度分析的时候，采用主应力来评价；但对于球墨铸铁，有比较高的延伸率，比如QT450-5、QT450-10，其延伸率分别为5%和10%，那么，对这样的材料作强度分析的时候，是不是应该采用Mises应力来评价？按说10%的延伸率比较高了，合金钢20Cr的延伸率也为10%，40Cr的仅为9%（刘鸿文的材料力学）。那么，在脆性材料和塑性材料之间的界限到底是多少呢？不能仅仅认为只要是铸铁，就得用主应力来评价吧！
GZSI2x2yiSimWe仿真论坛
问：在ABAQUS中,对于理想弹塑性材料,在计算时,ABAQUS中,材料屈服有很多准则.怎么在ABAQUS中定义我所要采用的Mises屈服准则?
在帮助里搜索一下就可以找到： ABAQUS Analysis User's Manual 11.2.1 Classical
metal plasticity：
The Mises yield surface is used to define isotropic yielding. It is
defined by giving the value of the uniaxial yield stress as a
function of uniaxial equivalent plastic strain, temperature, and/or
field variables. In ABAQUS/Standard the yield stress can
alternatively be defined in user subroutine UHARD.
Input File Usage: *PLASTIC
ABAQUS/CAE Usage: Property module: material editor: Mechanical /
Plasticity / Plastic
理想弹塑性容易有收敛问题，最好加入硬化特性，见Getting Started with
ABAQUS&& 10.4.4 Adding hardening
to the material model。
你的弹性模量是不是太小了，以钢为例，如果单位是N, mm，MPa，则弹性模量是210000, 屈服应力是420。
判断材料是否屈服
问：DOCUMENT上说PEEQ是等效塑性应变（Equivalent plastic
strain),PEMAG是塑性应变值（Plastic strain
magintude),感觉前者主要是判断材料是否屈服，而后者呢？？
在比例加载时（即加载过程中主应力方向和比值不变），大多数材料的PEMAG和PEEQ相等。这两个量的区别在于，PEMAG描述的是变形过程中某一时刻的塑性应变，与加载历史无关，而PEEQ是整个变形过程中塑性应变的累积结果。例如，单向拉伸一个圆柱体，使其发生塑性应变，再通过单向压缩使其恢复初始长度，则最终的PEMAG为0，而PEEQ是拉伸和压缩过程中塑性应变之和。
问：AC YIELD是材料进入塑性的标志，1是发生了塑性变形，0是没有，这样理解对么？
PE是塑性应变，可为什么AC YIELD为零的区域（没有发生塑性变形）也会有PE值呢？？
积分点（而不是节点）上的 AC YIELD 和 PEEQ (而不是PE)是完全吻合的： PEEQ
& 0 则AC YIELD = 1, 发生了塑性变形; PEEQ = 0 则AC YIELD = 0,
没有发生塑性变形。
云纹图是经过插值平均的，看上去AC YIELD 和 PEEQ
会显得不一致。如果逐个察看积分点上的结果，就可以看出它们是吻合的。
PEEQ 和 PE 的关系： PEEQ 类似于mises应力，PE类似于各个应力分量。
问：比如计算地基的极限承载力，m-c理想弹塑性，用计算结果不收敛作为破坏评判准则。
现在我做的时候，荷载是手动调整，每次调整完再提交计算，感觉比较麻烦
能否实现这样的功能：给出初始荷载以及增量步长，要是能自动调整就更好了。每一级计算完成后保存结果，并增加荷载继续计算，直至不收敛
ansys好像有循环语句可以实现，不知道在abaqus中如何实现。谢谢了~
不需要特别的设置，直接创建一个step, 设置initial increment（例如0.1）,
在这个step里施加一个载荷（比如1000N）, 于是step time = 0.1时，载荷就是100N; step time =
0.2时，载荷就是200N；超过极限载荷时，计算不收敛，ABAQUS会自动减小increment。
aba-aba,你有没有碰到过这种问题,在理想塑性状态(切线模量设为零)下,算出来的结果最大应力超过屈服应力?而并不是产生不收敛?
席勒同志好，关于这个问题，我的理解是这个样子地：
如果一个单元所受的载荷是单调增加的，则其积分点上的Mises应力和等效塑性应变PEEQ与塑性材料曲线中的真实应力和塑性应变相吻合。注意是积分点，而不是外插后得到的节点应力。
另外多说两句应变方面的问题，PEEQ描述的是整个变形过程中塑性应变的累积，，无论是拉应变还是压应变，都取绝对值，然后加在一起。如果分析过程中出现了卸载，则Mises应力和PEEQ与塑性材料数据中的真实应力和塑性应变不再吻合。这里所说的“卸载”不一定是模型外载荷的减小，在加载的过程中，摩擦或模型各部分的刚度差异有可能会改变应力场的分布，从而造成模型局部的卸载。
PEEQ和PEMAG的具体含义见
http://www.simwe.com/forum/post/view?bid=31&id=691447&tpg=1&ppg=1&sty=1#692219
无论有没有卸载，积分点上的Mises应力都不会超过理想塑性材料的屈服极限。卸载只是对PEEQ有影响。Getting
Started with ABAQUS&& 10.4.6
Postprocessing the results 有解释：
"Do the values listed in the contour legend surprise you? The
maximum stress is greater than 580 MPa, which should not be
possible since the material was assumed to be perfectly plastic at
this stress magnitude. "
"The Mises stress at an integration point can never exceed the
current yield stress of the element' however, the
extrapolated nodal values reported in a contour plot may do so. In
addition, the individual stress components may have magnitudes that
exceed the value of the
only the Mises stress
is required to have a magnitude less than or equal to the value of
the current yield stress."
另外，帮助中的这个例子还说了，理想塑性不容易收敛，最好加入硬化特性，让材料曲线中最大的应力应变大于模型中可能出现的应力应变。
作切削时候遇见这种情况怎么办？？？
刀具最后有点浸入工件里了|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent0@6d
但是还在切削 不知道怎么回事
模型做得真不错。工件的节点没有进入刀具，只是由于单元变形太大，所以单元位于刀具里面，可以试试：
- 尝试其他的接触算法，比如kinematic contact method, Balanced master-slave
contact (即把weight factor设为0.5). 见Getting Started with
ABAQUS& 12.7 Defining contact in
ABAQUS/Explicit，ABAQUS Analysis User's Manual&
21.4.4 Contact formulation for ABAQUS/Explicit contact pairs。
- 裂开的地方网格细化，比如工件沿刀尖圆弧分布10个单元。
- 刀具用analytical rigid，如果用discrete rigid,
刀尖圆弧处也要增加单元（比如10个）。见Getting Started with
ABAQUS&&&&
12.8 Modeling considerations in ABAQUS/Explicit。
- 减小刀尖圆弧的半径。
问：关于切削加工仿真的，
你的模型里网格没划分，接触对Int-1也没定义，根本就运行不了。
接触面不能只是刀头的一个点吧，应该是刀的两条折线。explicit接触问题参见Getting Started with
ABAQUS&&&&
12.7，12.8
想模拟材料破坏要定义失效判据，并让切下来的单元脱离模型，见ABAQUS Analysis User’s Manual
第11.6节Progressive damage and
failure。切削问题很复杂，建议你先做一些简单的explicit接触问题。
问：真诚求助高手指点－切削 网格变形
abaqus 6.5& 的 adaptive mesh 并不会细化网格，见
http://www.simwe.com/forum/thread--1.html
可以不用adaptive mesh ，自己把网格划细些，像
http://www.simwe.com/forum/thread--1.html 一样。 你的网格太粗了。
可以先把坯料长度改短些，以缩短计算时间，调试成功后再改为真实的坯料长度。
刀尖上要有圆角。关于接触的定义见&
http://www.simwe.com/forum/thread--2.html
ABAQUS常见问题汇总（1.0版）
Explicit中的接触问题
材料非线性实验与FEM的比较
**请教**材料非线性实验与模拟的差异问题!
我做了一些复合材料的层间动态压缩实验和模拟，实验测得的应力应变曲线表现出明显的非线性，见附图，按这条应力应变曲线作为材料参数建立有限元模型后，计算出的应力应变曲线和实验曲线对比，在失效应力前基本吻合还算可以，可是超过失效应力后，实验测得的仍是非线性的，应力下降应变继续增加，可是模拟出来的结果却是突然卸载，想请问一下为什么？是不是ABAQUS软件本身认为只要出现应力下降就是完全失效然后不再承载呀？看到的一些材料模型有理想弹塑性、弹性+塑性强化段，对于这种应力应变曲线模拟不出来吗？
iy1R‑Yw|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
5f9h2U1MD仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent另外问一下，材料弹性模量是不是不能太小呀，因为我没测准它的弹性模量，就想先选一下大概的值算一下，取在100GPa左右都可以算下去，取到20-30GPa的话就出错，原因是变形速率是波速的一点几倍，然后就出错算不下去了，网络分的感觉已经足够细了，还是有问题，在这也一块请教一下大家。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/Iq'?+B"Up
8Qjw)~6v0{D+ql
谢谢大家参于讨论！
估计是这样：
试验曲线是整个板的结果，板的一部分失效了，另一部分还没实效，能继续承载，所以应力不会一下子降下来；
你图中的FEM结果是一个单元的，这个单元失效后就一点也不能承载了，它的应力一下子降下来。
可以比较一下试验和FEM的载荷或支反力结果，它们都是基于整个板的。
支反力RF是“支承处的反作用力”，有边界条件的地方才有RF. 如果边界条件定义在整个面上, 可以用distributing
coupling把加载面和一个参考点耦合起来，然后把边界条件加在参考点上，计算后看参考点的RF。
透射杆在子弹射入方向（假定是方向2）上是固定的吧，施加了此边界条件的节点上的RF2之和，就是试样受到的载荷，对称轴上的RF没用。不用coupling也行。
接触问题技巧整理
(cafe) 【分享】ABAQUS
培训资料 － 接触单元部分－Lecture1
如何加入传统的接触单元，如：goodman单元、desai单元？
答：不像ANSYS，
ABAQUS不须定义接触元。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM5Hw2cnk‑q*t(@[1]{u
用*SURFACE在3D单元的面上定义SURFACE，然后用*CONTACT PAIR定义哪个面或哪个面接触。
请教aba_aba，如何设置库仑摩擦接触面属性，不知道在什么地方设置该接触面切向摩擦属性中的内聚力c值？
Getting Started with ABAQUS&
&附录 C.9.1 Creating an interaction property
其中用的是无摩擦的接触属性。库仑摩擦则需把 Mechanical / Tangential Behavior
中的frictionless改为penalty, 输入摩擦系数 friction coeff .
内聚力不知道，没用过。
[讨论] 关于接触分析中的收敛问题。
我做了一个圆柱压一个平面的接触分析，采用了两种方案
现在的问题，我发现在无摩擦计算中， 最后一步的残余位移和力都非常小： LARGEST RESIDUAL
&& & 3.245E-09，
LARGEST CORRECTION TO DISP.&
-4.941E-15www.simwe.com*x&xL[1]^7M.M]L&
www.simwe.com9T‑K4L
而在有摩擦的计算中，最后一步的残余位移和力都很大啊：仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMi#DB^L
LARGEST RESIDUAL FORCE&
|Simwe.com|仿真|设计|有限元|虚拟仪器,mKU
L.N?&XA2G$V‑N
LARGEST CORRECTION TO DISP.&
&&&3.786E-04www.simwe.comJ
t;a*lchU(W)w6`
请问这个是不是表明在有摩擦的计算中，目前还没有收敛啊？
收敛准则用默认的就可以。根据Getting Started With
ABAQUS&&8.2.2
节&&“ABAQUS/Standard中的平衡迭代和收敛”：
“CRITERION FOR RESIDUAL
FORCE&&FOR A NONLINEAR
PROBLEM：&&5.000E-03“ 的含义是
RESIDUAL FORCE 和 TIME AVG. FORCE 的比值小于 5.000E-03
即可，而不是&&RESIDUAL FORCE 的绝对值小于
5.000E-03 。
你的计算中 TIME AVG. FORCE =
2.302E+04，&&RESIDUAL FORCE =
做接触分析时，应观察msg文件中的以下信息
（接触状态的变化：1个节点由不接触变为接触，10个节点由接触变为不接触）
CONTACT CHANGE
SUMMARY:&&&&
1 CLOSURES&&&
10 OPENINGS.
SEVERE DISCONTINUITY
ITERATION&& 2 ENDS
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
4 OPENINGS.
如果分析能够收敛，CLOSURES和OPENINGS的数目会逐渐减少，否则一般是接触定义不当，可以考虑以下解决方法：
●& &
检查所定义的接触面、接触参数和边界条件是否正确；
●& &
在静力分析中，必须定义足够的约束条件，以消除各个平移和转动自由度上的刚体位移；
●& &
避免过约束（overconstraint）；
●& & 合理地定义接触面、接触参数和过盈配合；
●& & 使用足够细化的网格；
●& &
不要在接触面上使用C3D20、C3D20R和C3D10等单元；
●& &
在接触对上设置微小的过盈量，以保证在分析的一开始就已经建立起接触关系；
●& &
施加临时边界条件，以保证在接触关系建立之前，模型也不会出现刚体位移；
●& &
将分析过程分解为多个分析步来完成，让各个载荷分别在不同的分析步中逐步施加到模型上，避免使接触状态发生剧烈的改变。
下面按照应用领域对回帖分了类，但其实定义接触的原则是相同的，例如岩土类问题的接触问题的解决方法，往往同样适用于金属成型中的接触问题。
问：此类问题,如何在接触分析中约束?请指点!SimWe仿真论坛 U]I,wC%Y
,q$S*g2Yo.O@E
J/{^fT h-j1u-kTO9x
另外,在用壳模拟分析时,abaqus默认的是不是中面?能否更改?即除了材料方向外,如何更改厚度方向?www.simwe.comF$mn‑}&`&y*]Se
''接触分析中约束'': 不明白你的问题。
对称模型，只需取1/4分析，这样X和Z方向就有了足够的约束。
在加载荷之前，增加一个分析步，使part-1-1向上移动0.001，不加载荷，从而使接触关系平稳地建立起来。在下一个step中，去掉此边界条件，加载荷。
壳模拟分析： 见&
http://www.simwe.com/forum/thread--2.html
ABAQUS常见问题汇总（1.0版）
16.2& 用OFFSET参数偏置shell的中性面
面上的网格密度
问：如何得知接触表面网格太密或太疏而造成不收敛的问题呢..仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent-U&
www.simwe.comu3K Sk*Z2T
掌握 master 面网格密度比slave面网格密度疏，此点足够吗？SimWe仿真论坛7OZ}z8Q'[;v6S
[1x`仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
M.V(@%` Z& & &
}J&H+}接触分析中若有大变形情形，因为很不容易收敛
e3_zy1oZSimWe仿真论坛仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM"]{.I2j)y&HAw
改变 increment size 中的 minimum
值，将其值设小（default：1E-005）SimWe仿真论坛7f)YN:Ch&NN*X}
]#?m0l!`:U#S4}& &
c6b仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM能提高分析成功的机率吗？
网格疏密要靠经验，比如接触面在某个方向上只有一两个单元，就是太疏了，如果接触面在某个方向上有上百个单元，就可能是太密了。
大变形时应该减小increment size 中的 minimum
值，但如果接触或塑性材料参数定义得不正确，还是不能收敛。
为什么计算中断
求救：为什么我的模型提交JOB后计算一半就中断，提示错误：too many attempts made for this increment 可能是不收敛，减小最小增量步的值只能再多计算几步，但还是算不完。各位大侠请指教啊！
Tie的两个面之间的距离大于你定义的position tolerance = 0.0001, 所以tie没有完全起作用，骨头可以乱动，ABAQUS找不到合适的接触状态。
用display group来单独显示surface，例如 c4-1."Surf-c4 i" 和 i45-1."Surf-i45 s"，可以看到它们的距离大于0.0001。
在后处理中用view cut, 也可以查看各个面之间的距离。
有两种解决办法：
- 把 i45和i56 改薄些；
- 移动红烧排骨之间的相互位置，但这同时会改变接触面之间的距离。
Tie的 position tolerance 应该定义得大些，例如一个单元的大小，同时不要选Adjust slave surface initial position, 以免调整节点位置造成畸形单元。
接触面定义得太大，也会造成收敛困难。应该只把可能接触的位置定义为接触面。
c5-1."Surf-c5 p" 既是主面，又是从面，可能会有问题。应该分别定义两个面。
C3D4单元精度很差，应该用二阶单元 C3D10M。接触问题不能用二阶HEX单元，但可以用C3D10M。
在step 1中不要施加实际的载荷，而是施加一个3方向的微小的力，使接触关系建立起来，在step 2中再施加实际的载荷。
接触部位的网格划细些，否则只有几个单元发生接触，就不容易收敛。在后处理中察看接触面之间的距离（field
output的COPEN），可以看到有几个单元发生接触了。
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
如果位移远远小于一个单元的长度，可以用nlgeom = off.
如果还不收敛，可以增大 “允许的最大迭代次数”，见
http://www.simwe.com/forum/viewthread.php?tid=731009
http://www.simwe.com/forum/viewthread.php?tid=693237
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
关于接触产生节点应力集中的问题...请大家帮帮忙~拜托
我使用一只脚的模型
+zKLu& & &
v‑ywww.simwe.com使用静态分析方式
F&r!hSimWe仿真论坛将脚模型固定,地面向上撞击脚底SimWe仿真论坛Q'EGp
将地面设为master surface,而脚底设为slave
surface(网格较细).SimWe仿真论坛3z#|j& &
脚底与地板之间的接触条件是用surface to surface contact
7D,{s4&nV%OSimWe仿真论坛之间摩擦力设为0.11www.simwe.com1%u9K2VE
其它为预设
V;D3j|Simwe.com|仿真|设计|有限元|虚拟仪器可是结果如图
奇怪的应力看来是因为你的脚面本身就是凸凹不平的，凸的地方接触压力大，凹的地方接触压力小。
在后处理时无法正常察看CPRESS,&
说明接触面定义得不对（你的模型是用其它软件建的吧），你在ABAQUS/CAE里重新定义一下contact的salve
surface和master surface.
没必要把整个脚的表面定义为slave surface, 只定义足底可能接触的部分即可。
下图中，正中间的节点只属于下面两个单元，不属于上面的单元，这样的网格精度很差。每个节点都应该属于周围所有的单元。
单元类型C3D8的精度很差，应该用C3D8I或者C3D8R。
Step的increment size不要用fixed，这样出现收敛问题时，无法自动减小increment size。
脚后跟处应力很大，如果关心此处的应力，应该细化单元。
足底和地板接触前的应力很小，没必要模拟接触前的过程。可以一开始就让足底和地板接触。
如果想模拟瞬间的撞击，要用explicit作动态分析，
在后处理时，用display group 来只显示所有surface, 查看cpress和copen, 看到可能导致不收敛的原因是：
有3个接触对的master surface定义得太小，变形后接触点将要移出master surface之外 （你加4方向正转矩时能收敛，可以看看此时是不是没有这种问题）。
另外，各个面上都只有一个节点上发生接触，这可能会有问题。应该在此处细化单元，并适当增大adjust的值，使更多单元发生接触。
任务提交后提示以下错误：“Too many attempts made for this increment”
在msg文件中看到
CONTACT PAIR (ASSEMBLY_BLANKBOT,ASSEMBLY_TIE-1_DIEDURF) NODE
BLANK-1.5 IS& OVERCLOSED BY 0.0512228 WHICH IS TOO
这往往是因为接触面的法线方向定义反了。定义刚体和shell的surface时，要注意选择外侧。在后处理中可以显示法线方向：ABAQUS/CAE
Manual&&&&
37.7 Displaying element and surface normals。
问：请问， 过盈配合怎么做， 用哪个菜单？
----------------------------------------
*CONTACT INTERFERENCE
----------------------------------------
此关键词的使用方法为：
*CONTACT PAIR, INTERACTION = &接触属性的名称
&, ADJUST = &位置误差限度
&从面名称 &，&主面名称
*AMPLITUDE，NAME = &幅值曲线的名称 &
*CONTACT INTERFERENCE，AMPLITUDE = &幅值曲线的名称
&从面名称 &，&主面名称
&，&过盈量或间隙量 &
其中，参数&过盈量或间隙量 &为负值表示过盈量，正值表示间隙量。使用
INTERFERENCE类似于施加载荷，不能在initial分析步中对其进行定义，而只能在后续分析步中定义。
☆ ABAQUS/CAE操作： Interaction模块，主菜单Interaction → Create，点击Edit
Interaction对话框底部的Interference fit。
如果在分析结果中看到，使用 *CONTACT
INTERFERENCE所定义的过盈接触没有在模型中起作用，则有可能是以下原因引起的：
● *CONTACT
INTERFERENCE不能使用ABAQUS默认的幅值曲线Ramp（从1降至0），而要使用自定义的幅值曲线，使过盈接触的幅值在整个分析步中从0到1逐渐增大。
● 如果在ABAQUS/CAE模型中两个接触面之间有宽度为t的缝隙，而事实上它们之间应该存在接触关系，则关键词 *CONTACT
PAIR中的参数 ADJUST = &位置误差限度 &
必须略大于此缝隙的宽度t，否则ABAQUS会认为这两个接触面没有接触，*CONTACT INTERFERENCE 中的参数
&过盈量或间隙量 &不会起作用。
● 参数&过盈量或间隙量
&为负值才表示过盈接触，正值表示两个面之间存在缝隙。
问：使用clearance和使用CONTACT INTERFERENCE算出来的结果有什么区别么？
*CLEARANCE也可以，但有一些缺点，它只能用于small
sliding，不能在ABAQUS/CAE中定义，不能在step中修改过盈量，而且没有使用ramp,
也就是说，在分析的一开始全部过盈量就会被施加在模型上，如果此过盈量太大，ABAQUS就可能无法达到收敛，而且即使减小时间增量步也没用。
问：我的两个接触面是经过运动后才到达相互接触的状态的，但这时并不能达到理想的接触（两个面有部分gap,很小但的确没有接触，用contact
print 打印出来slave nodes相关项，都是zero)。我想在这个时候来adjust the contact
gap,让abaqus把接触面间的gap调成没有。
因为一开始时两个接触面的间距非常大，所以不可能在一开始就调整，而且我还需要对那个运动进行分析。
使用 *CONTACT INTERFERENCE可以在分析步中定义接触对之间的距离或过盈量，你可以试试。
使用 *CONTACT INTERFERENCE时要注意：
1) 其默认的幅值Ramp是从1降至0，不能使用，而要使用自定义的幅值曲线，使幅值在整个分析步中从0到1逐渐增大；
2)& 正值是缝隙，负值是过盈量。
问：我所关注的只是螺钉拧入之活的应力分布，而对于拧入的过程则不关心。我想不放把问题说的更细一些。
把螺钉拧入骨头可分为如下几种情况：
A、自攻型，就是硬生生的把螺钉拧入骨头。
B、助攻型，先在骨头上打一个比螺钉的直径小一些的洞，之后再把螺钉拧进去。这样做的好处是螺钉不容易折断。打不同直径的洞，在螺钉拧入后的应力分布应该是不同的。
C、用一个同样的螺丝钉（不锈钢）在骨头中钻个洞，再把需要的螺钉（钛）拧进去。
D、或者有时用一个稍小一点的不锈钢螺钉钻洞，再拧入钛螺钉。这样可以比C更稳固。
E、我不小心把洞钻大了，这样螺钉就容易松动。当然这在临床中是不允许的，但也存在。
各位大侠能否就这几种情况分别讨论一下，我觉得是不是可以用过盈装配来模拟？但具体A、B和D有什么样的区别呢？想不明白。等待着各位大侠解惑~~~
这么说我就明白了，弄明白了就感觉很爽。
一定是兽医，把骨头这么搞来搞去的，让病人求生不得，求死不能，关云长来了也受不了啊。
A（怒射破门）就是过盈配合非常大。
B（后卫助攻）就是过盈配合小一些。
C（换上替补）就是几乎没有过盈配合。
D（影子前锋）就是过盈配合也比较小。
E（自摆乌龙）就是接触面之间存在缝隙。
具体过盈配合该多大，估计该在螺钉上装些传感器，得到试验结果，再调整你的过赢配合量，与实验相吻合。过盈配合的模拟方法见。
zero4:& 用ABAQUS软件处理管土相互作用中的接触面问题
给有兴趣的朋友一个参考。
H%{:M(|&vX
附件: 1.part1.rar ( 13:36, 195.31 K)
问：计算中止，.msg文件里显示
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING
INSTANCE PART-1-1 D.O.F. 8 RATIO = 6.49141E+11
检查约束没有错误，而且节点346位于模型中部，不知道是什么原因引起的？望高手指点。附上文件，太大分了2部分。
桩上没有边界条件，完全靠接触面上的摩擦来消除刚体位移，但接触关系是在迭代过程中慢慢建立起来的，不是一开始就有摩擦力，所以会出现刚体位移。见
http://www.simwe.com/forum/post/view?bid=31&id=680814&tpg=1&ppg=1&sty=1&age=0#680814
接触对的主面上不能有尖角，桩的两侧要各定义一个接触对，底部用tie. 见
http://www.simwe.com/forum/viewthread.php?tid=699614
http://www.simwe.com/forum/viewthread.php?tid=679059
http://www.simwe.com/forum/viewthread.php?tid=684042
http://www.simwe.com/forum/viewthread.php?tid=679059
中谈了过约束的问题：
"msg文件中出现zero
pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact,
出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region,
然后选择区域时不要包含这一圈节点。"
你的桩底部只有一个单元，就没法用上述方法。要在底部细化网格才行。
不必整个模型都用那么细的网格，离接触面较远的地方可以粗一些。
新手请教：The plasticity/creep/connector friction algorithm
did not converge at 12 points
模型做得不错啊，改为弹性材料就可以收敛，说明模型没有太大的错误。我认为需要改进以下方面：
应力应变变化很剧烈的部位（桩底面上方和下方的的soil）,单元不能太狭长。原则是在后处理时不能看到单元之间的应力和PEEQ不连续（颜色发生跳跃）。
桩侧面接触状态变化很大, soil和桩的网格沿圆弧方向和纵向都要再细化一倍。
对称模型取1/4就够了。用C3D8R单元可以缩短计算时间。
rebar是truss单元，不能用*Solid Section,要用
&&*Solid Section, elset=***,
material=REBAR
&&截面面积
桩底部的从面网格还是比主面粗糙,会有问题。桩底部可以不用接触，改为tie，收敛就更容易：
&&*Tie, name=Tie-Bot, position
tolerance=0.1
&&soil-1.pilebot,
mpile-1.bot
我做了个轴对称模型，rebar相当于shell, 没有*initial condition,
见附件。从结果中看到，虽然应变很大，但位移很小，因此可以不用在step中设 NLGEOM=yes，这样收敛就更容易些。
塑性问题要把step 2的初始时间增量步和允许的下限设小些，比如：
&&0.01, 1., 1e-06, 1.
下面的警告信息都是说塑性问题不收敛：
ELEMENT SOIL-1.709 IS DISTORTING EXCESSIVELY.
***WARNING: THE PLASTICITY/CREEP/CONNECTOR FRICTION ALGORITHM DID
NOT CONVERGE AT 5 POINTS
Soil-Axisymm.cae(67.54k)
问：大侠们看看这个初始应力为什么总是施加不成功了？？
建议在step 1中只保留桩底部的接触，去掉其他接触。step 2中再加上所有接触。先做弹性材料，收敛了再做塑性。
问：求教：为什么基坑开挖中围护结构的作用没有体现？（附有inp文件）
简单介绍一下：通过读入初始应力的方法计算地应力平衡；step1中杀死所有的围护结构（包括3个支撑梁Graider1,Graider2,Graider3、一个连续墙Wall）；step2中激活wall、开挖第一步所要开挖的土体kw1、激活Graider1。后面的暂时没有计算。可是计算结果是wall和Graider1好像没有起作用。那位高手能帮我指点一下，先谢谢了!
围护结构和土体之间要定义接触或*rebar或*embeded，否则它们仅仅是空间位置重合，不会自动发生相互作用。
问：隧道开挖与管土相互作用中的接触问题求教
本人做了隧道开挖中的管土相互作用的模型。有些管－土接触问题求教。|Simwe.com|仿真|设计|有限元|虚拟仪器h&o},B0}(r%m:NQ
模型说明：开挖隧道位于管道正下方，由于是研究接触算法，仅做三步开挖，且没有进行支护。土体使用实体单元c3d8r，管道采用三维壳单元s4.模型1管土接触采用法向www.simwe.comkXP%z
硬接触，切向罚函数（文件名exva-h1-fenbu）。模型2采用tie约束管道与土体节点。（文件名tie-ex-fenbu)
&IfZD~'ev2.
计算后出现的问题SimWe仿真论坛vp2W&W8~E
（1）开挖隧道纵向中间部位土体时，两个模型的管道位移基本一致，接触算法作用不大。www.simwe.comY:K6Q
AQ[1]i8z,r
（2）隧道全部开挖的情况下：使用tie约束的模型2计算不收敛，通过查看塑性变形，发现和模型1的差不多，不知道为何不收敛。请问这种隧道开挖引起土体位移对管道的影响问题，设置接触的必要性以及模型2为何不收敛。
#k*B&k,~3Mwww.simwe.com仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMZe
[ 本帖最后由 doudout 于
或许是定义tie的时候应该加上 position tolerance=0.1, 以确保两个面tie在一起。
你的单元是c3d8, 精度很差，应该用c3d8r。
单位对吗？如果长度单位使用mm，那么质量的单位必须是“吨”，相应的密度单位就是吨/mm3，应力的单位是N/mm2（即MPa）。
见 http://www.simwe.com/forum/viewthread.php?tid=684273
问：做基桩的水平静载试验模拟，MSG文件提示迭代增量步太大，不能进行
取1/2原型，近场土C3D20R，远场CIN3D8，对称面上Y向对称约束。请赐教！
答：我认为模型中有以下问题：
1。塑性材料和接触面上都不能用C3D20R和C3D20单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用C3D8I
(在所关心的部位要让单元角度尽量接近90度)，如果只关心应变和位移，可以使用C3D8R,
几何形状复杂时，可以使用C3D10M.
我尝试把你模型中的接触关系以及材料参数中的*Mohr Coulomb 和 *Mohr Coulomb
Hardening都去掉，变为弹性材料, 仍使用C3D20R就收敛了。
2。接触对中的slave surface应该是材料较软，网格较细的面，你正好颠倒了。
接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于接触面之间的距离，否则ABAQUS会认为两个面没有接触：
*Contact Pair, interaction="SOIL PILE SIDE CONTACT", small sliding,
adjust=0.02
在ABAQUS/CAE中可以在edit interaction对话框中设定此参数。
4。定义tie时也应该设定类似的position tolerance:
*Tie, name=ShaftBottom, adjust=yes, position tolerance=0.1
Tie和contact pair都定义在step之前，你现在的模型是对的。
overconstraint警告信息如果只是出现在dat文件中，没有出现在msg文件中，就没问题，说明ABAQUS自动解决了过约束问题。
如果overconstraint警告信息也出现在msg文件中，说明ABAQUS无法自动解决此问题，这时分析往往不会收敛，在后处理时可以用display
group显示出现过约束的node set (ABAQUS会自动为这些节点建立node set),
需要你自己修改模型，避免过约束。
附上修改后的文件，运行时没出现零主元（zero pivot），我做的改动是：
- 去掉所有的 *Mohr Coulomb 和 *Mohr Coulomb Hardening
- 在*Contact Pair中添加 adjust=0.02
- 在*Tie 中添加 position tolerance=0.1
还应该改的是：
- 接触问题不能用C3D20R单元（所以目前施加载荷时很难收敛）。
- 接触对中的slave surface应该是材料较软，网格较细，即应细化SoilInnerSurface面。
msg文件中出现zero
pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact,
出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region,
然后选择区域时不要包含这一圈节点（我附上的文件中没有做这样的修改）。
无穷元我没用过，不太了解。
接触定义在哪个分析步取决于你模型的实际物理背景，如果从一开始两个面就是相接触的，就定义在initial或你的第一个分析步中；如果是后来才开始接触的，就定义在后面的分析步中。边界条件也是这样。
你现在把接触定义在后面的分析步中，应该是可以的。
你各个分析步的设置也没看出有什么问题，只要符合实际物理背景，并且利于收敛即可。
Data5.rar (94.68k)
岩土问题我没搞过，回答不了你，论坛上好像有这方面的资料，你搜一搜。
单元扭曲的模型传上来看看。
不太了解你的单向多循环加载法。原则上，如果是静载荷，用static分析步即可，如果加载时间很短（例如在地震中），或载荷本身的性质是动态的（例如来自旋转机械的载荷），在分析中不能忽略结构的惯性，就要采用动态分析，可以用ABAQUS/standard或explicit.
我在前面上传的文件里用*CONTROL设了允许的迭代次数18，意思是18次迭代不收敛时，才减小时间增量步（ABAQUS默认的值是12）。一般情况下不必设置此参数，如果在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛，就可以用*CONTROL来增大允许的迭代次数。
你的C3D20R单元导致了接触问题的不收敛，所以我用*CONTROL其实是没用的。
你的Test.inp也不能收敛啊。另外，如果你传上来两个模型，一个收敛，一个不收敛，应该写清楚这两个模型有什么区别，才好帮你分析。
你的桩用C3D15单元，是15-node quadratic triangular prism,
这还是一种二阶单元，不能用于塑性材料和接触问题。在TestModifiedContact.inp中把材料改为弹性，step
1就收敛了（step 1中没有接触, 但step 2有接触，又出现过度约束问题，前面已经讲过解决方法）。
你的桩还是应该像你以前的模型那样划分单元，用C3D8I. 另外，
尽管说主面的网格要比从面粗，但像你现在这样，主面上有很大尖角，还是可能出现收敛问题。原则上，每90度的圆弧上应该划分10个单元，适当少一些可能也行。
pair中的adjust=0.005，还是太小，在后处理时可以看到，接触面之间的距离大于0.005。把adjust设置大一些没关系，比如adjust=0.1。
在施加载荷前，应增加一个分析步，不施加载荷，只是让接触关系建立起来。原则是不要在一个分析步里让模型发生太多变化，那样就把ABAQUS搞晕了。
施加载荷时，如果是塑性材料，初始时间增量步应设小些，比如
0.1, 1., 1e-05, 1.
关于对称面上应力和位移，在后处理时用display group隐藏模型的一半即可。
问：仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentO&x&
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
整個被壓扁了& & 不如預期&
& [1]a3K]5fS‑q/M
照理&&是要&
&基礎&&沉入&&土壤內&
|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent9ZO"X‑uH@
或是&&對土壤&&造成&
&沉陷的效果& &^^
msg 文件中出现overconstraint警告，说明出现过约束。去掉rigid body约束就可以了。
foudation上要加对称边界条件。
长度的单位是什么？重力9.8太大了，和弹性模量的单位要统一。
土体不接触的地方单元可以划粗一些。
你的基础长度100，是100m吗？
数量级还是不对，用国际单位的话，以钢为例，弹性模量是2.1e+11 N/m2, 密度是7800
kg/m3.& 弹性模量如果太小，变形会非常大，难以收敛。
如果想模拟把基础压入土中，需要把土定义为塑性材料，并在接触面附近划分非常细的单元（塑性问题很复杂，你先把弹性作对了，再做塑性）。见
http://www.simwe.com/forum/post/view?bid=31&id=683533&sty=1
http://www.simwe.com/forum/post/view?bid=31&id=683101&sty=1
局部细化网格，要先在mesh模块中用菜单 tools / partition来分割要细化的部位。
问：这个桩基础大致是这样的, 顶部是一个直径在8CM左右的微型混凝土桩, 桩长5M左右. 桩下部有二个螺旋状的钢板,
直径为30CM左右, 厚度为1CM, 深度在地面下6到7米. 钢板焊接在圆形钢柱上, 其直径大约为4CM.案有限元模拟时,
土体为圆柱, 直径6M, 深度取为10M. 有限元分析时, 需考虑桩与土, 钢板与土以及钢柱与土之间的接触.
根据ABA-BA 大侠的见解, 土体和桩可采用C3D8I 单元, 我想请教的问题是,
(1), 本例的薄钢板到底应该采用3-D单元模拟(C3D8I)还是2-D单元(SHELL)模拟好.
(2), 同样, 细长的钢柱应该用何种单元,
(3), 因为各结构之间尺寸相差很大(比如, 钢板的1CM厚度和整个计算域10M相比), 如果全部采用三维实体单元,
整个模型大概为20万个单元, 对如此大的模型进行三维材料非线性分析并考虑接触问题, 大概需要多长的计算时间呢,
(4), 如果薄钢板采用SHELL单元, 如何设置它与土之间的接触面呢, (我的意思是, 当把一个无厚度的钢板放进土中,
如何分离出和钢板接触的二个面),
xile1981答：
1,4.如果要考虑要减少计算时间的话,我认为也可以使用shell单元对钢板进行建模,当定义接触面时使用
*SURFACE, NAME=surface-name
element-set-name, SNEG(or SPOS)来定义顶面或底面.
2.不知道你的钢柱有多长,个人建议还是用实体来建模吧,用beam的话感觉做接触效果不好.
3.计算时间是要根据你的机器来看的,没有一定的结果,以我以前计算过的模型来看,15万到20万个单元考虑接触(应该接触对要比你多,还有些材料非线性)需要一天左右,使用机器为(2-cpu,64bit,8G内存,linux系统)
问：非常感谢XILE1981大侠的回复. 关于设置接触面问题, 可能是我表达的不清楚, 我的意思是, 当把钢板放入土中的时候,
如果用无厚度的SHELL单元来模拟, 钢板上下顶面为MASTER SURFACE, 和钢板相接触的土面为从面,
这二个土体的从面如何得到呢.
是不是用钢板的厚度(1CM)在土中CUT出个空洞, 然后在空洞表面上设置从面, 但钢板SHELL
单元是通过抽取钢板的中性面得到的呢.确如果是这样, 那么从面和主面之间有间隙(钢板的一半厚度), 这就意味着主从面没有接触啊,
这如何解决呢.
xile1981答：哦,是这样啊,那你可以用
*CONTACT PAIR, INTERACTION=interaction_name, SMALL SLIDING,
来调整,使从面上单元节点调整到主面上.但可能会产生单元变形过大,但我觉得钢板厚度与网格尺寸相比较小的话,应该问题不大.
我没有做过土木,只是给你个参考吧,呵呵.
地基和土体的接触
我&&作出以後&
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3iT-EVI
試著& &加載&
&與9000N/m2&
hG:WJo仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent但&&跑出來的圖&
&不如預期& &
圖形如下:W7a+fR
P7@4P9~[1]X
0WGD&f.?|Simwe.com|仿真|设计|有限元|虚拟仪器仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'hH2ME
照理說&&應該&&是&
&基礎&&沉入&&土壤&
&如下圖仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMB;[jr5?
ol2|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMSimwe.com|仿真|设计|有限元|虚拟仪器!{[]/dN$z.~B
你那张成功的图里，应力最大值是10e7,
而失败的图里，应力最大值只有10e5。如果材料相同，应力小，变形当然就小。应力小的原因是载荷小，把载荷变为900000试试。
你那张成功的图里的deformation scale
factor是1吗？如果是的话，这样的大变形要在step里设置几何非线性nlgeom 为 on.
你用了*Solution Technique, type=CONTACT ITERATIONS，不知是否可行。
一般都是先作地应力平衡，再加载荷，不知你这样一步完成对不对。
既有隧道模拟中出现的问题
你的问题在 “ABAQUS常见问题汇总（1.0版）
8.5&&管土/桩土接触 ”
中都提到过，你仔细看过吗：
-&&***WARNING: SOLVER PROBLEM.
NUMERICAL SINGULARITY WHEN PROCESSING NODE
&& & LINER-1.24
D.O.F. 2 RATIO = 3.9.
- 网格太粗，90度圆弧上至少要有10个单元。
- 隧道在厚度方向不能只用一层单元，至少要有3层，用CPE4I单元。
- 塑性问题不能用二阶quad单元，soil单元要用CPE4I.
- contact要增大adjust的值，可以用 adjust = 0.5 .
板料成形的接触问题
请问厚度变化怎么看
在板料冲压过程中，冲压过后，我想看看板料的厚度变化情况，怎么看？|Simwe.com|仿真|设计|有限元|虚拟仪器]b}&V1yJ
我在输出变量设置里面看到有thickness，可是设置后，好像那个不是指板料的厚度
Hc^‑Za‑O1U仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&谢谢大侠指教
体单元没有法子直接输出厚度。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM"t&
Ex0P5i3A+A
其实体单元已经很直观了。
&A5E$~.x&_
y5M仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM可以查看厚度方向上下两个节点间的DISTANCE,www.simwe.comO'_/H*g
当然可能对不正，要采用一定的策略配准。
-y'Es5i&N({h厚度随时间变化的曲线某一点处的还是可以输出的。www.simwe.com&Nr'i#J
很复杂。。。。
板料成形还是壳单元比较划算些，
MtV4BXz;{[1]f‑K仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent一来可省大量的计算成本：如果要得到比较精确的结果，实体单元要划得比较细才行。
h3|[1]ksI9N|Simwe.com|仿真|设计|有限元|虚拟仪器如果用减缩单元，厚度方向要至少划四层才行，而别的方向上单元的尺寸要相对差不多。
h}-J"`就是差不多的一块板子，划实体单元要多用十几倍的单元，再加上最小单元长度变小了好多，又是十几倍的时间，如果用EXPLICIT，相同的问题，要多耗掉100多倍的时间，10多倍的磁盘空间和内存，是很不现实的。
0}(q;u._6Cf1y|Simwe.com|仿真|设计|有限元|虚拟仪器二来用实体单元好象模拟板材起皱有困难，理由见上。|Simwe.com|仿真|设计|有限元|虚拟仪器)e[1]r
三模拟回弹的时候约束不好加，也不容易收敛，至少在我看来是这样|Simwe.com|仿真|设计|有限元|虚拟仪器E)s7n+qa"x#elhN
除非你的板料相对过厚，或者有特别的目的。
问：我在模拟一板料冲压成形过程中，采用sc8r单元，但模拟第一步的0.9999时就出现不收敛的情况，这个问题怎么解决，*.inp文件在附件里面
在msg文件中看到：
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING
NODE SHEET-1.49 D.O.F. 3 RATIO = 4.5.
这说明SHEET-1在方向3上出现了无限大的刚体位移。
你的SHEET-1在方向3上没有边界条件，完全靠接触面上的摩擦来消除刚体位移，但接触关系是在迭代过程中慢慢建立起来的，不是一开始就有摩擦力，所以会出现刚体位移。
可以好好研究一下Getting Started with
ABAQUS中板料冲压的例子，其分析步和边界条件设置得很繁琐，但有其一定的道理。原则是，在第一个分析步中用临时的边界条件来消除所有刚体位移，让接触面上出现一定的压下量或过盈配合，使摩擦力建立起来，在后面的分析步中再去掉临时边界条件。
我试了试，确实很麻烦，不过还是搞出来了，见附件。看来主要还不是sc4r单元有什么特殊性，而是你模型中的其他问题。我做了以下修改：
- 在STEP 1中，固定坯料的整个左端面，这样就不会出现NUMERICAL SINGULARITY
警告信息。在后面的分析步中再去掉此边界条件。（如果只固定此面上的一条边，就不能收敛，很奇怪）.
- 你材料参数的屈服应力是1.37E+10，我想应该是1.37E+8吧。
- 在STEP 1中，让Holder下移0., 对坯料产生微小的压力，以便产生摩擦力，消除刚体位移。
- Holder刚性面上有尖角，在此角点上的法线方向是不唯一的，会造成收敛问题，将其改为圆角。
- 为die和坯料之间的*contact pair加上adjust=1e-07，否则不能收敛。
- 在STEP模块中，菜单other / General Solution Control,
把所有step的*control都恢复为默认值。点击菜单Model / Edit Keyword，在step
1中添加以下语句：
*CONTROLS, PARAMETERS= TIME INCREMENTATION
,,,,,,18,,,,
这是maximum number of severe discontinuity iterations allowed in an
increment. Default is 12, 现在改为18。
diagnostics中看到，接触状态不断变化的是坯料右端底部，可能是这里只有一侧发生接触，造成收敛困难。因此把die的接触面改短，这样还可以缩小接触区域，大大加快计算速度。
- *shell section用simpson积分，9个积分点。这是参考ABAQUS Example Problems Manual
1.3.3 Stretching of a thin sheet with a hemispherical
punch。你原来的Gauss积分我没试。
（经过以上改动，step 1能收敛了）
- Step 2中，不应把holder的所有边界条件都去掉，而应只去掉对U2的约束。
- BC和载荷的amplitude都用默认的Ramp.
（经过以上改动，step 2能收敛了）
- 在step "Contact 2"
中，对坯料施加一个向上的拉力，使其与punch贴在一起，否则无法收敛。下一个分析步中再去掉此拉力。
（经过以上改动，step "Contact 2" 能收敛了）
- 网格在长度方向太粗，大变形时只有很少的节点与punch相接触，导致收敛问题。在长度方向将网格细化，宽度方向改稀疏一些。
- 去掉了*potential, 不过或许有它也可以。
可能还应该改的是，摩擦系数0.162是否太大了？有润滑的话应该是0.1或更小，这样收敛就更容易些。
目前看到的结果是punch下压了分析步的0.05，后面的还在算. 等不及后面的结果了，明天再说吧。
问：在板料成形中出现的错误
在sta文件中看到
Maximum adjustment for initial overclosure is: 5.00001E-04 at
node 1 of instance BLANK-1.
***WARNING: For contact pairs that are not tied, unresolved
&&&&&&&&&&&
overclosures will result in initial accelerations on the order
&&&&&&&&&&&
the penetration distance divided by the time increment squared,
&&&&&&&&&&&
which can lead to significant mesh distortion. These
overclosures
&&&&&&&&&&&
can be avoided by user adjustment of the mesh.
后处理中看到increment 0时坯料就变形了，原因是上述overclosure.
出现这种过盈接触的原因是shell的厚度是0.001,
而模型中刚体和shell的距离是0，即过盈量为0.0005。我认为，解决方法是把刚体移开一些，使其和shell的初始距离为0.0005。见
ABAQUS Analysis User's Manual&
21.3.2 Surface properties for general contact
21.4.2 Surface properties for ABAQUS/Explicit contact pairs。
这样改后就可以运行了，但在2.8e-3秒时，看到坯料局部变形极大，原因可能是holder
force太大，材料不能自由流动。
另外，你用了刚塑性材料，没有硬化特性，不知是否会有问题，见Getting Started with
ABAQUS&& 10.4.4 Adding hardening
to the material model。
请教一个*control的问题
在用abaqus进行模拟时，如果不收敛，可以通过输入文件中用*CONTROLS关键字来更改默认的迭代及收敛准则。在看了好多遍keyword之后，对其中的参数还是高不清楚，希望高手给解释一下，每个参数都是什么意思啊？谢谢了！
用*CONTROL来增大 “允许的最大迭代次数” 不一定能解决问题，只有在以下情况下才有用：
在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛。
我一般只修改 “允许的最大迭代次数”，其他收敛参数修改后效果好像不大好。
具体方法见&
http://www.simwe.com/forum/viewthread.php?tid=693237
问：模型为二维弯曲的一半，板材为三层板。当将层与层中间的inter接触性质改为tie时，就可以收敛，但是我想计算一下在超过一定剪切强度下的摩擦情况。也就是说当界面剪切强度超过1E6Pa时，发生摩擦，摩擦系数为0.1的情况，大家帮我看看，我该怎么调整（Interaction
中的tolerance已调整为1E-8）
中间层为什么用三角形单元，那么小的尖角，单元形状太差。
尝试把材料都改为steel, 线弹性，摩擦属性不用taumax, 但仍然不收敛，说明不是这方面的问题。
问题是：接触面上单元非常多，需要增大允许的迭代次数。方法是：step模块，菜单other / Gerneral solution
controls, 点Edit，选择第一个分析步"establish contact I"，点continue,
选择specify，点Time Incrementation 标签页，点第一个more，把Is由默认的12改为30，点OK。
相应的inp文件：
*Controls, parameters=time incrementation
, , , , , , 30, , , ,
含义是：30次迭代不收敛时，才减小时间增量步。
这样改后就收敛了。另外，可以尝试把单元划得粗一些，就可以收敛得更快。
reduced integration单元要求在厚度上至少有4个单元。建议改为CPE4I单元。
问：这是个简单的拉伸成形模型，有几个问题没有解决：
（１）局部的单元变形太大，怎么控制让领近单元变形均匀点？
（２）我的回弹模型出了什么问题，一直没有成功，希望大侠能给指点一下，问题出在哪里！！
wy_anli答：
在成形过程中，那个应变比较大的单元已经超过了你定义材料的最大应力，也就是说当应力达到你定义的材料最高应力时，就会把材料当成超塑性体。
（见图中的5.856E8,你定义的最大应力为4。7429E8）
建议减小压边力
因为你采用的是显示算法，所有的加载应注意平滑加载，在计算完成后要看一下动能的历史输出，要平缓，另外比较动能与内能，动能/内能一般要小于8%
另外那个BLANK_UP_FORCE可以减少甚至去掉，因为BLANK_UP_FORCE会增大成形难度。
2.回弹要在隐式算法中做
你的那个图是变形完成后的应力情况，在成形的开始的部分应力已经超过了极限值，你可以从开始看一下应力变化情况，也可以看我传的那个图
你先把成形的部分模拟完毕（修改压边力，使成形合理），再进行回弹模拟
问：这个模型为什么不收敛（超弹和弹塑性接触）？二维模型为拉弯三层板，因为是对称的，所以只取其中一半。上左为凸模，上右为压料板，下为凹模，中间为三层板。
三层板中间为tie接触，上下两层为steel，中间那层板为hypereslatic超弹性材料。
附上inp文件
用线弹性材料试试，如果仍不收敛，就是模型参数定义得不对。见
http://www.simwe.com/forum/viewthread.php?tid=693237
http://www.simwe.com/forum/viewthread.php?tid=680814
http://www.simwe.com/forum/viewthread.php?tid=684042
http://www.simwe.com/forum/viewthread.php?tid=679059
问：（1）如果进行多道次的成形模拟，一般是怎么进行的？是不是：工序1（成形－回弹）－传递到－工序2（成形－回弹）－......这进行？
（2）第一道次的模具都可以在CAE里面完成设计等，哪后面工序用到的模具，有什么好办法添加进去：－直接在.INP文件里面修改？
先谢谢了！
可能可以在一个cae模型中把所有道次的模具都建在一起，然后在每个step中用*model
change来杀死/激活不同的单元和接触。
的接触问题
问：大家好，我现在做的是一个有压管道在外加作用物下出现凹坑的准静态分析
当然这里面存在接触（外加作用物和管外壁之间），不要接触直接加力我可以做出来，但是加了接触之后就不行了，并出现以下
警告：A slave surface definition has been defined on a face of element
1619 that contains no midface node. This may lead to convergence
difficulties. The element type will be converted and a midface node
generated automatically. The element number will be unchanged. See
the output from the preprocessor to confirm the new element and
node definitions. This message will not be repeated if subsequent
elements of this surface are converted
18 nodes may have incorrect normal definitions. The nodes have been
identified in node set WarnNodeIncorrectNormal.
The aspect ratio for 51 elements exceeds 100 to 1. The elements
have been identifiein element set WarnElemAspectRatio.
58 nodes have dof on which incorrect boundary conditions may have
been specified. The nodes have been identified in node set
WarnNodeBCIncorrectDof.
58 nodes have dof on which incorrect boundary conditions may have
been specified. The nodes have been identified in node set
WarnNodeBCIncorrectDof.
错误：118 elements have missing property definitions. The elements
have been identified in element set ErrElemMissingSection.
Analysis Input File Processor exited with an error.
在后处理中可以用display
group来显示错误信息中提到的那些集合WarnElemAspectRatio，WarnNodeIncorrectNormal，就可以看到问题在哪里。
- 接触问题和塑性材料不能用二阶quad和hex单元，圆管可以用C3D8R，但在厚度方向至少要有4个单元.
- 刚体建1/4模型没错，你定义刚体的基本方法也是对的，但*rigid body 的单元集合PICKEDSET10不存在,
需要重新定义一下.
- 施加了刚体约束就要把力加在参考点上。
- 参考点有旋转自由度，应定义边界条件来约束U1,U3,UR1,UR2,UR3。
- “118 elements have missing property
definitions”：施加刚体约束的单元也需要定义section，材料可以是任意的（比如steel）,类型是shell.
- “18 nodes may have incorrect normal definitions”：接触面的master
surface上不能有尖角，否则尖角处的法线方向不唯一。因此，主面应只选择可能发生接触的面，即indentor的底面和1/4模型的两个侧面，且三个面之间都要有过渡圆角，不能有尖角。另外90度圆弧上要至少划分10个单元。见附件。
- “The aspect ratio for 51 elements exceeds 100 to
1”：圆管底部网格太稀疏，单元过薄。不改可能也能算下去，这里不是大变形的关键部位。
- 你的塑性参数中的塑性应变怎么到了1360？到2就足够了。
- “58 nodes have dof on which incorrect boundary conditions may
have been specified”：这只是提示信息，一般不用理它。
问：刚性体与弹塑性性的接触问题
我做的是个管道凹坑成型问题，
/A&S'iX6d&t-]怎么算不了
a仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM这是我的cae文件仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent8r8cmLxo
麻烦大家给我看看仿真分析,有限元,模拟,计算
msg文件中有
&***WARNING: SOLVER PROBLEM. NUMERICAL
SINGULARITY WHEN PROCESSING NODE
&&&&&&&&&&&&
ASSEMBLY.1 D.O.F. 3 RATIO = 7..
应该约束刚体参考点的U1, U3, UR1, UR2, UR3.
管道大变形处的网格太粗。slave surface的网格要比master surface细。
C3D8单元精度很差，应该用C3D8R.
问：我做的一个素混凝土梁,四点弯曲试验.用的是smeared cracking模型,各项参数详见附件文件.
前半段计算正常,可以看见上升趋势以及峰值.问题是过了之后出现典型的突降,计算开始不收敛,最后因为步长太小而中止.
改用更小的步长,划分更小的单元,都不奏效.我自己的分析是因为中部的单元被拉破坏,从而导致计算不收敛?可是从11方向的应力看不出破坏的迹象(否则应该11方向为0),梁底部依旧有拉应力,这是怎么回事?
请参考ABAQUS Example Problems Manual 1.1.4 Indentation of an
elastomeric foam specimen with a hemispherical punch。
你的模型有几个问题：
- 一般的大变形可以用ABAQUS/Standard, 如果想模拟材料破坏就要用ABAQUS/Explicit，见ABAQUS
Analysis User’s Manual 第11.6节Progressive damage and
failure。可以查看一下PEEQ是否达到了破坏应变。
- 接触问题和塑性材料都不能用二阶quad或hex单元, 应该用CPS4R.
- 大位移大变形问题，要在step中把nlgeom设为on.
接触区域的网格要非常细化，沿着压头的180度圆弧应至少有20个单元，原则是变形后单元不能出现过度扭曲变形。可以在此区域partition一下，作局部细化。
- 接触从面应再扩大至少一倍，因为下压后需要和整个压头半圆面接触。
- 压头上不能有尖角，否则在尖角处的法线方向不唯一，会造成收敛问题。可以把压头做成360度的圆或类似于试管底部的形状。
- 你的 contact property 用SOFT/AUG. LAGRANGE,
不知是否有必要，一般用默认的hard就行。
问：阿巴，针对该模型，我调整如下，
单元换用CPS4R.
nlgeom设为on.
你建议接触区域的网格要非常细化，我开始也是这么做的，不过这个模型很奇怪，我采用非常细的网格计算还不如粗网格模型计算的好，反而因为穿透过多早早就中止了计算。
另外，我的接触相互作用之所以选择contact property 用SOFT/AUG.
LAGRANGE，主要是因为硬接触带来不少问题，所以为了计算能进行下去，采用在一定距离内发生压应力。
对于压头的形状，我也是尝试了很长时间，因为是理想刚体，不可能做成你建议的形状（假如你有办法请告诉一声，以便学习），我改来改去只有这个形状是软件可以接受的。不过我试过很多不同的网格划分，好像不是关键，因为结果都不是很好，计算全部都提前中断。
最后问一下你建议的ABAQUS/Explicit问题，我的第一个步骤是预压，采用ABAQUS/Standard，正式施压是第二个步骤，我想删除该步骤，重新建立一个，怎么找不到Explicit？6。5换了方法？
网格肯定要细化，不收敛是其他方面的问题造成的。
二维刚体的形状可以是任意的，想画圆的话，分别画四段90度的弧即可。
静态分析步后面不能接着explicit，要先删掉所有静态分析步。
问：请问一下这样的问题能不能做：下图所示为一实体的中截面，打算这样计算，先让物体在1 位置，计算完了之后移动一下到 2
位置再计算！模拟重叠宽度，由于模型的接触力不好简化，想直接利用接触来做，不知道可不可以！请高手指点一下！
移动到 &2 &位置:
&设定位移边界条件, &使x方向位移为某个值即可。
定义接触时要用finite &sliding.
&要把step中的NLGEOM（几何非线性）设为on。
可分为几个Step:
Step 1: 在 1 位置，在其质心施加一个很小的力，从而使接触平稳地建立起来。
Step 2: 在 1 位置，在其质心用一个随时间变化的力，作用一段时间（如 5s）;
Step 3: 载荷变为0，去掉接触（否则迭代过程会很慢），设定位移边界条件, 令x方向位移为某个值（即使其移动到 2
Step 4: 在2 位置，恢复接触，在其质心施加一个很小的力，从而使接触平稳地建立起来。
Step 5: 在质心作用同样的力。
Step 6: 载荷变为0，去掉接触。
如果用不考虑冲击效应和运动的惯性（即缓慢加载），就用静力分析；否则用动态分析，implicit 或 explicit
都可以，但加载时间5s的话，用动态分析可能计算时间很长。
可以重叠建立实体并划分网格.
去掉接触: 在abaqus/CAE的模型树中，可以把接触在某个step中 suppress.
上方实体除了X方向的移动，还可以向下移动一点点，以保证和下方实体相接触。下方实体不动。
问：刚体穿透，aba大侠看过来 ~
模型很简单，就是一个2d的圆筒受压。采用shell单元。麻烦大家看看，多谢了仿真分析,
删除interaction property中的normal behavier,
就是默认的hard接触，不允许穿透。
后处理时比例设为1, 就看不到穿透。
多谢，Aba,和前面一位大侠，问题解决了。都怪我没有注意deformation scale
factor是否为1，误以为穿透了。分享一点
W-Tk1^-VpSimWe仿真论坛www.simwe.com7V$z@$Sc8ax
刚体穿透解决办法SimWe仿真论坛:y6ZMUwO Y.z
;bX+d^,]v,s1.边界条件和模型,材料是首先需要检查的地方,通产都是最容易发生问题的
A:Q?4cK@nh0f%coSimWe仿真论坛2.尽量使用analytical rigid,而不是 discrete
oQ'f;av|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent3.检查刚体和变形体的法向方向是否正确|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|FluentD5W)y#q4W9IR
4.网格划分是需要fine meshed,从面要比主面细& &
& & D;h.f%oB$j0?6qQ
5.接触时主从面务必设置好,从面尽量使用bode-based element 而不是element-based
K.l3J仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMyccNB
[ 本帖最后由 impactboy 于
问：做得是三维接触分析,主要想知道接触区接触力的具体数值,以便进行下一步单个零件的受力分析.
Step模块，设置History Output，把Domain设为Interaction，选择接触对，选中CFN。
用 *CONTACT PRINT 可以将接触力CFN输出至DAT文件：
*CONTACT PRINT, SLAVE=从面名称
问：我用二维模型，椭圆物压平面物，计算得到了几个接触部分的history
output结果，发现跟实际情况数量级相差太远，一个是由接触压力得到的综合接触力（CFTM), 接触压力（cpress),
接触面积变化（CAREA)，而且接触面积变化不连续。再就是field output
中S子目录下pressure是什么压力，如果取接触面一点看的话，和接触压力是什么关系。
我主要是想得到接触体之间的相互作用合力大小，接触压力，单位为mm.
xile1981答：
你要得到相互作用力大小的话看CFN或者CFT(带摩擦的)就可以,如果跟你实际情况差的太远的话可能是你的模型有问题,检查一下接触的定义吧.
s子目录下的pressure是(equivalent pressure stress),
和接触压力(cpress)并不是一个概念,具体介绍可见手册.
答：法向接触力要看CFN(注意根据方向来选择CFN1,CFN2,CFN3),
你说的CFTM是接触压力和摩擦应力的合力，CPRESS是单位面积上的法向接触力（单位是 力/面积,
和接触力的数量级当然不一样）}
您的浏览器Javascript被禁用，需开启后体验完整功能，
享专业文档下载特权
&赠共享文档下载特权
&10W篇文档免费专享
&每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
数据库实验用表
&&数据库数据
阅读已结束，下载本文需要
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
加入VIP
还剩1页未读，
定制HR最喜欢的简历
你可能喜欢(window.slotbydup=window.slotbydup || []).push({
id: '2014386',
container: s,
size: '234,60',
display: 'inlay-fix'
&&|&&1次下载&&|&&总176页&&|
您的计算机尚未安装Flash，点击安装&
阅读已结束，如需下载到电脑，请使用积分（）
下载：60积分
相关分类推荐
0人评价3页
0人评价3页
0人评价3页
0人评价9页
0人评价9页
所需积分：（友情提示：大部分文档均可免费预览！下载之前请务必先预览阅读，以免误下载造成积分浪费！）
（多个标签用逗号分隔）
文不对题，内容与标题介绍不符
广告内容或内容过于简单
文档乱码或无法正常显示
文档内容侵权
已存在相同文档
不属于经济管理类文档
源文档损坏或加密
若此文档涉嫌侵害了您的权利，请参照说明。
我要评价：
价格：60积分VIP价：Tiddo de Ruiter
SATURNUSSTRAAT 91
2516 AG DEN HAAG
W www.tiddoderuiter.nl[转载]ABAQUS常见问题汇总&-&2.0版（7-9）
good．弹塑性分析
塑性问题不收敛时，msg文件中的常见现象是
1）&& 出现很多EQUILIBRIUM
ITERATION，且TIME INCREMENT 不断减小；
始终出现&& ***WARNING: THE STRAIN
INCREMENT HAS EXCEEDED FIFTY TIMES THE STRAIN TO
CAUSE& FIRST YIELD AT 1 POINTS
&***WARNING: THE STRAIN INCREMENT IS SO LARGE
THAT THE PROGRAM WILL NOT ATTEMPT& THE PLASTICITY
CALCULATION AT 1 POINTS
3）&& 在msg文件的结尾显示
***NOTE: THE SOLUTION APPEARS TO BE DIVERGING. CONVERGENCE IS
JUDGED UNLIKELY.
***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT
解决方法见下面回帖实例。
问：aba_aba老大，你说仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent+j4d&Cd-
用C3D15单元，是15-node quadratic triangular prism,
这还是一种二阶单元，不能用于塑性材料和接触问题。|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent(X&@(p/G$?:J0w
这个是什么原因？www.simwe.com)Te)E#|*@M
-- Getting Started with ABAQUS 12.4.4 Element selection:
When selecting elements for contact analyses in ABAQUS/Standard,
it is better, in general, to use first-order elements for those
parts of a model that will form a slave surface. Second-order
elements can sometimes cause problems in contact simulations
because of the way these elements calculate consistent nodal loads
for a constant pressure.
-- Getting Started with
ABAQUS&&&&
10.3 Selecting elements for elastic-plastic problems:
The fully integrated, second-order, solid elements available in
ABAQUS/Standard are very susceptible to volumetric locking when
modeling incompressible material behavior and, therefore, should
not be used in elastic-plastic simulations.
The second-order reduced-integration elements in ABAQUS/Standard
should be used with caution if the strains exceed 20&40% because at
this magnitude they can suffer from volumetric locking.
问：在我的模拟一个桶土接触问题中，我想加一个和坐标轴成一定角度的集中力，但是直接加集中力会导致不收敛，因此我改用加位移荷载的方法，求出反力，再反推合力，以及合力与坐标轴的夹角，请问这个合力方向是不是与所加的合位移方向相同，本人认为是相同，不知道大家有什其他的看法，请教一个？
有些情况下一样，有些情况下不一样，你试试就知道一样不一样了。施加位移荷载时的反力就等于力载荷。
“直接集中力会导致不收敛”：
一是加载之前要有一个只施加很小载荷的分析步，使各个接触关系都建立起来。
二是如果受集中力的材料是塑性，会因为局部变形过大而不收敛，应该将力施加在一个小区域上，让多个节点共同受力。
当设定stress 和 strain时总不能计算出结果
只有elastic的时候就没问题
当设定时注意什么
要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ. 见附件。
载荷不能太大，改为-1500。
问：二维模型为鋼板接頭實驗，因为是对称的，所以只取其中一半。上、下圓為螺栓（設定為rigid body），Step
1：先讓上；下螺栓接觸鋼板，並施加一微小壓力避免接觸時發生表面震顫。Step
2：下方螺栓固定不動，上方螺栓向上位移直到位移結束。
~7Z-S1_V dwww.simwe.com
qv‑V(y*k8E8oF仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM我的問題是：|Simwe.com|仿真|设计|有限元|虚拟仪器3L%Tfn.L[1]E:T/S
(1)我下方圓孔接觸的地方附近再割密一點，就會發生”too many attemps made for this
increment”，是不是指元素扭曲過大，收斂困難???（因為我把網格切大一點又就可以跑，我有試過aba-aba在別人的回文中的方法，將菜单other
/ Gerneral solution controls / Edit,把Is由默认的12改为30，可結果還是一樣。
t7p仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
y9Ae)bzPQ8B仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM(2)另外，因為我感興趣的部分是在位移過程中（step
2）我想知道鋼板所受的外力，因此我在上方螺栓作reaction force的輸出（REF-1），但是我發現分析出reaction
force來的結果其前面那一部份還蠻符合，但到了後面塑性的部分，力就開始發生不規則的震盪，不知道是哪裡有錯，仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1?
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMe&h5E
2裡，為什麼在一開始位移時卻沒有受力，是因為mesh不夠密嗎（平滑）？，導致一開始並沒有完全接觸到，所以沒受力。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM_I?
|q|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent針對鋼板下方螺栓圓孔的上部頂點位移與上方螺栓的反力，畫出load-displacement關係圖，與真實實驗結果作比較，供大家參考。希望知道問題的人，能不吝嗇指教我。SimWe仿真论坛l`ze&f‑U
www.simwe.comFJ3spY
謝謝了！！
附件: steel-plate.rar ( 20:35, 34.91 K)
msg文件中有
&***WARNING: SOLVER PROBLEM. NUMERICAL
SINGULARITY WHEN PROCESSING NODE STEEL PLATE-1.52 D.O.F. 2 RATIO =
意思是STEEL PLATE-1在2方向上缺乏约束。比如应固定STEEL PLATE-1的底边。
螺栓位移为17.7125，相当于圆孔的半径，圆孔处变形相当大，大变形处的网格需要相当细化，以避免大变形造成畸形单元。细化网格后可能load-displacement關係圖就不震荡了。
螺栓单元也太少，导致主面上有尖角，不容易收敛。每90度圆弧上应至少有10个单元，最好用analytical rigid.
要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ. 见
http://www.simwe.com/forum/viewthread.php?tid=717593
是塑性材料曲线的问题：
- 曲线开始的部分应力上下变化，使ABAQUS很难找到收敛的位置, 应该为一条光滑的曲线；
- 给出的塑性数据点太多，好像也会造成收敛困难，建议改为
&343.964,&&&&&&&
&391.794, 0.0178509
&409.007, 0.0247034
&423.079, 0.0308124
&437.593, 0.0373512
&450.328, 0.0411674
&470.377, 0.0512111
&490.794, 0.0612664
&507.881, 0.0714607
&522.588, 0.0812225
&538.149, 0.0909967
&550.917,& 0.100796
&595.463,& 0.138654
&628.621,& 0.174038
& 654.65,& 0.206945
&675.526,& 0.237475
&695.278,& 0.270728
&706.411,& 0.304198
&700.732,& 0.320962
&643.526,& 0.355134
& 571.81,& 0.374608
& 492.35,&&
470.,&&&&&&&
这样就收敛了（没有最后一行470, 1 ，也可以收敛 ）。在结果中可以看到：
- 大变形处单元扭曲非常厉害，应该在此区域 partition, 使用更加细化的网格。如果有ABAQUS
6.6，可以做自适应网格重划。其他应变小的区域，网格可以粗一些。
- PEEQ最大值为1.719，高于*Plastic最后一行中的塑性应变。对于这样大的变形，难以得到精确的塑性材料曲线：
拉伸试验到一定程度就断了，而且颈缩后的实验结果就开始不精确；压缩试验可以压到很大的应变，但这时摩擦力的影响变大，试样会出现鼓形，结果也不精确。好像没什么太好的办法。
-&&一般的大变形可以用ABAQUS/Standard,
如果想模拟材料破坏就要用ABAQUS/Explicit，见ABAQUS Analysis User’s Manual
第11.6节Progressive damage and failure。
- 弹簧刚度1000太大了。钢板底边的最大y向位移是12，所以弹簧分担的载荷是12*,
而支反力最大值为52070，与弹簧分担的载荷在同一个数量级上。弹簧的作用只是帮助收敛，不应该影响各个部件的受力和变形。可以试试减小弹簧刚度，例如刚度为1，不收敛的话再增大一些。
问：问题 自适应网格重划
之前我用二维模拟钢板接头实验，但发现模拟结果不好，所以改用3D模擬。上、下圆为螺栓（设定为rigid body），Step
1：先让上；下螺栓接触钢板，并施加一微小压力避免接触时发生表面震颤。Step
2：下方螺栓固定不动，上方螺栓向上位移直到位移结束。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMNR
B i&p5oza"H#A
5W~eF2i5a,h-k|Simwe.com|仿真|设计|有限元|虚拟仪器问题是：|Simwe.com|仿真|设计|有限元|虚拟仪器4T.h!o.?Z)]+H"VM:_
(1)还是老问题，我在上方螺栓作reaction force的输出（REF-1），但是我发现分析出reaction
force来的结果其前面那一部份还蛮符合，但到了后面塑性的部分，力发生不规则的震荡SimWe仿真论坛S
Q‑K/C9M%r&M0bT
V4yz&I:b/P"`r(2)我的是ABAQUS6.5，因为下方变形很大，为了避免元素扭曲过大，所以我用自适应网格重化下方大变形螺栓处，但却一开始就发生Abored，不知为什么？
t&J]仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
2D模型的结果怎么不好？你的应力应变在厚度方向几乎没有变化，所以2D和3D的结果应该是一样的，而3D的计算时间长得多。
ABAQUS 6.6才能真正实现自适应网格重划，ABAQUS
6.5的自适应网格主要用于Explicit，以及Standard分析中的表面磨损过程模拟，在ABAQUS
6.5一般的Standard分析中尽管也可以设定自适应网格，但不会起到明显的作用，只是大大增加计算时间。
下列警告都是非线性问题迭代过程中的正常现象，是ABAQUS正在尝试找到正确的解：
&***WARNING: THE SYSTEM MATRIX HAS 8 NEGATIVE
EIGENVALUES.
&***WARNING: THE STRAIN INCREMENT HAS EXCEEDED
FIFTY TIMES THE STRAIN TO CAUSE
&&&&&&&&&&&&
FIRST YIELD AT 34 POINTS
&***WARNING: EXCESSIVE DISTORTION AT A TOTAL OF
2 INTEGRATION POINTS IN SOLID
&&&&&&&&&&&&
(CONTINUUM) ELEMENTS
&***NOTE: ELEMENTS ARE DISTORTING EXCESSIVELY.
CONVERGENCE IS JUDGED UNLIKELY
（以当前的increment不能收敛，自动减小increment，重新迭代）.
在后处理时可以看到大变形而严重扭曲的单元，应在这些地方进行网格细化。
在你的实验中，钢板能够发生这么大的应变而不拉断吗？如果出现了材料破坏，用ABAQUS/standard得到的这种大变形结果就会和实际相差比较大。
螺栓可以用analytical rigid, 以缩短计算时间。
你说的“力发生不规则的震荡”，是说RF出现了下降/上升的变化吗？这好像是正常的，所谓
“局部卸载”。你的塑性材料应力应变曲线有下降段（从数据点
706.445,&& 0.30186
开始），就是说应变达到一定程度时，应力会下降，自然RF也就下降了。
你的试验和FEM已经吻合得不错了，不可能完全一致的。
ABAQUS/Standard不能模拟材料破坏，要用explicit, 见 ABAQUS Analysis User’s
Manual 第11.6节Progressive damage and failure。
ABAQUS 6.6 的自适应网格很不错，可以用于standard分析.
RF震荡，可以用排除法来找原因：
- 怀疑是摩擦造成的： 去掉摩擦，看是否还震荡；
- 怀疑是网格造成的： 螺栓用analytical
rigid，钢材大变形部位用非常非常细的网格，看是否还震荡（你现在大变形部位的单元发生严重扭曲，肯定是太粗了，得到的结果肯定很差）；
-&&&&&&&&&
怀疑是材料塑性曲线的下降段造成的：将材料曲线的下降段改为单调上升，看是否还震荡。
在msg文件中看到反复出现
&&&&&&&&&&&&&&
SEVERE DISCONTINUITY ITERATION&&
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
5 OPENINGS.
&&&&&&&&&&&&&&
SEVERE DISCONTINUITY ITERATION&&
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
0 OPENINGS.
即ABAQUS无法找到正确的接触状态。这是因为x,y方向没有边界条件，只有几个软弹簧，板还是会出现难以确定的位移。
应该尽量让每个部件在每个自由度上都有足够的边界条件：
- x方向： 你的模型是对称的，可以只取1/2，约束x方向位移。
- y方向： 板的顶部(大螺栓的一端）变形很小，可以约束此处的y方向位移，让大螺栓不动，小螺栓向下移动。
increment可以设为
0.001, 1., 1e-8, 1.
请教各位大侠这是怎么回事：
Solver problem. Zero pivot when processing D.O.F. 1 of 1 nodes. The
nodes have been identified in node set
WarnNodeSolvProbZeroPiv_1_1_83_1_2.
OVERCONSTRAINT CHECKS: The model is overconstrained in this
increment due to the interactions among nodes in node set
WarnNodeOverconZeroPivotStep1Inc83Iter1Cut1.
我没有看到OVERCONSTRAINT 信息。如果只是dat文件中有OVERCONSTRAINT
信息，msg文件中没有，就没问题。
刚塑性材料不容易收敛，见Getting Started with ABAQUS 10.4.4 Adding hardening to
the material model。另外要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ.
*Material, name=Material-1
10.0, 100.0
*Material, name=Material-2
10.0, 100.0
把缩放系数设为1，可以看到单元扭曲太大，所以无法收敛。要在大变形部位细化网格，变形小的部位的网格可以很粗。
&&&&&&&&&&&&&&&
问：桩土接触（群桩） 塑性为什么计算不出结果？（弹性可以）
本人一直在研究群桩的计算&&今天计算了一个七根桩的模型（三维）。www.simwe.com#Al%o
模型具体为：
&i8q5z3f[1]L&d3{ @n2}
q!{仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&&&桩：由圆桩（直径2.8m）等效为八面体桩&&桩长100m
全部入土& &梅花形排列www.simwe.com$_
a@**x3Y0r'C
&&&土：共七层土体&&有五层和桩接触&&土体40*40*200m
SimWe仿真论坛
&&&接触：每根桩和五层土体接触（共7*5＝35个接触对
node based surface），桩底采用tie约束&
&共27000左右个单元（c3d8r）
BA@ q-QN)d?F5b^&
&&&我分别用D-P模型和M-C模型都无法算出结果，后来我改为弹性可以算了而且结果还可以
并且平衡了地应力 I|
G|(x7& & &
&& & 讲了这么多
我就是想知道
为什么材料为塑性就算不出结果（经常是单元变形过大&&材料点发生屈服）
是我的塑性参数取的不合理吗？还是别的什么原因？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/_K7p
&& &一般土体yield
stress取多少？第一个step 平衡地应力可以用塑性材料吗？如果用弹性材料 后面计算怎么该为塑性材料？请各位给点建议
我是把yield stress
该大了&&一直不清楚这个值到底该怎么取&&m－c取为c的话
，土体也就几十Kpa&&是不是太小了？ 有没有人知道具体应该怎么取
？ 我200m的土&&如果yield
stress为30Kpa那么在自重作用下不是已经屈服了吗？
问：只在四节点施加一竖向位移，而水平方向怎么会出现应变，且同竖向的应变在同一量级？把弹性单元编为abaqus的uel子程序计算结果也不一致
边界条件对吗？三个节点的位移都固定死了，无法保持体积不变。水平方向应变是因为有泊松比。参见
http://www.simwe.com/forum/post/view?bid=31&id=717759&tpg=1&ppg=1&sty=1&age=0#717759
不同软件要比较单元积分点上的应力，节点应力是插值得到的，不准确。
先做普通弹性材料，正常了再做塑性。
另外，dat文件中有warning:
&***WARNING: A DILATION ANGLE OF 20.000 MAY
RESULT IN DECREASING PLASTIC
WORK&&&&&&&&&&&&&
OR UNSTABLE MATERIAL BEHAVIOR AT HIGH CONFINING STRESS STATES.
SET&&&&&&&&&&&&&
THE DILATION ANGLE LESS THAN 17.831 TO ENSURE STABLE
MATERIAL&&&&&&&&&&&&&
BEHAVIOR UNDER ALL LOADING CONDITIONS
问：在塑性加工中，不同的材料用如下变形抗力模型表示。请问：这个公式指的是整个弹塑性过程还是忽略弹性的？如果包括弹性，那么弹性模量又与该公式有何关系？屈服极限应力在该式中有体现吗？该公式反应的是真实应力应变关系吗？它的适用范围是什么？
不胜感谢！！！！！！
wy_anli答：塑性加工中的具体材料模型有成千上万个，大类可以分为理想刚塑性（忽略弹性、忽略硬化，主要应用于弹性变形相对于塑性变形较小，可以忽略的情况下，如：锻造等）、理想弹塑性（忽略硬化）、弹塑性（在材料成形时其弹性变形不可忽略，例如板料成形）、等等。
你的这个模型属于弹塑性模型，没有忽略弹性。
屈服极限是材料发生塑性变形时的临界应力，是材料的本身性质，在本公式中不能体现。
纵向应力与纵向应变的比例常数就是材料的弹性模量，也叫杨氏模量。也不能在此公式中反映出来。
任何材料模型都不可能完全反映真实的应力应变关系，当然次模型也不例外，只是最大近似。
一般只有根据材料提出材料模型，其适应范围也只能适用于特定材料以及特定场合。
wpengw123 wrote:
非常感谢。还是有些不解，就直接问了，不吝赐教！
那么该公式是否就是整个弹塑性应力应变关系的简化式呢？还是只描述塑性变形阶段？
是塑性应力应变关系的简化式，但还缺少抗拉强度或极限应变（当材料达到这个临界状态时就会发生断裂，继续算下去，就没有了意义）。弹性阶段在塑性阶段前考虑。
wpengw123 wrote:
如果包含弹性，那么杨氏模量是不是和公式有冲突呢？
杨式模量不会和公式有冲突，杨式模量是材料本身的参数，是材料在弹性阶段时应力应变线性关系的系数。Abqus中定义材料时，弹性和塑性是分开来的。弹性时只有杨式模量和波松比两个参数。
wpengw123 wrote:
而且它似乎常用于金属板料成形，优点是什么呢？
板料成形要考虑弹性变形，因为板料成形是大变形小应变情况，弹性应变不能忽略不计，并且在塑性成型结束后弹性变形要完全恢复（也就是回弹）。相比较来说，例如锻造等成形，弹性应变相对于塑性应变非常小，可以忽略不计，也就是要采用刚塑性材料模型。
wpengw123 wrote:
式中常数和变量的具体意义是什么呢？
参见塑性成形原理。
wpengw123 wrote:
怎样和单向拉深的实验数据联系，又怎样在Abaqus中使用呢？
单向拉深的实验数据通常是名义应力和名义应变。要通过名义值和真实值之间的关系式将其转化为真实值也就是等效应力和等效应变。Abaqus中的材料采用真实应力和真实应变。
wpengw123 wrote:
看到很多文章在模拟板料塑性成形时用到此模型，但是在给出杨氏模量及强度系数和硬化指数的取值时却都没有给出屈服极限，屈服极限不需要吗？
屈服极限也就是材料在屈服点上的应力应变状况。材料模型中应变为0时的应力就是屈服极限，大多数金属的屈服应力为材料弹性模量的0.05%-0.1%。
qyuan_13: &关于对球墨铸铁的强度评价标准
通常情况下，铸铁的脆性比较强，延伸率很低，一般以延伸率为0.2%作为屈服极限。对于脆性材料作强度分析的时候，采用主应力来评价；但对于球墨铸铁，有比较高的延伸率，比如QT450-5、QT450-10，其延伸率分别为5%和10%，那么，对这样的材料作强度分析的时候，是不是应该采用Mises应力来评价？按说10%的延伸率比较高了，合金钢20Cr的延伸率也为10%，40Cr的仅为9%（刘鸿文的材料力学）。那么，在脆性材料和塑性材料之间的界限到底是多少呢？不能仅仅认为只要是铸铁，就得用主应力来评价吧！
GZSI2x2yiSimWe仿真论坛
问：在ABAQUS中,对于理想弹塑性材料,在计算时,ABAQUS中,材料屈服有很多准则.怎么在ABAQUS中定义我所要采用的Mises屈服准则?
在帮助里搜索一下就可以找到： ABAQUS Analysis User's Manual 11.2.1 Classical
metal plasticity：
The Mises yield surface is used to define isotropic yielding. It is
defined by giving the value of the uniaxial yield stress as a
function of uniaxial equivalent plastic strain, temperature, and/or
field variables. In ABAQUS/Standard the yield stress can
alternatively be defined in user subroutine UHARD.
Input File Usage: *PLASTIC
ABAQUS/CAE Usage: Property module: material editor: Mechanical /
Plasticity / Plastic
理想弹塑性容易有收敛问题，最好加入硬化特性，见Getting Started with
ABAQUS&& 10.4.4 Adding hardening
to the material model。
你的弹性模量是不是太小了，以钢为例，如果单位是N, mm，MPa，则弹性模量是210000, 屈服应力是420。
判断材料是否屈服
问：DOCUMENT上说PEEQ是等效塑性应变（Equivalent plastic
strain),PEMAG是塑性应变值（Plastic strain
magintude),感觉前者主要是判断材料是否屈服，而后者呢？？
在比例加载时（即加载过程中主应力方向和比值不变），大多数材料的PEMAG和PEEQ相等。这两个量的区别在于，PEMAG描述的是变形过程中某一时刻的塑性应变，与加载历史无关，而PEEQ是整个变形过程中塑性应变的累积结果。例如，单向拉伸一个圆柱体，使其发生塑性应变，再通过单向压缩使其恢复初始长度，则最终的PEMAG为0，而PEEQ是拉伸和压缩过程中塑性应变之和。
问：AC YIELD是材料进入塑性的标志，1是发生了塑性变形，0是没有，这样理解对么？
PE是塑性应变，可为什么AC YIELD为零的区域（没有发生塑性变形）也会有PE值呢？？
积分点（而不是节点）上的 AC YIELD 和 PEEQ (而不是PE)是完全吻合的： PEEQ
& 0 则AC YIELD = 1, 发生了塑性变形; PEEQ = 0 则AC YIELD = 0,
没有发生塑性变形。
云纹图是经过插值平均的，看上去AC YIELD 和 PEEQ
会显得不一致。如果逐个察看积分点上的结果，就可以看出它们是吻合的。
PEEQ 和 PE 的关系： PEEQ 类似于mises应力，PE类似于各个应力分量。
问：比如计算地基的极限承载力，m-c理想弹塑性，用计算结果不收敛作为破坏评判准则。
现在我做的时候，荷载是手动调整，每次调整完再提交计算，感觉比较麻烦
能否实现这样的功能：给出初始荷载以及增量步长，要是能自动调整就更好了。每一级计算完成后保存结果，并增加荷载继续计算，直至不收敛
ansys好像有循环语句可以实现，不知道在abaqus中如何实现。谢谢了~
不需要特别的设置，直接创建一个step, 设置initial increment（例如0.1）,
在这个step里施加一个载荷（比如1000N）, 于是step time = 0.1时，载荷就是100N; step time =
0.2时，载荷就是200N；超过极限载荷时，计算不收敛，ABAQUS会自动减小increment。
aba-aba,你有没有碰到过这种问题,在理想塑性状态(切线模量设为零)下,算出来的结果最大应力超过屈服应力?而并不是产生不收敛?
席勒同志好，关于这个问题，我的理解是这个样子地：
如果一个单元所受的载荷是单调增加的，则其积分点上的Mises应力和等效塑性应变PEEQ与塑性材料曲线中的真实应力和塑性应变相吻合。注意是积分点，而不是外插后得到的节点应力。
另外多说两句应变方面的问题，PEEQ描述的是整个变形过程中塑性应变的累积，，无论是拉应变还是压应变，都取绝对值，然后加在一起。如果分析过程中出现了卸载，则Mises应力和PEEQ与塑性材料数据中的真实应力和塑性应变不再吻合。这里所说的“卸载”不一定是模型外载荷的减小，在加载的过程中，摩擦或模型各部分的刚度差异有可能会改变应力场的分布，从而造成模型局部的卸载。
PEEQ和PEMAG的具体含义见
http://www.simwe.com/forum/post/view?bid=31&id=691447&tpg=1&ppg=1&sty=1#692219
无论有没有卸载，积分点上的Mises应力都不会超过理想塑性材料的屈服极限。卸载只是对PEEQ有影响。Getting
Started with ABAQUS&& 10.4.6
Postprocessing the results 有解释：
"Do the values listed in the contour legend surprise you? The
maximum stress is greater than 580 MPa, which should not be
possible since the material was assumed to be perfectly plastic at
this stress magnitude. "
"The Mises stress at an integration point can never exceed the
current yield stress of the element' however, the
extrapolated nodal values reported in a contour plot may do so. In
addition, the individual stress components may have magnitudes that
exceed the value of the
only the Mises stress
is required to have a magnitude less than or equal to the value of
the current yield stress."
另外，帮助中的这个例子还说了，理想塑性不容易收敛，最好加入硬化特性，让材料曲线中最大的应力应变大于模型中可能出现的应力应变。
作切削时候遇见这种情况怎么办？？？
刀具最后有点浸入工件里了|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent0@6d
但是还在切削 不知道怎么回事
模型做得真不错。工件的节点没有进入刀具，只是由于单元变形太大，所以单元位于刀具里面，可以试试：
- 尝试其他的接触算法，比如kinematic contact method, Balanced master-slave
contact (即把weight factor设为0.5). 见Getting Started with
ABAQUS& 12.7 Defining contact in
ABAQUS/Explicit，ABAQUS Analysis User's Manual&
21.4.4 Contact formulation for ABAQUS/Explicit contact pairs。
- 裂开的地方网格细化，比如工件沿刀尖圆弧分布10个单元。
- 刀具用analytical rigid，如果用discrete rigid,
刀尖圆弧处也要增加单元（比如10个）。见Getting Started with
ABAQUS&&&&
12.8 Modeling considerations in ABAQUS/Explicit。
- 减小刀尖圆弧的半径。
问：关于切削加工仿真的，
你的模型里网格没划分，接触对Int-1也没定义，根本就运行不了。
接触面不能只是刀头的一个点吧，应该是刀的两条折线。explicit接触问题参见Getting Started with
ABAQUS&&&&
12.7，12.8
想模拟材料破坏要定义失效判据，并让切下来的单元脱离模型，见ABAQUS Analysis User’s Manual
第11.6节Progressive damage and
failure。切削问题很复杂，建议你先做一些简单的explicit接触问题。
问：真诚求助高手指点－切削 网格变形
abaqus 6.5& 的 adaptive mesh 并不会细化网格，见
http://www.simwe.com/forum/thread--1.html
可以不用adaptive mesh ，自己把网格划细些，像
http://www.simwe.com/forum/thread--1.html 一样。 你的网格太粗了。
可以先把坯料长度改短些，以缩短计算时间，调试成功后再改为真实的坯料长度。
刀尖上要有圆角。关于接触的定义见&
http://www.simwe.com/forum/thread--2.html
ABAQUS常见问题汇总（1.0版）
Explicit中的接触问题
材料非线性实验与FEM的比较
**请教**材料非线性实验与模拟的差异问题!
我做了一些复合材料的层间动态压缩实验和模拟，实验测得的应力应变曲线表现出明显的非线性，见附图，按这条应力应变曲线作为材料参数建立有限元模型后，计算出的应力应变曲线和实验曲线对比，在失效应力前基本吻合还算可以，可是超过失效应力后，实验测得的仍是非线性的，应力下降应变继续增加，可是模拟出来的结果却是突然卸载，想请问一下为什么？是不是ABAQUS软件本身认为只要出现应力下降就是完全失效然后不再承载呀？看到的一些材料模型有理想弹塑性、弹性+塑性强化段，对于这种应力应变曲线模拟不出来吗？
iy1R‑Yw|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
5f9h2U1MD仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent另外问一下，材料弹性模量是不是不能太小呀，因为我没测准它的弹性模量，就想先选一下大概的值算一下，取在100GPa左右都可以算下去，取到20-30GPa的话就出错，原因是变形速率是波速的一点几倍，然后就出错算不下去了，网络分的感觉已经足够细了，还是有问题，在这也一块请教一下大家。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/Iq'?+B"Up
8Qjw)~6v0{D+ql
谢谢大家参于讨论！
估计是这样：
试验曲线是整个板的结果，板的一部分失效了，另一部分还没实效，能继续承载，所以应力不会一下子降下来；
你图中的FEM结果是一个单元的，这个单元失效后就一点也不能承载了，它的应力一下子降下来。
可以比较一下试验和FEM的载荷或支反力结果，它们都是基于整个板的。
支反力RF是“支承处的反作用力”，有边界条件的地方才有RF. 如果边界条件定义在整个面上, 可以用distributing
coupling把加载面和一个参考点耦合起来，然后把边界条件加在参考点上，计算后看参考点的RF。
透射杆在子弹射入方向（假定是方向2）上是固定的吧，施加了此边界条件的节点上的RF2之和，就是试样受到的载荷，对称轴上的RF没用。不用coupling也行。
接触问题技巧整理
(cafe) 【分享】ABAQUS
培训资料 － 接触单元部分－Lecture1
如何加入传统的接触单元，如：goodman单元、desai单元？
答：不像ANSYS，
ABAQUS不须定义接触元。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM5Hw2cnk‑q*t(@[1]{u
用*SURFACE在3D单元的面上定义SURFACE，然后用*CONTACT PAIR定义哪个面或哪个面接触。
请教aba_aba，如何设置库仑摩擦接触面属性，不知道在什么地方设置该接触面切向摩擦属性中的内聚力c值？
Getting Started with ABAQUS&
&附录 C.9.1 Creating an interaction property
其中用的是无摩擦的接触属性。库仑摩擦则需把 Mechanical / Tangential Behavior
中的frictionless改为penalty, 输入摩擦系数 friction coeff .
内聚力不知道，没用过。
[讨论] 关于接触分析中的收敛问题。
我做了一个圆柱压一个平面的接触分析，采用了两种方案
现在的问题，我发现在无摩擦计算中， 最后一步的残余位移和力都非常小： LARGEST RESIDUAL
&& & 3.245E-09，
LARGEST CORRECTION TO DISP.&
-4.941E-15www.simwe.com*x&xL[1]^7M.M]L&
www.simwe.com9T‑K4L
而在有摩擦的计算中，最后一步的残余位移和力都很大啊：仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMi#DB^L
LARGEST RESIDUAL FORCE&
|Simwe.com|仿真|设计|有限元|虚拟仪器,mKU
L.N?&XA2G$V‑N
LARGEST CORRECTION TO DISP.&
&&&3.786E-04www.simwe.comJ
t;a*lchU(W)w6`
请问这个是不是表明在有摩擦的计算中，目前还没有收敛啊？
收敛准则用默认的就可以。根据Getting Started With
ABAQUS&&8.2.2
节&&“ABAQUS/Standard中的平衡迭代和收敛”：
“CRITERION FOR RESIDUAL
FORCE&&FOR A NONLINEAR
PROBLEM：&&5.000E-03“ 的含义是
RESIDUAL FORCE 和 TIME AVG. FORCE 的比值小于 5.000E-03
即可，而不是&&RESIDUAL FORCE 的绝对值小于
5.000E-03 。
你的计算中 TIME AVG. FORCE =
2.302E+04，&&RESIDUAL FORCE =
做接触分析时，应观察msg文件中的以下信息
（接触状态的变化：1个节点由不接触变为接触，10个节点由接触变为不接触）
CONTACT CHANGE
SUMMARY:&&&&
1 CLOSURES&&&
10 OPENINGS.
SEVERE DISCONTINUITY
ITERATION&& 2 ENDS
CONTACT CHANGE
SUMMARY:&&&&
CLOSURES&&&&
4 OPENINGS.
如果分析能够收敛，CLOSURES和OPENINGS的数目会逐渐减少，否则一般是接触定义不当，可以考虑以下解决方法：
●& &
检查所定义的接触面、接触参数和边界条件是否正确；
●& &
在静力分析中，必须定义足够的约束条件，以消除各个平移和转动自由度上的刚体位移；
●& &
避免过约束（overconstraint）；
●& & 合理地定义接触面、接触参数和过盈配合；
●& & 使用足够细化的网格；
●& &
不要在接触面上使用C3D20、C3D20R和C3D10等单元；
●& &
在接触对上设置微小的过盈量，以保证在分析的一开始就已经建立起接触关系；
●& &
施加临时边界条件，以保证在接触关系建立之前，模型也不会出现刚体位移；
●& &
将分析过程分解为多个分析步来完成，让各个载荷分别在不同的分析步中逐步施加到模型上，避免使接触状态发生剧烈的改变。
下面按照应用领域对回帖分了类，但其实定义接触的原则是相同的，例如岩土类问题的接触问题的解决方法，往往同样适用于金属成型中的接触问题。
问：此类问题,如何在接触分析中约束?请指点!SimWe仿真论坛 U]I,wC%Y
,q$S*g2Yo.O@E
J/{^fT h-j1u-kTO9x
另外,在用壳模拟分析时,abaqus默认的是不是中面?能否更改?即除了材料方向外,如何更改厚度方向?www.simwe.comF$mn‑}&`&y*]Se
''接触分析中约束'': 不明白你的问题。
对称模型，只需取1/4分析，这样X和Z方向就有了足够的约束。
在加载荷之前，增加一个分析步，使part-1-1向上移动0.001，不加载荷，从而使接触关系平稳地建立起来。在下一个step中，去掉此边界条件，加载荷。
壳模拟分析： 见&
http://www.simwe.com/forum/thread--2.html
ABAQUS常见问题汇总（1.0版）
16.2& 用OFFSET参数偏置shell的中性面
面上的网格密度
问：如何得知接触表面网格太密或太疏而造成不收敛的问题呢..仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent-U&
www.simwe.comu3K Sk*Z2T
掌握 master 面网格密度比slave面网格密度疏，此点足够吗？SimWe仿真论坛7OZ}z8Q'[;v6S
[1x`仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
M.V(@%` Z& & &
}J&H+}接触分析中若有大变形情形，因为很不容易收敛
e3_zy1oZSimWe仿真论坛仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM"]{.I2j)y&HAw
改变 increment size 中的 minimum
值，将其值设小（default：1E-005）SimWe仿真论坛7f)YN:Ch&NN*X}
]#?m0l!`:U#S4}& &
c6b仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM能提高分析成功的机率吗？
网格疏密要靠经验，比如接触面在某个方向上只有一两个单元，就是太疏了，如果接触面在某个方向上有上百个单元，就可能是太密了。
大变形时应该减小increment size 中的 minimum
值，但如果接触或塑性材料参数定义得不正确，还是不能收敛。
为什么计算中断
求救：为什么我的模型提交JOB后计算一半就中断，提示错误：too many attempts made for this increment 可能是不收敛，减小最小增量步的值只能再多计算几步，但还是算不完。各位大侠请指教啊！
Tie的两个面之间的距离大于你定义的position tolerance = 0.0001, 所以tie没有完全起作用，骨头可以乱动，ABAQUS找不到合适的接触状态。
用display group来单独显示surface，例如 c4-1."Surf-c4 i" 和 i45-1."Surf-i45 s"，可以看到它们的距离大于0.0001。
在后处理中用view cut, 也可以查看各个面之间的距离。
有两种解决办法：
- 把 i45和i56 改薄些；
- 移动红烧排骨之间的相互位置，但这同时会改变接触面之间的距离。
Tie的 position tolerance 应该定义得大些，例如一个单元的大小，同时不要选Adjust slave surface initial position, 以免调整节点位置造成畸形单元。
接触面定义得太大，也会造成收敛困难。应该只把可能接触的位置定义为接触面。
c5-1."Surf-c5 p" 既是主面，又是从面，可能会有问题。应该分别定义两个面。
C3D4单元精度很差，应该用二阶单元 C3D10M。接触问题不能用二阶HEX单元，但可以用C3D10M。
在step 1中不要施加实际的载荷，而是施加一个3方向的微小的力，使接触关系建立起来，在step 2中再施加实际的载荷。
接触部位的网格划细些，否则只有几个单元发生接触，就不容易收敛。在后处理中察看接触面之间的距离（field
output的COPEN），可以看到有几个单元发生接触了。
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
如果位移远远小于一个单元的长度，可以用nlgeom = off.
如果还不收敛，可以增大 “允许的最大迭代次数”，见
http://www.simwe.com/forum/viewthread.php?tid=731009
http://www.simwe.com/forum/viewthread.php?tid=693237
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
关于接触产生节点应力集中的问题...请大家帮帮忙~拜托
我使用一只脚的模型
+zKLu& & &
v‑ywww.simwe.com使用静态分析方式
F&r!hSimWe仿真论坛将脚模型固定,地面向上撞击脚底SimWe仿真论坛Q'EGp
将地面设为master surface,而脚底设为slave
surface(网格较细).SimWe仿真论坛3z#|j& &
脚底与地板之间的接触条件是用surface to surface contact
7D,{s4&nV%OSimWe仿真论坛之间摩擦力设为0.11www.simwe.com1%u9K2VE
其它为预设
V;D3j|Simwe.com|仿真|设计|有限元|虚拟仪器可是结果如图
奇怪的应力看来是因为你的脚面本身就是凸凹不平的，凸的地方接触压力大，凹的地方接触压力小。
在后处理时无法正常察看CPRESS,&
说明接触面定义得不对（你的模型是用其它软件建的吧），你在ABAQUS/CAE里重新定义一下contact的salve
surface和master surface.
没必要把整个脚的表面定义为slave surface, 只定义足底可能接触的部分即可。
下图中，正中间的节点只属于下面两个单元，不属于上面的单元，这样的网格精度很差。每个节点都应该属于周围所有的单元。
单元类型C3D8的精度很差，应该用C3D8I或者C3D8R。
Step的increment size不要用fixed，这样出现收敛问题时，无法自动减小increment size。
脚后跟处应力很大，如果关心此处的应力，应该细化单元。
足底和地板接触前的应力很小，没必要模拟接触前的过程。可以一开始就让足底和地板接触。
如果想模拟瞬间的撞击，要用explicit作动态分析，
在后处理时，用display group 来只显示所有surface, 查看cpress和copen, 看到可能导致不收敛的原因是：
有3个接触对的master surface定义得太小，变形后接触点将要移出master surface之外 （你加4方向正转矩时能收敛，可以看看此时是不是没有这种问题）。
另外，各个面上都只有一个节点上发生接触，这可能会有问题。应该在此处细化单元，并适当增大adjust的值，使更多单元发生接触。
任务提交后提示以下错误：“Too many attempts made for this increment”
在msg文件中看到
CONTACT PAIR (ASSEMBLY_BLANKBOT,ASSEMBLY_TIE-1_DIEDURF) NODE
BLANK-1.5 IS& OVERCLOSED BY 0.0512228 WHICH IS TOO
这往往是因为接触面的法线方向定义反了。定义刚体和shell的surface时，要注意选择外侧。在后处理中可以显示法线方向：ABAQUS/CAE
Manual&&&&
37.7 Displaying element and surface normals。
问：请问， 过盈配合怎么做， 用哪个菜单？
----------------------------------------
*CONTACT INTERFERENCE
----------------------------------------
此关键词的使用方法为：
*CONTACT PAIR, INTERACTION = &接触属性的名称
&, ADJUST = &位置误差限度
&从面名称 &，&主面名称
*AMPLITUDE，NAME = &幅值曲线的名称 &
*CONTACT INTERFERENCE，AMPLITUDE = &幅值曲线的名称
&从面名称 &，&主面名称
&，&过盈量或间隙量 &
其中，参数&过盈量或间隙量 &为负值表示过盈量，正值表示间隙量。使用
INTERFERENCE类似于施加载荷，不能在initial分析步中对其进行定义，而只能在后续分析步中定义。
☆ ABAQUS/CAE操作： Interaction模块，主菜单Interaction → Create，点击Edit
Interaction对话框底部的Interference fit。
如果在分析结果中看到，使用 *CONTACT
INTERFERENCE所定义的过盈接触没有在模型中起作用，则有可能是以下原因引起的：
● *CONTACT
INTERFERENCE不能使用ABAQUS默认的幅值曲线Ramp（从1降至0），而要使用自定义的幅值曲线，使过盈接触的幅值在整个分析步中从0到1逐渐增大。
● 如果在ABAQUS/CAE模型中两个接触面之间有宽度为t的缝隙，而事实上它们之间应该存在接触关系，则关键词 *CONTACT
PAIR中的参数 ADJUST = &位置误差限度 &
必须略大于此缝隙的宽度t，否则ABAQUS会认为这两个接触面没有接触，*CONTACT INTERFERENCE 中的参数
&过盈量或间隙量 &不会起作用。
● 参数&过盈量或间隙量
&为负值才表示过盈接触，正值表示两个面之间存在缝隙。
问：使用clearance和使用CONTACT INTERFERENCE算出来的结果有什么区别么？
*CLEARANCE也可以，但有一些缺点，它只能用于small
sliding，不能在ABAQUS/CAE中定义，不能在step中修改过盈量，而且没有使用ramp,
也就是说，在分析的一开始全部过盈量就会被施加在模型上，如果此过盈量太大，ABAQUS就可能无法达到收敛，而且即使减小时间增量步也没用。
问：我的两个接触面是经过运动后才到达相互接触的状态的，但这时并不能达到理想的接触（两个面有部分gap,很小但的确没有接触，用contact
print 打印出来slave nodes相关项，都是zero)。我想在这个时候来adjust the contact
gap,让abaqus把接触面间的gap调成没有。
因为一开始时两个接触面的间距非常大，所以不可能在一开始就调整，而且我还需要对那个运动进行分析。
使用 *CONTACT INTERFERENCE可以在分析步中定义接触对之间的距离或过盈量，你可以试试。
使用 *CONTACT INTERFERENCE时要注意：
1) 其默认的幅值Ramp是从1降至0，不能使用，而要使用自定义的幅值曲线，使幅值在整个分析步中从0到1逐渐增大；
2)& 正值是缝隙，负值是过盈量。
问：我所关注的只是螺钉拧入之活的应力分布，而对于拧入的过程则不关心。我想不放把问题说的更细一些。
把螺钉拧入骨头可分为如下几种情况：
A、自攻型，就是硬生生的把螺钉拧入骨头。
B、助攻型，先在骨头上打一个比螺钉的直径小一些的洞，之后再把螺钉拧进去。这样做的好处是螺钉不容易折断。打不同直径的洞，在螺钉拧入后的应力分布应该是不同的。
C、用一个同样的螺丝钉（不锈钢）在骨头中钻个洞，再把需要的螺钉（钛）拧进去。
D、或者有时用一个稍小一点的不锈钢螺钉钻洞，再拧入钛螺钉。这样可以比C更稳固。
E、我不小心把洞钻大了，这样螺钉就容易松动。当然这在临床中是不允许的，但也存在。
各位大侠能否就这几种情况分别讨论一下，我觉得是不是可以用过盈装配来模拟？但具体A、B和D有什么样的区别呢？想不明白。等待着各位大侠解惑~~~
这么说我就明白了，弄明白了就感觉很爽。
一定是兽医，把骨头这么搞来搞去的，让病人求生不得，求死不能，关云长来了也受不了啊。
A（怒射破门）就是过盈配合非常大。
B（后卫助攻）就是过盈配合小一些。
C（换上替补）就是几乎没有过盈配合。
D（影子前锋）就是过盈配合也比较小。
E（自摆乌龙）就是接触面之间存在缝隙。
具体过盈配合该多大，估计该在螺钉上装些传感器，得到试验结果，再调整你的过赢配合量，与实验相吻合。过盈配合的模拟方法见。
zero4:& 用ABAQUS软件处理管土相互作用中的接触面问题
给有兴趣的朋友一个参考。
H%{:M(|&vX
附件: 1.part1.rar ( 13:36, 195.31 K)
问：计算中止，.msg文件里显示
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING
INSTANCE PART-1-1 D.O.F. 8 RATIO = 6.49141E+11
检查约束没有错误，而且节点346位于模型中部，不知道是什么原因引起的？望高手指点。附上文件，太大分了2部分。
桩上没有边界条件，完全靠接触面上的摩擦来消除刚体位移，但接触关系是在迭代过程中慢慢建立起来的，不是一开始就有摩擦力，所以会出现刚体位移。见
http://www.simwe.com/forum/post/view?bid=31&id=680814&tpg=1&ppg=1&sty=1&age=0#680814
接触对的主面上不能有尖角，桩的两侧要各定义一个接触对，底部用tie. 见
http://www.simwe.com/forum/viewthread.php?tid=699614
http://www.simwe.com/forum/viewthread.php?tid=679059
http://www.simwe.com/forum/viewthread.php?tid=684042
http://www.simwe.com/forum/viewthread.php?tid=679059
中谈了过约束的问题：
"msg文件中出现zero
pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact,
出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region,
然后选择区域时不要包含这一圈节点。"
你的桩底部只有一个单元，就没法用上述方法。要在底部细化网格才行。
不必整个模型都用那么细的网格，离接触面较远的地方可以粗一些。
新手请教：The plasticity/creep/connector friction algorithm
did not converge at 12 points
模型做得不错啊，改为弹性材料就可以收敛，说明模型没有太大的错误。我认为需要改进以下方面：
应力应变变化很剧烈的部位（桩底面上方和下方的的soil）,单元不能太狭长。原则是在后处理时不能看到单元之间的应力和PEEQ不连续（颜色发生跳跃）。
桩侧面接触状态变化很大, soil和桩的网格沿圆弧方向和纵向都要再细化一倍。
对称模型取1/4就够了。用C3D8R单元可以缩短计算时间。
rebar是truss单元，不能用*Solid Section,要用
&&*Solid Section, elset=***,
material=REBAR
&&截面面积
桩底部的从面网格还是比主面粗糙,会有问题。桩底部可以不用接触，改为tie，收敛就更容易：
&&*Tie, name=Tie-Bot, position
tolerance=0.1
&&soil-1.pilebot,
mpile-1.bot
我做了个轴对称模型，rebar相当于shell, 没有*initial condition,
见附件。从结果中看到，虽然应变很大，但位移很小，因此可以不用在step中设 NLGEOM=yes，这样收敛就更容易些。
塑性问题要把step 2的初始时间增量步和允许的下限设小些，比如：
&&0.01, 1., 1e-06, 1.
下面的警告信息都是说塑性问题不收敛：
ELEMENT SOIL-1.709 IS DISTORTING EXCESSIVELY.
***WARNING: THE PLASTICITY/CREEP/CONNECTOR FRICTION ALGORITHM DID
NOT CONVERGE AT 5 POINTS
Soil-Axisymm.cae(67.54k)
问：大侠们看看这个初始应力为什么总是施加不成功了？？
建议在step 1中只保留桩底部的接触，去掉其他接触。step 2中再加上所有接触。先做弹性材料，收敛了再做塑性。
问：求教：为什么基坑开挖中围护结构的作用没有体现？（附有inp文件）
简单介绍一下：通过读入初始应力的方法计算地应力平衡；step1中杀死所有的围护结构（包括3个支撑梁Graider1,Graider2,Graider3、一个连续墙Wall）；step2中激活wall、开挖第一步所要开挖的土体kw1、激活Graider1。后面的暂时没有计算。可是计算结果是wall和Graider1好像没有起作用。那位高手能帮我指点一下，先谢谢了!
围护结构和土体之间要定义接触或*rebar或*embeded，否则它们仅仅是空间位置重合，不会自动发生相互作用。
问：隧道开挖与管土相互作用中的接触问题求教
本人做了隧道开挖中的管土相互作用的模型。有些管－土接触问题求教。|Simwe.com|仿真|设计|有限元|虚拟仪器h&o},B0}(r%m:NQ
模型说明：开挖隧道位于管道正下方，由于是研究接触算法，仅做三步开挖，且没有进行支护。土体使用实体单元c3d8r，管道采用三维壳单元s4.模型1管土接触采用法向www.simwe.comkXP%z
硬接触，切向罚函数（文件名exva-h1-fenbu）。模型2采用tie约束管道与土体节点。（文件名tie-ex-fenbu)
&IfZD~'ev2.
计算后出现的问题SimWe仿真论坛vp2W&W8~E
（1）开挖隧道纵向中间部位土体时，两个模型的管道位移基本一致，接触算法作用不大。www.simwe.comY:K6Q
AQ[1]i8z,r
（2）隧道全部开挖的情况下：使用tie约束的模型2计算不收敛，通过查看塑性变形，发现和模型1的差不多，不知道为何不收敛。请问这种隧道开挖引起土体位移对管道的影响问题，设置接触的必要性以及模型2为何不收敛。
#k*B&k,~3Mwww.simwe.com仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMZe
[ 本帖最后由 doudout 于
或许是定义tie的时候应该加上 position tolerance=0.1, 以确保两个面tie在一起。
你的单元是c3d8, 精度很差，应该用c3d8r。
单位对吗？如果长度单位使用mm，那么质量的单位必须是“吨”，相应的密度单位就是吨/mm3，应力的单位是N/mm2（即MPa）。
见 http://www.simwe.com/forum/viewthread.php?tid=684273
问：做基桩的水平静载试验模拟，MSG文件提示迭代增量步太大，不能进行
取1/2原型，近场土C3D20R，远场CIN3D8，对称面上Y向对称约束。请赐教！
答：我认为模型中有以下问题：
1。塑性材料和接触面上都不能用C3D20R和C3D20单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用C3D8I
(在所关心的部位要让单元角度尽量接近90度)，如果只关心应变和位移，可以使用C3D8R,
几何形状复杂时，可以使用C3D10M.
我尝试把你模型中的接触关系以及材料参数中的*Mohr Coulomb 和 *Mohr Coulomb
Hardening都去掉，变为弹性材料, 仍使用C3D20R就收敛了。
2。接触对中的slave surface应该是材料较软，网格较细的面，你正好颠倒了。
接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于接触面之间的距离，否则ABAQUS会认为两个面没有接触：
*Contact Pair, interaction="SOIL PILE SIDE CONTACT", small sliding,
adjust=0.02
在ABAQUS/CAE中可以在edit interaction对话框中设定此参数。
4。定义tie时也应该设定类似的position tolerance:
*Tie, name=ShaftBottom, adjust=yes, position tolerance=0.1
Tie和contact pair都定义在step之前，你现在的模型是对的。
overconstraint警告信息如果只是出现在dat文件中，没有出现在msg文件中，就没问题，说明ABAQUS自动解决了过约束问题。
如果overconstraint警告信息也出现在msg文件中，说明ABAQUS无法自动解决此问题，这时分析往往不会收敛，在后处理时可以用display
group显示出现过约束的node set (ABAQUS会自动为这些节点建立node set),
需要你自己修改模型，避免过约束。
附上修改后的文件，运行时没出现零主元（zero pivot），我做的改动是：
- 去掉所有的 *Mohr Coulomb 和 *Mohr Coulomb Hardening
- 在*Contact Pair中添加 adjust=0.02
- 在*Tie 中添加 position tolerance=0.1
还应该改的是：
- 接触问题不能用C3D20R单元（所以目前施加载荷时很难收敛）。
- 接触对中的slave surface应该是材料较软，网格较细，即应细化SoilInnerSurface面。
msg文件中出现zero
pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact,
出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region,
然后选择区域时不要包含这一圈节点（我附上的文件中没有做这样的修改）。
无穷元我没用过，不太了解。
接触定义在哪个分析步取决于你模型的实际物理背景，如果从一开始两个面就是相接触的，就定义在initial或你的第一个分析步中；如果是后来才开始接触的，就定义在后面的分析步中。边界条件也是这样。
你现在把接触定义在后面的分析步中，应该是可以的。
你各个分析步的设置也没看出有什么问题，只要符合实际物理背景，并且利于收敛即可。
Data5.rar (94.68k)
岩土问题我没搞过，回答不了你，论坛上好像有这方面的资料，你搜一搜。
单元扭曲的模型传上来看看。
不太了解你的单向多循环加载法。原则上，如果是静载荷，用static分析步即可，如果加载时间很短（例如在地震中），或载荷本身的性质是动态的（例如来自旋转机械的载荷），在分析中不能忽略结构的惯性，就要采用动态分析，可以用ABAQUS/standard或explicit.
我在前面上传的文件里用*CONTROL设了允许的迭代次数18，意思是18次迭代不收敛时，才减小时间增量步（ABAQUS默认的值是12）。一般情况下不必设置此参数，如果在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛，就可以用*CONTROL来增大允许的迭代次数。
你的C3D20R单元导致了接触问题的不收敛，所以我用*CONTROL其实是没用的。
你的Test.inp也不能收敛啊。另外，如果你传上来两个模型，一个收敛，一个不收敛，应该写清楚这两个模型有什么区别，才好帮你分析。
你的桩用C3D15单元，是15-node quadratic triangular prism,
这还是一种二阶单元，不能用于塑性材料和接触问题。在TestModifiedContact.inp中把材料改为弹性，step
1就收敛了（step 1中没有接触, 但step 2有接触，又出现过度约束问题，前面已经讲过解决方法）。
你的桩还是应该像你以前的模型那样划分单元，用C3D8I. 另外，
尽管说主面的网格要比从面粗，但像你现在这样，主面上有很大尖角，还是可能出现收敛问题。原则上，每90度的圆弧上应该划分10个单元，适当少一些可能也行。
pair中的adjust=0.005，还是太小，在后处理时可以看到，接触面之间的距离大于0.005。把adjust设置大一些没关系，比如adjust=0.1。
在施加载荷前，应增加一个分析步，不施加载荷，只是让接触关系建立起来。原则是不要在一个分析步里让模型发生太多变化，那样就把ABAQUS搞晕了。
施加载荷时，如果是塑性材料，初始时间增量步应设小些，比如
0.1, 1., 1e-05, 1.
关于对称面上应力和位移，在后处理时用display group隐藏模型的一半即可。
问：仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentO&x&
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
整個被壓扁了& & 不如預期&
& [1]a3K]5fS‑q/M
照理&&是要&
&基礎&&沉入&&土壤內&
|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent9ZO"X‑uH@
或是&&對土壤&&造成&
&沉陷的效果& &^^
msg 文件中出现overconstraint警告，说明出现过约束。去掉rigid body约束就可以了。
foudation上要加对称边界条件。
长度的单位是什么？重力9.8太大了，和弹性模量的单位要统一。
土体不接触的地方单元可以划粗一些。
你的基础长度100，是100m吗？
数量级还是不对，用国际单位的话，以钢为例，弹性模量是2.1e+11 N/m2, 密度是7800
kg/m3.& 弹性模量如果太小，变形会非常大，难以收敛。
如果想模拟把基础压入土中，需要把土定义为塑性材料，并在接触面附近划分非常细的单元（塑性问题很复杂，你先把弹性作对了，再做塑性）。见
http://www.simwe.com/forum/post/view?bid=31&id=683533&sty=1
http://www.simwe.com/forum/post/view?bid=31&id=683101&sty=1
局部细化网格，要先在mesh模块中用菜单 tools / partition来分割要细化的部位。
问：这个桩基础大致是这样的, 顶部是一个直径在8CM左右的微型混凝土桩, 桩长5M左右. 桩下部有二个螺旋状的钢板,
直径为30CM左右, 厚度为1CM, 深度在地面下6到7米. 钢板焊接在圆形钢柱上, 其直径大约为4CM.案有限元模拟时,
土体为圆柱, 直径6M, 深度取为10M. 有限元分析时, 需考虑桩与土, 钢板与土以及钢柱与土之间的接触.
根据ABA-BA 大侠的见解, 土体和桩可采用C3D8I 单元, 我想请教的问题是,
(1), 本例的薄钢板到底应该采用3-D单元模拟(C3D8I)还是2-D单元(SHELL)模拟好.
(2), 同样, 细长的钢柱应该用何种单元,
(3), 因为各结构之间尺寸相差很大(比如, 钢板的1CM厚度和整个计算域10M相比), 如果全部采用三维实体单元,
整个模型大概为20万个单元, 对如此大的模型进行三维材料非线性分析并考虑接触问题, 大概需要多长的计算时间呢,
(4), 如果薄钢板采用SHELL单元, 如何设置它与土之间的接触面呢, (我的意思是, 当把一个无厚度的钢板放进土中,
如何分离出和钢板接触的二个面),
xile1981答：
1,4.如果要考虑要减少计算时间的话,我认为也可以使用shell单元对钢板进行建模,当定义接触面时使用
*SURFACE, NAME=surface-name
element-set-name, SNEG(or SPOS)来定义顶面或底面.
2.不知道你的钢柱有多长,个人建议还是用实体来建模吧,用beam的话感觉做接触效果不好.
3.计算时间是要根据你的机器来看的,没有一定的结果,以我以前计算过的模型来看,15万到20万个单元考虑接触(应该接触对要比你多,还有些材料非线性)需要一天左右,使用机器为(2-cpu,64bit,8G内存,linux系统)
问：非常感谢XILE1981大侠的回复. 关于设置接触面问题, 可能是我表达的不清楚, 我的意思是, 当把钢板放入土中的时候,
如果用无厚度的SHELL单元来模拟, 钢板上下顶面为MASTER SURFACE, 和钢板相接触的土面为从面,
这二个土体的从面如何得到呢.
是不是用钢板的厚度(1CM)在土中CUT出个空洞, 然后在空洞表面上设置从面, 但钢板SHELL
单元是通过抽取钢板的中性面得到的呢.确如果是这样, 那么从面和主面之间有间隙(钢板的一半厚度), 这就意味着主从面没有接触啊,
这如何解决呢.
xile1981答：哦,是这样啊,那你可以用
*CONTACT PAIR, INTERACTION=interaction_name, SMALL SLIDING,
来调整,使从面上单元节点调整到主面上.但可能会产生单元变形过大,但我觉得钢板厚度与网格尺寸相比较小的话,应该问题不大.
我没有做过土木,只是给你个参考吧,呵呵.
地基和土体的接触
我&&作出以後&
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3iT-EVI
試著& &加載&
&與9000N/m2&
hG:WJo仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent但&&跑出來的圖&
&不如預期& &
圖形如下:W7a+fR
P7@4P9~[1]X
0WGD&f.?|Simwe.com|仿真|设计|有限元|虚拟仪器仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'hH2ME
照理說&&應該&&是&
&基礎&&沉入&&土壤&
&如下圖仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMB;[jr5?
ol2|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMSimwe.com|仿真|设计|有限元|虚拟仪器!{[]/dN$z.~B
你那张成功的图里，应力最大值是10e7,
而失败的图里，应力最大值只有10e5。如果材料相同，应力小，变形当然就小。应力小的原因是载荷小，把载荷变为900000试试。
你那张成功的图里的deformation scale
factor是1吗？如果是的话，这样的大变形要在step里设置几何非线性nlgeom 为 on.
你用了*Solution Technique, type=CONTACT ITERATIONS，不知是否可行。
一般都是先作地应力平衡，再加载荷，不知你这样一步完成对不对。
既有隧道模拟中出现的问题
你的问题在 “ABAQUS常见问题汇总（1.0版）
8.5&&管土/桩土接触 ”
中都提到过，你仔细看过吗：
-&&***WARNING: SOLVER PROBLEM.
NUMERICAL SINGULARITY WHEN PROCESSING NODE
&& & LINER-1.24
D.O.F. 2 RATIO = 3.9.
- 网格太粗，90度圆弧上至少要有10个单元。
- 隧道在厚度方向不能只用一层单元，至少要有3层，用CPE4I单元。
- 塑性问题不能用二阶quad单元，soil单元要用CPE4I.
- contact要增大adjust的值，可以用 adjust = 0.5 .
板料成形的接触问题
请问厚度变化怎么看
在板料冲压过程中，冲压过后，我想看看板料的厚度变化情况，怎么看？|Simwe.com|仿真|设计|有限元|虚拟仪器]b}&V1yJ
我在输出变量设置里面看到有thickness，可是设置后，好像那个不是指板料的厚度
Hc^‑Za‑O1U仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&谢谢大侠指教
体单元没有法子直接输出厚度。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM"t&
Ex0P5i3A+A
其实体单元已经很直观了。
&A5E$~.x&_
y5M仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM可以查看厚度方向上下两个节点间的DISTANCE,www.simwe.comO'_/H*g
当然可能对不正，要采用一定的策略配准。
-y'Es5i&N({h厚度随时间变化的曲线某一点处的还是可以输出的。www.simwe.com&Nr'i#J
很复杂。。。。
板料成形还是壳单元比较划算些，
MtV4BXz;{[1]f‑K仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent一来可省大量的计算成本：如果要得到比较精确的结果，实体单元要划得比较细才行。
h3|[1]ksI9N|Simwe.com|仿真|设计|有限元|虚拟仪器如果用减缩单元，厚度方向要至少划四层才行，而别的方向上单元的尺寸要相对差不多。
h}-J"`就是差不多的一块板子，划实体单元要多用十几倍的单元，再加上最小单元长度变小了好多，又是十几倍的时间，如果用EXPLICIT，相同的问题，要多耗掉100多倍的时间，10多倍的磁盘空间和内存，是很不现实的。
0}(q;u._6Cf1y|Simwe.com|仿真|设计|有限元|虚拟仪器二来用实体单元好象模拟板材起皱有困难，理由见上。|Simwe.com|仿真|设计|有限元|虚拟仪器)e[1]r
三模拟回弹的时候约束不好加，也不容易收敛，至少在我看来是这样|Simwe.com|仿真|设计|有限元|虚拟仪器E)s7n+qa"x#elhN
除非你的板料相对过厚，或者有特别的目的。
问：我在模拟一板料冲压成形过程中，采用sc8r单元，但模拟第一步的0.9999时就出现不收敛的情况，这个问题怎么解决，*.inp文件在附件里面
在msg文件中看到：
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING
NODE SHEET-1.49 D.O.F. 3 RATIO = 4.5.
这说明SHEET-1在方向3上出现了无限大的刚体位移。
你的SHEET-1在方向3上没有边界条件，完全靠接触面上的摩擦来消除刚体位移，但接触关系是在迭代过程中慢慢建立起来的，不是一开始就有摩擦力，所以会出现刚体位移。
可以好好研究一下Getting Started with
ABAQUS中板料冲压的例子，其分析步和边界条件设置得很繁琐，但有其一定的道理。原则是，在第一个分析步中用临时的边界条件来消除所有刚体位移，让接触面上出现一定的压下量或过盈配合，使摩擦力建立起来，在后面的分析步中再去掉临时边界条件。
我试了试，确实很麻烦，不过还是搞出来了，见附件。看来主要还不是sc4r单元有什么特殊性，而是你模型中的其他问题。我做了以下修改：
- 在STEP 1中，固定坯料的整个左端面，这样就不会出现NUMERICAL SINGULARITY
警告信息。在后面的分析步中再去掉此边界条件。（如果只固定此面上的一条边，就不能收敛，很奇怪）.
- 你材料参数的屈服应力是1.37E+10，我想应该是1.37E+8吧。
- 在STEP 1中，让Holder下移0., 对坯料产生微小的压力，以便产生摩擦力，消除刚体位移。
- Holder刚性面上有尖角，在此角点上的法线方向是不唯一的，会造成收敛问题，将其改为圆角。
- 为die和坯料之间的*contact pair加上adjust=1e-07，否则不能收敛。
- 在STEP模块中，菜单other / General Solution Control,
把所有step的*control都恢复为默认值。点击菜单Model / Edit Keyword，在step
1中添加以下语句：
*CONTROLS, PARAMETERS= TIME INCREMENTATION
,,,,,,18,,,,
这是maximum number of severe discontinuity iterations allowed in an
increment. Default is 12, 现在改为18。
diagnostics中看到，接触状态不断变化的是坯料右端底部，可能是这里只有一侧发生接触，造成收敛困难。因此把die的接触面改短，这样还可以缩小接触区域，大大加快计算速度。
- *shell section用simpson积分，9个积分点。这是参考ABAQUS Example Problems Manual
1.3.3 Stretching of a thin sheet with a hemispherical
punch。你原来的Gauss积分我没试。
（经过以上改动，step 1能收敛了）
- Step 2中，不应把holder的所有边界条件都去掉，而应只去掉对U2的约束。
- BC和载荷的amplitude都用默认的Ramp.
（经过以上改动，step 2能收敛了）
- 在step "Contact 2"
中，对坯料施加一个向上的拉力，使其与punch贴在一起，否则无法收敛。下一个分析步中再去掉此拉力。
（经过以上改动，step "Contact 2" 能收敛了）
- 网格在长度方向太粗，大变形时只有很少的节点与punch相接触，导致收敛问题。在长度方向将网格细化，宽度方向改稀疏一些。
- 去掉了*potential, 不过或许有它也可以。
可能还应该改的是，摩擦系数0.162是否太大了？有润滑的话应该是0.1或更小，这样收敛就更容易些。
目前看到的结果是punch下压了分析步的0.05，后面的还在算. 等不及后面的结果了，明天再说吧。
问：在板料成形中出现的错误
在sta文件中看到
Maximum adjustment for initial overclosure is: 5.00001E-04 at
node 1 of instance BLANK-1.
***WARNING: For contact pairs that are not tied, unresolved
&&&&&&&&&&&
overclosures will result in initial accelerations on the order
&&&&&&&&&&&
the penetration distance divided by the time increment squared,
&&&&&&&&&&&
which can lead to significant mesh distortion. These
overclosures
&&&&&&&&&&&
can be avoided by user adjustment of the mesh.
后处理中看到increment 0时坯料就变形了，原因是上述overclosure.
出现这种过盈接触的原因是shell的厚度是0.001,
而模型中刚体和shell的距离是0，即过盈量为0.0005。我认为，解决方法是把刚体移开一些，使其和shell的初始距离为0.0005。见
ABAQUS Analysis User's Manual&
21.3.2 Surface properties for general contact
21.4.2 Surface properties for ABAQUS/Explicit contact pairs。
这样改后就可以运行了，但在2.8e-3秒时，看到坯料局部变形极大，原因可能是holder
force太大，材料不能自由流动。
另外，你用了刚塑性材料，没有硬化特性，不知是否会有问题，见Getting Started with
ABAQUS&& 10.4.4 Adding hardening
to the material model。
请教一个*control的问题
在用abaqus进行模拟时，如果不收敛，可以通过输入文件中用*CONTROLS关键字来更改默认的迭代及收敛准则。在看了好多遍keyword之后，对其中的参数还是高不清楚，希望高手给解释一下，每个参数都是什么意思啊？谢谢了！
用*CONTROL来增大 “允许的最大迭代次数” 不一定能解决问题，只有在以下情况下才有用：
在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛。
我一般只修改 “允许的最大迭代次数”，其他收敛参数修改后效果好像不大好。
具体方法见&
http://www.simwe.com/forum/viewthread.php?tid=693237
问：模型为二维弯曲的一半，板材为三层板。当将层与层中间的inter接触性质改为tie时，就可以收敛，但是我想计算一下在超过一定剪切强度下的摩擦情况。也就是说当界面剪切强度超过1E6Pa时，发生摩擦，摩擦系数为0.1的情况，大家帮我看看，我该怎么调整（Interaction
中的tolerance已调整为1E-8）
中间层为什么用三角形单元，那么小的尖角，单元形状太差。
尝试把材料都改为steel, 线弹性，摩擦属性不用taumax, 但仍然不收敛，说明不是这方面的问题。
问题是：接触面上单元非常多，需要增大允许的迭代次数。方法是：step模块，菜单other / Gerneral solution
controls, 点Edit，选择第一个分析步"establish contact I"，点continue,
选择specify，点Time Incrementation 标签页，点第一个more，把Is由默认的12改为30，点OK。
相应的inp文件：
*Controls, parameters=time incrementation
, , , , , , 30, , , ,
含义是：30次迭代不收敛时，才减小时间增量步。
这样改后就收敛了。另外，可以尝试把单元划得粗一些，就可以收敛得更快。
reduced integration单元要求在厚度上至少有4个单元。建议改为CPE4I单元。
问：这是个简单的拉伸成形模型，有几个问题没有解决：
（１）局部的单元变形太大，怎么控制让领近单元变形均匀点？
（２）我的回弹模型出了什么问题，一直没有成功，希望大侠能给指点一下，问题出在哪里！！
wy_anli答：
在成形过程中，那个应变比较大的单元已经超过了你定义材料的最大应力，也就是说当应力达到你定义的材料最高应力时，就会把材料当成超塑性体。
（见图中的5.856E8,你定义的最大应力为4。7429E8）
建议减小压边力
因为你采用的是显示算法，所有的加载应注意平滑加载，在计算完成后要看一下动能的历史输出，要平缓，另外比较动能与内能，动能/内能一般要小于8%
另外那个BLANK_UP_FORCE可以减少甚至去掉，因为BLANK_UP_FORCE会增大成形难度。
2.回弹要在隐式算法中做
你的那个图是变形完成后的应力情况，在成形的开始的部分应力已经超过了极限值，你可以从开始看一下应力变化情况，也可以看我传的那个图
你先把成形的部分模拟完毕（修改压边力，使成形合理），再进行回弹模拟
问：这个模型为什么不收敛（超弹和弹塑性接触）？二维模型为拉弯三层板，因为是对称的，所以只取其中一半。上左为凸模，上右为压料板，下为凹模，中间为三层板。
三层板中间为tie接触，上下两层为steel，中间那层板为hypereslatic超弹性材料。
附上inp文件
用线弹性材料试试，如果仍不收敛，就是模型参数定义得不对。见
http://www.simwe.com/forum/viewthread.php?tid=693237
http://www.simwe.com/forum/viewthread.php?tid=680814
http://www.simwe.com/forum/viewthread.php?tid=684042
http://www.simwe.com/forum/viewthread.php?tid=679059
问：（1）如果进行多道次的成形模拟，一般是怎么进行的？是不是：工序1（成形－回弹）－传递到－工序2（成形－回弹）－......这进行？
（2）第一道次的模具都可以在CAE里面完成设计等，哪后面工序用到的模具，有什么好办法添加进去：－直接在.INP文件里面修改？
先谢谢了！
可能可以在一个cae模型中把所有道次的模具都建在一起，然后在每个step中用*model
change来杀死/激活不同的单元和接触。
的接触问题
问：大家好，我现在做的是一个有压管道在外加作用物下出现凹坑的准静态分析
当然这里面存在接触（外加作用物和管外壁之间），不要接触直接加力我可以做出来，但是加了接触之后就不行了，并出现以下
警告：A slave surface definition has been defined on a face of element
1619 that contains no midface node. This may lead to convergence
difficulties. The element type will be converted and a midface node
generated automatically. The element number will be unchanged. See
the output from the preprocessor to confirm the new element and
node definitions. This message will not be repeated if subsequent
elements of this surface are converted
18 nodes may have incorrect normal definitions. The nodes have been
identified in node set WarnNodeIncorrectNormal.
The aspect ratio for 51 elements exceeds 100 to 1. The elements
have been identifiein element set WarnElemAspectRatio.
58 nodes have dof on which incorrect boundary conditions may have
been specified. The nodes have been identified in node set
WarnNodeBCIncorrectDof.
58 nodes have dof on which incorrect boundary conditions may have
been specified. The nodes have been identified in node set
WarnNodeBCIncorrectDof.
错误：118 elements have missing property definitions. The elements
have been identified in element set ErrElemMissingSection.
Analysis Input File Processor exited with an error.
在后处理中可以用display
group来显示错误信息中提到的那些集合WarnElemAspectRatio，WarnNodeIncorrectNormal，就可以看到问题在哪里。
- 接触问题和塑性材料不能用二阶quad和hex单元，圆管可以用C3D8R，但在厚度方向至少要有4个单元.
- 刚体建1/4模型没错，你定义刚体的基本方法也是对的，但*rigid body 的单元集合PICKEDSET10不存在,
需要重新定义一下.
- 施加了刚体约束就要把力加在参考点上。
- 参考点有旋转自由度，应定义边界条件来约束U1,U3,UR1,UR2,UR3。
- “118 elements have missing property
definitions”：施加刚体约束的单元也需要定义section，材料可以是任意的（比如steel）,类型是shell.
- “18 nodes may have incorrect normal definitions”：接触面的master
surface上不能有尖角，否则尖角处的法线方向不唯一。因此，主面应只选择可能发生接触的面，即indentor的底面和1/4模型的两个侧面，且三个面之间都要有过渡圆角，不能有尖角。另外90度圆弧上要至少划分10个单元。见附件。
- “The aspect ratio for 51 elements exceeds 100 to
1”：圆管底部网格太稀疏，单元过薄。不改可能也能算下去，这里不是大变形的关键部位。
- 你的塑性参数中的塑性应变怎么到了1360？到2就足够了。
- “58 nodes have dof on which incorrect boundary conditions may
have been specified”：这只是提示信息，一般不用理它。
问：刚性体与弹塑性性的接触问题
我做的是个管道凹坑成型问题，
/A&S'iX6d&t-]怎么算不了
a仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM这是我的cae文件仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent8r8cmLxo
麻烦大家给我看看仿真分析,有限元,模拟,计算
msg文件中有
&***WARNING: SOLVER PROBLEM. NUMERICAL
SINGULARITY WHEN PROCESSING NODE
&&&&&&&&&&&&
ASSEMBLY.1 D.O.F. 3 RATIO = 7..
应该约束刚体参考点的U1, U3, UR1, UR2, UR3.
管道大变形处的网格太粗。slave surface的网格要比master surface细。
C3D8单元精度很差，应该用C3D8R.
问：我做的一个素混凝土梁,四点弯曲试验.用的是smeared cracking模型,各项参数详见附件文件.
前半段计算正常,可以看见上升趋势以及峰值.问题是过了之后出现典型的突降,计算开始不收敛,最后因为步长太小而中止.
改用更小的步长,划分更小的单元,都不奏效.我自己的分析是因为中部的单元被拉破坏,从而导致计算不收敛?可是从11方向的应力看不出破坏的迹象(否则应该11方向为0),梁底部依旧有拉应力,这是怎么回事?
请参考ABAQUS Example Problems Manual 1.1.4 Indentation of an
elastomeric foam specimen with a hemispherical punch。
你的模型有几个问题：
- 一般的大变形可以用ABAQUS/Standard, 如果想模拟材料破坏就要用ABAQUS/Explicit，见ABAQUS
Analysis User’s Manual 第11.6节Progressive damage and
failure。可以查看一下PEEQ是否达到了破坏应变。
- 接触问题和塑性材料都不能用二阶quad或hex单元, 应该用CPS4R.
- 大位移大变形问题，要在step中把nlgeom设为on.
接触区域的网格要非常细化，沿着压头的180度圆弧应至少有20个单元，原则是变形后单元不能出现过度扭曲变形。可以在此区域partition一下，作局部细化。
- 接触从面应再扩大至少一倍，因为下压后需要和整个压头半圆面接触。
- 压头上不能有尖角，否则在尖角处的法线方向不唯一，会造成收敛问题。可以把压头做成360度的圆或类似于试管底部的形状。
- 你的 contact property 用SOFT/AUG. LAGRANGE,
不知是否有必要，一般用默认的hard就行。
问：阿巴，针对该模型，我调整如下，
单元换用CPS4R.
nlgeom设为on.
你建议接触区域的网格要非常细化，我开始也是这么做的，不过这个模型很奇怪，我采用非常细的网格计算还不如粗网格模型计算的好，反而因为穿透过多早早就中止了计算。
另外，我的接触相互作用之所以选择contact property 用SOFT/AUG.
LAGRANGE，主要是因为硬接触带来不少问题，所以为了计算能进行下去，采用在一定距离内发生压应力。
对于压头的形状，我也是尝试了很长时间，因为是理想刚体，不可能做成你建议的形状（假如你有办法请告诉一声，以便学习），我改来改去只有这个形状是软件可以接受的。不过我试过很多不同的网格划分，好像不是关键，因为结果都不是很好，计算全部都提前中断。
最后问一下你建议的ABAQUS/Explicit问题，我的第一个步骤是预压，采用ABAQUS/Standard，正式施压是第二个步骤，我想删除该步骤，重新建立一个，怎么找不到Explicit？6。5换了方法？
网格肯定要细化，不收敛是其他方面的问题造成的。
二维刚体的形状可以是任意的，想画圆的话，分别画四段90度的弧即可。
静态分析步后面不能接着explicit，要先删掉所有静态分析步。
问：请问一下这样的问题能不能做：下图所示为一实体的中截面，打算这样计算，先让物体在1 位置，计算完了之后移动一下到 2
位置再计算！模拟重叠宽度，由于模型的接触力不好简化，想直接利用接触来做，不知道可不可以！请高手指点一下！
移动到 &2 &位置:
&设定位移边界条件, &使x方向位移为某个值即可。
定义接触时要用finite &sliding.
&要把step中的NLGEOM（几何非线性）设为on。
可分为几个Step:
Step 1: 在 1 位置，在其质心施加一个很小的力，从而使接触平稳地建立起来。
Step 2: 在 1 位置，在其质心用一个随时间变化的力，作用一段时间（如 5s）;
Step 3: 载荷变为0，去掉接触（否则迭代过程会很慢），设定位移边界条件, 令x方向位移为某个值（即使其移动到 2
Step 4: 在2 位置，恢复接触，在其质心施加一个很小的力，从而使接触平稳地建立起来。
Step 5: 在质心作用同样的力。
Step 6: 载荷变为0，去掉接触。
如果用不考虑冲击效应和运动的惯性（即缓慢加载），就用静力分析；否则用动态分析，implicit 或 explicit
都可以，但加载时间5s的话，用动态分析可能计算时间很长。
可以重叠建立实体并划分网格.
去掉接触: 在abaqus/CAE的模型树中，可以把接触在某个step中 suppress.
上方实体除了X方向的移动，还可以向下移动一点点，以保证和下方实体相接触。下方实体不动。
问：刚体穿透，aba大侠看过来 ~
模型很简单，就是一个2d的圆筒受压。采用shell单元。麻烦大家看看，多谢了仿真分析,
删除interaction property中的normal behavier,
就是默认的hard接触，不允许穿透。
后处理时比例设为1, 就看不到穿透。
多谢，Aba,和前面一位大侠，问题解决了。都怪我没有注意deformation scale
factor是否为1，误以为穿透了。分享一点
W-Tk1^-VpSimWe仿真论坛www.simwe.com7V$z@$Sc8ax
刚体穿透解决办法SimWe仿真论坛:y6ZMUwO Y.z
;bX+d^,]v,s1.边界条件和模型,材料是首先需要检查的地方,通产都是最容易发生问题的
A:Q?4cK@nh0f%coSimWe仿真论坛2.尽量使用analytical rigid,而不是 discrete
oQ'f;av|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent3.检查刚体和变形体的法向方向是否正确|www.Simwe.com|有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|FluentD5W)y#q4W9IR
4.网格划分是需要fine meshed,从面要比主面细& &
& & D;h.f%oB$j0?6qQ
5.接触时主从面务必设置好,从面尽量使用bode-based element 而不是element-based
K.l3J仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMyccNB
[ 本帖最后由 impactboy 于
问：做得是三维接触分析,主要想知道接触区接触力的具体数值,以便进行下一步单个零件的受力分析.
Step模块，设置History Output，把Domain设为Interaction，选择接触对，选中CFN。
用 *CONTACT PRINT 可以将接触力CFN输出至DAT文件：
*CONTACT PRINT, SLAVE=从面名称
问：我用二维模型，椭圆物压平面物，计算得到了几个接触部分的history
output结果，发现跟实际情况数量级相差太远，一个是由接触压力得到的综合接触力（CFTM), 接触压力（cpress),
接触面积变化（CAREA)，而且接触面积变化不连续。再就是field output
中S子目录下pressure是什么压力，如果取接触面一点看的话，和接触压力是什么关系。
我主要是想得到接触体之间的相互作用合力大小，接触压力，单位为mm.
xile1981答：
你要得到相互作用力大小的话看CFN或者CFT(带摩擦的)就可以,如果跟你实际情况差的太远的话可能是你的模型有问题,检查一下接触的定义吧.
s子目录下的pressure是(equivalent pressure stress),
和接触压力(cpress)并不是一个概念,具体介绍可见手册.
答：法向接触力要看CFN(注意根据方向来选择CFN1,CFN2,CFN3),
你说的CFTM是接触压力和摩擦应力的合力，CPRESS是单位面积上的法向接触力（单位是 力/面积,
和接触力的数量级当然不一样）}
我要回帖
更多关于 tid口服 的文章
更多推荐
- ·我像韩红点评过的一首歌歌名像是个人名吗？
- ·升级两个人玩的游戏有哪些怎么玩
- ·因为比较胖，多 囊 多囊卵巢最快受孕方法不孕吗？
- ·有哪些值得推荐的王者代打电脑软件推荐？
- ·有哪些值得推荐的王者代打电脑软件推荐？
- ·碰到这种没眼睛的简述行政主管职责我是不是要去找老板好好谈谈
- ·求解啊各位世界杯大神预测6！
- ·亚预赛菲律宾vs菲律宾澳大利亚第三节节打群架，蓝球打了三节打完结束，第四节没打中断外围怎么结算？
- ·一部美国老的探讨成年轻人影片免费播美国电影的电影
- ·电路图中虚线表示什么有一条虚线
- ·平安车主信用贷易贷对信用的要求高吗
- ·10Kw能30kw用多大的电缆电缆
- ·梦见往玉米地里浇水昨天下午刚打了百草枯今天就浇水了，对玉米有影响吗
- ·每天薄薄一薄生平片黄瓜片会胖吗
- ·贵阳护照办理地点7个办理护照的地点在哪里？
- ·女友家里没人，她女朋友让我去玩吧她家，她平时很清纯很可爱，我们今年都高考过，应该能去一座城市
- ·你喜欢钟林，最喜欢的用英语怎么说说？
- ·关于这个日本搜查官系列图解的图，叫什么，在哪儿搜
- ·痔疮痔疮便血吃什么食物好好？
- ·童心宝腹泻贴使用的秋季腹泻饮食注意事项项是什么呀？
- ·平稳，是tid，还是did？
- ·国民经济基础是什么问题求详细解答～
- ·求请老师指点怎么说合适?点
- ·求missin' you like crazy 结局
- ·305匹压缩机工作电流电流是多少
- ·谁知道田布扁鹊仓公列传译文的翻译 谢谢
- ·怎么让如何哄睡小月龄宝宝宝仰睡
- ·反流性食道炎胸骨疼会导致胸骨疼痛吗，怎么根治
- ·低烧就是说免疫力低下会低烧吗.现在又开始发烧了我该怎么办
- ·没有流鼻血可颈部有肿块是鼻咽癌颈部肿块位置图吗
- ·中考考什么27号考什么时候
- ·整天活在恐惧里当中，非常恐惧自己的未来，担心自己活不了多久，该如何的自我调节，自我救赎呢？
- ·机在工地干活挖掘机，我挖掘机的话干完了，要不到钱怎么办
- ·职工集资款利息怎么算承诺的利息不还怎么办
