abaqus一个从初学到精通粘弹性的分析的经验积累.docx
《abaqus一个从初学到精通粘弹性的分析的经验积累.docx》由会员分享，可在线阅读，更多相关《abaqus一个从初学到精通粘弹性的分析的经验积累.docx（19页珍藏版）》请在冰点文库上搜索。
abaqus一个从初学到精通粘弹性的分析的经验积累
问题积累〔待续〕
1.abaqus如何调整图例的大小，就是云图左上角那个图框，字太小了看不清！
！
直接设置图例的字体大小就可以：
工具栏viewport>viewportannotationoptions>legend（选项卡）>text（选项）>setfont（按钮）>size，修改size选项中的数字，就可以修改图例大小了。
2.cohesiveelementＡＢＡＱＵＳ
在6.11使用cohesiveelement，定义cohesive材料属性的时候主要步骤：
1.定义一个材料的名字，比方cohesive，不要去定义任何属性〔弹性，弹塑性等等〕。
2.翻开工具栏model--editkeywords，在inp中手动添加材料的各种属性。
PS:
定义section的时候选cohesive，elementcontrol选sweep，elementtype选cohesive，这些是使用cohesiveelement的根本步骤。
zerothickness的cohesivesection设定abaqus所谓的zero-thickness，其实就是定义cohesivesection的initialthickness=1.0。
你可以在定义section的时候定义〔specify〕，也可以用系统默认的thickness〔也是1.0〕，这样有关cohesiveelement的计算当中，就有displacement〔位移〕=strain〔应变〕*thickness（1.0）=strain的数值。
我们知道从1914年Ingless和1921年Griffith提出断裂力学开始，一直到60年代都停留在线弹性断裂力学（LEFM）的层次。
后来由於发现在裂纹尖端进入塑性区后用LEFM仍然无法解决stresssingularity的问题。
1960年由Barenblatt和Dugdale率先提出了nonlinear/plasticfracturemechnics的概念，在裂纹前端引入了plasticzone,这也就是我们现在用的cohesivefracturemechnics的前身。
当时这个概念还没引起学术界的轰动。
直到1966年Rice发现J-integral及随后发现在LEFM中J-integral是等于energyreleaserate的关系（赞大牛Rice一个，在读Ph.D期间就推导出对断裂力学有重要意义的J-integral了。
。
。
，实令我辈汗颜）。
随后在工程中发现了越来越多的LEFM无法解释的问题。
cohesivefracturemechnics开始引起更多的关注。
在研究以混凝土为代表的quassi-brittlematerial时，cohesivefracturemechnics提供了非常好的结果，所以在70年代到90年代，cohesivefracturemechnics被大量应用于混凝土研究中。
目前比拟常用的方法主要是fictitiouscrackapproach和effective-elasticcrackapproach或是称为equivalent-elasticcrackapproach.其中fictitiouscrackapproach只考虑了Dugdale-Barenblattenergymechanism而effective-elasticcrackapproach只考虑了基於LEFM的Griffith-Irwinenergydissipationmechanism，但作了一些修正。
简单来说，我读了一下ABAQUScohesiveelement的理论帮助，个人觉得ABAQUS的cohesiveelement采用的是广泛应用于混凝土的类似fictitiouscrack的方法。
只考虑了Dugdale-Barenblattenergymechanism。
这其中softeninglaw的影响是非常重要的。
但ABAUQS似乎只提供了linear或者exponential的softeninglaw，复杂的本构关系还需要另想方法。
至於基於Griffith-Irwinenergydissipationmechanism的J-integral值可以在LEFM分析中单独算。
（ABAQUS用的是SuoZhigang和Hutchinson在1990一篇论文中提出的方法）目前cohesivefracturemechnics已经被应用于各种材料。
不过在使用到纳米或者更小数量级的研究中碰到了不少问题，可能需要结合位错和分子动力学的一些理论，我现在的研究中也碰到类似的一些问题，希望和大家一起讨论。
做裂纹ABAQUS有几种常见方法。
最简单的是用debond命令,定义*FRACTURECRITERION,TYPE=XXX,参数。
。
。
***DEBOND,SLAVE=XXX,MASTER=XXX,timeincrement=XX0,1,……......time,0要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset，然后在*INITIALCONDITIONS,TYPE=CONTACT中定义master,slave,及指定的Nset这种方法用途其实较为有限。
另一种方法，在interaction模块，special,定义crackseam,网格最好细化，用collapseelement模拟singularity.这种方法可以计算J积分，应力强度因子等常用的断裂力学参数.裂尖及奇异性定义:
在interaction-special,先定义crack,定义好裂尖及方向,然后在singularity选择：
midsidenodeparameter:
输入0.25,然后选Collapsedelementside,duplicatenodes，8节点单元对应（1/r）+（1/r^1/2）奇异性。
这里midsidenodeparameter选0.25对应裂尖collapse成1/4节点单元。
如果midsidenodes不移动到1/4处,那么对应（1/r）奇异性,适合perfectplasticity的情况.网格划分:
裂尖网格划分有一些技巧需要注意，partition后先处理最外面的正方形，先在对角线和边上布点，记住要点constraint,然后选第三个选项donotallowthenumberofelementstochange不准seed变化，密度可以自己调整.最里面靠近圆的正方形可以只在对角线上布点.也可以进一步分割内圆及在圆周上布点.里面裂尖周围的内圆选freemesh,elementtype选cps6或者cpe6，外面四边形选sweepmesh,elementtype选cps8或者cpe8,记住把quad下那个缩减积分的勾去掉。
补充一下这种方法的几个常见问题，见不少朋友问过。
主要是对断裂力学的理解问题，顺便一起答复一下。
1．为什么我设置理想弹塑性（epp）分析的时候得到的xx,yy方向或者最大应力值Sxx,Syy会超过材料的屈服强度Sy呢,这分析结果可能吗？
这是因为在ABAQUS中对应等于材料的屈服强度的是vonMises等效应力Se=Sy，因此在平面应变的条件下，xx方向的应力Sxx=Sy*pi/SRQT（3）>Sy,而Syy=Sy*（2+pi）/SRQT（3）,大概是3倍的屈服应力。
所以得到大于材料的屈服强度的xx及yy方向应力是正常的。
2．为什么设置collapseelement的时候对弹性分析在中间就一个点而要把单元边上的中点移到1/4处，但弹塑性分析却要在中间设置一圈点并且保持单元边上的中点位置不变呢？
这个其实不是随便定的，在有限元中分析裂纹时，对弹性分析需要模拟裂尖1/SQRT（r）的奇异性，这样在把单元边上的中点移到1/4处后计算出来的等参单元拉格郎日型函数对应的ufield正好包含1/SQRT（r）项，事实上这一方法在断裂力学的数值模拟开展史上是很巧妙的一个发现，至今仍然被广泛采用。
至于理想弹塑性分析需要模拟裂尖1/r的奇异性,这样大家都知道在把单元边上的点放在到1/2处后计算出来的正常的等参单元拉格郎日型函数对应的ufield包含1/r项,可以模拟弹塑性分析需要的裂尖1/r的奇异性。
所以在看似动手点几下就能实现的分析模式后面有很清楚漂亮的理论作支持。
还有就是比拟新的cohesiveelement单元。
需要定义damageinitiation和evolution的准那么，softening准那么目前只有linear和exponential，但对一般材料也够用了。
然后通过设置后处理displaygroup可以看到裂纹扩展情况。
裂纹扩展不是ABAQUS的强项，目前比拟方便的只能用cohesiveelement，我做过几个模型效果还可以，但对应的参数需要一定的实验数据支持，否那么做出来了也不知道对不对。
或者使用python控制seam尺寸,然后移动partition和网格，比拟麻烦。
3.分析步
分析步主要有初始分析步和后续分析步，每个分析步可以用来描述一个分析过程，例如在后续分析步中施加不同荷载，在初始分析步中施加边界条件等。
4.增量步
增量步是在分析步里面根据模型计算收敛情况设置的，简单模型可以设置较少的增量步，并可使初始增量为1；复杂模型设置多一点增量步，并减少初始增量值。
超过设置的允许增量步数，那么计算停止。
补充；静力分析时，分析总时间（Timeperiod），是个相对概念，没有单位的。
尽管静力分析，你的荷载可能是一下子定义施加上了，但是aba计算时，数值积分是一步一步迭代计算的，需要N个迭代步，每个迭代步都需要求解联立耦合的方程组，程序计算需要一定的时间，这个时间就是时间增量increments，假假设一共迭代N步，这N个increments之和就是分析总时间（Timeperiod）。
一般情况下，程序一开始按照你设定初始时间增量步进行计算，
1〕假设是在这个初始时间增量步下不能收敛，就会减小时间增量直到减小到你设定的最小时间增量为止，假设还是不能收敛，那么程序就会给你报错。
2〕假设是按照初始时间增量步可以顺利计算，程序会自动调大时间增量，一边加快计算进度，时间增量步调大的上限就是你设定的最大增量〔Maximumnumberofincrements〕。
故而，最大增量步，初始时间增量，最小时间增量，最大时间增量等你的预先设置，在很大程度上影响着程序能不能顺利计算和计算速度的快慢。
5接触问题
接触问题，如果一开始并没接触上，是有间隙的，那么第一步就是要接触上。
可以用一个小载荷，但用一个位移，比初始间隙稍微大一点点，或许更好。
然后在后继步中，去掉这个位移，并加上载荷。
甚至可以加一个大一点的位移，用比拟小的增量步，多输出几个增量步的结果，查看接触力的增加。
在快要到达所应加的力的地方，做一个重启动，去掉位移，把最后的力加上去。
加位移要比加力更容易收敛。
6.abaqus中的收敛
要取得一个精确的结果需要以下这些：
a,网格的收敛
b,瞬态问题的精确的时间整合
c,非线性求解过程适当的收敛
另外，一个精确的结果还需要在创立一个适宜的有限元模型时好的工程评价，包括材料，载荷，边界条件以及求解过程。
补充：
如果不收敛，一般使计算陷于死循环，为了防止这种情况，程序里都参加是否收敛的判断.
有限元法也是近似解法，收敛是作为计算结束，结果无限接近真实结果的一种判定标准。
7.在abaqus配置完成后，点击运行，如果出现ServerStartFailed.TheServerMayAlreadyBeRunning!
解决方法是先点击startserver右边的stopserver，然后再点击startserver，如果配置没错的话，这里会ServerStartsuccessful的。
8如何显示梁截面〔如何显示三维梁模型〕
a〕无论是运算还是默认显示，ABA中的梁都是一条线，很多人想看梁截面〔一般一个星期有人问一次〕。
显示梁截面：
view->assemblydisplayoption->renderbeamprofiles，自己调节系数
:
//forum.simwe/viewthread.php?
tid=835478&page=1#pid1531086
b〕后处理到底能不能显示梁截面？
在deformedshape和undeformedshape都能用上面的方法显示梁截面，在应力云图〔contour〕不能显示。
c〕也经常有人问起如何显示壳单元的厚度
:
//forum.simwe/thread-865887-1-1.html
9怎么在局部坐标系下建立参考点
在前处理中，已经建立了局部坐标系，如何在局部坐标系中建立参考点？
这个有点麻烦，看看konadoul图文并茂的例如吧。
:
//forum.simwe/viewthread.php?
tid=863389&highlight=%D7%F8%B1%EA%CF%B5
10Documentaion〔help文件〕不能搜索
首先保证你准确的安装了Documentaion〔先安装Documentaion再安装程序〕，其次有问题你可以重新安装
一次Documentaion。
如果你和我一样比拟懒不想安装，看看下面的方法是否管用吧。
1〕控制面板---效劳找到texis双击查看是不是automatic，如果不是就设置为automatic
2〕你可以用这个：
:
//name:
2080/v6.8/注：
name是你的计算机名；6.8是版本号，
比方你用6.6的就改为6.6.〔我的在自从不能搜索之后我就一直这么用的〕
:
//forum.simwe/viewthread.php?
tid=861085&extra=page%3D12%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
:
//forum.simwe/viewthread.php?
tid=698879
11建立几何模型草绘sketch的时候，发现画布尺寸太小了
1〕这个在createpart的时候就有approximatesize，你可以定义适宜的〔比你的定性尺寸大一倍〕；
如果你已经在sketch了，可以在edit菜单－－sketchoption－－grid更改
2〕这里如果你选择constriant标签，还能更改尺寸精度
12想输出几何模型
part步，file，outport－－part
想导入几何模型？
part步，file，import－－part
13如何定义局部坐标系
Tool－CreateDatum－CSYS－－建立坐标系方式－－选择直角坐标系or柱坐标系or球坐标
14如何在局部坐标系定义载荷
laod－－Editload－－CSYS－Edit〔在BC中同理〕选用你定义的局部坐标系
15如何定义随变载荷
amplitude这个不多说了，强烈建议看?
常见问题2.0?
小康大侠图示空间变载：
:
//forum.simwe/thread-867236-1-5.html〔强烈推荐照此演示操作一回〕
:
//forum.simwe/viewthread.php?
tid=861727&highlight=%CB%E6%CA%B1%BC%E4%B1%E4%BB%AF%B5%C4%D4%D8%BA%C9
設定於空間中變化的負載：
A;v%|!
:
//forum.simwe/viewthrea...p;page=1#pid1556636
加载梯形载荷：
:
//forum.simwe/viewthread.php?
tid=870350&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
16怎么知道模型单元数目〔一共有多少个单元〕
在mesh步，meshverify可以查到单元类型，数目以及单元质量一目了然
Query---element也可以查询的。
:
//forum.simwe/viewthread.php?
tid=813860&highlight=%B5%A5%D4%AA%CA%FD%C4%BF
:
//forum.simwe/viewthread.php?
tid=858700&highlight=%B6%E0%C9%D9%B8%F6%B5%A5%D4%AA
17怎样查询某个面的面积
下面是好心人给出的几种方法〔有图文版〕，看看你喜欢那种了
:
//forum.simwe/viewthread.php?
tid=855289&highlight=%C3%E6%BB%FD
:
//forum.simwe/viewthread.php?
tid=859665&highlight=%C3%E6%BB%FD
18分享一个Documentation搜索文件小技巧：
1在操作过程中有疑问需要搜索Documentation文件，如果你在相关CAE界面直接按F1，这样
Documentation会自动链接到相关领域，省去搜索、查找的麻烦〔你搜索一个关键词，会出来
一堆相关的东西，到底哪一个是你最想要的呢？
〕
2在搜索的时候，搜索一个词组〔比方hardcontact〕就采用"hardcontact",如果你输入hardcontact，
系统就默认为是你想搜索hard和contact，会出来很多你用不着的东西，反倒把你想要的给淹没了。
加不加引号就是这个区别〔注意，必须是英文输入法的引号""哦〕
19模型显示dependentpart不能mesh，怎么办？
1〕你可以点模型树在part步mesh，也可以
2〕在模型树assembly下点相应的part，右键，makeindependent
20想隐藏一些part以便更清楚的看见其他part，edge等
view－AssemblyDisplayOptions，打勾
21想打印或者保存图片
file－print－OK就打印
file－print－把Destination改为file并指定路径，format可选，就可以把当前模型视图保存为图形文件〔我一般用png〕
22我建立了实体模型，实体单元只有3个自由度，如何施加弯矩〔or扭矩〕？
常见问题2.0有说。
:
//forum.simwe/viewthread.php?
tid=677582&extra=page%3D8%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
如何施加扭矩〔例如〕
:
//forum.simwe/thread-864413-1-1.html
23已经建好2D模型，想得到3D模型，怎么实现？
:
//forum.simwe/viewthread.php?
tid=854372&extra=page%3D10%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68&page=1
24CAE文件太大，上传不便，可否缩小？
保存种子删除网格，另存为新文件并关闭之，他就变小很多了。
另外，看看高人如何文档瘦身的：
:
//forum.simwe/viewthread.php?
tid=861908&extra=page%3D10%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
25你喜欢用Hypermesh划网格，用aba提交运算？
看看老向的好东西吧
很好的HM转aba的视频〔不过建议高手就这个问题收集一个专题，呵呵，呼唤下〕
:
//forum.simwe/viewthread.php?
tid=860878&extra=page%3D13%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
在Hyperworks版置顶的这个由binling_love大侠总结的专题
:
//forum.simwe/thread-862441-1-1.html〔HyperMesh与ABAQUS的接口模板（汇总〕〕非常详细！
26在CAE中定义节点集，如何将外表所有节点全部选中？
請看echo_shanghai及garethz的方法.
:
//forum.simwe/thread-774697-1-1.html
27如何控制加载的数值〔如何按照你的加载步长来加载而不是按照aba默认的
步长加载〕
:
//forum.simwe/thread-864942-1-1.html
28如何更改CAE界面默认颜色
view->Grahphicoptions->viewportBackground->Solid->choosethewitecolour!
然后在file->saveoptions...
29如何增大允许的迭代次数
这个一般情况下不需要更改，除非你有特殊需要才改的。
方法：
在step模块，菜单other---Gerneralsolutioncontrols--Edit，选择第一
分析步"establishcontactI"，点continue,选择specify，点TimeIncrementation
标签，点第一个more，把Is由默认的12改为30，点OK。
30提示“磁盘空间缺乏〞怎么办？
UNABLETOCOMPLETEFILEWRITE.CHECKTHATSUFFICIENTDISKSPACE31AVAILABLE
看看这2帖的方法
:
//forum.simwe/thread-866003-1-1.html
:
//forum.simwe/thread-862120-1-1.html
32abaqus有没有镜像〔mirror〕功能，abaqus如何缩放模型大小
镜像：
模型树右键点copy，改下面相应设置：
:
//forum.simwe/viewthread.php?
tid=261578&extra=page%3D3%26amp%3Bfilter%3Ddigest
缩放：
模型树copypart>scale
33abaqus如何加载重力：
可以用CAE操作，你也可以编辑keywords实现：
:
//forum.simwe/viewthread.php?
tid=45661&extra=page%3D4%26amp%3Bfilter%3Ddigest〔老帖子〕
:
//forum.simwe/thread-867007-1-2.html〔新帖子〕
:
//forum.simwe/viewthread.php?
tid=867007&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D68
加载体力〔惯性力〕呢？
和加重力很类似的〔重力是体力中常见的一种嘛〕：
load---bodyload
34ProE.导入ABAQUS.CAE
konadoul大侠图文并茂：
:
//forum.simwe/viewthread.php?
tid=865114&extra=page%3D1%26amp%3Bfilter%3Ddigest
怎样把Autocad导入abaqus：
AUTOCAD必须保存为dxf格式，再file
升级会员 