【技术实现步骤摘要】
本专利技术属于力学建模与分析领域，具体涉及一种基于里兹法的功能梯度材料的半解析分析方法。
技术介绍
1、随着飞行器设计技术的发展，飞行器的飞行速度的极限也逐渐被突破，很多飞行器都达到了高超声速的范围。而高速飞行给飞行器设计又带来了很多新的挑战。在高速飞行的过程中，飞行器的表面与空气会产生剧烈的摩擦，产生大量的气动热，使得飞行器表面可达到上千度的高温。为了抵御如此高的温度，通常都需要在表面铺设热防护材料，而不同材料之间由于物理属性存在突变，极易导致应力集中、剥离、分层等现象。针对这一现象，功能梯度材料(functionally graded materials,fgm)的思想便应运而生，即两种不同的材料通过指定的连续函数的分布规律进行混合，使得物理属性可以实现连续变化。而功能梯度材料这种新的材料形式，与传统的均匀材料相比，有着很大不同：
2、首先，在功能梯度材料中，材料属性不再是一个常值，而是随着空间位置而发生变化，需要先将功能梯度材料的材料属性通过两种组分材料的属性表示出来。
3、其次，由于材料属性在厚度方向上不均匀，导致在受载下，对称面的应变并不为0，即中性面并不在对称面上，经典板理论中，横向位移为0的假设并不成立，不能直接使用现有的均匀板理论。
4、此外，功能梯度材料的热传导系数也同样随空间位置而变化，稳态热传导微分方程中的热传导系数不再为常数，微分方程的性质从线性微分方程变化为非线性微分方程，需要对功能梯度材料的温度场进行表征。
5、还有，功能梯度材料由于材料属性并不是均匀的
6、目前国内外对功能梯度材料的力学性能的研究主要可分为解析/半解析法和数值解法。其中，数值解法包括有限元法，无网格法，广义微分求积法等。数值解法中，由于功能梯度材料沿厚度方向材料属性的变化很大，为了更精确地模拟功能梯度材料的属性，需要在厚度方向上离散的节点很多，网格量很大，计算效率很低。而解析/半解析法中，目前大多数研究都是通过平衡方程的形式推导出偏微分方程，再通过galerkin法离散得到常微分方程组或代数方程组进行求解，而平衡方程的建立需要对该力学系统更为熟悉，且不容易快速扩展，galerkin法属于加权余量法，得到的未知系数物理意义不清晰。
技术实现思路
1、为了解决上述问题，本专利技术提供一种基于里兹法的功能梯度材料的半解析分析方法。
2、本专利技术的技术方案如下：
3、一种基于里兹法的功能梯度材料的半解析分析方法，所述方法包括以下步骤：
4、s1获取功能梯度材料的位移边界条件，基于所述位移边界条件获得位移试函数；其中位移试函数由位移系数和位移基函数的乘积构成；
5、s2通过计算模块计算所述功能梯度材料的各种能量；本专利技术将不同种类的能量进行了模块法处理，可以分门别类地单独计算各种能量，并在各个能量计算模块中传递共用的物理量信息，并根据功能梯度材料的特点，将每一种能量形式进行特殊处理，并将各能量计算模块有机组合在一起，形成信息交互。而在力学系统中，与结构自身关系最密切的就是应变能，而对于功能梯度材料，由于其工作环境经常在高温中，热对结构的影响也很重要。
6、s3将s2计算的各种能量通过哈密顿原理联系起来，得到总的能量泛函，实际的可存在的物理系统的能量泛函的变分为0；
7、s4根据变分法原理，总的能量泛函的变分为0，等价于，对各个位移系数求得的偏导数为0，由此可得到关于位移系数的常微分方程组或代数方程组，直接求解该方程组可得到位移系数；
8、s5基于s4得到的位移系数求解位移试函数，得到所述功能梯度材料的所有响应物理量。
9、进一步地，所述s2中计算模块包括热计算模块和应变能计算模块。
10、进一步地，所述热计算模块具体为：
11、在给定的温度边界条件下，通过求解非线性常微分方程得到内部温度场与空间位置的表达式；
12、通过新的弹性模量和热膨胀系数，计算热应力；
13、将热应力视为压缩载荷，在功能梯度材料总的计算域内进行积分，得到热应力带来的压缩载荷所作的功，即为热应力对功能梯度材料的能量影响。
14、进一步的，所述热计算模块还包括：
15、如果需要考虑材料属性与温度相关的特点，则需要将求得的温度场的表达式，代入材料属性随温度变化的函数关系式中，得到新的材料属性，并将新的与温度相关的材料属性输出出来。
16、进一步的，所述应变能计算模块具体为：
17、首先根据位移试函数计算功能梯度材料的几何参数，所述几何参数包括斜率和曲率；
18、根据计算得到的功能梯度材料内部的材料属性与空间位置的关系，得到物理中面的位置；
19、根据几何参数和物理中面位置，得到所述功能梯度材料内部的应变；
20、根据功能梯度材料的物理属性，得到所述功能梯度材料内部的应力；
21、根据上两个步骤得到的应力与应变的乘积，在功能梯度材料的计算域内进行积分，得到所述功能梯度材料的应变能。
22、进一步的，所述物理中面的位置为应变为0处的平面的位置。
23、本专利技术的技术效果：
24、(1)流程化，模块化，将各项因素对力学系统的影响通过能量的形式表示，可以方便得到功能梯度材料中；
25、(2)建模简单，通过基本的物理概念便可建立得到能量表达式，进而得到常微分方程组或代数方程组，完成求解；
26、(3)得到的结果形式简单，求解速度快，可以更方便得到各物理量对力学性能的影响；
27、(4)可以考虑热环境对功能梯度材料的影响，包括热应力带来的载荷和温度对材料属性的影响。
本文档来自技高网...【技术保护点】
1.一种基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述方法包括以下步骤：
2.根据权利要求1所述的基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述S2中计算模块包括热计算模块和应变能计算模块。
3.根据权利要求2所述的基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述热计算模块具体为：
4.根据权利要求2所述的基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述热计算模块还包括：
5.根据权利要求2所述的基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述应变能计算模块具体为：
【技术特征摘要】
1.一种基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述方法包括以下步骤：
2.根据权利要求1所述的基于里兹法的功能梯度材料的半解析分析方法，其特征在于，所述s2中计算模块包括热计算模块和应变能计算模块。
3.根据权利要求2所述的基于里兹法的功能梯度材...
【专利技术属性】
技术研发人员：王晓喆，李畅，万志强，杨超，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。