- 博客(20)
- 资源 (2)
- 问答 (3)
- 收藏
- 关注
RSS订阅原创 有限元中的高斯积分
对于复杂的积分怎么计算呢？这就提出了数值积分，他就是通过将连续区间进行离散，对离散量求和便是积分的近似值。1D积分 ∫fdξ 在-1到1上求积分，数值积分就是，将-1到1离散为n个点，ΣAk*f(ξk) 其中ξk为-1，-1+2/n，-1+4/n，....，1。Ak为对应的ξk在每个点处的权系数。=========================================高斯积分：...
2020-03-31 23:50:56
3591
原创 高等数学中的积分
在说高斯积分之前，先说一下自己对积分的理解。1. 1D积分∫dx 结果就是一个长度∫f(x)dx f可以看作高度 则结果就是一个面积∫f(x)dx f为线密度 则结果就是一根铁丝的质量(长度方向远大于其余两个方向)注：不一定是铁丝，只要物体的长度方向远大于其余两个方向即可，比如杆。2. 2D积分∫∫dxdy 结果就是一个面积∫∫f(x,y)dxdy f可看作高度 则结果...
2020-03-31 22:36:50
1664
2
原创 CFG and D-U Paths
最近在写软件源代码课程的作业，作业要求：输入任意一段c语言，输出CFG和D-U Paths。下面说一下概念。变量的定义与使用(du of vars)：如 i = i+1 先对 i 使用再对 i 定义。变量的使用：变量被定义以后可分为两种用途，一是用于判断(predicate, p-use) 如x>10，二是用于计算(computation, c-use)如y=x+1。定义清晰路径：(...
2020-03-30 20:20:40
534
1
原创 使用 graphviz dot 绘制有向图
1.下载，解压2.将dot添加到环境变量3.编写dot文件4.将dot文件转化为PNG打开cmdcd到dot文件的目录下dot -Tpng example1.dot -o example1.png5.附代码digraph G{node[color="red"];A[label="傅里叶变换",fontname="FangSong"];node [color="yello...
2020-03-22 13:57:21
1507
原创 信号的概率密度曲线
刚刚研究了一下信号的概率密度曲线，还挺有趣的。所谓概率密度曲线，横轴为信号幅值，纵轴为幅值出现的频率。只要次数足够大，频率可以表示概率。1.构造曲线F1 与 F2 和 F3F1 = np.sin(w*t+phi1) 简谐信号F2 = F1+np.sin(2*w*t+phi2) 一倍频与二倍频信号的叠加F3 = np.exp(-0.05*t)*F1 衰减的简谐信号2. 概率密度分...
2020-03-21 20:18:20
3459
原创 欧拉角(转子动力学)
1.提出目的？确定绕定点转动时刚体的位置2.第一类欧拉角如何得到？xyz坐标系xy平面绕z轴转动ψ 得到x1y1z坐标系 ψ称为进动角/旋进角y1z平面绕x1轴转动θ 得到x1y2z1坐标系 θ称为方位角/挠曲角/章动角x1y2平面绕z1轴转动φ 得到x2y3z1坐标系 φ称为自转角3.第二类欧拉角/广义欧拉角/卡尔丹角 如何得到？xyz坐标系xz平面绕y轴转动α...
2020-03-20 19:38:03
2207
原创 子模型分析
1.创建模型并分析求解。相较于划分更细的子模型，可以称原模型为粗糙模型。保存coarse.db。求解得到coarse.rst。==========================================================/COM,ANSYS MEDIA REL. 150 (11/8/2013) REF. VERIF. MANUAL: REL. 150/VERIFY,V...
2020-03-17 23:53:27
1146
1
原创 如何添加环境变量
我之前下载了mysql，但是每次启动mysql只能先cd到bin文件夹下面，很不方便。还有就是如果windows如果已经安装python3，想要安装python2，为了不引起冲突，安装时先暂时不加入环境变量，之后自己手动加入。=============mysql========================bin文件夹路径 D:\MySQL\mysql-8.0.18-winx64\bin...
2020-03-16 21:07:12
1399
原创 用三角形单元划分平面结构后的求解
上次我写过一篇，求三角形单元划分平面梁后的单元矩阵，今天写的包括求解单元刚度矩阵、总体刚度矩阵、单元弹性系数矩阵。同时列出方程组，带入已知量，求解未知支反力与位移。这个程序用于求解的问题需满足：1.满足下面的四种假设2.平面问题附代码：# -*- coding: utf-8 -*-"""Created on Tue Mar 10 18:49:35 2020@auth...
2020-03-15 13:16:27
882
原创 exp(j*t) = cost+j*sint
1.exp(jt) = cost+j*sint将sin画到纵轴 cos画到横轴上将sin画到x轴上，cos画到y轴上，t画到z轴上振动力学中，转轴上一点处振动位移Asinwt，与其方向垂直的一点处振动位移Bcoswt(当然两点在同一半径处)（而且因为无阻尼 所以垂直方向 相位差为90度），然后用Bcoswt+j*Asinwt就可以表示整个转轴上的受力。实数表示的一个平面，虚数表示的...
2020-03-15 12:41:09
1773
原创 numpy.arange 与 numpy.linspace
自己在计算机中画的函数都是离散的，都是经过采样的。采样间隔为dt 采样频率为fs=1/dtt = np.linspace(0,t总,t总/dt,endpoint=False)t = np.arange(0,t总,dt) （当然也表明了两个函数怎么使用）linspace 开始 结束 个数 生成的序列默认包含结束，如果endpoint=False则不包含arange 开始 结...
2020-03-12 21:54:14
720
原创 非周期信号 进行傅里叶变换
概念：非周期序列包括：准周期与瞬态信号。前者是指 组成信号的正弦函数频率之比不是有理数。后者指 一段时间内信号的频率，幅值是变化的。==============准周期信号======================# 采样间隔dt = 0.001# 采样时间t1与t2t1 = np.linspace(0,100,int(100/dt)) #因为100/dt是一个浮点数 需要int转...
2020-03-12 14:09:02
6282
原创 整数个周期与非整数个周期信号 进行傅里叶变换的差异
#采样间隔 dt = 0.001#采样时间t1 0到100s 采样时间t2 1到100st1 = np.linspace(0,100,int(100/dt)) #因为100/dt是一个浮点数 需要int转化一下t2 = np.linspace(1,100,int(100/dt))=================================================...
2020-03-12 12:29:58
1550
原创 刚度与强度
1.四大强度理论所谓强度，就是指抵抗断裂与屈曲的能力。四大强度理论是用来解释强度失效的，1和2解释断裂，3和4解释屈曲。四大强度理论：1最大拉应力理论2最大伸长线应变理论3最大切应力理论4畸变能密度理论/von mises理论2.质量与转动惯量，转动惯量与面积惯性矩面积惯性矩：截面抵抗弯曲的性质质量惯性矩/转动惯量：转动惯性的度量质量：平动惯性的度量3.刚度所谓刚度...
2020-03-10 13:34:40
3731
原创 异形梁的子结构分析
1.子结构分析普通有限元模型是由很多小单元构成的。而子结构分析，就是将这些小单元组装成一个超级单元，然后用超级单元进行其他分析，诸如静力分析，动力分析等等。步骤：建模定义主自由度生成超级单元(子结构分析)使用超级单元(进行其他分析)扩展超级单元后处理2.下面是一个异形梁进行静力分析的例子finish/clear/filname,beam !文件名/prep7...
2020-03-10 13:14:16
898
1
原创 模态分析与屈曲分析
1.两者都是求解广义特征方程：KΦ=λMΦ 模态分析特征方程KΦ=λSΦ 屈曲分析特征方程其中S为应力刚度矩阵也称几何刚度矩阵,与载荷有关。K为刚度矩阵，与结构自身属性有关。2.模态分析的方法(一般屈曲分析也适用)2.1缩减法将结构的自由度定义为主自由度与从自由度，用主自由度来表达原方程。主自由度至少为感兴趣频率的两倍。现在ansys的版本中已没有该方法，已用子结构分析...
2020-03-10 13:06:16
4523
原创 NewMark法求解动力学响应
问题描述：两个质量块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 ...
2020-03-07 18:40:18
11298
1
原创 N-P求解非线性方程组
对于那些看别人博客，自己却从来不写的人，我希望你们也多写。对于那些只会说别人:''你写的这个别人已经写过了''的人，请你躝开。==========================================用牛顿迭代法(Newton-Raphson)求解非线性方程组x1**2-10*x1+x2**2+8=0x1*x2**2+x1-10*x2+8=0图1是随着求解次数的x1与x...
2020-03-07 18:22:40
465
原创 用三角形单元划分平面梁后的单元矩阵
这是明天的作业 我想着发出来算了问题描述：平面梁 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) 通过这个公式便可计...
2020-03-03 19:32:39
1127
原创 APDL制作动画
针对前面的例子 《坎贝尔图分析》==================================================1.在求解完之后，查看Results Summary可看出：施加了9次转速，每次施加求解一次，共9个载荷步。扩展了30阶模态，也就是每个载荷步内有30个子步。最后一个载荷步便是加载完全的，他的30个子步便是自己要的30阶振型。2.FREQ(Damped...
2020-03-01 20:21:47
2874
masm link debug dosbox
2018-06-03
深度学习英伟达显卡配置
2023-09-07
问一道关于信号处理的问题，帮我解答一下？bai'tuo
2019-10-27
pyinstaller打包后exe不能运行
2019-10-06
不知道如何通过浏览器访问ubuntu sever
2017-08-16
TA创建的收藏夹 TA关注的收藏夹
TA关注的人