I have programmed a heating cooling PID Controller for DL 06 PLC(Koyo / Direct Logic). It is not easy because using ladder logic to make a PID loop. And finally i have finished it.
Due to efficiency we have bought a Direct Logic PLC. DL 06 type is a cheap PLC. It is my first experiences.
I Selected the PID function in directsoft 5.3 and starting entering my addressing.
Loop 1
Table start address V1600
Setpoint Variable V1602
Process Variable V1603
Output V1605
I don't see anything in ladder logic. The PID loop operates outside of the PLC ladder logic. You can use ladder logic to interface with the PID loop. Most setups will use some ladder logic.
I have just copied from http://en.wikipedia.org/wiki/PID_controller. It is very useful for me to know how PID works.
A proportional–integral–derivative controller (PID controller) is a generic control loop feedback mechanism widely used in industrial control systems. A PID controller attempts to correct the error between a measured process variable and a desired setpoint by calculating and then outputting a corrective action that can adjust the process accordingly and rapidly, to keep the error minimal.
[edit] General
The PID controller calculation (algorithm) involves three separate parameters; the proportional, the integral and >derivative values. The proportional value determines the reaction to the current error, the integral value determines the reaction based on the sum of recent errors, and the derivative value determines the reaction based on the rate at which the error has been changing. The weighted sum of these three actions is used to adjust the process via a control element such as the position of a control valve or the power supply of a heating element.
By tuning the three constants in the PID controller algorithm, the controller can provide control action designed for specific process requirements. The response of the controller can be described in terms of the responsiveness of the controller to an error, the degree to which the controller overshoots the setpoint and the degree of system oscillation. Note that the use of the PID algorithm for control does not guarantee optimal control of the system or system stability.
Some applications may require using only one or two modes to provide the appropriate system control. This is achieved by setting the gain of undesired control outputs to zero. A PID controller will be called a PI, PD, P or I controller in the absence of the respective control actions. PI controllers are particularly common, since derivative action is very sensitive to measurement noise, and the absence of an integral value may prevent the system from reaching its target value due to the control action.
Note: Due to the diversity of the field of control theory and application, many naming conventions for the relevant variables are in common use.
[edit] Control loop basics
A familiar example of a control loop is the action taken to keep one's shower water at the ideal temperature, which typically involves the mixing of two process streams, cold and hot water. The person feels the water to estimate its temperature. Based on this measurement they perform a control action: use the cold water tap to adjust the process. The person would repeat this input-output control loop, adjusting the hot water flow until the process temperature stabilized at the desired value.
Feeling the water temperature is taking a measurement of the process value or process variable (PV). The desired temperature is called the setpoint (SP). The output from the controller and input to the process (the tap position) is called the manipulated variable (MV). The difference between the measurement and the setpoint is the error (e), too hot or too cold and by how much.
As a controller, one decides roughly how much to change the tap position (MV) after one determines the temperature (PV), and therefore the error. This first estimate is the equivalent of the proportional action of a PID controller. The integral action of a PID controller can be thought of as gradually adjusting the temperature when it is almost right. Derivative action can be thought of as noticing the water temperature is getting hotter or colder, and how fast, anticipating further change and tempering adjustments for a soft landing at the desired temperature (SP).
Making a change that is too large when the error is small is equivalent to a high gain controller and will lead to overshoot. If the controller were to repeatedly make changes that were too large and repeatedly overshoot the target, the output would oscillate around the setpoint in either a constant, growing, or decaying sinusoid. If the oscillations increase with time then the system is unstable, whereas if they decay the system is stable. If the oscillations remain at a constant magnitude the system is marginally stable. A human would not do this because we are adaptive controllers, learning from the process history, but PID controllers do not have the ability to learn and must be set up correctly. Selecting the correct gains for effective control is known as tuning the controller.
If a controller starts from a stable state at zero error (PV = SP), then further changes by the controller will be in response to changes in other measured or unmeasured inputs to the process that impact on the process, and hence on the PV. Variables that impact on the process other than the MV are known as disturbances. Generally controllers are used to reject disturbances and/or implement setpoint changes. Changes in feed water temperature constitute a disturbance to the shower process.
In theory, a controller can be used to control any process which has a measurable output (PV), a known ideal value for that output (SP) and an input to the process (MV) that will affect the relevant PV. Controllers are used in industry to regulate temperature, pressure, flow rate, chemical composition, speed and practically every other variable for which a measurement exists. Automobile cruise control is an example of a process which utilizes automated control.
Due to their long history, simplicity, well grounded theory and simple setup and maintenance requirements, PID controllers are the controllers of choice for many of these applications.
PID controller theory
This section describes the parallel or non-interacting form of the PID controller. For other forms please see the Section "Alternative notation and PID forms".
The PID control scheme is named after its three correcting terms, whose sum constitutes the manipulated variable (MV). Hence:
where
- P_{out}, I_{out}, and D_{out} are the contributions to the output from the PID controller from each of the three terms, as defined below.
[edit] Proportional term
The proportional term (sometimes called gain) makes a change to the output that is proportional to the current error value. The proportional response can be adjusted by multiplying the error by a constant K_{p}, called the proportional gain.
The proportional term is given by:
where
- P_{out}: Proportional term of output
- K_{p}: Proportional gain, a tuning parameter
- e: Error = SP − PV
- t: Time or instantaneous time (the present)
A high proportional gain results in a large change in the output for a given change in the error. If the proportional gain is too high, the system can become unstable (See the section on loop tuning). In contrast, a small gain results in a small output response to a large input error, and a less responsive (or sensitive) controller. If the proportional gain is too low, the control action may be too small when responding to system disturbances.
In the absence of disturbances, pure proportional control will not settle at its target value, but will retain a steady state error that is a function of the proportional gain and the process gain. Despite the steady-state offset, both tuning theory and industrial practice indicate that it is the proportional term that should contribute the bulk of the output change.
[edit] Integral term
The contribution from the integral term (sometimes called reset) is proportional to both the magnitude of the error and the duration of the error. Summing the instantaneous error over time (integrating the error) gives the accumulated offset that should have been corrected previously. The accumulated error is then multiplied by the integral gain and added to the controller output. The magnitude of the contribution of the integral term to the overall control action is determined by the integral gain, K_{i}.
The integral term is given by:
where
- I_{out}: Integral term of output
- K_{i}: Integral gain, a tuning parameter
- e: Error = SP − PV
- t: Time or instantaneous time (the present)
- τ: a dummy integration variable
The integral term (when added to the proportional term) accelerates the movement of the process towards setpoint and eliminates the residual steady-state error that occurs with a proportional only controller. However, since the integral term is responding to accumulated errors from the past, it can cause the present value to overshoot the setpoint value (cross over the setpoint and then create a deviation in the other direction). For further notes regarding integral gain tuning and controller stability, see the section on loop tuning.
Derivative term
The rate of change of the process error is calculated by determining the slope of the error over time (i.e., its first derivative with respect to time) and multiplying this rate of change by the derivative gain K_{d}. The magnitude of the contribution of the derivative term (sometimes called rate) to the overall control action is termed the derivative gain, K_{d}.
The derivative term is given by:
where
- D_{out}: Derivative term of output
- K_{d}: Derivative gain, a tuning parameter
- e: Error = SP − PV
- t: Time or instantaneous time (the present)
The derivative term slows the rate of change of the controller output and this effect is most noticeable close to the controller setpoint. Hence, derivative control is used to reduce the magnitude of the overshoot produced by the integral component and improve the combined controller-process stability. However, differentiation of a signal amplifies noise and thus this term in the controller is highly sensitive to noise in the error term, and can cause a process to become unstable if the noise and the derivative gain are sufficiently large.
[edit] Summary
The proportional, integral, and derivative terms are summed to calculate the output of the PID controller. Defining u(t) as the controller output, the final form of the PID algorithm is:
where the tuning parameters are:
- Proportional gain, K_{p}
- larger values typically mean faster response since the larger the error, the larger the Proportional term compensation. An excessively large proportional gain will lead to process instability and oscillation.
- Integral gain, K_{i}
- larger values imply steady state errors are eliminated more quickly. The trade-off is larger overshoot: any negative error integrated during transient response must be integrated away by positive error before we reach steady state.
- Derivative gain, K_{d}
- larger values decrease overshoot, but slows down transient response and may lead to instability due to signal noise amplification in the differentiation of the error.
[edit] Loop tuning
If the PID controller parameters (the gains of the proportional, integral and derivative terms) are chosen incorrectly, the controlled process input can be unstable, i.e. its output diverges, with or without oscillation, and is limited only by saturation or mechanical breakage. Tuning a control loop is the adjustment of its control parameters (gain/proportional band, integral gain/reset, derivative gain/rate) to the optimum values for the desired control response.
The optimum behavior on a process change or setpoint change varies depending on the application. Some processes must not allow an overshoot of the process variable beyond the setpoint if, for example, this would be unsafe. Other processes must minimize the energy expended in reaching a new setpoint. Generally, stability of response (the reverse of instability) is required and the process must not oscillate for any combination of process conditions and setpoints. Some processes have a degree of non-linearity and so parameters that work well at full-load conditions don't work when the process is starting up from no-load. This section describes some traditional manual methods for loop tuning.
There are several methods for tuning a PID loop. The most effective methods generally involve the development of some form of process model, then choosing P, I, and D based on the dynamic model parameters. Manual tuning methods can be relatively inefficient.
The choice of method will depend largely on whether or not the loop can be taken "offline" for tuning, and the response time of the system. If the system can be taken offline, the best tuning method often involves subjecting the system to a step change in input, measuring the output as a function of time, and using this response to determine the control parameters.
Choosing a Tuning Method | ||
Method | Advantages | Disadvantages |
Manual Tuning | No math required. Online method. | Requires experienced personnel. |
Ziegler–Nichols | Proven Method. Online method. | Process upset, some trial-and-error, very aggressive tuning. |
Software Tools | Consistent tuning. Online or offline method. May include valve and sensor analysis. Allow simulation before downloading. | Some cost and training involved. |
Cohen-Coon | Good process models. | Some math. Offline method. Only good for first-order processes. |
[edit] Manual tuning
If the system must remain online, one tuning method is to first set K_{i} and K_{d} values to zero. Increase the K_{p} until the output of the loop oscillates, then the K_{p} should be left set to be approximately half of that value for a "quarter amplitude decay" type response. Then increase K_{i} until any offset is correct in sufficient time for the process. However, too much K_{i} will cause instability. Finally, increase K_{d}, if required, until the loop is acceptably quick to reach its reference after a load disturbance. However, too much K_{d} will cause excessive response and overshoot. A fast PID loop tuning usually overshoots slightly to reach the setpoint more quickly; however, some systems cannot accept overshoot, in which case an "over-damped" closed-loop system is required, which will require a K_{p} setting significantly less than half that of the K_{p} setting causing oscillation.
Effects of increasing parameters | ||||
Parameter | Rise time | Overshoot | Settling time | Error at equilibrium |
K_{p} | Decrease | Increase | Small change | Decrease |
K_{i} | Decrease | Increase | Increase | Eliminate |
K_{d} | Indefinite (small decrease or increase)^{[1]} | Decrease | Decrease | None |
[edit] Ziegler–Nichols method
Another tuning method is formally known as the Ziegler–Nichols method, introduced by John G. Ziegler and Nathaniel B. Nichols. As in the method above, the K_{i} and K_{d} gains are first set to zero. The P gain is increased until it reaches the critical gain, K_{c}, at which the output of the loop starts to oscillate. K_{c} and the oscillation period P_{c} are used to set the gains as shown:
Ziegler–Nichols method | ||||
Control Type | K_{p} | K_{i} | K_{d} | |
P | 0.50K_{c} | - | - | |
PI | 0.45K_{c} | 1.2K_{p} / P_{c} | - | |
PID | 0.60K_{c} | 2K_{p} / P_{c} | K_{p}P_{c} / 8 |
[edit] PID tuning software
Most modern industrial facilities no longer tune loops using the manual calculation methods shown above. Instead, PID tuning and loop optimization software are used to ensure consistent results. These software packages will gather the data, develop process models, and suggest optimal tuning. Some software packages can even develop tuning by gathering data from reference changes.
Mathematical PID loop tuning induces an impulse in the system, and then uses the controlled system's frequency response to design the PID loop values. In loops with response times of several minutes, mathematical loop tuning is recommended, because trial and error can literally take days just to find a stable set of loop values. Optimal values are harder to find. Some digital loop controllers offer a self-tuning feature in which very small setpoint changes are sent to the process, allowing the controller itself to calculate optimal tuning values.
Other formulas are available to tune the loop according to different performance criteria.
[edit] Modifications to the PID algorithm
The basic PID algorithm presents some challenges in control applications that have been addressed by minor modifications to the PID form.
One common problem resulting from the ideal PID implementations is integral windup. This problem can be addressed by:
- Initializing the controller integral to a desired value
- Increasing the setpoint in a suitable ramp
- Disabling the integral function until the PV has entered the controllable region
- Limiting the time period over which the integral error is calculated
- Preventing the integral term from accumulating above or below pre-determined bounds
- Freezing the integral function in case of disturbances
- If a PID loop is used to control the temperature of an electric resistance furnace, the system has stabilized and then the door is opened and something cold is put into the furnace the temperature drops below the setpoint. The integral function of the controller tends to compensate this error by introducing another error in the positive direction. This can be avoided by freezing of the integral function after the opening of the door for the time the control loop typically needs to reheat the furnace.
- Replacing the integral function by a model based part
- Often the time-response of the system is approximately known. Then it is an advantage to simulate this time-response with a model and to calculate some unknown parameter from the actual response of the system. If for instance the system is an electrical furnace the response of the difference between furnace temperature and ambient temperature to changes of the electrical power will be similar to that of a simple RC low-pass filter multiplied by an unknown proportional coefficient. The actual electrical power supplied to the furnace is delayed by a low-pass filter to simulate the response of the temperature of the furnace and then the actual temperature minus the ambient temperature is divided by this low-pass filtered electrical power. Then, the result is stabilized by another low-pass filter leading to an estimation of the proportional coefficient. With this estimation it is possible to calculate the required electrical power by dividing the set-point of the temperature minus the ambient temperature by this coefficient. The result can then be used instead of the integral function. This also achieves a control error of zero in the steady-state but avoids integral windup and can give a significantly improved control action compared to an optimized PID controller. This type of controller does work properly in an open loop situation which causes integral windup with an integral function. This is an advantage if for example the heating of a furnace has to be reduced for some time because of the failure of a heating element or if the controller is used as an advisory system to a human operator who may or may not switch it to closed-loop operation or if the controller is used inside of a branch of a complex control system where this branch may be temporarily inactive.
Many PID loops control a mechanical device (for example, a valve). Mechanical maintenance can be a major cost and wear leads to control degradation in the form of either stiction or a deadband in the mechanical response to an input signal. The rate of mechanical wear is mainly a function of how often a device is activated to make a change. Where wear is a significant concern, the PID loop may have an output deadband to reduce the frequency of activation of the output (valve). This is accomplished by modifying the controller to hold its output steady if the change would be small (within the defined deadband range). The calculated output must leave the deadband before the actual output will change.
The proportional and derivative terms can produce excessive movement in the output when a system is subjected to an instantaneous step increase in the error, such as a large setpoint change. In the case of the derivative term, this is due to taking the derivative of the error, which is very large in the case of an instantaneous step change. As a result, some PID algorithms incorporate the following modifications:
- Derivative of output
- In this case the PID controller measures the derivative of the output quantity, rather than the derivative of the error. The output is always continuous (i.e., never has a step change). For this to be effective, the derivative of the output must have the same sign as the derivative of the error.
- Setpoint ramping
- In this modification, the setpoint is gradually moved from its old value to a newly specified value using a linear or first order differential ramp function. This avoids the discontinuity present in a simple step change.
- Setpoint weighting
- Setpoint weighting uses different multipliers for the error depending on which element of the controller it is used in. The error in the integral term must be the true control error to avoid steady-state control errors. This affects the controller's setpoint response. These parameters do not affect the response to load disturbances and measurement noise.
[edit] Limitations of PID control
While PID controllers are applicable to many control problems, they can perform poorly in some applications.
PID controllers, when used alone, can give poor performance when the PID loop gains must be reduced so that the control system does not overshoot, oscillate or hunt about the control setpoint value. The control system performance can be improved by combining the feedback (or closed-loop) control of a PID controller with feed-forward (or open-loop) control. Knowledge about the system (such as the desired acceleration and inertia) can be fed forward and combined with the PID output to improve the overall system performance. The feed-forward value alone can often provide the major portion of the controller output. The PID controller can then be used primarily to respond to whatever difference or error remains between the setpoint (SP) and the actual value of the process variable (PV). Since the feed-forward output is not affected by the process feedback, it can never cause the control system to oscillate, thus improving the system response and stability.
For example, in most motion control systems, in order to accelerate a mechanical load under control, more force or torque is required from the prime mover, motor, or actuator. If a velocity loop PID controller is being used to control the speed of the load and command the force or torque being applied by the prime mover, then it is beneficial to take the instantaneous acceleration desired for the load, scale that value appropriately and add it to the output of the PID velocity loop controller. This means that whenever the load is being accelerated or decelerated, a proportional amount of force is commanded from the prime mover regardless of the feedback value. The PID loop in this situation uses the feedback information to effect any increase or decrease of the combined output in order to reduce the remaining difference between the process setpoint and the feedback value. Working together, the combined open-loop feed-forward controller and closed-loop PID controller can provide a more responsive, stable and reliable control system.
Another problem faced with PID controllers is that they are linear. Thus, performance of PID controllers in non-linear systems (such as HVAC systems) is variable. Often PID controllers are enhanced through methods such as PID gain scheduling or fuzzy logic. Further practical application issues can arise from instrumentation connected to the controller. A high enough sampling rate, measurement precision, and measurement accuracy are required to achieve adequate control performance.
A problem with the Derivative term is that small amounts of measurement or process noise can cause large amounts of change in the output. It is often helpful to filter the measurements with a low-pass filter in order to remove higher-frequency noise components. However, low-pass filtering and derivative control can cancel each other out, so reducing noise by instrumentation means is a much better choice. Alternatively, the differential band can be turned off in many systems with little loss of control. This is equivalent to using the PID controller as a PI controller.
[edit] Cascade control
One distinctive advantage of PID controllers is that two PID controllers can be used together to yield better dynamic performance. This is called cascaded PID control. In cascade control there are two PIDs arranged with one PID controlling the set point of another. A PID controller acts as outer loop controller, which controls the primary physical parameter, such as fluid level or velocity. The other controller acts as inner loop controller, which reads the output of outer loop controller as set point, usually controlling a more rapid changing parameter, flowrate or acceleration. It can be mathematically proven^{[citation needed]} that the working frequency of the controller is increased and the time constant of the object is reduced by using cascaded PID controller.^{[vague]}.
[edit] Physical implementation of PID control
In the early history of automatic process control the PID controller was implemented as a mechanical device. These mechanical controllers used a lever, spring and a mass and were often energized by compressed air. These pneumatic controllers were once the industry standard.
Electronic analog controllers can be made from a solid-state or tube amplifier, a capacitor and a resistance. Electronic analog PID control loops were often found within more complex electronic systems, for example, the head positioning of a disk drive, the power conditioning of a power supply, or even the movement-detection circuit of a modern seismometer. Nowadays, electronic controllers have largely been replaced by digital controllers implemented with microcontrollers or FPGAs.
Most modern PID controllers in industry are implemented in programmable logic controllers (PLCs) or as a panel-mounted digital controller. Software implementations have the advantages that they are relatively cheap and are flexible with respect to the implementation of the PID algorithm.
[edit] Alternative nomenclature and PID forms
[edit] Ideal versus standard PID form
The form of the PID controller most often encountered in industry, and the one most relevant to tuning algorithms is the standard form. In this form the K_{p} gain is applied to the I_{out}, and D_{out} terms, yielding:
where
- T_{i} is the integral time
- T_{d} is the derivative time
In the ideal parallel form, shown in the controller theory section
the gain parameters are related to the parameters of the standard form through and . This parallel form, where the parameters are treated as simple gains, is the most general and flexible form. However, it is also the form where the parameters have the least physical interpretation and is generally reserved for theoretical treatment of the PID controller. The standard form, despite being slightly more complex mathematically, is more common in industry.
[edit] Laplace form of the PID controller
Sometimes it is useful to write the PID regulator in Laplace transform form:
Having the PID controller written in Laplace form and having the transfer function of the controlled system, makes it easy to determine the closed-loop transfer function of the system.
[edit] Series/interacting form
Another representation of the PID controller is the series, or interacting form
where the parameters are related to the parameters of the standard form through
- , , and
with
- .
This form essentially consists of a PD and PI controller in series, and it made early (analog) controllers easier to build. When the controllers later became digital, many kept using the interacting form.
[edit] Discrete implementation
The analysis for designing a digital implementation of a PID controller in a Microcontroller (MCU) or FPGA device requires the standard form of the PID controller to be discretised ^{[2]}. Approximations for first-order derivatives are made by backward finite differences. The integral term is discretised, with a sampling time Δt,as follows,
The derivative term is approximated as,
Thus, a velocity algorithm for implementation of the discretised PID controller in a MCU is obtained,
[edit] Pseudocode
Here is a simple software loop that implements the PID algorithm:
previous_error = 0 integral = 0 start: error = setpoint - actual_position integral = integral + (error*dt) derivative = (error - previous_error)/dt output = (Kp*error) + (Ki*integral) + (Kd*derivative) previous_error = error wait(dt) goto start
[edit] External links
[edit] PID tutorials
- PID Tutorial
- P.I.D. Without a PhD: a beginner's guide to PID loop theory with sample programming code
- What's All This P-I-D Stuff, Anyhow? Article in Electronic Design
- Shows how to build a PID controller with basic electronic components (pg. 22)
[edit] Simulations
- Free, real-time PID simulator for Windows
- PID controller using MatLab and Simulink
- PID controller laboratory, Java applets for PID tuning
[edit] Special topics and PID control applications
- Proven Methods and Best Practices for PID Control
- PID Control Primer Article in Embedded Systems Programming
[edit] References
- ^ http://saba.kntu.ac.ir/eecd/pcl/download/PIDtutorial.pdf
- ^ http://www.scribd.com/doc/19070283/Discrete-PI-and-PID-Controller-Design-and-Analysis-for-Digital-Implementation
- Liptak, Bela (1995). Instrument Engineers' Handbook: Process Control. Radnor, Pennsylvania: Chilton Book Company. pp. 20–29. ISBN 0-8019-8242-1.
- Tan, Kok Kiong; Wang Qing-Guo, Hang Chang Chieh (1999). Advances in PID Control. London, UK: Springer-Verlag. ISBN 1-85233-138-0.
- Van, Doren, Vance J. (July 1, 2003). "Loop Tuning Fundamentals". Control Engineering (Red Business Information). http://www.controleng.com/article/CA307745.html.
- Sellers, David. "An Overview of Proportional plus Integral plus Derivative Control and Suggestions for Its Successful Application and Implementation" (PDF). http://www.peci.org/library/PECI_ControlOverview1_1002.pdf. Retrieved 2007-05-05.
- Graham, Ron (10/03/2005). "FAQ on PID controller tuning". http://web.archive.org/web/20050206113949/www.tcnj.edu/~rgraham/PID-tuning.html. Retrieved 2009-01-05.
Interesting post and thanks for sharing. Some things in here I have not thought about before. Thanks for making such a cool post.
that's will be good for my reference on my website.
thanks for the post, yes you're right about PID
good... job
Wow.., kontennya bagus.., buat perkembangan pendidikan :)
Kunjungi juga Hon Book Store ya.., silahkan pilih2 bukunya, mau liat2 aja juga boleh :)
Kalo' beli dapatkan diskon hingga 15%.
http://www.honbookstore.com/
thank you very much for this information so valuable webcam sexy
I read this article and I loved it, congratulations
lineas eroticas
It's a great pleasure reading your post.It's full of information I am looking for and I love to post a comment that "The content of your post is awesome" Great work!
www.google.com
Hi there There. I discovered your weblog using msn. This is a rather well written article. I'll be certain toCara Menaikan berat bada
Pelangsing badan
Menu diet menurunkan berat badan
Makanan untuk diet
bookmark it and return to learn more of your helpful information. Thanks for the post. I will positively return
This website is really a walk-through for all of the details you desired about this and did not know who to ask.Certainly a amazing item of content ... It has appropriate details. Thanks for publishing this.
http://social24h.com/
The clarity in your post is easily spectacular and I can assume you are a professional on this field. Well with your consent allow me to catch your increases feed to hold up to designated day with incoming post.
http://www.hknightlifeescort.com/.
This website is excellent and so is how the subject matter was explained. I also like a number of the comments too.Ready for
subsequent post.
resep makanan diet
makanan sehat untuk anak
susu untuk diet
cara sehat menurunkan berat badan
resep smoothies untuk diet
beras merah untuk diet
manfaat buah pepaya untuk diet
teh hijau untuk diet
cara untuk memutikan wajah
thanks for sharing...
obat pelangsing alami
kurang mengerti memang
obat kurus
PIDthe best thing. Nice post