24小时热门版块排行榜    

查看: 1315  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +5 yunziaaaaa 2026-03-01 6/300 2026-03-01 23:57 by ccp273206157
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 材料类求调剂 +10 wana_kiko 2026-02-28 12/600 2026-03-01 22:10 by 海嵙Y
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 12/600 2026-03-01 19:48 by 好好好1233
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 295求调剂 +7 19171856320 2026-02-28 7/350 2026-03-01 18:54 by 18137688336
[考研] 材料学调剂 +9 提神豆沙包 2026-02-28 11/550 2026-03-01 18:15 by ms629
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
请填处理意见