CN105787199A - 一种叶轮机转子叶片从热态到冷态的叶型转换方法 - Google Patents
一种叶轮机转子叶片从热态到冷态的叶型转换方法 Download PDFInfo
- Publication number
- CN105787199A CN105787199A CN201610153913.6A CN201610153913A CN105787199A CN 105787199 A CN105787199 A CN 105787199A CN 201610153913 A CN201610153913 A CN 201610153913A CN 105787199 A CN105787199 A CN 105787199A
- Authority
- CN
- China
- Prior art keywords
- blade
- hot
- blade profile
- profile
- cold conditions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明涉及一种叶轮机转子的叶型转换方法,特别涉及一种叶轮机转子叶片从热态到冷态的叶型转换方法。所述的方法在已知叶轮机转子叶片热态叶型的前提下,基于有限元逆向迭代的叶型转换方法,获取叶片的冷态叶型。所述的方法提高了计算精度,并且节省了人力,提高了效率。
Description
技术领域
本发明涉及一种叶轮机转子的叶型转换方法,特别涉及一种叶轮机转子叶片从热态到冷态的叶型转换方法。
背景技术
当前叶轮机转子叶片设计中,经常需要在热态和冷态两种叶型之间进行转换。所谓冷态叶型,即工程图上给出的机械加工叶型,提供给叶片制造加工方和结构强度振动分析人员;而热态叶型指的是气动设计人员在压气机设计阶段给出的工作状态叶型,即高速旋转的柔性叶片在离心载荷、表面气动载荷和温度载荷的多重作用下形成的变形后的叶型。由于叶轮机转子叶片的冷态和热态叶型之间差别较大,存在弹性恢复角、径向拉伸、轴向偏移和周向挠度,因此作为叶轮机完整设计链中必不可少的一环,热态和冷态叶型之间的转换精度尤其需要得到设计者的重视。当前主流方法是采用近似线性化方法处理叶型从热态到冷态的转换问题,这种方法将叶片冷热态之间的变形视为完全线性过程,通过对热态叶型实施强度分析获得近似的变形,再将近似得到的变形从热态叶型中减去来获得冷态叶型。现有的方法至少存在以下三个缺点:
首先,现有方法无法考虑叶片的非线性作用,这对于过去刚性较强的压气机叶片是没有问题的,但对于新型的高负荷高性能叶片而言,其气动载荷更高而刚性更差,不能忽视由几何非线性导致的叶片大变形。因此,压气机尤其是低压压气机叶片的变形应视为非线性,不适宜进行简单的线性化处理。
其次,现有方法通过试探获得了“伪”冷态叶片,使用“伪”冷态叶片变形而获得的最终热态叶片方案与初始的气动方案存在一定差别,导致最终方案无法完全忠实于设计者的原始设计思路。
最后,现有方法的操作流程繁琐,需要多次在气动设计和结构设计间反复修正,必须进行大量的人工干预,实施效率较低。
发明内容
本发明解决的技术问题为:克服常规的从热态到冷态的叶型转换方法的不足,提供一种基于有限元逆向迭代的叶型转换方法,以准确获得叶片的冷态叶型,使工厂根据本方法提供的冷态叶型加工得到的叶片在投入实际使用时,其工作状态的热态叶型能够与气动设计人员给定的原始气动热态叶型一致,从而保证实际叶片能够忠实于设计者的意图。
本发明的技术方案为:所述的方法在已知叶轮机转子叶片热态叶型的前提下,基于有限元逆向迭代的叶型转换方法,获取叶片的冷态叶型。
作为本技术方案的一种改进,所述的方法包括如下步骤:
步骤一,获取叶片热态叶型坐标,并计算叶片的容差ε1、ε2和ε3;
步骤二,由流体计算软件对热态叶型进行三维计算流体动力学分析,获得叶片表面静压;
步骤三,将叶片热态叶型的流体网格和叶片表面静压转化为叶片热态叶型的结构网格上的气动力;
步骤四,由已知的热态叶型坐标和气动力,将叶片热态叶型结构网格转化为叶片冷态叶型的结构网格;
步骤五:将叶片由热态叶型坐标到模拟热态叶型坐标的有限元网格的变形,通过三维贴体插值得到模拟热态叶型的流体网格;
步骤六:由流体计算软件对所获得的模拟热态叶型的流体网格进行三维计算流体动力学分析;
步骤七:如果步骤六和步骤二两次流体动力学分析获得的压气机特性间的差异不满足容差,则用模拟热态叶型坐标代替热态叶型坐标,回到步骤三,否则计算结束,此时得到的叶片冷态叶型的结构网格即为所需要求的叶轮机转子叶片的冷态叶型的坐标。
作为本技术方案的一种改进,通过三维贴体插值,将叶片热态叶型的流体网格和叶片表面静压转化为叶片热态叶型的结构网格上的气动力。
作为本技术方案的一种改进,通过迭代法将叶片热态叶型结构网格转化为叶片冷态叶型的结构网格。
本发明的有益效果为:(1)常规方法都采用直接近似法处理叶型从热态到冷态的转换问题,对热态叶型实施强度分析获得变形,再将其从热态叶型中减去来获得近似的冷态叶型,这大幅降低了计算精度。本发明通过反复逆向迭代的方法随时修改结构刚度矩阵和气动载荷,克服了这一缺陷。
(2)常规方法采用了完全线性化假设,这对于弯扭掠复合的先进压气机叶片是不适合的,针对此类型叶片计算出的弹性恢复角的误差极大,这将导致所生产出的叶片在工作状态的热态叶型严重偏离原始气动设计。本发明在计算中引入了几何非线性效应,采用了非线性方法描述叶型从热态到冷态的转换过程,避免了大量的线性化假设对计算结果精度的影响。
(3)本发明采用了三维贴体形函数插值方法来传递气动载荷和结构变形,相对于常规的近似距离权函数插值法,可以提高插值精度和插值效率。
(4)常规方法的操作流程繁琐,需要多次在气动设计和结构设计间反复修正,必须进行大量的人工干预,实施效率较低,而本发明采用Fortran程序对整个转换过程进行了自动化操作,大大提高了计算效率,更加简便易行。
附图说明
图1是本发明的实现流程图;
图2是本发明的三维贴体插值流程图;
图3是本发明的三维贴体插值示意图;
图4是本发明的热态到冷态叶型迭代流程图;
图5是本发明实施例中的流体表面网格;
图6是本发明实施例中的热态到冷态迭代历程图;
图7是本发明实施例中的叶片变形曲线。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是叶轮机转子叶片从热态到冷态的叶型转换方法,方法的实现流程如图1所示。所述的方法在已知叶轮机转子叶片热态叶型的前提下,基于有限元逆向迭代的叶型转换方法,获取叶片的冷态叶型。
步骤一:给定叶片热态叶型(其坐标为XH)、计算容差ε1、ε2和ε3。
步骤二:由流体计算软件对热态叶型XH进行三维计算流体动力学分析,获得叶片表面静压。
步骤三:将叶片热态叶型的流体网格和叶片表面静压通过三维贴体插值转化为叶片热态叶型的结构网格上的气动力Fp,实施流程如图2所示,其具体步骤如下:
步骤3.1:遍历所有耦合界面上的流体节点,找到距离结构节点s最近的流体节点f。
步骤3.2:找到所有包含流体节点f的流体表面网格单元c1,c2,…,cm,并将这些流体表面网格单元沿垂直网格面方向进行拉伸,构造成为拥有8个节点的六面体临时单元h1,h2,…,hm。
步骤3.3:对每个临时单元hi,求出s相对于单元hi的局部坐标(ri,si,ti),并求出s与单元hi中心的相对距离
步骤3.4:找到与s相对距离最小的单元ht,即dt=min(d1,d2,…,dm),单元ht即为结构节点s所对应的临时单元。
步骤3.5:假设结构节点s在流体表面网格单元上的投影点为s',如图3所示。投影点s'在整体坐标系下的坐标为(xs',ys',zs'),忽略流体节点与其投影点之间的距离,近似有xs'≈xs,ys'≈ys,zs'≈zs,基于有限元形函数理论,投影点s'相对于临时单元ht的局部坐标(r,s,t)可以通过形函数插值得到
其中(xi,yi,zi),1≤i≤8为临时单元ht中8个节点在整体坐标系下的坐标,Ni为投影点处的形函数
这个方程为三元三次方程,直接求解比较困难,可以通过近似法求解,如拟牛顿法和梯度法等。
步骤3.6:假设结构节点s上的流体静压为ps,投影点s'的位移为ps',近似有ps≈ps',投影点s'的静压可以通过形函数插值得到,即式中pi是临时单元ht的8个节点i上的静压,由于构造临时单元ht的8个节点中的4个临时节点上的形函数N=0,可得由此就通过流体表面网格上各节点的静压插值获得了结构节点s的静压ps,并将其转换为结构节点上的气动力Fp。
步骤四:由已知的热态叶型坐标XH和气动力Fp,通过迭代将叶片热态叶型结构网格转化为叶片冷态叶型的结构网格,假设所计算得到的冷态叶型的坐标为XC,实施流程如图4所示,具体步骤如下:
步骤4.1:设叶片的转速为Ω,则引入矩阵
步骤4.2:令i=0,并假设当前冷态叶型坐标
步骤4.3:根据当前冷态叶型坐标求得结构的小位移刚度矩阵K0、离心刚度矩阵KC和所受的离心力FC,分别为
其中所有的积分都是在结构单元e上进行体积积分,式中B0代表小位移刚度矩阵,D为弹性矩阵,ρ为密度,N为形函数矩阵,[x0y0z0]T为初始坐标。
步骤4.4:令n=0,并假设初始位移
步骤4.5:根据an求得应力刚化矩阵Kσ、大变形刚度矩阵KL和反力Fnr,分别为:
其中:
σ=[σxσyσzτxyτyzτxz]T
BL=AG
步骤4.6:根据(K0-KC+Kσ+KL)dan+1=FC+Fp-Fnr求出增量位移dan+1。
步骤4.7:计算新位移an+1=an+dan+1。
步骤4.8:如果满足收敛条件||dan+1||2<ε1||an+1||2,则获得总增量刚度矩阵和位移a(i)=an+1,否则令n=n+1,前往步骤4.5。
步骤4.9:计算模拟热态叶型坐标和下一步新的冷态叶型坐标
步骤4.10:如果满足收敛条件则获得冷态叶型坐标和模拟热态叶型坐标否则令i=i+1,前往步骤4.3。
步骤五:将叶片由热态叶型坐标XH到模拟热态叶型坐标X'H的有限元网格的变形通过三维贴体插值得到模拟热态叶型的流体网格,其方法与步骤三相似。
步骤六:由流体计算软件对所获得的模拟热态叶型X'H的流体网格进行三维计算流体动力学分析。
步骤七:如果步骤六和步骤二两次流体动力学分析获得的压气机特性间的差异不满足容差ε3,则用X'H代替XH,回到步骤三,否则计算结束,此时得到的XC即为所需要求的叶轮机转子叶片的冷态叶型的坐标。
现采用某航空发动机压气机转子叶片为实例,对本发明作进一步的详细说明。本例采用自编程序的形式实现发明内容,自编程序的开发平台是工程计算语言Fortran90,编译平台是IntelVisualFortran。例中流体动力学计算采用了商用软件NUMECA。
步骤1:由气动设计人员给出该叶片热态叶型的叶型坐标表。设定计算容差ε1=1e-6、ε2=1e-6和ε3=1e-3。同时给定叶片材料参数和设计转速。采用自编程序的trans_NUMECA函数将叶型坐标表转化为geomTurbo几何文件。
步骤2:由商用软件NUMECA读入步骤1给出的geomTurbo几何文件,采用NUMECA的AutoGrid5模块进行网格划分得到流体网格,然后采用NUMECA的Fine模块进行参数设定和流体动力学计算,得到设计转速下的叶片表面静压分布和压气机特性(流量-总压比关系曲线、流量-效率关系曲线等)。由NUMECA的CFView模块给出叶片流体表面网格的网格坐标和表面静压,输出至文件surfpre.dat中。叶片流体表面网格如图5所示。
步骤3:调用自编程序的trans_fmesh_smesh函数,读入surfpre.dat文件,通过插值获得叶片的8节点六面体结构网格,存入数组coord(num_node,3)和elem(num_elem,8)中,其中num_node为结构网格的节点数,每个节点存入x、y和z三个方向的坐标值;num_elem为结构网格的单元数,每个单元存入8个节点的编号。然后采用三维贴体形函数插值将流体网格的表面静压插值到结构网格的节点上,形成结构节点力,存入列阵gfp中。
步骤4:采用自编程序的hot_cold函数,经过迭代获得了冷态叶型和模拟热态叶型,本例中经过了8次迭代达到收敛标准,收敛历程如图6所示。
步骤5:调用自编程序的trans_smesh_fmesh函数,采用三维贴体形函数插值,将模拟热态叶型的结构网格转化为气动造型输入的geomTurbo几何文件。
步骤6:由商用软件NUMECA读入步骤5给出的geomTurbo几何文件,采用NUMECA的AutoGrid5模块进行网格划分得到流体网格,然后采用NUMECA的Fine模块进行参数设定和流体动力学计算,得到设计转速下的叶片表面静压分布和压气机特性。由NUMECA的CFView模块给出叶片流体表面网格的网格坐标和表面静压,输出至文件surfpre_new.dat中。
步骤7:将步骤6计算获得的压气机特性与步骤2的计算结果进行比较。如果不满足容差,则用步骤4获得的模拟热态叶型代替给定热态叶型,返回步骤3继续进行计算,直至满足收敛条件为止。本例中经过2次迭代达到收敛标准。此时所输出的gongyi.dat文件即为所需要输出的冷态叶型的叶型坐标表,经过一定数据标准化处理后则成为提供给工厂加工使用的工艺文件。同时还输出了hot_cold.dat文件,在文件中给出了叶片从热态到冷态的几何变形,包括叶片的径向伸长量、周向挠度变化、轴向位置变化和弹性恢复角。该文件可以输出给Excel或Origin进行后处理,获得叶片的变形曲线,如图7所示。图中横坐标为变形量,纵坐标为叶高。
Claims (4)
1.一种叶轮机转子叶片从热态到冷态的叶型转换方法,其特征为:所述的方法在已知叶轮机转子叶片热态叶型的前提下,基于有限元逆向迭代的叶型转换方法,获取叶片的冷态叶型。
2.根据权利要求1所述的一种叶轮机转子叶片从热态到冷态的叶型转换方法,其特征为,所述的方法包括如下步骤:
步骤一,获取叶片热态叶型坐标,并计算叶片的容差ε1、ε2和ε3;
步骤二,由流体计算软件对热态叶型进行三维计算流体动力学分析,获得叶片表面静压;
步骤三,将叶片热态叶型的流体网格和叶片表面静压转化为叶片热态叶型的结构网格上的气动力;
步骤四,由已知的热态叶型坐标和气动力,将叶片热态叶型结构网格转化为叶片冷态叶型的结构网格;
步骤五:将叶片由热态叶型坐标到模拟热态叶型坐标的有限元网格的变形,通过三维贴体插值得到模拟热态叶型的流体网格;
步骤六:由流体计算软件对所获得的模拟热态叶型的流体网格进行三维计算流体动力学分析;
步骤七:如果步骤六和步骤二两次流体动力学分析获得的压气机特性间的差异不满足容差,则用模拟热态叶型坐标代替热态叶型坐标,回到步骤三,否则计算结束,此时得到的叶片冷态叶型的结构网格即为所需要求的叶轮机转子叶片的冷态叶型的坐标。
3.根据权利要求2所述的一种叶轮机转子叶片从热态到冷态的叶型转换方法,其特征为:通过三维贴体插值,将叶片热态叶型的流体网格和叶片表面静压转化为叶片热态叶型的结构网格上的气动力。
4.根据权利要求2所述的一种叶轮机转子叶片从热态到冷态的叶型转换方法,其特征为:通过迭代法将叶片热态叶型结构网格转化为叶片冷态叶型的结构网格。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510387123X | 2015-07-06 | ||
CN201510387123 | 2015-07-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105787199A true CN105787199A (zh) | 2016-07-20 |
CN105787199B CN105787199B (zh) | 2019-05-07 |
Family
ID=56393965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610153913.6A Active CN105787199B (zh) | 2015-07-06 | 2016-03-17 | 一种叶轮机转子叶片从热态到冷态的叶型转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105787199B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339557A (zh) * | 2016-08-29 | 2017-01-18 | 张明 | 旋转机械叶片反问题设计的计算方法 |
CN106709179A (zh) * | 2016-12-21 | 2017-05-24 | 中国燃气涡轮研究院 | 一种适用于有限元自由网格的叶片气动力分布数据处理方法 |
CN109556765A (zh) * | 2018-12-10 | 2019-04-02 | 中国航发四川燃气涡轮研究院 | 一种叶片非接触振动应力测量值换算方法 |
CN109595040A (zh) * | 2018-12-21 | 2019-04-09 | 杭州汽轮动力集团有限公司 | 一种燃气轮机叶片扭曲规律设计方法 |
CN109684674A (zh) * | 2018-12-04 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种舱门气动载荷处理方法 |
CN111368373A (zh) * | 2020-03-01 | 2020-07-03 | 南京航空航天大学 | 一种冷态叶型与热态叶型的迭代分析方法 |
CN112131681A (zh) * | 2020-09-15 | 2020-12-25 | 中国航发沈阳发动机研究所 | 一种冷热态结构变换方法 |
CN112597613A (zh) * | 2021-03-08 | 2021-04-02 | 中国航发上海商用航空发动机制造有限责任公司 | 确定压气机叶片所受气动合力及其作用点的方法 |
CN113094963A (zh) * | 2021-06-09 | 2021-07-09 | 中国航发上海商用航空发动机制造有限责任公司 | 一种复合材料叶片冷态叶型的计算方法 |
CN114818436A (zh) * | 2021-12-14 | 2022-07-29 | 泸州懋威科技有限公司 | 一种离心叶轮和涡轮叶片叶型冷热态快速转换方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN102799730A (zh) * | 2012-07-13 | 2012-11-28 | 北京航空航天大学 | 一种燃气轮机风扇叶片反扭过程的预估方法 |
CN105382489A (zh) * | 2015-11-19 | 2016-03-09 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种钛合金板料成型环形件热定型装置及工艺 |
-
2016
- 2016-03-17 CN CN201610153913.6A patent/CN105787199B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN102799730A (zh) * | 2012-07-13 | 2012-11-28 | 北京航空航天大学 | 一种燃气轮机风扇叶片反扭过程的预估方法 |
CN105382489A (zh) * | 2015-11-19 | 2016-03-09 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种钛合金板料成型环形件热定型装置及工艺 |
Non-Patent Citations (4)
Title |
---|
孙树立等: "基于节点坐标修正和网格修匀的涡轮叶片", 《机械工程学报》 * |
牛冬生: "涡轮叶尖间隙数值分析和冷热态尺寸换算方法研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
王延荣等: "基于热态叶型的小轮毂比弯掠叶片", 《航空动力学报》 * |
郑赟等: "跨声速风扇叶片的反扭设计研究", 《机械工程学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339557A (zh) * | 2016-08-29 | 2017-01-18 | 张明 | 旋转机械叶片反问题设计的计算方法 |
CN106709179A (zh) * | 2016-12-21 | 2017-05-24 | 中国燃气涡轮研究院 | 一种适用于有限元自由网格的叶片气动力分布数据处理方法 |
CN106709179B (zh) * | 2016-12-21 | 2020-04-07 | 中国燃气涡轮研究院 | 一种叶片气动力分布数据处理方法 |
CN109684674B (zh) * | 2018-12-04 | 2023-05-05 | 中国航空工业集团公司西安飞机设计研究所 | 一种舱门气动载荷处理方法 |
CN109684674A (zh) * | 2018-12-04 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种舱门气动载荷处理方法 |
CN109556765A (zh) * | 2018-12-10 | 2019-04-02 | 中国航发四川燃气涡轮研究院 | 一种叶片非接触振动应力测量值换算方法 |
CN109595040A (zh) * | 2018-12-21 | 2019-04-09 | 杭州汽轮动力集团有限公司 | 一种燃气轮机叶片扭曲规律设计方法 |
CN111368373A (zh) * | 2020-03-01 | 2020-07-03 | 南京航空航天大学 | 一种冷态叶型与热态叶型的迭代分析方法 |
CN111368373B (zh) * | 2020-03-01 | 2022-06-03 | 南京航空航天大学 | 一种冷态叶型与热态叶型的迭代分析方法 |
CN112131681A (zh) * | 2020-09-15 | 2020-12-25 | 中国航发沈阳发动机研究所 | 一种冷热态结构变换方法 |
CN112131681B (zh) * | 2020-09-15 | 2024-02-02 | 中国航发沈阳发动机研究所 | 一种冷热态结构变换方法 |
CN112597613A (zh) * | 2021-03-08 | 2021-04-02 | 中国航发上海商用航空发动机制造有限责任公司 | 确定压气机叶片所受气动合力及其作用点的方法 |
CN113094963A (zh) * | 2021-06-09 | 2021-07-09 | 中国航发上海商用航空发动机制造有限责任公司 | 一种复合材料叶片冷态叶型的计算方法 |
CN114818436A (zh) * | 2021-12-14 | 2022-07-29 | 泸州懋威科技有限公司 | 一种离心叶轮和涡轮叶片叶型冷热态快速转换方法 |
CN114818436B (zh) * | 2021-12-14 | 2024-03-29 | 泸州懋威科技有限公司 | 一种离心叶轮和涡轮叶片叶型冷热态快速转换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105787199B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105787199A (zh) | 一种叶轮机转子叶片从热态到冷态的叶型转换方法 | |
CN111159636B (zh) | 一种柔性多体系统动力学半解析灵敏度分析方法 | |
CN102799730A (zh) | 一种燃气轮机风扇叶片反扭过程的预估方法 | |
CN109726437B (zh) | 一种舱门气动载荷等效节点力处理方法 | |
Gang et al. | Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
CN110717275A (zh) | 一种针对压水堆堆芯的三维中子通量数值模拟方法 | |
CN115630542B (zh) | 一种薄壁加筋结构的加筋布局优化方法 | |
CN112414668B (zh) | 一种风洞试验数据静气弹修正方法、装置、设备及介质 | |
CN105205293A (zh) | 用于获得飞机部件气动载荷的方法和系统 | |
CN111079327B (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
CN104573178B (zh) | 一种整体叶轮强度有限元计算方法 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Jonsson et al. | Development of flutter constraints for high-fidelity aerostructural optimization | |
CN117171894B (zh) | 一种考虑静稳定裕度约束的飞行器布局气动优化设计方法 | |
CN104504227A (zh) | 锻造液压机多维虚拟样机协同仿真方法及系统 | |
CN101334803B (zh) | 基于变形修正的动力减缩算法及超单元构造方法 | |
CN106126871A (zh) | 一种pscad电磁暂态仿真用的调速器模型建模方法 | |
CN106339557A (zh) | 旋转机械叶片反问题设计的计算方法 | |
CN100587696C (zh) | 一种参数化的静态超单元构造方法 | |
CN105550424A (zh) | 一种基于rbf网格变形插值序列的筛选方法 | |
CN113486512A (zh) | 一种功能梯度变厚度叶片模型的颤振分析方法 | |
Pałys et al. | Using metamodeling and fluid-structure interaction analysis in multi-objective optimization of a butterfly valve | |
CN109684674A (zh) | 一种舱门气动载荷处理方法 | |
CN116090274B (zh) | 基于量子计算的材料变形模拟方法、装置、终端及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |