Another Image Won the FWW

作者spanzhang使用最新版FMF软件创作了一幅名为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、付费专栏及课程。

余额充值