## General Mechanical

#### Import Bearing coefficients via APDL command

• langea
Subscriber

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 = 3TYPE,2numstr,elem,100001*if,bearing_opt,eq,3,then*do,i,1,n,1n,100000+i,0,5e-3,L1+((L2-L1)/(n-1))*(i-1)!nodes for the bearing statorREAL,i!assigning real constantse,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 frequencybearing_opt = 3*if,bearing_opt,eq,3,then	/input,'ansys_bearing_manual_%current rotational velocity in Hz%','inp','P:Simulation (aktuell)ansysMagnetansys_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

• Sandeep Medikonda
Ansys Employee

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

• langea
Subscriber

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)ansysMagnetansys_bearings_rillen',2,0
*endif```

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

Best regards

• jpasquerell
Ansys Employee

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

• langea
Subscriber

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
Subscriber

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

• peteroznewman
Subscriber

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

• Sankar93
Subscriber

How can i assign radial stiffness in this element

• peteroznewman
Subscriber

K11 and K22 are values for radial stiffness.

• Sankar93
Subscriber

Then k12 and k21 for which stiffness

• peteroznewman
Subscriber

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
Subscriber

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

• peteroznewman
Subscriber

Read the help I pointed out above.