clear all
syms L H E P
x=[L,H,E,P];
mu_x=[9,0.2,2e5,1000];
sigma_x=[0.9,0.02,2e4,100];
g=beam1(x);
g1=0.1125-g
[m, n]=size(x);
for i=1:m
fidout=fopen('modeling.txt','w'); %建立并打开一个txt文件用以写入ANSYS命令流
if fidout==-1
error('Can not rewrite the APDL file');
else
% fprintf(fidout,'L=%6.4f\n',x(i,1));
% fprintf(fidout,'H=%8.6f\n',x(i,2));
% fprintf(fidout,'E=%13.1f\n',x(i,3));
% fprintf(fidout,'P=%8.6f\n',x(i,4));
%ANSYS命令流
fprintf(fidout,'/PREP7 \n'); %进入前处理器
fprintf(fidout,'ET,1,BEAM188 \n'); %选取单元类型
fprintf(fidout,'MPTEMP,,,,,,,, \n'); %定义材料属性
fprintf(fidout,'MPTEMP,1,0 \n');