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

为什么用函数编辑高斯函数无法加载轨迹为圆的热源,请大家指导!

[复制链接]
发表于 2016-9-6 10:33:59 | 显示全部楼层 |阅读模式 来自 天津
/FILNAME,weldingstress,0
/PREP7  
!定义焊接参数
L=1E-1      !焊件的长度
H=6E-3      !焊件的高度
W1=360    !整个圆360度
w=5        !焊接角速度
R=0.007       !电弧有效加热半径
Qm=4e7 !加热斑点中心最大热流密度
!*  
ET,1,PLANE55
ET,2,SOLID70


MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000   
MPTEMP,6,1500   
MPTEMP,7,1700   
MPTEMP,8,2500
   
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,40
MPDATA,KXX,1,,55
  
MPDATA,DENS,1,,7820
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090

MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,780
MPDATA,C,1,,820   

MPDATA,EX,1,,2.05E11
MPDATA,EX,1,,1.87E11
MPDATA,EX,1,,1.5E11
MPDATA,EX,1,,0.7E11
MPDATA,EX,1,,0.2E11
MPDATA,EX,1,,0.19E2
MPDATA,EX,1,,0.18E2
MPDATA,EX,1,,0.12e2
  
MPDATA,PRXY,1,,0.28
MPDATA,PRXY,1,,0.29
MPDATA,PRXY,1,,0.31
MPDATA,PRXY,1,,0.35
MPDATA,PRXY,1,,0.4  
MPDATA,PRXY,1,,0.45
MPDATA,PRXY,1,,0.48
MPDATA,PRXY,1,,0.5

UIMP,1,REFT,,,20
MPDATA,ALPX,1,,1.1e-5   
MPDATA,ALPX,1,,1.22e-5  
MPDATA,ALPX,1,,1.39e-5  
MPDATA,ALPX,1,,1.48e-5  
MPDATA,ALPX,1,,1.34e-5  
MPDATA,ALPX,1,,1.33e-5  
MPDATA,ALPX,1,,1.32e-5  
MPDATA,ALPX,1,,1.31e-5  


TB,BISO,1,6,2,  
TBTEMP,20   
TBDATA,,220e6,0,,,,
TBTEMP,250  
TBDATA,,175e6,0,,,,
TBTEMP,500  
TBDATA,,80e6,0,,,,  
TBTEMP,750
TBDATA,,40E6,0,,,,   
TBTEMP,1000
TBDATA,,10E6,0,,,,   
TBTEMP,1500
TBDATA,,1E-5,0,,,,   


K,1,L,L,0
K,2,L,-L,0
K,3,-L,-L,0
K,4,-L,L,0


A,1,2,3,4

ESIZE,0.005
AMESH,1

TYPE,   2   
EXTOPT,ESIZE,2,0,   
EXTOPT,ACLEAR,1
!*  
EXTOPT,ATTR,1,0,0   
REAL,_Z4
ESYS,0  
!*  
VOFFST,1,H, ,  

NUMMRG,ALL, , , ,LOW

/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS   
*SET,_FNCNAME,'zhu1yzzb'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),5   
*SET,_FNC_C2(1),0.007   
*SET,_FNCCSYS,1
! /INPUT,H:\Desktop\csys1zhugui.func,,,1
*DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS%  
! Begin of equation: 4e7*exp(-3*((0.5)^2+(w*{time})^2)/r^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   1, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.5, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 17, -3   
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, 17, 3, 1
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -1, 17, -3  
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 18, 17, -1  
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -3, 4, -2   
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 4e7, 0, 0, -1  
*SET,%_FNCNAME%(0,19,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 4e7*exp(-3*((0.5)^2+(w*{time})^2)/r^2)   
TUNIF,20, !定义初始温度

SFA,2,1,CONV,30,20
SFA,3,1,CONV,30,20

SFA,3,1,CONV,30,20
SFA,4,1,CONV,30,20
SFA,5,1,CONV,30,20
!施加高斯热源  
SFA,1,1,HFLUX, %zhu1yzzb%   

OUTRES,ALL,ALL,

TIME,W1/W !设置求解时间
AUTOTS,-1   
NSUBST,50,50,50
KBC,0   
TSRES,ERASE
LSWRITE,1,  !写入载荷文件为1
!*  
TIME,80
AUTOTS,1
NSUBST,20,20,20
KBC,0   
!*  
TSRES,ERASE
LSWRITE,2,

TIME,100
AUTOTS,1
NSUBST,30,30,30
KBC,0   
!*  
TSRES,ERASE
LSWRITE,3,

TIME,1100
AUTOTS,1
NSUBST,105,105,105
KBC,0   
!*  
TSRES,ERASE
LSWRITE,4,
  
LSSOLVE,1,4,1,  !开始求解

我是用柱坐标csys1,w为角速度,w1是360圆
 楼主| 发表于 2016-9-7 16:56:07 | 显示全部楼层 来自 天津
Simdroid开发平台
打开ANSYS经典界面,在函数编辑器下创建如下函数:

高斯函数=Qm*exp(-3* {r^2/R^2})     R为加热半径   r为点到电弧加热的中心

回复 不支持

使用道具 举报

发表于 2016-9-13 08:29:47 | 显示全部楼层 来自 上海闵行区
  你的函数是时间和 r 的函数,问题是:
  1  r 是如何随着时间变化的?当你施加该函数时,没有定义过 r,因此 r 默认为 0;
  2  你只有一个载荷步,且是在其它设置之前施加了载荷,那么它将以时刻 0 时的 r 施加载荷,以后 r 是不会变化的。因此你需要将整个时间过程分解为若干载荷步,在每个载荷步中重新计算 r 和施加载荷,这样,载荷才会随着时间改变。
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-14 21:15:15 | 显示全部楼层 来自 天津
htbbzzg 发表于 2016-9-13 08:29
你的函数是时间和 r 的函数,问题是:
  1  r 是如何随着时间变化的?当你施加该函数时,没有定义过 r, ...

前辈您好,这个高斯热源函数是  定义半径和时间的关系,我现在想用柱坐标来计算,这样的话 在csys1中,X是半径  Y表示角度   所以我想用 角速度来确定时间
回复 不支持

使用道具 举报

发表于 2016-9-15 07:58:45 | 显示全部楼层 来自 上海闵行区
本帖最后由 htbbzzg 于 2016-9-15 08:00 编辑

  定义函数时,可以使用自变量 时间,但不能使用角速度作为自变量,所以你需要改为使用时间来表示角速度。  另外,上次写的第二点你考虑过没有?
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-15 10:17:48 | 显示全部楼层 来自 天津
htbbzzg 发表于 2016-9-15 07:58
定义函数时,可以使用自变量 时间,但不能使用角速度作为自变量,所以你需要改为使用时间来表示角速度。  ...

前辈,祝您中秋节快乐,多谢指点,目前第二点我还不会改造,这个是我找了一个例子,原来是施加的高斯热源走直线,是可以直接运行的,我就在例子基础上面直接改为走 路径为圆的  曲线,所以下载的载荷步这块,我还没有整
回复 不支持

使用道具 举报

发表于 2016-9-15 11:03:52 | 显示全部楼层 来自 上海闵行区
  定义函数时,在函数定义对话框中可以选择使用直角坐标系,也可以使用圆柱坐标系。至于后面走直线,还是走圆弧,与函数定义没有直接关系,是由你自行定义的。关键是准确确定函数加载的位置。
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-15 19:25:45 | 显示全部楼层 来自 天津
htbbzzg 发表于 2016-9-15 11:03
定义函数时,在函数定义对话框中可以选择使用直角坐标系,也可以使用圆柱坐标系。至于后面走直线,还是走 ...

前辈,我一直在迷惑,在函数编辑器中得x y 代表的是什么,难道不是函数施加的区域吗?
这个高斯热源函数中,r^2  代表点到电弧的加热中心距离,不应该控制  r  来确定路径吗?

您说的我自行定义,函数准确加载的位置 指的是什么意思?  我在参考  原来走直线的例子时,那上面并没有选择节点,而是直接  控制 函数中得  x  y  time   ,直接施加在面上去
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-15 19:27:29 | 显示全部楼层 来自 天津
liu19930001 发表于 2016-9-15 19:25
前辈,我一直在迷惑,在函数编辑器中得x y 代表的是什么,难道不是函数施加的区域吗?
这个高斯热源函数 ...

他在确定  x  y  time  的时候,就是根据  距离原心的距离确定的,前辈,你有没有关于函数加载这一块的资料推荐,我现在找不到相关的资料,手头就这一份资料
回复 不支持

使用道具 举报

发表于 2016-9-16 10:09:07 | 显示全部楼层 来自 上海闵行区
  你的关键问题是:施加高斯热源的命令实在所有载荷步之前,因此它是按照当时的 x,y 和 time 来施加的,后面不再变化:
    !施加高斯热源  
    SFA,1,1,HFLUX, %zhu1yzzb%   
    因此,你需要将整个加载过程切分为若干小的时间步,然后使用循环语句,逐步进行计算,每改变一次时间,就施加一次高斯热源,这样高斯热源才会动起来。
    至于函数加载的资料,网上很多,这里给你几份看看吧:
   
   
   
   
   
   

本帖子中包含更多资源

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

×
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-16 11:12:34 | 显示全部楼层 来自 天津
htbbzzg 发表于 2016-9-16 10:09
你的关键问题是:施加高斯热源的命令实在所有载荷步之前,因此它是按照当时的 x,y 和 time 来施加的,后 ...

嗯,多谢您的指点,我就沿着您的思路  检查下这个命令流,还有前辈,您的邮箱方便留一下吗?这两天我没有用函数加载,直接采用施加热流密度,让热源动了起来,我想把命令流给你发一下
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-16 11:25:11 | 显示全部楼层 来自 天津
liu19930001 发表于 2016-9-16 11:12
嗯,多谢您的指点,我就沿着您的思路  检查下这个命令流,还有前辈,您的邮箱方便留一下吗?这两天我没有 ...

前辈,我就是根据这个附件  改写的命令流,他原来是走的直线,也是在  载荷步之前就是时间了高斯函数,为什么这个例子就可以移动

还有,我的邮箱是1039049732@qq.com  

本帖子中包含更多资源

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

×
回复 不支持

使用道具 举报

发表于 2016-9-17 07:25:34 | 显示全部楼层 来自 上海闵行区
  我的邮箱:htbbzzg@163.com
  不过,我最近家里事情很多,未必能及时给你回信。
回复 不支持

使用道具 举报

 楼主| 发表于 2016-9-17 08:59:18 | 显示全部楼层 来自 天津
htbbzzg 发表于 2016-9-17 07:25
我的邮箱:
  不过,我最近家里事情很多,未必能及时给你回信。

好的,那前辈您先忙,我这几天好好整理一下,然后把原文件文件什么的,一起发过去
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-3-29 09:43 , Processed in 0.037935 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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