- 博客(79)
- 收藏
- 关注
原创 传热学一些“数”和意义
Nu与Bi的表达式相同,但是意义是无量纲的h。它们表达式里的长度取值不同,比如同样一个平板,Bi的L是厚度,Nu是长度,因为Bi面向固体,λ为固体的,Nu的λ是气体的。管子的特征长度是直径,板子如果单面受热特征长度就是一个厚度,如果双面受热特征长度就是半个厚度。速度边界层厚度/温度边界层厚度 (粘性扩散能力/热扩散能力)物体单位面积上的导热热阻/单位表面积上的对流换热热阻。
2024-10-08 16:00:00 200
原创 关于直肋的效率
这个的分子的导数小于0,分子的零点明显为0,因此在x>0时该导数为负,因此随mH的减小,效率提升。可以如此理解,由效率表达式上面的分子在x趋近于∞时为1,可以看做一个反比例函数。肋片越矮胖,肋片各个部分的温度越趋近于肋根温度,效率越高。对于截面为圆柱形的,效率随半径的增大而增大。假如固定其他条件,当λ变大,mH减小,对。效率随材料的导热系数的增大而增大。效率随对流换热系数的增大而减小。效率随肋片的增高而减小。进行求导,得到导数为。
2024-10-08 15:50:24 675
原创 k=1.4,R=287.06. 总温和临界声速的关系
T*c_cr100183.003643101183.916385102103104105106107108109110191.93584111112113114115116117118119120121122123124125126127128129130131132133134135136137138139215.758113140141217.304786142143144145146147。
2024-10-08 08:46:18 117
原创 关于复变函数解析与可导
一般情况下,有lm(z)或|z|存在的函数不可导,但是也有例外,比如。一般看f(z)在定义域D内解析只看是否满足柯西-黎曼方程,即。两个式子同时成立,但是还要满足u和v在D内可微。
2024-10-06 18:24:56 992
原创 关于abaqus里一些问题的记录
为了帮助您直观地了解偏差系数的影响,Abaqus/CAE 会显示它将围绕与您输入的设置相对应的圆圈创建的元素数。指定最小尺寸控件有助于防止 Abaqus/CAE 在高曲率区域创建不必要的精细网格(如果这些区域对分析意图并不重要)。最后输出的结果没有力,可能是因为材料的参数给的不对,忘记了压强的单位是MPa,尤其注意密度单位t/mm^2。在进行布种时,会遇到最大偏离因子和最小尺寸因子,在帮助文档里,是这么解释。这里有一个恢复模块,不是回弹,是重新提交之前被中断的explicit任务。输入绝对最小元素大小。
2024-10-04 08:30:07 305
原创 传热学测温套管相关
公式的推导过程见传热学的书籍,假设套管长为H,外径为d,管厚δ。套管伸入大管子,大管子壁上温度为t0,温度计示数为tw,则外温是。,如果给定测量误差α和来流温度tf,则θ=αtf。对于一般的肋片,如果顶端绝热,则顶端的过余温度为。,如果已知m,可以由此解出H,代码如下。
2024-09-30 21:20:13 666
原创 关于特征方程的想法
解常微分方程所用到的特征方程或许和拉氏变换有关,比如解。,这分母正是特征方程,可以解出通解为。,这和由特征方程解出的一样,其中。
2024-09-18 14:24:13 954
原创 matlab实现简单的保角变换
用虚数的思想可以在虚坐标系内绘制圆,同样可以用虚数的想法将这个圆进行变换。用MATLAB绘制一个坐标在(1,1)的圆代码如下。
2024-09-04 14:21:10 400
原创 古典显示格式解一偏微分方程并绘制结果的彩色图
最后画图时如果用mesh(x,t,T),会显示一个三维彩色网格图如下。如果想要实心的格子,就用surf(x,t,T),效果如下。通过转动也可以显示出x-T或t-T二维图。以上公式的Latex代码。MATLAB代码如下。
2024-09-02 17:00:35 240
原创 离散的点进行积分
这里的代码表明,拟合后的结果可以直接转化为句柄,但并不是标准的形式,因此不能用integral直接积分。2,可以先进行多项式拟合,再积分。1,可以直接用梯形公式。
2024-09-01 09:15:47 272
原创 穿越沙漠小游戏的Python代码
本游戏以命令行的形式进行,玩家可以通过改变代码的方式改相关的参数,游戏背景请参考。这里假设在一条线上移动,天气随机决定。
2024-08-28 10:35:45 143
原创 logistic方程求解和模拟
其中N(t)为t时刻种群的数量,N为环境最大容纳的的种群数。k的取值通常与N有关,大约在1/N和4/N之间。logistic规律可以用以下式子表示。假设一些参数后的代码如下。
2024-08-25 13:53:23 246
原创 采用不高于3次的勒让德多项式拟合原函数
利用勒让德多项式进行拟合的区域是[-1,1],如果不是这个区域,比如是[a,b],利用。integrate为求积分的函数,见下。如果用MATLAB,则可以使用函数句柄。转化到[-1,1]。参考以下例题计算系数。
2024-08-24 14:41:11 145
原创 一个可以搜索页面内超链接的HTML页面
页面效果如上,含有标题,搜索框和一些超链接,在搜索框输入一些文字可以在下方显示含有这些文字的超链接。这里的搜索不仅可以中文,还可以是英文,数字。
2024-08-24 12:32:19 335
原创 用C#写一个随机音乐播放器
链接:https://pan.baidu.com/s/105w5fij6kkrUoWcjMuiBWA?这里的文件夹地址是直接Ctrl+Shift+C复制过来的,注意要去掉引号。实际是随机打开一个文件,鉴于VS内置的库不能播放除了wav以外的音频。form1中namespce里的代码如下。
2024-08-19 13:34:15 635
原创 两种用MATLAB绘制色块的方法
这种方式本质上是创造一个图片,因此显示时没有坐标轴,不过由于仅仅填充像素,当尺寸设置过小,图像也过小。绘制色块首先可以想到用填充像素的方式。另外可以想到填充用坐标线围成的网格。
2024-08-19 12:58:54 366
原创 y=λsin(πx)分岔的研究
加密步长,取r区间为[0.4,0.9],绘制feigenbaum图。用如下代码对收敛的λ的值进行探究,这里的r代表λ。使用如下的迭代格式,λ为可变的参数。可以观察到从0.9开始就开始分岔。
2024-07-25 22:21:27 224
原创 简单的差分格式解一个一元二阶常微分方程
将a到b的区域平均分成m份,解以上的式子可以看做解一个三对角方程组。这里也包含了追赶法求三对角矩阵的函数。对于一个如下的一元二阶常微分方程。其中二阶导数可以换成差分格式。之后展开成的矩阵形式,取。这里是左右两边同时乘以。
2024-07-06 11:46:41 749 1
原创 C++创建一个可以在不同源文件之间访问的二维数组
这里定义了数组的大小,定义二维数组时要在前面加extern。首先需要创建一个源文件(如Moudle.h),代码如下。之后,可以在一个源文件里给u赋值。注意这里u的声明要放函数外面。在主文件里,代码如下。
2024-07-05 13:53:05 116
原创 MATLAB将两个折线图画在一个图里
输入行数和列数,点击开始填入数据,其中第一列为x值,后面几列,每一列都是y坐标值,填好后点击画在同一张图里即可。点击置零就把所有数变成0,另外也可以选择节点样式。.mlapp格式的文件如下。
2024-07-02 16:12:57 600 1
原创 关于一维,二维正态分布的绘制
完全正相关或完全负相关,即一个变量的增加(减少)总是伴随着另一个变量按确定比例的增加(减少)。,可以在保持各自分布形态不变的情况下,描述两个随机变量之间不同的依赖关系。二维正态分布的概率密度函数(PDF)可以用两个一维正态分布的随机变量。是独立的,其联合分布可以分解为各自边际分布的乘积;和两个随机变量标准差的乘积的比值,即。之间的线性关系强度和方向,由相关系数。对于二维正态分布,其联合概率密度函数。若要算某一点的概率密度,代码如下。,以及它们之间的相关系数。绘制二维正态分布代码如下。本文采用AI辅助创作。
2024-07-02 08:04:53 1014
原创 强度理论的相当应力在坐标图里的表示
以下的图均只考虑二维,且没有切应力,有切应力的情况可以转换成只有两个正应力的情况,另外,这里画出的图均认为相当应力就是许用应力。两种强度理论的图均明显有一个转角,这是因为主应力的转变,当。从图形上看,在极坐标系下的和在直角坐标系下的相同。这里绘制畸变能理论的图,公式为。第一,二强度理论用于脆性材料。绘制最大拉应变理论的代码如下。第三,四强度理论用于脆性材料。
2024-07-01 13:00:04 181
原创 学习数值方法解常微分方程的笔记
比如y''' + 2y'' + 2y' + y = 80,化为y1' = y2, y2' = y3, y3' = 80 - 2y3 - 2y2 - y1,这个80 - 2y3 - 2y2 - y1相当于一阶里的f(t,y)工程上在满足精度要求的前提下,多用4阶龙格库塔法解常微分方程。取步长0.001,步数10,C语言代码如下。对于多阶常微分方程,可以逐级化为一阶的。与精确值的误差几乎没有。
2024-04-27 11:30:47 276 1
原创 单片机小项目——直流电机+按键
不过在未按下按键时,电动机仍然有一定的运转,而且在按下前15个按键的瞬间电动机仍要加速转动,猜想是矩阵键盘的扫描与电机共用一个IO口,互相影响了电平。使用4*4键盘控制,第一个到第15按键分别对应数码管显示1到F,计划在按下第16的按键时电动机运转数码管显示的秒数,代码如下。利用普中单片机的代码实现在按下第k个独立按键时,直流电机运作k秒。
2024-04-20 22:32:46 616
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人