顶部banner
Logo 折叠
全网首发仿真工程师项目实训
下载app
返回旧版
首页 成长助手 小邻学院 社区 发现
职业认证 企业服务 行业会议

登录解锁更多功能

还没有账号?立即注册

技术邻
技术引领职场价值
电话
0571-86682823
商务合作
service@jishulink.com
  • 全部  > 
  • CAE仿真

新人求问大神一些问题

浏览: 2019 回答: 5

当我把宏文件导入到ANSYS后出现下面问题

CXTKJ(LOIB}UC)]E$AU2LHP.png

附上宏文件命令:

/prep7

*del,all

cmsel,s,longgu_line

pai=3.1415926

e_tan=2.1e5

f_y=235

lmd_y=1.05

lmd_z=1.05

*get,lmin,line,0,num,min

*get,lnum,line,0,count

*get,lmax,line,0,num,max

ll=lmin

*dim,changxibi,,lunm,4

esll,s

*get,enum,elem,,count

*dim,result,,enum,2

*dim,result_la,,enum,2

*dim,result_ya,,enum,2

num_of_line=0

num_of_line_la=0

num_of_line_ya=0

num_of_ele=0

num_of_ele_la=0

num_of_ele_ya=0

*do,m,1,lnum,1

num_of_line=num_of_line+1

lsel,s,line,,ll,ll,1

*get,length_of_line,line,ll,leng

changxibi(num_of_line,4)=length_of_line

esll,s

/post1

etable,sdi,smisc,31

etable,sbi_yt,smisc,32

etable,sbi_yb,smisc,33

etable,sbi_zt,smisc,34

etable,sbi_zb,smisc,35

etable,axial_f,smisc,1

smax,sbi_y_max,sbi_yt,sbi_yb,1,1

smin,sbi_y_min,sbi_yt,sbi_yb,1,1

smax,sbi_z_max,sbi_zt,sbi_zb,1,1

smin,sbi_z_min,sbi_zt,sbi_zb,1,1

*get,emin,elem,0,num,min

*get,enum,elem,,count

*get,emax,elem,0,num,max

ee=emin

*get,sect_num,elem,ee,attr,secn

*get,area_sec,secp,sect_num,prop,area

*get,iz_sec,secp,sect_num,prop,izz

*get,iy_sec,secp,sect_num,prop,iyy

iyy=sqrt(iy_sec/area_sec)

izz=sqrt(iz_sec/area_sec)

lmd_cxb_y=length_of_line/iyy

lmd_cxb_z=length_of_line/izz

ne_y=(pai**2*e_tan*area_sec)/(1.1*lmd_cxb_y**2)

ne_z=(pai**2*e_tan*area_sec)/(1.1*lmd_cxb_z**2)

changxibi(num_of_line)=ll

lmd_cxb=lmd_cxb_y

*if,lmd_cxb_y,le,lmd_cxb_z,then

*endif

changxibi(num_of_line,2)=lmd_cxb

*get,str_zhou,elem,ee,etab,sdi,,,

*if,str_zhou,gt,1e-8,then

        mark=1

        changxibi(num_of_line,3)=1

        num_of_line_la=num_of_line_la+1

                   fai_wd=0

*elseif,str_zhou,lt,1e-8,then

        mark=2

        changxibi(num_of_line,3)=-1

        num_of_line_ya=num_of_line_ya+1

        lmd_n_y=(lmd_cxb_y/pai)*sqrt(f_y/e_tan)

        *if,lmd_n_y,le,0.215,then

                  fai_wd_y=1-0.65*lmd_n_y**2

        *elseif,lmd_n_y,gt,0.215,then

                  changshu=0.965+0.3*lmd_n_y+lmd_n_y**2

                  fai_wd_y=(changshu-sqrt(changshu**2-4*lmd_n_y**2))/(2*lmd_n_y**2)

        *endif

        lmd_n_z=(lmd_cxb_z/pai)*sqrt(f_y/e_tan)

        *if,lmd_n_z,le,0.215,then

                  fai_wd_z=1-0.65*lmd_n_z**2

        *elseif,lmd_n_z,gt,0.215,then

                  changshu=0.965+0.3*lmd_n_z+lmd_n_z**2

                  fai_wd_z=(changshu-sqrt(changshu**2-4*lmd_n_z**2))/(2*lmd_n_z**2)

        *endif

*endif

fai_wd=fai_wd_y

*if,fai_wd_z,le,fai_wd_y,then

fai_wd=fai_wd_z

*endif

changxibi(num_of_line,4)=fai_wd

*if,mark,eq,1,then

*do,n,1,enum,1

*get,stress_zhou,elem,ee,etab,sdi,,,

*get,stress_wq_y,elem,ee,etab,sbi_y_max,,,

*get,stress_wq_z,elem,ee,etab,sbi_z_max,,,

num_of_ele=num_of_ele+1

num_of_ele_la=num_of_ele_la+1

result(num_of_ele,1)=ee

result(num_of_ele,2)=(stress_zhou+sress_wq_y/lmd_y+sress_wq_z/lmd_z)

result_la(num_of_ele_la,1)=ee

result_la(num_of_ele_la,2)=(stress_zhou+sress_wq_y/lmd_y+sress_wq_z/lmd_z)

ee=elnext(ee)

*enddo

*elseif,mark,eq,2,then

*do,n,1,enum,1

*get,stress_zhou,elem,ee,etab,sdi,,,

*get,stress_wq_y,elem,ee,etab,sbi_y_min,,,

*get,stress_wq_z,elem,ee,etab,sbi_z_min,,,

*get,zhouli,elem,ee,etab,axial_f,,,

num_of_ele=num_of_ele+1

num_of_ele_ya=num_of_ele_ya+1

result(num_of_ele,1)=ee

result(num_of_ele,2)=(stress_zhou/fai_wd+stress_wq_y/(lmd_y*(1.0-0.8*zhouli/ne_y))+stress_wq_z/(lmd_z*(1.0-0.8*zhouli/ne_z)))

result_ya(num_of_ele_ya,1)=ee

result_ya(num_of_ele_ya,2)=(stress_zhou/fai_wd+stress_wq_y/(lmd_y*(1.0-0.8*zhouli/ne_y))+stress_wq_z/(lmd_z*(1.0-0.8*zhouli/ne_z)))

ee=elnext(ee)

*enddo

*endif

cmsel,s,longgu_line

ll=lsnext(ll)

*enddo


大神们能帮我看看错在哪吗?改了好几天,依然毫无头绪!!!不胜感激啊!!

ANSYS CAE 结构CAE

全部回答 (1)

默认 最新
力争ANSYS大神 2016年9月15日
没有大神吗@ANSYS专家
2016年9月15日
评论 4 点赞

相似问题

查看全部
  • ANSYS APDL语言编程 循环? 2个回答

    我正在编写一个APDL命令,程序准备实现的功能如下: 1,由于是大变形问题,所以采用自适应时间步长 2,在每个substep停下,将发生位移的接触单元上的载荷去除 3,继续下一个substep,直到所求的时间为止。 但是我根据这个想法遍了一段命令,却并不能实现我想的功能,想请问各位专家,问题出在哪里了? 命令如下: *DO,I,1,3000 *IF,I,EQ,1 /SOLU ANTYPE,0 TI

  • 求解释ANSYS命令流? 7个回答

    finish /prep7 *set,strd,90e6 esel,none cm,bldam,elem cm,bl1dam,elem cm,bl2dam,elem cm,bl3dam,elem cm,bl4dam,elem cm,bl5dam,elem cm,bl6dam,elem cm,bl7dam,elem cm,bl8dam,elem cm,bl9dam,elem cm,bl0dam,el

  • Ansys apdl 结合生死单元仿真激光烧蚀纤维复合材料,表面温度过高? 1个回答

    考虑相变(直接汽化),手动添加了大概的焓值,表面基质相汽化温度3460°C,纤维相汽化温度4000°C,模型如图 加载激光热流密度,生死单元求解,网格最细划到了50等分(100等分不求解,不知道为啥) 热流密度大的话,基质和纤维烧了一部分,但温度最高达到了10000°C 热流密度小的话,只烧了基质,但表面温度最高还是有5000°C,按理来说高于汽化温度的不是应该都杀死了吗? 是我两种不同材料的生死

推荐阅读

用户自定义场变量子程序USDFLD从入门到高级 (如何通过USDFLD实现本构模型参数随状态改变)

用户自定义场变量子程序USDFLD从入门到高级 (如何通过USDFLD实现本构模型参数随状态改变)

长江学渣 长江学渣
¥48
0#块箱梁托架法施工结构模拟

0#块箱梁托架法施工结构模拟

yudachuan1105 yudachuan1105
¥200
ABAQUS钢管混凝土柱温度场及耐火性能分析(未完)

ABAQUS钢管混凝土柱温度场及耐火性能分析(未完)

地下结构设计
¥40
前沿技术!大数据分析及人工智能在优化软件中的应用

前沿技术!大数据分析及人工智能在优化软件中的应用

IDAJ中国 IDAJ中国
¥9.99
新一代智能头灯的动态设计评估与仿真

新一代智能头灯的动态设计评估与仿真

Ansys中国 Ansys中国
免费
UG培训第九课:自由曲面构造法

UG培训第九课:自由曲面构造法

luffy8610 luffy8610
¥20
应用ANSYS瞬态动力学法模拟啮合齿轮的高速转动

应用ANSYS瞬态动力学法模拟啮合齿轮的高速转动

夏日星空 夏日星空
¥35
一线科技工作者接受仿真咨询服务的全过程经验分享

一线科技工作者接受仿真咨询服务的全过程经验分享

技术邻直播 技术邻直播
免费
混凝土材性试块拉压数值模拟(ABAQUS通法建模初级案例3)

混凝土材性试块拉压数值模拟(ABAQUS通法建模初级案例3)

大平-结构工程 大平-结构工程
¥299
汽车仪表模具模流分析的实战讲解

汽车仪表模具模流分析的实战讲解

北卡 北卡
免费
UG有限元基础教程

UG有限元基础教程

moonshine🤓 moonshine🤓
免费
ABAQUS桁架结构强度分析

ABAQUS桁架结构强度分析

wj_2704 wj_2704
¥500
ADAS功能软件基础介绍

ADAS功能软件基础介绍

Alex王 Alex王
免费
abaqus lamb波传播分析

abaqus lamb波传播分析

abaquser abaquser
¥25
abaqus模拟桩(三维)全过程

abaqus模拟桩(三维)全过程

红日当空q3120210076 红日当空q3120210076
¥62
TurboTides 2024R2全新版本发布会--智能算法驱动的全自动优化透平机械集成设计平台

TurboTides 2024R2全新版本发布会--智能算法驱动的全自动优化透平机械集成设计平台

TurboTides TurboTides
免费
场景仿真加速智能网联开发测试进程

场景仿真加速智能网联开发测试进程

海克斯康设计与仿真 海克斯康设计与仿真
免费
ANSYS speos汽车信号灯案例实操教程

ANSYS speos汽车信号灯案例实操教程

IDAJ中国 IDAJ中国
¥9.9
ANSYS SI/PI/EMI&TI 2019 R3 新功能介绍

ANSYS SI/PI/EMI&TI 2019 R3 新功能介绍

Ansys中国 Ansys中国
免费
复杂电驱动系统动态实时功率测试

复杂电驱动系统动态实时功率测试

HBK测试与测量 HBK测试与测量
免费