TPMS曲面的matlab生成及可视化-导出stl文件-导出abaqus的inp文件-周期性边界条件施加

自编写函数包含的功能:

tpms曲面生成

tpms可视化

tpms曲面采用拉普拉斯方法进行平滑(为了可以在abaqus中通过单元质量检测),函数为SmoothIsosurface

tpms导出stl文件,函数为WriteSTL

tpms导出inp文件(施加周期性边界条件)函数为WriteINP

matlab 复制代码
clc
clear
Numberofgridpoints=[20,20,20];
Numberofunitcells=[1,1,1];
[x, y, z]=meshgrid(linspace(0, Numberofunitcells(1),Numberofgridpoints(1)), ...
                linspace(0, Numberofunitcells(2),Numberofgridpoints(2)), ...
                linspace(0, Numberofunitcells(3),Numberofgridpoints(3)));

cx = cos(2*pi*x); cy = cos(2*pi*y); cz = cos(2*pi*z); 
sx = sin(2*pi*x); sy = sin(2*pi*y); sz = sin(2*pi*z);
c2x = cos(2*2*pi*x); c2y = cos(2*2*pi*y); c2z = cos(2*2*pi*z); 
s2x = sin(2*2*pi*x); s2y = sin(2*2*pi*y); s2z = sin(2*2*pi*z);
c3x = cos(3*2*pi*x); c3y = cos(3*2*pi*y); c3z = cos(3*2*pi*z); 
s3x = sin(3*2*pi*x); s3y = sin(3*2*pi*y); s3z = sin(3*2*pi*z);

f= (sx.*sy).*sz+(sx.*cy).*cz+(cx.*sy).*cz+(cx.*cy).*sz;
%f = cx+cy+cz;

FV=isosurface(x,y,z,f,0); 
%FV3=reducepatch(FV3,30^3);%缩减FV3的补片的面数,同时尝试保留原始的整体形状
% FV=reducepatch(FV,0.9);

newFV=ClearIsosurface(FV,1e-3);
newFV = SmoothIsosurface(newFV,1e-3);
WriteINP(newFV, 'main.inp')
相关推荐
Python×CATIA工业智造1 小时前
Frida RPC高级应用:动态模拟执行Android so文件实战指南
开发语言·python·pycharm
我叫小白菜2 小时前
【Java_EE】单例模式、阻塞队列、线程池、定时器
java·开发语言
狐凄2 小时前
Python实例题:基于 Python 的简单聊天机器人
开发语言·python
weixin_446122463 小时前
JAVA内存区域划分
java·开发语言·redis
悦悦子a啊3 小时前
Python之--基本知识
开发语言·前端·python
QuantumStack4 小时前
【C++ 真题】P1104 生日
开发语言·c++·算法
whoarethenext4 小时前
使用 C++/OpenCV 和 MFCC 构建双重认证智能门禁系统
开发语言·c++·opencv·mfcc
代码的奴隶(艾伦·耶格尔)5 小时前
后端快捷代码
java·开发语言
hie988946 小时前
MATLAB锂离子电池伪二维(P2D)模型实现
人工智能·算法·matlab
Jay_5156 小时前
C++多态与虚函数详解:从入门到精通
开发语言·c++