IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> Python知识库 -> 基于python的掺杂介质六面体nastran网格生成脚本 -> 正文阅读

[Python知识库]基于python的掺杂介质六面体nastran网格生成脚本

1.效果图

1.1 网格

在这里插入图片描述在这里插入图片描述
M e s h 1 Mesh_1 Mesh1? M e s h 2 Mesh_2 Mesh2?

1.2 应力

(使用COMSOL求解)

在这里插入图片描述在这里插入图片描述
m i s s e s misses misses应力切应力

3.程序

import numpy as np
import scipy.interpolate as si
import matplotlib.pyplot as plt

N_origin= 10 # 初始矩阵尺寸
N_inter= 50 # 插值之后矩阵尺寸
t_0= np.linspace(0,1,N_origin)
t_1= np.linspace(0,1,N_inter)
originMat= np.random.random(( N_origin, N_origin, N_origin))
interMat= np.zeros(( N_inter, N_inter, N_inter))
# mesh_01, mesh_02, mesh_03 = np.meshgrid( t_0, t_0, t_0)

# originMat[originMat> 0.3]= 1.0
# originMat[originMat<= 0.3]= 0.0
# plt.imshow(originMat[0,:,:], cmap= "gray")
# plt.colorbar()
# plt.show()

mesh_1 = np.meshgrid( t_1, t_1, t_1)
xi_list= np.zeros((N_inter*N_inter*N_inter, 3))
# [[x1, y1, z1], [x2, y2, z2],...,[xn, yn, zn]]
for i in range(N_inter):
    for j in range(N_inter):
        for k in range(N_inter):
            for m in range(3):
                xi_list[((i*N_inter)+j)*N_inter+ k, m]= mesh_1[m][i, j, k]
interList= si.interpn(points= ( t_0, t_0, t_0),values= originMat,xi= xi_list, method= "linear")
for i in range(N_inter):
    for j in range(N_inter):
        for k in range(N_inter):
            interMat[i, j, k]= interList[((i*N_inter)+j)*N_inter+ k]

# 二值化
interMat[interMat> 0.4]= 1.0
interMat[interMat<= 0.4]= 0.0

# plt.imshow(interMat[0,:,:], cmap= "gray")
# plt.colorbar()
# plt.show()
# interMat= "[读取到的矩阵]"
# N_inter = "[读取到的矩阵尺寸]"


# 开始写入网格文件
file= open("mesh.nas", mode= "w")
file.write("$ Generated by COMSOL 5.6.0.401\n$ Jul 15 2021, 15:23\nBEGIN BULK\n")
file.write("$ Grid data section\n")
for i in range(N_inter +1):
    for j in range(N_inter +1):
        for k in range(N_inter +1):
            file.write("GRID{:12d}        ".format(((i*(N_inter+1))+j)*(N_inter+1)+ k+1))
            for xi in [k/N_inter*1.0, j/N_inter*1.0, i/N_inter*1.0]:
                if abs(xi- int(xi))< 1e-6:
                    file.write("{:8.5f}".format(xi))
                else:
                    file.write("{:8.6f}".format(xi))
            file.write("\n")
file.write("$ Element data section\n")
"""
CHEXA          1       1       1       2       7       3       4       6+CONT
+CONT         15       9
"""
for i in range(N_inter):
    for j in range(N_inter):
        for k in range(N_inter):
            num_ele= (i*N_inter+j)*N_inter+ k+1
            node_1= num_ele+ (i*N_inter+ j)+(i*(N_inter+1))
            node_2= node_1+ 1
            node_3= node_1+ N_inter+1+1
            node_4= node_3- 1
            node_5= node_1 + (N_inter+1)*(N_inter+1)
            node_6= node_5 + 1
            node_7= node_5+ N_inter+1+1
            node_8= node_7- 1
            file.write("CHEXA{:11d}{:8.0f}{:8d}{:8d}{:8d}{:8d}{:8d}{:8d}+CONT\n+CONT{:11d}{:8d}\n".format(
                num_ele,interMat[i,j,k]+1,node_1, node_2, node_3, node_4, node_5, node_6, node_7, node_8 ))
file.write("$ Property data section\n")
file.write("PSOLID         1\n")
file.write("PSOLID         2\n")
file.write("ENDDATA\n")
file.close()
  Python知识库 最新文章
Python中String模块
【Python】 14-CVS文件操作
python的panda库读写文件
使用Nordic的nrf52840实现蓝牙DFU过程
【Python学习记录】numpy数组用法整理
Python学习笔记
python字符串和列表
python如何从txt文件中解析出有效的数据
Python编程从入门到实践自学/3.1-3.2
python变量
上一篇文章      下一篇文章      查看所有文章
加:2021-08-02 10:46:25  更:2021-08-02 10:48:29 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年12日历 -2024/12/25 14:57:15-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码
数据统计