学习笔记（3）——（改进）使用ANSYS调用MATLAB解决大量不同几何参数和荷载样本下连续求解结点位移

本篇文章是对“学习笔记（1）——使用MATLAB和ANSYS解决大量不同几何参数和荷载样本下连续求解结点位移”的改进”：
（1）将代码“fprintf(fid,[repmat(’%.4f\t’, 1, size(a,2)), ‘\n’], a’);”改为“fprintf(fid,[repmat(’%.4f ‘, 1, size(a,2)), ‘\n’], a’);”
（2）不需要使用UltraEdit对样本进行处理数据整理，直接将MATLAB生成数据文件导入ANSYS。

问题描述

本文的主要目的是已知模型的几何参数（弹性模量 E、泊松比 $\mu$、截面宽 b 和高 h 以及长度尺寸 L ）和荷载 P 这6个自变量 X_i 的分布类型，解决在大量的不同自变量样本下，如何利用ANSYS进行参数化建模求解相应的最大位移 Y ，最后将位移数据导出成 txt 文档。
以厚悬臂梁为例（如下图），进行样本处理和建模分析。

厚悬臂梁的几何参数和荷载的服从正态分布。

求：厚悬臂梁的顶端位移。

求解思路

1、根据自变量的分布类型，使用MATLAB生成100个样本数据；
2、在ANSYS中建立一个二维数组 beam ，将样本数据txt文件读入到ANSYS中；
3、在ANSYS中建立一个一维数组 weiyi ，采用循环语句进行建模，将求解的位移数据储存在数组中；
4、将一维数组 weiyi 中的100个位移数据导出成 txt 文件。

具体步骤

1、利用MATLAB生成100个样本数据

% 用matlab生成样本数据并导出成txt文件，数据表示不是科学计数法
clear
clc
format long g
% the samples
Mu = [ 2.5,200,0.225,1.0,3,3.5 ];                              % mean value
Std = Mu*0.1;                                                 % standard deviation
Sigma = diag(Std.^2);                                         % standard deviation;
N = 10000;
input = lhsnorm(Mu, Sigma, N);
a = input;
fid = fopen('beam.txt', 'wt');                                % -t模式按照文本而非二进制模式读写 
fprintf(fid,[repmat('%.4f ', 1, size(a,2)), '\n'], a');      % 转置并适时(写完一列后)添加换行
fclose(fid); 

利用ANSYS读取数据的过程解释
由于在读入时，采用的格式是(6f9.4)。其中：“6”：自变量个数，数据的列数；“f9.4”：一次读取9个字符，保留小数点后4位。
按照行将数据读写进在ANSYS中创建的100 x 6 的数组beam中，首先是将数据第一行读入到数组的第一行中，从右向左，一次读取9个字符，（“*”表示空格）。
第一次读入“***3.9041”，第二次读入“***2.7998”，第三次读入“***1.0221”，
第四次读入“***0.2573”，第五次读入“*187.4173”，第六次读入“***2.3509”。
第一行读取结束后，然后按照同样的规则依次读取下面每一行。

ANSYS的txt文件数据读入格式：

*Create, Dataread, mac                                     ! 创建一个宏文件，名为Dataread
*DIM, beam, array, 100, 6                                  ! 创建样本数据存储的数组（100 x 6）
*VREAD, beam(1,1), beam, TXT, , JIK, 6, 100                !  读入beam.txt文件中的数据
(6f9.4)                ! 读入格式，其中“6”表示自变量个数；“f9.4”表示每次从每一行从右往左读入9个字符，数据保留小数点后4位。
*End                                                       ! 结束宏文件的创建
Dataread                                                   ! 运行宏文件

ANSYS的txt文件数据写出格式：

*create, datawrite, mac                                     ! 创建一个宏文件，名为Datawrite
*cfopen, weiyi, txt                                         ! 创建一个写出文件，名为位移，格式为txt
*do, j, 1, 100                                              ! 利用*do命令逐个写出数组 weiyi 中的数据
  *vwrite, weiyi(j,1)
  (f9.6)                                 ! 写出格式，其中“f9.9”表示每次从每一行从右往左读入9个字符，数据保留小数点后9位。
*enddo                                                      ! 结束循环
*cfclose                                                    ! 关闭文件
*end                                                        ! 宏文件创建完毕
Datawrite                                                   ! 运行宏文件

3、使用ANSYS进行参数化编程建模
将txt文件读入ANSYS后，采用循环语句进行重复建模分析，并将每一次计算的位移储存在创建的 weiyi 数组中。当模型计算100次之后，将存入数组中的100个位移数据写入txt文件中。需要注意在建模过程中单位的统一。

finish
/clear

! 文件写入
*Create, dataread, mac
*DIM, beam, array, 100, 6                      ! 创建样本数据存储的数组
*VREAD, beam(1,1), beam, TXT, , JIK, 6, 100
(6f9.4)
*End
Dataread
*DIM, weiyi, array, 100, 1                     ! 创建位移数据存储的数组

/prep7
! 参数定义
*do, i, 1, 100

 !************************
  et, 1, beam189
  mp, ex, 1, beam(i,2)*1000
  mp, dens, 1, 10000e-12
  mp, prxy, 1, beam(i,3)

 !************************
  sectype, 1, beam, rect
  secdata, beam(i,4), beam(i,5)
 !************************
  k, 1, 0, 0, 0
  k, 2, beam(i,6), 0, 0
  k, 3, beam(i,6)/2, beam(i,6)/2, 0
  
  l, 1, 2
  latt, 1, , 1, , , 3, 1
  lesize, 1, , , 10
  lmesh, all
 !************************
  /solu
  dk, 1, all, 0
  fk, 2, fy, -beam(i,1)*1000
  allsel, all
  solve
  fdele,all
  
  /post1
  weiyi(i) = uy(2)
  
  /prep7
  lclear, all
  ldele, all
  kdele, 1, 2

*enddo

! 将位移数据导出
*create, datawrite, mac
*cfopen, weiyi, txt
*do, j, 1, 100
  *vwrite, weiyi(j,1)
  (f9.6)
*enddo
*cfclose
*end
Datawrite

结束语

下图显示了beam.txt文件中的数据与读入ANSYS中数组的数据，可以看出，读入的结果是正确的。

在ANSYS的运行路径中可找到 weiyi.txt ，即为我们所求的位移，共有100个位移数据。如下图所示：

可以看出，数组中的位移数据正确写入txt文件中了。
得到位移数据的txt文件之后，我们就可以使用 MATLAB 将位移数据导入，结合自变量样本点进行后续的分析，如 代理模型（Meta-model或Surrogate model） 的构建， 全局灵敏性分析（Global sensitivity analysis） 等。