【数学建模】高温作业专用服装设计(2018A)隐式差分推导

为方便计算,对区域进行离散化处理,采用隐式差分格式进行离散计算。隐式差分格式如图:
请添加图片描述

每层材料内部

对第 j j j层材料:

其中, λ j \lambda_j λj表示第 j j j层的热扩散率, c j c_j cj表示第 j j j层的比热容, ρ j \rho_j ρj表示第 j j j层的材料密度。

∂ T j ∂ t = λ j c j ρ j ∂ 2 T j ∂ x 2 \dfrac{\partial T_j}{\partial t} =\dfrac{\lambda_j}{c_j \rho_j}\dfrac{\partial^2 T_j}{\partial x^2} tTj=cjρjλjx22Tj

对方程左端采用向后差分格式,方程右端采用中心差分格式,有:

T j ( x i , t k + 1 ) − T j ( x i , t k ) Δ t = λ j c j ρ j T j ( x i + 1 , t k + 1 ) − 2 T j ( x i , t k + 1 ) + T j ( x i − 1 , t k + 1 ) ( Δ x ) 2 \dfrac{T_j(x_i,t_{k+1})-T_j(x_i,t_{k})}{\Delta t} =\dfrac{\lambda_j}{c_j \rho_j} \dfrac{T_j(x_{i+1},t_{k+1})-2T_j(x_i,t_{k+1}) + T_j(x_{i-1},t_{k+1})}{(\Delta x)^2} ΔtTj(xi,tk+1)Tj(xi,tk)=cjρjλj(Δx)2Tj(xi+1,tk+1)2Tj(xi,tk+1)+Tj(xi1,tk+1)

整理得:
− T j ( x i , t k ) = λ j Δ t c j ρ j ( Δ x ) 2 F ( X ) − T j ( x i , t k + 1 ) -T_j(x_i,t_{k}) = \dfrac{\lambda_j \Delta t}{c_j \rho_j (\Delta x)^2} F(X) - T_j(x_i,t_{k+1}) Tj(xi,tk)=cjρj(Δx)2λjΔtF(X)Tj(xi,tk+1)

F ( X ) = T j ( x i + 1 , t k + 1 ) − 2 T j ( x i , t k + 1 ) + T j ( x i − 1 , t k + 1 ) F(X) = T_j(x_{i+1},t_{k+1})-2T_j(x_i,t_{k+1}) + T_j(x_{i-1},t_{k+1}) F(X)=Tj(xi+1,tk+1)2Tj(xi,tk+1)+Tj(xi1,tk+1)
F ( X ) = r j F(X) = r_j F(X)=rj

整理得:

T j ( x i , t k ) = − r j T j ( x i + 1 , t k + 1 ) − r j T j ( x i − 1 , t k + 1 ) + ( 1 + 2 r j ) T j ( x i , t k + 1 ) T_j(x_i,t_k) = -r_jT_j(x_{i+1},t_{k+1})-r_jT_j(x_{i-1},t_{k+1})+(1+2r_j)T_j(x_i,t_{k+1}) Tj(xi,tk)=rjTj(xi+1,tk+1)rjTj(xi1,tk+1)+(1+2rj)Tj(xi,tk+1)

左边界条件

T j ( x , 0 ) = 37 ° T_j(x,0) = 37\degree Tj(x,0)=37°
此前求得的左边界条件:

− k 1 ∂ T 1 ∂ x + k o u t T 1 ∣ x = 0 = k o u t T o u t -k_1 \dfrac{\partial T_1}{\partial x} + k_{out}T_1|_{x=0} = k_{out}T_{out} k1xT1+koutT1x=0=koutTout

对其进行离散化处理:

− k 1 T ( 2 , k ) − T ( 0 , k ) 2 Δ x = k o u t ( T o u t − T ( 1 , k ) ) -k_1 \dfrac{T_(2,k) - T(0,k)}{2\Delta x} = k_{out}(T_{out} - T(1,k)) k1xT(2,k)T(0,k)=kout(ToutT(1,k))

整理可得:

T ( 0 , k ) = 2 Δ x k o u t T o u t k 1 − 2 Δ x k o u t T ( 1 , k ) k 1 + T ( 2 , k ) T(0,k) = \dfrac{2\Delta xk_{out}T_{out}}{k_1} - \dfrac{2\Delta xk_{out}T(1,k)}{k_1} + T(2,k) T(0,k)=k1xkoutToutk1xkoutT(1,k)+T(2,k)

根据每层材料内部,有

T j ( 1 , k − 1 ) = − r 1 T j ( 2 , k ) − r 1 T j ( 0 , k ) + ( 1 + 2 r 1 ) T j ( 1 , k ) T_j(1,k-1) = -r_1T_j(2,k)-r_1T_j(0,k)+(1+2r_1)T_j(1,k) Tj(1,k1)=r1Tj(2,k)r1Tj(0,k)+(1+2r1)Tj(1,k)

联立得:

T j ( 1 , k − 1 ) = − r j T j ( 2 , k ) − r j ( 2 Δ x k o u t T o u t k 1 − 2 Δ x k o u t T ( 1 , k ) k 1 + T ( 2 , k ) ) + ( 1 + 2 r j ) T j ( 1 , k ) T_j(1,k-1) = -r_jT_j(2,k)-r_j(\dfrac{2\Delta xk_{out}T_{out}}{k_1} - \dfrac{2\Delta xk_{out}T(1,k)}{k_1} + T(2,k))+(1+2r_j)T_j(1,k) Tj(1,k1)=rjTj(2,k)rj(k1xkoutToutk1xkoutT(1,k)+T(2,k))+(1+2rj)Tj(1,k)

T j ( 1 , k − 1 ) + r 1 2 Δ x k o u t T o u t k 1 = − 2 r 1 T j ( 2 , k ) + ( 1 + 2 r 1 + r 1 2 Δ x k o u t k 1 ) T j ( 1 , k ) T_j(1,k-1)+r_1\dfrac{2\Delta xk_{out}T_{out}}{k_1} = -2r_1T_j(2,k)+(1+2r_1+r_1\dfrac{2\Delta xk_{out}}{k_1})T_j(1,k) Tj(1,k1)+r1k1xkoutTout=2r1Tj(2,k)+(1+2r1+r1k1xkout)Tj(1,k)

右边界条件

此前求得的右边界条件:

k 4 ∂ T 1 ∂ x + k o u t T 4 ∣ x = x s k i n = k s k i n T s k i n k_4 \dfrac{\partial T_1}{\partial x} + k_{out}T_4|_{x=x_{skin}} = k_{skin}T_{skin} k4xT1+koutT4x=xskin=kskinTskin

对其进行离散化处理:

k 4 T ( N + 1 , k ) − T ( N − 1 , k ) 2 Δ x = k s k i n ( T s k i n − T ( N , k ) ) k_4 \dfrac{T_(N+1,k) - T(N-1,k)}{2\Delta x} = k_{skin}(T_{skin} - T(N,k)) k4xT(N+1,k)T(N1,k)=kskin(TskinT(N,k))

整理可得:

T ( N − 1 , k ) = T ( N + 1 , k ) − 2 Δ x k s k i n T s k i n k 4 + 2 Δ x k s k i n T ( N , k ) k 4 T(N-1,k) = T(N+1,k) -\dfrac{2\Delta xk_{skin}T_{skin}}{k_4} + \dfrac{2\Delta xk_{skin}T(N,k)}{k_4} T(N1,k)=T(N+1,k)k4xkskinTskin+k4xkskinT(N,k)
根据每层材料内部,有

T j ( N , k − 1 ) = − r 1 T j ( N + 1 , k ) − r 1 T j ( N − 1 , k ) + ( 1 + 2 r 1 ) T j ( N , k ) T_j(N,k-1) = -r_1T_j(N+1,k)-r_1T_j(N-1,k)+(1+2r_1)T_j(N,k) Tj(N,k1)=r1Tj(N+1,k)r1Tj(N1,k)+(1+2r1)Tj(N,k)

联立得:

T j ( N , k − 1 ) = − r j T j ( N + 1 , k ) − r j ( T ( N + 1 , k ) − 2 Δ x k s k i n T s k i n k 4 + 2 Δ x k s k i n T ( N , k ) k 4 ) + ( 1 + 2 r ) T j ( N , k ) T_j(N,k-1) = -r_jT_j(N+1,k)-r_j(T(N+1,k) -\dfrac{2\Delta xk_{skin}T_{skin}}{k_4} + \dfrac{2\Delta xk_{skin}T(N,k)}{k_4})+(1+2r)T_j(N,k) Tj(N,k1)=rjTj(N+1,k)rj(T(N+1,k)k4xkskinTskin+k4xkskinT(N,k))+(1+2r)Tj(N,k)

T j ( N , k − 1 ) − 2 Δ x k s k i n T s k i n k 4 = − 2 r j T j ( N + 1 , k ) + ( 1 + 2 r j − r j 2 Δ x k s k i n k 4 ) T j ( N , k ) T_j(N,k-1)-\dfrac{2\Delta xk_{skin}T_{skin}}{k_4} = -2r_jT_j(N+1,k)+(1+2r_j-r_j\dfrac{2\Delta xk_{skin}}{k_4})T_j(N,k) Tj(N,k1)k4xkskinTskin=2rjTj(N+1,k)+(1+2rjrjk4xkskin)Tj(N,k)

N = L N=L N=L

本文参考自【数学建模】2018年国赛A题详解(一)

更多建模相关信息看Github仓库(包括真题练习)

更多数学建模信息

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://xiahunao.cn/news/3246642.html

如若内容造成侵权/违法违规/事实不符,请联系瞎胡闹网进行投诉反馈,一经查实,立即删除!

相关文章

RFID(NFC) CLRC663非接触读取芯片GD32/STM32 SPI读取

文章目录 基本介绍硬件配置连接硬件连接详解程序代码代码解释 基本介绍 CLRC663 是高度集成的收发器芯片,用于 13.56 兆赫兹的非接触式通讯。CLRC663 收发器芯片支 持下列操作模式 • 读写模式支持 ISO/IEC 14443A/MIFARE • 读写模式支持 SO/IEC 14443IB • JIS X…

全网超详细Redis主从部署(附出现bug原因)

主从部署 整体架构图 需要再建两个CentOs7,过程重复单机部署 http://t.csdnimg.cn/zkpBE http://t.csdnimg.cn/lUU5gLinux环境下配置redis 查看自己ip地址命令 ifconfig 192.168.187.137 进入redis所在目录 cd /opt/software/redis cd redis-stable 进入配置文件 vim redi…

【JavaEE精炼宝库】 初识网络原理——网络通信基础 | 协议

文章目录 一、网络发展史1.1 独立模式:1.2 网络互连:1.3 局域网(LAN):1.4 广域网(WAN): 二、网络通信基础2.1 IP地址:2.2 端口号: 三、协议3.1 协议的概念&am…

Python基础语法篇(上)

Python基础语法(上) 一、基知二、基本数据类型(一)标准数据类型(二)数据类型转换 三、字符串基本操作(一)字符串的索引和切片(二)字符串的拼接 三、运算符四、…

在golang中Sprintf和Printf 的区别

最近一直在学习golang这个编程语言,我们这里做一个笔记就是 Sprintf和Printf 的区别 fmt.Sprintf 根据格式化参数生成格式化的字符串并返回该字符串。 fmt.Printf 根据格式化参数生成格式化的字符串并写入标准输出。由上面就可以知道,fmt.Sprintf返回的…

AI第二课堂第一次笔记

conda的使用 在输入cmd进入终端后,使用命令 conda create -n env_name python3.10 创建环境 命令 conda activate env_name 打开环境,如:使用 conda deactivate退出指令 2.python一些常见操作 python中的文件打开与关闭 调开源的库 p…

C++--lambda表达式

介绍 一个lambda表达式表示一个可调用的代码单元。我们可以将其理解为一个未命名的内联函数。和函数类型,lambda有一个返回值,一个参数列表和一个函数体,但比函数多一个捕获列表。具体形式如下: [捕获列表](参数列表) ->返回值类型 {函数体}其中:捕获列表:可以捕获定义lam…

NineData全面支持PostgreSQL可视化表结构设计

“PostgreSQL 是最像 Oracle 的开源关系型数据库“,也正因为如此,很多企业都青睐 PostgreSQL,拿它当成 Oracle 的替代品。所以毫无疑问,目前 PostgreSQL 在企业中非常常见。 对于直接接触 PostgreSQL 的开发人员而言,…

60K起?“软件安全岗”比“网络安全岗”薪资高在哪里?

在网络世界的江湖中,“软件安全”与“网络安全”这两大“武林高手”都肩负着守护数字领域和平的重任。不过,眼尖的小伙伴们可能发现了,软件安全岗位的薪资待遇往往比网络安全岗位要丰厚那么一些,这到底是为啥呢?今天&a…

使用ETLCloud实现MySQL数据库与StarRocks数据库同步

在现代数据架构中,数据同步是保证数据一致性和分析准确性的关键步骤之一。本文将介绍如何利用ETLCloud技术实现MySQL数据库与StarRocks数仓数据库的高效数据同步,以及其在数据管理和分析中的重要性。 数据同步的重要性 在数据驱动的时代,企…

【紫光同创盘古PGX-Nano教程】——(盘古PGX-Nano开发板/PG2L50H_MBG324第十三章)蓝牙透传实验例程说明

本原创教程由深圳市小眼睛科技有限公司创作,版权归本公司所有,如需转载,需授权并注明出处(www.meyesemi.com) 适用于板卡型号: 紫光同创PG2L50H_MBG324开发平台(盘古PGX-Nano) 一:…

OZON成本低卖价高产品,OZON单价高产品

Top1 Crocs鞋扣配件 Джиббитсы набор украшения для обуви крокс значки на crocs клипсы аксессуары украшение для сабо бутылка вина бокал 商品id:1409545850…

第一百六十七节 Java IO教程 - Java Zip字节数组

Java IO教程 - Java Zip字节数组 校验和 Java在java.util.zip包中提供了一个Adler32类来计算数据字节的Adler-32校验和。 我们需要调用这个类的update()方法将字节传递给它。 在同一个包中还有另一个名为CRC32的类,它允许您使用CRC32算法计算校验和。 以下代码…

使用geoipupdate自动更新GeoIP数据库

一、 什么是 GeoIP? 通过在数据库中将地理位置和 IP 地址相互映射,软件程序便可以使用 IP 地址来确定其对应的地理位置,其中包括国家/地区、州/省、城市、邮政编码、纬度/经度、ISP、区号和其他信息。 很多软件都使用 MaxMind 的数据库对 IP…

选择Maya进行3D动画制作与渲染的理由

如果你对3D动画充满热情并追求成为专业3D动画师的梦想,你一定听说过Maya——近年来3D动画的行业标准。Maya被3D艺术家广泛使用,你是否想知道为什么Maya总是他们的首选?下面一起来了解下。 一、什么是Maya? 由Autodesk开发的Maya是…

昇思25天学习打卡营第24天|基于MindSpore的Diffusion扩散模型

Diffusion扩散模型 本文基于Hugging Face:The Annotated Diffusion Model一文翻译迁移而来,同时参考了由浅入深了解Diffusion Model一文。 关于扩散模型(Diffusion Models)有很多种理解,本文的介绍是基于denoising di…

浅谈Git

一:什么是 git git一种开源的分布式版本控制系统,可以有效、高速地处理从很小到非常大的项目版本管理。 下图是 git 的一个工作流程简图 二:什么是 分布/集中式版本控制系统 软件开发过程中,要解决多人协作的问题,需要…

FPGA CFGBVS 管脚接法

说明 新设计了1个KU040 FPGA板子,回来之后接上JTAG FPGA不识别。做如下检查: 1、电源测试点均正常; 2、查看贴片是否有漏焊,检查无异常,设计上NC的才NC; 3、反复检查JTAG接线是否异常,贴片是…

PostgreSQL 的 generic xlog 模块代码解读

为了描述方便,定义一些名词: wal/WAL -- Write Ahead Log,指PG中事务日志模块,这里所指的其实比较模糊,相关的函数,内存都算,自定义WAL就是指,用户可以决定自己存储在shared buffer …

数据库的约束条件和用户管理

约束条件: 主键:主键约束 primary key 用于标识表中的主键列的值,而且这个值是全表当中唯一的,而且只不能为null 一个表只能有一个主键。 外键:用来建立表与表之间的关系。确保外键中的值于另一个表的主键值匹配&a…