为了正常的体验网站,请在浏览器设置里面开启Javascript功能!

将ANSYS导入ABAQUS

2019-06-14 10页 doc 25KB 46阅读

用户头像

is_219945

暂无简介

举报
将ANSYS导入ABAQUSSYS软件的参数化建模(APDL)极其方便,而ABAQUS卓越的非线性计算功能使其成为有限元软件中的贵族,如果能整合两者的优势,有限元模拟计算就会相当高效。下面就将ANSYS模型如何导入ABAQUS简单交流一下。(其中部分内容参考了SIMWE上的精华贴) 首先在ANSYS编写APDL语言输出模型的节点和单元数据文件。如果直接用ANSYS中的WRITE NODE FILE和WRITE ELEM FILE得到的节点和单元数据文件,数据之间只有空格,没有逗号,不符合ABAQUS的INPUT文件格式要求。输出节点信息的APDL如下...
将ANSYS导入ABAQUS
SYS软件的参数化建模(APDL)极其方便,而ABAQUS卓越的非线性计算功能使其成为有限元软件中的贵族,如果能整合两者的优势,有限元模拟计算就会相当高效。下面就将ANSYS模型如何导入ABAQUS简单交流一下。(其中部分内容参考了SIMWE上的精华贴) 首先在ANSYS编写APDL语言输出模型的节点和单元数据文件。如果直接用ANSYS中的WRITE NODE FILE和WRITE ELEM FILE得到的节点和单元数据文件,数据之间只有空格,没有逗号,不符合ABAQUS的INPUT文件格式。输出节点信息的APDL如下: Allsel,all !选中所有 ! 输出节点 *GET,NNode,NODE,,COUNT, , , , ! 得到当前模型中的总节点数 *CFOPEN,ansystoabaqus,inp *DO,I,0,NNode *VWRITE,Chrval(i),',',NX(I),',',NY(I),',',NZ(I) (A8,A1,F10.5,A1,F10.5,A1,F10.5) *ENDDO Allsel,all 输出单元信息时要弄清楚单元类型,编写APDL对每种类型的每个单元输出单元编号,各节点的编号。常见的SOLID65和SOLID45、SOLID95单元的APDL输出命令如下: *GET,NumberofNode,NODE,,COUNT, , , , ESEL,S,ENAME,,65 *GET,NElem,ELEM,,COUNT, , , , ! 得到当前模型中的总单元数 !对单元集进行循环 *GET,nd,ELEM,,NUM,MIN, , , , *DO,I,1,NElem !得到当前单元的类型 *GET,ENAME,ELEM,nd,ATTR,ENAM !如果是65号(65号单元) *IF,ENAME,EQ,65,THEN !得到该单元的节点编号 *GET,EN1,ELEM,nd,NODE,1 *GET,EN2,ELEM,nd,NODE,2 *GET,EN3,ELEM,nd,NODE,3 *GET,EN4,ELEM,nd,NODE,4 *GET,EN5,ELEM,nd,NODE,5 *GET,EN6,ELEM,nd,NODE,6 *GET,EN7,ELEM,nd,NODE,7 *GET,EN8,ELEM,nd,NODE,8 *VWRITE,Chrval(nd),',',Chrval(EN1),',',Chrval(EN2),',',Chrval(EN3),',',Chrval(EN4),',',Chrval(EN5),', ',Chrval(EN6),',',Chrval(EN7),',',Chrval(EN8) (A8,8(A1,A8)) nd = ELnext(nd) *END IF *ENDDO Allsel,all *GET,NumberofNode,NODE,,COUNT, , , , ESEL,S,ENAME,,45 !solid45单元8节点六面体 *GET,NElem,ELEM,,COUNT, , , , ! 得到当前模型中的总单元数 !对单元集进行循环 *GET,nd,ELEM,,NUM,MIN, , , , *DO,I,1,NElem *GET,ENAME,ELEM,nd,ATTR,ENAM !得到当前单元的类型 *IF,ENAME,EQ,45,THEN !如果是45号(45号单元) !得到该单元的节点编号 *GET,EN1,ELEM,nd,NODE,1 *GET,EN2,ELEM,nd,NODE,2 *GET,EN3,ELEM,nd,NODE,3 *GET,EN4,ELEM,nd,NODE,4 *GET,EN5,ELEM,nd,NODE,5 *GET,EN6,ELEM,nd,NODE,6 *GET,EN7,ELEM,nd,NODE,7 *GET,EN8,ELEM,nd,NODE,8 *VWRITE,Chrval(nd),',',Chrval(EN1),',',Chrval(EN2),',',Chrval(EN3),',',Chrval(EN4),',',Chrval(EN5),', ',Chrval(EN6),',',Chrval(EN7),',',Chrval(EN8) (A8,8(A1,A8)) nd = ELnext(nd) *END IF *ENDDO *VWRITE ('*ELEMENT,TYPE=C3D15 ,ELSET=Esolid2') Allsel,all *GET,NumberofNode,NODE,,COUNT, , , , ESEL,S,ENAME,,95 !solid95单元20节点六面体 *GET,NElem,ELEM,,COUNT, , , , !得到当前模型中的总单元数 !对单元集进行循环 *GET,nd,ELEM,,NUM,MIN, , , , ! 得到当前模型中的最小单元号 *DO,I,1,NElem *GET,ENAME,ELEM,nd,ATTR,ENAM !得到当前单元的类型 *IF,ENAME,EQ,95,THEN !如果是95号(95号单元) !得到该单元的节点编号 *GET,EN1,ELEM,nd,NODE,1 *GET,EN2,ELEM,nd,NODE,2 *GET,EN3,ELEM,nd,NODE,3 *GET,EN4,ELEM,nd,NODE,4 *GET,EN5,ELEM,nd,NODE,5 *GET,EN6,ELEM,nd,NODE,6 *GET,EN7,ELEM,nd,NODE,7 *GET,EN8,ELEM,nd,NODE,8 *GET,EN9,ELEM,nd,NODE,9 *GET,EN10,ELEM,nd,NODE,10 *GET,EN11,ELEM,nd,NODE,11 *GET,EN12,ELEM,nd,NODE,12 *GET,EN13,ELEM,nd,NODE,13 *GET,EN14,ELEM,nd,NODE,14 *GET,EN15,ELEM,nd,NODE,15 *GET,EN16,ELEM,nd,NODE,16 *GET,EN17,ELEM,nd,NODE,17 *GET,EN18,ELEM,nd,NODE,18 *GET,EN19,ELEM,nd,NODE,19 *GET,EN20,ELEM,nd,NODE,20 ENN1=EN1 ENN2=EN2 ENN3=EN3 ENN4=EN5 ENN5=EN6 ENN6=EN7 x1=Nx(EN1) y1=Ny(EN1) Z1=Nz(EN1) x2=Nx(EN2) y2=Ny(EN2) Z2=Nz(EN2) ENN7=NumberofNode+(I-1)*5+1 ENN8=EN10 ENN9=EN12 x5=Nx(EN5) y5=Ny(EN5) Z5=Nz(EN5) x6=Nx(EN6) y6=Ny(EN6) Z6=Nz(EN6) ENN10=NumberofNode+(I-1)*5+2 ENN11=EN14 ENN12=EN16 ENN13=NumberofNode+(I-1)*5+3 ENN14=NumberofNode+(I-1)*5+4 ENN15=NumberofNode+(I-1)*5+5 x3=Nx(EN3) y3=Ny(EN3) Z3=Nz(EN3) x7=Nx(EN7) y7=Ny(EN7) Z7=Nz(EN7) *VWRITE,Chrval(ENN7),',',(X1+X2)/2,',',(y1+y2)/2,',',(z1+z2)/2 (A8,3(A1,F10.5)) *VWRITE,Chrval(ENN10),',',(X5+X6)/2,',',(y5+y6)/2,',',(z5+z6)/2 (A8,3(A1,F10.5)) *VWRITE,Chrval(ENN13),',',(X5+X1)/2,',',(y5+y1)/2,',',(z5+z1)/2 (A8,3(A1,F10.5)) *VWRITE,Chrval(ENN14),',',(X6+X2)/2,',',(y6+y2)/2,',',(z6+z2)/2 (A8,3(A1,F10.5)) *VWRITE,Chrval(ENN15),',',(X3+X7)/2,',',(y3+y7)/2,',',(z3+z7)/2 (A8,3(A1,F10.5)) nd =ELnext(nd) *END IF *ENDDO !输出实体单元 *VWRITE ('*ELEMENT,TYPE=C3D8,ELSET=Esolid1') Allsel,all 运行以上命令后会在ANSYS的工作目录下得到名为ansystoabaqus.inp的文件,里面包含了模型的节点和单元信息。接下来就是编写ABAQUS的inp文件了。格式如下 *Heading ** Job name: imput Model name: ^^^^ ** Generated by: Abaqus/CAE 6.9-1 *Preprint, echo=NO, model=NO, history=NO, contact=NO ** ** PARTS ** *Part, name=^^^^^^^ *End Part ** ** ** ASSEMBLY ** *Assembly, name=Assembly ** *Instance, name=^^^^^, part=^^^^^^ *Node !!!!!!!!!!!!!!!!!!!!!!!!!ansystoabaqus.inp中的节点文件!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *Element,type=^^^^^ !!!!!!!!!!!!!!!!!!!!!!!!!ansystoabaqus.inp中的单元文件!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 最后在ABAQUS中点主菜单FILE\IMPORT\MODEL,选择要导入的INP文件,在窗口顶部环境栏的MODEL下拉列表中,就会出现与此INP文件同名的模型。
/
本文档为【将ANSYS导入ABAQUS】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。 本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。 网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。

历史搜索

    清空历史搜索