应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 机械 » 流体机械 » 小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢
51/1返回列表
查看: 1617  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Dr_Kevin

木虫 (小有名气)

[求助] 小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢 已有2人参与

小弟求滑动轴承非线性油膜力matlab程序，最好是差分法的，谢谢
回复此楼

» 收录本帖的淘帖专辑推荐

摩擦疲劳振动与可靠性 滑动轴承及铸造专题

» 猜你喜欢
1楼 2015-11-02 11:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

熊友平

新虫 (初入文坛)
为什么我运行出错了嘞

发自小木虫Android客户端
赞一下 回复此楼
6楼2018-07-19 23:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

雪之星愿

新虫 (初入文坛)

【答案】应助回帖

% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m;                                    %沿圆周方向均匀划分m格   
delta_n=2/n;                                        %沿轴向均匀划分n格
for j=1:n+1                                         %对各节点赋压力初值为0
    for i=1:m+1
        P(i,j)=0;
        Q(i,j)=0;
    end
end
S=0;
T=0;
for j=2:n                                    %第一次计算各内部节点的压力值
    for i=2:m
       A=(1+e*cos((i+1/2-1)*delta_m))^3;
                B=(1+e*cos((i-1/2-1)*delta_m))^3;
                C=((1/ratio)*delta_m/delta_n)^2*(1+e*cos((i-1)*delta_m))^3;
                D=C;
                E=A+B+C+D;
                F=delta_m*e*(cos((i+1/2-1)*delta_m)-cos((i-1/2-1)*delta_m));
                P(i,j)=(A*P(i+1,j)+B*P(i-1,j)+C*P(i,j+1)+D*P(i,j-1)-F)/E;
                if (P(i,j)<=0)                                   %引入雷诺边界条件   
            for k=i:m
                P(k,j)=0;
                S=S+abs(Q(k,j));
                T=T+abs(P(i,j));
                Q(k,j)=P(k,j);
            end
            break;
                end
                S=S+abs(P(i,j)-Q(i,j));
                T=T+abs(P(i,j));
                Q(i,j)=P(i,j);
    end
end
while(S/T>0.001)                          %循环计算各内部节点的压力值
    S=0;
    T=0;
    for j=2:n
                for i=2:m
                        A=(1+e*cos((i+1/2-1)*delta_m))^3;
                        B=(1+e*cos((i-1/2-1)*delta_m))^3;
                        C=((1/ratio)*delta_m/delta_n)^2*(1+e*cos((i-1)*delta_m))^3;
                        D=C;
                        E=A+B+C+D;
                        F=delta_m*e*(cos((i+1/2-1)*delta_m)-cos((i-1/2-1)*delta_m));
                        P(i,j)=(A*P(i+1,j)+B*P(i-1,j)+C*P(i,j+1)+D*P(i,j-1)-F)/E;
                        if (P(i,j)<=0)
                for k=i:m
                    P(k,j)=0;
                    S=S+abs(Q(k,j));
                    T=T+abs(P(k,j));
                    Q(k,j)=P(k,j);
                end
                break;
                        end
            S=S+abs(P(i,j)-Q(i,j));
            T=T+abs(P(i,j));
            Q(i,j)=P(i,j);
        end
        end
end
%P;                    %此部分是当本函数单独使用时，用于绘制离散压力分布
%for i=1:m+1
%   x(i)=(i-1)*delta_m;
%end
%for j=1:n+1
%   y(j)=(j-1)*delta_n;
%end
%figure
赞一下(2人) 回复此楼
2楼2016-01-14 15:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dr_Kevin

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 雪之星愿 at 2016-01-14 15:28:03
% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m;                                    %沿圆周方向均匀划分m格 ...

谢谢~~只是。。。您这个程序中的F表达式对么。。。
赞一下 回复此楼
3楼2016-01-15 09:54:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dr_Kevin

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 雪之星愿 at 2016-01-14 15:28:03
% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m;                                    %沿圆周方向均匀划分m格 ...

你是把F里做减法的1直接运算掉了~~看明白了~谢谢
赞一下 回复此楼
4楼2016-01-15 10:23:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭