CN109766613A - 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 - Google Patents
基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 Download PDFInfo
- Publication number
- CN109766613A CN109766613A CN201811637938.9A CN201811637938A CN109766613A CN 109766613 A CN109766613 A CN 109766613A CN 201811637938 A CN201811637938 A CN 201811637938A CN 109766613 A CN109766613 A CN 109766613A
- Authority
- CN
- China
- Prior art keywords
- jump
- probability
- thin
- sound pressure
- equation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000010276 construction Methods 0.000 title abstract 4
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 239000002131 composite material Substances 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims abstract description 17
- 230000014509 gene expression Effects 0.000 claims abstract description 16
- 239000010410 layer Substances 0.000 claims description 23
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 230000007704 transition Effects 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 10
- 238000005452 bending Methods 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 5
- 239000002356 single layer Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000009191 jumping Effects 0.000 abstract description 7
- 238000009795 derivation Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000004672 jump response Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开一种基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,包括:步骤1:获取待预测复合材料板铺层信息、材料和几何参数;步骤2:计算待预测复合材料的总体刚度矩阵;步骤3:热噪声载荷下薄壁结构动力学方程推导,获取无量纲单模态方程;步骤4:获取待预测复合材料板的刚度曲线,利用刚度曲线极值点获得跳变临界力fcr关于屈曲系数s的函数表达式;步骤5:列出给定声压级SPL的噪声载荷概率密度函数g(p);步骤6:比较噪声载荷的概率密度函数和跳变临界力对应的临界声压pcr,推导跳变临界声压级表达式,并给出对应的跳变概率。本发明对给定的薄壁结构以及温度噪声载荷,能基于理论计算立即给出跳变概率,具有高效快速的优点。
Description
技术领域
本发明属于动力学响应预测领域,特别涉及一种基于概率理论进行飞行器薄壁结构在热噪声载荷下的突弹跳变边界及跳变概率预测的方法。
背景技术
飞行器壁板结构在服役过程中将受到严酷的气动热和强噪声激励,其响应将表现出强非线性特征。高温使薄壁结构膨胀延伸导致热屈曲,在强噪声激励下,结构将在两个屈曲后平衡位置间跳动,即产生突弹跳变,频繁的应力反转容易造成结构的疲劳破坏。不同温度下,引起突弹跳变的临界载荷也各不相同,他们构成跳变边界。准确预测跳变边界有助于把握结构的响应规律,为结构的安全性评价和优化设计提供指导。
目前已有的研究突弹跳变的方法包括等效线性法[1]、Fokker-Planck-Kolmogorov方法[2](FPK法)、有限元法[3]、实验法等[3]。
文献[2]中的方法基于单模态方程,采用FPK方法获得位移响应的概率密度函数,根据概率密度函数随载荷的演化规律与跳变响应规律之间的关联,采用概率密度的极大值与极小值的比值构建判别参数Δ,当Δ靠近0时不发生跳变,当Δ靠近1时,响应为持续性跳变,在0和1之间,Δ越大,跳变频次越高。
等效线性法仅适用于稳态的无跳变和持续性跳变响应,无法预测非稳态的间歇性跳变。FPK法不能定量给出跳变概率,不能直观评估安全裕度。有限元法、实验法作为常规方法,需要通过大量的计算和实验,效率不高,且无法揭示跳变的内在机理。
参考文献:
[1]Ng,C.F."The analysis of non-linear dynamic behavior(includingsnap-through)of postbuckled plates by simple analytical solution,"NASATechnical Memorandum.1988.
[2]Liu,L.,Lv,B.Y.,and He,T.R."The stochastic dynamic snap-throughresponse of thermally buckled composite panels,"Composite Structures Vol.131,2015,pp.344-355.
[3]Goley,G.,Zappia,B.,Beberniss,T.,and Shukla,A."Effect of Loading onthe Snap-Through Response of a Post-Buckled Beam,"49th AIAA/ASME/ASCE/AHS/ASCStructures,Structural Dynamics,and Materials Conference,16th AIAA/ASME/AHSAdaptive Structures Conference,10th AIAA Non-Deterministic ApproachesConference,9th AIAA Gossamer Spacecraft Forum,4th AIAA MultidisciplinaryDesign Optimization Specialists Conference.American Institute of Aeronauticsand Astronautics,Schaumburg,IL,2008.
发明内容
本发明的目的在于提供一种基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,以解决上述技术问题。本发明从突弹跳变的内在机理出发,能覆盖包括间歇性跳变在内的所有响应类型,能定量给出跳变概率,便于评估安全裕度。对给定的薄壁结构以及温度噪声载荷,能基于理论计算立即给出跳变概率,具有高效快速的优点。
为实现以上的目的,本发明采取以下技术方案:
基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,包括以下步骤:
步骤1:获取待预测复合材料板铺层信息、材料和几何参数;
步骤2:计算待预测复合材料的总体刚度矩阵;
步骤3:热噪声载荷下薄壁结构动力学方程推导,获取无量纲单模态方程;
步骤4:获取不同温度下结构的刚度曲线,利用刚度曲线极值点获得跳变临界力fcr关于屈曲系数s的函数表达式;
步骤5:列出给定声压级SPL的噪声载荷概率密度函数g(p);
步骤6:比较噪声载荷的概率密度函数和跳变临界力对应的临界声压pcr,推导跳变临界声压级表达式,并给出对应的跳变概率。
进一步的,步骤1中:
获取复合材料板的铺层信息,包括铺设层数N、每层厚度h、每层铺设方向θ;
获取每层的材料参数,包括弹性模量E1,E2,泊松比ν12,ν21,剪切模量G12,G13,G23,密度ρ,热膨胀系数α1,α2;
获取几何参数,包括长Lx、宽Ly、厚h。
进一步的,步骤2中,首先确定每个单一铺层刚度矩阵,第k层的刚度矩阵为:
其中,
然后计算总体刚度矩阵
其中,Aij,Bij,Dij分别为拉伸、拉弯耦合、弯曲刚度矩阵元素,z为厚度方向坐标,zk+1和zk分别为第k层上下表面坐标,T(k)为第k层坐标转换矩阵。
θ(k)为第k层铺设角。
进一步的,步骤3具体包括:
采用假设模态法,将响应表示成双重傅里叶级数并截取第一项,应用伽辽金法将板的偏微分动力学方程转化为常微分方程,并将方程无量纲化得:
其中q=w/h为无量纲位移,w为板中点位移,h为板厚;
τ=t/γ为无量纲时间;
为无量纲基频,β=Ly/Lx,Lx和Ly分别为板的长和宽,β小于1;
ξ为名义粘性阻尼;
s=T/T0为屈曲系数,T为板的温升,T0为板发生热屈曲的临界温升,板内各处温升相等,
其中,分别为单位温升引起的x,y方向的热内力以及xy平面内的剪切力,
κ为立方项系数
其中Lz=zN+1-z1为板的总厚度,为时间缩放系数,[A*]=[A]-1为拉伸刚度矩阵的逆矩阵,
f=pcoefficientp(t)为激励项,p(t)为板上施加的声压,为载荷转换系数。
进一步的,步骤4具体包括:
考虑式(4)对应的静平衡方程,即令加速度项和速度项为零,得到力与位移的关系:
绘制s取不同值时式(8)对应的刚度曲线,其中力f转换为声压p,无量纲位移q转换为有量纲位移w;在s>1情况下,取极值点B和C的纵坐标即为跳变临界声压pcr:
pcr=fcr/pcoefficient (9)
其中
进一步的,步骤5具体包括:
假设声压幅值为高斯白噪声,其幅值p为零均值的正态分布,其概率密度函数表示为:
对给定声压级SPL的高斯白噪声,其均方根值为
Prms=2×10(SPL/20-5) (13)
由于均值为零,均方根值Prms就是标准差σ,因此上式转换为
σ=2×10(SPL/20-5) (14)。
进一步的,步骤6具体包括:
绘制概率密度函数式(12)对应的钟形曲线,并将跳变临界声压pcr标记在横坐标上,横坐标绝对值大于pcr的曲线围成的面积为跳变发生的概率。
进一步的,钟形曲线的宽度正比于标准差σ,pcr=nσ时,跳变刚开始发生,联立(9)(10)(11)(14)式得跳变临界声压级:
跳变概率的表达式为:
相对于现有技术,本发明具有以下有益效果:
1)从跳变机理出发结合噪声载荷概率密度函数直观地揭示了突弹跳变的本质;
2)从概率角度出发,能全面覆盖无跳变、间歇性跳变和持续性跳变响应类型,并给出跳变发生的概率,所得跳变边界可用于直观评估安全裕度;
3)本方法基于理论分析,推导出跳变临界声压级的表达式,能快速给出不同尺寸和材料参数复合材料板在一定温度和噪声声压级范围内的跳变边界,有助于把握结构的全局特性,加速优化设计迭代进程,相比于有限元和实验方法有巨大的时间成本和人力物力成本优势。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明方法的流程图;
图2为刚度曲线;
图3为声压概率密度分布与跳变临界声压;
图4为理论与有限元跳变边界对比。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明,旨在对本发明提供进一步的详细说明。除非另有指明,本发明所采用的所有技术术语与本申请所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式,而并非意图限制根据本发明的示例性实施方式。
请参阅图1所示,本发明提供一种基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,包括以下步骤:
步骤1:获取待预测复合材料板铺层信息、材料和几何参数;
步骤2:计算待预测复合材料的总体刚度矩阵;
步骤3:热噪声载荷下薄壁结构动力学方程推导,获取无量纲单模态方程;
步骤4:获取不同温度下结构的刚度曲线,利用刚度曲线极值点获得跳变临界力fcr关于屈曲系数s的函数表达式;
步骤5:列出给定声压级SPL(Sound pressure level)的噪声载荷概率密度函数g(p);
步骤6:比较噪声载荷的概率密度函数和跳变临界力对应的临界声压pcr,推导跳变临界声压级表达式,即跳变边界,并给出对应的跳变概率。
步骤1中,获取复合材料板的铺层信息,包括铺设层数N、每层厚度h、每层铺设方向θ。获取每层的材料参数,包括弹性模量E1,E2,泊松比ν12,ν21,剪切模量G12,G13,G23,密度ρ,热膨胀系数α1,α2。几何参数为长、宽、厚:Lx,Ly,h。
步骤2中,首先确定每个单一铺层刚度矩阵,第k层的刚度矩阵为:
其中,
然后计算总体刚度矩阵
其中,Aij,Bij,Dij分别为拉伸、拉弯耦合、弯曲刚度矩阵元素,z为厚度方向坐标,zk+1和zk分别为第k层上下表面坐标,T(k)为第k层坐标转换矩阵。
θ(k)为第k层铺设角。
需要说明的是,单层板及其他均质材料板可视为复合材料板的特例,只需计算单层刚度矩阵即可,本方法同样适用。
步骤3的方程推导过程如下:基于矩形板横向振动微分方程及变形协调方程,考虑温度引起的热应力,在本构关系中加入热内力;基于冯·卡门大挠度板理论,考虑非线性应变位移关系。采用假设模态法,将响应表示成双重傅里叶级数并截取第一项,应用伽辽金法将板的偏微分动力学方程转化为常微分方程,并将方程无量纲化得
其中q=w/h为无量纲位移,w为板中点位移,h为板厚;
τ=t/γ为无量纲时间;
为无量纲基频,β=Ly/Lx,Lx和Ly分别为板的长和宽,β须小于1;
ξ为名义粘性阻尼;
s=T/T0为屈曲系数,T为板的温升,T0为板发生热屈曲的临界温升,板内各处温升相等,
其中,分别为单位温升引起的x,y方向的热内力以及xy平面内的剪切力,
κ为立方项系数
其中Lz=zN+1-z1为板的总厚度,为时间缩放系数,[A*]=[A]-1为拉伸刚度矩阵的逆矩阵,
f=pcoefficientp(t)为激励项,p(t)为板上施加的声压,为载荷转换系数。
步骤中4,考虑式(4)对应的静平衡方程,即令加速度项和速度项为零,得到力与位移的关系
绘制s取不同值时式(8)对应的曲线如图2所示,此即刚度曲线,其中力f转换为声压p,无量纲位移q转换为有量纲位移w。在s>1情况下,取极值点B和C的纵坐标即为跳变临界声压pcr。
pcr=fcr/pcoefficient(9)
其中
进一步,所述步骤5,假设声压幅值为高斯白噪声,其幅值p为零均值的正态分布,因此其概率密度函数可表示为
对给定声压级SPL的高斯白噪声,其均方根值为
Prms=2×10(SPL/20-5) (13)
由于均值为零,均方根值Prms就是标准差σ,因此上式转换为
σ=2×10(SPL/20-5) (14)
步骤6中,绘制概率密度函数式(12)对应的钟形曲线,并将跳变临界声压pcr标记在横坐标上,如图3所示,横坐标绝对值大于pcr的曲线围成的面积(图中阴影)即为跳变发生的概率。钟形曲线的宽度正比于标准差σ,假设pcr=nσ时,跳变刚开始发生,则联立(9)(10)(11)(14)式可得跳变临界声压级。
这就是跳变边界表达式。
跳变概率的表达式为
本实施例中,预测该复合材料板的跳变边界及特定声压级下的跳变概率的方法,包括以下步骤:
步骤1:获取该复合材料的几何和材料参数。该复合材料尺寸为300×200×3mm,材料参数如表1所示。该材料采用整体成型,且厚度极薄,因此其材料特性采用工程常数表示,不涉及具体铺层信息,在后续计算中视作单一铺层处理,板的材料主方向与坐标系平行,即沿矩形板的长宽方向。
表1材料参数
步骤2:将表1参数代入(1)至(3)式,可得总体刚度矩阵[K],式中空白处为0。
步骤3:所述方程详细推导过程见于文献[2],此处将直接根据步骤3对各变量的说明及具体表达式,获得无量纲单模态方程的各项系数,各项参数如下:ω0=1.3777,β=0.6667,ξ取0.0019,T0=252.76,κ=0.8173,γ=6.24e-4,pcoefficient=2.60e-5。
步骤4:将步骤3获得的各项参数代入式(8),并绘制刚度曲线(如图2所示),其横纵坐标已转化为有量纲的声压和位移。在s>1(即热屈曲)情况下,如s=2,取刚度曲线极值点B和C的纵坐标即为跳变临界声压pcr,pcr由(9)至(11)式求出,图中零点A和A’为热屈曲平衡位置。为计算简洁,中间推导过程暂不带入具体数值。
步骤5:由式(12)至(14)确定给定声压级SPL的高斯白噪声声压的概率密度函数。为计算简洁,中间推导过程暂不带入具体数值。
步骤6:以屈曲系数s=1.2,声压级SPL=165dB为例,此时pcr=3.8kPa,σ=3.56kPa绘制概率密度函数式(12)对应的钟形曲线,并将跳变临界声压pcr标记在横坐标上,如图3所示,横坐标绝对值大于pcr的曲线围成的面积(图中阴影)即为跳变发生的概率。根据3σ准则,当pcr=3σ时,非跳变区刚好占据钟形曲线的大部分面积,在此基础上若增大声压级SPL,则钟形曲线向两侧拓宽,跳变区面积从几乎为零开始逐渐增加,因而跳变概率也逐渐增加,因此3σ准则可作为确定跳变边界的依据。将n=3及步骤3中各项参数代入式(15)可得理论跳变边界表达式。
为验证所提方法的准确性,建立相应的有限元模型,对s=1.2、1.4、1.6、1.8、2.0及SPL=151~180dB,间隔1dB,共计150个组合工况进行扫描计算,观察获得每个屈曲系数下发生跳变的最小声压级,构成跳变边界。比较本方法所得理论跳变边界与有限元计算的跳变边界如图4所示,两者最大误差1.2dB,最大相对误差0.69%,因此本发明所提方法具有较高的精度。对给定的屈曲系数和声压级,由式(16)可得跳变概率。同样以s=1.2,SPL=165dB为例,代入(16)式,可得跳变概率Psnap=0.286。
对任意给定的屈曲系数和声压级组合工况(s,SPL),从图4可以立即判断该工况下是否发生跳变,如果发生跳变则可以进一步给出跳变概率,如无跳变发生,也可以根据该工况点(s,SPL)离跳变边界的距离远近直观判断安全裕度。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
Claims (8)
1.基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,包括以下步骤:
步骤1:获取待预测复合材料板铺层信息、材料和几何参数;
步骤2:计算待预测复合材料的总体刚度矩阵;
步骤3:热噪声载荷下薄壁结构动力学方程推导,获取无量纲单模态方程;
步骤4:获取待预测复合材料板的刚度曲线,利用刚度曲线极值点获得跳变临界力fcr关于屈曲系数s的函数表达式;
步骤5:列出给定声压级SPL的噪声载荷概率密度函数g(p);
步骤6:比较噪声载荷的概率密度函数和跳变临界力对应的临界声压pcr,推导跳变临界声压级表达式,并给出对应的跳变概率。
2.根据权利要求1所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤1中:
获取复合材料板的铺层信息,包括铺设层数N、每层厚度h、每层铺设方向θ;
获取每层的材料参数,包括弹性模量E1,E2,泊松比ν12,ν21,剪切模量G12,G13,G23,密度ρ,热膨胀系数α1,α2;
获取几何参数,包括长Lx、宽Ly、厚h。
3.根据权利要求1所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤2中,首先确定每个单一铺层刚度矩阵,第k层的刚度矩阵为:
其中,
然后计算总体刚度矩阵
其中,Aij,Bij,Dij分别为拉伸、拉弯耦合、弯曲刚度矩阵元素,z为厚度方向坐标,zk+1和zk分别为第k层上下表面坐标,T(k)为第k层坐标转换矩阵;
θ(k)为第k层铺设角。
4.根据权利要求1所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤3具体包括:
采用假设模态法,将响应表示成双重傅里叶级数并截取第一项,应用伽辽金法将板的偏微分动力学方程转化为常微分方程,并将方程无量纲化得:
其中q=w/h为无量纲位移,w为板中点位移,h为板厚;
τ=t/γ为无量纲时间;
为无量纲基频,β=Ly/Lx,Lx和Ly分别为板的长和宽,β小于1;
ξ为名义粘性阻尼;
s=T/T0为屈曲系数,T为板的温升,T0为板发生热屈曲的临界温升,板内各处温升相等,
其中,分别为单位温升引起的x,y方向的热内力以及xy平面内的剪切力,
κ为立方项系数
其中Lz=zN+1-z1为板的总厚度,为时间缩放系数,[A*]=[A]-1为拉伸刚度矩阵的逆矩阵,
f=pcoefficientp(t)为激励项,p(t)为板上施加的声压,为载荷转换系数。
5.根据权利要求1所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤4具体包括:
考虑式(4)对应的静平衡方程,即令加速度项和速度项为零,得到力与位移的关系:
绘制s取不同值时式(8)对应的刚度曲线,其中力f转换为声压p,无量纲位移q转换为有量纲位移w;在s>1情况下,取极值点B和C的纵坐标即为跳变临界声压pcr:
pcr=fcr/pcoefficient (9)
其中
6.根据权利要求1所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤5具体包括:
假设声压幅值为高斯白噪声,其幅值p为零均值的正态分布,其概率密度函数表示为:
对给定声压级SPL的高斯白噪声,其均方根值为
Prms=2×10(SPL/20-5) (13)
由于均值为零,均方根值Prms就是标准差σ,因此上式转换为
σ=2×10(SPL/20-5) (14)。
7.根据权利要求6所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,步骤6具体包括:
绘制概率密度函数式(12)对应的钟形曲线,并将跳变临界声压pcr标记在横坐标上,横坐标绝对值大于pcr的曲线围成的面积为跳变发生的概率。
8.根据权利要求7所述的基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法,其特征在于,
钟形曲线的宽度正比于标准差σ,pcr=nσ时,跳变刚开始发生,联立(9)(10)(11)(14)式得跳变临界声压级:
跳变概率的表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811637938.9A CN109766613B (zh) | 2018-12-29 | 2018-12-29 | 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811637938.9A CN109766613B (zh) | 2018-12-29 | 2018-12-29 | 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109766613A true CN109766613A (zh) | 2019-05-17 |
CN109766613B CN109766613B (zh) | 2020-10-27 |
Family
ID=66453214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811637938.9A Active CN109766613B (zh) | 2018-12-29 | 2018-12-29 | 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109766613B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539107A (zh) * | 2020-04-23 | 2020-08-14 | 北京玻钢院复合材料有限公司 | 一种格栅增强复合材料圆柱壳体临界轴压载荷快速预测方法 |
CN112861329A (zh) * | 2021-01-22 | 2021-05-28 | 郑州轻工业大学 | 热管理下电热微驱动器的输出特性快速计算与设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
EP2450172B1 (en) * | 2010-11-04 | 2014-12-03 | Rolls-Royce plc | Method of manufacturing a composite material hollow axisymmetrical body |
CN105760577A (zh) * | 2016-01-28 | 2016-07-13 | 北京航空航天大学 | 一种含不确定性金属结构声振疲劳寿命的估算方法 |
CN106529090A (zh) * | 2016-12-12 | 2017-03-22 | 中国航天标准化研究所 | 一种航天电子类产品可靠性评估方法 |
CN107220461A (zh) * | 2017-06-26 | 2017-09-29 | 大连理工大学 | 一种变刚度复合材料板壳结构高效优化方法 |
-
2018
- 2018-12-29 CN CN201811637938.9A patent/CN109766613B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2450172B1 (en) * | 2010-11-04 | 2014-12-03 | Rolls-Royce plc | Method of manufacturing a composite material hollow axisymmetrical body |
CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
CN105760577A (zh) * | 2016-01-28 | 2016-07-13 | 北京航空航天大学 | 一种含不确定性金属结构声振疲劳寿命的估算方法 |
CN106529090A (zh) * | 2016-12-12 | 2017-03-22 | 中国航天标准化研究所 | 一种航天电子类产品可靠性评估方法 |
CN107220461A (zh) * | 2017-06-26 | 2017-09-29 | 大连理工大学 | 一种变刚度复合材料板壳结构高效优化方法 |
Non-Patent Citations (3)
Title |
---|
Y. SHA等: "Nonlinear response characteristics of thinwalled structures under thermo-acoustic loadings", 《ACTA AERONAUTICA ASTRONAUTICA SINICA》 * |
吕冰洋: "热环境下热防护结构动响应行为研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
李纪永等: "热-噪声载荷作用下薄壁结构非线性动态响应分析", 《机械科学与技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539107A (zh) * | 2020-04-23 | 2020-08-14 | 北京玻钢院复合材料有限公司 | 一种格栅增强复合材料圆柱壳体临界轴压载荷快速预测方法 |
CN111539107B (zh) * | 2020-04-23 | 2023-05-09 | 北京玻钢院复合材料有限公司 | 一种格栅增强复合材料圆柱壳体临界轴压载荷快速预测方法 |
CN112861329A (zh) * | 2021-01-22 | 2021-05-28 | 郑州轻工业大学 | 热管理下电热微驱动器的输出特性快速计算与设计方法 |
CN112861329B (zh) * | 2021-01-22 | 2022-10-25 | 郑州轻工业大学 | 热管理下电热微驱动器的输出特性快速计算与设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109766613B (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Przekop et al. | Dynamic snap-through of thin-walled structures by a reduced-order method | |
CN109033705B (zh) | 一种空间褶皱薄膜弹性应变能和主应力有限元分析方法 | |
Khorshidi et al. | Flutter analysis of sandwich plates with functionally graded face sheets in thermal environment | |
CN109766613B (zh) | 基于概率理论的热噪声载荷下薄壁结构突弹跳变预测方法 | |
CN112580236B (zh) | 一种热防护连接结构非线性动态响应的快速分析方法 | |
Chen et al. | Periodic response of blades having three-dimensional nonlinear shroud constraints | |
Liu et al. | Nonlinear dynamic response and structural evaluation of container ship in large freak waves | |
Coaquira et al. | Parametric instability and nonlinear oscillations of an FRP channel section column under axial load | |
Perraud et al. | Laminar-turbulent transition on aerodynamic surfaces with imperfections | |
Golmakani et al. | Nonlinear bending analysis of ring-stiffened circular and annular general angle-ply laminated plates with various boundary conditions | |
Viau et al. | Dynamic analysis methods for modelling timber assemblies subjected to blast loading | |
Beisheim et al. | On the three-dimensional finite element analysis of dovetail attachments | |
AsadiGorgi et al. | Effects of all-over part-through cracks on the aeroelastic characteristics of rectangular panels | |
Locke et al. | Finite element, large-deflection random response of thermally buckled beams | |
Zhou et al. | Grouping response method for equivalent static wind loads based on a modified LRC method | |
Gru¨ ber et al. | The impact of viscous effects on the aerodynamic damping of vibrating transonic compressor blades—a numerical study | |
Mao et al. | Comparative study of flange-to-seal contact couplings with bolt relaxation under creep condition | |
Li et al. | Finite element analysis of the ultimate strength of aluminum-stiffened panels with fixed and floating transverse frames | |
Spottswood et al. | Reduced-order models for a shallow curved beam under combined loading | |
Wang et al. | A compact cfd-based reduced order modeling for gust analysis | |
Nayak et al. | Linear transient dynamic analysis of plates with and without cutout | |
Alanbay et al. | Lowest twelve frequencies of sandwich plates using third-order shear-normal deformation theory | |
Hoshi et al. | Numerical simulations of transonic flutter on a three-dimensional wing | |
Bhagat et al. | Vibrational characteristics of truss core sandwich panel under thermal environment: effect of core topology | |
Alder | Development and validation of a partitioned fluid-structure solver for transonic panel flutter with focus on boundary layer effects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |