24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1331  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344 材料专业 求调剂211 无地域要求 +5 hualkop 2026-04-11 5/250 2026-04-11 23:13 by 852137818
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 268求调剂 +3 你好tg 2026-04-05 5/250 2026-04-11 16:45 by Delta2012
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 326分,一志愿沪9,求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 22408 一志愿双一流人工智能300分 四六级,数据分析国奖 +4 zzfeng123 2026-04-06 6/300 2026-04-07 21:02 by zzfeng123
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
信息提示
请填处理意见