【转载】ANSYS的APDL与C语言混合编程(实例)
原文地址:http://www.cnblogs.com/lyq105/archive/2010/05/04/1727557.html
本文讨论的不是利用C语言为ANSYS写扩展(或者说是用户子程序),而是讨论结合APDL对ANSYS的现有的计算程序进行复用,例如对ANSYS强大的前后处理程序,包括网格剖分,计算结果显示等,或者利用有限元计算部分的模块,使用ANSYS方便的边界条件加载,方程组求解等。
0. 为什么进行混合编程
既然ANSYS已经提供了APDL这么强大的编程语言,为什么还要用C语言来与APDL混合编程呢?下面的几个原因,首先是前处理建模的要求,APDL并不能应对实际中非常复杂的模型的构建,本文的出发点之一就是复合材料分析中的几何模型构建问题,由于材料中的增强体为随机的颗粒,每一次模拟中,颗粒都是随机生成的,使用APDL并不方便;其次就是利用ANSYS的计算资源的考虑,比如上述问题,就可以利用C语言来生成模型的数据,利用APDL来构建几何模型,以及剖分有限元网格,用于自己的计算程序;还有在一些问题中,其边界条件,载荷等也是通过程序算出来的话,也可以通过混合编程对问题求解;另外,可以发挥C语言的优势,比如说计算效率,现有的库等,总而言之,混合编程能更充分利用ANSYS的计算程序以及C语言的优势,并且可以实现计算机上的自动化计算。
1. 混合编程的原理
ANSYS与C语言混合编程的原理很简单,实际上就是用C语言生成ANSYS的apdl命令流文件,然后通过命令行调用ANSYS进行计算。这个过程同样可以使用其他的编程语言(比如FORTRAN,Python等)来进行,可以发挥各种语言的优势以及对现有的程序库的使用。
2. 如何使用C语言的实现
知道了混合编程的原理,那么就很简单了,其过程无非就是利用C语言打开一个文件,写入命令流,关闭文件,调用ANSYS进行计算。
文件操作需要使用的函数有
1 fp=fopen("mac.inp","w"); //打开文件
2 fprintf(fp, "APDL_Commands,%f\n",SomeData); //写入文件
3 fclose(fp); //关闭文件
调用ANSYS的命令行计算则使用
system("D:\\Progra~1\\AnsysI~1\\v100\\ANSYS\\bin\\intel\\ansys100 -b -i temp.inp -o mesh.log");
3. 实例
下面给几个APDL和C语言混合编程的实例
3.1 复合材料中的圆性颗粒网格剖分
图中的颗粒是C语言中随机生成的圆形,结果如图所示:
相应的生成上述网格的函数为
1 /*调用ANSYS生成网格*/
2 void GenMesh(struct circle0 *circle,double xmax,double xmin,double ymax,double ymin)
3 {
4
5 FILE* ansysfile;
6 ansysfile=fopen("temp.inp","w");
7
8 fprintf(ansysfile,"/CLEAR,NOSTART\n");
9 fprintf(ansysfile,"/PREP7\n");
10
11 //创建基体
12 fprintf(ansysfile,"RECTNG,%f,%f,%f,%f,\n",xmin,xmax,ymin,ymax);
13
14 //创建颗粒
15 for (int i=0; i<N0; i++)
16 {
17 if (circle[i].r<0.00000001)
18 {
19 continue;
20 }
21 fprintf(ansysfile,"K,%d,%f,%f\n",i+1000,
22 circle[i].centre_x,circle[i].centre_y);
23 fprintf(ansysfile,"CIRCLE,%d,%f\n",i+1000,circle[i].r);
24 }
25 fprintf(ansysfile,"ASBL,1,ALL\n");
26 fprintf(ansysfile,"AGLUE,All\n");
27
28 //创建单元类型
29 fprintf(ansysfile,"ET,1,PLANE42\n");
30 //创建材料
31 fprintf(ansysfile,"MP,DENS,1,0.1\n");//材料一
32 fprintf(ansysfile,"MP,DENS,2,10\n");//材料二
33 //为基体赋一号材料值
34 fprintf(ansysfile,"ASEL,S,AREA,,ARINQR(0,14)\n");
35 fprintf(ansysfile,"AATT,1\n");
36
37 //为颗粒赋二号材料值
38 fprintf(ansysfile,"ASEL,ALL\n");
39 fprintf(ansysfile,"ASEL,U,AREA,,ARINQR(0,14)\n");
40 fprintf(ansysfile,"AATT,2\n");
41 fprintf(ansysfile,"ALLSEL,ALL\n");
42
43 //生成网格
44 //网格尺寸
45 fprintf(ansysfile,"SMRT,3\n");
46 //fprintf(ansysfile,"AESIZE,ALL,0.6,\n");
47
48 fprintf(ansysfile,"MSHAPE,1,2D\n");
49 fprintf(ansysfile,"MSHKEY,0\n");
50 //开始划分网格
51 fprintf(ansysfile,"AMESH,ALL\n");
52
53 //输出网格
54 // fprintf(ansysfile,"/INPUT,ExportMesh,inp\n");
55
56 fclose(ansysfile);
57 //---------------调用ANSYS生成网格
58 system("D:\\Progra~1\\AnsysI~1\\v100\\ANSYS\\bin\\intel\\ansys100 -b -i temp.inp -o mesh.log");
59 // system("tec360 tecplot.plt");
60 }
3.2 复合材料中的椭圆形颗粒网格剖分
对于椭圆形颗粒,网格剖分结果如图所示
生成上述网格的代码为
1
2 //生成颗粒是椭圆时的网格
3 //椭圆的数据结构是下面的struct
4 //Genmesh函数是调用ANSYS生成网格的函数
5 void GenMesh(struct elliptic0 *ellipse,int ellipsenum,double xmax,double xmin,double ymax,double ymin)
6 {
7
8 FILE* ansysfile;
9 ansysfile=fopen("temp.inp","w");
10 fprintf(ansysfile,"FINISH\n");
11 fprintf(ansysfile,"/CLEAR,NOSTART\n");
12 fprintf(ansysfile,"/PREP7\n");
13
14 //创建基体
15 fprintf(ansysfile,"RECTNG,%f,%f,%f,%f,\n",xmin,xmax,ymin,ymax);
16
17 //创建椭圆颗粒
18 for (int i=0; i<ellipsenum; i++)
19 {
20 //把工作平面移到椭圆圆心。
21 fprintf(ansysfile,"wpave,%f,%f\n",ellipse[i].center_x,ellipse[i].center_y);
22 //旋转工作平面
23 fprintf(ansysfile,"wprota,%f\n",180.0/pi*ellipse[i].theta);
24 //激活当前活动坐标系为工作平面
25 fprintf(ansysfile,"csys,4\n");
26 //生成上半个椭圆
27 fprintf(ansysfile,"k,%d,%f\n",2*i+1000,ellipse[i].a);
28 fprintf(ansysfile,"k,%d,%f\n",2*i+1001,-ellipse[i].a);
29 fprintf(ansysfile,"cswpla,%d,2,%f\n",2000+2*i,(ellipse[i].b/ellipse[i].a));
30 fprintf(ansysfile,"l,%d,%d\n",2*i+1000,2*i+1001);
31 fprintf(ansysfile,"reflect_line=_return\n");//记录线的编号
32 fprintf(ansysfile,"csys,4\n"); //切换活动坐标系为工作平面
33 fprintf(ansysfile,"lsymm,y,reflect_line\n");//反射生成椭圆
34 //调整工作平面使其与整体坐标系重合
35 fprintf(ansysfile,"wprota,%f\n",-180.0/pi*ellipse[i].theta);
36 fprintf(ansysfile,"CSYS,0\n");
37 fprintf(ansysfile,"WPAVE,0,0,0\n");
38 }
39
40 fprintf(ansysfile,"ASBL,1,ALL\n");
41 fprintf(ansysfile,"AGLUE,All\n");
42
43 //创建单元类型
44 fprintf(ansysfile,"ET,1,PLANE42\n");
45 //创建材料
46 fprintf(ansysfile,"MP,DENS,1,0.1\n");//材料一
47 fprintf(ansysfile,"MP,DENS,2,10\n");//材料二
48 //为基体赋一号材料值
49 fprintf(ansysfile,"ASEL,S,AREA,,ARINQR(0,14)\n");
50 fprintf(ansysfile,"AATT,1\n");
51 //为颗粒赋二号材料值
52 fprintf(ansysfile,"ASEL,ALL\n");
53 fprintf(ansysfile,"ASEL,U,AREA,,ARINQR(0,14)\n");
54 fprintf(ansysfile,"AATT,2\n");
55 fprintf(ansysfile,"ALLSEL,ALL\n");
56 //生成网格
57 //网格尺寸
58 fprintf(ansysfile,"SMRT,3\n");
59 //fprintf(ansysfile,"AESIZE,ALL,0.6,\n");
60
61 fprintf(ansysfile,"MSHAPE,1,2D\n");
62 fprintf(ansysfile,"MSHKEY,0\n");
63 //划分网格
64 fprintf(ansysfile,"AMESH,ALL\n");
65
66 //输出网格
67 fprintf(ansysfile,"/INPUT,ExportMesh,inp\n");//输出网格到Tecplot
68 fprintf(ansysfile,"EWRITE,elemlist1,txt\n");//输出单元
69 fprintf(ansysfile,"NWRITE,nodelist1,txt\n");//输出节点
70
71 fclose(ansysfile);
72 //---------------调用ANSYS生成网格
73 system("D:\\Progra~1\\AnsysI~1\\v100\\ANSYS\\bin\\intel\\ansys100 -b -i temp.inp -o mesh.log");
74 system("tec360 tecplot.plt");
75 }
4.上述方法的缺点
--------先想想在写。
【转载】ANSYS的APDL与C语言混合编程(实例)的更多相关文章
- Java语言与C语言混合编程(2)--在Java中调用C语言本地库
在上一篇文章中介绍了Java语言中的native关键字,以及Java语言调用C语言的编译生成本地动态链接库(DLL)实现加法运算的小例子,本文通过一个更加详细的例子,深入讲解Java语言调用C语言的函 ...
- Java语言与C语言混合编程(1)--Java native 关键字
一. 什么是 native Method 简单地讲,一个 native Method 就是一个java调用非java代码的接口.一个 native Method 是这样一个java的方法:该方法的实现 ...
- Android程序中,内嵌ELF可执行文件-- Android开发C语言混合编程总结
前言 都知道的,Android基于Linux系统,然后覆盖了一层由Java虚拟机为核心的壳系统.跟一般常见的Linux+Java系统不同的,是其中有对硬件驱动进行支持,以避开GPL开源协议限制的HAL ...
- linux下C语言多线程编程实例
用一个实例.来学习linux下C语言多线程编程实例. 代码目的:通过创建两个线程来实现对一个数的递加.代码: //包含的头文件 #include <pthread.h> #include ...
- Swift 与 C 语言混合编程
前言 作为一种可与 Objective-C 相互调用的语言,Swift 也具有一些与 C 语言的类型和特性,如果你的代码有需要,Swift 也提供了和常见的 C 代码结构混合编程的编程方式. 1.基本 ...
- Dart 调用C语言混合编程
Dart 调用C语言本篇博客研究Dart语言如何调用C语言代码混合编程,最后我们实现一个简单示例,在C语言中编写简单加解密函数,使用dart调用并传入字符串,返回加密结果,调用解密函数,恢复字符串内容 ...
- SQL+C#:一次多语言混合编程的经验总结
1.用JAVA做,采取轮询策略: 2.用sql语言+C#混合编程,采取触发策略
- 5种语言混合编程:C++、JS、python、Lisp、汇编
/* 混合C++.JS.python.Lisp.汇编 1种语言,5种语法 */ main { //C++ vector<int> v; v.push(2); putsl(v.size()) ...
- Swift语言与Objective-C语言混合编程
首先创建一个Swift的Single View工程 然后直接在工程中新建OC文件: 然后选择OC语言之后会问你是否自动创建OC和Swift的中间文件: 然后工程文件夹里就有了三个文件: 现在OC头文件 ...
随机推荐
- Linux_命令_积累
1.ps 查看进程状态 ZC: "ps -a" 和 "ps a" 有区别... (具体查看 "man ps") 1.1.ps aux 1.2 ...
- Codeforces Round #263 (Div. 1)
B 树形dp 组合的思想. Z队长的思路. dp[i][1]表示以i为跟结点的子树向上贡献1个的方案,dp[i][0]表示以i为跟结点的子树向上贡献0个的方案. 如果当前为叶子节点,dp[i][0] ...
- 【转】linux命令详解:md5sum命令
[转]linux命令详解:md5sum命令 转自:http://blog.itpub.net/29320885/viewspace-1710218/ 前言 在网络传输.设备之间转存.复制大文件等时,可 ...
- Shell重定向文件描述符
#!/bin/bash 最近在看shell,各种困惑,不过解决困惑的感觉还是很不错的.废话少说,linux中使用文件描述符来标识每个文件对象.文件描述符为一个非负整数,可以唯一标识会话中打开 ...
- 《大道至简》之第一章:编程的精义读后感(JAVA伪代码)
——大道至简之编程的精义读后感(JAVA伪代码) import.java.大道至简.*; import.java.愚公移山.*; public class YuGongYiShan{ 愚公 = {项目 ...
- iOS7中的ViewController切换
转自:https://onevcat.com/2013/10/vc-transition-in-ios7/ iOS 7 SDK之前的VC切换解决方案 在深入iOS 7的VC切换效果的新API实现之前, ...
- svn用法:摘自http://wh-vip-126-com200907043220.iteye.com/blog/424642
下载最新版本svn(http://www.visualsvn.com/server/ ) 1.安装SVN 2.创建user与repository 3.Eclipse下配置SVN插件 3.1帮助- ...
- 介绍开源的.net通信框架NetworkComms框架之九 合并DLL
原文网址: http://www.cnblogs.com/csdev Networkcomms 是一款C# 语言编写的TCP/UDP通信框架 作者是英国人 以前是收费的 目前作者已经开源 许可是 ...
- iOS开发拓展篇—CoreLocation定位服务
iOS开发拓展篇—CoreLocation定位服务 一.简单说明 1.CLLocationManager CLLocationManager的常用操作和属性 开始用户定位- (void)startUp ...
- Python基础教程-第一章-变量、函数、字符串
1.1变量 变量基本上就是代表(或者引用)某个值的名字,举例来说,如果希望用x代表3,只需要执行下面的语句即可: >>>x = 3 这样的操作称为赋值(assignment),值3赋 ...