【隐式动态求解器】使用非线性Newmark方法的隐式动态求解器（Matlab代码实现）

  💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

 ⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

使用非线性Newmark方法的隐式动态求解器及示例

函数调用： Result=Newmark_Nonlinear(Elements,Material,Support,Free,M,C,f,fs,delta)

输入参数： Elements：一个结构体，包含各元素的信息。Elements{i}.DOFs=[j k]表示元素i连接节点j与k的自由度；Elements{i}.Material=m表示元素i采用材料m的属性。

Material：一个结构体，包含双线性弹簧的材料属性。其中，Material{m}.k1为弹簧初始刚度；Material{m}.x1为弹簧刚度开始减小的变形阈值；Material{m}.k2为减小后的弹簧刚度。

Support：一个向量，表示固定支座的自由度，大小为(nSupport,1)。

Free：一个向量，表示自由度，大小为(nFree,1)。

M：质量矩阵，大小为(nFree*nFree)。

C：阻尼矩阵，大小为(nFree*nFree)。

f：外部力矩阵，大小为(nFree,N)，其中N为动态力数据点的数量。

fs：采样频率。

delta：残余力的收敛标准。

输出： Result：一个结构体，包含以下内容： Result.Displacement：位移矩阵，大小为(nFreeN)。 Result.Velocity：速度矩阵，大小为(nFreeN)。 Result.Acceleration：加速度矩阵，大小为(nFree*N)。

注释：假设元素为连接节点的双线性刚度弹簧（无迟滞效应）。

📚2 运行结果

部分代码：

clear; clc; close all;

%Input

M=5*diag(ones(3,1));               %Mass matrix
C=0.01*M;                           %Damping matrix
f=randn(3,10000);                  %force
fs=100;                            %Sampling frequency
delta=10e-4;                       %Convergance criterion for residual force

%Material types definition
Material{1}.k1=1000;               %Spring stiffness
Material{1}.x1=0.01;               %Spring deformation beyond which the stiffness decreases
Material{1}.k2=10;                 %Reduced stiffness

%Elements connectivity and material used
Elements{1}.DOFs=[1 2];    Elements{1}.Material=1;     
Elements{2}.DOFs=[2 3];    Elements{2}.Material=1;    
Elements{3}.DOFs=[3 4];    Elements{3}.Material=1;    
%Support and Free DOFs
Support=[1];               Free=[2 3 4];

Result=Newmark_Nonlinear(Elements,Material,Support,Free,M,C,f,fs,delta);

t=[0:1/fs:9999/fs];
figure;
subplot(3,1,1);
plot(t,Result.Displacement(1,:)); xlabel('Time'); ylabel('DSP2');
subplot(3,1,2);
plot(t,Result.Displacement(2,:)); xlabel('Time'); ylabel('DSP3');
subplot(3,1,3);
plot(t,Result.Displacement(3,:)); xlabel('Time'); ylabel('DSP4');

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]苟英旗.基于Newmark法的线性隐式数值积分方法及其在阻尼器减震结构动力分析中的应用[J].[2024-07-27].

[2]范玉川,黄清云,鲁艳,等.基于Newmark-β法的非线性体系动载荷识别[J].东北大学学报：自然科学版, 2019, 40(12):5.DOI:CNKI:SUN:DBDX.0.2019-12-016.

🌈4 Matlab代码实现

资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取