Displaying publications 1 - 20 of 140 in total

Abstract:
Sort:
  1. Syahrom A, Abdul Kadir MR, Abdullah J, Öchsner A
    Med Biol Eng Comput, 2011 Dec;49(12):1393-403.
    PMID: 21947767 DOI: 10.1007/s11517-011-0833-0
    The relationship between microarchitecture to the failure mechanism and mechanical properties can be assessed through experimental and computational methods. In this study, both methods were utilised using bovine cadavers. Twenty four samples of cancellous bone were extracted from fresh bovine and the samples were cleaned from excessive marrow. Uniaxial compression testing was performed with displacement control. After mechanical testing, each specimen was ashed in a furnace. Four of the samples were exemplarily scanned using micro-computed tomography (μCT) and three dimensional models of the cancellous bones were reconstructed for finite element simulation. The mechanical properties and the failure modes obtained from numerical simulations were then compared to the experiments. Correlations between microarchitectural parameters to the mechanical properties and failure modes were then made. The Young's modulus correlates well with the bone volume fraction with R² = 0.615 and P value 0.013. Three different types of failure modes of cancellous bone were observed: oblique fracture (21.7%), perpendicular global fracture (47.8%), and scattered localised fracture (30.4%). However, no correlations were found between the failure modes to the morphological parameters. The percentage of error between computer predictions and the actual experimental test was from 6 to 12%. These mechanical properties and information on failure modes can be used for the development of synthetic cancellous bone.
    Matched MeSH terms: Finite Element Analysis
  2. Norliza Mohd. Zain, Zuhaila Ismail
    MATEMATIKA, 2019;35(2):213-227.
    MyJurnal
    Blood flow through a bifurcated artery with the presence of an overlapping stenosis located at parent’s arterial lumen under the action of a uniform external magnetic field is studied in this paper. Blood is treated as an electrically conducting fluid which exhibits the Magnetohydrodynamics principle and it is characterized by a Newtonian fluid model. The governing equations are discretized using a stabilization technique of finite element known as Galerkin least-squares. The maximum velocity and pressure drop evaluated in this present study are compared with the results found in previous literature and COMSOL Multiphysics. The solutions found in a satisfactory agreement, thus verify the source code is working properly. The effects of dimensionless parameters of Hartmann and Reynolds numbers in the fluid’s velocity and pressure are examined in details with further scientific discussions.
    Matched MeSH terms: Finite Element Analysis
  3. Ismail, A.Y., Noerpamoengkas, A., Zakaria, S.I.F.S.
    MyJurnal
    In this paper, an extensive work on the natural frequency of perforated plate has
    been made by introducing micro-holes on the plate. The micro-holes, which is known
    to have a remarkable performance in several applications, were arranged diagonally
    among the perforated holes in order to make a new combination of micro-macro
    perforated plate. A 3D geometrical model of the plate was made in the Autodesk
    Inventor and the Finite Element (FE) simulation was employed to calculate the
    natural frequency and visualize the mode shape. Four models were made with
    various micro-holes diameter starting from 0,25, 0,5, 0,75 and 0,9 mm, respectively.
    The macro holes diameter, however, was kept constant to be 1,5 mm purposely to
    know the effect of micro holes in particular. The results from the models were then
    compared to a single perforated plate (PP) for clarification. It is found that the micro
    holes gives considerable effect to the perforated plate natural frequency. Similar to
    the author previous findings, the diameter of the micro holes is proportional to the
    natural frequency reduction.
    Matched MeSH terms: Finite Element Analysis
  4. Mutafi A, Yidris N, Ishak MR, Zahari R
    Heliyon, 2018 Nov;4(11):e00937.
    PMID: 30839801 DOI: 10.1016/j.heliyon.2018.e00937
    Steel sections are normally shaped via cold work manufacturing processes. The extent of cold work to shape the steel sections might induce residual stresses in the region of bending. Previously, researchers had performed studies on the influences of local buckling on the failure behavior of steel compression members which shown that failure will happen when most of the yielding has extended to the middle surface in the bend region of the sections. Therefore, these cold work methods may have major effect on the behavior of the steel section and also its load-bearing capability. In addition, another factor may play significant role in formed section's load-bearing capacity which is the longitudinal residual strain. The longitudinal residual strain raised during forming procedure can be used to define the section imperfection of the formed section and its relation to the existence of defects. Therefore, the main motivation of this research paper is to perform three-dimensional finite element (3D-FE) to investigate peak longitudinal residual strains of a thin-walled steel plate with large bending angle along member length. A 3D finite element simulation in ABAQUS has been employed to simulate this forming process. The study concluded that the longitudinal residual strain at the section corner edge was higher than those at the rest of the corner region. These strains at the edge were higher than the yield strain


    (



    ε


    y



    )


    of the formed section which occurred due to the lack of transverse restraint. This made the plate edge tended to bend toward the normal direction when it was under a high transverse bending. This causes a significant difference in longitudinal strain at the plate edge.
    Matched MeSH terms: Finite Element Analysis
  5. Chuah HG, Abd Rahim I, Yusof MI
    Comput Methods Biomech Biomed Engin, 2010 Jun;13(3):319-26.
    PMID: 19927241 DOI: 10.1080/10255840903208189
    The stress shielding effect is an event in which the replacement implant limits the load transferred to bone and the ineffective stress in the vertebrae causes bony growth to cease. In the present study, a 3D finite element L4-L5 model was developed and subjected to a 1200 N compression preload. Five groups of muscle forces were applied on L4 under flexion-extension, lateral bending and axial rotation. Topology optimisation was employed for reducing the stress shielding effect by removing the ineffective material from the design domain. The optimised design was designed with polyaryletheretherketone (PEEK) titanium and cortical materials to encounter the shielding response. The stress responses show that the new design increased the stress magnitude by at least 17.10, 18.11 and 18.43% in 4 Nm of flexion-extension, lateral bending and axial rotation, respectively. In conclusion, the material factor did not significantly alter the stress magnitude, but volume was the key factor in reducing the stress shielding effect.
    Matched MeSH terms: Finite Element Analysis
  6. Mohd Razip Wee MF, Jaafar MM, Faiz MS, Dee CF, Yeop Majlis B
    Biosensors (Basel), 2018 Dec 05;8(4).
    PMID: 30563159 DOI: 10.3390/bios8040124
    Gallium Nitride (GaN) is considered as the second most popular semiconductor material in industry after silicon. This is due to its wide applications encompassing Light Emitting Diode (LED) and power electronics. In addition, its piezoelectric properties are fascinating to be explored as electromechanical material for the development of diverse microelectromechanical systems (MEMS) application. In this article, we conducted a theoretical study concerning surface mode propagation, especially Rayleigh and Sezawa mode in the layered GaN/sapphire structure with the presence of various guiding layers. It is demonstrated that the increase in thickness of guiding layer will decrease the phase velocities of surface mode depending on the material properties of the layer. In addition, the Q-factor value indicating the resonance properties of surface mode appeared to be affected with the presence of fluid domain, particularly in the Rayleigh mode. Meanwhile, the peak for Sezawa mode shows the highest Q factor and is not altered by the presence of fluid. Based on these theoretical results using the finite element method, it could contribute to the development of a GaN-based device to generate surface acoustic wave, especially in Sezawa mode which could be useful in acoustophoresis, lab on-chip and microfluidics applications.
    Matched MeSH terms: Finite Element Analysis
  7. An X, Chong PL, Zohourkari I, Roy S, Merdji A, Linda Gnanasagaran C, et al.
    Proc Inst Mech Eng H, 2023 Aug;237(8):1008-1016.
    PMID: 37477395 DOI: 10.1177/09544119231187685
    The mechanical properties of tissue scaffolds are essential in providing stability for tissue repair and growth. Thus, the ability of scaffolds to withstand specific loads is crucial for scaffold design. Most research on scaffold pores focuses on grids with pore size and gradient structure, and many research models are based on scaffolding with vertically arranged holes. However, little attention is paid to the influence of the distribution of holes on the mechanical properties of the scaffold. To address this gap, this research investigates the effect of pore distribution on the mechanical properties of tissue scaffolds. The study involves four types of scaffold designs with regular and staggered pore arrangements and porosity ranging from 30% to 80%. Finite element analysis (FEA) was used to compare the mechanical properties of different scaffold designs, with von-Mises stress distribution maps generated for each scaffold. The results show that scaffolds with regular vertical holes exhibit a more uniform stress distribution and better mechanical performance than those with irregular holes. In contrast, the scaffold with a staggered arrangement of holes had a higher probability of stress concentration. The study emphasized the importance of balancing porosity and strength in scaffold design.
    Matched MeSH terms: Finite Element Analysis
  8. Harith H, Schmutz B, Malekani J, Schuetz MA, Yarlagadda PK
    Med Eng Phys, 2016 Mar;38(3):280-5.
    PMID: 26739124 DOI: 10.1016/j.medengphy.2015.11.012
    Anatomically precontoured plates are commonly used to treat periarticular fractures. A well-fitting plate can be used as a tool for anatomical reduction of the fractured bone. Recent studies highlighted that some plates fit poorly for many patients due to considerable shape variations between bones of the same anatomical site. While it is impossible to design one shape that fits all, it is also burdensome for the manufacturers and hospitals to produce, store and manage multiple plate shapes without the certainty of utilization by a patient population. In this study, we investigated the number of shapes required for maximum fit within a given dataset, and if they could be obtained by manually deforming the original plate. A distal medial tibial plate was automatically positioned on 45 individual tibiae, and the optimal deformation was determined iteratively using finite element analysis simulation. Within the studied dataset, we found that: (i) 89% fit could be achieved with four shapes, (ii) 100% fit was impossible through mechanical deformation, and (iii) the deformations required to obtain the four plate shapes were safe for the stainless steel plate for further clinical use. The proposed framework is easily transferable to other orthopaedic plates.
    Matched MeSH terms: Finite Element Analysis
  9. Karimzadeh A, R Koloor SS, Ayatollahi MR, Bushroa AR, Yahya MY
    Sci Rep, 2019 10 31;9(1):15763.
    PMID: 31673118 DOI: 10.1038/s41598-019-51904-4
    This study investigates the capacity of the nano-indentation method in the mechanical characterization of a heterogeneous dental restorative nanocomposite using experimental and computational approaches. In this respect, Filtek Z350 XT was selected as a nano-particle reinforced polymer nanocomposite with a specific range of the particle size (50 nm to 4 µm), within the range of indenter contact area of the nano-indentation experiment. A Sufficient number of nano-indentation tests were performed in various locations of the nanocomposite to extract the hardness and elastic modulus properties. A hybrid computational-experimental approach was developed to examine the extracted properties by linking the internal behaviour and the global response of the nanocomposite. In the computational part, several representative models of the nanocomposite were created in a finite element environment to simulate the mechanism of elastic-plastic deformation of the nanocomposite under Berkovich indenter. Dispersed values of hardness and elastic modulus were obtained through the experiment with 26.8 and 48.5 percent average errors, respectively, in comparison to the nanocomposite properties, respectively. A disordered shape was predicted for plastic deformation of the equilateral indentation mark, representing the interaction of the particles and matrix, which caused the experiment results reflect the local behaviour of the nanocomposite instead of the real material properties.
    Matched MeSH terms: Finite Element Analysis
  10. Seow LL, Toh CG, Fok AS, Wilson NH
    Am J Dent, 2008 Oct;21(5):331-6.
    PMID: 19024261
    PURPOSE: To investigate the level and distribution of stresses in endodontically treated maxillary premolar teeth restored using various cavity designs of bonded all-ceramic restorations. The hypothesis tested was that the various all-ceramic approaches, including incorporating a pulp chamber extension in the restoration, had no influence on the stresses in the restored tooth unit.
    METHODS: Finite element packages Patran and Abaqus were used for the stress analysis. The cavity designs investigated include: (1) inlay (I); (2) inlay with palatal cusp coverage (IPC); (3) onlay (O); (4) inlay with pulp chamber extension (IPE); (5) inlay with palatal cusp coverage and pulp chamber extension (IPCPE); and (6) onlay with pulp chamber extension (OPE).
    RESULTS: In each case, tensile stresses were found to be concentrated subjacent to the occlusal fossa. Peak tensile stress and peak shear stress values along the tooth/restoration interface for IPC, O IPCPE and OPE cavity designs were found to be associated with the axiogingival line angle. Overall, the order of the various forms of restoration investigated in terms of the maximum principal stress (from greatest to lowest) was as follows: IPE > IPCPE > OPE > I > IPC > O.
    Matched MeSH terms: Finite Element Analysis
  11. Pirouzi G, Abu Osman NA, Eshraghi A, Ali S, Gholizadeh H, Wan Abas WA
    ScientificWorldJournal, 2014;2014:849073.
    PMID: 25197716 DOI: 10.1155/2014/849073
    Socket is an important part of every prosthetic limb as an interface between the residual limb and prosthetic components. Biomechanics of socket-residual limb interface, especially the pressure and force distribution, have effect on patient satisfaction and function. This paper aimed to review and evaluate studies conducted in the last decades on the design of socket, in-socket interface pressure measurement, and socket biomechanics. Literature was searched to find related keywords with transtibial amputation, socket-residual limb interface, socket measurement, socket design, modeling, computational modeling, and suspension system. In accordance with the selection criteria, 19 articles were selected for further analysis. It was revealed that pressure and stress have been studied in the last decaeds, but quantitative evaluations remain inapplicable in clinical settings. This study also illustrates prevailing systems, which may facilitate improvements in socket design for improved quality of life for individuals ambulating with transtibial prosthesis. It is hoped that the review will better facilitate the understanding and determine the clinical relevance of quantitative evaluations.
    Matched MeSH terms: Finite Element Analysis
  12. Panjehpour M, Chai HK, Voo YL
    PLoS One, 2015;10(6):e0130734.
    PMID: 26110268 DOI: 10.1371/journal.pone.0130734
    Deep beams are commonly used in tall buildings, offshore structures, and foundations. According to many codes and standards, strut-and-tie model (STM) is recommended as a rational approach for deep beam analyses. This research focuses on the STM recommended by ACI 318-11 and AASHTO LRFD and uses experimental results to modify the strut effectiveness factor in STM for reinforced concrete (RC) deep beams. This study aims to refine STM through the strut effectiveness factor and increase result accuracy. Six RC deep beams with different shear span to effective-depth ratios (a/d) of 0.75, 1.00, 1.25, 1.50, 1.75, and 2.00 were experimentally tested under a four-point bending set-up. The ultimate shear strength of deep beams obtained from non-linear finite element modeling and STM recommended by ACI 318-11 as well as AASHTO LRFD (2012) were compared with the experimental results. An empirical equation was proposed to modify the principal tensile strain value in the bottle-shaped strut of deep beams. The equation of the strut effectiveness factor from AASHTTO LRFD was then modified through the aforementioned empirical equation. An investigation on the failure mode and crack propagation in RC deep beams subjected to load was also conducted.
    Matched MeSH terms: Finite Element Analysis
  13. Mat Daud NI, Viswanathan KK
    PLoS One, 2019;14(7):e0219089.
    PMID: 31269073 DOI: 10.1371/journal.pone.0219089
    Vibrational behaviour of symmetric angle-ply layered circular cylindrical shell filled with quiescent fluid is presented. The equations of motion of cylindrical shell in terms of stress and moment resultants are derived from the first order shear deformation theory. Irrotational of inviscid fluid are expressed as the wave equation. These two equations are coupled. Strain-displacement relations and stress-strain relations are adopted into the equations of motion to obtain the differential equations with displacements and rotational functions. A system of ordinary differential equation is obtained in one variable by assuming the functions in separable form. Spline of order three is applied to approximate the displacement and rotational functions, together with boundary conditions, to get a generalised eigenvalue problem. The eigenvalue problem is solved for eigen frequency parameter and associate eigenvectors of spline coefficients. The study of frequency parameters are analysed using the parameters the thickness ratio, length ratio, angle-ply, properties of material and number of layers under different boundary conditions.
    Matched MeSH terms: Finite Element Analysis
  14. Patil PG, Seow LL, Uddanwadikar R, Ukey PD
    J Prosthet Dent, 2021 Jan;125(1):138.e1-138.e8.
    PMID: 33393474 DOI: 10.1016/j.prosdent.2020.09.015
    STATEMENT OF PROBLEM: Mini implants (<3 mm in diameter) are being used as an alternative to standard implants for implant-retained mandibular overdentures; however, they may exhibit higher stresses at the crestal level.

    PURPOSE: The purpose of this finite element analysis study was to evaluate the biomechanical behavior (stress distribution pattern) in the mandibular overdenture, mucosa, bone, and implants when retained with 2 standard implants or 2 mini implants under unilateral or bilateral loading conditions.

    MATERIAL AND METHODS: A patient with edentulous mandible and his denture was scanned with cone beam computed tomography (CBCT), and a 3D mandibular model was created in the Mimics software program by using the CBCT digital imaging and communications in medicine (DICOM) images. The model was transferred to the 3Matics software program to form a 2-mm-thick mucosal layer and to assemble the denture DICOM file. A 12-mm-long standard implant (Ø3.5 mm) and a mini dental implant (Ø2.5 mm) along with the LOCATOR male attachments (height 4 mm) were designed by using the SOLIDWORKS software program. Two standard or 2 mini implants in the canine region were embedded separately in the 3D assembled model. The base of the mandible was fixed, and vertical compressive loads of 100 N were applied unilaterally and bilaterally in the first molar region. The material properties for acrylic resin (denture), titanium (implants), mucosa (tissue), and bone (mandible) were allocated. Maximum von Mises stress and strain values were obtained and analyzed.

    RESULTS: Maximum stresses of 9.78 MPa (bilaterally) and 11.98 MPa (unilaterally) were observed in 2 mini implants as compared with 3.12 MPa (bilaterally) and 3.81 MPa (unilaterally) in 2 standard implants. The stress values in the mandible were observed to be almost double the mini implants as compared with the standard implants. The stresses in the denture were in the range of 3.21 MPa and 3.83 MPa and in the mucosa of 0.68 MPa and 0.7 MPa for 2 implants under unilateral and bilateral loading conditions. The strain values shown similar trends with both implant types under bilateral and unilateral loading.

    CONCLUSIONS: Two mini implants generated an average of 68.15% more stress than standard implants. The 2 standard implant-retained overdenture showed less stress concentration in and around implants than mini implant-retained overdentures.

    Matched MeSH terms: Finite Element Analysis
  15. Yanliang Shang, Shouji Du, Biao Shao, Tongyin Han
    Sains Malaysiana, 2017;46:2091-2099.
    A large number of shallow buried tunnels are built in the city nowadays and the special strata such as large upper-soft and lower-hard ground often encountered. Deformation control of strata is the focus issue related to the construction safety. Based on Dalian metro Hing Street station with the classical geological condition of upper-soft and lower-hard ground, this paper fully used a combined control method including six different support measures to control the deformation of surrounding rock. 3D finite element model was setup to analyze the construction effect of combined control measures and the monitoring in-site was carried out to verify the deformation control effect of combined control method. It shows that the maximum surface subsidence value is gradually reduced with the support measures gradually increasing. In the case of various supports the maximum sedimentation value is 2.67 cm, which is 42. 1% lower than that of not using control method and the control effect is obvious. In addition, it can be seen that the two-layer initial support and additional large arch foot have the best effect on controlling the ground surface settlement with reduction of 11.7% and 20.2%, respectively. The research results can provide practical experience for the construction of such tunnels, and guide the design and construction of the tunnel in the future.
    Matched MeSH terms: Finite Element Analysis
  16. Bajuri MN, Isaksson H, Eliasson P, Thompson MS
    Biomech Model Mechanobiol, 2016 12;15(6):1457-1466.
    PMID: 26951049
    The healing process of ruptured tendons is problematic due to scar tissue formation and deteriorated material properties, and in some cases, it may take nearly a year to complete. Mechanical loading has been shown to positively influence tendon healing; however, the mechanisms remain unclear. Computational mechanobiology methods employed extensively to model bone healing have achieved high fidelity. This study aimed to investigate whether an established hyperelastic fibre-reinforced continuum model introduced by Gasser, Ogden and Holzapfel (GOH) can be used to capture the mechanical behaviour of the Achilles tendon under loading during discrete timepoints of the healing process and to assess the model's sensitivity to its microstructural parameters. Curve fitting of the GOH model against experimental tensile testing data of rat Achilles tendons at four timepoints during the tendon repair was used and achieved excellent fits ([Formula: see text]). A parametric sensitivity study using a three-level central composite design, which is a fractional factorial design method, showed that the collagen-fibre-related parameters in the GOH model-[Formula: see text] and [Formula: see text]-had almost equal influence on the fitting. This study demonstrates that the GOH hyperelastic fibre-reinforced model is capable of describing the mechanical behaviour of healing tendons and that further experiments should focus on establishing the structural and material parameters of collagen fibres in the healing tissue.
    Matched MeSH terms: Finite Element Analysis
  17. Jaafar, M.S., Wong, N.J., Noorzaei, J., Thanoon, W.A.
    MyJurnal
    This study presents an efficient finite element analysis technique which shows great versatility in
    modelling of precast composite flooring system subjected to static loadings. The method incorporates sliding and opening in the analysis of composite structures using the interface element which was specifically designed to simulate the actual behaviour at the interfaces between contacting materials. A three-dimensional finite element model of the precast composite slab which exhibits discontinuous behaviour was performed to demonstrate the potential and applicability of the proposed method of analysis. The results of the analysis demonstrate that the overall response of a discontinuous system to external loading is significantly affected by the bonding condition at the interfaces between the contacting materials.
    Matched MeSH terms: Finite Element Analysis
  18. Lim LG, Pao WK, Hamid NH, Tang TB
    Sensors (Basel), 2016 Jul 04;16(7).
    PMID: 27384567 DOI: 10.3390/s16071032
    A 360° twisted helical capacitance sensor was developed for holdup measurement in horizontal two-phase stratified flow. Instead of suppressing nonlinear response, the sensor was optimized in such a way that a 'sine-like' function was displayed on top of the linear function. This concept of design had been implemented and verified in both software and hardware. A good agreement was achieved between the finite element model of proposed design and the approximation model (pure sinusoidal function), with a maximum difference of ±1.2%. In addition, the design parameters of the sensor were analysed and investigated. It was found that the error in symmetry of the sinusoidal function could be minimized by adjusting the pitch of helix. The experiments of air-water and oil-water stratified flows were carried out and validated the sinusoidal relationship with a maximum difference of ±1.2% and ±1.3% for the range of water holdup from 0.15 to 0.85. The proposed design concept therefore may pose a promising alternative for the optimization of capacitance sensor design.
    Matched MeSH terms: Finite Element Analysis
  19. Masni-Azian, Tanaka M
    Comput Biol Med, 2018 07 01;98:26-38.
    PMID: 29758454 DOI: 10.1016/j.compbiomed.2018.05.010
    Intervertebral disc degeneration involves changes in its material properties that affect the mechanical functions of the spinal system. However, the alteration of the biomechanics of a spinal segment through specific material degradation in a specific region is poorly understood. In this study, the influence of the constitutive material degeneration of disc tissues on the mechanics of a lower lumbar spinal unit was examined using a three-dimensional nonlinear finite element model of the L4-L5 functional spinal unit. Different grades of disc degeneration were simulated by introducing a degeneration factor to the corresponding material properties to represent fibrous nucleus, increased fibre and ground substance laxity, increased fibre stiffness and total annular fracture along posterior and posterolateral regions. The model was loaded with an axial compression of 500 N and pure moments of up to 10 Nm to simulate extension, flexion, lateral bending and axial rotation. To validate the model, the spinal motion and intradiscal pressure of healthy and degenerated discs with existing in vitro data were compared. The disc with a fibrous nucleus and the presence of intradiscal pressure increase the spinal instability during flexion and axial rotation, and the absence of intradiscal pressure increases the spinal instability in all directions. Bulging displacement and shear strains in the disc with total fracture and ground substance laxity are high in all of the loading cases. Our study could provide useful information to enhance our understanding of the influence of each constitutive component of the intervertebral disc on the mechanics of the spinal segment.
    Matched MeSH terms: Finite Element Analysis
  20. Masni-Azian, Tanaka M
    Comput Methods Biomech Biomed Engin, 2017 Aug;20(10):1066-1076.
    PMID: 28532164 DOI: 10.1080/10255842.2017.1331345
    In the biomechanics field, material parameters calibration is significant for finite element (FE) model to ensure a legit estimation of biomechanical response. Determining an appropriate combination of calibration factors is challenging as each constitutive component responds differently. This study proposes a statistical factorial analysis approach using L16(4(5)) orthogonal array to evaluate material nonlinearity and applicable calibration factor of the intervertebral disc FE model in pure moment. The calibrated model exhibits improved agreement to the experimental findings for all directions. Appropriate combination of calibration parameter reduces the estimation gap to the experimental findings, ensuring agreeable biomechanical responses.
    Matched MeSH terms: Finite Element Analysis*
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links