实干、实践、积累、思考、创新！ Tag: 结构工程博士 工程师+程序员 | 结构动力学与抗震 | 复杂结构分析与设计 | 地铁振动减振 | 人行&风致振动控制 | 海上风电CAE研发 | TMD&TLD阻尼器 | 固体力学与有限元 | 程序开发&软件定制
程序作者 ( 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中的命名规则，即CPS8。 问题描述( Problem Description) XZ平面内的悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa，材料泊松比为0.3。 MATLAB编程实现该悬臂梁的弹性静力分析，采用Q4平面单元进行模拟，并将MATLAB的计算结果与Abaqus分析结果进行对比。 MATLAB 编程 disp stress Abaqus Model …
(业余时间 和小伙伴一起写写有限元程序 ） 程序作者 ( 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.华南理工大学建筑设计研究院 基本概念 ( Concept ) 欧拉梁单元基于一定的假设（Kirchhoff假设），在梁的高度远小于其跨度的时候，可以忽略梁的横向剪切变形，此时采用欧拉梁单元进行模拟，能够得到较为满意的结果。但对于跨高比较小的深梁，梁的剪切变形将引起附加挠度，使得原来垂直于轴线的截面在变形后将不再与轴线垂直，且发生翘曲。此时需采用能够考虑横向剪切变形的梁单元进行模拟。 考虑剪切变形修正的经典梁单元和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编程计算的结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程 SAP2000 Midas Gen 注释 ( …
（ 有空和小伙伴一起写写有限元程序 ） 程序作者 ( 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 …
【编程笔记】：异步 I/O 的基本概念与Ajax的异步刷新机制类似，如下图。 简而言之就是，发出请求后无需等待服务器处理结果便可继续进行其他函数调用，发出端无法准确知道服务器返回结果的时间。注重结果，不关心过程。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
申请了GitHub，地址: https://github.com/JdCui 以后在这里分享自己写的一些源代码。
（1）阻塞函数 指的是，函数调用成后、返回结果前（处理完该函数的任务前），当前线程会挂起，直至调用函数结束返回结果，线程才能再次运行。 （2）非阻塞函数 概念与阻塞函数相对，指的是，函数调用后会立即返回（该函数的任务有没有完成没有关系），因此线程不不会挂起和等待，可继续往下运行。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
Python 出现”ValueError:invalid literal for int() with base 10: \xef\xbb\xbf0″ 错误
做个备忘，写程序的时候用到版本号。从这个表也可以看出，1997年颁布的 R14 是一个AutoCAD的一个转折。R14版本以后，AutoCAD功能得到了极大的增强。 Official Name Version DWG Tag Release AutoCAD Version 1.0 1 1 AutoCAD Version 1.2 1.2 2 AutoCAD Version 1.3 1.3 3 AutoCAD Version 1.4 1.4 4 AutoCAD …
Microsoft Office Version Number [微软Office软件的版本编号]
做个记录备忘~~ 🙂 🙂 🙂 🙂 Visual C++ _MSC_VER _MSC_FULL_VER 1.0 800 3.0 900 4.0 1000 4.2 1020 5.0 1100 6.0 1200 6.0 SP6 1200 12008804 7.0 1300 …
利用MATLAB提供的API，将EXCEL数据转换为MATLAB的.mat文件
3D Truss by MATLAB用MATLAB实现3D桁架静力分析。
元包自动机（Cellular Automata）这里就不介绍了，可以在我的其他博文中了解到，有兴趣的朋友也可以去百度一切。这里利用元包自动机进行森林火灾蔓延模拟，当然这肯定不是基于物理的模拟了。最简单的基于元包自动机的森林火灾模型（Cellular Automata Forest Fire Model）大家可以在这个网站 Cellular Automata Forest Fire Model 了解到，作者用Java写了小程序进行模拟。 森林火灾模型介绍（The rules of Cellular Automata Forest Fire Model）： 燃烧的树（红色格子）会在下一步变成空地（白色格子）。A burning tree (red cell) becomes an empty site (white cell)。 如果一颗非燃烧的树（绿色格子）周围有树燃烧，那么它会在下一步变成燃烧的树（白色格子）。A tree (green cell) becomes …
最近在微信上有一个很火的小游戏——Core Ball，又名见缝插针游戏。游戏有很多个不同的版本，不同人编写。游戏玩法很加单，就是用一个轮盘不断的转，轮盘上插满了球，然后点击鼠标，发射球（针头）到轮盘上，但是球不能和轮盘上已有的球发生碰撞，碰撞则重新玩该关卡。游戏难点是：寻找合适的时机发射球，考验你的耐性，考验你是否眼疾手快。趁着休息时间，我也来编写一个：CJD Core Ball。大家快来玩。
生命游戏其实是一个零玩家游戏（百度百科），这个小游戏是一个简单的二维元包自动机（Cellular Automaton）,1970年由英国的数学家John Horton Conway提出，关于生命游戏的介绍参考Wikipedia的Conway’s Game of Life。
[01] 著:《PERFORM-3D原理与实例》
[02] 著:《有限单元法-编程与软件应用》
[03] 著:《结构地震反应分析-编程与软件应用》
[04] 著:《有限单元法 Python编程》
[05] 著:《结构地震动力响应Python编程》
[06] 著:《Grasshopper建筑结构参数化建模应用实例》
[07] 土木工程试验数据处理软件汇总
[08] 自编程序 [Software Box]
