找回密码
 注册
Simdroid-非首页
楼主: blueshell

[INP文件] *Model change 的应用技巧初步

[复制链接]
发表于 2008-9-11 18:05:52 | 显示全部楼层 来自 四川成都

请问blueshell大哥

blueshell大哥好
    老师给了一片文章让再现一下,可做起来遇到很多问题,还请blueshell大哥帮忙看一下:
    文章主要是为了模拟摩擦磨损过程中表面有微凸体作用时的力学行为,所以自己考虑着采用把一些表面单元给删除掉来实现。
    主要分析步骤:
    step1:将一个半球向下压向一平板,位移为0.5mm,
    step2:将半球上边界固定,给平板一固定频率和唯一幅值的三角位移
    step3:删除接触表面的一些单元,并让平板继续滑动,来模拟微凸体的影响
    遇到问题;
    1:通过在step3中添加*MODEL CHANGE, TYPE=ELEMENT, REMOVE来实现时,step3根本没有运算,并有
     ***WARNING: THE STRAIN INCREMENT HAS EXCEEDED FIFTY TIMES THE STRAIN TO CAUSE FIRST YIELD AT
     1POINTS出错信息,请blueshell大哥帮忙用钝化单元帮忙修改一下。
    2:不知表面删除一些单元后,刚开始采用的面面接触方式是否还起作用,删除单元后,接触面间的接触关系是不是因删除了 一 些单元需重新建立。
    3:删除单元后,微凸体之间间隔交错出现的关系实现了没有。
   麻烦blueshell大哥了,这里小弟谢谢先

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

 楼主| 发表于 2008-9-11 19:08:39 | 显示全部楼层 来自 韩国
Simdroid开发平台
首先声明,涉及接触的问题一般都是很复杂的,我对这方面还不了解。从运行后提供的信息来看,你主要的错误来源,不是你上边列出的那个warning,而是下边的这些warning,也给出了一些提示。一般情况下,*Model change, remove 掉单元后,与单元相关的性质都钝化了。也许问题正是出在这里,也就是说,钝化掉相应单元后,你原先有的接触就跟着钝化单元消失了,此时的计算就会出现多次数值迭代不成功的情况。你看看信息文件,也许对你有帮助,并相信你能找到解决办法。
***WARNING: THE SYSTEM MATRIX HAS 4 NEGATIVE EIGENVALUES.
             EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: SOLVER PROBLEM.  ZERO PIVOT WHEN PROCESSING NODE 92 INSTANCE
             PART-1-1 D.O.F. 1
***NOTE: A ZERO PIVOT WARNING (IF THE *MODEL  CHANGE OPTION WAS USED) MIGHT
          INDICATE EITHER AN OVERCONSTRAINT OR THAT AN ELEMENT REMOVAL HAS
          REMOVED ALL OF THE STIFFNESS ASSOCIATED WITH ONE OR MORE
          DEGREES-OF-FREEDOM FROM AN ACTIVE NODE.
***NOTE: FOR EXAMPLE, IF A SHELL ELEMENT AND A BRICK ELEMENT SHARE A COMMON
          NODE, AND THE SHELL ELEMENT IS REMOVED, THE BENDING STIFFNESS AT THE
          NODE WOULD BE SET TO ZERO BUT THE NODE WOULD STILL BE ACTIVE SINCE
          THE BRICK ELEMENT IS STILL PART OF THE MODEL.  ONE POSSIBLE FIX IN
          THIS CASE WOULD BE TO CONSTRAIN THE INACTIVE ROTATION
          DEGREES-OF-FREEDOM USING A BOUNDARY CONDITION.
***NOTE: IF AN OVERCONSTRAINT OCCURS, THE CHAINS OF CONSTRAINTS CAUSING THE
          OVERCONSTRAINT CAN BE PRINTED OUT USING *CONSTRAINT
          CONTROLS,PRINT=YES.
OVERCONSTRAINT CHECKS:
An overconstraint was detected at node 92 INSTANCE PART-1-1. There are multiple constraints
applied directly or chained constraints that are applied indirectly at this node.
A list of nodes and chained constraints between these nodes that most likely lead
to the detected overconstraint can be printed out to the message file by including
*constraint controls, print=yes in the input file.
回复 不支持

使用道具 举报

发表于 2008-9-13 09:06:34 | 显示全部楼层 来自 四川成都

谢谢blueshell大哥

blueshell大哥,恩啊。我再把接触改动后试一下,还有blueshell大哥知道:
如何把一个面上的节点应力和位移信息输出到一个文件里,而在修改完的模型中作为初始边界条件加进去吗?感觉和子模型不一样,但又不知道怎么弄,自己刚接触这方面的东西还有很多不懂的地方,还望大哥指点啊!
回复 不支持

使用道具 举报

 楼主| 发表于 2008-9-13 15:28:40 | 显示全部楼层 来自 韩国
为实现这个功能,可能要用到子程序USDFLD。由GETVRM函数获得你关心的节点应力和位移,添加语句,可将你关心的内容输出到文件或存到状态变量(State variables)中,存在状态变量中的数据可以被修改,也可以被其他子程序调用,也许还能在后续的计算中调用该数据(这点还没去了解)。由于涉及到的问题比较多,建议你看看Abaqus Keywords Reference Manual   中的 *USER DEFINED FIELD及相关的知识点。

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2008-9-13 16:10:05 | 显示全部楼层 来自 四川成都
恩啊,谢谢blueshell大哥。我好好看一下,不知道大哥有没有这方面的例子什么的,有的话感觉着对比着效果会好一点。再次谢谢blueshell大哥,好人啊,呵呵
回复 不支持

使用道具 举报

发表于 2008-9-16 15:26:39 | 显示全部楼层 来自 河北秦皇岛

blueshell

大哥,
你那天给我那个我算出来了,现在做了一个大件,为什么算不下去,出现下面的错误,我就是按以前那样修改的INP文件。还有就是那个复制单元集在CAE里怎么设置啊。

***ERROR: in keyword *MODELCHANGE, file "1.inp", line 21060: Unknown assembly
           set SET-1,
***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO
          INTERPRET SOME DATA.  SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION

[ 本帖最后由 张晓敏 于 2008-9-16 15:29 编辑 ]
回复 不支持

使用道具 举报

 楼主| 发表于 2008-9-16 20:22:05 | 显示全部楼层 来自 韩国
复制单元因为与钝化单元完全重合,所以,最好在INP文件中进行处理。所以,一开始就得对要钝化和非钝化的单元进行区分,最好在assembly时设置不同的单元集,以便在打开inp文件的时候,可以很方便的进行处理。CAE中不支持Model change,只能在inp文件中添加。提示的错误信息是set-1出错,可以进行针对性查找。

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2009-1-22 16:02:34 | 显示全部楼层 来自 内蒙古兴安盟
谢谢楼主!!!
回复 不支持

使用道具 举报

发表于 2009-1-23 20:46:08 | 显示全部楼层 来自 广东梅州
如果钝化和非钝化单元上还作用了相同的表面均布荷载时,model change时,怎么加呢?
回复 不支持

使用道具 举报

 楼主| 发表于 2009-1-24 09:57:52 | 显示全部楼层 来自 韩国
Any distributed loads, fluxes, flows, and foundations specified for inactive elements are also inactive. However, unless you explicitly remove them, records of these loads are still kept and are listed in the data (.dat) file as though the elements were still present. Continuation of loads across steps is not affected by removal; on element reactivation unremoved distributed loads are also reactivated.

By default, if a distributed load is applied to an element that is being reactivated in a step, the distributed load magnitude is scaled up linearly from zero to its end-of-step value during the step. If such a load is applied with an amplitude reference, the magnitude value given by the amplitude reference is scaled again by a value that ramps from zero to one throughout the step. This scheme ensures that reactivation has a smooth effect on the solution, even in cases where a distributed load with an amplitude reference on a reactivated element is carried over from a previous step.
为了避免出错,可以考虑将remove和reactive分别建一个分析步,其上的荷载在随后的分析步中加上。
回复 不支持

使用道具 举报

发表于 2009-2-12 11:20:11 | 显示全部楼层 来自 广东广州
如果不采用备份钝化单元的方法,直接杀生,激活单元。连接部分的单元会有什么误差呢?
因为当部件很多时,组装也较复杂时,采用你说的这种方法会很麻烦。
回复 不支持

使用道具 举报

 楼主| 发表于 2009-2-12 14:16:42 | 显示全部楼层 来自 韩国
你说的是事实,在结构很复杂的时候,这种方法在建模以及后处理的时候都很复杂。至于是否存在其他方法,我不是很清楚。但根据个人经验,如果没有共享节点的“完全弹性单元”,则钝化单元会在其原始位置被激活,以至于在后处理时会出现钝化单元与原始结构(非钝化部分)之间出现变形差异。要获得全结构,尤其是钝化部分的真实变形,同样需要对相应的结果进行进一步的分析处理。
回复 不支持

使用道具 举报

发表于 2009-2-16 20:30:39 | 显示全部楼层 来自 福建福州
好东西
回复 不支持

使用道具 举报

发表于 2009-2-17 00:59:50 | 显示全部楼层 来自 北京
多谢lz分享
回复 不支持

使用道具 举报

发表于 2009-2-17 01:26:02 | 显示全部楼层 来自 重庆北碚区
多谢lz分享
回复 不支持

使用道具 举报

发表于 2009-2-17 01:47:15 | 显示全部楼层 来自 重庆北碚区
好像曲线有问题?
回复 不支持

使用道具 举报

发表于 2009-2-17 04:33:44 | 显示全部楼层 来自 美国
很有启发,谢谢楼主
回复 不支持

使用道具 举报

 楼主| 发表于 2009-2-17 09:04:03 | 显示全部楼层 来自 韩国
难得大家如此捧场,特奉献一个小工具。由于采用了“完全弹性”的共享节点的重复单元,因此,在钝化单元对应的节点上的应力、变形或其它结果都会由两种单元(钝化单元本身和附加的重复单元)换算过来的。而重复单元主要是用来跟踪定位的,其结果一般都很小,并不是我们需要的,因此需要将这部分结果过滤掉。特用excel的宏命令来实现这个功能(附件excel文件),只需将需要整理的数据拷贝到第一第二列,点一下计算按钮就可以得到需要的结果。enjoy

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

发表于 2009-2-23 14:39:32 | 显示全部楼层 来自 广东广州
不知道楼主有没有分析过 with strain 和strain free的区别? 用with strain时,对解决被杀死单元与激活单元共同变形是否会好一点?
回复 不支持

使用道具 举报

 楼主| 发表于 2009-2-23 14:59:40 | 显示全部楼层 来自 韩国
我没有做过 WITH STRAIN 方面的分析。而帮助文件中
ADD=WITH STRAIN to specify that stress/displacement elements are reactivated with strain.
此处的with strain 我也没搞明白是代表什么意思,究竟是额外指定的strain,还是直接就是跟未钝化单元共同变形?你倒可以一试。
我前面提供的这种方法本身就能很好的解决共同变形的问题,只是在选取所谓的“完全弹性单元”时,要注意选取其材料参数,应避免两点:一尽量减少对未钝化的单元的影响;二尽量使其自身的应力足够小。以达到跟踪未钝化单元节点变形后坐标的目的。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 02:58 , Processed in 0.044407 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表