CN106055772A - 基于波长的微振动分析有限元建模单元尺度确定方法 - Google Patents
基于波长的微振动分析有限元建模单元尺度确定方法 Download PDFInfo
- Publication number
- CN106055772A CN106055772A CN201610362822.3A CN201610362822A CN106055772A CN 106055772 A CN106055772 A CN 106055772A CN 201610362822 A CN201610362822 A CN 201610362822A CN 106055772 A CN106055772 A CN 106055772A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- wavelength
- formula
- platform structure
- 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 30
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 239000000463 material Substances 0.000 claims description 10
- 238000005452 bending Methods 0.000 claims description 9
- 238000010008 shearing Methods 0.000 claims description 7
- 239000007787 solid Substances 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 abstract description 2
- 230000005284 excitation Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 229920000049 Carbon (fiber) Polymers 0.000 description 5
- 239000004917 carbon fiber Substances 0.000 description 5
- 230000001413 cellular effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical group C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 210000000744 eyelid Anatomy 0.000 description 4
- 230000002787 reinforcement Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种基于波长的微振动分析有限元建模单元尺度确定方法,步骤S1:确定平台结构所受激励的频率范围,并以最高激励频率的1.5~2倍作为基准频率。步骤S2:分析卫星平台结构系统的组成,选择平台结构系统中的典型结构作为基本的建模对象,并计算各种典型结构的等效力学参数。步骤S3:确定各种典型结构中振动波在其中的传递速度。步骤S4:根据波速计算平台结构中各种典型结构中振动波的最短波长。步骤S5:在最短波长内取离散点数为20~25,进而得到平台结构中各种典型结构在建模时的单元尺度。本发明得到的单元尺度可直接用于平台结构的有限元建模,在保证模型计算精度收敛的同时,可以有效减小计算代价。
Description
技术领域
本发明涉及动力学分析技术领域,具体的涉及一种基于波长的微振动分析有限元建模单元尺度确定方法。
背景技术
航天器上运动部件的工作会导致局部或整体的微幅、宽频、不造成结构破坏的微振动,其影响高精度航天器上有效载荷的正常工作。通常采用有限元方法建立卫星平台结构的动力学模型,之后对微振动及其在航天器平台中的传递进行研究。然而,由于整星结构的复杂性以及微振动分析的高精度需求,对有限元建模时的单元尺度选择提出了较高的要求。
目前对于有限元建模时单元尺度的选择大多是通过建模者的经验进行选取,而有限元建模时单元尺度的选择,关系到模型复杂度、计算代价、计算精度。单元尺度的选择通常需要满足计算精度收敛性的要求,即单元尺度的减小不再显著改善计算精度。因此,现有技术中为了获得合适的单元尺度,需要选择不同的单元尺度多次重复建立有限元模型,不断重复尝试对比多次计算精度,然后从中选择合适的单元尺度。但是对于复杂的航天器平台结构而言,重复建模所带来的工作量巨大,几乎成为不可能完成的任务。
发明内容
为了解决上述技术问题,本发明提供了一种基于波长的微振动分析有限元建模单元尺度确定方法。
本发明提供的一种基于波长的微振动分析有限元建模单元尺度确定方法,包括以下步骤:
步骤S1:根据平台结构上安装的振源,确定平台结构所受激励的频率范围,并以最高激励频率的1.5~2倍作为基准频率ω;
步骤S2:分析卫星平台结构系统的组成,选择平台结构中的各典型结构作为基本建模对象,将各典型结构等效为各向同性结构,并计算各典型结构的等效弹性模量与等效密度;
步骤S3:对所得基本建模对象确定其动力学建模时采用的有限单元类型,有限单元类型包括实体单元、板单元、梁单元,并分别计算各有限单元类型的各项波速;
步骤S4:取步骤S3中所得的各项波速中各种弹性波的最小波速cmin,根据公式(8)计算平台结构中各种典型结构中振动波的最短波长λ,
其中,ω为激励频率,此处的各种弹性波包括步骤S3中的纵向拉伸波、剪切横向波、弯曲波、扭转剪切波;
步骤S5:在最短波长λ内取离散点数n为20~25,得到平台结构中各种典型结构在建模时的单元尺度为λ/n。
进一步地,根据公式(1)计算得出实体单元中的纵向拉伸波的波速clS:
其中,E为材料弹性模量,ρ为材料密度,μ为泊松比;
根据公式(2)计算实体单元中的剪切横向波的波速cSS:
进一步地,根据公式(3)计算板单元中的纵向拉伸波的波速clp:
根据公式(4)计算得出板单元中弯曲波的波速cfp:
其中,ω为激励频率,h为板的厚度。
进一步地,根据公式(5)计算得出梁单元中的纵向拉伸波的波速clb:
根据公式(6)计算得出梁单元中的扭转剪切波的波速ctb:
其中,G为剪切模量,J为横截面的极惯性矩,I为横截面的惯性矩。
进一步地,根据公式(7)计算得到梁部件中的弯曲波的波速cfb:
其中,S为梁截面面积。
本发明的技术效果:
本发明提供的基于波长的微振动分析有限元建模单元尺度确定方法,可直接确定单元尺度,且能够保证计算精度收敛。
本发明提供的基于波长的微振动分析有限元建模单元尺度确定方法,针对航天器平台结构微振动分析要求,对有限元建模时的单元尺度得出量化的结果,以便运用于对微振动分析时有限元建模的单元尺度规定中,从而实现在保证计算精度收敛的情况下,直接快速的确定单元尺度,减少反复尝试单元尺度带来的巨大工作量。
具体请参考根据本发明的基于波长的微振动分析有限元建模单元尺度确定方法提出的各种实施例的如下描述,将使得本发明的上述和其他方面显而易见。
附图说明
图1是本发明提供的基于波长的微振动分析有限元建模单元尺度确定方法的流程示意图;
图2是本发明优选实施例中所处理某卫星平台部分舱段结构外形示意图;
图3是本发明优选实施例中所处理某卫星平台部分舱段结构内部示意图(舱段的上、下盖板未在图中示出);
图4是对比例中单元最大尺度为0.1m时舱段的有限元模型一示意图;
图5是本发明优选实施例中单元最大尺度为0.02m时舱段的有限元模型二示意图;
图6是对比例中单元最大尺度为0.01m时舱段的有限元模型三示意图;
图7是模型一~三的有限元模型频响计算结果对比曲线。
具体实施方式
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
参见图1,本发明提供的基于波长的微振动分析有限元建模单元尺度确定方法,包括以下步骤:
步骤S1:根据平台结构上安装的振源,确定平台结构所受激励的频率范围,并以最高激励频率的1.5~2倍作为基准频率ω。此处的平台结构可以为卫星等常用航天器作为平台,针对具体平台的具体问题进行参数的确定。
步骤S2:分析卫星平台结构系统的组成,选择平台结构中的各典型结构作为基本建模对象,将各典型结构等效为各向同性结构,并计算各典型结构的等效弹性模量与等效密度。此处的典型结构是指各种航天器在轨使用的不同材料、不同截面参数的梁、板或实体结构,例如3cm厚的M40蒙皮/铝蜂窝夹芯板。此处所用的等效弹性模量和等效密度均按现有方法进行计算。
步骤S3:对所得基本建模对象确定其动力学建模时采用的有限单元类型,有限单元类型包括实体单元、板单元、梁单元,并分别计算各有限单元类型的各项波速;此处对于各项波速的计算可以按常规方法进行。
步骤S4:取步骤S3中所得的各项波速中各种弹性波的最小波速cmin,根据公式(8)计算平台结构中各种典型结构中振动波的最短波长λ,
其中,ω为激励频率,此处的各种弹性波包括步骤S3中的纵向拉伸波、剪切横向波、弯曲波、扭转剪切波。
步骤S5:在最短波长λ周期内取离散点数n为20~25,得到平台结构中各种典型结构在建模时的单元尺度为λ/n。
利用弹性波在结构中传递的波长来确定有限元建模时的单元尺度,并在每个周期波长内取离散点数为20~25。此时仿真计算结果能在较短时间内快速收敛,从而减少了计算量,保证仿真结果的精度。
优选的,对于实体单元,根据公式(1)计算得出实体单元中的纵向拉伸波的波速clS:
其中,E为材料弹性模量,ρ为材料密度,μ为泊松比;
根据公式(2)计算实体单元中的剪切横向波的波速cSS:
按此方法计算能提高仿真结果的精度。
优选的,对于板单元,根据公式(3)计算板单元中的纵向拉伸波的波速clp:
根据公式(4)计算得出板单元中弯曲波的波速cfp:
其中,ω为激励频率,h为板的厚度。
按此方法计算能提高仿真结果的精度。
优选的,针对梁单元,根据公式(5)计算得出梁单元中的纵向拉伸波的波速clb:
根据公式(6)计算得出梁单元中的扭转剪切波的波速ctb:
其中,G为剪切模量,J为横截面的极惯性矩,I为横截面的惯性矩。
按此方法计算能提高仿真结果的精度。
优选的,根据公式(7)计算得到梁部件中的弯曲波的波速cfb:
其中,S为梁截面面积。
按此方法计算能提高仿真结果的精度。
以下结合具体实例对本发明提供的基于波长的微振动分析有限元建模单元尺度确定方法仅需详细的说明。
以某卫星平台结构的部分舱段为对象,如图2和图3所示,具体描述基于波长的微振动分析有限元建模单元尺度确定方法,下面结合附图说明具体实施方式。实例中,平台结构由含加强筋的舱板组成,舱板均为碳纤维面板/铝蜂窝夹芯材料,加强筋为碳纤维方管。坐标轴X、Y、Z的定义如图2所示。
步骤S1:某卫星平台结构上的激励源包括SADA、飞轮、控制力矩陀螺等,其微振动激励频率位于5Hz~200Hz之间,因此选择基准频率ω=300Hz。
步骤S2:分析平台结构系统的组成,选择平台结构系统中的典型结构作为基本的建模对象,平台结构的整体尺寸为825mm×2680mm×1150mm,平台结构主要由含加强筋的板件构成,其中板件分别命名为东板、西板、南板、北板、中隔板、上板和下板,东板为3cm厚的M40蒙皮蜂窝板,其他均为25.6cm厚的铝蒙皮蜂窝板。因此以M40碳纤维方管、3cm厚的M40蒙皮蜂窝板、25.6cm厚的铝蒙皮蜂窝板作为典型结构。计算得到各种结构的等效力学参数如下表所示。
表1本实施例中所处理典型结构的等效力学参数表
名称 | 等效弹性模量(GPa) | 等效密度(Kg/m3) | 泊松比 |
M40碳纤维方管 | 180 | 1600 | 0.34 |
M40蒙皮蜂窝夹芯板 | 2.34 | 38.5 | 0.34 |
铝蒙皮蜂窝夹芯板 | 1.64 | 66.1 | 0.34 |
步骤S3:对各建模对象确定动力学建模时采用的有限单元类型。舱板厚度方向上的尺寸远小于其他两个方向上的尺寸,可选取板单元进行结构离散;对于其中的加强筋可以采用梁单元进行建模。由于平台结构中仅采用板单元和梁单元进行建模,因此仅计算板单元和梁单元中振动波的波长。
针对板单元,根据下式计算得出板单元中的纵向拉伸波的波速clp
其中,E为材料弹性模量,ρ为材料密度,μ为泊松比。根据下式计算得出板单元中弯曲波的波速cfp
其中,ω为激励频率,h为板的厚度。计算可得到M40蒙皮蜂窝夹芯板中300Hz拉伸弹性波的波速为8290m/S,弯曲弹性波的波速为135m/S。铝蒙皮蜂窝夹芯板中300Hz拉伸弹性波的波速为5297m/S,弯曲弹性波的波速为108m/S。
针对梁单元,根据下式计算得出梁单元中的纵向拉伸波的波速clb
根据下式计算得出梁单元中的扭转剪切波的波速为ctb
其中,G为剪切模量,J为梁横截面的极惯性矩,I为梁横截面的惯性矩。梁部件中的弯曲波的波速为cfb
其中,S为梁截面面积。计算得到M40碳纤维方管中,300Hz拉伸弹性波的波速为3354m/S,300Hz扭转弹性波的波速为7855m/S。300Hz弯曲弹性波的波速为169m/S。
步骤S4:根据步骤S3得到的波速,得到各种弹性波波速的最小值为108m/S,计算得到平台结构中各种典型结构中振动波的最短波长为0.36m,取整约为0.4m。
步骤S5:通过大量仿真显示,在每个周期内选择离散点数为20~25时,计算结果已经收敛,即继续增加离散点数不会明显提高求解精度,因此选择单元尺度为0.4m/20,即0.02m。
对如图2所示的部分平台结构系统,利用Patran/Nastran软件进行有限元建模,为了对比分析,分别选择单元尺度为0.1m、0.02m和0.01m,所得到的模型分别如图4(称为模型一)、图5(称为模型二)和图6(称为模型三)所示。为便于对图5(本发明提供方法得到的单元尺度)与其他对比例进行对比分析,人为设定图4和图6中模型一、三中的的单元尺度。图4中模型一总共1505个单元;图5中模型二的单元尺度由本发明方法给出,总共36854个单元;图6中模型三总共14750个单元。可以看出随着单元尺度的降低,模型规模迅速增加,由此必然带来计算代价的迅速增加。
在图2中平台结构底板的A点作用5~200Hz的单位简谐加速度激励,并且将图2中平台结构的C、D、E、F点设置为固支边界条件,分析得到各个模型中图2中平台结构上板的B点加速度响应输出如图7所示。从图7可以观察到,模型二计算得到的B点加速度响应曲线与模型三的几乎完全重合,而模型一得到的加速度响应曲线与模型三差别较大,不具有对比意义。由此可以判定模型二计算精度已经收敛。三个模型在同一台计算机上进行计算,模型一的计算时间为7s,模型二的计算时间为168s,模型三的计算时间为2150s。综上可以看出,相对于模型三,根据本发明方法给出的单元尺度建立的模型二计算代价显著降低,且计算精度已经收敛。从而验证了通过本发明计算得到的单元尺度在保证模型计算收敛的同时,可以减小计算代价。
本领域技术人员将清楚本发明的范围不限制于以上讨论的示例,有可能对其进行若干改变和修改,而不脱离所附权利要求书限定的本发明的范围。尽管己经在附图和说明书中详细图示和描述了本发明,但这样的说明和描述仅是说明或示意性的,而非限制性的。本发明并不限于所公开的实施例。
通过对附图,说明书和权利要求书的研究,在实施本发明时本领域技术人员可以理解和实现所公开的实施例的变形。在权利要求书中,术语“包括”不排除其他步骤或元素,而不定冠词“一个”或“一种”不排除多个。在彼此不同的从属权利要求中引用的某些措施的事实不意味着这些措施的组合不能被有利地使用。权利要求书中的任何参考标记不构成对本发明的范围的限制。
Claims (5)
1.一种基于波长的微振动分析有限元建模单元尺度确定方法,其特征在于,包括以下步骤:
步骤S1:根据平台结构上安装的振源,确定平台结构所受激励的频率范围,并以最高激励频率的1.5~2倍作为基准频率ω;
步骤S2:分析卫星平台结构系统的组成,选择平台结构中的各典型结构作为基本建模对象,将各典型结构等效为各向同性结构,并计算各典型结构的等效弹性模量与等效密度;
步骤S3:对所得基本建模对象确定其动力学建模时采用的有限单元类型,有限单元类型包括实体单元、板单元、梁单元,并分别计算各有限单元类型的各项波速;
步骤S4:取步骤S3中所得的各项波速中各种弹性波的最小波速cmin,根据公式(8)计算平台结构中各种典型结构中振动波的最短波长λ,
其中,ω为激励频率,此处的各种弹性波包括步骤S3中的纵向拉伸波、剪切横向波、弯曲波、扭转剪切波;
步骤S5:在最短波长λ内取离散点数n为20~25,得到平台结构中各种典型结构在建模时的单元尺度为λ/n。
2.根据权利要求1所述的基于波长的微振动分析有限元建模单元尺度确定方法,其特征在于,根据公式(1)计算得出实体单元中的纵向拉伸波的波速clS:
其中,E为材料弹性模量,ρ为材料密度,μ为泊松比;
根据公式(2)计算实体单元中的剪切横向波的波速cSS:
3.根据权利要求1所述的基于波长的微振动分析有限元建模单元尺度确定方法,其特征在于,根据公式(3)计算板单元中的纵向拉伸波的波速clp:
根据公式(4)计算得出板单元中弯曲波的波速cfp:
其中,ω为激励频率,h为板的厚度。
4.根据权利要求1所述的基于波长的微振动分析有限元建模单元尺度确定方法,其特征在于,根据公式(5)计算得出梁单元中的纵向拉伸波的波速clb:
根据公式(6)计算得出梁单元中的扭转剪切波的波速ctb:
其中,G为剪切模量,J为横截面的极惯性矩,I为横截面的惯性矩。
5.根据权利要求1所述的基于波长的微振动分析有限元建模单元尺度确定方法,其特征在于,根据公式(7)计算得到梁部件中的弯曲波的波速cfb:
其中,S为梁截面面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362822.3A CN106055772B (zh) | 2016-05-27 | 2016-05-27 | 基于波长的微振动分析有限元建模单元尺度确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362822.3A CN106055772B (zh) | 2016-05-27 | 2016-05-27 | 基于波长的微振动分析有限元建模单元尺度确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106055772A true CN106055772A (zh) | 2016-10-26 |
CN106055772B CN106055772B (zh) | 2019-02-26 |
Family
ID=57175552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362822.3A Expired - Fee Related CN106055772B (zh) | 2016-05-27 | 2016-05-27 | 基于波长的微振动分析有限元建模单元尺度确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106055772B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488812A (zh) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | 卫星在轨微振动仿真建模修正方法 |
CN103605834A (zh) * | 2013-10-31 | 2014-02-26 | 上海卫星工程研究所 | 卫星飞轮被动隔振系统设计方法 |
CN103712759A (zh) * | 2014-01-06 | 2014-04-09 | 北京卫星环境工程研究所 | 航天器整星微振动力学环境地面试验方法 |
-
2016
- 2016-05-27 CN CN201610362822.3A patent/CN106055772B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488812A (zh) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | 卫星在轨微振动仿真建模修正方法 |
CN103605834A (zh) * | 2013-10-31 | 2014-02-26 | 上海卫星工程研究所 | 卫星飞轮被动隔振系统设计方法 |
CN103712759A (zh) * | 2014-01-06 | 2014-04-09 | 北京卫星环境工程研究所 | 航天器整星微振动力学环境地面试验方法 |
Non-Patent Citations (5)
Title |
---|
LI D, JIANG J, LIU W, ET AL.: "A new mechanism for the vibration control of large flexible space structures with embedded smart devices", 《 IEEE/ASME TRANS. MECHATRONICS》 * |
YIHUAN L, DAOKUI L I, GUOJIN T.: "Motion Planning for Vibration Reducing of Free-floating Redundant Manipulators Based on Hybrid Optimization Approach", 《CHINESE JOURNAL OF AERONAUTICS》 * |
孟光,周徐斌,苗军: "航天重大工程中的力学问题", 《力学进展》 * |
廖一寰, 李道奎, 唐国金: "自由漂浮柔性机械臂系统的动力学建模与仿真", 《国防科技大学学报》 * |
马兴瑞,韩增尧,邹元杰,丁继锋: "航天器力学环境分析与条件设计研究进展", 《宇航学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106055772B (zh) | 2019-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nguyen et al. | An analytical solution for buckling and vibration analysis of functionally graded sandwich beams using a quasi-3D shear deformation theory | |
Ghafoori et al. | Dynamic analysis of laminated composite plates traversed by a moving mass based on a first-order theory | |
Nayak et al. | Process for design optimization of honeycomb core sandwich panels for blast load mitigation | |
Quinn et al. | Non-prismatic sub-stiffening for stiffened panel plates—Stability behaviour and performance gains | |
CN103336871A (zh) | 一种确定复合材料厚度方向平面内等效工程常数的方法 | |
Wennberg et al. | Orthotropic models of corrugated sheets in finite element analysis | |
Kuo | Flutter of thermally buckled angle-ply laminates with variable fiber spacing | |
Cho et al. | Ultimate strength formulation for stiffened plates subjected to combined axial compression, transverse compression, shear force and lateral pressure loadings | |
Li et al. | A piecewise shear deformation theory for free vibration of composite and sandwich panels | |
CN104573372A (zh) | 一种网状可展开天线展开过程索力分析方法 | |
Mondal et al. | Dynamic instability of damped composite plates with embedded delaminations | |
Rodrigues et al. | Experimental and finite element modal analysis of variable stiffness composite laminated plates | |
CN106055772A (zh) | 基于波长的微振动分析有限元建模单元尺度确定方法 | |
Marjanović et al. | Object-oriented framework for 3D bending and free vibration analysis of multilayer plates: Application to cross-laminated timber and soft-core sandwich panels | |
Schreiber et al. | A holistic approach for local buckling of composite laminated beams under compressive load | |
Nasihatgozar et al. | General equations for free vibrations of thick doubly curved sandwich panels with compressible and incompressible core using higher order shear deformation theory | |
Zhou et al. | Analysis of dynamic fracture parameters in functionally graded material plates with cracks by graded finite element method and virtual crack closure technique | |
Lopatin et al. | Fundamental frequency and design of the CFCF composite sandwich plate | |
Panahandeh-Shahraki et al. | Nonlinear buckling analysis of laminated composite curved panels constrained by Winkler tensionless foundation | |
Rodriguez et al. | Static aeroelastic analysis with an inviscid Cartesian method | |
Tokura et al. | Forming process simulation of truss core panel | |
Sorrenti et al. | Free vibration analysis of angle-ply laminated and sandwich plates using enhanced refined zigzag theory | |
Liu et al. | On the prediction improvement of transverse stress distributions in cross-ply laminated beams: advanced versus conventional beam modelling | |
Ranaivomiarana et al. | Optimal orthotropy and density distribution of two-dimensional structures | |
Adámek et al. | A two-layer elastic strip under transverse impact loading: Analytical solution, finite element, and finite volume simulations |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190226 |