一种梯形直波纹板生产方法技术

一种梯形直波纹板生产方法，将梯形直波纹板分为多个计算单元；建立计算单元的应变能方程；求解计算单元的周期函数；基于计算单元的周期函数计算计算单元的等效板拉伸刚度、等效板弯曲刚度及等效板拉弯耦合；利用分析软件基于计算单元的等效板拉伸刚度、等效板弯曲刚度及等效板拉弯耦合计算梯形直波纹板的拉伸刚度、拉弯刚度及拉伸耦合。本发明专利技术通过具有相应等效板刚度的平板代替梯形直波纹板，从而大大减少单元数量和全局分析时有限元模型的总自由度，通过梯形直波纹板截面分析得到等效刚度。通过有限元分析或封闭形式解析解来计算梯形直波纹板的整体性能，极大地提高了梯形直波纹板性能计算的效率，进而提高了梯形直波纹板的生产效率。

【技术实现步骤摘要】
一种梯形直波纹板生产方法
本专利技术涉及生产方法
，尤其涉及一种梯形直波纹板生产方法。
技术介绍
板壳是由两个接近的曲面或平面构成的薄三维体。在工程应用中，薄板缺乏抗剪、抗压和抗弯曲载荷的能力，可通过波纹结构进行弥补。与平板相比，波纹板在平面上沿某个方向存在重复曲率，其宏观行为主要受此曲率控制。一般来说，波纹板的弯曲不能脱离其在一定方向上的拉伸，这使得其结构分析比平板更复杂。在梯形直波纹板的生产过程中，需要先对梯形直波纹板进行结构设计，并对梯形直波纹板的各项性能进行计算，当性能满足要求后再进行生产。现有技术中，通常将设计好的梯形直波纹板结构直接进行有限元分析，然而，这种方法需要规模较大的有限元模型，且计算量大，效率低。因此，本专利技术公开了一种梯形直波纹板生产方法，通过具有相应等效板刚度的平板代替梯形直波纹板，从而大大减少单元数量和全局分析时有限元模型的总自由度，通过梯形直波纹板截面分析得到等效刚度。然后，通过有限元分析或封闭形式解析解来计算梯形直波纹板的整体性能，极大地提高了梯形直波纹板性能计算的效率。
技术实现思路
针对现有技术存在的上述不足，本专利技术要解决的技术问题是：如何减少单元数量和全局分析时有限元模型的总自由度，从而提高梯形直波纹板性能计算的效率，进而提高梯形直波纹板的生产效率。为解决上述技术问题，本专利技术采用了如下的技术方案：一种梯形直波纹板生产方法，包括根据性能需求设计梯形直波纹板结构，计算所述梯形直波纹板结构的拉伸刚度及弯曲刚度，当所述梯形直波纹板的拉伸刚度及拉弯刚度满足所述梯形直波纹板的性能需求时，按照设计的梯形直波纹板结构制造所述梯形直波纹板，其中，计算所述梯形直波纹板的拉伸刚度、弯曲刚度及拉弯耦包括如下步骤：将所述梯形直波纹板分为多个计算单元；建立所述计算单元的应变能方程；基于所述计算单元的应变能方程求解所述计算单元的周期函数；基于所述计算单元的周期函数计算所述计算单元的等效板拉伸刚度及等效板弯曲刚度；利用分析软件基于所述计算单元的等效板拉伸刚度及等效板弯曲刚度计算所述梯形直波纹板的全局行为。梯形直波纹板的全局应变包括梯形直波纹板的应力和应变。本专利技术通过具有相应等效板刚度的平板代替梯形直波纹板，从而大大减少单元数量和全局分析时有限元模型的总自由度，通过梯形直波纹板截面分析得到等效刚度。然后，通过有限元分析或封闭形式解析解来计算梯形直波纹板的整体性能，极大地提高了梯形直波纹板性能计算的的效率。优选地，建立所述计算单元的应变能方程包括以下步骤：建立每个所述计算单元的应变能方程其中J表示所述计算单元的应变能、G表示剪切模量、t表示厚度、a表示表面度量张量、μ为第一计算参数，μ＝v/(1-v)、γ表示拉伸应变、ρ表示弯曲应变、v表示泊松比，在本专利技术中，字符的下标为1表示此字符代表的量的方向为梯形波波纹的正方向，下标为3表示此字符代表的量的方向为垂直于梯形波波纹的正方向向上，下标为2表示此字符代表的量的方向为梯形波波纹正方向和垂直于梯形波波纹正方向向上两个方向的夹角方向，下标为连续的两个数字时，带有前一个下标数字的此字符与带有后一个下标数字的此字符相乘，<>表示求<>内的量关于X在-1/2至1/2间的积分，X为无量纲计算单元坐标，X＝x1/L，L表示一个完整的梯形波结构的长度，x表示在计算单元中心建立的坐标系的坐标。在本专利技术中，每个计算单元，包括了梯形波纹板中一个角以及构成这个角的两个边的一半，如图2所示，波纹的形状可由x3(X)＝Lφ(X)这一周期函数描述。斜率k和表面度量张量a的行列式分别可表示为：显然，对于梯形直波纹板，k和|a|为梯形常数。x3(X)可分解为式中：P为折点位置，上标内的括号表示梯形号，梯形号为1表示计算单元内左边的直线，为2表示右边的直线，如图3所示。描述波纹板位移的函数u可认为是计算单元坐标、坐标x和y的函数，即：式中：υ为有效板位移，为X的周期性函数，下标α为自然数，取值为1或2。根据周期性结构的一般理论，可以定义为在计算单元上的平均值：υi(x1,x2)＝<ui(X,x1,x2)>(4)在本专利技术中i为自然数，i＝1、2或3，显然：波纹板拉伸应变的主导项可表示为：弯曲应变的主导项变为：计算单元的应变能可写为：式中：μ＝v/(1-v)，v为泊松比，G为剪切模量。材料参数G,μ和厚度t是X的函数，通过观察式(6)和式(7)，不包含项，仅与相关，涉及为此可将总应变能分为三部分，其中J1对应于J2对应于J3包含常量部分。优选地，基于所述计算单元的应变能方程求解所述计算单元的周期函数包括：基于所述计算单元的应变能方程求解所述计算单元的第一周期函数：υ表示有效板位移、φ表示描述波纹形状的函数、c1表示第一积分系数、k表示斜率、Z表示第二积分系数，在本专利技术中，字符的下标为中间用逗号隔开的两个数字时，表示带有后一个下标数字的X对带有前一个下标数字的此字符求导；基于所述计算单元的应变能方程求解所述计算单元的第二周期函数：t表示梯形直波纹板的厚度；基于所述计算单元的应变能方程求解所述计算单元的第三周期函数：首先求解与有关的主导应变能为：周期函数可通过式(5)约束下最小化式(9)中的得到。通过引入Lagrange乘子λ2处理约束，相应的Euler-Lagrange方程为。以及边界条件：式中：方括号表示整个计算单元域两端值之间的差异，例如表示两个条件，一个是周期性边界条件另一个是连续性在本专利技术中，公式间用逗号隔开，则表示为不同的条件，例如式(11)中第二个条件即为由式(11)的第二个条件得到λ2＝0，因此将式(13)代入式(6)第二个方程，得到对式(14)两侧相对X积分，得到的显式表达式:类似地，使用Lagrange乘子λ1,λ3来处理式(5)中和的约束。相应的Euler-Lagrange方程为以及边界条件为建立式(17)中的第五个条件，需要两段之间角度变化的几何关系，如图4所示。在经典板理论(Kirchhoff-Love板理论)下构建波纹板模型时，剪切效应不会导致角度θ的变化，即Δθ＝0。δ为变分符号。一个计算单元内两块直板的角度变化定义为：式中：K为位置矢量矩阵，u为位移矢量矩阵，s(α)分别是计算单元内从段(1)和段(2)连接测量的长度，段(1)表示计算单元内左边的直板，段(2)表示计算单元内右边的直板。忽略高阶项的小位移，式(18)改写为其中：将式(3)代入式(19)和式(20)，得到边界条件在不连续点处用式(21)修改式(17)中的第五个条件得到式(17)的第二和第四个条件得到λ1＝λ3＝0和第一积分系数因此考虑式(22)的条件，在计算单元长度上积分(k×(23)-(24))。得到c3＝0(25)式(23)和式(24)可以简化为对式(26)积分得到考虑式(7)，重写式(28)如表示分段(1)的零阶近似应变，上标为数字加梯形号时，表示梯形号对应的分段处求括号外的数字的阶数的近似应变，边界条件式(21)下沿计算单元长度进行积分，得到不连续点处对式(28)进行评估，得到以及注意对式(29)积分并考虑c4，得到因此其中：k×(27)+(31)得到沿计算单元长度积分，得到注意到对式(34)两侧相对本文档来自技高网...

【技术保护点】
一种梯形直波纹板生产方法，其特征在于，包括根据性能需求设计梯形直波纹板结构，计算所述梯形直波纹板结构的拉伸刚度及弯曲刚度，当所述梯形直波纹板的拉伸刚度及拉弯刚度满足所述梯形直波纹板的性能需求时，按照设计的梯形直波纹板结构制造所述梯形直波纹板，其中，计算所述梯形直波纹板的拉伸刚度、弯曲刚度及拉弯耦包括如下步骤：将所述梯形直波纹板分为多个计算单元；建立所述计算单元的应变能方程；基于所述计算单元的应变能方程求解所述计算单元的周期函数；基于所述计算单元的周期函数计算所述计算单元的等效板拉伸刚度及等效板弯曲刚度；利用分析软件基于所述计算单元的等效板拉伸刚度及等效板弯曲刚度计算所述梯形直波纹板的全局行为。

【技术特征摘要】
1.一种梯形直波纹板生产方法，其特征在于，包括根据性能需求设计梯形直波纹板结构，计算所述梯形直波纹板结构的拉伸刚度及弯曲刚度，当所述梯形直波纹板的拉伸刚度及拉弯刚度满足所述梯形直波纹板的性能需求时，按照设计的梯形直波纹板结构制造所述梯形直波纹板，其中，计算所述梯形直波纹板的拉伸刚度、弯曲刚度及拉弯耦包括如下步骤：将所述梯形直波纹板分为多个计算单元；建立所述计算单元的应变能方程；基于所述计算单元的应变能方程求解所述计算单元的周期函数；基于所述计算单元的周期函数计算所述计算单元的等效板拉伸刚度及等效板弯曲刚度；利用分析软件基于所述计算单元的等效板拉伸刚度及等效板弯曲刚度计算所述梯形直波纹板的全局行为。2.如权利要求1所述的梯形直波纹板生产方法，其特征在于，建立所述计算单元的应变能方程包括以下步骤：建立每个所述计算单元的应变能方程

【专利技术属性】
技术研发人员：钟轶峰，罗丹，李帛书，邓兵，李奕豪，刘一钉，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆,50