蓝羽浅葱-CSDN博客

原创共轭传热拓扑优化算例分享

一个简单的共轭传热拓扑优化算例分享

2022-01-15 18:39:35 2865 4

文章目录MMA算法的推导及代码详解问题描述算法推导结果展示代码及注释参考文献MMA算法的推导及代码详解对于变密度的参数化方法，设计变量x为材料相对密度，在已知材料的物性，包括弹性模型、密度以及给定载荷的条件下，我们希望简支梁的柔度最小，或者说使得结构势能最小（结构在力作用下的位移，该力做的功也就是势能）。那么当材料总体积保持为常数不变，给定载荷不变的情况下柔度最小可以理解为结构刚度最强。问题描述在结构力学相关设计中，通常会出现结构冗余的情况，例如受到一定载荷作用下的悬臂梁和简支梁其结构质量可以进一步

2022-01-08 17:06:59 7015 15

原创标准oc算法的推导与99行代码详解

标准OC算法的数学推导与经典99行代码的详解

2022-01-07 17:41:19 2996 3

原创翻译翻译什么是拓扑优化算法

翻译翻译什么是拓扑优化算法广义地讲，优化问题就是在某种条件下寻找最优解。在问题与条件已经其关系明了的条件下，我们需要寻找某种能够在有限步内由已知条件收敛到未知变量的方法，这种方法被称为优化算法。对于优化算法来说最重要的要素有两个，一个是方向另一个是步长，怎样才能找到最优点多在的方向，用怎样的速度（步长）向最优点趋近是优化理论的主要研究对象。根据方向的选择方法优化算法可以被分成两类，第一种是基于梯度的优化算法，比如在某个连续的函数上求梯度，则最优点会在梯度或者梯度的负方向上。另一种就是无梯度算法，本质上说都

2022-01-07 16:46:39 2194

原创 COMSOL进行拓扑优化时遇到的问题及解决方案

一、变量未定义：CFLCMP通过全局搜索我们可以发现，在层流的方程视图中用CFLCMP来定义CFL数。这么看来CFL数是内置变量。查阅用户手册我们会找到以下描述。If the automatic expression for CFLloc is set to the built-in variable CFLCMP, then theautomatic setting suggests a PID regulator for the pseudo time step in the defaults

2021-07-22 16:58:13 6183

原创 HITRAN数据库的使用及考虑辐射时喷管烧蚀的仿真（附代码）

目录编译环境的搭建安装python系统变量的配置anaconda的安装与使用pycharm的安装与使用pycharm与anaconda的联动HITRAN数据库的使用HAPI第三方库的调用HAPI中函数的用法普朗克平均吸收系数的计算网格划分Fluent喷管烧蚀仿真辐射处理代码编译环境的搭建Python编译环境的搭建需要花费一定的精力，包括解释器的下载与配置，集成开发环境的选择，系统变量和用户变量的配置，除此之外还需要不断完善本地第三方库的配置，流行的做法是使用anaconda对众多第三方库进行管理。安装

2021-07-02 21:12:43 1963 1

原创喷管烧蚀仿真过程中的常见问题

1.计算过程中报错，浮点数溢出一般来说，浮点数溢出和边界条件的设定关系较大，可以查看boundary conditions2.稳态计算的结果与物理实际不符边界上没有问题，但是总体不符合物理常识，需要检查燃气性质。如果将密度设置为constant，则会出现如图所示的问题。因该将fluid定义为理想气体，这样在喷管中流动可压就可以得到符合物理常识的结果。3.计算结果总体符合常识但是细节处存在问题，或者计算过程中出现waning一般来说初场的设置会影响计算结果的准确性，例如出现图中所示的问题是由

2021-05-08 15:32:20 1255

原创关于在ICEM中导入CAD几何模型文件的问题

关于在ICEM中导入CAD几何模型文件的问题本例中所用到的bates-71b喷管模型已经上传到资源用传统方法导入几何模型（import gem）会出现很多问题，大多时候是warnning还是可以正常导入但是在导出网格导入FLUENT的时候会出现问题，比如FLUENT检查网格质量的时候会提示存在纵坐标为负数的网格，这有可能是一开始导入网格文件时出现的问题。最常见的报错是Sorry: Unable to guess the file content child process exited abnorm

2021-04-28 19:53:32 6208

原创 HITRAN数据库的使用及平均吸收系数的计算(二)

目前已经发展出多个python第三方库可以对HITRAN数据库进行处理，比较有名的是RADIS库。详情请移步https://github.com/radis/radis大家可以看RADIS开发者给我的回复。https://github.com/radis/radis/issues/249Q:I want to calculate Planck’s mean absorption coefficient using radis someday.A:This would be very easy to

2021-04-22 13:22:31 1541

原创 HITRAN数据库的使用以及普朗克平均吸收系数的计算

HITRAN数据库的使用一、HAPI第三方库的调用参考文章列表：1.https://blog.csdn.net/qq_38307826/article/details/1075560562.https://blog.csdn.net/Gou_Hailong/article/details/105586880总的来说就是从官网上下载HAPI第三方库，然后从pycharm中import，再然后就是下载数据，最后调用第三方库对数据进行必要的处理。当然，这个处理过程可能还会涉及到其他用于科学计算的第三方

2021-04-20 17:24:16 2290

原创热力学第二定律的局限性与对热寂说的反驳

摘自冯青《工程热力学》开尔文和克劳修斯在热力学第二定律的建立中做出了重大的贡献，但他们却把热力学第二定律做了不适当的推广，先后在1862年和1867年提出了错误的“热寂说”。他们把整个宇宙当做一个经历一系列连续变化的孤立系统来研究，并根据熵增原理得出结论：认为热总是不断地从高温物体传向低温物体，功总是自发地转变为热。这样，经过足够长的时间后，宇宙中的能量将由于不断贬值终于失去转化为其他能量的能力；这时，宇宙的温度将趋于一致，宇宙的熵趋于最大值，一切运动都将停息，宇宙将变成静止的热死亡的状态。这即是所谓的“

2021-04-04 17:02:58 1565

原创 ICEM/FLUENT学习笔记（三）

带边界层的非结构化网格导入fluent，选择二维求解器。初始化之后曲线收敛压力云图，与物理分析相符速度矢量图速度云图网格文件https://download.csdn.net/download/qq_42183549/15706516

2021-03-10 13:37:55 1123

原创 ICEM/FULENT学习笔记（二）

在导入网格文件时FLUENT报错，日志如下：Error at host: Build Grid: Aborted due to critical error.0: Cell is missing face: Cell 19, zone 9, (ID: 20): at location ( 1.89242e+01, 6.32246e+02, -8.55000e+01).Error at Node 0: Build Grid: Aborted due to critical error.Error:

2021-02-23 21:10:42 9917 1

原创 ICEM/FLUENT学习笔记（一）

ICEM为什么无法导出msh文件？ICEM在导出网格文件时无法找到已经保存的uns文件但事实上如果选中all file还是能看到UNS点打开后没有反应，同时日志出现：You may either select a superdomain or a set of structured domains. If you select structured domains, they will be converted into a superdomain using the “Structured -

2021-02-20 16:25:05 3883 2

原创 Python学习笔记（十七）

对应17章内容import requestsurl='https://api.github.com/search/repositories?q=language:python&sort=stars'headers={'Accept':'application/vnd.github.v3+json'}r=requests.get(url,headers=headers)print(f"Status code:{r.status_code}")response_dict=r.json()

2021-02-18 17:53:26 316 1

原创 Python学习笔记（十六）

前十六章要点总结方法：assertEqual(v,str)\setUp()\subplots()生成图片\plot(x,y,c=,alpha=)生成图表\set_title(‘str’,fontsize=)\set_xlabel(‘str’,fontsize=)\set_ylabel(‘str’,fontsize=)\tick_params(axis=‘both’,labelsize=)\use(seaborn)\scatter(x,y,c=,cmap=,edgecolors=,s=)\show()\s

2021-01-30 23:05:50 112 1

原创 Python学习笔记（十五）

十六章拓展plotly express 的渐变https://blog.csdn.net/tMb8Z9Vdm66wH68VX1/article/details/107421627?ops_request_misc=%7B%22request%5Fid%22%3A%22161175139116780264014887%22%2C%22scm%22%3A%2220140713.130102334.pc%5Fall.%22%7D&request_id=161175139116780264014887&

2021-01-30 23:04:49 274

原创 Python学习笔记（十四）

十六章 JSON过长的地址名不能使用分行filename='D:\\pycharmprograms\\codes\\chapter_16\\mapping_global_data_sets\\data\\eq_data_1_day_m1.json'导致：SyntaxError: EOL while scanning string literalpython非字符串中出现中文字符：data=pd.DataFrame( data=zip(lons,lats,titles,mags),col

2021-01-30 23:03:08 329 1

原创 Python学习笔记（十三）

十六章CSVax.set_title等函数中的标题应该使用英文单引号而不是双引号，否则继续输入中文时程序会误以为汉字字体错误。RuntimeWarning: Glyph 24180 missing from current font. font.set_text(s, 0.0, flags=flags)其他部分程序运行比较成功import csvimport matplotlib.pyplot as pltimport matplotlibfrom datetime import dateti

2021-01-25 19:08:49 467

原创 Python学习笔记（十二）

十五章直方图plotly并不支持.html之外的格式的输出Your filename d6.png didn’t end with .html. Adding .html to the end of your file.使用plotly时并不要求在使用时首先设置汉字字体。plotly 相关的博客Python数据可视化教程：基于Plotly的动态可视化绘图_aiwen2100-CSDN博客...

2021-01-25 19:07:34 203

原创 Python学习笔记（十一）

十五章拓展<ax.plot(x, y1, ‘-rv’,label=‘y1’)>这段代码里面的’-rv’包含了三个重要细节: <->表示实线表示红色表示下三角形。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-3tiyjMlP-1611572769259)(https://s3-us-west-2.amazonaws.com/secure.notion-static.com/19fb111b-7793-48d7-a525-9fbb2919a6f8/

2021-01-25 19:06:27 76

原创 Python学习笔记（十）

对应第十五章折线图部分import matplotlib.pyplot as pltsquares= [1,4,9,16,25]fig,ax = plt.subplots()ax.plot(squares)ax.set_title("平方数",fontsize=24)ax.set_xlabel("值",fontsize=14)ax.set_ylabel("平方",fontsize=14)ax.tick_params(axis='both',labelsize=14)plt.show(

2021-01-25 19:05:35 2752

原创 Python学习笔记（九）

前十章要点整理方法：init()\read()\readlines()\replace(’’)\split()函数：input()\int()\float()\remove()\super()\open(path,file)\close()\json.dump()\json.load(value,file)语句：break\continue\def*(元组)**(字典)\try-except\类的命名使用驼峰命名法，类中的方法、属性应使用小写字母加下划线的形式...

2021-01-25 19:04:14 82

原创 Python学习笔记（八）

对应第十章如果文件内容不是逐行输出则会输出列表形式。try: with open('D:\\pycharmprograms\\pi.txt') as f: contents=f.readlines()except FileNotFoundError: print("FileNotFoundError")else: print(contents)>>>['3.1415926535\n', ' 8979323846']json数据格式的列

2021-01-16 20:48:02 106

原创 Python学习笔记（七）

对应第九章面向对象编程语言：类: 一个模板, (人类)—是一个抽象的，没有实体的对象: (eg: 张三，李四)属性: (表示这类东西的特征, 眼睛，嘴巴，鼻子)方法： (表示这类物体可以做的事情， eg: 吃饭，睡觉，学习)类的意义可以参考：init前时双下划线__而不是_。class Restaurant: def _init_(self,restaurant_name,cuisine_type): self.restaurant_name=rest

2021-01-16 20:46:24 131

原创 Python学习笔记（六）

对应第八章内容当函数有默认值时未必要使用关键字实参def name(first_name,last_name='biden'): print(f"{first_name} {last_name}")name('joe','wisle')>>>joe wisle函数内变量的生存期定义域只在函数内，主程序中可有同名变量def name(first_name,last_name='biden'): full_name=f"{first_name} {last_nam

2021-01-16 20:45:21 129

原创 Python学习笔记（五）

对应第七章name=['joe','jojo','joe','jojo']name2=set(name)print(name)print(name2)name3=list(name2)print(name3)print(name3[0])['joe', 'jojo', 'joe', 'jojo']{'joe', 'jojo'}['joe', 'jojo']joe在条件判据中active与active==Ture等价active=Trueif active==True:

2021-01-16 20:44:14 51

原创 Python学习笔记（四）

前五章总结方法：title()\upper\lower\rstrip\lstrip\stripappend\insert\pop\remove\sort\reverse\range\items\key\value\函数：sorted\min\max\sum\set语句：del

2021-01-08 21:52:02 59

原创 Python学习笔记（三）

name={'j1':'john','j2':'joe','j3':'johnson','j4':'jack'}a=name.items()print(a)dict_items([('j1', 'john'), ('j2', 'joe'), ('j3', 'johnson'), ('j4', 'jack')])但是items返回的值并不支持输出name={'j1':'john','j2':'joe','j3':'johnson','j4':'jack'}a=name.items()print

2021-01-08 21:51:06 79

原创 Python学习笔记（二）

car='bmw'print(car=='bmw')“==、and、or”运算是布尔运算，有返回值True/Falsea=[]if a: print("1")else: print("2")可以用if a 做判断，但是if a 本身并不返回值a=[]print(if a)File "<ipython-input-9-90ac662d1df0>", line 2 print(if a) ^SyntaxError: invalid

2021-01-08 21:50:22 233

原创 Python学习笔记（一）

sorted（）函数用法详解：https://www.notion.so/python-5e85d913a5c94c4cb31d9250af6ef6b9#89280ae9db8443fa868857dcf97fd972print(first_name)File "<stdin>", line 1print(first_name)^IndentationError: unexpected indentpython的缩进是有语法意义的，首行的空格会引起歧义。full_name=f"

2021-01-08 21:49:10 226