CN108614922B - 一种基于应力杂交化后处理的板壳边缘效应的计算方法 - Google Patents
一种基于应力杂交化后处理的板壳边缘效应的计算方法 Download PDFInfo
- Publication number
- CN108614922B CN108614922B CN201810300367.3A CN201810300367A CN108614922B CN 108614922 B CN108614922 B CN 108614922B CN 201810300367 A CN201810300367 A CN 201810300367A CN 108614922 B CN108614922 B CN 108614922B
- Authority
- CN
- China
- Prior art keywords
- stress
- plate shell
- boundary
- edge effect
- displacement
- 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.)
- Active
Links
- 230000000694 effects Effects 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 title claims abstract description 36
- 238000009396 hybridization Methods 0.000 title claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000005452 bending Methods 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 10
- 239000000463 material Substances 0.000 claims description 9
- 238000012805 post-processing Methods 0.000 claims description 8
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 4
- 238000007670 refining Methods 0.000 claims description 4
- 230000000295 complement effect Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- NTSBMKIZRSBFTA-AIDOXSFESA-N Digoxigenin bisdigitoxoside Chemical compound C1[C@H](O)[C@H](O)[C@@H](C)O[C@H]1O[C@@H]1[C@@H](C)O[C@@H](O[C@@H]2C[C@@H]3[C@]([C@@H]4[C@H]([C@]5(CC[C@@H]([C@@]5(C)[C@H](O)C4)C=4COC(=O)C=4)O)CC3)(C)CC2)C[C@@H]1O NTSBMKIZRSBFTA-AIDOXSFESA-N 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000005034 decoration Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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]
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于应力杂交化后处理的板壳边缘效应的计算方法，属于CAE仿真、工程结构力学计算领域，具体为：1）基于Mindlin‑Reissner理论建立板壳结构的力学模型；2）使用有限元软件进行求解，得到结构整体的位移结果；3）对结构自由边界临近区域进行应力杂交化后处理，得到应力分布并提取出应力峰值以及峰值所出现的位置作为参考指标；4）将有限元网格模型精细化处理后重复步骤2与步骤3的操作，并判断参考指标是否满足收敛性条件以或精度要求。本发明解决了工程分析中板壳结构边缘效应问题难以准确高效计算的缺陷，该方计算效率更高，能够求解更广泛的工程结构，具有更好的灵活性与适用性。
Description
技术领域
本发明属于CAE仿真、工程结构力学计算领域，具体是指一种基于应力杂交化后处理的板壳边缘效应的计算方法。
背景技术
板壳结构是一类非常重要的结构形式，被广泛地应用于航空航天、土建、船舶等工程领域。板壳结构边缘效应是一种典型的局部应力现象，是指在某些工况下自由边界邻近非常狭小的区域内，应力分布会出现非常剧烈的梯度变化。而这一狭小的区域称为板壳的边界层，尺度与板壳厚度同量级。在一些工程应用领域中，例如复合材料层合板(Vidal P,Polit O,D′Ottavio M,et al.Assessment of the refined sinus plate finiteelement:Free edge effect and Meyer-Piening sandwich test.Finite Elements inAnalysis and Design,2014,92(0):60-71.)、超大型浮式海洋结构(Wang C M,Xiang Y,Utsunomiya T,et al.Evaluation of modal stress resultants in freely vibratingplates.International Journal of Solids and Structures,2001,38(36-37):6525-6558.)，该区域内的应力峰值极有可能会导致板壳结构的破坏或失效，因此准确的求解边缘效应现象非常重要。
在绝大多数的工程应用中，板壳结构分析均会使用有限元方法进行求解。但是常规的有限元方法是一种基于分片插值的纯数值方法。受限于单元的插值阶次不足，其在处理边缘效应这类局部应力问题时表现十分乏力。针对边缘效应问题的现有技术大致可以分为两类。第一类是通过网格自适应技术来对板壳结构的边界层区域进行处理，通过减少边界层区域的网格尺寸来提高有限元计算的精度(Schwab C,Suri M.The p and hpversions of the finite element method for problems with boundarylayers.Mathematics of Computation,1996,65(216):1403-1429.)。这是比较常用的方式。但是由于边缘效应的作用范围远小于结构的整体尺寸，自适应计算的收敛速度会受到很大的影响。特别是在分析显式动力学问题时，边界层内过小的单元尺寸会严重降低整体计算的稳定时间步长，从而导致计算成本的显著增加。第二类则是使用特殊构造的有限元模型来处理边界层区域。例如Ramesh等(Ramesh S S,Wang C M,Reddy J N,etal.Computation of stress resultants in plate bending problems using higher-order triangular elements.Engineering Structures,2008,30(10):2687-2706.)通过大幅度增加节点自由度数量来提高插值阶次。但是这种做法同样会导致高昂的计算成本。
本发明人也曾通过构造特殊的有限元模型来有效地求解该类问题(Shang Y,CenS,Li C F,et al.An effective hybrid displacement function element method forsolving the edge effect of Mindlin-Reissner plate.International Journal forNumerical Methods in Engineering,2015,102(8):1449-1487.)(Shang Y,Cen S,Li Z,et al.Improved Hybrid Displacement Function(IHDF)element scheme for analysisof Mindlin-Reissner plate with edge effect.International Journal forNumerical Methods in Engineering,2017,111(12):1120-1169.)，但是其依赖于特定的有限元分析程序，其在通用性以及求解效率上均有明显不足。此外，该方法只能够用于分析各向同性材料的平板结构边缘效应问题。
发明内容
本发明针对现有技术中存在的问题，公开了一种基于应力杂交化后处理的板壳边缘效应的计算方法，该方法是首先基于有限元软件强大的求解能力，利用求解结果对有限元网格模型进行应力杂交化后处理，本发明方法计算效率更高，能够求解更广泛的工程结构，具有更好的灵活性与适用性。解决了现有技术中存在的问题。
本发明是这样实现的：
一种基于应力杂交化后处理的板壳边缘效应的计算方法，其特征在于，步骤如下：
步骤一：建立的板壳结构模型：基于Mindlin-Reissner理论对实际的板壳结构进行简化，建立相应的力学模型；M-R理论考虑了板壳结构的横向剪切变形，能够有效地解释边缘效应存在的原因；因此在本方法中，具体选用Mindlin-Reissner板壳理论。
步骤二：针对步骤一建立的板壳结构力学模型；使用有限元软件对其进行求解，得到结构整体的位移、应变能、应力、应变；本发明中使用的有限元软件中商业的有限元软件均可，为方便也可以指定为具体的软件，如Abaqus。)商业有限元软件具备强大的求解能力，该方法充分利用其强大的求解能力，计算效率更高，能够求解更广泛的工程结构，具有更好的灵活性与适用性。
步骤三：针对板壳结构发生边缘效应现象的部分，即自由边界的临近区域，提取出由步骤二计算得到的该区域的位移；然后对该区域进行应力杂交化后处理，从而获得边界层区域内的应力分布，并提取应力峰值与所出现的位置作为评估参数；
步骤四：对板壳结构的有限元网格模型减少网格尺寸，(可以将网格的特征尺寸减到一半或者根据实际的用户需要减少网格尺寸)增加网格密度，然后重复步骤二以及步骤三，得到一组新的应力峰值与位置；评估前后两次计算所得的结果是否满足用户预先设定的收敛准则以及或所需的计算精度；是否满足的判据是两次计算的相对误差低于某一个值。这个值通常由用户根据自己的要求确定，例如可取1％或更低。值越低，计算成本越高。
4.1，如若不满足，则进一步通过减少网格尺寸、增加网格密度来细化有限元网格模型，并重复步骤二到步骤四，然后再次进行收敛性判断，直至满足收敛准则以及所需的计算精度；
4.2，如若满足，即得到符合要求的板壳结构边缘效应的内力结果。
进一步，所述的步骤一具体为：
1.1，M-R理论考虑了板壳结构的横向剪切变形，能够有效地解释边缘效应存在的原因。使用Mindlin-Reissner板壳理论可以对板壳结构的厚度方向给予几何建模，从而将结构简化为一个二维结构，使用其中性层的二维拓扑来进行表征。
在Mindlin-Reissner理论中，板壳结构每一个点有6个自由度，包括3个平动自由度和3个转动自由度。将这些自由度可以表示为下述矢量形式：
q＝[u v w θx θy θz]T (1)
式中，q为自由度矢量，u,v,w为三个坐标方向的平动自由度，θx，θy，θz为三个方向的转动自由度。
根据变形模式特征，这6个自由度又可以分为两组，分别对应着面内变形部分：
qin＝[u v θz]T (2)
以及横向弯曲变形部分：
qout＝[w θx θy]T (3)
1.2，通过理论分析可知边缘效应局部应力现象主要需要考虑其弯曲变形模式，Mindlin-Reissner理论中与弯曲变形有关的广义应变分量和应力分量为别为
Ε＝[κx κy κxy γxz γyz]T (4)
R＝[Mx My Mxy Tx Ty]T (5)
式中，E为广义应变分量，κx,κx,κxy为弯曲曲率，γxz,γyz为横向剪切应变；R为应力分量，Mx,Mx,Mxy为弯矩和扭矩，Tx,Ty为横向剪力。
其中，扭矩与剪力分量因为边缘效应会在很小的范围内发生剧烈的梯度变化，应变与位移之间的得到：
κx＝-θx,x,κy＝-θy,y,κxy＝-θx,y-θy,x (6)
γx＝w,x-θx,γy＝w,y-θy (7)
应力需满足平衡方程：
Mx,x+Mxy,y-Tx＝0,Mxy,x+My,y-Ty＝0,Tx,x+Ty,y+q＝0 (8)
式中，q是横向分布载荷；
应变与应力之间的本构关系可以表示为
D是材料的弹性矩阵，具体表达需要根据板壳结构的具体材料决定,Db是与弯曲变形相关的子矩阵，Ds是与横向剪切变形相关的子矩阵。
进一步，所述的步骤二具体为：
2.1，使用有限元软件对板壳结构模型进行分析求解，得到整个结构的位移与应力结果。
2.2，根据结构的边界约束类型，甄别出其自由边界部分，确定需要进行应力杂交化后处理的区域；
2.3，提取2.2中所述区域内的网格信息以及节点位移结果，并从中分离出与离面弯曲变形有关的位移分量。
进一步，为了求解边缘效应，应力杂交化后处理主要是对弯曲变形部分进行。所述的步骤三具体为对弯曲变形部分进行应力杂交化后处理，过程如下：
3.1，应力杂交化后处理是对弯曲变形部分进行的，其核心思想是由杂交元技术衍生而来。首先设定一组弯曲应力场，即：
R＝Sβ (10)
β是一组待定参数，S是一组应力试探函数；
3.2，根据步骤二得到的结构整体的位移qout设定边界位移模式：
d＝Nqout (11)
选用梁的变形模式来近似模拟板壳结构的边界变形，由此可得N的具体表达；
3.3，上述弯曲应力场与边界位移在每个单元中需要满足最小余能原理，即
式中第一项为域内积分，Ω表示控制域；第二项为边界积分，Γ表示控制域的边界；
3.4，为了准确反应出自由边界处部分内力为零的特殊性，需要使用Lagrange乘子法对式(13)进行修正，即将其变换为：
式中，
表示在自由边界s上的内力矢量；λ是Lagrange乘子；通过推导，可以得到待定参数β与步骤二所得位移qout之间的关系：
式中M、
是计算过程中所会用到的中间矩阵，其可以由下述的步骤计算得到：
M＝∫∫ΩSTD-1SdΩ (16)
Hλ＝(ΓM-1ΓT)-1ΓM-1H (18)
H＝∫ΓSTLTNdΓ (19)
Γ＝∫SLSdS (20)
L是方向余弦矩阵；
3.5，在求解式(15)后再将β代入到式(10)中即可得到后处理之后的应力结果；此外，对式(10)求导可以进一步得到各应力分量的峰值以及其所出现的位置；
3.6，边缘效应是一种典型的局部应力现象，在很小的尺度范围内，应力会出现非常显著的梯度变化。因此式(10)中的试探函数S必须具备有效模拟该现象的能力。本方法的一个关键点就是将具有下述指数函数形式的应力试函数扩充到式(10)中：
式中m，n是在切平面内该自由边界的方向余弦；K是与结构厚度以及材料属性相关的参数；
是多项式函数，可取不同的阶次，具体选取视需求而定。
本发明与现有技术相比的有益效果在于：
相比较于现有技术中使用纯数值的有限元方法以及依赖于特定的有限元模型、有限元分析程序，本发明首先利用成熟的有限元分析软件进行整体结构的求解，再针对板壳结构发生边缘效应现象的部分进行杂交化后处理，根据是否满足收敛性需求进行精细化处理，最终得到符合要求的板壳结构边缘效应的内力结果。首先在分析程序求解边缘效应问题上，本发明的方法计算效率更高，能够求解更广泛的工程结构，具有更好的灵活性与适用性；此外，与通过网格细分求解边缘效应问题的常规做法相比，本发明的方法有效的克服了细分网格导致计算规模过大的缺点，显著地提高了计算效率与计算精度。
附图说明
图1为本发明的方法原理流程图；
图2为运用本发明方法求解边缘效应的实例示意图；
图3为运用本发明方法求解实例所得到的扭矩分量和剪力分量的分布曲线示意图。
具体实施方式
为使本发明的目的、技术方案及效果更加清楚，明确，以下参照附图并举实例对本发明进一步详细说明。应当指出此处所描述的具体实施仅用以解释本发明，并不用于限定本发明。
如图1所示，图1为一种基于应力杂交化后处理的板壳边缘效应的计算方法，主要分为4个部分：
第一部分：基于Mindlin-Reissner理论对实际的板壳结构进行简化，建立相应的力学模型用于求解；具体步骤为：
使用Mindlin-Reissner板壳理论可以对板壳结构的厚度方向不予几何建模，从而将结构简化为一个二维结构，使用其中性层的二维拓扑来进行表征。在M-R理论中，每一个节点有6个自由度，包括3个平动自由度和3个转动自由度。将这些自由度可以表示为下述矢量形式：
q＝[u v w θx θy θz]T (1)
式中，q为自由度矢量，u,v,w为三个坐标方向的平动自由度，θx，θy，θz为三个方向的转动自由度；
根据变形模式，这6个自由度又可以分为两组，分别对应着面内变形部分与离面弯曲变形部分，即
qin＝[u v θz]T (2)
qout＝[w θx θy]T (3)
边缘效应局部应力现象主要需要考虑其弯曲变形模式。根据M-R理论，与弯曲变形有关的广义应变分量和应力分量为别为
Ε＝[κx κy κxy γxz γyz]T (4)
R＝[Mx My Mxy Tx Ty]T (5)
式中，E为广义应变分量，κx,κx,κxy为弯曲曲率，γxz,γyz为横向剪切应变；R为应力分量，Mx,Mx,Mxy为弯矩和扭矩，Tx,Ty为横向剪力；
其中，扭矩与剪力分量因为边缘效应会在很小的范围内发生剧烈的梯度变化。应变与位移之间的得到：
κx＝-θx,x,κy＝-θy,y,κxy＝-θx,y-θy,x (6)
γx＝w,x-θx,γy＝w,y-θy (7)
应力需满足平衡方程：
Mx,x+Mxy,y-Tx＝0,Mxy,x+My,y-Ty＝0,Tx,x+Ty,y+q＝0 (8)
应变与应力之间的本构关系可以表示为
式中弹性矩阵的具体表达需要根据板壳结构具体材料与铺层来决定。
第二部分：对于步骤一中建立的结构模型，使用有限元软件软件，例如Abaqus，对其进行整体求解，得到整个结构的位移。根据结构的边界约束类型，甄别出其自由边界部分，确定需要应力杂交化后处理的区域。提取该部分区域内的网格信息，以及从节点位移中提取出与离面弯曲变形有关的位移分量，即式(3)，以用于进一步的后处理。
第三部分：针对板壳结构发生边缘效应现象的部分，即自由边界的临近区域，提取出由步骤二计算得到的该区域的位移；然后对该区域进行应力杂交化后处理，从而获得边界层区域内的应力分布，并提取应力峰值与所出现的位置作为评估参数。具体步骤为：
应力杂交化后处理是对弯曲变形部分进行的，其核心思想是由杂交元技术衍生而来。首先假设一组符合要求的弯曲应力场，即
R＝Sβ (10)
β是一组待定参数，S是一组应力试探函数。然后根据步骤二得到的位移qout假设一个合理的边界位移模式
d＝Nqout (11)
具体选用梁的变形模式来近似模拟板壳结构的边界变形，由此可以的N的具体表达。上述假设的应力场与边界位移在每个单元中需要满足最小余能原理，即
为了准确反应出自由边界处部分内力为零的特殊性，需要使用Lagrange乘子法对式(13)进行修正，即将其变换为：
通过推导可以得到待定参数β与步骤二所得位移qout之间的关系
式中M、
可由下述运算得到：
M＝∫∫ΩSTD-1SdΩ (16)
Hλ＝(ΓM-1ΓT)-1ΓM-1H (18)
H＝∫ΓSTLTNdΓ (19)
Γ＝∫SLSdS (20)
式中，L是方向余弦矩阵。在求解式(15)后再将β代入到式(10)中即可得到后处理之后的应力结果。此外，对式(10)求导可以进一步得到各应力分量的峰值以及其所出现的位置。
边缘效应是一种典型的局部应力现象，具有强烈的梯度变化。因此式(10)中的试探函数S必须具备有效模拟该现象的能力。对此需要将具有下述指数函数形式的应力试函数扩充到式(10)中：
式中m，n是在切平面内该自由边界的方向余弦；K是与结构厚度以及材料属性相关的参数；
是多项式函数，可取不同的阶次。
第四部分具体为：
将板壳结构的有限元网格模型进行细化，减少网格尺寸(将网格的特征尺寸减到一半，或者根据用户实际需要减少网格尺寸)，增加网格密度，然后重复步骤二与步骤三，获得边缘效应一组新的内力峰值与峰值所在位置。评估前后两次计算结果的收敛性是否满足精度的要求。
若不满足要求，则再次对有限元网格模型进行精细化处理，并重复步骤二到步骤四的分析计算。若满足要求，则输出最终次的计算结果。
图2是一个边缘效应简单算例的示意图：边长为a的一个矩形方板承受横向均布载荷作用。左右边界为自由边，上下边界为简支约束。利用对称性，可以只取四分之一结构，即阴影部分，进行建模并分析。
如图3所示，由应力杂交化后处理得到的边缘效应扭矩分量和剪力分量分布曲线示意图，3a所示的是扭矩分量沿AB边分布；3b所示的是剪力分量沿AB边分布。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于应力杂交化后处理的板壳边缘效应的计算方法，其特征在于，步骤如下：
步骤一：基于Mindlin-Reissner理论对实际的板壳结构进行简化，建立相应的力学模型；
步骤二：针对步骤一建立的板壳结构力学模型；使用有限元软件对其进行求解，得到结构整体的位移、应变能、应力、应变结果；
步骤三：针对板壳结构发生边缘效应现象的部分，即自由边界的临近区域，提取出由步骤二计算得到的该区域的位移；然后对该区域进行应力的杂交化后处理，从而获得自由边界的临近区域内的应力分布，并提取应力峰值与所出现的位置作为评估参数，具体为：
对横向弯曲变形部分进行应力杂交化后处理，过程如下：
3.1，首先设定一组弯曲应力场，即
R＝Sβ (10)
R是应力矢量，S是一组应力试探函数，β是相应的待定参数；
3.2，根据步骤二得到的横向弯曲变形部分qout设定边界位移：
d＝Nqout (11)
选用梁的变形模式来近似模拟板壳结构的边界变形，由此可得N的具体表达；
3.3，上述弯曲应力场与边界位移在每个单元中需要满足最小余能原理，即
式中第一项为域内积分，Ω表示控制域；第二项为边界积分，Γ表示控制域的边界；D是材料的弹性矩阵；
3.4，使用Lagrange乘子法对式(13)进行修正，将自由边界处的零值力边界条件加以考虑，即将其变换为：
式中，
表示在自由边界s上的内力矢量；λ是Lagrange乘子；通过推导，可以得到待定参数β与步骤二所得位移qout之间的关系：
式中M、
是计算过程所会用到的中间矩阵，其可以由下述的步骤计算得到：
M＝∫∫ΩSTD-1SdΩ (16)
Hλ＝(ΓM-1ΓT)-1ΓM-1H (18)
H＝∫ΓSTLTNdΓ (19)
Γ＝∫SLSdS (20)
L是方向余弦矩阵；
3.5，在求解式(15)后再将β代入到式(10)中即可得到后处理之后的应力结果；对式(10)求导可以进一步得到各应力分量的峰值以及其所出现的位置；
3.6，将下述指数函数形式的应力试函数扩充到式(10)中：
式中m，n是在切平面内该自由边界s的方向余弦；K是与结构厚度以及材料属性相关的参数；
是多项式函数；
步骤四：对板壳结构的有限元网格模型减少网格尺寸，增加网格密度，然后重复步骤二以及步骤三，得到一组新的应力峰值与位置；评估前后两次计算所得的结果是否满足用户预先设定的收敛准则或设定的计算精度；
4.1，如若不满足，则进一步通过减少网格尺寸、增加网格密度来细化有限元网格模型，并重复步骤二到步骤四，然后再次进行收敛性判断，直至满足收敛准则以及所需的计算精度；
4.2，如若满足，即得到符合要求的板壳结构边缘效应的内力结果。
2.根据权利要求1所述的一种基于应力杂交化后处理的板壳边缘效应的计算方法，其特征在于，所述的步骤一具体为：
1.1，在Mindlin-Reissner理论中，板壳结构每一个点有6个自由度，包括3个平动自由度和3个转动自由度；将这些自由度可以表示为下述矢量形式
q＝[u v w θx θy θz]T (1)
式中，q为自由度矢量，u,v,w为三个坐标方向的平动自由度，θx，θy，θz为三个方向的转动自由度；
根据变形模式特征，这6个自由度又可以分为两组，分别对应着面内变形部分：
qin＝[u v θz]T (2)
以及横向弯曲变形部分：
qout＝[w θx θy]T (3)
1.2，通过理论分析可知边缘效应局部应力现象主要需要考虑其弯曲变形模式，Mindlin-Reissner理论中与弯曲变形有关的广义应变分量和应力分量为别为
Ε＝[κx κy κxy γxz γyz]T (4)
R＝[Mx My Mxy Tx Ty]T (5)
式中，E为广义应变分量，κx,κx,κxy为弯曲曲率，γxz,γyz为横向剪切应变；R为应力矢量，Mx,Mx,Mxy为弯矩和扭矩，Tx,Ty为横向剪力；
其中，扭矩与剪力分量因边缘效应发生梯度变化，应变与位移之间的得到：
κx＝-θx,x,κy＝-θy,y,κxy＝-θx,y-θy,x (6)
γx＝w,x-θx,γy＝w,y-θy (7)
应力需满足平衡方程：
Mx,x+Mxy,y-Tx＝0,Mxy,x+My,y-Ty＝0,Tx,x+Ty,y+q＝0 (8)
式中，q是横向分布载荷；
应变与应力之间的本构关系可以表示为：
式中，D是材料的弹性矩阵,Db是与弯曲变形相关的子矩阵，Ds是与横向剪切变形相关的子矩阵。
3.根据权利要求1所述的一种基于应力杂交化后处理的板壳边缘效应的计算方法，其特征在于，所述的步骤二具体为：
2.1，使用有限元软件对板壳结构模型进行分析求解，得到整个结构的位移、应变能、应力、应变结果；
2.2，根据板壳结构的边界约束类型，甄别出其自由边界部分，确定需要进行应力杂交化后处理的区域；
2.3，提取2.2中所述区域内的网格信息以及节点位移结果，并从中分离出与离面弯曲变形相关的位移分量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810300367.3A CN108614922B (zh) | 2018-04-04 | 2018-04-04 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810300367.3A CN108614922B (zh) | 2018-04-04 | 2018-04-04 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108614922A CN108614922A (zh) | 2018-10-02 |
| CN108614922B true CN108614922B (zh) | 2021-08-06 |
Family
ID=63659606
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810300367.3A Active CN108614922B (zh) | 2018-04-04 | 2018-04-04 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108614922B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112613148B (zh) * | 2020-12-30 | 2024-05-28 | 一重集团大连工程技术有限公司 | 一种基于数值分析变形数据的核电装备设计方法 |
| CN113722824B (zh) * | 2021-08-30 | 2024-01-12 | 江南造船(集团)有限责任公司 | 适用于有限元分析的船舶板结构简化方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1399750A (zh) * | 1999-10-15 | 2003-02-26 | 莫尔德弗洛有限公司 | 用于结构分析的设备和方法 |
| JP2004046380A (ja) * | 2002-07-09 | 2004-02-12 | Tadahiko Kawai | 近似解の収束手段 |
| CN106547999A (zh) * | 2016-12-08 | 2017-03-29 | 湖南大学 | 基于显式动力学有限元方法的拉延模结构分析 |
| CN106844907A (zh) * | 2017-01-04 | 2017-06-13 | 中南大学 | 回转窑筒体支承位置及筒体厚度的协同优化方法 |
| CN107292047A (zh) * | 2017-06-21 | 2017-10-24 | 贵州财经大学 | 一种基于最大位移的复杂板壳厚度快速优化设计方法 |
-
2018
- 2018-04-04 CN CN201810300367.3A patent/CN108614922B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1399750A (zh) * | 1999-10-15 | 2003-02-26 | 莫尔德弗洛有限公司 | 用于结构分析的设备和方法 |
| JP2004046380A (ja) * | 2002-07-09 | 2004-02-12 | Tadahiko Kawai | 近似解の収束手段 |
| CN106547999A (zh) * | 2016-12-08 | 2017-03-29 | 湖南大学 | 基于显式动力学有限元方法的拉延模结构分析 |
| CN106844907A (zh) * | 2017-01-04 | 2017-06-13 | 中南大学 | 回转窑筒体支承位置及筒体厚度的协同优化方法 |
| CN107292047A (zh) * | 2017-06-21 | 2017-10-24 | 贵州财经大学 | 一种基于最大位移的复杂板壳厚度快速优化设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| An edge-based smoothed finite element method (ES-FEM) with stabilized discrete shear gap technique for analysis of Reissner–Mindlin plates;H. Nguyen-Xuan;《Comput. Methods Appl. Mech. Engrg》;20090624;全文 * |
| New hybrid-Trefftz Mindlin-Reissner plate elements for efficiently modeling the edge zones near free/SS1 edges;Yan Shang;《Engineering Computations》;20180202;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108614922A (zh) | 2018-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Do et al. | Material optimization of tri-directional functionally graded plates by using deep neural network and isogeometric multimesh design approach | |
| JP7058902B2 (ja) | ハイブリッド繊維複合材料の板巻きシェル構造に対する高速協調最適化方法 | |
| Zhou et al. | Vibration and flutter analysis of supersonic porous functionally graded material plates with temperature gradient and resting on elastic foundation | |
| Cirak et al. | Integrated modeling, finite-element analysis, and engineering design for thin-shell structures using subdivision | |
| Ngo-Cong et al. | Free vibration analysis of laminated composite plates based on FSDT using one-dimensional IRBFN method | |
| CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
| Khoei et al. | Crack growth modeling via 3D automatic adaptive mesh refinement based on modified-SPR technique | |
| CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
| Lu et al. | Parallel mesh adaptation for high-order finite element methods with curved element geometry | |
| CN110955941B (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
| Ahn et al. | Response surface method for airfoil design in transonic flow | |
| CN108614922B (zh) | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 | |
| Nefovska-Danilovic et al. | Shear deformable dynamic stiffness elements for a free vibration analysis of composite plate assemblies–Part I: Theory | |
| CN116011301B (zh) | B样条等几何状态空间有限元方法 | |
| Li et al. | Modeling crack propagation with the extended scaled boundary finite element method based on the level set method | |
| Ma et al. | Generative design of stiffened plates based on homogenization method | |
| CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
| CN106650073B (zh) | 一种平衡力系边界条件子模型分析方法 | |
| Deng et al. | Prediction of transonic flow over supercritical airfoils using geometric-encoding and deep-learning strategies | |
| Kourepinis et al. | Higher-order discontinuous modeling of fracturing in concrete using the numerical manifold method | |
| CN116757026B (zh) | 一种基于等几何分析的板架结构极限强度分析实现方法 | |
| Park et al. | Multidisciplinary design optimization of a structurally nonlinear aircraft wing via parametric modeling | |
| Siegler et al. | Supersonic airfoil shape optimization by variable-fidelity models and manifold mapping | |
| De Barcellos et al. | AC k continuous generalized finite element formulation applied to laminated Kirchhoff plate model | |
| Yang et al. | Numerical analysis of large-amplitude ship motions using FV-based Cartesian grid method |
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 |