| 查看: 5614 | 回复: 3 | ||
[求助]
comsol中怎么生成三维随机粗糙表面 已有1人参与
|
| 粗糙表面类似于正态分布 |
回复此楼» 猜你喜欢
为什么到现在没收到通知？
已经有3人回复
-
有没有仍没收到信息的
已经有5人回复
-
面上意见出来了
已经有10人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
【答案】应助回帖
感谢参与，应助指数 +1
|
您好： 我这里有MATLAB生成的粗糙裂纹程序，你可以导入COMSOL试试。还有COMSOL种随机函数，好像可以构建。希望有所帮助。 clear clc N=512; L=0.05; delta=0.05;%表面均方根粗糙度 betax=10;%x方向的相关长度 betay=10; C=1;%功率谱密度 dx=L/N;dy=dx; NN=-N/2:N/2-1; [Nx,Ny]=meshgrid(NN,NN); taux=dx.*Nx;tauy=dy.*Ny; %%生成具有指定自相关函数的粗糙表面 eta=randn(N,N);%高斯分布白噪声 A=fft2(eta);%傅里叶变换 R=zeros(N,N); R=delta^2*exp(-2.3*((taux/betax).^2+(tauy/betay).^2).^0.5);%自相关函数 Gz=1/(2*pi^2).*fft2(R);%功率谱密度函数 H=(Gz/C).^0.5;%传递函数 Z=H.*A;%表面高度的傅里叶变换 z=ifft2(Z);%表面高度分布 figure(1) meshz(abs(z)); title('rough surface'); axis square |
2楼2016-04-25 17:14:59
3楼2016-05-11 11:05:00
leslie1996
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2.5
- 散金: 5
- 帖子: 19
- 在线: 12.9小时
- 虫号: 7334571
- 注册: 2017-10-20
- 性别: MM
- 专业: 金属结构材料
4楼2018-10-31 20:34:01
30