Load Calculation of Hydrodynamic Journal Bearing

18 次查看(过去 30 天)
I have the pressure field generated in circumferencial and axial direction. But can't crack the double integration part. Can anybody help me with a code for that? It would be big help.
  2 个评论
Dyuman Joshi
Dyuman Joshi 2023-12-30
Is theta_prime related to theta? If yes, then how?
Are p and R constants?
Adwait Mahajan
Adwait Mahajan 2023-12-30
编辑:Adwait Mahajan 2023-12-30
R is a constant. p is a function of theta_prime and z.
theta_prime = theta - alpha
where, alpha is the attitude angle of the bearing.

请先登录,再进行评论。

采纳的回答

Hassaan
Hassaan 2023-12-30
编辑:Hassaan 2023-12-30
you could use integral2 or dblquad function:
% Define the pressure function, replace this with your actual function
pressure = @(theta_prime, z) cos(theta_prime); % Example function
% Integration bounds
theta_prime_min = 0;
theta_prime_max = 2 * pi;
z_min = 0;
z_max = 1; % Replace with the actual length L of the bearing
% Radius R is a constant, replace with the actual radius of the bearing
R = 1;
% Radial load Wr calculation using double integration
Wr = integral2(@(theta_prime, z) pressure(theta_prime, z) .* R .* cos(theta_prime), ...
z_min, z_max, theta_prime_min, theta_prime_max);
% Tangential load Wt calculation using double integration
Wt = integral2(@(theta_prime, z) pressure(theta_prime, z) .* R .* sin(theta_prime), ...
z_min, z_max, theta_prime_min, theta_prime_max);
disp(Wr);
4.5699
disp(Wt);
2.2245
Replace the pressure function and the bounds with your actual function and bearing dimensions.
------------------------------------------------------------------------------------------------------------------------------------------------
If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.
  1 个评论
Adwait Mahajan
Adwait Mahajan 2023-12-30
Thanks a lot. I was trying to implement Simpson's integral rule. But this makes my work way easier.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by