实干、实践、积累、思考、创新！ Tag: 结构工程博士 结构工程师 伪程序员 建筑结构 抗震设计 振动控制 弹塑性 编程 软件开发 有限元 超高层 超限设计
实干、实践、积累、思考、创新。 # 傅里叶变换公式的几种形式 研究振动控制、随机振动等，离不开傅里叶变换（Fourier transforms）。对于傅里叶变换，不同文献、书籍有时候会采用不同形式的公式，刚开始看的时候有点凌乱，后面才理清楚，其实不同公式形式本质上都是等价。为了便于后续学习，以下总结几种常见的傅里叶变换公式形式。PS.这里傅里叶变换，指的是非周期函数（周期可以理解为无限长）的傅里叶变换，不是傅里叶级数（Fourier series），傅里叶级数是针对周期函数的说法。 ## 第一种形式 一种通用的形式，\(\Phi \left( p \right)\)和\(F\left( x \right)\)互为傅里叶变换对： $$\Phi \left( p \right) \Leftrightarrow F\left( x \right)$$ $$\Phi \left( w \right) = {1 \over {2\pi …
实干、实践、积累、思考、创新。 这几天更新了FOUR_TRAN （ [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] ），顺便做个风振响应的频谱分析例子。 下图所示，为某超高层整体分析得到的结构的X向顶点风振加速度时程，将其导入FOUR_TRAN，点击Run Fourier Analysis，可在第二个图中看出频谱。 由频谱图可看出，X向风振响应主要受结构的X向平动振型控制，其中主要受X向一阶平动振型控制。 进一步可在FOUR_TRAN中将X轴换为周期的形式，可看到各阶控制振型的对应的周期。 同理，将Y向风振响应时程导入FOUR_TRAN，也可看到Y向风振响应主要受结构的Y一阶平动振型控制。 PS. 一般情况，高耸结构或超高层结构，风振响应主要受一阶平动振型控制。 这也是风振响应与地震响应差别的一个地方。 案例 ( Examples ) [00] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] …
实干、实践、积累、思考、创新。 程序图标 ( Icon ) 程序介绍 ( Introduction) 小伙伴建议下编写的一个小工具，其主要功能是求解平均与标准偏差谱曲线。用户可以导入一组曲线，程序自动求解平均曲线与指定倍数标准偏差的偏差曲线。 程序可考虑正态分布和对数正态分布模型进行评估，并可将分析结果保存为文本文件或者直接输出EXCEL图表。 The basic function of the program is to calculate the average and standard deviation curves. The program can consider Normal distribution and …
实干、实践、积累、思考、创新。 2018写的题目，今天填坑吧。。。。 相关博文 ( Related Topics) [01] [编程][算法][几何] 编程训练: 凸包问题 [Convex-Hull Problem] [02] [编程][算法][几何] 编程训练: 离散点的闭合路径 [Closed Path for points] [03] [编程][算法][几何] 编程训练: 3次B样条曲线 …
实干、实践、积累、思考、创新。 2020 年写的题目，现在已经是2022年了…… 这两天补上笔记 案例 ( Examples ) [01] [数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波] [02] …
实干、实践、积累、思考、创新。 不考虑梁本身的伸长、弯曲等自身变形，单纯考虑水平放置的梁跨中发生挠度 的情况下，水平伸缩量Δ的大小。这其实是一个纯数学的推导，已经和物理无关了（PS. 题目说梁其实不对），其实说的是一根不可伸缩的直线，跨中发生挠度 ，变为两根直线后，水平伸缩量是多大。 其中 \(L\)为直线的总长度，\(\Delta \)为水平伸缩量的大小，\(X\)为直线发生倾斜变形后，水平投影长度的一半。 由总长不变，可得水平伸缩量 \(\Delta = L – 2X\) 同时 \(L\)、\(\Delta \) 及\(X\) 之间满足以下三角函数关系 \(2X = 2\sqrt {{{\left( {\frac{L}{2}} \right)}^2} – {\omega ^2}} \) 将上式代入第一个公式，可得 …
实干、实践、积累、思考，创新。 相关博文 ( Related Topics) [01] [编程][算法][几何] 编程训练: 凸包问题 [Convex-Hull Problem] [02] [编程][算法][几何] 编程训练: 离散点的闭合路径 [Closed Path for points] [03] [编程][算法][几何] 编程训练: 3次B样条曲线 [Cubic …
实干、实践、积累、思考、创新。 阿基米德螺线（亦称等速螺线）(Archimedean spiral)，得名于公元前三世纪希腊数学家阿基米德。阿基米德螺线是一个点匀速离开一个固定点的同时又以固定的角速度绕该固定点转动而产生的轨迹。更多的介绍可以参见：阿基米德螺旋线百度百科词条。 简单的说阿基米德螺线可以用极坐标表示为： r = a+b*θ； 其中其中 a 和 b 均为实数，θ为转角。当θ=0时，a为起点到极坐标原点的距离,b为螺旋线每增加单位角度r随之对应增加的数值。 我们可以通过以下公式将极坐标转换为笛卡尔坐标： x = (a+b*θ)*cos(θ)； y = (a+b*θ)*sin(θ)； 于是就可以非常简单画出图形了，废话不多说，我们用EXCEL来作图。 以下是b=0情况下， a=0及a=1000 时的螺旋线: 以下是a=0情况下， b=10及b=15时的螺旋线: 最后我们可以把曲线导入之前的编写的 GAAGM软件 [软件][工具][编程] GAAGM: Graph …
实干、实践、积累、思考，创新。 1976 年 , 美国数学家梅 (May .R ) 在美国 《 自然 》 杂志上发表的题为 “ 具有复杂的动力学的简单模型 ” 文章中指出 , 在生态学中一些非常简单的确定性的数学模型却能产生看似随机的行为。如虫口模型：xn+ 1= μxn ( 1 -x n )，其中 xn 是第n 年的虫口数，xn+ 1 …
实干、实践、积累、思考，创新！ 直接看图。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 最小二乘法曲线拟合( Least square curve fitting )，又叫最小二乘法多项式曲线拟合，根据给定的m个点,并不要求这条曲线精确地经过这些点，而是曲线y=f(x)的近似曲线y= φ(x)。 使得近似曲线与y=f(x)的偏差最小。按偏差平方和最小的原则选取拟合曲线，并且采取二项式方程为拟合曲线的方法,称为最小二乘法。在某些情况有应用。 以下是编程测试的结果。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 一个数学题，红酒杯与白酒杯的一道数学题。从结构求真群听来的。主要是这么一个问题，有一杯白酒和一杯红酒，从白酒杯里蒯一勺白酒倒入红酒杯里，然后再从红酒杯里蒯一勺倒入白酒杯里，问：白酒杯里的红酒多，还是红酒杯里的白酒多? 于是写这个小软件测试一下，结果也十分有趣，看完结果也恍然大悟。 首先，相互蒯一勺后，白酒杯里的红酒和红酒杯里的白酒是一样的多的。 而且，只要不断相互倒酒，那么白酒杯里的红酒和红酒杯里的白酒就会一直是一样多。 其实与每次倒酒的量无关，不管倒几次，不管每次倒的量是否相等，只要最终状态还是两杯一样多，那么就是一样。最终状态控制，以白酒杯为例，出去多少白酒，就会进来多少红酒。 而且，不断相互倒酒，当倒酒次数无限多厚，白酒杯里面的红酒和白酒，红酒杯里面的红酒和白酒就一样多了。 其实也很好理解，相互不断倒酒，无限多次地倒，相当于酒杯倒均匀了。 十分有意思。感兴趣的朋友可以下载来玩玩： Download: 红酒与白酒 White_RedWine.rar 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新 编程训练。有空的时候做些简单项目训练。 …… 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、实践、积累、思考，创新。 编程训练，Delaunay三角剖分算法。初步测算，这一块还有很多问题需要解决。 微信公众号 ( Wechat Subscription)Delaunay三角剖分算 欢迎关注 “结构之旅” 微信公众号
Derivation of the ε-δ identity [ε-δ等式的推导]。from 崔济东，崔济东的博客，www.jdcui.com, CJD, JidongCui
Orthogonality of Trigonometric Functions [三角函数的正交性]
悬链线（Catenary）是一种曲线，它的形状与悬在两端的绳子因均匀引力作用下掉下来之形相似。适当选择坐标系后，悬链线的方程是一个双曲余弦函数。（1） 悬链线的基本概念（2） 悬链线的公式推导（3） 悬链线的参数确定思路及数值求解算法
康威常数 (Conway’s Constant) 是 Look-and-say sequence 相邻两项数字长度的比值的极限，常用希腊字母λ表示，约等于1.303577，由数学家 John Conway 发现。
初等元包自动机；Elementary Cellular Automata
著名的几何概率问题 —— 蒲丰投针问题（Buffon’s Needle problem ），最初由数学家Georges-Louis Leclerc, Comte de Buffon于18世纪提出。问题可表述为：假定长度为L的的针，随机投到画满间距为T的平行线的纸上，求针和平行线相交的概率。同时有趣的是，该概率值和圆周率(PI)有关系，因此，我们可以利用投针试验来计算圆周率(PI)的值。实际上，这个试验有价值的地方在于：我们可以设计一个试验，它的概率与我们感兴趣的一个变量有关，然后通过大量试验来估算这个量，这其实就是目前使用得十分广泛的蒙塔卡罗法（Monte Carlo method）的主要思想。
前两天在网上无意看到了很多关于分形的东西。感觉十分漂亮。分形体现了一种有序和无序的统一。有兴趣的可以上网搜些资料看看，还有很多专门的分形软件、网站、论坛。 最出名的分形我想就是MandelbrotSet和JuliaSet了。因为分形几何学的创始人就是IBM研究室的数学家B.B.Mandelbrot。下面简单介绍一下MandelbrotSet和JuliaSet的算法。其实算法原理还是比较简单。 MandelbrotSet由一个二次迭代方程：zn+1 = zn2 +c，z0=c所得到。也就是对所有的复平面上的点c都进行迭代，不发散的点都属于Mandelbrot集的范围。 JuliaSet由一个二次映射：zn+1 = zn2 +c 产生。c为一定值。也就是选定一个常复数，然后对所有的复平面上的点z进行迭代，不发散的点都属于Julia集的范围。 其中z和c都是复数。迭代的次数越多，图形显示就越精细，可以无限细分，同时分形还有自相似性，就是无限放大一个细小的部分，看起来和原来的部分一样。 由于分形的无限性，对图形的显示和算法都是极大的挑战。一般编程的书都会用这个例子来介绍多核计算、并行计算。哈哈。手痒了，还是编个程序来看一下效果吧。截图如下： 1、先是Mandelbrot Set（迭代次数50）： （1）来个黑白调色方案（黑色的就是Mandelbrot集的范围） （2）黑白的调色方案太单调了，来看看彩色的效果如何 2、JuliaSet是什么样的呢，我们接下来看看，这次选择迭代次为500次，从下面的图也可以看到，显示的内容更加细致和漂亮。由于Julia集和选择的复常数C有关，下面选择了几个典型的例子来看看效果。。 （1）c=-0.8+0.17i （2）c=-0.7-0.38i (3) c=0.3 哈哈，真的很漂亮，配色很重要，配色没配好，在网上随便找了几个配色方案。算法也没有仔细去考虑。很多书有专门的分形算法介绍。 小结：数学真的很美，很好玩。
[01] Ground Motion Selection (选波) 服务
[02] 书:《PERFORM-3D原理与实例》
[03] 书:《有限单元法-编程与软件应用》
[04] 书:《结构地震反应分析-编程与软件应用》
[05] 博士论文:《RC梁、柱及剪力墙变形性能指标限值研究与试验验证》（Ph.D. Paper）
[06] Software Notes [软件笔记汇总]
[07] 土木工程试验数据处理软件汇总（New!!!）
[08] 自编程序 [Software Box]（New!!!）
[09] 手绘大样 [Detail Drawing]（New!!!）