实干、实践、积累、思考、创新！ Tag: 结构工程博士 结构工程师 程序员 结构抗震设计 振动控制 光伏支架 加固改造 软件研发 编程 弹塑性 有限元 超高层 超限设计
实干、实践、积累、思考，创新。 最近闲暇的时候翻翻 随机振动的书，又搞到傅里叶分析，功率谱等。 在测试网上的MATLAB功率谱分析程序的时候，发现功率谱居然是负数，(又暴露了知识面的短缺， 😥 ) 后仔思考发现，主要是因为功率谱使用了 分贝 为计量单位，比如通过公式 10*log10(Pxx) 将原本的功率谱Pxx进行转换，可见当Pxx处于0~1之间的时候，通过log10(Pxx)转换出来的公式就是负数了。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、实践、积累、思考，创新。 小伙伴让写个EEEP(Equivalent Energy Elastic Plastic) 等效屈服点及延性分析程序，用MATLAB测试编写。改天找时间添加到 YPD（[软件][科研][试验] 2020版 YPD: Yield Point and Ductility [2020版 YPD: 等效屈服点和延性分析程序]）软件。 相关博文( Related Topics) [01] [软件][科研][试验] 2020版 YPD: Yield Point and Ductility [2020版 YPD: 等效屈服点和延性分析程序] …
实干、实践、积累、思考、创新。 如题，一看是一个莫名其妙的想法。实际上也是一个错误的想法，不过，当时脑子一热，就测算一下。具体看看结果。 模型： 20层的剪切层模型。 阻尼矩阵： （1）模态阻尼，20阶振型计算 （2）模态阻尼，20阶振型计算，把非对角线元素取为0。 采用MATLAB编程，采用Newmark-β积分法进行弹性时程分析，两种阻尼模型的计算结果对比如下。 其中，参考阻尼为完整的模态阻尼，对比阻尼为去掉阻尼矩阵对角线元素后的矩阵。 由图可见，采用仅保留对角线元素的阻尼矩阵，结构的位移及剪力响应远小于完整的阻尼矩阵，而楼层加速度响应似乎相差不大！！十分诧异！！ 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
目前市面上，关于FEM理论(有限元理论)及FEM编程的书有许多，关于FEM软件应用的书也有许多，但关于两者之间结合的书比较少，这是目前市面书籍的gap，也是很多有限元初学者学习过程中存在的gap，因此我们写作了本书。本书将FEM基础理论、编程、及软件应用结合起来，讲完FEM基础理论，然后编程，最后采用结构工程师常用的几个软件进行同一个案例的分析并将分析结果与编程结果进行对比，当看到自己编写的FEM代码与商业软件的计算结果5个小数点内重合时，理论、编程、应用之间的gap自然而然就消除了。这正是我们编书的初衷，也是作者本人学习FEM基础理论，学习软件所采用的思路及方法。
因此这不仅仅是一本讲FEM，讲FEM编程或者讲FEM软件应用的书，而背后更重要的是，讲述作者推荐的一种学习FEM，学习软件的方法及思考过程，并不是纯粹的编著或者某些已有书籍的翻版。对于结构工程师及其他FEM理论应用者，十分重要的一点是在理解FEM基础理论的情况下灵活地运用FEM软件于实际问题。因此作者认为该书对于有限单元法的初学者及应用FEM解决实际工程问题的工程师有一定的用处。商品时代，给客户多一点选择。
花了点时间，整理了一下 www.jdcui.com的软件笔记，并汇总到以下这个页面，有需要的小伙伴可以收藏一下，持续更新。链接：http://www.jdcui.com/?page_id=2850。Abauqs，Ansys,ETABS,SAP2000,MIDAS,YJK,PKPM,XTRACT,MATLAB,ATUOCAD,PERFORM-3D,OPENSEES,HYPERMESH
实干、实践、积累、思考、创新。 接着博文《[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础]》继续介绍经典材料力学或结构力学课本上介绍的压杆稳定问题。该部分内容也是 书本 《有限单元法：编程与软件应用》屈曲分析章节的部分内容节选。 11.5 屈曲分析3：压杆稳定 作为屈曲分析的补充，本节讨论一下压杆稳定问题。 图 11‑10 压杆支座情况 算例结构为一根等截面轴心受压直杆，直杆材料为钢管，钢管外径100mm，管厚5mm，高5m。一共考虑了5种约束情况，分别为：1.两端铰接；2.一端铰接、一端嵌固；3.两端嵌固；4.一端嵌固、一端滑动；5.一端嵌固、一端自由。针对每种约束情况，分别将压杆划分为1个、2个、5个、10个、20个梁单元，进行屈曲分析。 11.5.1 MATLAB代码与注释 本节以底端嵌固、上端自由、划分20个单元的情况为例，给出进行屈曲分析所需的主要代码。 % Pressed Bar Buckling Analysis % Author ： JiDong Cui（崔济东），Xuelong Shen（沈雪龙） % Website : …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 今天是20181231,2018年的最后一天，希望新的一年会更好。
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。
实干、实践、积累、思考、创新。 《有限单元法：编程与软件应用》 曲屈分析 章节的内容节选： 11.1 稳定问题分类 结构失稳（屈曲）是指在外力作用下结构的平衡状态开始丧失，稍有扰动变形便迅速增大，最后使结构发生破坏。稳定问题一般分为两类，第一类是理想化的情况，即达到某种荷载时，除结构原来的平衡状态存在外，还可能出现第二个平衡状态，所以又称平衡分岔失稳或分支点失稳 (Buckling, Bifuraction)（图1中OAB曲线），对应于数学中是求解特征值问题，故又称特征值屈曲，此类结构失稳时相应的荷载称为屈曲荷载。第二类是结构失稳时，变形将迅速增大，而不会出现新的变形形式，又称极值点失稳(Instability)（图1中OCD曲线），结构失稳时相应的荷载称为极限荷载。此外还有一种跳跃失稳(Snap-through)（图1中OEF曲线），当荷载达到某值时（E点），结构平衡状态发生一个明显的跳跃，突然过渡到非临近的另一个具有较大位移的平衡状态（F点），由于在跳跃时结构通常已经破坏，其后的状态不能被利用，所以可归入第二类失稳。 图 11‑1 稳定问题 11.2 最小势能原理 11.2.1 系统的平衡 若系统处于平衡状态，则一定存在一种位移状态，使得系统总势能为驻值，即 (11.2‑1) 公式中 …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 接着博文《[动力学] 模态分析基础 (Basics of modal analysis)》继续介绍框架单元用于模态分析的方法。该部分内容也是 书本 《有限单元法：编程与软件应用》模态分析章节的部分内容节选。 10.3.1 梁单元质量矩阵 集中质量矩阵（局部坐标） 设杆材料密度为 ，单元长度为 ，截面积为 ，每个节点分担单元1/2的平动质量，无转动惯量，则单元质量矩阵 (10.3‑1) 10.3.2 算例：2D框架模态分析 算例采用与2D框架结构静力分析中相同的结构，采用欧拉梁单元，材料密度为2.5493e-9t/mm3，采用集中质量矩阵，将单元质量集中于两端节点，且不考虑集中质量后的节点转动质量。 由于算例模态分析的Matlab代码与前面章节中静力分析的Matlab代码大体相同，因此这里仅给进行模态分析所需的新增代码。 % Truss 2D modal analysis …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 接着博文《[动力学] 模态分析基础 (Basics of modal analysis)》继续介绍桁架单元用于模态分析的方法。该部分内容也是 书本 《有限单元法：编程与软件应用》模态分析章节的部分内容节选。 10.2 2D桁架结构模态分析 模态分析需要的参数主要有结构的刚度矩阵和质量矩阵，其中的刚度矩阵与静力分析中的刚度矩阵相同，这里仅给出进行模态分析所需的质量矩阵等其他条件。 10.2.1 桁架单元质量矩阵 集中质量矩阵（局部坐标） 设杆材料密度为 ，单元长度为 ，截面积为 ，每个节点分担单元1/2质量，则单元质量矩阵 (10.2‑1) 整体坐标下的质量矩阵转换 整体坐标下单元一致质量矩阵为 …
(业余时间 和小伙伴一起写写有限元程序 ） 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 4节点四面体线性“完全积分”单元，用Abaqus中的命名规则，即 C3D4。 问题描述( Problem Description) 悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，受重力作用。材料弹性模量E=200000MPa，材料泊松比为0.3。 基于MATLAB编程实现该悬臂梁的弹性静力分析，采用C3D4单元进行模拟，并将MATLAB的计算结果与Abaqus、midas Gen分析结果进行对比。 MATLAB 编程 model disp stress …
程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所； 2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 平面6节点二次“完全积分”单元，用Abaqus中的命名规则，即CPS6。 问题描述( Problem Description) XZ平面内的悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa，材料泊松比为0.3。 MATLAB编程实现该悬臂梁的弹性静力分析，采用CPS6单元进行模拟，并将MATLAB的计算结果与Abaqus分析结果进行对比。 MATLAB 编程 disp stress Abaqus Model stress 单元应力对比 (MATLAB …
(业余时间 和小伙伴一起写写有限元程序 ） 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 8节点六面体线性“完全积分”单元，用Abaqus中的命名规则，即 C3D8。 问题描述( Problem Description) 悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，受重力作用。材料弹性模量E=200000MPa，材料泊松比为0.3。 基于MATLAB编程实现该悬臂梁的弹性静力分析，采用C3D8单元进行模拟，并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程 disp stress SAP2000 …
(业余时间 和小伙伴一起写写有限元程序 ） 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所； 2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 平面4节点线性“完全积分”单元，用Abaqus中的命名规则，即CPS4。 问题描述( Problem Description) XZ平面内的悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa，材料泊松比为0.3。 MATLAB编程实现该悬臂梁的弹性静力分析，采用Q4平面单元进行模拟，并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程 disp stress SAP2000 …
（ 业余时间 和小伙伴一起写写有限元程序 ） 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 三节点三角形单元有三个节点和三个直边，单元内应力为常数，因此称为常应变三角形单元，即Constant Stress Triangle Element（常应力三角形单元），简称CST单元。 问题描述( Problem Description) 平面内的悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa，材料泊松比为0.3。 基于MATLAB编程实现该悬臂梁的弹性静力分析，采用CST平面单元进行模拟，并将MATLAB的计算结果与SAP2000、Midas Gen分析结果进行对比。 MATLAB 编程 mesh …
（ 有空和小伙伴一起写写有限元程序 ） 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 欧拉梁单元基于一定的假设（Kirchhoff假设），在梁的高度远小于其跨度的时候，可以忽略梁的横向剪切变形，此时采用欧拉梁单元进行模拟，能够得到较为满意的结果。但对于跨高比较小的深梁，梁的剪切变形将引起附加挠度，使得原来垂直于轴线的截面在变形后将不再与轴线垂直，且发生翘曲。此时需采用能够考虑横向剪切变形的梁单元进行模拟。 考虑剪切变形修正的经典梁单元和Timoshenko梁单元是两种较为常用的能够考虑梁剪切变形的梁单元，但这两种梁单元仍假定原来垂直于中面的截面在梁变形后仍保持为平面。 在列式方面，Timoshenko梁的基本特点是将挠度和截面转角分别插值。Timoshenko梁单元应用广泛，也容易据此推广到板壳单元 问题描述( Problem Description) 一榀XZ平面内的刚架结构，结构几何信息如图所示；节点1、5处为固定支座，节点4处受到+x方向P=200kN的集中力作用。结构中各杆件采用相同的材料，弹性模量E=30000MPa，梁、柱截面面积分别为0.08m2和0.16m2，梁、柱截面惯性矩分别为0.0128/12m4和0.0256/12m4，梁柱抗剪面积分别为0.0667m2和0.1333m2，材料泊松比为0.2。 基于MATLAB编程实现该框架结构的弹性静力分析，所有构件采用剪切修正梁单元进行模拟，并将基于MATLAB编程计算的结果与Abaqus分析结果进行对比。 MATLAB 编程 Abaqus 位移结果对比： 支座反力对比： 可以看出，MATLAB编程计算结果是和Abaqus计算结果是完全一致的。 …
（和小伙伴，一起给大家带来有限元编程案例。） 程序作者 ( Author) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 问题描述( Problem Description) 节点1处为固定铰支座，节点4处为滑动铰支座，节点5、6、7处分别受到-y方向P=100000N的集中力作用；结构中各杆件采用相同的规格，其中弹性模量E=200000MPa，截面积A=4532mm2。 基于MATLAB编程实现该桁架结构的弹性静力分析，并将MATLAB的计算结果与SAP2000、Midas Gen分析结果进行对比。 MATLAB 编程 SAP2000 Midas Gen ABAQUS 可以看出，MATLAB编程计算结果是和 SAP2000,MidasGen 及 Abaqus计算结果是完全一致的。 注释 ( Comments ) …
（ 有空和小伙伴一起写写有限元程序 ） 程序作者 ( Author) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院 问题描述( Problem Description) 节点1、5处为固定支座，节点4处受到+x方向P=200kN的集中力作用。结构中各杆件采用相同的材料，弹性模量E=30000MPa，梁、柱截面面积分别为0.08m2和0.16m2，梁、柱截面惯性矩分别为0.0128/12 m4和0.0256/12 m4。 基于MATLAB编程实现该框架结构的弹性静力分析，所有构件采用欧拉梁单元进行模拟，并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程结果 SAP2000 计算结果 Midas Gen计算结果 注释 ( Comments ) 源代码已经收录到新书《有限单元法：编程与软件应用》，详见链接：http://www.jdcui.com/?page_id=9731 ( 如果您发现有错误，欢迎批评指正。邮箱：jidong_cui@163.com …
Dynamic Analysis of SDOF System by MATLAB [MATLAB单自由度系统动力响应分析]
利用MATLAB提供的API，将EXCEL数据转换为MATLAB的.mat文件
利用 XTRACT 和 MATLAB 制作 PMM 3D屈服面
3D Truss by MATLAB用MATLAB实现3D桁架静力分析。
MATLAB 2D Truss 用MATLAB实现2D桁架
[01] Ground Motion Selection (选波) 服务
[02] 著:《PERFORM-3D原理与实例》
[03] 著:《有限单元法-编程与软件应用》
[04] 著:《结构地震反应分析-编程与软件应用》
[05] 著:《有限单元法 Python编程》（New!!!）
[06] Software Notes [软件笔记汇总]
[07] 土木工程试验数据处理软件汇总（New!!!）
[08] 自编程序 [Software Box]（New!!!）
[09] 手绘大样 [Detail Drawing]（New!!!）
