close all
clear all
clc
%% 离散曲率值矩阵k,弧长矩阵ss,(包含边界点k=0,s=0)求出对应弧段的弧长-曲率线性方程
%Strain=[0 1.75044E-04 1.53376E-04 1.32575E-04 1.12096E-04 9.16930E-05 7.12955E-05 5.08464E-05 3.01552E-05]%应变
%Dspl=[0 -8.57915E-04 -2.36228E-03 -4.51358E-03 -7.22375E-03 -1.04064E-02 -1.39754E-02 -1.78449E-02 -2.19286E-02]%位移
Data= load('F:\ANSYS working directory\modal2\simulation.txt'); %从数据采集文件中提取
Datas=Data(:,1);
Strain=[zeros(1,4);reshape(Datas,8,4)];
s(1)=0;i=2:9;s(i)=0.15+0.1*(i-2); %弧长:单位m
ss=repmat(s',1,4);
h=0.005;%厚度
k=2*Strain/h;%曲率
for i=1:1:8
for j=1:4
a(i,j)=(k(i+1,j)-k(i,j))/(ss(i+1,j)-ss(i,j)); %线性方程的参数
b(i,j)=(k(i,j)*ss(i+1,j)-k(i+1,j)*ss(i,j))/(ss(i+1,j)-ss(i,j));
end
end