Another Image Won the FWW

本文作者spanzhang(张友邦)介绍了一幅名为“Burn”的分形艺术作品,该作品被选入本周的FWW。通过使用FMF软件的最新版本及其新特性“ImageBasedTrap”,作者分享了创作复杂且美观的分形图案的经验。

本文属spanzhang(张友邦)原创,发布地址为:http://blog.youkuaiyun.com/spanzhang。转载或引用请注明原文之出处,谢谢!

It really means an honor to me. One of my images, named Burn, was selected into the FWW this week. Here's the thumbnail of it.

I cooked this image with my newest version of FMF. It has a coloring component named 'Image Based Trap' which allows the user to select an outside image as the trap basis and another optional image as the weight coefficient provider. But as all can see, to produce a well-looking fractal is not an easy thing, especially to freshman.

At the forum, I found this reply to a thread. A number of people here have worked with Ferryman Fractals and have enjoyed that software very much.  Spanzhang is the author of the program, and he is very good about answering questions.  I hope that helps.  You could try searching (ferryman) in the galleries to find people who have used it.

clear syms xi theta t real; load('ModeFUVWn_0_m=2.mat'); FY = 1; % 周向模态函数 NX=length(ModeFU); NY=length(FY); % % 参数赋值 % L = 0.1; % 圆柱壳长度 (m) % h = 5e-4; % 厚度 (m) % R = 0.05; % 半径 (m) % rho = 2600; % 密度 (kg/m³) % zeta_val = 0.05; % 阻尼比 % mu = 0.3; % 泊松比 % E = 7e10; % 弹性模量 (Pa) % Omega_bar = 0; % 无量纲转速 % L_r = L / R; % 无量纲长度 % H_r = h / R; % 无量纲厚度 %参数2 %圆柱壳参数设置2 L = 520e-3; % 长度,单位:m (从520 mm转换) R = 149.4e-3; % 半径,单位:m (从149.4 mm转换) h = 0.519e-3; % 厚度,单位:m (从0.519 mm转换) E = 1.95e11; % 弹性模量,单位:Pa rho = 7800; % 密度,单位:kg/m³ mu = 0.3; % 泊松比 (代码中使用mu表示) Omega1=0; N_theta_0_1=0; %矩阵无量纲化需在能量方程统一除以rho * h * R * L * h/(gamma^2) % 计算 Q11, Q22, Q12, G Q11 = E / (1 - mu^2); Q22 = E / (1 - mu^2); Q12 = (mu * E) / (1 - mu^2); G = E / (2 * (1 + mu)); gamma = R*sqrt(rho*(1-mu^2)/E); % 计算积分矩阵 H_i [Fuu00x,Fuu11x] = Integraluux(ModeFU); [Fww00x,Fww11x,Fww22x,Fww02x] = Integralwwx(ModeFW); [Fuw10x] = Integraluwx(ModeFU,ModeFW); [F00y,F01y,F02y,F11y,F22y] = Integraly(FY,0,2*pi); [H1, H2, H3, H6, H13, H14, H15, H16, H17, H18] = ... ComputeAllIntegralMatricest_n_0(Fuu00x, Fuu11x, Fuw10x,... Fww00x, Fww22x, Fww02x, Fww11x, F00y, F01y, F02y, F11y, F22y); % --- 构造无量纲矩阵 --- % (1) 质量矩阵 M_bar M_bar = (L* R* h* rho)* blkdiag(H1, H13); % (4) 刚度矩阵 K_bar K11 = (R* h* Q11)/ L* H2 + (N_theta_0_1 + h* G)/ R* L* H3; K13 = h* Q12* H6; K33 = h* Q22/ R* L* H13 + (R* h^3* Q11/12)/L^3 * H14 + ... (h^3* Q12/ (12*R))/L* H15 + (h^3*Q22/(12*R^3))*L* H16 + ... (h^3*G/(3*R))/L * H17 + N_theta_0_1/ R*L* H18; K_bar = [K11, K13; K13',K33]; % 合并边界条件 K_global = K_bar; A=M_bar\K_global; [Te3,Omega3]=eigs(A); Omega=diag(Omega3); f3=sort(real(sqrt(Omega)./(2*pi))); 2*2矩阵,求解出来2个频率,怎么区分哪个是u对应的,哪个是w对应的频率
最新发布
09-07
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值