Summary Page of Programming Training [编程训练-汇总页]

将一些编程训练的资料汇总到【编程】->【编程训练】目录,可以通过路径【编程】->【编程训练】进行访问,如下图所示。 也可以通过连接直接访问:http://www.jdcui.com/?page_id=13144 目前编程训练放了以前训练的一些笔记。以后有空做编程训练都汇总到这个页面。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号  

[编程笔记] 关于字符编码(Character encoding)

实干、实践、积累、思考、创新。 简单的笔记。 ASCII:最早的字符编码,1个字节,最早是用7个二进制位(128个字符),后面又扩展到8个二进制位(256字符,扩展ASCII),主要用来表示英文字母。 Unicode: 由于ASCII编码不够用,不适合考虑各国语言,于是人们就把全世界的语言字符都放进来,通盘考虑,形成Unicode编码 (统一码、万国码、单一码),解决传统的字符编码方案的局限而产生的,它为每种语言中的每个字符设定了统一并且唯一的二进制编码,以满足跨语言、跨平台进行文本转换、处理的要求。 Unicode字符集又有多种不同的编码方式,如常见的UTF-8,UTF-16,UTF-32等等。 其中UTF-8是互联网上使用最广的一种Unicode实现方法。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[FEM][MATLAB][有限元][编程] 压杆稳定问题MATLAB有限元编程 (《有限单元法-编程与软件应用》章节节选)

实干、实践、积累、思考、创新。 接着博文《[力学][有限元][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 : …

[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础] (《有限单元法-编程与软件应用》章节节选)

实干、实践、积累、思考、创新。 《有限单元法:编程与软件应用》 曲屈分析 章节的内容节选: 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) 公式中               …

[FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Frame Elements) (框架单元模态分析编程)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 接着博文《[动力学] 模态分析基础 (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 …

[FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Truss Element) (桁架单元模态分析编程)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 接着博文《[动力学] 模态分析基础 (Basics of modal analysis)》继续介绍桁架单元用于模态分析的方法。该部分内容也是 书本 《有限单元法:编程与软件应用》模态分析章节的部分内容节选。 10.2 2D桁架结构模态分析 模态分析需要的参数主要有结构的刚度矩阵和质量矩阵,其中的刚度矩阵与静力分析中的刚度矩阵相同,这里仅给出进行模态分析所需的质量矩阵等其他条件。 10.2.1 桁架单元质量矩阵 集中质量矩阵(局部坐标) 设杆材料密度为 ,单元长度为 ,截面积为 ,每个节点分担单元1/2质量,则单元质量矩阵                 (10.2‑1) 整体坐标下的质量矩阵转换 整体坐标下单元一致质量矩阵为    …

[编程][算法][游戏] 编程训练: 3 5 7 游戏开发

程序图标 ( Program Icon )      程序介绍 ( Program Introduction) 3 5 7扑克牌游戏,基本规则是,每次从一行中抽取任意张牌,最后一张牌,谁抽到谁输。 前一排年会上,容总发起了 3 5 7游戏挑战,我连续下了两盘,都输了,所以回去研究了一下,游戏十分有意思,也是很有规律的,并编写了这个小程序。现在看这个游戏,还是很简单的。感兴趣的小伙伴可以试试。 程序界面 ( Program Interface ) 下载 ( Download ) GAME 357: Game357_byCJD_20180929.7z 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[VBA][编程] VBA使用“TextFrame2”设置坐标轴属性出错

用VBA开发EXCEL,为修改坐标轴刻度线字体的大写,通过录制宏看到如下代码: ActiveChart.Axes(xlValue).Select With Selection.Format.TextFrame2.TextRange.Font   .Size = 12   .Bold = True End With 可以从 宏 中看见,代码是通过 TextFrame2 属性对刻度线姿态进行控制的。 但将上述代码在程序中编译,则发现报错,提示  TextFrame2 无法识别,不存在。 不知道是Bug还是什么,最后通过 TickLabels 属性解决。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[FEM][有限元][编程][Matlab][Code by myself] 4节点四面体单元 (C3D4)

  (业余时间 和小伙伴一起写写有限元程序 ) 程序作者 ( 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 …

[FEM][有限元][编程][Matlab][Code by myself] 平面6节点二次“完全积分”单元(CPS6)

程序作者 ( 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 …