my blogs include programming, math, physics.回首过去，放眼未来。

非牛顿流体非定常流动非牛顿流体：粘度与剪切速率有关非定常：参数随时间变化问题连续性方程、运动方程、本构方程、粘度方程（Bird-Carreau模型），运动方程不考虑惯性项运动方程多了一项，表明速度随时间变化求解四边形二次单元划分单元方程M、B、C、F、D为形函数矩阵与形函数矩阵导数构成的量进行积分得到的。D还包括粘度，而粘度需要剪切速率来计算。整体方程也就是，进行迭代，直到x收敛。你也可以计算一下，考虑惯性项的非牛顿流体非定常流动。补充整体刚度方程求解时出现的

一维一次常微分方程维：坐标系维数次：导数阶数方程及边界条件注：边界条件包括场值、场导数值求解用一维线性单元求解，对应固体力学就是杆单元。不同的书里面表达方式不同，插值函数，形函数，插值函数矩阵，形函数矩阵Galerkin有限元法——权函数就是插值函数数值解残差单元内，权函数与残差的内积为0，得到加权余量方程，化简加权余量方程就能得到单元有限元方程。书中都是通过galerkin有限元法推导的。将所有单元的K、b进行集总，得到，u即为节点解。解析解为u=.

杆单元梁单元平面问题单元空间问题单元轴对称问题单元板单元壳单元ansys中的单元单元分类七大类单元：结构单元、热单元、电磁单元、耦合场单元、流体单元、网分单元、显示动力分析单元结构：点单元、线单元、梁单元、实体单元、壳单元、实体壳单元、界面单元、多点约束单元、层合单元热：点单元、线单元、实体单元、电单元、壳单元电磁：plane、solid、hf显示动力单元黏弹实体单元visco耦合场单元：solid、plane、rom接触单元conta、target流体

一般动力学方程 MU'' + CU' + KU = f转子动力学方程 MU'' + (G+C)*U' + (B+K)*U = f （静坐标系下)其中G是陀螺矩阵、B是旋转阻尼矩阵分析步骤：建立模型定义单元类型定义材料定义转速考虑是否引入陀螺效应划分网格求解查看结果...

对于复杂的积分怎么计算呢？这就提出了数值积分，他就是通过将连续区间进行离散，对离散量求和便是积分的近似值。1D积分 ∫fdξ 在-1到1上求积分，数值积分就是，将-1到1离散为n个点，ΣAk*f(ξk) 其中ξk为-1，-1+2/n，-1+4/n，....，1。Ak为对应的ξk在每个点处的权系数。=========================================高斯积分：...

在说高斯积分之前，先说一下自己对积分的理解。1. 1D积分∫dx 结果就是一个长度∫f(x)dx f可以看作高度 则结果就是一个面积∫f(x)dx f为线密度 则结果就是一根铁丝的质量(长度方向远大于其余两个方向)注：不一定是铁丝，只要物体的长度方向远大于其余两个方向即可，比如杆。2. 2D积分∫∫dxdy 结果就是一个面积∫∫f(x,y)dxdy f可看作高度 则结果...

上次我写过一篇，求三角形单元划分平面梁后的单元矩阵，今天写的包括求解单元刚度矩阵、总体刚度矩阵、单元弹性系数矩阵。同时列出方程组，带入已知量，求解未知支反力与位移。这个程序用于求解的问题需满足：1.满足下面的四种假设2.平面问题附代码：# -*- coding: utf-8 -*-"""Created on Tue Mar 10 18:49:35 2020@auth...

1.exp(jt) = cost+j*sint将sin画到纵轴 cos画到横轴上将sin画到x轴上，cos画到y轴上，t画到z轴上振动力学中，转轴上一点处振动位移Asinwt，与其方向垂直的一点处振动位移Bcoswt(当然两点在同一半径处)（而且因为无阻尼 所以垂直方向 相位差为90度），然后用Bcoswt+j*Asinwt就可以表示整个转轴上的受力。实数表示的一个平面，虚数表示的...

1.四大强度理论所谓强度，就是指抵抗断裂与屈曲的能力。四大强度理论是用来解释强度失效的，1和2解释断裂，3和4解释屈曲。四大强度理论：1最大拉应力理论2最大伸长线应变理论3最大切应力理论4畸变能密度理论/von mises理论2.质量与转动惯量，转动惯量与面积惯性矩面积惯性矩：截面抵抗弯曲的性质质量惯性矩/转动惯量：转动惯性的度量质量：平动惯性的度量3.刚度所谓刚度...

1.两者都是求解广义特征方程：KΦ=λMΦ 模态分析特征方程KΦ=λSΦ 屈曲分析特征方程其中S为应力刚度矩阵也称几何刚度矩阵,与载荷有关。K为刚度矩阵，与结构自身属性有关。2.模态分析的方法(一般屈曲分析也适用)2.1缩减法将结构的自由度定义为主自由度与从自由度，用主自由度来表达原方程。主自由度至少为感兴趣频率的两倍。现在ansys的版本中已没有该方法，已用子结构分析...

问题描述：两个质量块m1 m2由三个弹簧连接k1 k2 k3二自由度系统，m1为1kg m2为2kg k1为1000N/m k2为1000N/m k3为1000N/m刚度阵为[[k1+k2,-k2],[-K2,K2+K3]]质量阵用集中质量阵表示。阻尼采用比例阻尼。newmark法 有两个参数β与γ(有的地方称呼不一样)。该方法将连续时间进行离散，从初始条件x0 ...

对于那些看别人博客，自己却从来不写的人，我希望你们也多写。对于那些只会说别人:''你写的这个别人已经写过了''的人，请你躝开。==========================================用牛顿迭代法(Newton-Raphson)求解非线性方程组x1**2-10*x1+x2**2+8=0x1*x2**2+x1-10*x2+8=0图1是随着求解次数的x1与x...

这是明天的作业 我想着发出来算了问题描述：平面梁 t=1 E=1 u=0 正方形边长a=1单元节点编号 按照逆时针的顺序E1 1 3 2E2 2 3 4E3 4 3 5E4 5 6 4节点坐标 x y--------- 1 0 02 0 a3 a 04 a a5 2a 06 2a aKe = t*A*(B.T*D*B) 通过这个公式便可计...

利用ansys导出节点与单元信息 然后输入给自己用c fortran python写的程序进行计算ansys2019的输出命令与之前不一样了，自己看的曾攀老师的那本书给的还是之前的命令。FINISH/CLEAR,START/TITLE,Export Element and Node Info/PREP7ET,1,PLANE42N,1,0,2N,2,0,1N,3,1,1N,4...

循环对称结构分析问题描述finish/clear/filname,modal2/title,modal analysis of a gear/prep7et,1,shell63et,2,solid45mp,ex,1,2e8mp,nuxy,1,0.3mp,dens,1,7.8e-6r1=30r2=r1+10r3=r1+30r4=r1+50r5=r1+70...

问题描述finish/clear/filname,modal/title,modal analysis of a modal airplane wing/prep7et,1,plane42et,2,solid45mp,ex,1,38000mp,dens,1,1.033e-3mp,nuxy,1,0.3k,1 !k,1,0,0,0k,2,2 !k,2,...