CCFF薄板弯曲
3视图(30天)
显示旧的评论
你好,我有这段代码,测量误差和模式下薄板的形状点力边界条件的ccs(夹在两边,只是suppourted另两个)。我需要改变它的边界条件CCFF(夹在两边和自由在另外两个)。我也需要改变tho点区域力量。有人可以请帮助。
clc
清晰的所有
关闭所有
% %一般参数:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
%频率范围:
df0 = 1;
f0 = 1: df0:2000;
%的材料参数:
E0 = 5.5 e8;
rho0 = 7800;
h0 = 0.003;
nu0 = 0.33;
eta0 = 0.01;
m: %维度
L0 = [0.8, 0.4];
%的模式:
nmodes = 16;
%点力量和响应位置m:
F0 = 1;
xf0 = [0.6, 0.3];
xr0 = [0.6, 0.3];
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %求解特征方程:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
fun0 = @(β)cos(β)* cosh(β)1;
fund0 = @(β)sin(β)* cosh(β)+ cos(β)* sinh(β);
fun1 = @(β)beta-fun0(β)。/ fund0(β);
beta0 =((1:4) + 1/2) *π;
beta1 = 0(大小(beta0));
为j = 1:4
beta00 = beta0 (j);
为k = 1:1000
beta1 (j) = fun1 (beta00);
如果(abs (fun0 (beta1 (j))) < 1平台以及)
打破
其他的
beta00 = beta1 (j);
结束
结束
结束
如果(nmodes > 4)
beta1 = [beta1, ((5: nmodes) + 1/2) *π);
结束
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %振动模式的形状:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
om0 = 2 *π* f0;
在x %固支边界条件:
k1 = beta1 / 10 (1);
γ= @ (m) (cos (k1 (m) * 10 (1) cosh (k1 (m) * 10(1)))。/(罪(k1 (m) * 10 (1) sinh (k1 (m) * 10 (1)));
模型= @ (x, m) (cos (k1 (m) * x)γ(m) * sin (k1 (m) * x)) cosh (k1 (m) * x) * (1-gamma1 (m) *双曲正切(k1 (m) * x));
mode2 = @ (x, m) (cos (k1 (m) * x) sin (k1 (m) * x)) +(罪(k1 (m) * 10 (1) cos (k1 (m) * 10 (1)))。* exp (k1 (m) * (10 (1) - x)) exp (k1 (m) * x);
modedd1 = @ (x, m) k1 (m) ^ 2 * (- (cos (k1 (m) * x)γ(m) * sin (k1 (m) * x)) cosh (k1 (m) * x) * (1-gamma1 (m) *双曲正切(k1 (m) * x)));
modedd2 = @ (x, m) k1 (m) ^ 2 * (- (cos (k1 (m) * x) sin (k1 (m) * x)) +(罪(k1 (m) * 10 (1) cos (k1 (m) * 10 (1)))。* exp (k1 (m) * (10 (1) - x)) exp (k1 (m) * x));
%在y方向简支边界万博1manbetx条件:
k2 = (1: nmodes) *π/ 10 (2);
mode3 = @ (y, n)罪(k2 (n) * y);
modedd3 = @ (y, n) k2 (n) ^ 2 *罪(k2 (n) * y);
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %图形输出:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
x1 = linspace (0, 10 (1), 1000);
日元= linspace (0, 10 (2), 1000);
%夹在方向上振动模式:
图(1)
次要情节(1、2、1)
为k = 1:4
如果(k < = 5)
情节(x1,模型(x1, k),“b -”,“线宽”,2)
持有在
其他的
情节(x1, -mode2 (x1, k),“b -”,“线宽”,2)
持有在
结束
结束
持有从
ylim ([2])
网格(“上”)
%简支在y方向上振万博1manbetx动模式:
次要情节(1、2、2)
为k = 1:4
情节(mode3日元(y1, k),“b -”,“线宽”,2)
持有在
结束
持有从
ylim ([2])
网格(“上”)
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %刚度和质量:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
D0 = E0 * (1 + 1 j * eta0) * h0 ^ 3 / (12 * (1-nu0 ^ 2));
m0 = rho0 * h0;
%初始化所有可能的模式:
[钱数,j₀] = meshgrid (1: nmodes, 1: nmodes);
I0 = I0 (:);
j₀= j₀(:);
%质量矩阵
M0 = M0 *(刺激(10)/ 2)* speye (nmodes ^ 2);
%刚度矩阵
K01 =诊断接头(D0 * (k1(钱数)^ 4 + k2 (j₀)。^ 4) *(刺激(10)/ 2));
K02x = 0 (nmodes ^ 2, 1);
K02y = 0 (nmodes ^ 2, 1);
x0 = linspace (0, 10 (1), 1000);
y0 = linspace (0, 10 (2), 1000);
%固支边界条件:
k01 = @ (x, m, p) modedd1 (x,米)。*模型(x, p);% m < 5 & p < 5
k02 = @ (x, m, p) modedd1 (x,米)。* mode2 (x, p);% m < 5 & p > = 5
k03 = @ (x, m, p) modedd2 (x,米)。*模型(x, p);% m > = 5 & p < 5
k04 = @ (x, m, p) modedd2 (x,米)。* mode2 (x, p);% m > = 5 & p > = 5
%简支边界条件:万博1manbetx
k05 = @ (y, n, q modedd3 (y, n)。* mode3 (y, q);
在x %固支边界条件:
count0 = 1;
为k = 1: nmodes
为l = 1: nmodes
如果(k < 5 & & l < 5)
K02x (count0) = trapz (x0, k01 (x0, k, l));
count0 = count0 + 1;
elseif(k < 5 & & l > = 5)
K02x (count0) = trapz (x0, k02 (x0, k, l));
count0 = count0 + 1;
elseif(k > = 5 & & l < 5)
K02x (count0) = trapz (x0, k03 (x0, k, l));
count0 = count0 + 1;
elseif(k > = 5 & & l > = 5)
K02x (count0) = trapz (x0, k04 (x0, k, l));
count0 = count0 + 1;
结束
结束
结束
%在y方向简支边界万博1manbetx条件:
count0 = 1;
为k = 1: nmodes
为l = 1: nmodes
K02y (count0) = trapz (y0, k05 (y0, k, l));
count0 = count0 + 1;
结束
结束
% x和y方向上的融合:
K02x =重塑(K02x nmodes nmodes);
K02y =重塑(K02y nmodes nmodes);
K02x =重塑(K02x nmodes nmodes);
K02x = kron (K02x (nmodes));
K02y =重塑(K02y nmodes nmodes);
K02y = repmat (K02y nmodes nmodes);
%总刚度矩阵:
K00 =稀疏(K01 + 2 * D0 *稀疏(K02x * K02y));
K0 = (K00 + K00。') / 2;
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %点力和反应:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
load0 = 0 (nmodes ^ 2, 1);
resp0 = 0 (nmodes ^ 2, 1);
count0 = 1;
为k = 1:元素个数(钱数)
如果(I0 (k) < 5)
load0 (count0) = F0 *模型(xf0 (1), I0 (k))。* mode3 (xf0 (2), j₀(k));
resp0 (count0) =模型(xr0 (1), I0 (k))。* mode3 (xr0 (2), j₀(k));
count0 = count0 + 1;
其他的
load0 (count0) = F0 * mode2 (xf0 (1), I0 (k))。* mode3 (xf0 (2), j₀(k));
resp0 (count0) = mode2 (xr0 (1), I0 (k))。* mode3 (xr0 (2), j₀(k));
count0 = count0 + 1;
结束
结束
%方程解系统
w0 = 0(大小(f0));
为k = 1:元素个数(f0)
Mat0 = K0-om0 (k) ^ 2 * M0;
C0 = Mat0 \ load0;
w0 (k) = 1 j * om0 (k) * (C0。* resp0)总和;
结束
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
% %图形输出:
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
图(1)
次要情节(2,1,1)
情节(f0, 10 * log10 ((abs (w0) / 5 e-8)。^ 2/2),“b -”,“线宽”,2)
网格(“上”)
次要情节(2,1,2)
情节(f0 (180 / pi) *角(w0),“b -”,“线宽”,2)
网格(“上”)
% = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
