matlab简易编程,MATLAB简单编程
本帖最后由 wanggh 于 2016-12-5 14:47 编辑
用数值差分、SOR迭代法求雷诺方程和用牛顿迭代法求解轴向柱塞泵滑靴副压力场的算法,
%油膜厚度场、压力场迭代 MATLAB只认弧度制,不认角度制
clear all; % (60rpm=1rad/s)
wg=1000; %1000rpm缸体公转速度; 单位: 转每分钟 rpm
wz=60; %ws为滑靴自转速度 单位: 转每分钟 rpm
RCP=45*10^(-3); %柱塞分布圆半径,单位: m
pc=24000000; %柱塞腔压力,单位:Pa(单位不能变,就是Pa) Pa
aa=0.87; %压降系数,和阻尼孔直径、长度等有关,无量纲
pr=pc/aa; %支承室压力 单位: Pa
NIANDU=0.04; %介质动压粘度, 单位: Pa*S
beita=18*pi/180; %柱塞泵斜盘倾角 单位:rad 弧度制
nx=100; %沿x方向上区间个数
ny=100; %沿y方向上区间个数
fai=60*pi/180; %滑靴转动角度,(如果可以,fai为向量)单位: rad
R0=15.4*10^(-3); %滑靴外径 单位: m
r0=10.1*10^(-3); %支承室半径 单位 : m
TT=100; %时间
dt=TT/(2*nx); %时间间隔
deltax=2*pi/nx; %x方向上区间长度 (极角),单位: rad
deltay=(R0-r0)/ny; %y方向上区间长度,半径 单位: m
ms=50*10^(-3); %滑靴质量, 单位: kg
mh=56*10^(-3); %柱塞质量 单位: kg
lg=8*10^(-3); %滑靴质心与球头中心距离, 单位: m
ls=10*10^(-3); %滑靴底面到球头的距离, 单位: m
RP=12.4*10^(-3); %柱塞半径 单位: m
RFF=10.5; %柱塞孔的等效液阻,单位:Pa*s/mm^3(这个数值没用上额)
rou=RCP*sqrt(1+(tan(beita))^2*(cos(fai))^2); %公式(数组) m
vs=wg*rou*2*pi/60; %滑靴速度(数组) m/s
f0=0;
MOX=0;
MOY=0; %三个函数初始化
P=0.75*pr*ones(nx+1,ny+1); %压力矩阵,初始值 Pa
P2=ones(nx+1,ny+1); %压力矩阵,初始值
P3=ones(nx+1,ny+1); %压力矩阵,初始值
H=ones(nx+1,ny+1); %油膜迭代初始矩阵 油膜厚度
HL=ones(nx+1,ny+1); %油膜厚度对时间求导矩阵
W=ones(nx+1,ny+1); %系数初始矩阵,西
E=ones(nx+1,ny+1); %系数初始矩阵,东
S=ones(nx+1,ny+1); %系数初始矩阵,南
N=ones(nx+1,ny+1); %系数初始矩阵,北
G=ones(nx+1,ny+1); %系数初始矩阵,东西南北四个系数之和
B=ones(nx+1,ny+1); %原雷诺方程右边常系数,,矩阵初始化
FO1=pi*RP^2*pc; %柱塞腔作用力 N=m^2*Pa
FO2=(ms+mh)*(wg/60)^2*RCP*cos(fai)*tan(beita); %惯性力 (wg/60)为角速度 N
MX=ms*(wg/60)^2*lg*RCP*cos(beita); %滑靴油膜的离心力矩 N*m
con=0; %第二次迭代次数初值
count=0; %压力迭代次数初始值(第三层)
coo=0; %第一层迭代次数初始值
h0=10^(-3)*[0.008 0.01 0.009]' ; %油膜初值 单位: m
hl=10^(-3)*[0.001 0.002 0.003]'; %油膜变化率初值
2016-11-23 21:56 上传
2016-11-23 21:56 上传
matlab简易编程,MATLAB简单编程相关推荐
- matlab和C/C++混合编程--Mex
最近的项目需要matlab和C的混合编程,经过一番努力终于完成了项目要解决的问题.现在就将Mex的一些经验总结一下,当然只是刚刚开始,以后随着学习的深入继续添加.首先讲讲写Mex的一些常规规定,然后我 ...
- Matlab与C/C++混合编程接口及应用
http://www.cnblogs.com/lidabo/archive/2012/08/24/2654148.html 在参考文献基础上.补充和完善了. Matlab与C/C++混合编程接口及应用 ...
- matlab充分利用性能,Matlab高性能编程——代码优化和并行计算
Jeremy Lin @HQU Update: 2014/4/29 Matlab代码优化 Matlab是一种高级计算机语言,同时也是一个用于算法开发,数据可视化,数据分析和数值计算的交互式工作环境. ...
- Matlab与C/C++混合编程 (基于Opencv库)
之前用过基于VS2018 与MATLAB2018a 混合编程(C++特性)(见https://blog.csdn.net/wwwoowww/article/details/83013801),奈何后来 ...
- matlab usb cam,Matlab环境下使用USB2.0摄像头编程image acquisition tools(读书笔记) | 学步园...
(1)Matlab环境下使用USB2.0摄像头编程 Matlab中的图像获取工具箱给我们提供了必要的函数,我们直接调用就可以了,主要就是简单的介绍如何使用该工具箱进行对USB2.0摄像头的编程. TA ...
- matlab编程范例_编程范例到底是什么?
matlab编程范例 Any fool can write code that a computer can understand. Good programmers write code that ...
- matlab与c 接口与混合编程,Matlab与C/C++混合编程接口及应用方法解析
1 引言 Matlab 是当前应用最为广泛的数学软件,具有强大的数值计算.数据分析处理.系统 分析.图形显示甚至符号运算等功能[1].利用这一完整的数学平台,用户可以快速实现十分 复杂的功能,极大地提 ...
- matlab和C/C++混合编程--Mex (转载)
matlab和C/C++混合编程--Mex 最近的项目需要matlab和C的混合编程,经过一番努力终于完成了项目要解决的问题.现在就将Mex的一些经验总结一下,当然只是刚刚开始,以后随着学习的深入继续 ...
- MATLAB和C语言混合编程示例
MATLAB虽然很好用,尤其在数值运算,图像处理,仿真方面很简洁.方便,但有时间当处理比较大的数据时速度确实有点让人捉急,不过MATLAB提供了和C/C++混合编程的接口,以此兼具二者的优点来缓解这个 ...
- matlab和C/C++混合编程--Mex (六)参数传递
最近的项目需要matlab和C的混合编程,经过一番努力终于完成了项目要解决的问题.现在就将Mex的一些经验总结一下,当然只是刚刚开始,以后随着学习的深入继续添加.首先讲讲写Mex的一些常规规定,然后我 ...
最新文章
- Ubuntu 10.10 安装 libx11-dev
- Linux下C高手成长过程----经典书籍推荐
- python前端接口_Python接口自动化——Web接口
- php实现页面雪花效果,CSS3如何实现雪花飘落的效果
- 戴AirPods pro走路会传来咚咚的脚步声,有什么办法能消除这个声音?
- Hyper-V实战:高可用性-配置Hyper-V群集
- SHELL脚本“syntax error: unexpected end of file”解决方案[半原创]
- 第十一次立会暨第四次集体开发
- vue 项目中使用阿里巴巴矢量图标库iconfont
- [VC] 【Visual Studio】2005~2015中文完整旗舰版(附序列号)
- 2023系统分析师综合知识必备知识点
- OpenHarmony hdc_std无反应的解决办法
- Cannot find current proxy: Set 'exposeProxy' property on Advised to 'true' to 以及Spring事务失效的原因和解决方案
- POJ:1182 食物链(带权并查集)
- java native函数库_Java 层调用 Native 层函数的两种方式
- conda 使用清华大学开源软件镜像
- 【C和指针】const指针
- npm--踩坑--npm audit fix 解决方法
- 【字符串】字符串长度与字节长度
- 通过root用户与hadoop用户分别启动zookeeper,两种启动结果为什么会截然不同?