24小时热门版块排行榜    

查看: 1275  |  回复: 5

vicly2013

铁虫 (小有名气)

[求助] 在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪? 已有2人参与

如题!!程序见附件!
我建立的模型如下图所示:

在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪?
模型2.jpg


在施加移动旋转高斯热源时,通过*get命令获取各单元的质心坐标。运行后,各坐标获取成功且不为0,但热源密度youbu=0(见下图). 不知哪里不对,请各位大神帮忙看看,非常感谢!!!
错误如下:在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪?-1
Error4.jpg


在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪?-2
Error5.jpg




在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪?-3
0.1s状态.jpg


在施加旋转高斯热源时总提示错误,但感觉程序没问题,不知错在哪?-4
0.2s状态.jpg

[ Last edited by vicly2013 on 2013-12-14 at 17:07 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 29-2.inp
  • 2013-12-14 17:05:46, 5.85 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
1592203609: 金币+2, 谢谢回帖! 2013-12-14 21:43:18
你可以检查一下循环中用到的 y 的值是否可能为 0:
*DO,i,1,max_time,1
    nsel,s,loc,y,0,ly1     !  节点 y 坐标为 0
    nsel,r,loc,z,lz1-lz2,lz1
    nsel,r,loc,x,(i-1)*lsize2,i*lsize2
    cm,eee,elem
    cmsel,s,eee
    cmlist,eee,1
    *get,emax,elem,,count
    *get,elem1,elem,,num,min
    *status
    *do,j,1,emax,1
        x=nex(elem1)
        y=ney(elem1)         !  此处 y 值是否可能为 0?
        z=nez(elem1)
        distance=sqrt((x-i*lsize2)**2+z**2)   !!!距离
        r=sqrt(log(ly1/y)/3)*r0               !!!定义某一界面处高斯旋转半径的大小  ! 如果 y 为 0,则出现被零除的错误。
2楼2013-12-14 21:18:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vicly2013

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by htbbzzg at 2013-12-14 21:18:34
你可以检查一下循环中用到的 y 的值是否可能为 0:
*DO,i,1,max_time,1
    nsel,s,loc,y,0,ly1     !  节点 y 坐标为 0
    nsel,r,loc,z,lz1-lz2,lz1
    nsel,r,loc,x,(i-1)*lsize2,i*lsize2
    cm,eee,el ...

Y值也不为0 啊。上图有*status命令列出的数据
3楼2013-12-14 22:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

你可以在有 "/y" 的语句之前添加语句:
    *if, abs(y), le, 1e-10, then
          输出相应语句中的主要参数 (如单元编号、单元重心坐标、以及循环变量的值,等)
    *endif
可以确认是否有 y=0 的情况及其发生的位置。
计算时可以使用 /output 语句将输出转换到一个文件中,以便出错时查询。
4楼2013-12-15 08:04:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vicly2013

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by htbbzzg at 2013-12-15 08:04:15
你可以在有 "/y" 的语句之前添加语句:
    *if, abs(y), le, 1e-10, then
          输出相应语句中的主要参数 (如单元编号、单元重心坐标、以及循环变量的值,等)
    *endif
可以确认是否有 y=0 的 ...

嗯,谢谢老师的回复。我尝试下。
5楼2013-12-15 15:43:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luo_kuijin

银虫 (小有名气)

【答案】应助回帖

引用回帖:
5楼: Originally posted by vicly2013 at 2013-12-15 15:43:34
嗯,谢谢老师的回复。我尝试下。...

您好楼主,您这个环形热源加载的问题解决了么?我现在也在做这个关于环形加载的不会,请楼主把您之前做的案例分享一下,谢谢。luo.kuijin@163.com
6楼2015-12-16 09:49:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vicly2013 的主题更新
信息提示
请填处理意见