等效仿真方法、系统、计算机设备及计算机可读存储介质与流程

本发明涉及模态仿真
技术领域：
，尤其涉及一种考虑电机定子层叠结构的等效仿真方法、系统、计算机设备及计算机可读存储介质。
背景技术：
：近年来随着新能源汽车产业的飞速发展，新能源汽车的噪声问题也越来越引起人们的关注。与传统内燃机汽车不同，新能源汽车尤其是纯电动汽车的噪声主要来自驱动电机，其噪声总值虽然不大，但频率较高，容易引起驾驶员和乘客的烦恼。电机噪声主要由电磁力激励下的定子及机壳的振动产生，因此对电机定子进行准确的有限元建模和模态分析是研究电机噪声的关键。在电机设计过程中，为了减少定子的涡流损耗，定子铁芯通常由很薄的硅钢片沿轴向堆叠而成，放置在定子槽中的绕组则由一根根铜导线绕制而成，因此定子铁芯和绕组均是层叠结构而非实体结构，这样就给电机定子的模态仿真带来了很大的困难。现有的电机定子模态仿真方法中，大部分是将定子铁芯和绕组当成实体结构来建模，然后分别赋予硅钢片和铜的材料参数，仿真结果相比试验测试结果存在较大的偏差；还有一小部分也是将定子铁芯和绕组当成实体结构，分别对定子铁芯和绕组进行单独建模，然后分别赋予各向异性材料参数，这种方法需要分别确定定子铁芯和绕组的等效各向异性材料参数，需要确定的参数多，且要进行两次模态试验，工作量大，精度不易控制。此外，现有的方法均没有对定子有限元建模和网格划分过程进行控制，这样也会引入额外的误差。技术实现要素：本发明所要解决的技术问题在于，提供一种考虑电机定子层叠结构的等效仿真方法、系统、计算机设备及计算机可读存储介质，只需构建一个实体模型，可有效减少材料参数，易于控制，提高仿真精度。为了解决上述技术问题，本发明提供了一种考虑电机定子层叠结构的等效仿真方法，包括：建立定子铁芯三维有限元网格及定子绕组等效质量点，并关联所述定子铁芯三维有限元网格及定子绕组等效质量点以构成定子三维有限元模型；根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数，并采用所述定子三维有限元模型分析所述材料参数的取值对定子模态频率的影响规律；根据所述影响规律从所述材料参数中选取目标材料参数，并对所述目标材料参数的取值进行迭代修正，直到所述三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代；将修正后的目标材料参数作为定子的等效材料参数并赋予所述定子三维有限元模型以进行定子模态仿真。作为上述方案的改进，所述建立定子铁芯三维有限元网格及定子绕组等效质量点，并关联定子铁芯三维有限元网格及定子绕组等效质量点以构成定子三维有限元模型的步骤包括：构建定子铁芯三维有限元网格；在定子铁芯的几何中心位置构建定子绕组等效质量点；将所述定子绕组等效质量点与所述铁芯三维有限元网格中每个定子齿表面的网格节点相连接，以将定子绕组的质量平均分配到每个定子齿上；确定所述定子铁芯的等效密度及定子绕组等效质量点的质量。作为上述方案的改进，所述构建定子铁芯三维有限元网格的步骤包括：构建定子齿二维平面模型；对所述定子齿二维平面模型进行网格划分，生成定子齿二维平面网格；对所述定子齿二维平面网格进行周向旋转阵列，生成定子截面二维平面网格；将所述定子截面二维平面网格沿定子轴向拉伸，生成定子铁芯三维有限元网格。作为上述方案的改进，所述构建定子齿二维平面模型的步骤包括：根据电机定子的几何参数绘制定子齿二维平面模型，或者从预设定子三维有限元模型中切割出定子齿二维平面模型。作为上述方案的改进，所述确定定子铁芯的等效密度及定子绕组等效质量点的质量的步骤包括：根据公式ρ＝m/v，计算定子铁芯的等效密度，其中，ρ为定子铁芯的等效密度，m为实际定子铁芯的质量，v为定子铁芯三维有限元网格的体积；将实际定子绕组的质量作为所述定子绕组等效质量点的质量。作为上述方案的改进，所述根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数并采用定子三维有限元模型分析材料参数的取值对定子模态频率的影响规律的步骤包括：根据定子层叠结构的特点确定各向异性材料的独立材料参数；给所述独立材料参数赋予初始值；改变任一独立材料参数的取值并保持其他独立材料参数固定不变，采用所述定子三维有限元模型仿真分析定子模态频率随所述独立材料参数的变化规律；采用同样的方法分析定子模态频率随其他独立材料参数的变化规律。作为上述方案的改进，所述对目标材料参数的取值进行迭代修正直到三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代的步骤包括：给所述目标材料参数赋值，并采用所述三维有限元模型以对定子进行模态仿真，生成定子模态仿真结果；将所述定子模态仿真结果与试验结果进行对比，若定子模态仿真结果与试验结果之间的偏差超出目标范围，则重新对所述目标材料参数进行赋值，直到所述定子模态仿真结果与试验结果与试验结果之间的偏差在目标范围内时停止迭代。相应地，本发明还提供了一种考虑电机定子层叠结构的等效仿真系统，包括：模型构建模块，用于建立定子铁芯三维有限元网格及定子绕组等效质量点，并关联所述定子铁芯三维有限元网格及定子绕组等效质量点以构成定子三维有限元模型；规律分析模块，用于根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数，并采用所述定子三维有限元模型分析所述材料参数的取值对定子模态频率的影响规律；迭代修正模块，用于根据所述影响规律从所述材料参数中选取目标材料参数，并对所述目标材料参数的取值进行迭代修正，直到所述三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代；模态仿真模块，用于将修正后的目标材料参数作为定子的等效材料参数并赋予所述定子三维有限元模型以进行定子模态仿真。作为上述方案的改进，所述模型构建模块包括：网格构建单元，用于构建定子铁芯三维有限元网格；质量点构建单元，用于在定子铁芯的几何中心位置构建定子绕组等效质量点；关联单元，用于将所述定子绕组等效质量点与所述铁芯三维有限元网格中每个定子齿表面的网格节点相连接，以将定子绕组的质量平均分配到每个定子齿上；参数确定单元，用于确定所述定子铁芯的等效密度及定子绕组等效质量点的质量。相应地，本发明还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行上述等效仿真方法的步骤。相应地，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现等效仿真方法的步骤。实施本发明，具有如下有益效果：本发明将定子绕组等效为质量点而不建立定子绕组的实体模型，并将定子绕组等效为质量点与定子铁芯三维有限元网格进行关联以构成定子三维有限元模型，因此，本发明只需要进行一次模态试验，大大地减少了仿真方法的工作量。同时，本发明用具有与定子模态仿真结果相关的各向异性材料参数来等效考虑定子铁芯和定子绕组的层叠结构，材料参数少，精度易于控制，仿真精度高。附图说明图1是本发明考虑电机定子层叠结构的等效仿真方法的第一实施例流程图；图2是本发明考虑电机定子层叠结构的等效仿真方法的第二实施例流程图；图3是本发明中构建定子铁芯三维有限元网格的示意图；图4是本发明考虑电机定子层叠结构的等效仿真系统的结构示意图；图5是本发明中模型构建模块的结构示意图。具体实施方式为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。参见图1，图1显示了本发明考虑电机定子层叠结构的等效仿真方法的第一实施例流程图，包括：s101，建立定子铁芯三维有限元网格及定子绕组等效质量点，并关联所述定子铁芯三维有限元网格及定子绕组等效质量点以构成定子三维有限元模型。现有的仿真方法需要对定子铁芯和定子绕组分别建立实体模型，因此需要进行两次模态试验；与现有技术不同的是，本发明中将定子绕组等效为质量点而不需要建立绕组的实体模型，因此只需进行一次模态试验。综上，本发明提供的方法工作量少，需要确定的参数少，精度易于控制，仿真精度高。s102，根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数，并采用所述定子三维有限元模型分析所述材料参数的取值对定子模态频率的影响规律。具体地，所述根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数并采用定子三维有限元模型分析材料参数的取值对定子模态频率的影响规律的步骤包括：(1)根据定子层叠结构的特点确定各向异性材料的独立材料参数。(2)给所述独立材料参数赋予初始值。(3)改变任一独立材料参数的取值并保持其他独立材料参数固定不变，采用所述定子三维有限元模型仿真分析定子模态频率随所述独立材料参数的变化规律。(4)采用同样的方法分析定子模态频率随其他独立材料参数的变化规律。需要说明的是，本发明的思路是用虚拟的各向异性材料来等效考虑定子铁芯和定子绕组的层叠结构，因此对电机定子进行模态仿真首先需要确定各向异性材料的参数，而要确定这些参数的取值则需要先研究各个参数对定子模态的影响规律。由于各向异性材料的三个弹性常数(杨氏模量、剪切模量和泊松比)共有9个材料参数，而这些材料参数并不全部都是相互独立，因此，需要根据定子的层叠结构特点确定这9个材料参数中的独立材料参数(即，与定子模态仿真结果相关的材料参数)。具体地，各向异性材料的三个弹性常数(杨氏模量e、剪切模量g和泊松比v)共有9个材料参数，分别是三个方向的杨氏模量(ex、ey、ez)，三个平面的剪切模量(gxy、gyz、gzx)及三个泊松比(vxy、vyz、vzx)。如图3所示，由于定子铁芯由硅钢片沿轴向(z向)堆叠而成，所以硅钢片所在的平面(xy平面)是一个各向同性的平面，因此有ex＝ey，gyz＝gzx。同时，由于金属材料的泊松比的变化范围很小，而且泊松比对模态频率的影响不大，所以直接假设vxy＝vyz＝vzx＝0.3。因此，定子铁芯等效的各向异性材料的9个材料参数中只有4个独立材料参数，分别是ex(ey)、ez、gxy和gyz(gzx)。另外，分析各独立材料参数(ex(ey)、ez、gxy和gyz(gzx))对定子模态频率的影响规律，具体就是改变其中某一个独立材料参数的取值，同时保持其他独立材料参数固定不变，然后采用所述定子三维有限元模型对定子进行模态仿真，观察定子模态频率随每一个独立材料参数的变化规律。例如，针对某新能源车用驱动电机的48槽定子。令各向异性材料的三个独立材料参数为ez＝210gpa，gxy＝gyz(gzx)＝80.77gpa(硅钢片的材料参数为ex＝ey＝ez＝210gpa，gxy＝gyz＝gzx＝80.77gpa)，然后将ex(ey)的取值从210gpa分十步均匀减小到21gpa，通过模态仿真即可得到每一步的定子模态频率和定子模态振型。同时，采用同样的方法研究其他三个独立材料参数对定子模态频率的影响规律。结果表明：x(y)方向的杨氏模量ex(ey)以及xy平面的剪切模量gxy对定子径向模态的模态频率影响很大，定子径向模态的模态频率随ex(ey)和gxy的减小而降低，而z方向的杨氏模量ez以及yz(zx)平面的剪切模量gyz(gzx)对定子径向模态的模态频率几乎没有影响；x(y)方向的杨氏模量ex(ey)，xy平面的剪切模量gxy以及yz(zx)平面的剪切模量gyz(gzx)对定子反对称模态的模态频率影响很大，而z方向的杨氏模量ez对定子反对称模态的模态频率几乎没有影响。因此，杨氏模量ex(ey)、剪切模量gxy及剪切模量gyz(gzx)对定子模态频率的影响较大。s103，根据所述影响规律从所述材料参数中选取目标材料参数，并对所述目标材料参数的取值进行迭代修正，直到所述三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代。为了快速确定定子的等效各向异性材料参数，使定子模态仿真结果接近试验结果，需要重点关注那些对定子模态频率影响大的材料参数，因此，一般情况下，将对定子模态频率影响较大的材料参数作为目标材料参数。例如，针对某新能源车用驱动电机的48槽定子，杨氏模量ex(ey)、剪切模量gxy及剪切模量gyz(gzx)对定子模态频率的影响较大，可作为目标材料参数。具体地，所述对目标材料参数的取值进行迭代修正直到三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代的步骤包括：(1)给所述目标材料参数赋值，并采用所述三维有限元模型以对定子进行模态仿真，生成定子模态仿真结果。(2)将所述定子模态仿真结果与试验结果进行对比，若定子模态仿真结果与试验结果之间的偏差超出目标范围，则返回步骤(1)重新对所述目标材料参数进行赋值，直到所述定子模态仿真结果与试验结果与试验结果之间的偏差在目标范围内时停止迭代。需要说明的是，所述定子模态仿真结果包括定子模态频率及对应的定子模态振型，同时，所述试验结果包括试验模态频率及试验模态振型。比对时，根据定子模态振型查找匹配的试验模态振型，再将定子模态振型所对应的定子模态频率与试验模态振型所对应的试验模态频率进行比对，若定子模态频率与试验模态频率之间的偏差超出目标范围，则重新对目标材料参数进行赋值。其中，所示目标范围优选为5％，但不以此为限制，可根据实际情况进行预先设置。下表为某新能源车用驱动电机的目标材料参数进行迭代修正后，得出的定子铁芯各向异性等效材料参数：s104，将修正后的目标材料参数作为定子的等效材料参数并赋予所述定子三维有限元模型以进行定子模态仿真。将步骤s103中反复迭代得到目标材料参数作为定子最终的等效材料参数，再将所述等效材料参数赋予步骤s101中建立的定子三维有限元模型进行定子模态仿真，最后仿真得到的模态仿真结果能较准确地反映实际定子的振动特性。下表为某新能源车用驱动电机的模态仿真结果与试验模态结果的比对情况：模态仿真结果试验结果相对误差二阶模态645hz634hz1.74％三阶模态1781hz1765hz0.91％四阶模态3179hz3327hz4.45％五阶模态4617hz4797hz3.75％呼吸模态5635hz5546hz1.60％因此，从模态仿真结果与试验模态结果的对比可以看出采用本发明提供的等效仿真方法得到的模态仿真结果能较准确地反映实际定子的振动特性。由上可知，本发明通过对电机定子进行准确的模态分析，工作量少，仿真精度高，仿真结果可用于后续的振动噪声分析和电机结构优化，为电机的设计和优化奠定了基础。参见图2，图2显示了本发明考虑电机定子层叠结构的等效仿真方法的第二实施例流程图，包括：s201，构建定子铁芯三维有限元网格。需要说明的是，实际定子铁芯的每个定子齿都是相同的，所以在有限元建模和网格划分过程中需要保证定子结构的对称性，减少由结构不对称带来的误差。如图3所示，所述构建定子铁芯三维有限元网格的步骤包括：(1)构建定子齿二维平面模型。具体地，可根据电机定子的几何参数绘制定子齿二维平面模型，或者从预设定子三维有限元模型中切割出定子齿二维平面模型。其中，所述定子齿二维平面模型用于表示一个定子齿的二维平面模型。(2)对所述定子齿二维平面模型进行网格划分，生成定子齿二维平面网格。(3)对所述定子齿二维平面网格进行周向旋转阵列，生成定子截面二维平面网格。需要说明的是，对一个定子齿二维平面网格进行周向旋转阵列后，可得到整个定子截面的二维平面网格，从而保证每个定子齿的网格一致性。(4)将所述定子截面二维平面网格沿定子轴向拉伸，生成定子铁芯三维有限元网格。优选地，拉伸长度等于实际定子的轴向长度。s202，在定子铁芯的几何中心位置构建定子绕组等效质量点。考虑到实际的定子绕组是缠绕在定子齿上，其质量主要分布在定子齿上，且定子绕组刚度相对较小，对电机定子模态的影响相对较小，此外也为了减少控制变量的个数，本发明不对定子绕组建立实体模型，而是用质量点来代替。本发明中，在定子铁芯的几何中心位置处(即，定子轴向中间位置横截面的圆心位置处)建立定子绕组等效质量点，用于存放定子绕组的质量，而不对定子绕组建立实体模型。s203，将所述定子绕组等效质量点与所述铁芯三维有限元网格中每个定子齿表面的网格节点相连接，以将定子绕组的质量平均分配到每个定子齿上。根据实际定子绕组是缠绕在定子齿上特点，本发明将定子绕组等效质量点与定子齿表面的网格节点连接起来，从而将定子绕组的质量平均分配到每个定子齿上。因此，与现有技术不同的是，本发明不需要单独建立定子绕组的实体模型，定子绕组的质量通过等效质量点来考虑并将其质量平均分配到每个定子齿上，而定子绕组的刚度则等效到定子铁芯中，通过给定子铁芯赋予各向异性材料来等效考虑实际定子总成(定子铁芯和定子绕组)的刚度。s204，确定所述定子铁芯的等效密度及定子绕组等效质量点的质量。实际的定子铁芯是由很薄的硅钢片沿轴向堆叠而成，但是在建模过程中无法为每一片硅钢片单独建模，只能将整个定子铁芯当成一个实体来建模，所以定子有限元模型不能直接赋予硅钢片的密度，而是要赋予一个等效密度。具体地，可根据公式ρ＝m/v，计算定子铁芯的等效密度，其中，ρ为定子铁芯的等效密度，m为实际定子铁芯的质量，v为定子铁芯三维有限元网格的体积。同时，可将实际定子绕组的质量作为所述定子绕组等效质量点的质量。例如，针对某新能源车用驱动电机的48槽定子，其定子绕组质量为4.1kg，定子铁芯等效密度为7458kg/m3，明显小于硅钢片的密度(7850kg/m3)。需要说明的是，通过步骤s201-204后，即可形成完整的定子三维有限元模型，以便进行模态仿真处理。s205，根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数，并采用所述定子三维有限元模型分析所述材料参数的取值对定子模态频率的影响规律。具体地，所述根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数并采用定子三维有限元模型分析材料参数的取值对定子模态频率的影响规律的步骤包括：(1)根据定子层叠结构的特点确定各向异性材料的独立材料参数；(2)给所述独立材料参数赋予初始值；(3)改变任一独立材料参数的取值并保持其他独立材料参数固定不变，采用所述定子三维有限元模型仿真分析定子模态频率随所述独立材料参数的变化规律；(4)采用同样的方法分析定子模态频率随其他独立材料参数的变化规律。s206，根据所述影响规律从所述材料参数中选取目标材料参数，并对所述目标材料参数的取值进行迭代修正，直到所述三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代。具体地，所述对目标材料参数的取值进行迭代修正直到三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代的步骤包括：(1)给所述目标材料参数赋值，并采用所述三维有限元模型以对定子进行模态仿真，生成定子模态仿真结果。(2)将所述定子模态仿真结果与试验结果进行对比，若定子模态仿真结果与试验结果之间的偏差超出目标范围，则返回步骤(1)重新对所述目标材料参数进行赋值，直到所述定子模态仿真结果与试验结果与试验结果之间的偏差在目标范围内时停止迭代。s207，将修正后的目标材料参数作为定子的等效材料参数并赋予所述定子三维有限元模型以进行定子模态仿真。由上可知，本发明将定子绕组等效为质量点而不建立定子绕组的实体模型，然后根据实际定子绕组缠绕在定子齿上的特点，将定子绕组等效质量点的质量平均分配到每个定子齿上，接着再对定子铁芯进行等效处理，用具有与定子模态仿真结果相关的各向异性材料参数来等效考虑定子铁芯和定子绕组的层叠结构，仿真过程只需要进行一次模态试验。相比之下，现有的仿真方法需要分别对定子铁芯和定子绕组建立实体模型，然后分别确定定子铁芯和定子绕组的各向异性材料参数，且要进行两次模态试验，也就是说现有方仿真法需要确定的各向异性材料参数个数和试验次数比本发明提供的仿真方法多一倍，因此，本发明提供的方法工作量少，需要确定的参数少，精度易于控制，仿真精度高。参见图4，图4显示了本发明考虑电机定子层叠结构的等效仿真系统100的具体结构，其包括：模型构建模块1，用于建立定子铁芯三维有限元网格及定子绕组等效质量点，并关联所述定子铁芯三维有限元网格及定子绕组等效质量点以构成定子三维有限元模型。模型构建模块1仅需对定子铁芯建立实体模型，并将定子绕组等效为质量点而不需要建立绕组的实体模型，可大大提高仿真效率。规律分析模块2，用于根据各向异性材料的特点确定与定子模态仿真结果相关的材料参数，并采用所述定子三维有限元模型分析所述材料参数的取值对定子模态频率的影响规律。具体地，规律分析模块2根据定子层叠结构的特点确定定子等效各向异性材料的9个材料参数(三个方向的杨氏模量、三个平面的剪切模量和三个泊松比)中的独立材料参数；给所述独立材料参数赋予初始值；改变其中某一个独立材料参数的取值并保持其他独立材料参数固定不变，采用所述定子三维有限元模型仿真分析定子模态频率随该独立材料参数的变化规律；最后，采用同样的方法研究定子模态频率随其他独立材料参数的变化规律。迭代修正模块3，用于根据所述影响规律从所述材料参数中选取目标材料参数，并对所述目标材料参数的取值进行迭代修正，直到所述三维有限元模型仿真得到的定子模态仿真结果与定子模态试验结果之间的偏差在目标范围内时停止迭代。具体地，迭代修正模块3根据所述影响规律从所述材料参数中提取目标材料参数后，分别给所述目标材料参数赋值，然后采用所述三维有限元模型以对定子进行模态仿真，生成定子模态仿真结果，再将所述定子模态仿真结果与试验结果进行对比，若定子模态仿真结果与试验结果之间的偏差超出目标范围，则重新对所述目标材料参数进行赋值，直到所述定子模态仿真结果与试验结果与试验结果之间的偏差在目标范围内时停止迭代，从而实现了目标材料参数的取值的最优化。模态仿真模块4，用于将修正后的目标材料参数作为定子的等效材料参数并赋予所述定子三维有限元模型以进行定子模态仿真。因此，本发明通过对电机定子进行准确的模态分析，工作量少，仿真精度高，仿真结果可用于后续的振动噪声分析和电机结构优化，为电机的设计和优化奠定了基础。如图5所示，所述模型构建模块1包括：网格构建单元11，用于构建定子铁芯三维有限元网格。具体地，网格构建单元11先构建定子齿二维平面模型，再对所述定子齿二维平面模型进行网格划分以生成定子齿二维平面网格，然后对所述定子齿二维平面网格进行周向旋转阵列以生成定子截面二维平面网格，最后将所述定子截面二维平面网格沿定子轴向拉伸，生成定子铁芯三维有限元网格。质量点构建单元12，用于在定子铁芯的几何中心位置构建定子绕组等效质量点。关联单元13，用于将所述定子绕组等效质量点与所述铁芯三维有限元网格中每个定子齿表面的网格节点相连接，以将定子绕组的质量平均分配到每个定子齿上。参数确定单元14，用于确定所述定子铁芯的等效密度及定子绕组等效质量点的质量。具体地，可根据公式ρ＝m/v，计算定子铁芯的等效密度，其中，ρ为定子铁芯的等效密度，m为实际定子铁芯的质量，v为定子铁芯三维有限元网格的体积；同时，可将实际定子绕组的质量作为所述定子绕组等效质量点的质量。因此，与现有技术不同的是，本发明不需要单独建立定子绕组的实体模型，定子绕组的质量通过等效质量点来考虑并将其质量平均分配到每个定子齿上，而定子绕组的刚度则等效到定子铁芯中，通过给定子铁芯赋予各向异性材料来等效考虑实际定子总成(定子铁芯和定子绕组)的刚度。相应地，本发明还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述等效仿真方法的步骤。同时，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述等效仿真方法的步骤。以上所述是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。当前第1页12