lijil168的专栏

Excite Power Unit 2020的有限元子结构文件转换只支持Abaqus2018、2019,Excite Power Unit 2016只支持Abaqus6.13、6.14，目前使用Abaqus2017，无法进行转换。

要得到大象的身高的样子，用高h(x,y)作为目标函数是显而易见的，x,y是空间的坐标，即设计变量，平面（x,y）构成设计空间。但是很多时候，评价问题的指标（目标函数）并不是显而易见的，比如，要评估齿轮系统的振动强度，是用齿轮间的啮合力呢，还是齿轮的角加速度，还是所有齿轮的惯量加权角加速度？这需要好好的思考一下。另外，影响目标问题的因素可能有很多，有一些是可控的，有一些是不可控的，有一些对目标问题是敏感的，有一些是不敏感的，有一些影响大些，有一些影响小些，有一些是可观测的，有一些是不可观测的。

【代码】hypermesh二次开发tcl脚本。

发动机悬置的动刚度和工作时的频率、预载荷、动态载荷的幅值有关，要确定准确的动刚度，需要对发动机工作时的悬置动刚度进行测试，对怠速时悬置动刚度，AVL Excite 悬置数学模型有几种，比较简单的一种为：用FTAB模拟静刚度，用SLS模拟动刚度，计算模型如下图，动刚度为复数，有模和相位角（滞后角，与阻尼有关）

1.2、轴瓦节点计算并建立2.4、镜像选择的面

Entities>

2、绘制弹簧2.1、绘制弹簧截面

打开ABAQUS CAE后，rpy脚本日志（abaqus.rpy）如下，一个Viewport对象’Viewport: 1’已经建立。

本文是对《简明微积分》的学习总结，从微分到积分、重积分、线积分、面积分，从Green公式到Guass公式、Stokes公式，从梯度到散度、旋度，这可能是力学中最重要的几个数学概念。从微分到外微分，则更是让人眼界大开。“Stokes公式是微积分的顶峰。从理论上讲，这是微积分的终点，也是微积分从古典走向现代的入口处。”自《简明微积分》。

键合图-多体及多物理系统建模利器1、 键合图的广义状态变量1.1势是流的原因1.2流是势的原因1.3动量是势对时间的累积效应1.4位移是流对时间的累积效应1.5功率是势和流的相干效应1.6动能（磁能）是流对冲量（磁通量）的累积效应1.7势能（电能）是势对位移（电荷）的累积效应2、键合图到方框图的转换3、键合图列写系统状态方程4、Matlab Simulink方框图建模5、AMESIM键合图建模6、Adams微分方程建模1、 键合图的广义状态变量1.1势是流的原因1.2流是势的原因1.3动量是势对时

对系统求微分，传递函数乘s，对系统求积分 ，传递函数乘1/s。

[TOC]能量（Power、Energy、Activity)–从另一个角度看物理)

AMEsim用2d接触模型来模拟同步器的动态行为：

临界速度和临界位移的意义。

齿轮敲击不一定发生在驱动力矩为负时，也不一定发生在齿轮减速时，要具体问题具体分析，否则就容易犯经验主义的错误！

AMEsim齿轮噪声仿真1.模型由AMEsim例子改造而成2、爆压3、曲轴转速和一轴转速，曲轴转速被调制4、传动齿轮转速5、太阳轮、行星架、行星轮转速，太阳轮与行星架固定，行星轮空套在行星架上，转速波动很大6、传动齿轮敲击能量7.行星轮敲击能量1.模型由AMEsim例子改造而成2、爆压3、曲轴转速和一轴转速，曲轴转速被调制4、传动齿轮转速5、太阳轮、行星架、行星轮转速，太阳轮与行星架固定，行星轮空套在行星架上，转速波动很大6、传动齿轮敲击能量7.行星轮敲击能量

【代码】挠性航天器姿态机动动力学模型及PD鲁棒控制。

*当用鼠标选取与虚轴相交的点时，找到系统振荡开始的增益Km=14,对应的极点有三个（红色的点），振荡频率wm=14rad/s，根据整定公式可得Kp=8.8371;基于稳定性分析的频域响应PID整定方法：对于给定的被控对象传递函数，可以得到其根轨迹，对应穿越jw轴的点，增益Km为系统开始振荡时的增益K值，wn为振荡频率。整定前系统根轨迹图见Figure1，系统有3个极点，根轨迹图有3根，从极点开始，终止与无限远（因为系统无0点，否则有终止于0点的线），**整定后的系统根轨迹图见Figure3，

开始接触以前不敢碰的问题——传热及流体力学问题。

有限元学习笔记-虚功平衡方程的线性化及牛顿-拉普森迭代算法。张量微分是解连续介质物理问题的利器。

开始啃最硬的骨头，非线性连续介质力学问题的有限元解法及张量应用。《有限元法》对非线性问题的表示和《张量分析》对张量的表示一样酷。有了张量分析，连续介质物理问题就如鱼得水，否则寸步难行。

python tkinter UI设计

终于见到了《张量分析》中的曲线坐标、雅克比矩阵与雅克比行列式。有限元学习笔记——等参有限单元。

有限元的混合格式，利用虚位移原理和拉格朗日乘子算法推导胡—鹫津泛函。学习有限元法，对物理和数学都是较大的考验，不会物理，看不懂问题，不会数学，无法表示问题，求解问题和解释问题。

有限元法》求表面力产生的节点载荷，公式很简单，但要掌握其物理意义并不轻松，理解了其中含义，就可以举一反三，处理其他形式的载荷。

开始学习《有限元法》。一个简单的力学问题，从微分、变分到虚位移原理建模，从加权余量法、里茨法到伽辽金和有限元算法，问题很小，意义很大。

因为当系统出现故障时，故障通常会对系统信号产生调制作用，使系统的响应信号振幅产生周期性的变化（调幅），或者频率产生周期性的变化（调频），这种周期性变化的频率对应于故障发生的频率，变化的大小对应于故障的强弱，经过本软件希尔伯特变换得到瞬时频率和瞬时振幅，即可找到故障发生的频率和故障强弱，希尔伯特变换在故障诊断中具有广泛的运用。5、将调制后的信号进行FFT变换：频分析->频谱，可看出以30HZ的载波频率为中心，间隔5HZ及10HZ的边频（调制信号的频率），通过边频分析定位故障源也是NVH故障分析的一种方法。

难度升级，要想不迷路，一、掌握公式的内在逻辑和外在形式之间的关系，二、牢记四大算法：“分”“合”“缩”“变”。

张量分析之张量微分学习笔记：

从幂函数、幂级数、复变函数到张量函数。在变化中“看见”不变的东西是认识的重要一步，认识了张量，我们就可以在“平面”世界看到“立体”世界的样子。

终于敲开张量分析之门。从直角坐标系、到直线坐标系、到曲线坐标系，从标量、到矢量、到张量，抽象度越高，威力越大。雅可比矩阵和雅可比行列式的物理意义，了不起的雅可比！用爱因斯坦约定求和进行张量定义和张量算法推导，完美！

对抽象的东西，如果能找到具体的参照物，即使只是特例，也是很有意义的，因为我们可以“看到”以前所看不到的东西。矢量、矩阵、行列式计算的物理意义。

变分法求解两点间直线距离最短，问题很小，意义很大。

动力学系统刚度及模态有效质量的计算公式推导

用Abaqus进行仿真计算的加速度冲击响应的结果与SPECTR部分不一致，尝试找一下原因。结果与Abaqus仿真结果基本一致，加速度响应从14HZ到18Hz中途有下降，而。

norm_u=np.linalg.norm(u,ord=None,axis=0) #按列求模，对单位矢量模=1。

很明显，上部直线的长度增加，这说明1方向的应力，σ11，是拉伸的。类似地，下部直线的长度缩短，说明σ11是压缩的。而如果对载荷产生的位移类型有怀疑，则应采用不同的单元类型。剪力锁闭只影响受弯曲载荷的完全积分线性单元，这些单元的功能在受纵向或剪切荷载时并没有问题。释考虑一个受纯弯的结构中的一小块材料，材料将产生的弯曲如图4－4所示。谓“规则形状”是指单元的边相交成直角，而任何的节点位于边的中点。，l是长度，b是宽度，d是梁的高度。模拟采用了线性或二次的完全积分单元，并说明了。.........

求解的过程我就不展开了，让我们重点关注结果的形式和一些重要的参数，wn是系统的无阻尼固有频率，本例是1Hz，wd是系统的阻尼固有频率,w是激励力的频率，本例是2Hz,zeta是阻尼比。大家好，我今天的培训是“机械系统的模态和振动”，说是培训，其实是谈一谈自己多年对NVH学习的心得体会，谈起NVH，可能很多人都觉得它很“玄”，即使是在NVH领域工作多年的人。测它，这是“测”得的小球的位移图。除了看到小球振幅逐渐减小外，我们还能发现小球的振动频率不变，是1Hz，另外，它的振幅是按指数减小的。......

建立刚体动力学模型，不需要网格，只需在。

假设一自由度系统受到的激励力是均值为0、自谱为S0=1的理想白噪声平稳过程。求系统响应的均方值与自谱（功率谱密度函数，设m=1,k=(2*pi)**2,阻尼系数zeta=0.16,即构造固有频率为1Hz,快速衰减的系统。激励力随机时，和确定时类似，都需要利用系统的脉冲响应函数h(t)或者复频响应函数H(w)。............