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

登录解锁更多功能

还没有账号?立即注册

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

新人求问大神一些问题

浏览: 1979 回答: 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,按理来说高于汽化温度的不是应该都杀死了吗? 是我两种不同材料的生死

推荐阅读

转子旋转的周期性模型-水冷电机散热仿真

转子旋转的周期性模型-水冷电机散热仿真

技术邻小李 技术邻小李
¥100
ABAQUS热切割数值模拟

ABAQUS热切割数值模拟

hitliuyong hitliuyong
¥99
Altair SimSolid 求解精度控制和对标网络研讨会

Altair SimSolid 求解精度控制和对标网络研讨会

ALTAIR ALTAIR
免费
Altair基于仿真驱动的制造工艺解决方案

Altair基于仿真驱动的制造工艺解决方案

ALTAIR ALTAIR
免费
HBK声学测量基础

HBK声学测量基础

HBK声学与振动 HBK声学与振动
免费
Mathcad基本操作

Mathcad基本操作

HYT HYT
¥89
后排座椅和安全带固定点强度分析

后排座椅和安全带固定点强度分析

Crisby_Vectory_TrHo Crisby_Vectory_TrHo
¥80
workbench纯六面体网格划分

workbench纯六面体网格划分

八号楼停水 八号楼停水
¥100
BOX-3D变宽斜钢箱梁绘图及功能更新演示

BOX-3D变宽斜钢箱梁绘图及功能更新演示

敦樸DUNPU 敦樸DUNPU
免费
新四化浪潮下,传统汽车人如何积极转型

新四化浪潮下,传统汽车人如何积极转型

夸克 夸克
免费
Cadence Fidelity 水泵水轮机CFD模拟解决方案和应用

Cadence Fidelity 水泵水轮机CFD模拟解决方案和应用

Cadence楷登 Cadence楷登
免费
abaqus_python提取odb中的最大值结果3

abaqus_python提取odb中的最大值结果3

tanquery tanquery
免费
十一特惠--Adams Car 模型建立的相关硬点及运动副

十一特惠--Adams Car 模型建立的相关硬点及运动副

仿真分析 仿真分析
¥9.9
ANSYS新版本功能速递: HFSS微放电仿真

ANSYS新版本功能速递: HFSS微放电仿真

Ansys中国 Ansys中国
免费
0#块箱梁托架法施工结构模拟

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

yudachuan1105 yudachuan1105
¥200
基于Hypermesh与ABAQUS联合仿真的摄像头模组跌落分析——工作创建以及检查

基于Hypermesh与ABAQUS联合仿真的摄像头模组跌落分析——工作创建以及检查

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

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

地下结构设计
¥40
5G终端天线仿真设计方法及其应用

5G终端天线仿真设计方法及其应用

Ansys中国 Ansys中国
免费
abaqus三维切削数值模拟(sph法)

abaqus三维切削数值模拟(sph法)

abaquser abaquser
¥30
Python面向对象的有限元编程

Python面向对象的有限元编程

weilai00 weilai00
¥60