033_Structure_Static_In_Matlab求解结构静力学问题两套方法

🕗 发布于 2024-11-07 03:24 FEM 结构静力学 matlab 结构力学有限元

在这里插入图片描述

结构静力学问题

静力学问现在是已经很简单的问题，在材料各向同性的情况下，对于弹性固体材料，很容易通过有限元求解。特别是线弹性问题，方程的矩阵形式可以很容易的写出（准确得说是很容易通过有限元表达），然后通过求解线性方程组得到结果。

基本上，求解的就是如下形式的方程：

$\frac{\partial^2 u}{\partial t^2} + d \frac{\partial u}{ \partial t} - \nabla \cdot (c \nabla u) + a u = f$

其中， $m$ 是质量矩阵， $c$ 是阻尼矩阵， $a$ 是刚度矩阵， $f$ 是外力矩阵, $u$ 是不同问题中的因变量，结构问题中是位移，热传导问题中是温度，电磁场问题中是电势。而方程中的 $d$ 只有在时变问题中才会出现，是一个时间相关的粘性系数。

当然，当方程的变量有多个时，方程就可以写成矢量的形式：

$\mathbf{m} \frac{\partial^2 \mathbf{u}}{\partial t^2} + \mathbf{d} \frac{\partial\mathbf{u}}{\partial t} - \nabla \cdot (\mathbf{c} \otimes \nabla \mathbf{u}) + \mathbf{a} \mathbf{u} = \mathbf{f}$

相应的系数就变成矩阵或者向量，式中的 $\otimes$ 也就是张量积或者叫克罗内克积、Kronecker乘积。

接下来我们先介绍下传统方式的求解结构静力学问题，然后介绍新的统一模型的求解方式。

Matlab PDE Toolbox传统形式

这个工具箱最早的时候就实现了固体结构静力学、热传导、电磁场等问题的求解。

工具箱把不同领域的问题分成了不同的模型，采用同一个接口createpde来创建模型。这个接口的第一个参数是模型的名字，根据名称返回不同的模型对象。

structural：结构问题
thermal：热传导问题
electromagnetic：电磁场问题

每种问题，又包含了若干不同的分析类型，就拿结构问题来说，有：

"static-solid" | "static-planestress" | "static-planestrain" | "static-axisymmetric" | "transient-solid" | "transient-planestress" | "transient-planestrain" | "transient-axisymmetric" | "modal-solid" | "modal-planestress" | "modal-planestrain" | "modal-axisymmetric" | "frequency-solid" | "frequency-planestress" | "frequency-planestrain" | "frequency-axisymmetric"

首先，我们创建一个结构模型。

% 创建PDE模型
model = createpde('structural', 'static-solid');

model

model = 

  StructuralModel - 属性:

              AnalysisType: "static-solid"
                  Geometry: []
        MaterialProperties: []
                 BodyLoads: []
        BoundaryConditions: []
      ReferenceTemperature: []
    SuperelementInterfaces: []
                      Mesh: []
             SolverOptions: [1x1 pde.PDESolverOptions]

这个StructuralModel对象中，有各种属性，可以很容易的猜出其含义。比如MaterialProperties就是材料属性，BodyLoads就是体力，BoundaryLoads就是边界力，BoundaryConditions就是边界条件。

下面我们找一个最简单的问题，就是一个悬臂梁的问题，然后求解。

% 定义几何
gm = multicuboid(15, 1, 1);
model.Geometry = gm;
% 几何模型
figure
pdegplot(model, 'FaceLabels', 'on', 'FaceAlpha', 0.5);
title('Geometry with Face Labels - Domain Specific Modeling');

在这里插入图片描述

我们这里专门把面的编号写出来，方便后面的边界条件的设置。

材料特性都是瞎写的。

% 指定材料属性
structuralProperties(model, 'YoungsModulus', 210E9, 'PoissonsRatio', 0.3, 'MassDensity', 7800);

接下来就是产生网格，然后求解。

% 生成网格
generateMesh(model, 'Hmax', 0.2);

网格的最大边长设为0.1，这个值是可以调整的，这个值越小，网格越密，计算结果越精确，但是计算量也越大。

% 绘制网格
figure
pdeplot3D(model);
title('Mesh with Hmax = 0.2');

在这里插入图片描述

边界条件，我们把梁的右边固定，左边施加一个力。还记得加上重力的影响。看到我们这个力是一个向量，斜向下拉。

% 应用边界条件
structuralBC(model, 'Face', 3, 'Constraint', 'fixed');
structuralBoundaryLoad(model, 'Face', 5, 'SurfaceTraction', [0; -1E5; -1E5]);

% 重力
structuralBodyLoad(model, 'GravitationalAcceleration', [0; 0; -9.81]);

调用函数solve求解。

% 求解PDE
result = solve(model);

首先我们可以绘制整体的位移大小，基本上，就是经典的悬臂梁的感觉。不过，1米的悬臂梁还是很强很强的……

% 后处理结果
pdeplot3D(model, 'ColorMapData', result.Displacement.Magnitude);

title('Displacement Magnitude - Static Analysis - Domain Specific Modeling');

在这里插入图片描述

当然我们还可以计算应力，然后绘制应力的分布，应力集中的方式也很明显。

vonMisesStress = result.VonMisesStress;
pdeplot3D(model, 'ColorMapData', vonMisesStress);

title('VonMisesStress-Structural Static-Domain Specific Model');

在这里插入图片描述

统一求解框架

从2023a开始，Matlab开始引入了一个统一的求解框架来替代传统的各个专业的模型。这与工业软件的大的趋势也是一致的，多学科融合、数字化双创等等。

利用统一个数模、或者同一套网格，来求解不同的问题，这样可以减少重复工作，提高效率。

这个框架入口，就是femodel。就上面的悬臂梁问题，我们可以用这个框架来求解。

model = femodel(AnalysisType = "structuralStatic", Geometry = gm);

model

model = 

  1x1 femodel array

Properties for analysis type: structuralStatic

            AnalysisType: "structuralStatic"
                Geometry: [1x1 fegeometry]
      MaterialProperties: [0x1 materialProperties]

   Boundary Conditions
                  FaceBC: [0x6 faceBC]
                  EdgeBC: [0x12 edgeBC]
                VertexBC: [0x8 vertexBC]

   Loads
                CellLoad: [0x1 cellLoad]
                FaceLoad: [0x6 faceLoad]
                EdgeLoad: [0x12 edgeLoad]
              VertexLoad: [0x8 vertexLoad]

   Other Parameters
    ReferenceTemperature: []
           SolverOptions: [1x1 pde.PDESolverOptions]

  Show <a href="matlab:if exist('model','var'), displayAllProperties(model); else, disp('无法显示变量 model 的属性，因为该变量不再存在。'), end">所有属性</a>

这个对象也没有什么特别的，其属性也很容易从名字看出内容。

绘制模型的函数一点都没变……也能够认识这个对象

pdegplot(model, FaceLabels='on', FaceAlpha=0.5);
title('Geometry with Face Labels Displayed - femodel');

在这里插入图片描述

设置材料属性，采用了属性赋值的形式，并提供了一个新的函数materialProperties来设置材料属性。

划分网格的函数没有改变，唯一变化的时候必须把模型对象传递出去，不会直接在位修改输入参数的属性，这说明这个对象并不是handle对象。

最后，产生的网格记录在了model.Geometry.Mesh中，这个对象是一个FEMesh对象，这个对象也是一个新的对象，用来表示网格。

model.MaterialProperties = materialProperties('YoungsModulus', 210E9, 'PoissonsRatio', 0.3, 'MassDensity', 7800);
model = generateMesh(model);
pdeplot3D(model.Geometry.Mesh);

model.Geometry.Mesh

title('Mesh with default parameters - femodel');

ans = 

  FEMesh - 属性:

             Nodes: [3x1267 double]
          Elements: [10x600 double]
    MaxElementSize: 0.6027
    MinElementSize: 0.3013
     MeshGradation: 1.5000
    GeometricOrder: 'quadratic'

在这里插入图片描述

设置边界条件的方式也略有不同，从调用函数修改为了属性赋值的方式。我们可以看到，边界条件变成了一个数组。

FaceBC：面的边界条件，对应faceBC对象/函数
EdgeBC：边的边界条件，对应edgeBC对象/函数
VertexBC：点的边界条件，对应vertexBC对象/函数

载荷条件、初值条件等等也是类似的。

从这三个对象的数组，也可以看到模型包含面、线、点三种对象的数目。当然，设置边界条件的时候，也是按照这个顺序来设置的。

当然，载荷有一个体积载荷，对应CellLoad属性和cellLoad对象/函数。

model.FaceBC(3) = faceBC(Constraint='fixed');

model.FaceLoad(5) = faceLoad(SurfaceTraction=[0; -1E5; -1E5]);

model.CellLoad = cellLoad(Gravity=[0; 0; -9.81]);

当然，我们如果对于Matlab结构数组的操作比较熟悉，也可以对多个对象进行操作，比如设置多个边界条件。

比如面边界，经过我们设置之后，可以看到：

model.FaceBC

ans = 

  1x6 faceBC array

Properties for analysis type: structuralStatic

Index    Constraint    XDisplacement    YDisplacement    ZDisplacement
  1          []             []               []               []      
  2          []             []               []               []      
  3        fixed            []               []               []      
  4          []             []               []               []      
  5          []             []               []               []      
  6          []             []               []               []      

  Show <a href="matlab:if exist('ans','var'), displayAllProperties(ans); else, disp('无法显示变量 ans 的属性，因为该变量不再存在。'), end">所有属性</a>

model.FaceBC(2:4) = faceBC(Constraint='fixed');

model.FaceBC

ans = 

  1x6 faceBC array

Properties for analysis type: structuralStatic

Index    Constraint    XDisplacement    YDisplacement    ZDisplacement
  1          []             []               []               []      
  2        fixed            []               []               []      
  3        fixed            []               []               []      
  4        fixed            []               []               []      
  5          []             []               []               []      
  6          []             []               []               []      

  Show <a href="matlab:if exist('ans','var'), displayAllProperties(ans); else, disp('无法显示变量 ans 的属性，因为该变量不再存在。'), end">所有属性</a>

这样就很容易的设置了多个边界条件。

接下来就是基本一样的求解过程，调用solve函数，然后绘制结果。

% 先把面边界条件设回去，不然三个面被固定，结果会很好玩的
model.FaceBC([2, 4]) = faceBC();

ret = solve(model);

pdeplot3D(model.Geometry.Mesh, ColorMapData=ret.Displacement.Magnitude);
title('Displacement Magnitude-Structural Static-femodel');

在这里插入图片描述

vonMisesStress = ret.VonMisesStress;

pdeplot3D(model.Geometry.Mesh, ColorMapData=vonMisesStress);

title('VonMisesStress-Structural Static-femodel');

在这里插入图片描述

总的来说，统一接口的方式，也并没有什么大的不同之处。我感觉可能会让共享数模来分析不同学科的交叉问题会更加直观一些。当然，后面再学习的时候，建议还是直接从femodel开始，因为这个是未来的发展方向。

虽然，资料确实比较少，例子也还有好多没有调整过来。

总结

求解静力学问题，求解过程和结果都没有什么让人惊喜的地方。

下面，在讨论一个周期加载的问题，这个问题是一个动力学问题。

我们就直接用统一建模框架来做。

原文地址：https://blog.csdn.net/withstand/article/details/143449169

免责声明：本站文章内容转载自网络资源，如本站内容侵犯了原著者的合法权益，可联系本站删除。更多内容请关注自学内容网（zxcms.com）！

上一篇：写过的试卷怎么打印出新的？
下一篇：【字符串匹配算法】BF与KMP算法

量子力学（1）如何理解海森堡的不确定性原理
海森堡不确定性原理
阅读更多2024-11-09
5G-A如何与AI融合发展？华为MBBF2024给出解答
2024年，对于整个ICT产业来说，是极为重要的一年。一方面，随着3GPP R18标准的正式冻结，我们全面进入了5G-A时代，迎来了这项全新技术标准的规模商用。另一方面，在经历了“百模大战”的喧嚣之后
阅读更多2024-11-09
C++ Double Dispatch，即双重调度
对于初次接触或曾有所闻但未曾深入了解的开发人员而言，调度（Dispatching）这一概念值得我们深入探讨。调度，从字面意义上理解，即将某物发送至特定目的地。在C++编程中，调度同样指将控制权从一处转
阅读更多2024-11-09
编译器、IDE对C/C++新标准的支持
很遗憾的是，微软VC++对C99的支持直到VS2015才基本完全支持，在VS2012之前，仅仅支持C89. 也正是2015年之前，这么多人骂微软并纷纷转移到其它C语言编译器的原因。VS2019开始才正
阅读更多2024-11-09
css:基础
陶喆 {font-family: '宋体'.周杰伦 {可以看见部分字体改变了我们还可以设置多个字体，可以方便用户在查看网页的时候，如果没安装第一种字体，你们就执行第二种字体.eat {一般都使用系统默
阅读更多2024-11-09
FPGA课程内容升级【纯逻辑和图像开发】
手写axi lite控制器（master和slave）mipi-datamover-hdmi应用实例。mipi-axi-dma-hdmi应用实例。datamover读写 PL DDR。mipi-vdm
阅读更多2024-11-09
前端页面性能优化的常见问题与解决方案
在当今互联网高速发展的时代，前端页面的性能对于用户体验至关重要。一个加载缓慢、交互卡顿的页面很可能会导致用户流失。本文将深入探讨前端页面性能优化中常见的问题以及相应的解决方案。
阅读更多2024-11-09
linux常见命令
linux常见命令
阅读更多2024-11-09
Java后端请求想接收多个对象入参的数据方法
在Java后端开发中，如果我们希望接收多个对象作为HTTP请求的入参，可以使用Spring Boot框架来简化这一过程。Spring Boot提供了强大的RESTful API支持，能够方便地处理各种
阅读更多2024-11-09
使用 PageHelper 在 Spring Boot 项目中实现分页查询
在 Web 应用开发中，分页查询是非常常见的需求，特别是在涉及大量数据的应用场景中，通过分页可以减少数据加载压力，提升系统性能。然而，手动实现分页查询需要编写繁琐的 SQL 语句和逻辑代码，容易出现错
阅读更多2024-11-09

033_Structure_Static_In_Matlab求解结构静力学问题两套方法

结构静力学问题

Matlab PDE Toolbox传统形式

统一求解框架

总结

相关文章