ABAQUS常见问题汇总 - 1.0版
作者:aba_aba
2006.06.14
在SIMWE论坛注册5个月以来,在ABAQUS版回帖860篇,知无不言,言无不尽。在此整理一下,以方便大家查找。下面绝大多数是我的回帖,有少部分是其他网友的回答,都注明了作者。
下面只摘录了帖子中的一些主要内容,有些地方可能上下文不太连贯,完整的讨论请大家根据相应链接去论坛上察看。 如果这些内容对你有帮助,希望你也能常上论坛来回答一下别人的问题。众人拾柴,火焰才高。
目录 点击小节标题,可以跳到相应的内容(有些WORD版本可能需要按住ctrl键)
1. 论坛提问要诀................................................................................................................................... 4 2. 故障诊断基本方法 ........................................................................................................................... 4 3. 约束刚体位移 / NUMERICAL SINGULARITY ..................................................................................... 4 4.过约束(Overconstraint) .............................................................................................................. 7 5.其他常见error和warning信息....................................................................................................... 7 6.ABAQUS无法运行(安装、License、环境变量) ....................................................................... 10 7.弹塑性分析 .................................................................................................................................... 14 7.1 塑性问题不收敛的常见现象 ..................................................................................................... 14 7.2 接触问题和塑性材料不要用二阶单元 ....................................................................................... 14 7.3 不要在塑性材料上施加点载荷 ................................................................................................. 15 7.4 塑性材料参数 ........................................................................................................................... 15 7.5 屈服准则 .................................................................................................................................. 18 7.6 后处理 / 判断材料是否屈服 ...................................................................................................... 20
8. 接触分析 ........................................................................................................................................ 22 8.1 接触分析不收敛的常见现象和解决方法 ................................................................................... 22 8.2 接触面上的网格密度 ................................................................................................................ 22 8.3 接触面的法线方向 .................................................................................................................... 23 8.4 过盈接触 .................................................................................................................................. 23 8.5 管土/桩土接触 .......................................................................................................................... 25
8.6 板料成形的接触问题 ................................................................................................................ 30 8.7 凹坑成型的接触问题 ................................................................................................................ 34 8.8 刚体穿透 .................................................................................................................................. 37 8.9 接触力 ...................................................................................................................................... 37 8.10 Explicit中的接触问题............................................................................................................. 38 8.11 改变接触状态 ......................................................................................................................... 40 8.12 其他与接触有关的Error和Warning信息 .............................................................................. 41
9. 各个量的单位................................................................................................................................. 42 9.1 各个量的单位要相互一致 ......................................................................................................... 42 9.2 时间/增量步的含义 ................................................................................................................... 44
10.导入 / 修改 / 求解inp文件 .......................................................................................................... 45 10.1 把inp文件导入ABAQUS/CAE ............................................................................................. 45 10.2 运行inp文件.......................................................................................................................... 46 10.3 修改inp文件.......................................................................................................................... 47
11.用户子程序和FORTRAN ........................................................................................................... 49 12.创建/修改part ............................................................................................................................ 50 13.在Assembly中定位 .................................................................................................................... 51 14.划分网格 ..................................................................................................................................... 52 15.选择单元类型 .............................................................................................................................. 54 16.截面属性和材料 .......................................................................................................................... 56 16.1 定义截面属性 ......................................................................................................................... 56 16.2 用OFFSET参数偏置shell的中性面 ....................................................................................... 57 16.3 材料属性 ................................................................................................................................ 57
17.载荷 / 速度场 .............................................................................................................................. 58 17.1 集中力 .................................................................................................................................... 58 17.2 线载荷 .................................................................................................................................... 59 17.3 面载荷 .................................................................................................................................... 60 17.4 Body Force / Gravity .............................................................................................................. 60 17.5 随时间变化的载荷 / Amplitude .............................................................................................. 61 17.6 在实体单元上施加弯矩 / 旋转速度 ......................................................................................... 62
18.边界条件 ..................................................................................................................................... 67 19.定义Set / 选取模型的局部 .......................................................................................................... 67 20.弹簧 ............................................................................................................................................ 68 21.约束(coupling, tie) .................................................................................................................. 68 22.后处理 ......................................................................................................................................... 69 22.1 显示局部坐标系上的结果 ....................................................................................................... 69 22.2 绘制曲线(X–Y data) .......................................................................................................... 69 22.3 field output / history output ................................................................................................. 70 22.4 显示结果 ................................................................................................................................ 71 22.5 输出shell积分点上的应力和应变 .......................................................................................... 73
23.多步骤分析.................................................................................................................................. 74 23.1 重启动分析(Restart).......................................................................................................... 74
23.2 子模型(Submodel) ............................................................................................................ 74 23.3 用 *import 把计算结果传递到下一个过程中 .......................................................................... 74
24.岩土分析 ..................................................................................................................................... 75 24.1 初始地应力平衡 ..................................................................................................................... 75 24.2 基坑开挖 ................................................................................................................................ 77 24.3 钢筋混凝土(Rebar Layer和Embedded Element) .............................................................. 77 24.4 实体单元的截面力/弯矩/转角 ................................................................................................. 78 24.5 软土固结 ................................................................................................................................ 80
25.热分析 ......................................................................................................................................... 81 26.动态分析 ..................................................................................................................................... 83 26.1 Explicit分析 ........................................................................................................................... 83 26.2 分析结果与amplitude不一致 ................................................................................................ 84 26.3 固有频率 ................................................................................................................................ 84 26.4 *DYNAMIC分析 .................................................................................................................... 85
27.屈曲分析 ..................................................................................................................................... 86 28.裂纹 ............................................................................................................................................ 87 29.连接单元(CONNECTOR) ......................................................................................................... 88 30.惯性释放(INERTIA RELIEF) .................................................................................................... 88 31.非线性问题的算法 ....................................................................................................................... 89 32.ABAQUS和其它软件的比较 ....................................................................................................... 89 33.使用ABAQUS帮助手册 ............................................................................................................. 90 34.论坛的“全文搜索”功能 ............................................................................................................ 90 35.关于找枪手.................................................................................................................................. 91
1. 论坛提问要诀
在论坛上的提问是否能够得到解答,在很大程度上取决于提问的方式。提问时只是笼统地说“运行不了”、“不收敛”,大家就无从入手分析,再牛的高手也不能未卜先知。应该说清大致的操作过程,看到的现象,出现了什么提示信息。
Msg文件结尾的“ ERROR: Too Many Attempts Made for This Increment” 只是说无法收敛,而不收敛的原因则可能是多种多样的,往往需要看到模型才能找出具体原因。
应该尽可能把你的模型传上来,最好是ABAQUS/CAE模型,它比INP文件更容易看清楚边界条件,并且可以修改网格。论坛要求附件不得大于256K, 如果你的模型太大,可以考虑以下解决方法:
- 在ABAQUS/CAE中用save as, 换个文件名,文件会变小很多;
- 把网格删掉,然后save as, CAE模型就会很小,但应该保留你的网格种子,以便大家复原你的网格(不收敛的原因往往是网格太粗)。
- 申请一个email邮箱,上传你的模型,然后向大家公布邮箱密码。
如果你的问题解决了,希望你能抽空把你的解决方法简单写两句贴出来,对大家会有很大帮助。
2. 故障诊断基本方法
提交分析后,应该
1)察看dat文件,出现的warning一般可以不去理会,但如果出现error, 则分析无法进行;
2)如果分析能够进行,则察看msg文件(ABAQUS/Standard分析)或sel文件(ABAQUS/Explicit 分析)。如果msg文件出现error, 则分析同样无法进行;msg文件中的大部分warning都可以不去理会,但以下warning意味着模型有问题:
● Numerical Singularity (见下面3节);
● Zero Pivot (见下面??????????节);
● Too Many Attempts Made for This Increment:分析无法收敛,而不收敛的原因则可能是多种多样的,常见原因是接触或非线性材料定义不当,或网格太粗,需要具体问题具体分析
解决办法:
1)根据上述error或warning信息找出问题。
2)简化模型,去掉模型中的复杂参数(接触、塑性材料等),做一个简单模型,成功了再做复杂的。
3. 约束刚体位移 / NUMERICAL SINGULARITY
这是论坛上最常见的问题!!!
在静力分析中,必须对模型中所有实体的所有平移和转动自由度上定义足够的约束条件,以避免它们出现不确定的刚体位移。否则就会看到msg文件中出现
WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY
这时分析一般无法收敛。
下面是论坛上回帖的实例:
http://www.simwe.com/forum/viewthread.php?tid=419588
问:有警告说一个点自由度1,自由度2值都很大,但是检查发现这个点没问题,边界条件什么的也没问题,这怎么回事呢?
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE 15294 D.O.F. 2 RATIO = 2.48305E+11
答:\说明NODE 15294 所在的实体在方向2出现无限大的刚体位移。
可以在此实体上的任意一点和地面之间定义一个很软的弹簧,以消除刚体位移。方法是:interaction模块,菜单special / springs-dashpots / create, 选connect points to ground, 选节点,Degree of freedom 设为出现了刚体位移的自由度,spring stiffness为一个较小的值(太小则不足以消除刚体位移,太大则会影响变形)。
如果多个方向上出现了刚体位移,就要分别在相应的方向上各定义一个弹簧。
spring 所在的节点在弹簧方向的位移乘以spring stiffness,就是弹簧所分担的载荷,它应该远远小于在此方向上的外载荷。如果模型位移很小,我常常把spring stiffness设为1.
http://www.simwe.com/forum/viewthread.php?tid=686203
问:我做了一椭圆物接触平面的模型,接触受力后,椭圆物(约束为固定)平移了,不知是啥原因
答:椭圆没平移啊,只看到平板平移了。平板在u1上没有边界条件,接触面又不是水平的,当然受力后会平移的。
另外,椭圆顶部网格形状太差,接触区域的网格要细化。
这是基本概念啊,静态分析中,每个实体的每个自由度上都要有足够的约束,否则它会有不确定的刚体位移或旋转.
每次计算时,都应该看msg文件中有没有warning或error。你的msg文件中有
warning: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE FLAT-1.14 D.O.F. 1 RATIO = 2.50278E+010.
说明 FLAT-1的D.O.F. 1 上出现了无限大的刚体位移。消除的方式可以是边界条件,约束,接触,摩擦,弹簧等。你如果不希望平板在x方向运动,需要施加边界条件u1=0.
对称问题,对1/2建模,约束整个对称面上的U1. 不对称的话,也要想其他办法施加足够的约束。
如果只有负特征值警告,没有numberical singularity, 计算能收敛,就没问题。
http://www.simwe.com/forum/viewthread.php?tid=690585&highlight=?ùò??a???£Díà?μ?????2??t1?
答:已经说得没法再明白了。你加位移载荷的模型spring_contact.inp里:
- Flat-1: U1=U2=0
- Flat-2: U1=0, U2受弹簧约束 - 椭圆: U1=0, U2有给定的位移
所以这个模型里的每个部件关于x,y轴的刚体平移和关于z轴的刚体转动都得到了足够的约束。
你加力载荷的模型sphere_flat.inp里:
- 椭圆: U1=U2=0, 没问题。
- 平板: y方向受力,上面和椭圆接触,所以y方向也没问题。但x方向没有施加任何约束,尽管这个方向没有力,对你来说平板不会发生x方向的刚体平移,但这不符合有限元的计算要求,x方向的刚体位移是不确定的。这并不是说,你会在后处理中看到x方向发生很大的刚体位移,而是ABAQUS根本就无法计算,得不到任何结果。这和y方向的载荷大小没有关系。
所以要根据对称性对1/2建模,约束整个对称面上的U1.
不用看job diagnostics ,错误和警告信息都在分析生成的文件sphere_flat.msg里。
你要是还不明白,我就只能去跳楼了。
http://www.simwe.com/forum/viewthread.php?tid=695202
问:我用constrain_equation 将梁和板系在一起,进行均布力受压分析,但是发现,当我将1,3自由度限制时,板梁都弯曲了,2方向有偏移
但是我将1,2,3自由度一起限制时,发现板梁都不弯了,可能原理上有点问题,附上文件,哪位大虾给与指点!!谢谢
后来用tie 试了一下,效果还挺好的,就是不知道对不对,因为帮助文件里说tie用于定义面和面的,梁和面这样定义不知道行不??
答:每次计算时,都应该看msg文件中有没有warning或error. 两个模型的在msg文件中都有
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE 104 D.O.F. 4 RATIO = 3.80284E+015.
说明 D.O.F. 4 上缺乏约束。你只约束U1,U2,U3, 对实体单元是可以的,但shell和beam单元都有旋转自由度,要施加足够的约束来消除刚体位移 (虽然对你来说,约束了两端的U1,U2,U3就没有刚体位移了,但这不符合ABAQUS的要求)。你的第一个模型虽然算出来了,但结果不一定正确。
解决方法:对称模型可取一半,在对称面上约束UR1,UR2,UR3.
梁和面tie应该可以。
4.过约束(Overconstraint)
http://www.simwe.com/forum/viewthread.php?tid=718246
问:做了个基坑开挖的题,出现了很多过约束问题,怎么解决过约束的问题呀?
答:参见
http://www.simwe.com/forum/post/view?bid=31&id=705685&tpg=1&ppg=1&sty=1#705754
接触对的主面上不能有尖角,桩的两个侧面要分别定义接触对,底部可能可以用tie.
slave surface的网格要比master surface细。
过约束可能是因为被挖的土上既定义了接触,又要被杀死,这二者相矛盾。可以试试为每段被挖的土单独定义一个接触,挖土时先deactivate这个接触,再杀死单元。
http://www.simwe.com/forum/viewthread.php?tid=724602&pid=738324&page=1&extra=page=1#pid738324
问:
在job-monitor中有个警告是这样
Solver problem. Zero pivot when processing D.O.F. 2 of 1 nodes. The nodes have been identified in node set WarnNodeSolvProbZeroPiv_2_1_1_1_1.
是不是刚度矩阵的问题?我的程序是照着一个讲座的材料上抄下来的,应该没什么问题的哪位老兄能指点下!
SimWe仿真论坛5o-p:td/k+Kò[)T
仿真分析有限元模拟
答:
Zero pivot 往往意味着 OVERCONSTRAINT。此警告信息如果只是出现在dat文件中,没有出现在msg文件中,就没问题,说明ABAQUS自动解决了过约束问题。
如果overconstraint警告信息也出现在msg文件中,说明ABAQUS无法自动解决此问题,这时分析往往不会收敛,在后处理时可以用display group显示出现过约束的node set WarnNodeSolvProbZeroPiv_2_1_1_1_1. 这时需要你自己修改模型,避免过约束。
5.其他常见error和warning信息
http://www.simwe.com/forum/viewthread.php?tid=692452
问:There is zero MOMENT everywhere in the model based on the default criterion. please check the value of the average MOMENT during the current iteration to verify that the MOMENT is small enough to be treated as zero. if not, please use the solution controls to reset the criterion for zero MOMENT.
该怎么办啊,谢谢大家了!
答:这只是提醒你模型中只有位移边界条件,没有施加载荷。如果能收敛,结果也正常,就没关系。很多警告都只是提示信息,不意味着模型有错误。
http://www.simwe.com/forum/viewthread.php?tid=683255
问:1 ***WARNING: DEGREE OF FREEDOM 4 IS NOT ACTIVE IN THIS MODEL AND CAN NOT BE RESTRAINED
2 ***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT: ANALYSIS TERMINATED
3 ***ERROR: NODE SET ASSEMBLY_CENTER HAS NOT BEEN DEFINED
4 ***ERROR: A BOUNDARY CONDITION HAS BEEN SPECIFIED ON NODE SET ASSEMBLY_CENTER BUT THIS NODE SET IS NOT ACTIVE IN THE MODEL
对于以上错误,哪位大侠指点一下,是哪里出现了错误,如何解决它, 在此表示感谢。
答:
1.―DEGREE OF FREEDOM 4 IS NOT ACTIVE‖:很多警告信息都不意味着有错误,可以不管它。
2.―TOO MANY ATTEMPTS‖:如果以当前的时间增量步无法在规定的迭代次数内达到收敛,ABAQUS会自动减小时间增量步,重新开始迭代,即所谓的―Cutback‖。如果这样仍不能收敛,则会继续减小时间增量步。如果达到了规定的Cutback最大次数(默认值为5次),或时间增量步长减小到所规定的最小限度(默认值为10-5),ABAQUS就会中止分析,并在MSG文件的结尾处显示上述错误信息。
不收敛的可能原因很多,就像天上的星星,数也数不清。所谓 ―正确的模型都是相似的,不正确的模型则各有各的不正确‖。
3.―NODE SET ASSEMBLY_CENTER HAS NOT BEEN DEFINED‖:模型中用到了名为center的集合,但没有定义它。
4.―***ERROR: A BOUNDARY CONDITION HAS BEEN SPECIFIED ON NODE SET
ASSEMBLY_CENTER。BUT THIS NODE SET IS NOT ACTIVE IN THE MODEL‖:同上。
运行INP文件时会调用子程序,二者不能集成在一起。
http://www.simwe.com/forum/viewthread.php?tid=690862
问:大家好,我是一个AQBAQUS的初学者,最近在计算软土固结问题时,总是出现 FIXED TIME INCREMENT IS TOO LARGE
这样的问题,屡试屡败,所以非常想知道如何才能正确设置time increment选项及其中的窍门,不知道哪位高手肯指导一二,本人不胜感激。
答:岩土我没搞过,不一定是increment的问题,你检查一下单元类型,材料参数,单位,初始地应力等是否正确。建议在帮助文件或论坛里找一个正确的例子,把它一点点改成你的参数,改到哪里出错了,就找到了问题所在。
http://www.simwe.com/forum/viewthread.php?tid=697404
问:我在运行abaqus时出现错误:
abaqus error:the executable C:\\abaqus\\6.5-1\\exec\\standard.exe aborted with system error \拒绝访问\
答:先看dat和msg文件中有没有warning或error,如果没有,常见的有三种原因:
- abaqus_v6.env中的参数pre_memory设置得过大
- 用户子程序中出现问题。检查在INP文件中调用用户子程序时,是否使用了正确的参数。可以尝试用一个简单的例子来调试用户子程序,确保其正确性。另外在编写用户子程序时,应注意尽可能多设置一些出错提示。
- ABAQUS本身的bug。可以尝试去掉模型中不必要的特殊设置,各参数尽量使用ABAQUS的默认值。如果初始的模型没问题,而在修改了某些参数之后出现上述错误信息,可以尝试逐个将这些参数恢复至初始状态。
在本版搜索\拒绝访问\可找到同类问题的解答。
(自己的经验:)
问: *contact pair, interaction=INTPROP-NOFRICTION, type=NODE TO SURFACE, no thickness, small sliding, adjust=0.02
错误信息:
***ERROR: THE AVERAGE NORMAL VECTOR CANNOT BE COMPUTED AT NODE 493 INSTANCE ERSATZRING-1 ON THE SURFACE ASSEMBLY_SURF-ERSATZ-AXIAL-MIT-NABE. PLEASE CHECK THE SURFACE DEFINITION OR USE THE *NORMAL KEYWORD TO SPECIFY THE NORMAL VECTOR.
答: 圆角处种子太稀疏,接触面的slave surface被adjust之后造成畸形单元。
http://www.simwe.com/forum/viewthread.php?tid=728204
问:我又重新写了一次input文件,这次没有参数,可以再cae中直观地看到图形,但运行后,在.dat文件中总是出现这个错
误:
***ERROR: The area of 54 elements is zero, small, or negative. Check
coordinates or node numbering, or modify the mesh seed. The elements have been identified in element set ErrElemAreaSmallNegZero.
***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO INTERPRET SOME DATA. SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent-t6|#g*Eo*kZE,c A
www.simwe.comh'?0kQ0HQ*[5D!c
SimWe仿真论坛s4P'[3N&O|Simwe.com|仿真设计有限元虚拟仪器w1v2w%|Simwe.com|仿真设计
答:
你可能是单元的节点编号不对,应该是逆时针方向。
在后处理中用display group可以看到上述错误信息提到的集合ErrElemAreaSmallNegZero. 在assembly模块中菜单view / Assembly Display Options / Mesh,可以显示单元和节点编号。菜单tools / query / element 可以查询单元编号。
http://www.simwe.com/forum/viewthread.php?tid=703681 问:
i-w I%K I
我是初学者,用abaqus计算用户手册带的例子.
www.simwe.com~;i版本是6.5.1的.出现下列问题: ***ERROR: Issue cannot be deleted Not all data Released
答:根据你的错误信息 ―Issue cannot be deleted Not all data Released‖,在网上搜了一下,找到一个德文版的解答:
okay, wurde gekl?rt mittlererweile. In Systemsteuerung->System->Erweitert->Systemleistung->Datenausführungsverhinderung den pre und den standard hinzufügen, und voila, l?uft wie geschmiert.
大意是:windows控制面板 / 系统 / 高级 / 系统性能 / 数据运行连接 (我这里没有中文系统,不知中文windows下用的是哪个词),把 pre.exe 和 standard.exe 添加进去。
不过我在windows xp下没找到这个―数据实现连接‖。
http://www.simwe.com/forum/viewthread.php?tid=704712
问:在abaqus中可以将摩尔库伦模型中的凝聚力和摩擦角定义为塑性应变的函数,从而模拟应变软化。请问:有人用过这样做过?请给点经验。
本人在计算过程中发现很那收敛,常出现两个典型的错误提示:
(1) time incrementation is exceeded and (2) convergence is not possible.
答:先不要用这些复杂的设置,用最普通的参数算,收敛了再一点点改成复杂的。
http://www.simwe.com/forum/viewthread.php?tid=718464
问:在Assembly中通过Instance Part命令将Punch,Specimen和Die分别添加到组合图。发现刚性体和柔性体不能同时装进去。老是出现:part dimensionality does not match the dimensonalty of the existing instances。
答:create各个part时,modeling space要相同。
http://www.simwe.com/forum/viewthread.php?tid=716573
问:墙的重力加不上去,提示有mixed dimentionaries?老大怎么办?
答:墙的单元是CPS4R, 应该是 CPE4R. 土的单元CPE4精度很差,要用CPE4R。
杆应该改成梁吧。
另外,dat文件中有warning:
***WARNING: IN THIS TYPE OF ANALYSIS ONE WOULD NORMALLY EXPECT THAT A DENSITY FOR MATERIAL CONCRETE2 WOULD BE DEFINED.
6.ABAQUS无法运行(安装、License、环境变量)
http://www.simwe.com/forum/viewthread.php?tid=681671
问:今天换了个杀毒软件和防火墙,运行cae时,出现:cannot connect to license sever~~~~,请问为什么?!!
答:ABAQUS的多个文件都要访问网络,这时防火墙一般会弹出提示,应选择 ―允许此文件访问网络‖。可以检查防火墙的设置,不要让它屏蔽那些ABAQUS安装目录下的文件。实在不行就关掉防火墙。
http://www.simwe.com/forum/viewthread.php?tid=685885
问:我是按照例子的步骤做下来的,最后submit时,打开monitor监控 monitor中显示如下信息
The executable d:\\abaqus\\6.4-1\\exec\\pre.exeaborted with system error code 29539. 不知道何故? 请高手解答一下! 3x!!
答:在本版搜索 error code 29539,可以找到下面的帖子:
http://www.simwe.com/forum/viewthread.php?tid=151256
我的abaqus安装的时候全都通过了,可是在我做练习的时候,到job模块的时候我create job后submit的时候就出现下面的问题了,
ABAQUS Error: The executable g:\\ABAQUS\\6.4-1\\exec\\pre.exe aborted with system error code 29539.
但是我在command里面运行abaqus j=deform还是可以的,msg文件中也没有错误,结果用view也可以看。
以上是我这几天安装abaqus以来遇到的一个一直没有解决的问题,问了无数的人都没有结果,今天终于凭借我自己的脑袋解决了这个问题。
其实很简单,原因就在计算机名上, 我原来的计算机名子是6-33
就是这个名子惹的事,后来我改成全是字母的就搞定了,我相信把中间的短横去掉就可以了,只要没有短横就可以了。
哈哈,想知道我是怎么想起来可能是这个原因的吗,因为我以前用CFX的时候就遇到过相似的问题,今天突然想到有可能是这个原因就试了一下,没想到就通过了,这回斑主是不是应该给我加分了呀。
http://www.simwe.com/forum/viewthread.php?tid=687097
问:我安装一路均没问题但最后在verification时其它均pass,但老是ACE failed,运行ACE时出现的错误信息如下
cannot find server hostname in network database the lookup for the hostname on the server line in the
license file failed this ofen happens when nis or dns or the hosts file is incorrect. workaround:use Use IP-Address (e.g., 123.456.789.123) instead of hostname feature:cae-teaching
Hostname:
License path: 27000@USER-A6E7EA6B9E.;27000@
FLEXlm error: -14,7. System Error: 11004 \but no record
abaqus error:abaqus/cae kernel exited with an error
上述是什么意思 我该怎么做才能正确安装呢 我都困惑了几天了 请各位abaqus高人多多指点 谢谢 谢谢
答:在本版搜索\,有很多解决方法的提示.
你的计算机名是不是 ―USER-A6E7EA6B9E. ‖?或许是结尾处的点或者中间的横线导致的问题。
关掉防火墙。
http://www.simwe.com/forum/thread-727593-1-2.html
问:我的cae只有把显卡驱动卸了才能用,可能是两者冲突,怎么设置 exceed才能解决? 我的显卡是nvidia的。 hanma答:打开你的安装目录,site文件夹里 找到abaqus.env的文件,或者是*env具体名字我忘了,它可以用文本文件打开,然后托到最低下,大概倒数第四行左右有一个选项, hardwareaccelaration=hardware.......之类的一个参数
让他=OFF,注意一定要大写,就是把硬件加速关掉,就可以正常用了(6.4,6。5没有这个问题),显卡----不成问题!
|Simwe.com仿真设计有限元虚拟仪器-k$P)YD9^A-b:}www.simwe.com%h4Z,\\-a#e+v p0Q
http://www.simwe.com/forum/viewthread.php?tid=709422
问:我的ABAQUS原先安装好了,也可以用,可后来不知道什么原因出现问题。现象是建模画图时除了画点以外的任何操作都会使程序没有响应,重装过几次都不行。实在没办法了,希望各位帮忙解决一下。(注明:原先是可以用的)
答:在本版搜索\硬件加速“。
在本版搜索\显卡“。
http://www.simwe.com/forum/viewthread.php?tid=697681
问:小弟在用ABAQUS的GUI模式建模时发现每次在LOAD步定义约束时,abaqus就突然跳出,重复了好几次都是如此,正不知道是怎么回事?请各位大侠指点一下,谢了先, <98 out of 99 licenses remain available>.
GUI detected error while waiting for ipc connection to close. ABAQUS Error: ABAQUS/CAE Kernel exited with an error. 请按任意键继续. . .
答:在本版搜索 “GUI detected error”
http://www.simwe.com/forum/viewthread.php?tid=698879
问:求救!document无法搜索
答:要先安装document的cd, 再安装ABAQUS软件的CD。
在本版查找 ―documentation 搜索‖,以及 ―帮助 搜索‖。
在abaqus_v6.env中把
doc_root=\
改为 doc_root=http://winr:2080.v6.5
http://www.simwe.com/forum/viewthread.php?tid=681116
磁盘空间不足
问:我在进行接触分析遇到了一个问题,错误提示为:ERROR: UNABLE TO COMPLETE FILE WRITE. CHECK THAT SUFFICIENT DISK SPACE IS
AVAILABLE. FILE IN USE AT FAILURE IS subcasejob-1_metsp.uft.
可我计算机C盘剩余8G,D盘45G,abaqus安装在D盘,内存为8G 单元数量为26万体单元
可为什么还说磁盘空间不够呢。请指教
答:分析过程中除了在当前工作目录下生成odb, msg, dat等文件外,还会在默认的临时文件目录下生成很大的临时文件(此目录好像是c:/temp),分析中止后,这些临时文件会被自动删掉,所以尽管你看到c:有很大空间,但在分析过程中可能会空间不够。
你的d盘很大,可以把临时文件的目录修改到d:, 方法是:
在ABAQUS安装目录下的文件夹site中可以找到环境文件abaqus_v6.env,在其开头添加
scratch='D:\\scratch'
保存此文件,在D盘创建目录scratch. 重新启动ABAQUS/CAE.
另外,在abaqus_v6.env中还可以看到关于内存的设置:
# 对INP文件进行预处理所允许使用的最大内存 pre_memory = \#
# 使用ABAQUS/Standard求解时所允许使用的最大内存 standard_memory = \
可以把这两个值改大些(但不要超过你计算机内存的实际大小).
melxc2003:察看你的临时文件所在的硬盘,如果硬盘是fat32类型的,那么存放文件不能大于4G,否则会抱错。
Jacobi:这么大的模型,建议用linux version
http://www.simwe.com/forum/viewthread.php?tid=701624
问:abaqus在linux系统上要比winxp节省很多时间…………不知有没有哪位也试验过?
答:linux要比winxp的分析速度快15%, 且linux能够利用的内存上限要比winxp大。
磁盘阵列没用过。如果分析时内存的使用量超出了物理内存的大小,买一条内存肯定是有用的。
http://www.simwe.com/forum/viewthread.php?tid=719091
问:电脑自动关机,abaqus自动保存了文件,以前,都询问是否恢复。这次,cae不能读了,多了一个rec文件,错误提示: cae file can[t be opened: the database is corrupted. ren abaqus cae recover=
在command里运行abaqus cae recover=
答:好像没什么办法。自己要注意经常保存,关掉CAE才是真正保存了模型。
7.弹塑性分析
7.1 塑性问题不收敛的常见现象
塑性问题不收敛时,msg文件中的常见现象是
1) 出现很多EQUILIBRIUM ITERATION,且TIME INCREMENT 不断减小;
2) 始终出现 ***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
解决方法见下面回帖实例。
7.2 接触问题和塑性材料不要用二阶单元
http://www.simwe.com/forum/viewthread.php?tid=647602&page=1&extra=#pid729426
问:aba_aba老大,你说
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent+j4dCd-\\ vI2Fh用C3D15单元,是15-node quadratic triangular prism, 这还是一种二阶单元,不能用于塑性材料和接触问题。这个是什么原因?
www.simwe.com)Te)E#|*@M|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent(X&@(p/G$?:J0w
答:
-- 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.
7.3 不要在塑性材料上施加点载荷
http://www.simwe.com/forum/viewthread.php?tid=725926
问:在我的模拟一个桶土接触问题中,我想加一个和坐标轴成一定角度的集中力,但是直接加集中力会导致不收敛,因此我改用
加位移荷载的方法,求出反力,再反推合力,以及合力与坐标轴的夹角,请问这个合力方向是不是与所加的合位移方向相同,本人认为是相同,不知道大家有什其他的看法,请教一个?
答:有些情况下一样,有些情况下不一样,你试试就知道一样不一样了。施加位移荷载时的反力就等于力载荷。
“直接集中力会导致不收敛”:
一是加载之前要有一个只施加很小载荷的分析步,使各个接触关系都建立起来。
二是如果受集中力的材料是塑性,会因为局部变形过大而不收敛,应该将力施加在一个小区域上,让多个节点共同受力。
7.4 塑性材料参数
http://www.simwe.com/forum/viewthread.php?tid=717593
问:
当设定stress 和 strain时总不能计算出结果 只有elastic的时候就没问题 谢谢了
当设定时注意什么
答:
要让*Plastic最后一行的数据大于模型中可能出现的Mises应力和PEEQ. 见附件。
载荷不能太大,改为-1500。
http://www.simwe.com/forum/viewthread.php?tid=729810
问:二维模型为鋼板接頭實驗,因为是对称的,所以只取其中一半。上、下圓為螺栓(設定為rigid body),Step 1:先讓上;
下螺栓接觸鋼板,並施加一微小壓力避免接觸時發生表面震顫。Step 2:下方螺栓固定不動,上方螺栓向上位移直到位移結束。
我的問題是:
(1)我下方圓孔接觸的地方附近再割密一點,就會發生‖too many attemps made for this increment‖,是不是指元素扭曲過大,收斂困難???(因為我把網格切大一點又就可以跑,我有試過aba-aba在別人的回文中的方法,將菜单other / Gerneral solution controls / Edit,把Is由默认的12改为30,可結果還是一樣。
(2)另外,因為我感興趣的部分是在位移過程中(step 2)我想知道鋼板所受的外力,因此我在上方螺栓作reaction force的輸出(REF-1),但是我發現分析出reaction force來的結果其前面那一部份還蠻符合,但到了後面塑性的部分,力就開始發生不規則的震盪,不知道是哪裡有錯,
(3)step 2裡,為什麼在一開始位移時卻沒有受力,是因為mesh不夠密嗎(平滑)?,導致一開始並沒有完全接觸到,所以沒受力。
|Simwe.com|仿真设计有限元虚拟仪器3L%Tn.LE:T/S
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1?~z7OD ~
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMe&h5E S$N
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM_Ii-}針對鋼板下方螺栓圓孔的上部頂點位移與上方螺栓的反力,畫出load-displacement關係圖,與真實實驗結果作比較,供大家參考。希望知道問題的人,能不吝嗇指教我。
謝謝了!!
SimWe仿真论坛e&f-U
www.simwe.comJ3s附件: steel-plate.rar (2006-6-9 20:35, 34.91 K)
答:
msg文件中有
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE STEEL PLATE-1.52 D.O.F. 2 RATIO = 5.51056E+014.
意思是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很难找到收敛的位置, 应该为一条光滑的曲线; - 给出的塑性数据点太多,好像也会造成收敛困难,建议改为
*Plastic
343.964, 0. 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, 0.38856 470., 1.
这样就收敛了(没有最后一行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*1000=12000, 而支反力最大值为52070,与弹簧分担的载荷在同一个数量级上。弹簧的作用只是帮助收敛,不应该影响各个部件的受力和变形。可以试试减小弹簧刚度,例如刚度为1,不收敛的话再增大一些。
http://www.simwe.com/forum/viewthread.php?tid=714029
问:
请教各位大侠这是怎么回事:
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 *Elastic 225., 0.3 *Plastic 0.35,0. 10.0, 100.0 **
*Material, name=Material-2 *Elastic 170., 0.3 *Plastic 0.158, 0.
0.22838, 0.0016 0.45333, 0.02239 0.57685, 0.04982 0.63916, 0.06956 0.69071, 0.08936 0.72243, 0.10325 0.77079, 0.1271 0.8234, 0.15695 0.86677, 0.18485 0.88119, 0.19482 10.0, 100.0
把缩放系数设为1,可以看到单元扭曲太大,所以无法收敛。要在大变形部位细化网格,变形小的部位的网格可以很粗。
http://www.simwe.com/forum/thread-728444-1-1.html
问:桩土接触(群桩) 塑性为什么计算不出结果?(弹性可以)
本人一直在研究群桩的计算 今天计算了一个七根桩的模型(三维)。 模型具体为:
桩:由圆桩(直径2.8m)等效为八面体桩 桩长100m 全部入土 梅花形排列 土:共七层土体 有五层和桩接触 土体40*40*200m
接触:每根桩和五层土体接触(共7*5=35个接触对 node based surface),桩底采用tie约束 共27000左右个单元(c3d8r)
我分别用D-P模型和M-C模型都无法算出结果,后来我改为弹性可以算了而且结果还可以 并且平衡了地应力
讲了这么多 我就是想知道 为什么材料为塑性就算不出结果(经常是单元变形过大 材料点发生屈服) 是我的塑性参数取的不合理吗?还是别的什么原因?
一般土体yield stress取多少?第一个step 平衡地应力可以用塑性材料吗?如果用弹性材料 后面计算怎么该为塑性材料?请各位给点建议
www.simwe.com#Al%o
www.simwe.com$_ a@*\\*x3\\Y0r'C
SimWe仿真论坛{:WF&D4| I|(x7\\ X
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/_K7pm%l(_!v\
lyj198171答:我是把yield stress 该大了 一直不清楚这个值到底该怎么取 m-c取为c的话 ,土体也就几十Kpa 是不是太小了? 有没有人知道具体应该怎么取 ? 我200m的土 如果yield stress为30Kpa那么在自重作用下不是已经屈服了吗?
http://www.simwe.com/forum/viewthread.php?tid=717089
问:只在四节点施加一竖向位移,而水平方向怎么会出现应变,且同竖向的应变在同一量级?把弹性单元编为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
7.5 屈服准则
http://www.simwe.com/forum/viewthread.php?tid=688718
问:在塑性加工中,不同的材料用如下变形抗力模型表示。请问:这个公式指的是整个弹塑性过程还是忽略弹性的?如果包括弹性,那么弹性模量又与该公式有何关系?屈服极限应力在该式中有体现吗?该公式反应的是真实应力应变关系吗?它的适用范围是什么?
不胜感谢!!!!!!
wy_anli答:塑性加工中的具体材料模型有成千上万个,大类可以分为理想刚塑性(忽略弹性、忽略硬化,主要应用于弹性变形相对于塑性变形较小,可以忽略的情况下,如:锻造等)、理想弹塑性(忽略硬化)、弹塑性(在材料成形时其弹性变形不可忽略,例如板料成形)、等等。 你的这个模型属于弹塑性模型,没有忽略弹性。
屈服极限是材料发生塑性变形时的临界应力,是材料的本身性质,在本公式中不能体现。
纵向应力与纵向应变的比例常数就是材料的弹性模量,也叫杨氏模量。也不能在此公式中反映出来。 任何材料模型都不可能完全反映真实的应力应变关系,当然次模型也不例外,只是最大近似。 一般只有根据材料提出材料模型,其适应范围也只能适用于特定材料以及特定场合。
wpengw123 wrote:
非常感谢。还是有些不解,就直接问了,不吝赐教!
那么该公式是否就是整个弹塑性应力应变关系的简化式呢?还是只描述塑性变形阶段?
是塑性应力应变关系的简化式,但还缺少抗拉强度或极限应变(当材料达到这个临界状态时就会发生断裂,继续算下去,就没有了意义)。弹性阶段在塑性阶段前考虑。
wpengw123 wrote:
如果包含弹性,那么杨氏模量是不是和公式有冲突呢?
杨式模量不会和公式有冲突,杨式模量是材料本身的参数,是材料在弹性阶段时应力应变线性关系的系数。Abqus中定义材料时,弹性和塑性是分开来的。弹性时只有杨式模量和波松比两个参数。 wpengw123 wrote:
而且它似乎常用于金属板料成形,优点是什么呢?
板料成形要考虑弹性变形,因为板料成形是大变形小应变情况,弹性应变不能忽略不计,并且在塑性成型结束后弹性变形要完全恢复(也就是回弹)。相比较来说,例如锻造等成形,弹性应变相对于塑性应变非常小,可以忽略不计,也就是要采用刚塑性材料模型。 wpengw123 wrote:
式中常数和变量的具体意义是什么呢?
参见塑性成形原理。
wpengw123 wrote:
怎样和单向拉深的实验数据联系,又怎样在Abaqus中使用呢?
单向拉深的实验数据通常是名义应力和名义应变。要通过名义值和真实值之间的关系式将其转化为真实值也就是等效应力和等效应变。Abaqus中的材料采用真实应力和真实应变。
wpengw123 wrote:
看到很多文章在模拟板料塑性成形时用到此模型,但是在给出杨氏模量及强度系数和硬化指数的取值时却都没有给出屈服极限,屈服极限不需要吗?
屈服极限也就是材料在屈服点上的应力应变状况。材料模型中应变为0时的应力就是屈服极限,大多数金属的屈服应力为材料弹性模量的0.05%-0.1%。
http://www.simwe.com/forum/viewthread.php?tid=724557
qyuan_13: 关于对球墨铸铁的强度评价标准
通常情况下,铸铁的脆性比较强,延伸率很低,一般以延伸率为0.2%作为屈服极限。对于脆性材料作强度分析的时候,采用主应力来评价;但对于球墨铸铁,有比较高的延伸率,比如QT450-5、QT450-10,其延伸率分别为5%和10%,那么,对这样的材料作强度分析的时候,是不是应该采用Mises应力来评价?按说10%的延伸率比较高了,合金钢20Cr的延伸率也为10%,40Cr的仅为9%(刘鸿文的材料力学)。那么,在脆性材料和塑性材料之间的界限到底是多少呢?不能仅仅认为只要是铸铁,就得用主应力来评价吧!
http://www.simwe.com/forum/viewthread.php?tid=702428
问:在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。
7.6 后处理 / 判断材料是否屈服
http://www.simwe.com/forum/viewthread.php?tid=691447
问:DOCUMENT上说PEEQ是等效塑性应变(Equivalent plastic strain),PEMAG是塑性应变值(Plastic strain magintude),感觉前者主要是判断材料是否屈服,而后者呢??
答:在比例加载时(即加载过程中主应力方向和比值不变),大多数材料的PEMAG和PEEQ相等。这两个量的区别在于,PEMAG描述的是变形过程中某一时刻的塑性应变,与加载历史无关,而PEEQ是整个变形过程中塑性应变的累积结果。例如,单向拉伸一个圆柱体,使其发生塑性应变,再通过单向压缩使其恢复初始长度,则最终的PEMAG为0,而PEEQ是拉伸和压缩过程中塑性应变之和。
http://www.simwe.com/forum/viewthread.php?tid=566054
问: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类似于各个应力分量。
http://www.simwe.com/forum/viewthread.php?tid=696024
问:比如计算地基的极限承载力,m-c理想弹塑性,用计算结果不收敛作为破坏评判准则。
现在我做的时候,荷载是手动调整,每次调整完再提交计算,感觉比较麻烦
能否实现这样的功能:给出初始荷载以及增量步长,要是能自动调整就更好了。每一级计算完成后保存结果,并增加荷载继续计算,直至不收敛
ansys好像有循环语句可以实现,不知道在abaqus中如何实现。谢谢了~
答:不需要特别的设置,直接创建一个step, 设置initial increment(例如0.1), 在这个step里施加一个载荷(比如1000N), 于是step time = 0.1时,载荷就是100N; step time = 0.2时,载荷就是200N;超过极限载荷时,计算不收敛,ABAQUS会自动减小increment。
问:to 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 有解释:
\
MPa, which should not be possible since the material was assumed to be perfectly plastic at this stress magnitude. \
\material; 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 current yield stress; only the Mises stress is required to have a magnitude less than or equal to the value of the current yield stress.\
另外,帮助中的这个例子还说了,理想塑性不容易收敛,最好加入硬化特性,让材料曲线中最大的应力应变大于模型中可能出现的应力应变。
8. 接触分析
8.1 接触分析不收敛的常见现象和解决方法
做接触分析时,应观察msg文件中的以下信息
(接触状态的变化:1个节点闭合,10个节点开放)
CONTACT CHANGE SUMMARY: 1 CLOSURES 10 OPENINGS.
SEVERE DISCONTINUITY ITERATION 2 ENDS
CONTACT CHANGE SUMMARY: 0 CLOSURES 4 OPENINGS.
如果分析能够收敛,CLOSURES和OPENINGS的数目会逐渐减少,否则一般是接触定义不当,可以考虑以下解决方法:
● 检查所定义的接触面、接触参数和边界条件是否正确;
● 在静力分析中,必须定义足够的约束条件,以消除各个平移和转动自由度上的刚体位移; ● 避免过约束(overconstraint);
● 合理地定义接触面、接触参数和过盈配合; ● 使用足够细化的网格;
● 不要在接触面上使用C3D20、C3D20R和C3D10等单元;
● 在接触对上设置微小的过盈量,以保证在分析的一开始就已经建立起接触关系; ● 施加临时边界条件,以保证在接触关系建立之前,模型也不会出现刚体位移;
● 将分析过程分解为多个分析步来完成,让各个载荷分别在不同的分析步中逐步施加到模型上,避免使接触状态发生剧烈的改变。
下面按照应用领域对回帖分了类,但其实定义接触的原则是相同的,例如岩土类问题的接触问题的解决方法,往往同样适用于金属成型中的接触问题。
8.2 接触面上的网格密度
http://www.simwe.com/forum/viewthread.php?tid=723292 问: 如何得知接触表面网格太密或太疏而造成不收敛的问题呢..
www.simwe.comu3K S*Z2T
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent-U P vB q-A$kW
掌握 master 面网格密度比slave面网格密度疏,此点足够吗?
接触分析中若有大变形情形,因为很不容易收敛
改变 increment size 中的 minimum 值,将其值设小(default:1E-005)
SimWe仿真论坛7Oz8Q'[;v6S
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM\{.I2j)yw
SimWe仿真论坛7f)YN:Ch&NN*X
能提高分析成功的机率吗?
答:网格疏密要靠经验,比如接触面在某个方向上只有一两个单元,就是太疏了,如果接触面在某个方向上有上百个单元,就可能是太密了。
大变形时应该减小increment size 中的 minimum 值,但如果接触或塑性材料参数定义得不正确,还是不能收敛。
8.3 接触面的法线方向
http://www.simwe.com/forum/viewthread.php?tid=707677
问:任务提交后提示以下错误:―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 SEVERE.
这往往是因为接触面的法线方向定义反了。定义刚体和shell的surface时,要注意选择外侧。在后处理中可以显示法线方向:ABAQUS/CAE User's Manual 37.7 Displaying element and surface normals。
8.4 过盈接触
http://www.simwe.com/forum/viewthread.php?tid=672658
问:请问, 过盈配合怎么做, 用哪个菜单?
答:---------------------------------------- *CONTACT INTERFERENCE
---------------------------------------- 此关键词的使用方法为:
*CONTACT PAIR, INTERACTION = <接触属性的名称 >, ADJUST = <位置误差限度 > <从面名称 >,<主面名称 > ……
*AMPLITUDE,NAME = <幅值曲线的名称 > …… *STEP ……
*CONTACT INTERFERENCE,AMPLITUDE = <幅值曲线的名称 > <从面名称 >,<主面名称 >,<过盈量或间隙量 > ……
*END STEP
其中,参数<过盈量或间隙量 >为负值表示过盈量,正值表示间隙量。使用 *CONTACT 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 中的参数 <过盈量或间隙量 >
不会起作用。
● 参数<过盈量或间隙量 >为负值才表示过盈接触,正值表示两个面之间存在缝隙。
http://www.simwe.com/forum/viewthread.php?tid=692864
问:使用clearance和使用CONTACT INTERFERENCE算出来的结果有什么区别么?
答:*CLEARANCE也可以,但有一些缺点,它只能用于small sliding,不能在ABAQUS/CAE中定义,不能在step中修改过盈量,而且没有使用ramp, 也就是说,在分析的一开始全部过盈量就会被施加在模型上,如果此过盈量太大,ABAQUS就可能无法达到收敛,而且即使减小时间增量步也没用。
http://www.simwe.com/forum/viewthread.php?tid=687881
问:我的两个接触面是经过运动后才到达相互接触的状态的,但这时并不能达到理想的接触(两个面有部分gap,很小但的确没有接触,用contact print 打印出来slave nodes相关项,都是zero)。我想在这个时候来adjust the contact gap,让abaqus把接触面间的gap调成没有。
因为一开始时两个接触面的间距非常大,所以不可能在一开始就调整,而且我还需要对那个运动进行分析。
答:使用 *CONTACT INTERFERENCE可以在分析步中定义接触对之间的距离或过盈量,你可以试试。
使用 *CONTACT INTERFERENCE时要注意:
1) 其默认的幅值Ramp是从1降至0,不能使用,而要使用自定义的幅值曲线,使幅值在整个分析步中从0到1逐渐增大;
2) 正值是缝隙,负值是过盈量。
http://www.simwe.com/forum/viewthread.php?tid=691170
问:我所关注的只是螺钉拧入之活的应力分布,而对于拧入的过程则不关心。我想不放把问题说的更细一些。
把螺钉拧入骨头可分为如下几种情况:
A、自攻型,就是硬生生的把螺钉拧入骨头。
B、助攻型,先在骨头上打一个比螺钉的直径小一些的洞,之后再把螺钉拧进去。这样做的好处是螺钉不容易折断。打不同直径的洞,在螺钉拧入后的应力分布应该是不同的。
C、用一个同样的螺丝钉(不锈钢)在骨头中钻个洞,再把需要的螺钉(钛)拧进去。 D、或者有时用一个稍小一点的不锈钢螺钉钻洞,再拧入钛螺钉。这样可以比C更稳固。 E、我不小心把洞钻大了,这样螺钉就容易松动。当然这在临床中是不允许的,但也存在。
各位大侠能否就这几种情况分别讨论一下,我觉得是不是可以用过盈装配来模拟?但具体A、B和D有什么样的区别呢?想不明白。等待着各位大侠解惑~~~
答:这么说我就明白了,弄明白了就感觉很爽。
一定是兽医,把骨头这么搞来搞去的,让病人求生不得,求死不能,关云长来了也受不了啊。
A(怒射破门)就是过盈配合非常大。 B(后卫助攻)就是过盈配合小一些。
C(换上替补)就是几乎没有过盈配合。 D(影子前锋)就是过盈配合也比较小。 E(自摆乌龙)就是接触面之间存在缝隙。
具体过盈配合该多大,估计该在螺钉上装些传感器,得到试验结果,再调整你的过赢配合量,与实验相吻合。过盈配合的模拟方法见。
8.5 管土/桩土接触
http://www.simwe.com/forum/viewthread.php?tid=725275&page=1&extra=page=1#pid744132
问:隧道开挖与管土相互作用中的接触问题求教
本人做了隧道开挖中的管土相互作用的模型。有些管-土接触问题求教。
1. 模型说明:开挖隧道位于管道正下方,由于是研究接触算法,仅做三步开挖,且没有进行支护。土体使用实体单元c3d8r,管道采用三维壳单元s4.模型1管土接触采用法向
硬接触,切向罚函数(文件名exva-h1-fenbu)。模型2采用tie约束管道与土体节点。(文件名tie-ex-fenbu) 2. 计算后出现的问题 ( 1)开挖隧道纵向中间部位土体时,两个模型的管道位移基本一致,接触算法作用不大。
(2)隧道全部开挖的情况下:使用tie约束的模型2计算不收敛,通过查看塑性变形,发现和模型1的差不多,不知道为何不收敛。请问这种隧道开挖引起土体位移对管道的影响问题,设置接触的必要性以及模型2为何不收敛。
[ 本帖最后由 doudout 于 2006-5-12 09:30 编辑 ]
|Simwe.com|仿真设计有限元虚拟仪器h&o,B0}(r%m:NQ
www.simwe.comkP%z
~@WK
SimWe仿真论坛p2W&W8~www.simwe.comY:K6QQi8z,rAP
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMn(x:D图片附件: 模型1管道位移.jpg (2006-5-12 09:30, 68.48 K)
图片附件: 模型1竖向位移.jpg (2006-5-12 09:30, 98.26 K)
答:
或许是定义tie的时候应该加上 position tolerance=0.1, 以确保两个面tie在一起。
你的单元是c3d8, 精度很差,应该用c3d8r。
单位对吗?如果长度单位使用mm,那么质量的单位必须是“吨”,相应的密度单位就是吨/mm3,应力的单位是N/mm2(即MPa)。
见 http://www.simwe.com/forum/viewthread.php?tid=684273
http://www.simwe.com/forum/viewthread.php?tid=705685
问:计算中止,.msg文件里显示
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE 346 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 中谈了过约束的问题:
\文件中出现zero pivot说明ABAQUS无法自动解决过约束问题,例如在桩底部的最外一圈节点上即定义了tie,又定义了contact, 出现过约束。解决方法是在选择tie或contact的slave surface时,将类型设为node region, 然后选择区域时不要包含这一圈节点。\
你的桩底部只有一个单元,就没法用上述方法。要在底部细化网格才行。
不必整个模型都用那么细的网格,离接触面较远的地方可以粗一些。
http://www.simwe.com/forum/viewthread.php?tid=684042
问:
新手请教: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的初始时间增量步和允许的下限设小些,比如:
*Static
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)
http://www.simwe.com/forum/viewthread.php?tid=699614
问:大侠们看看这个初始应力为什么总是施加不成功了??
答:建议在step 1中只保留桩底部的接触,去掉其他接触。step 2中再加上所有接触。先做弹性材料,收敛了再做塑性。
http://www.simwe.com/forum/thread-499681-1-1.html
zero4: 用ABAQUS软件处理管土相互作用中的接触面问题
给有兴趣的朋友一个参考。 1/3
附件: 1.part1.rar (2005-6-10 13:36, 195.31 K)
http://www.simwe.com/forum/viewthread.php?tid=678974
问:做基桩的水平静载试验模拟,MSG文件提示迭代增量步太大,不能进行
取1/2原型,近场土C3D20R,远场CIN3D8,对称面上Y向对称约束。请赐教!
答:我认为模型中有以下问题:
1。塑性材料和接触面上都不能用C3D20R和C3D20单元,这可能是你收敛问题的主要原因。如果需要得到应力,可以使用C3D8I (在所关心的部位要让单元角度尽量接近90度),如果只关心应变和位移,可以使用C3D8R, 几何形状复杂时,可以使用C3D10M.
我尝试把你模型中的接触关系以及材料参数中的*Mohr Coulomb 和 *Mohr Coulomb Hardening都去掉,变为弹性材料, 仍使用C3D20R就收敛了。
2。接触对中的slave surface应该是材料较软,网格较细的面,你正好颠倒了。
3。 接触面之间有微小的距离,定义接触时要设定―Adjust=位置误差限度‖,此误差限度要大于接触面之间的距离,否则ABAQUS会认为两个面没有接触:
*Contact Pair, interaction=\
在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)
岩土问题我没搞过,回答不了你,论坛上好像有这方面的资料,你搜一搜。
单元扭曲的模型传上来看看。
不太了解你的单向多循环加载法。原则上,如果是静载荷,用standard的static分析步即可,如果加载时间很短(例如在地震中),或载荷本身的性质是动态的(例如来自旋转机械的载荷),在分析中不能忽略结构的惯性,就要采用动态分析,可以用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个单元,适当少一些可能也行。
*contact pair中的adjust=0.005,还是太小,在后处理时可以看到,接触面之间的距离大于0.005。把adjust设置大一些没关系,比如adjust=0.1。
在施加载荷前,应增加一个分析步,不施加载荷,只是让接触关系建立起来。原则是不要在一个分析步里让模型发生太多变化,那样就把ABAQUS搞晕了。
施加载荷时,如果是塑性材料,初始时间增量步应设小些,比如
*Static
0.1, 1., 1e-05, 1.
关于对称面上应力和位移,在后处理时用display group隐藏模型的一半即可。
http://www.simwe.com/forum/viewthread.php?tid=721435&pid=726515&page=1&extra=page=1&filter=type&typeid=2#pid726515 问:
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentOx S'yn7V仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM Y]'T)}整個被壓扁了 不如預期 照理 是要 基礎 沉入 土壤內 或是 對土壤 造成 沉陷的效果 ^^
a3K]5fS-q/M|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent9ZO\@
答: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来分割要细化的部位。
http://www.simwe.com/forum/viewthread.php?tid=696477
问:这个桩基础大致是这样的, 顶部是一个直径在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, ADJUST=?
来调整,使从面上单元节点调整到主面上.但可能会产生单元变形过大,但我觉得钢板厚度与网格尺寸相比较小的话,应该问题不大.
我没有做过土木,只是给你个参考吧,呵呵.
http://www.simwe.com/forum/viewthread.php?tid=725476
问:地基和土体的接触
我 作出以後
試著 加載 900N/m2 與9000N/m2 比較 但 跑出來的圖 不如預期 圖形如下
照理說 應該 是 基礎 沉入 土壤 如下圖
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3iT-仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'hH2ME `Simwe.com|仿真设计有限元虚拟仪器!{]/dN$z.~B仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMB;[5? k
答:你那张成功的图里,应力最大值是10e7, 而失败的图里,应力最大值只有10e5。如果材料相同,应力小,变形当然就小。
应力小的原因是载荷小,把载荷变为900000试试。
你那张成功的图里的deformation scale factor是1吗?如果是的话,这样的大变形要在step里设置几何非线性nlgeom 为 on.
你用了*Solution Technique, type=CONTACT ITERATIONS,不知是否可行。
一般都是先作地应力平衡,再加载荷,不知你这样一步完成对不对。
8.6 板料成形的接触问题
http://www.simwe.com/forum/viewthread.php?tid=680567
问:我在模拟一板料冲压成形过程中,采用sc8r单元,但模拟第一步的0.9999时就出现不收敛的情况,这个问题怎么解决,*.inp文件在附件里面
答:在msg文件中看到:
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE SHEET-1.49 D.O.F. 3 RATIO = 4.54101E+013.
这说明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.00010001, 对坯料产生微小的压力,以便产生摩擦力,消除刚体位移。
- 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。
- 在job 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 \中,对坯料施加一个向上的拉力,使其与punch贴在一起,否则无法收敛。下一个分析步中再去掉此拉力。
(经过以上改动,step \能收敛了)
- 网格在长度方向太粗,大变形时只有很少的节点与punch相接触,导致收敛问题。在长度方向将网格细化,宽度方向改稀疏一些。
- 去掉了*potential, 不过或许有它也可以。
可能还应该改的是,摩擦系数0.162是否太大了?有润滑的话应该是0.1或更小,这样收敛就更容易些。
目前看到的结果是punch下压了分析步的0.05,后面的还在算. 等不及后面的结果了,明天再说吧。
Ubend3.zip (106.84k)
http://www.simwe.com/forum/viewthread.php?tid=725202
问:在板料成形中出现的错误
答:在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 initial overclosures will result in initial accelerations on the order of 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。
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,选择第一个分析步\,点continue, 选择specify,点Time Incrementation 标签页,点第一个more,把Is由默认的12改为30,点OK。
相应的inp文件:
*Controls, parameters=time incrementation , , , , , , 30, , , ,
含义是:30次迭代不收敛时,才减小时间增量步。
这样改后就收敛了。另外,可以尝试把单元划得粗一些,就可以收敛得更快。
reduced integration单元要求在厚度上至少有4个单元。建议改为CPE4I单元。
http://www.simwe.com/forum/viewthread.php?tid=706374
问:这是个简单的拉伸成形模型,有几个问题没有解决:
(1)局部的单元变形太大,怎么控制让领近单元变形均匀点?
(2)我的回弹模型出了什么问题,一直没有成功,希望大侠能给指点一下,问题出在哪里!!
wy_anli答:在成形过程中,那个应变比较大的单元已经超过了你定义材料的最大应力,也就是说当应力达到你定义的材料最高应力时,就会把材料当成超塑性体。 (见图中的5.856E8,你定义的最大应力为4。7429E8) 建议减小压边力
1. 因为你采用的是显示算法,所有的加载应注意平滑加载,在计算完成后要看一下动能的历史输
出,要平缓,另外比较动能与内能,动能/内能一般要小于8%
另外那个BLANK_UP_FORCE可以减少甚至去掉,因为BLANK_UP_FORCE会增大成形难度。 2.回弹要在隐式算法中做
你的那个图是变形完成后的应力情况,在成形的开始的部分应力已经超过了极限值,你可以从开始看一下应力变化情况,也可以看我传的那个图
你先把成形的部分模拟完毕(修改压边力,使成形合理),再进行回弹模拟
http://www.simwe.com/forum/viewthread.php?tid=646079
问:这个模型为什么不收敛(超弹和弹塑性接触)?二维模型为拉弯三层板,因为是对称的,所以只取其中一半。上左为凸模,上右为压料板,下为凹模,中间为三层板。
三层板中间为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
http://www.simwe.com/forum/viewthread.php?tid=716122
问:(1)如果进行多道次的成形模拟,一般是怎么进行的?是不是:工序1(成形-回弹)-传递到-工序2(成形-回弹)-......这进行?
(2)第一道次的模具都可以在CAE里面完成设计等,哪后面工序用到的模具,有什么好办法添加进去:-直接在.INP文件里面修改?
先谢谢了!
答:可能可以在一个cae模型中把所有道次的模具都建在一起,然后在每个step中用*model change来杀死/激活不同的单元和接触。
8.7 凹坑成型的接触问题
http://www.simwe.com/forum/viewthread.php?tid=683533
问:大家好,我现在做的是一个有压管道在外加作用物下出现凹坑的准静态分析
当然这里面存在接触(外加作用物和管外壁之间),不要接触直接加力我可以做出来,但是加了接触之后就不行了,并出现以下
警告: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‖:这只是提示信息,一般不用理它。
rigid.cae (49.03k)
http://www.simwe.com/forum/viewthread.php?tid=730082
问:刚性体与弹塑性性的接触问题
我做的是个管道凹坑成型问题, 怎么算不了 这是我的cae文件
麻烦大家给我看看
M$u
仿真分析有限元模拟计算仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent8r8cmo
答:msg文件中有
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE ASSEMBLY.1 D.O.F. 3 RATIO = 7.3154E+013.
应该约束刚体参考点的U1, U3, UR1, UR2, UR3.
管道大变形处的网格太粗。slave surface的网格要比master surface细。
C3D8单元精度很差,应该用C3D8R.
http://www.simwe.com/forum/viewthread.php?tid=683101
问:我做的一个素混凝土梁,四点弯曲试验.用的是smeared cracking模型,各项参数详见附件文件.
前半段计算正常,可以看见上升趋势以及峰值.问题是过了之后出现典型的突降,计算开始不收敛,最后因为步长太小而中止.
改用更小的步长,划分更小的单元,都不奏效.我自己的分析是因为中部的单元被拉破坏,从而导致计算不收敛?可是从11方向的应力看不出破坏的迹象(否则应该11方向为0),梁底部依旧有拉应力,这是怎么回事?
Job-2.inp (179.09k)
答:请参考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,要先删掉所有静态分析步。
http://www.simwe.com/forum/viewthread.php?tid=726388
问:请问一下这样的问题能不能做:下图所示为一实体的中截面,打算这样计算,先让物体在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方向的移动,还可以向下移动一点点,以保证和下方实体相接触。下方实体不动。
8.8 刚体穿透
http://www.simwe.com/forum/viewthread.php?tid=722253&pid=726202&page=1&extra=page=1#pid726202
问:刚体穿透,aba大侠看过来 ~
模型很简单,就是一个2d的圆筒受压。采用shell单元。麻烦大家看看,多谢了
仿真分析
答:删除interaction property中的normal behavier, 就是默认的hard接触,不允许穿透。
后处理时比例设为1, 就看不到穿透。
impactboy 多谢,Aba,和前面一位大侠,问题解决了。都怪我没有注意deformation scale factor是否为1,误以为穿透了。分享一点
刚体穿透解决办法
1.边界条件和模型,材料是首先需要检查的地方,通产都是最容易发生问题的 2.尽量使用analytical rigid,而不是 discrete rigid 3.检查刚体和变形体的法向方向是否正确 4.网格划分是需要fine meshed,从面要比主面细 5.接触时主从面务必设置好,从面尽量使用bode-based element 而不是element-based 的surface
[ 本帖最后由 impactboy 于 2006-4-23 10:35 编辑 ]
www.simwe.com7V$z@$Sc8aSimWe仿真论坛:y6ZMwO Y.z|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|FluentD5W)y#q4W9I D;h.f%o$j0?6qQ仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMB A(t
8.9 接触力
http://www.simwe.com/forum/viewthread.php?tid=692527
问:做得是三维接触分析,主要想知道接触区接触力的具体数值,以便进行下一步单个零件的受力分析. 谢谢!!!!!
答:Step模块,设置History Output,把Domain设为Interaction,选择接触对,选中CFN。
用 *CONTACT PRINT 可以将接触力CFN输出至DAT文件:
*CONTACT PRINT, SLAVE=从面名称 CFN
http://www.simwe.com/forum/viewthread.php?tid=698091
问:我用二维模型,椭圆物压平面物,计算得到了几个接触部分的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是单位面积上的法向接触力(单位是 力/面积, 和接触力的数量级当然不一样)。
8.10 Explicit中的接触问题
http://www.simwe.com/forum/viewthread.php?tid=722822&page=1&extra=#pid727580
问:作切削时候遇见这种情况怎么办??? 刀具最后有点浸入工件里了 但是还在切削 不知道怎么回事
答:模型做得真不错。工件的节点没有进入刀具,只是由于单元变形太大,所以单元位于刀具里面,可以试试:
- 尝试其他的接触算法,比如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。
- 减小刀尖圆弧的半径。
http://www.simwe.com/forum/viewthread.php?tid=707106
问:关于切削加工仿真的,
|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent0@6d KN6f0R6r`.G U
答:你的模型里网格没划分,接触对Int-1也没定义,根本就运行不了。
接触面不能只是刀头的一个点吧,应该是刀的两条折线。explicit接触问题参见Getting Started with ABAQUS 12.7,12.8
想模拟材料破坏要定义失效判据,并让切下来的单元脱离模型,见ABAQUS Analysis User’s Manual 第11.6节Progressive damage and failure。切削问题很复杂,建议你先做一些简单的explicit接触问题。
http://www.simwe.com/forum/viewthread.php?tid=709083
问:我试过好几次,子弹就是不动,现传上模型CAE,请aba_aba帮我检查一下。非常感谢
strikemodel-6.4-1.rar (156.41k)
答:祝贺,你实现了传说中的隔山打牛。
别光写问我一个人,高手还是很多的,大家都发言论坛才有生命力,不能光靠一两个人。
察看子弹参考点的u,总是no value. 说明设置field out时没有输出子弹的计算结果,所以显示不出来。
改了以后,看到的是防弹衣把子弹挡了回来,真牛啊,刀枪不入。
http://www.simwe.com/forum/viewthread.php?tid=707952
问:模拟胀形过程,上下模具和压边都是解析刚体,工件为可变形体。 计算过程无提示错误,但是结果中工件和刚性体之间发生了穿透。 不知什么原因,附上inp文件。
答:刚性面必须是主面,刚性圆弧面的单元要足够细化,最好用analytical rigid.
ABAQUS Analysis User's Manual 21.4.1 Defining contact pairs in ABAQUS/Explicit
你用的
*Contact Pair, interaction=fric, mechanical constraint=PENALTY, cpset=Int-punch blank punch-1.Surf-punch, _PickedSet77_CNS_
我试了一下,把两个面换过来就不穿透了,但不知这样对不对。explicit我很少做。
Getting Started with ABAQUS 13.5 Example: forming a channel in ABAQUS/Explicit中用的是kinematic contact method,可以模仿一下。
ABAQUS Analysis User's Manual 21.4.4 Contact formulation for ABAQUS/Explicit contact pairs:
―As with the pure master-slave kinematic contact algorithm, there is no resistance to master surface nodes penetrating slave surface faces with the pure master-slave penalty contact algorithm. Using a sufficiently refined mesh on the slave surface will help correct this problem.―
http://www.simwe.com/forum/viewthread.php?tid=708886 的例子不错
8.11 改变接触状态
http://www.simwe.com/forum/viewthread.php?tid=674821
问:如何把两个不同位置的parts在它们经过运动后相遇时将相遇部分DOF锁在一起?
就是说,一开始两个parts在不相接触的位置, 然后一个part因受外力变型,使两个parts的某些elements相触. 这时我想把这些elements相关的地方(nodes or surfaces)的自由度锁在一起. 然后外力被撤去. 由于相关部分的锁定(tie, or bond, or what ever you wanna call it), 一个part企图恢复型变会使另一part也变形.
刚开始学abaqus几天,我能想到的一些方法如下,其中有几个已经试过不行了. 1)
试过将MPC放在非 1st step,
(n keyword *MPC, file \
It can be suboption for the following keyword: assembly, instance, part) 不行 2)
Abaqus Analysis User's manual 22.1.8 Breakable bonds 说:
*BOND can be defined only in the first step of a simulation 不行
2b ) Bond + remove/reactivating
Bond could be used only under ABAQUS/Explicit 3)
Abaqus keywords reference manual *fastener
Level: Part Part Instance Assembly 不行
4) tie的contact pair,先remove, 然后再Reactivating
Abaqus Analysis User's manual 21.2.4 Breakable bonds 说:
Removal and reactivation of contact pairs cannot be performed with ―tied‖ contact pairs
5) cohesive 不行(我试了,但试不出来)
6)connector element
答:try this:
*SURFACE BEHAVIOR, NO SEPARATION
Include this parameter to prevent any separation of the two surfaces once contact has been established.
注意要使用finite sliding,因为你的接触区域是变化的。如果使用small sliding, 开始不接触的区域在以后也不会再接触。
http://www.simwe.com/forum/viewthread.php?tid=684866
问:请问ABAQUS中能否将用*TIE 约束的接触面在STEP中解除邦定,如果能怎么实现这一功能,谢谢DX!!!
答:可能以下方法可以,但我没试过:
不用tie,用*CONTACT PAIR,TIED (在ABAQUS/CAE中无法定义,只能手写)
TIED : Include this parameter to indicate that the surfaces of this *CONTACT PAIR are to be “tied” together for the duration of the simulation. The ADJUST parameter is required when the TIED parameter is used.
然后在step中用 *model change, remove 来去掉接触。
问:Removal and reactivation of contact pairs in ABAQUS/Standard:
cannot be performed with ―tied‖ contact pairs
答:那就不用*CONTACT PAIR,TIED ,而是用普通的*CONTACT PAIR,但在CAE中定义接触属性的法向行为时,不选择Allow seperation after contact. 这样和*CONTACT PAIR,TIED的效果是一样的。以后再去掉接触即可。
还不行的话,再换一种 ―双面胶‖ 方法:在两个面之间定义一层很薄的单元,分别和两个面粘在一起(tie或者*CONTACT PAIR,TIED ),以后再把这层单元杀死,不知是否可行。
8.12 其他与接触有关的Error和Warning信息
http://www.simwe.com/forum/thread-724601-1-1.html
问:***NOTE: THE ROTATIONAL DEGREES OF FREEDOM OF THE NODES THAT CONSTITUTE THE
MASTER SURFACE WILL NOT BE CONSTRAINED BETWEEN THE SLAVE AND MASTER
PAIR (QWTSURF, QWQSURF).
***WARNING: NODE 2 ON SURFACE QLTSURF HAS FACETS WITH NORMAL VECTORS DIFFERING BY MORE THAN 30 DEGREES. CONVERGENCE DIFFICULTIES MAY OCCUR AT THIS NODE WITH FINITE-SLIDING CONTACT. THE NORMAL CONTACT DIRECTION AT THIS NODE WILL BE (-0.70711,0.70711,0.0000). 是什么意思,该如何改正?谢谢!
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMl)e g;c?8g\仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,a6}6yS
|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|FluentQ6U/x*UQbO6C
答:这个提示信息可能不用管它,只是说在接触中不考虑壳单元的旋转自由度。
http://www.simwe.com/forum/viewthread.php?tid=722397&pid=729331&page=1&extra=#pid729331
问:小弟使用如aba_aba前輩之教學的方式F2D2來用...
不過在Part2 的接觸面上卻出現錯誤...
\
ELEMENT TO DEFINE A SURFACE IN THIS CURRENT VERSION. THIS ELEMENT WILL NOT BE CONSIDERED AS PART OF THIS SURFACE.\
搞了幾天都不知要改哪... 可否請大大幫我看一下... 感謝~
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentw-}K!Rm6F仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM5}O1s0N-w5?仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMmm zM#R5[B$xwww.simwe.com/H lv2?/D6N-k(o-}
答:可能是这种单元不能定义基于单元的surface, 可以试试使用基于节点的surface,方法是: 不要直接定义surface, 而是定义set, 然后在定义contact时,将slave surface的类型设为node region, 选中此set.
9. 各个量的单位
9.1 各个量的单位要相互一致
http://www.simwe.com/forum/viewthread.php?tid=684273
问:请问,ABAQUS中有专门设置单位的地方吗? 我发现有两种单位组合,但如何确定种类呢?
答:
ABAQUS中的量都是没有单位的,但用户要自己保证量纲的统一,例如可以都使用国际单位,即m、kg、N、s等,相应的密度单位是kg/m3,应力的单位是N/m2(即Pa)。
需要注意的是,如果长度单位使用mm,那么质量的单位必须是―吨‖,相应的密度单位就是吨/mm3,应力的单位是N/mm2(即MPa)。
必须同时使用mm和―吨‖的原因是:根据―力 = 质量 加速度‖,各个量之间的换算关系为:N = kg ? m / s2 = 0.001吨 ? 1000 mm / s2 = 吨 ? mm / s2。
-----------------------------------------------------------------------------
焦 = N * m, 所以如果用N, mm,能量就要用mJ (0.001焦)。
可用的单位是 长度 m 力 N 质量 Kg 时间 s 应力 Pa 能量 J
密度 Kg/m3
例如:
*Density 7800, *Elastic 2.1e11., 0.3
*Plastic 4.2e8, 0
或者
长度 mm 力 N
质量 Tonne(1000kg) 时间 s 应力 MPa 能量 mJ
密度 Tonne/mm3
例如:
*Density
7.8e-09, *Elastic
210000., 0.3
问:我认为你说的不对
你使用长度单位为mm,质量单位可以随便定就用kg就可以,时间单位是另一个 基本单位,其他各个量都可以根据这三个来换算得到,当然热的不包括,而不是 你说的mm必须对吨,heihei
补充一下,因为你不必一定让力的单位是牛顿,
可以是达因,或者KN,MN等等
答:
你说的在理论上是可以的,但如果长度单位为mm, 质量单位用kg,根据 N = 吨 ? mm / s2,就有:
0.001N = kg ? mm / s2
也就是说力的单位是0.001N, 即 0.1 达因, 或者0.000001 KN, 0.000000001 MN, 这些用起来都很不方便。所以还是该用mm, 吨,N, 或者国际单位。
http://www.simwe.com/forum/thread-726603-1-1.html
问:求救,這樣的施力為何會出現錯誤...
答:单位不统一。你的力是1471N,材料数据是
*Uniaxial Test Data 20.1625, 5.473 ......
意思是名义应力是20.1625 N/mm2 还是 N/m2 ? 你单元的长度约为 0.058. 总之这样的材料属性和模型尺寸,会被1471N压得稀巴烂,所以不收敛。见
http://www.simwe.com/forum/viewthread.php?tid=684273
我把力改为-0.01471, 就收敛了。另外,施力点要约束U1和UR3 (基于你的局部坐标系)。见附件。
材料参数的数量级对吗?Getting Started with ABAQUS 10.7.2 中的Uniaxial test data 是
Stress (Pa) Strain 0.054E6 0.0380 0.152E6 0.1338 ......
它的单位是Pa, 你用MPa。
用MPa就意味着你的长度单位是mm。你模型的宽度是64mm,这样小的一块薄橡胶受到150kg的压力,肯定会被压成扁扁的一小团,按照你现在的边界条件,单元畸变会非常严重,无论怎么划分网格也很难收敛,除非是用ABAQUS 6.6的remesh功能。
你试验中的橡胶厚度是多少?模型中plane stress的单元厚度是1,承受的压力就应该不是150kg了吧。plane stress问题中的pressure含义见
http://www.simwe.com/forum/viewthread.php?tid=676775
另外重要的一点,你上下都用tie, 橡胶两侧在x方向无法自由移动(如同被放进了一个罐头盒里),所以压缩到一定程度就压不下去了,看到的现象就是不收敛。建议去掉下面的tie和两侧的边界条件(如果符合实际情况的话),改为橡胶底部边界条件U2=0,对称模型取一半,对称面上U1=0. 力载荷可以改为位移载荷(顶部刚体向下移动一段距离),就更容易收敛。
超弹性不知能不能用二阶单元,用一阶单元可以大大减少计算时间。超弹性是否该使用杂交单元(hybrid formulation),
但plane stress不能用hybrid formulation,plane strain才可以。
角部网格畸变也是不收敛的重要原因,可以通过partition来改善,见Getting Started with ABAQUS 10.7.2 图10-63,图10-64.附件是我做的修改。
http://www.simwe.com/forum/viewthread.php?tid=683628
问:我练习那个用户指南里的那个斜板的例题,线性的算出来了,可是非线性的一步都没有算就出错了, 我是按照书上做的,错误原因有1次是说那个attempt 太多了,还有几次就没有说明.请大虾们指点迷津.另外,错误信息到哪里去看啊,才学几天,问题很弱 ,望高手给予帮助。
答:
把标准答案中的CAE模型和INP文件与你的一点点比较,终于发现了几点区别:
最重要的是,你模型的几何尺寸是标准答案中的100倍,厚度却是一样的,也就是你的单元太薄,把shell section的厚度改为0.8就收敛了。或者把单元加密100倍可能也行。
另外,你的弹性模量是3E+9, 应该是3E+10。
你也没有为材料定义方向,但这可能问题不大,只是得到的结果会和标准答案不同。
所以,只要你有耐心,把标准答案和你模型的区别一点点找出来,总是能找出问题所在的。
9.2 时间/增量步的含义
http://www.simwe.com/forum/viewthread.php?tid=680199
问:Abaqus的History Output中的时间是真实的时间吗?还是时间增量步的概念?
答:
如果是和时间无关的静态分析就不是真实事件,如果是和时间相关的蠕变分析、动力学分析等就是真实时间
http://www.simwe.com/forum/viewthread.php?tid=725493&pid=743800&page=1&extra=page=1#pid743800 问:在稳态分析中,和时间无关,为何也要定义初始时间步,总时间?
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluento4GF z;op2[6P-k*SOILS
wf*W^
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentR-\\!g c_ j
Initial time increment,Time period
是不是不用定义Initial time increment,Time period 也可以?
SimWe仿真论坛u.s+K*H's#v5}7A&m
答:*SOILS我没用过,但估计和*static是类似的,Initial time increment和Time period 是相对的时间概念,increment小,则load之类的量就小,就容易收敛。
10.导入 / 修改 / 求解inp文件
10.1 把inp文件导入ABAQUS/CAE
http://www.simwe.com/forum/viewthread.php?tid=683623
问:用户手册里的例子里的各种.inp文件是否可以变成ABAQUS CAE 中的模型建立过程,这样我就能看懂了,很不好意思问这种幼稚的问题,请大虾们指点,谢谢.
答:
菜单File / Import / Model. 但要注意下面贴子中提到的问题:
http://www.simwe.com/forum/viewthread.php?tid=679644
问:
为什么我导入的help里的INP文件都不能直接进行求解呢?
如:abaqus example : 5.1.1 Thermal-stress analysis of a disc brake。里面的多个INP文件,我导入了―discbrake_3d.inp‖,进行求解,为什么出现错误?
哪位朋友告诉我一下,学习help时,help给出的INP文件,如何处理? 在此表示感谢。
答:
INP文件中的有些关键词是ABAQUS/CAE所不支持的,导入INP文件时会在窗口底部的信息区中看到警告信息:
WARNING: The following keywords/parameters are not yet supported by the input file reader.
这时导入ABAQUS/CAE的模型是不完整的,所以分析时出错。
如果在ABAQUS Command窗口中输入以下命令来提交分析,就没有这种问题:
abaqus job=INP文件的名称
问:我输入:abaqus job=inp2d.inp 按回车后出现如下: SyntaxError: invalid syntax
答:不是在ABAQUS/CAE底部的>>> 后面输入,而是:
在WINDOWS中点击 [开始] → [程序] → [ABAQUS 6.5-1] → [ABAQUS Command],然后在DOS窗口中输入:
abaqus job=2D
http://www.simwe.com/forum/viewthread.php?tid=725948
问:我在看有些帮助例子时,为什么有些inp文件输入以后什么也没有
是没有part吗,但是打开inp文件,节点单元都定义有阿!
abaqus 中的example 下面的Submodeling of a stacked sheet metal assembl 我将stackedassembly_s4r_global.inp打开,没有错误提示,就是什么也没有! 帮忙看看,谢谢!
答:CAE窗口底部有提示信息: AbaqusException: Include file \
需要把这个例子中的stackedassembly_s4r_global_mesh.inp也存在同一个目录下。
10.2 运行inp文件
http://www.simwe.com/forum/viewthread.php?tid=702294
问:如何在command中运行.inp文件? 我运行之后出现
ABAQUS warning:The .inp extension has been removed from the jod identifier Input file :
这种问题应该如何解决?
答:inp文件要放在command窗口所显示的目录下,“job=” 后面的文件名要和inp文件名完全相同。
http://www.simwe.com/forum/viewthread.php?tid=698349
问:只知道在command 中提交分析,可是计算途中想修改模型,可是不知道该怎么kill job ?
shock111答:abaqus suspend j=xxxx 这个任务是可以恢复的
abaqus terminate j=xxxx 这个就直接杀了
handful答:abaqus suspend j=xxxx
abaqus resume j=xxxx
从上次分析结束的地方重新开始分析
答:或者也可以在windows任务管理器里直接中止进程standard.exe或explicit.exe
http://www.simwe.com/forum/viewthread.php?tid=694281
问:我想在matlab中反复调用abaqus,
1)如何判断abaqus是否完成运算(因为要等它算完后再调用,要调用很多次>2k,同时执行好象也不可能执行这么多任务的)
2)如何从.dat的文件中读数据?
我本来在建模时都是手动打开.dat的文件,找相关数据的(我在inp中加入*node print).
但现在要调用abaqus数千次(因为要用mont carlo法),每一次都要读1-2个数据,并把数据赋给一数组,然后用matlab进行相关运算得出结论。
我怎样才能让matlab自动去读abaqus运算完后我所需要的特定的几个点的位移数据?
答:C 和 FORTRAN 都有函数来读写文件,MATLAB估计也有吧。
ABAQUS分析顺利结束后,会在msg和dat文件结尾显示\MATLAB读到这个,就知道ABAQUS分析结束了。
一些可供借鉴的:如果想在dos窗口中依次运行两个job, 可以自己写一个bat批处理文件,里面包括以下语句:
call abaqus job=*** interactive call abaqus job=*** interactive
其中的interactive可以保证分析完一个job, 再分析另一个。
10.3 修改inp文件
http://www.simwe.com/forum/viewthread.php?tid=689187
问:我在inp文件定义了三个节点集:bottop,toptop,toptop, 然后在最后一个step最后设置了所定义节点的坐标输出,
结果inp文件前处理错误。请对此了解的同仁帮我看一看,谢谢! 附上inp文件
答:对于在part中定义的集合,在后面引用时要加上instance名称:
*Node print, nset=blank-bot-1.bottop coord
http://www.simwe.com/forum/viewthread.php?tid=697544
问:CAE告诉我说,An initial condition has been specified on element 0 but this element has not been defined……
我tmd的哪来的一个element 0 ~~~~!!
答:inp文件中initial conditions后面的数据有错误,比如有空行,ABAQUS就把那个错误的数据理解为了Element 0。
不要急躁,一般都是你自己的错误,极少数情况才是ABAQUS的bug.
同一行数据之间的空格不要紧. inp文件中的常见错误有:
- 出现空行
- 数据之间没有用逗号分开,只用了空格 - 用了中文的逗号
- keyword拼写错误或位置错误
http://www.simwe.com/forum/viewthread.php?tid=723358&pid=730354&page=1&extra=#pid730354
问:我尝试使用一种在cae下有的单元类型,然后再修改inp文件,改成无限元单元,用来定义边界;
出现的问题:
1。在cae中定义的单元类型较少,假如使用与实体其它部分(partition过)相同的单元类型,在生成的inp文件中,将该类型相同的单元全部放在一起:
2。好像在cae下直接修改的关键字命令,似乎不太好用,而且你在cae下的修改,不能即时地在cae下的菜单中修改keywords,好像没有刷新。
自己觉得可能的解决方法:
1。定义单元集,然后,再在inp文件中,增加单元类型定义,重新定义成无限元,这样不知道对单元的二次定义,会影响吗?? 2。inp文件的导出似乎不太明确,自己无法直接导出,只能导入,只有在提交job的时候出现inp文件,不知道理解得对不对!!
希望哪位达人能指导一下!!
www.simwe.com-La;J7}u\q#Cwww.simwe.comN/l1w2^Q7lH7X
www.simwe.com2h
答:在inp文件中改了之后,直接在dos窗口中提交分析,不要再import到CAE里面了。
CAE中选择的单元形状、节点数、坐标要和无限元完全一样,在inp中只是把*ELEMENT,TYPE=*** 改为*ELEMENT,TYPE= $$$.
提交job的时候, 不要点submit, 而是点write input, 就输出了inp文件。
http://www.simwe.com/forum/viewthread.php?tid=724350&pid=737435&page=1&extra=#pid737435
问:在计算接触时有这个错误:
***ERROR: in keyword *CONFLICTS, file \ *Conflicts is generated by ABAQUS/CAE to highlight changes made via the Keywords Editor that conflict with changes made subsequently via the native GUI.
请高手指教一下这是哪里出错了
www.simwe.com;lD2h-y'Q:N5P
答:在edit keywords之后修改了模型,导致相互冲突的keywords。在edit keywords中可以找到*CONFLICTS, 修改正确即可。
你是否用过菜单Model / Edit Keyword,添加了什么语句?如果是,再次进入Edit Keyword,找到*CONFLICTS, 删掉它,把你添加的语句移至正确的位置。
http://www.simwe.com/forum/viewthread.php?tid=726921
问:经过十几天的努力,终于搞定了这个INPUT文件,可惜总是运行不了,检查了很多遍也不知道错误在哪里? 只好来请教高
手。希望哪位高手赐教!
ABA_ABA 请你帮忙看看,多谢了!
|Simwe.com|仿真设计有限元虚拟仪器d2`%K
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentJ$JE4g {3?
答:错误信息: ABAQUS Error: The following file(s) could not be located: concrete.inc
因为你的inp文件里有*INCLUDE,INPUT=concrete.inc。我暂时去掉include这一行,运行后在dat文件里看到大量错误信息,其中大多是关于自动生成的pes文件中的错误。
建议你先不要搞这些复杂的参数化输入,从简单的inp文件开始,把接触问题做出来就不容易了。
运行后在dat文件里看到错误信息,例如 ***ERROR: in keyword *ELEMENT, file \Non integer in element connectivity for element 51。
根据这个提示,看生成的文件 Bikont.pes 的第51行,下面出现很多问号,这里本应是节点编号的。Bikont.pes 里所有出现问号的地方恐怕都是不正常的。
解决方法是把你的inp文件简化,不要一下子就使用这么多参数化输入,先只做一两个单元,一两个参数,不要用接触,成功了再作复杂的。
11.用户子程序和FORTRAN
http://www.simwe.com/forum/viewthread.php?tid=683381
问:我是第一次调用subroutine.调试的*.inp和*.f也是从abaqus中的help找来的. 在调试过程中,已通过了fortran的编译,生成了*.obj文件,我已将以下的文件
dp.exe,aba_para_sp,aba_para_dp,并将aba_para_sp改为aba_para,但调试仍出现以下问题(如图),请大家看看是什么问题?谢谢了
答:
在本版搜索―valueerror‖,可以找到很多。
我运行没问题。ABAQUS要求的FORTRAN版本是Compaq Visual Fortran 6.0。
不用自己在Fortran里编译*.f。直接把帮助文件里的 uhetvalx.inp 和 uhetvalx.f 保存在默认的工作目录下,把 uhetvalx.f 改名为uhetvalx.for. 然后输入:
abaqus job=uhetvalx user=uhetvalx interactive
答:
我总算把这个问题解决了,是fortran软件忍得祸。因为我在另一个机器也装了abaqus,和compaq visual fortan 6.5,还是出现同样的问题。我就开始想,是不是软件的问题,找了另一个形式的(里的内容一样,只是serial不同),安装完后,verify 成功,真高兴呀!
http://www.simwe.com/forum/viewthread.php?tid=729244
问:INP文件和.for文件一起执行加入umat...
偶用document里的文件来运行:abaqus job=rotcyl2d user=rotcyl2dvumat interactive结果报错
roblem during compilation-df.exe not found in PATH.
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM6iOwx6Oj2k p
答:
首先要保证安装了正确的fortran版本,见ABAQUS Installation and Licensing Guide A.1 Software。 Windows下的ABAQUS 6.5 要求是Compaq Visual Fortran 6.0。
错误信息Problem during compilation- df.exe not found in PATH, 意思是找不到fortran的df.exe , 可以设定一下windows的path.
http://www.simwe.com/forum/viewthread.php?tid=723356
问:小弟新学ABAQUS(版本6.4),课题方向是结构抗火,现在有问题请教如下:
我想分析帮助文件里面的例子ABAQUS Example Problems Manual(5.1.4 Radiation analysis of a plane finned surface),里面给出了两个文件:Input file
radiationfinnedsurf.inp Fire transient problem.
radiationfinnedsurf.f User subroutine FILM used in radiationfinnedsurf.inp. 其中第二个是inp文件所调用的子程序subroutine film,用于定义与温度相关的热传导系数。 我是这样做的:
1.按照有人的说法将子程序的后缀改为.for(不改的也试过)
2.将这两个文件拷贝到我abaqus分析目录D:>Temp下,然后用命令运行: D:\\ Temp>abaqus job=radiationfinnedsurf user=radiationfinnedsurf
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluenta:Pe}
仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent1`%W#R8i%O Nw^
SimWe仿真论坛e)v0D-@\l8~5tMU5SlOJ-X~ Q|www.Simwe.com|有限元模拟计算力学航空航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|FluentV ma,~z,w
问题是: 电脑告诉我找不到子程序。
答:
1.将这两个文件拷贝到分析目录D:>Temp下 2.将子程序的后缀改为.for
3。 然后在D:\\Temp下面用命令运行
abaqus job=radiationfinnedsurf user=radiationfinnedsurf
这样就对了。注意如果用记事本之类的软件保存了子程序文件,它会自动在扩展名.for后面再加上.txt,可以在dos下用dir命令看看扩展名是否真的是.for
http://www.simwe.com/forum/viewthread.php?tid=700587
问:我装了ABAQUS6.5然后又装了Vc++.net和Intel fortran8.0怎么ABAQUS不能调用,请问它们的正确安装顺序和该怎样设定环境变量?
答:windows下,ABAQUS6.5要求的是Compaq Visual Fortran 6.0 和 Microsoft Visual C++ Version 6.0。
Linux/x86-32下,ABAQUS6.5要求的是Intel Visual Fortran 8.0 和 Intel Visual C++ 8.0。
我是先装fortran, 再装ABAQUS。你是ABAQUS不能用,还是用户子程序不能运行?
http://www.simwe.com/forum/viewthread.php?tid=705429
问:我在用UMAT计算过程中,算到中途就碰到这样的问题:在MSG、DAT和STA文件中都没有看到什么错误,但是计算却停止了,只在LOG文件中有这样一段话,请大侠们帮我分析一下。 ABAQUS Error: The executable D:\\abaqus\\6.5-1\\exec\\standard.exe aborted with system error \拒绝访问。 \
答:在帮助中找个UMAT的例子试一下,如果一切正常,就不是用户子程序安装的问题,而是在运行你自己的UMAT过程中出现了错误,原因可能是UMAT程序有问题或调用UMAT时所传递的参数错误。
12.创建/修改part