Import Bearing coefficients via APDL command

  • 249 Views
  • Last Post 3 weeks ago
langea posted this 19 September 2018

 

Hello.

I want to analyze rotordynamics with fluid bearings. I have a working APDL script, however I would like to perform the same analysis using Workbench. My problem is as follows:

Since my bearing characteristics are both speed-dependent and length-dependent, I must discretize the whole bearing with 51 bearing elements. While this works neatly in APDL, creating named sections and manually entering all the coefficients is a rather long process.

Is there a way to implement the APDL commands in Workbench?

In /prep7 I define the nodes for the bearings. In /solu I use a loop to import the bearing coefficents corresponding to the respective rotational velocity. So, I think I have to actually insert two APDL command blocks.

I have tried to insert a command block right after I define my connections and insert the commands to create the bearing nodes:

L1 = 7.5e-3 

L2 = 18e-3 + L1

 

ET,2,COMBI214,,0,1  !Unsymmetric bearing in XY-Plane

 

bearing_opt = 3

TYPE,2

numstr,elem,100001

*if,bearing_opt,eq,3,then

*do,i,1,n,1

n,100000+i,0,5e-3,L1+((L2-L1)/(n-1))*(i-1)!nodes for the bearing stator

REAL,i!assigning real constants

e,node(0,0,L1+((L2-L1)/(n-1))*(i-1)),100000+i!assigning bearing elements over 2 nodes each

*enddo

*endif

 In the solution command block I must get the current rotational velocity and then import the corresponding bearing coefficients. The bearing coefficient files are labeled by the rotational velocity in rounds per second. So a rotational velocity of 600 rpm (10 Hz) is labeled as “ansys_bearing_manual_10.

! Get current rotational velocity in rpm, devide by 60 to get rotational frequency

bearing_opt = 3

*if,bearing_opt,eq,3,then

/input,'ansys_bearing_manual_%current rotational velocity in Hz%','inp','P:\Simulation (aktuell)\ansys\Magnet\ansys_bearings_rillen',2,0

cmomega,Rotor,,,f(j)*2*pi

*endif

1.) I am not sure if my approach is correct and possible.

2.) I do not know how to get the rotational velocity.

I am grateful for help of any sort (code, command, example,&hellip. I have attached the APDL Input script, the bearing coefficient files and a dummy Workbench project to implement the commands. If you need further information, please tell me.

Thanks in advance and best regards

Attached Files

Order By: Standard | Newest | Votes
jpasquerell posted this 20 September 2018

 

This is not a documented command so you should independently confirm that it is working for your model. 

*GET,OMEGAX,common,,acelcm,,real,31
*GET,OMEGAY,common,,acelcm,,real,32
*GET,OMEGAZ,common,,acelcm,,real,33

  • Liked by
  • SandeepMedikonda
  • langea
SandeepMedikonda posted this 19 September 2018

Hi Iangea,

  You can write out the database using the CDWRITE command and bring it into Mechanical using External Model. Additionally, if everything you want isn't imported, then you can use command snippets to implement the code you need.

Regards,

Sandeep

  • Liked by
  • langea
langea posted this 20 September 2018

Hi SandeepMedikonda,

Using the External Model does not give me the ability to easily link analysis systems using a GUI. Even though I already looked at an input file for a prestressed modal analysis and used the External model for post-processing in Mechanical Workbench, my goal is to switch to Workbench and implement command snippets.

Regarding code snippets, I think I got the first code snippet to work by leaving out the element type number:

ET,,COMBI214,,0,1  !Unsymmetric bearing in XY-Plane

as it defaults to 1+ current maximum.

For the second code snippet, I thought about getting the current rotational velocity by using the *Get command, but I don't know how exactly.

*GET,frpm,??? ! Get current rotational velocity
fHz = frpm/60
bearing_opt = 3

*if,bearing_opt,eq,3,then
    /input,'ansys_bearing_manual_%fHz%','inp','P:\Simulation (aktuell)\ansys\Magnet\ansys_bearings_rillen',2,0
*endif

How can I get the rotational velocity of the current step?

Best regards

 

langea posted this 21 September 2018

Thanks to both of you. Although the command partially works, I could not get the comment and my inputs to work in a solution loop.

I have decided to stick to APDL as it gives me more flexibility and use the suggested external model method to post process my results in Workbench. Originally, I could not to that, since that procedure apparently does not work with SOLID272 elements. But since I want to incorporate pre-stress effects, I must switch from those elements regardless.

Sankar93 posted this 4 weeks ago

I want to know about the physical meaning of stiffness coefficients in COMBI214 elements. What does k11, k12, k21,k22  stand for?

peteroznewman posted this 4 weeks ago

You can find this diagram in the ANSYS help system for COMBI214.

Sankar93 posted this 4 weeks ago

How can i assign radial stiffness in this element 

peteroznewman posted this 4 weeks ago

K11 and K22 are values for radial stiffness.

Sankar93 posted this 3 weeks ago

Then k12 and k21 for which stiffness

 

peteroznewman posted this 3 weeks ago

Sankar93, did you look in the ANSYS Help system?

Use this URL for help in ANSYS 2019 R2, by following these instructions.

https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v194/ans_thry/thy_el214.html

I copy from that page...

Sankar93 posted this 3 weeks ago

how to assign the axial and torsional stiffness in the bearing connection COMBI214

peteroznewman posted this 3 weeks ago

Read the help I pointed out above.

Close