matlab newmark车路耦合代码_基于变量耦合与软件集成的CFRP导流罩一体化设计方法...

笔者近期忙于毕业应聘与毕业论文，9.21-10.15的24天内未正常更新公众号，故补充更新3篇文章，保持一周更新一篇的状态。

基于变量耦合与软件集成的CFRP导流罩一体化设计方法，涉及一体化设计领域。现在的设计，往往考虑零部件在单个学科领域下如何优化设计使得所需功能最佳，近几年流行的协同优化设计方法开始考虑零部件宏细观的协同优化。但是零部件设计往往涉及多个学科。例如：机械结构优化设计领域、材料内部细观结构优化设计领域、空气动力学领域以及流体力学领域等等，考虑一体化设计时，往往能够得到的实际功能效果最理想状态。

本文提出一种一体化设计的新方法，并且以导流罩为案例做进一步说明。首先对优化设计变量分为公有变量与私有变量，对学科私有变量与公有变量进行耦合操作，之后基于软件脚本语言传递的功能进行一体化设计，并且对目标函数进行学科加权操作，实现一体化设计在多学科上的平衡。

本文以碳纤维复合材料(CFRP)风力发电机导流罩为例，为这一问题提供一种新的思路及解决方案。风力发电机导流罩的一体化设计与导流罩结构设计学科、材料学科中的CFRP铺层、流体动力学的流固耦合相关。

    该设计框架使用编程/脚本功能集成了一组设计软件。在SIEMENS NX中进行几何建模，在ICEM CFD中进行气动网格划分，在ANSYS FLUENT中进行流动求解，APDL中进行结构有限元建模与分析。采用NURBS思想编写MATLAB程序，导入APDL中进行CFRP铺层建模。通过对导流罩形状和尺寸的描述来评估反映导流罩性能的结构和气动指标。为了克服仿真模型的高成本，使高保真度优化问题得到有效的解决，采用了基于代理模型的优化策略。通过在低速空气导流设计中的应用，研究了该方法的可靠性。优化目标是在保持结构安全的前提下，使风力发电导流罩的减阻汇流效果最优，从而最大程度的实现风力发电。

1.基于变量耦合与软件集成的CFRP导流罩一体化设计方法。其特征在于所述设计方法由以下步骤实现：

步骤1：在软件集成操作前，首先对设计变量进行变量耦合操作，提高优化效率，例如：CFRP导流罩铺层厚度参数可以与公有变量中的结构厚度参数耦合。

步骤1说明：不同的学科私有变量之间不存在耦合，然而，从数学的角度来看，私有变量可能与学科中的公共变量相结合，提高了接下来遗传算法寻优操作的效率，耦合方法如下所示：

传统的协同优化模型：

(a)

(b)

图a为变量可以耦合             图b为变量不可以耦合

dX-Y≤a是判定公有变量与私有变量是否可以耦合的判定依据，a为一个给定的小正数。X与Y为设计变量的各自的范数空间，通过数学方法变量是否可以判定耦合，即两个变量X和Y的耦合是两个空间的距离(即X和Y)。

    变量耦合结束之后，进行一体化设计操作，一体化设计分为系统级与子系统级，子系统级也称为学科级。零件轻量化设计为系统级，结构优化、CFRP铺层优化、流固耦合等为学科级。系统级采用遗传算法进行寻优，使用线性加权方法将各学科级分配权重系数，整合成系统级的优化模型。学科子系统级采用序列二次规划优化算法，在导流罩流固耦合学科中引入代理模型近似。

步骤2，在SIEMENS NX中进行几何建模，建好的模型通过软件后台语言grip传递到CFD中。在NX提供的各种自动化工具中，图形编程交互语言(GRIP)灵活性较好。脚本开发的三个步骤，即编译，链接和执行，使用一组在批处理文件上操作的系统命令实现自动化。

步骤3，APDL的结构建模主要是零件内部的碳纤维复合材料铺层建模，通过初始铺层设计，编写遗传算法优化铺层角度；

gripbatch-c-dev:osfile -name:-dir: filespec

说明：output_file_name指定输出输出的文件名，directory_name是脚本所在的目录，filespec是脚本名。

步骤4，在ICEM CFD中进行气动网格划分，批处理模式执行通过以下系统命令实现 -batch -script 脚本是使用工具命令语言/工具包(Tcl/Tk)编写的脚本执行的，编写脚本自动完成网格生成过程。包括以下步骤(1)导入CAD几何模型

(2)创建流体域并定义边界表面

(3)划分计算域

(4)设置间距参数

(5)生成网格

(6)导出网格到msh文件

步骤5，在ANSYS FLUENT中进行流动求解，该模块负责导流罩表面的流场计算。采用软件ANSYS FLUENT 15.0作为流动求解器，通过以下系统命令流实现自动化 3d -hidden -t4 -i 。

步骤6，在APDL中进行结构有限元建模；把 MATLAB和ANSYS 联合应用于结构优化设计中，其中关键技术是 MATLAB自动调用运行于后台的 ANSYS。

在MATLAB的主程序中，通过SYSTEM 函数调用ANSYS批处理(BATCH)模式。ANSYS的批处理模式基于APDL语言，Matlab软件中编写自动调用Ansys APDL的语句，语句如下：

SYSTEM(’”D:\ProgramFiles\ANSYSinc\vl50\a-nsys\bin\intel\ANSYS150"-b-p ane3fl-i input file-o out file.bat’),其中："D:\ProgramFiles\ANSYSinc\vl50\ansys\bin\int-el\ANSYS150"为ANSYS软件安装的路径；

—b为ANSYS的批处理运行模式；

—p ane3fl为ANSYS多物理场分析模块Multiphysics的产品输入代码(ANSYS的help文件中有详述)；

—i input file为APDL语言宏文件所在的路径；

—o out file.bat为ANSYS运行结果输出文件，内容包括调用APDL命令流、ANSYS运行日志及错误报告等。

在不打开ANSYS情况下，MATLAB自动调用ANSYS, 使ANSYS后台运行，具体方法为：调用MATLAB中的system函数，格式中的执行程序与ANSYS 安装目录保持一致，每一子路径祈使用"\"隔开，ANSYS被MATLAB调用，在计算机后台工作，包括创建分析模型、应力分析及后处理等；

基于变量耦合与软件集成的CFRP导流罩一体化设计方法

其特征在于，步骤5具体包括：

创建* .txt文件，打开此文本文件并写入结构自身初始参数及与结构相关联的外部参数；具体方法为：采用读写模式创建一个*.txt文本文件，运用写入语句把结构外形尺寸参数、单元属性参数及结构外部作用载荷写入所创建的*.txt文本文件中，采用关闭文件语句关闭*.txt文本文件。具体程序为：

fp=fopen(‘*.txt’,’a’);%’*.txt’为文件名，’a’为打开方式；

fprintf(fp,’%d ’,x);%fp为文件句柄，制定要写入数据的文件，%d后有空格；

fclose(fp);%关闭文件。

Matlab打开并读取Ansys保存的**.txt文本文件，读取的格式如下：

fid=fopen('C:\ATOM.txt','r');

tline=fgetl(fid);

y=str2num(tline);

fclose(fid);

其中，fgetl指从已经打开的文件中读取一行，并且丢掉末尾的换行符；tline是把字符串转换成数值而命名的一个中间变量；str2num是把字符串转化成数值；

基于变量耦合与软件集成的CFRP导流罩一体化设计方法。其特征在于该方法可以同时考虑零部件宏观结构学科优化、微观材料铺层学科优化与流体动力学等多学科的协同，对零部件进行一体化设计。过程中首先进行变量耦合操作，再利用软件脚本语言传递参数实现一体化设计，快速求解最优的模型，以此达到一体化设计的要求。