[zx145833]

Hello, I'm currently simulate the oscillation laser welding process. The waveform of oscillation path in y direction is triangle wave (From zero to amplitude value A, then from A to -A, eventually from -A to zero). The cycle duty is 1/f_osc (f_osc is the oscillation frequency). The program of oscillation path in y direction is as below. Simulation results shows the position in y direction at the initial time-step is -A mm, while 0 mm for the expected result. Anyone who encounters the similar problem. Could you give me some advice? Thanks for your help.

Best wishes, Xu Zhao.

// The program

time=RP_Get_Real("flow-time");
time1=((time*f_osc)-int(time*f_osc))/f_osc; //Time normalization in the
interval [0,T]
if(time1<1/(4*f_osc)) // As aforementioned, the position in y
direction can be regarded as a three-
stage linear functions
{
y=xc[1]-4*A*f_osc*time1; // A is the oscillation amplitude
}
else if(time1>=1/(4*f_osc)&&time1<3/(4*f_osc))
{
y=xc[1]-A+4*A*f_osc*(time1-1/(4*f_osc));
}
else
{
y=xc[1]+A-4*A*f_osc*(time1-3/(4*f_osc));
}

[pakk]

You don't show how you define f_osc, there might be a problem there.
If you do one of the following, you will have a problem with integer division that gives the result that you see :

int f_osc=10;
#define f_osc 10

This problem is avoided with these definitions:

float f_osc=10;
#define f_osc 10.0

[zx145833]

Originally Posted by pakk

You don't show how you define f_osc, there might be a problem there.
If you do one of the following, you will have a problem with integer division that gives the result that you see :

int f_osc=10;
#define f_osc 10

This problem is avoided with these definitions:

float f_osc=10;
#define f_osc 10.0

Thank you for the reply. I didn't pay attention to the details and would correct it and finally the testing is successful. The habit of making program is indeed important. Thanks again.

Xu Zhao