找回密码
 注册
Simdroid-非首页
查看: 5285|回复: 152

[热分析] [原创]:焊接模拟之圆管对接

[复制链接]
发表于 2011-12-20 21:59:07 | 显示全部楼层 |阅读模式 来自 北京
经过版内几位朋友的指导,做了一个圆管对接焊的模拟.效果差了点,主要是参数设置的比较随便,另外电脑硬件限制,模型不能分的很细,但是思想应该是可以说明白.


附件包括焊接的子程序和模型,建模的方法和平板的类似,不过子程序中需要对平板焊接的子程序进行修改(主要是坐标变换),具体可以看看附件.




欢迎讨论,谢绝灌水.


本帖子中包含更多资源

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

×

评分

2

查看全部评分

发表于 2015-11-24 21:33:56 | 显示全部楼层 来自 湖南长沙
Simdroid开发平台
版主,在子程序中如何控制焊接方向呢?
到底是顺时针还是逆时针?
回复 1 不支持 0

使用道具 举报

 楼主| 发表于 2011-12-21 08:38:40 | 显示全部楼层 来自 上海
smallmj 发表于 2011-12-21 08:04
这个很不错,问个问题,子程序能不能用python来写?

子程序只能fortran了,python是用来做前后处理的。
回复 1 不支持 0

使用道具 举报

发表于 2011-12-21 08:04:44 | 显示全部楼层 来自 广东深圳
这个很不错,问个问题,子程序能不能用python来写?
回复 1 不支持 0

使用道具 举报

发表于 2011-12-21 09:57:00 | 显示全部楼层 来自 重庆沙坪坝区
这个必须顶啊,谢谢
回复 不支持

使用道具 举报

发表于 2011-12-26 09:18:37 | 显示全部楼层 来自 安徽合肥
谢谢楼主的分享            
回复 不支持

使用道具 举报

发表于 2011-12-26 09:24:33 | 显示全部楼层 来自 安徽合肥
beta=0.0
      aacos=a/cos(beta)
      cccos=c*cos(beta)
      aaacos=aa/cos(beta)
      cccos=c*cos(beta)
这段程序有什么作用    十分感谢!
回复 不支持

使用道具 举报

 楼主| 发表于 2011-12-26 09:42:56 | 显示全部楼层 来自 上海
lidongxu2006 发表于 2011-12-26 09:24
beta=0.0
      aacos=a/cos(beta)
      cccos=c*cos(beta)

这个本来是为了计算如果焊接电弧的中心线和焊接工件表面的法线之间有偏角beta的时候用的,一般情况下两者是基本重合的,beta=0。不过这个例子中没有用到,呵呵
回复 不支持

使用道具 举报

发表于 2011-12-26 11:22:51 | 显示全部楼层 来自 安徽合肥
JingheSu 发表于 2011-12-26 09:42
这个本来是为了计算如果焊接电弧的中心线和焊接工件表面的法线之间有偏角beta的时候用的,一般情况下两者 ...

谢谢楼主
回复 不支持

使用道具 举报

发表于 2011-12-26 15:03:13 | 显示全部楼层 来自 安徽合肥
deta=theta-phi,为什么有这一步啊     此外,abaqus调用子程序的具体过程怎样啊?           谢谢楼主!
回复 不支持

使用道具 举报

 楼主| 发表于 2011-12-26 16:17:42 | 显示全部楼层 来自 上海
lidongxu2006 发表于 2011-12-26 15:03
deta=theta-phi,为什么有这一步啊     此外,abaqus调用子程序的具体过程怎样啊?           谢谢楼主! ...

具体过程,呵呵,这个我还真不懂。
deta=theta-phi
theta是当前单元的位置,phi是当前热源中心的位置,减一下就是为了确定当前单元到底处于双椭球的那一个半球。
回复 不支持

使用道具 举报

发表于 2011-12-26 18:48:18 | 显示全部楼层 来自 安徽合肥
JingheSu 发表于 2011-12-26 16:17
具体过程,呵呵,这个我还真不懂。
deta=theta-phi
theta是当前单元的位置,phi是当前热源中心的位置,减 ...

谢谢楼主    呵呵                        
回复 不支持

使用道具 举报

发表于 2011-12-26 19:10:11 | 显示全部楼层 来自 安徽合肥
JingheSu 发表于 2011-12-26 16:17
具体过程,呵呵,这个我还真不懂。
deta=theta-phi
theta是当前单元的位置,phi是当前热源中心的位置,减 ...

版主,还想请教您个问题  ,就是上面的例子中您是怎样定义的单元生死呢?具体怎么实现的啊    非常感谢!
回复 不支持

使用道具 举报

 楼主| 发表于 2011-12-26 21:12:59 | 显示全部楼层 来自 北京
lidongxu2006 发表于 2011-12-26 19:10
版主,还想请教您个问题  ,就是上面的例子中您是怎样定义的单元生死呢?具体怎么实现的啊    非常感谢! ...

版内搜索model change
回复 不支持

使用道具 举报

发表于 2011-12-26 23:34:12 | 显示全部楼层 来自 广东深圳
我做了测试,发现算不完,一直出现warning,说是应变已经达到了50倍。

点评

奇怪,我的算完了,6.10+fortran11+vs08环境  发表于 2011-12-27 08:45
回复 不支持

使用道具 举报

发表于 2011-12-27 09:02:08 | 显示全部楼层 来自 广东深圳
我的是win7 32 + 6.10+fortran11+vs2008
回去再试试看。
回复 不支持

使用道具 举报

发表于 2011-12-27 14:01:51 | 显示全部楼层 来自 广东深圳
在笔记本上算了3个小时了,还在算。大概要算多久?
我的配置:xp + 4G内存 + abaqus6.10+vs2008+fortran11

点评

我自己电脑上xp + 2G内存 + abaqus6.10+vs2008+fortran11大概算了有6个小时,你可以看看sta文件情况,不是一下就知道计算速度了?  发表于 2011-12-27 15:12
回复 不支持

使用道具 举报

发表于 2011-12-27 15:41:01 | 显示全部楼层 来自 广东深圳
本帖最后由 smallmj 于 2011-12-27 15:41 编辑

这是我的warning信息

OUTPUT AT EXACT, PREDEFINED TIME POINTS WAS REQUESTED IN THIS STEP. Abaqus MIGHT ADJUST THE TIME INCREMENTATION BUT WILL NOT VIOLATE THE MINIMUM TIME INCREMENT ALLOWED IN THE STEP IN ORDER TO WRITE OUTPUT AT THOSE POINTS. CONSEQUENTLY, IN SOME CASES Abaqus MIGHTWRITE OUTPUT IMMEDIATELY AFTER THE SPECIFIED TIME POINTS, INSTEAD OF AT THE TIME POINTS. IN ADDITION, THE NUMBER OF INCREMENTS REQUIRED TO COMPLETE THE STEP WILL IN GENERAL INCREASE.

There are 2 unconnected regions in the model.

There is zero FORCE everywhere in the model based on the default criterion. please check the value of the average FORCE during the current iteration to verify that the FORCE is small enough to be treated as zero. if not, please use the solution controls to reset the criterion for zero FORCE.

There is zero FORCE everywhere in the model based on the default criterion. please check the value of the average FORCE during the current iteration to verify that the FORCE is small enough to be treated as zero. if not, please use the solution controls to reset the criterion for zero FORCE.

The strain increment has exceeded fifty times the strain to cause first yield at 20 points

The strain increment has exceeded fifty times the strain to cause first yield at 54 points

The strain increment has exceeded fifty times the strain to cause first yield at 4 points

The strain increment has exceeded fifty times the strain to cause first yield at 6 points

The strain increment has exceeded fifty times the strain to cause first yield at 118 points

The strain increment has exceeded fifty times the strain to cause first yield at 10 points

The strain increment has exceeded fifty times the strain to cause first yield at 28 points

The strain increment has exceeded fifty times the strain to cause first yield at 2 points

The strain increment has exceeded fifty times the strain to cause first yield at 42 points

The strain increment has exceeded fifty times the strain to cause first yield at 14 points

The strain increment has exceeded fifty times the strain to cause first yield at 140 points

The strain increment has exceeded fifty times the strain to cause first yield at 28 points

The strain increment has exceeded fifty times the strain to cause first yield at 22 points
回复 不支持

使用道具 举报

发表于 2011-12-27 15:51:48 | 显示全部楼层 来自 广东深圳
刚开始用ABAQUS,还不怎么会用,呵呵。
回复 不支持

使用道具 举报

发表于 2011-12-27 15:57:18 | 显示全部楼层 来自 广东深圳
还有,我算到13.4秒的时候已经57个step了,我看楼主给的图中算完15秒才32个step。
我就是把楼主给的文件下载下来后计算的。不知道为啥不一样。
回复 不支持

使用道具 举报

 楼主| 发表于 2011-12-27 16:00:14 | 显示全部楼层 来自 上海
smallmj 发表于 2011-12-27 15:57
还有,我算到13.4秒的时候已经57个step了,我看楼主给的图中算完15秒才32个step。
我就是把楼主给的文件下 ...

总共应该有65个step左右,呵呵。
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-4-27 00:17 , Processed in 0.051760 second(s), 18 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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