原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
温度场与热变形问题 9-1 温度场与热变形问题 9-2 温度场问题的基本方程 9-3 平面稳态温度场的有限元法 9-4 热变形的计算;2.1 直梁
2.1.1梁有限元模型
2.1.2节点位移与节点载荷
2.1.3单元刚度矩阵
2.1.4单元刚度矩阵的叠加
2.1.5边界条件
2.1.6工程实例
2.2 平面刚架
2.2.1有限元法基本思想节点位移与节点载荷
2.2.2单元刚度矩阵
2.2.3单元刚度矩阵的坐标变换
2.2.4总的刚度矩阵叠加
2.2.5位移法基本方程
2.3工程实例
2.2.1有限元法基本思想节点位移与节点载荷
2.2.2单元刚度矩阵
2.2.3单元刚度矩阵的坐标变换
2.2.4总的刚度矩阵叠加
2.2.5位移;9-1 温度场与热变形问题;;设微元在dt内，温度升高为：
相应所积蓄的热量为：
同一时间内，微元体沿x方向传入和传出的热量之差，即净热量为：
类似，y，z方向的净热量：
即传入微元体的净热量为：
由热传导定律：热流密度与温度梯度成正比，而方向相反，即：
代入上式得传入微元体净热量为：;设微元体内有热源，其热源密度为Q(x,y,z,t)，则该热源在dt内所共给的热量为：
据热平衡得一般热传导微分方程：;整理得：
满足上述热传导方程的解有无限多个，为了确定真实的温度场，必须知道物体初始瞬态的温度分布，即初始条件，称为第一类边界条件
同时，还需知道物体表面与周围介质间进行热交换的规律，即边界条件，称为第二类边界条件。;1、三维瞬态热传导方程及边界条件
2、二维稳态热传导方程及边界条件
;9-3 平面稳态温度场的有限元法;2、平面稳态温度场的泛函
求满足平面温度场方程及边界条件的温度场T(x,y),设k为常数
据变分原理，此问题等价于求泛函J[T(x,y)]的极值函数，参考相关教材，可得上述热传导作为欧拉方程的相应泛函：
;3、温度场单元分析
图示求解域离散为若干三角形单元，含有边界的单元，称为边界单元，任取一个单元i,j,k,如图。
A、温度插值函数
在边界线(如ij)上的任一点的温度T，可用两个端点的节点温度线性插值表示：
;B、单元温度刚度矩阵
从温度场插值函数可知，温度场已离散到全部节??上，即求温度场实际是求节点的温度值。因而，泛函式实际已成为描述未知节点温度的多元函数，而不是温度场T(x,y)的函数，即问题转化为求多元函数的极值
设求解域有n个节点温度未知量，则泛函J[T(x,y)]转化为 的形式，极值条件为：
设单元只有三节点温度，jk为边界，将温度插值函数代入前述的泛函，并求导得极值条件：
;上式第一部分为内部单元的温度刚阵：
对于内部单元的温度刚阵，i,j,k三点轮换，记为矩阵形式：
第二部分：
记为矩阵形式：
两部分相加可得边界单元的温度刚阵：
;3、整体温度场方程
为n个线性方程组，对于每个方程而言，是对绕节点m的所有单元求和，如图，节点5，则绕节点5的单元为1，2，3，而其它单元不含节点5，即它们的泛函对 的偏导为0，可不考虑，即
如单元1，3为边界单元，则按边界单元刚阵计算；如单元2为内部单元，则按内部单元刚阵计算。
如此整理可得整体代数方程组：
对于其他带热源的稳态温度场或三维温度场计算其方法相似。;9-4 热变形的计算;设热变形引起的初应变：
则考虑初应变情况的弹性方程(如平面应力问题)：
应力方程：
对比不考虑初应变的应力方程：
刚度方程：
;
文档评论（0）