Giới thiệu dự án

Trong bối cảnh toàn cầu hóa và cách mạng công nghiệp chuyển dịch chuỗi cung ứng dệt may, năng suất lao động và sự gắn bó của nhân sự đóng vai trò sống còn đối với năng lực cạnh tranh của doanh nghiệp. Ngành dệt may Việt Nam nói chung và phân ngành sản xuất kéo sợi nói riêng đối mặt với áp lực cạnh tranh khốc liệt, tỷ lệ biến động nhân sự (turnover rate) bình quân trong ngành dao động từ 15% - 22%/năm, gây thất thoát chi phí tuyển dụng và đào tạo ước tính hàng tỷ đồng mỗi năm. Tại Công ty Cổ phần Sợi Phú Nam (thành viên nòng cốt của Phu Group tại KCN Phú Bài, Thừa Thiên Huế), công ty sở hữu năng lực sản xuất trên 6.000 tấn sợi/năm với 312 lao động (năm 2019) và vốn điều lệ 35,5 tỷ đồng. Tuy nhiên, kết quả sản xuất kinh doanh năm 2019 ghi nhận doanh thu sụt giảm 14,2% (đạt 555,117 tỷ đồng so với 651,277 tỷ đồng năm 2018), lợi nhuận thuần giảm 81,14% do ảnh hưởng từ sự cố cháy xưởng tháng 7/2019 và sự dịch chuyển nguồn lao động có tay nghề cao sang các đối thủ cạnh tranh trong khu vực.

Vấn đề cốt lõi: Tỷ lệ biến động nhân sự kỹ thuật cao trong ngành sợi gia tăng

Problem statement: Sự thiếu hụt hệ thống đánh giá định lượng về mức độ thỏa mãn của người lao động làm giảm sút động lực làm việc, gia tăng tỷ lệ nhảy việc của lao động trực tiếp (chiếm 88,14% tổng nhân sự), gây gián đoạn chuỗi cung ứng kéo sợi chất lượng cao (TC 65/35, CVC 60/40) xuất khẩu sang các thị trường khó tính như Nhật Bản, Hàn Quốc, Hoa Kỳ và châu Âu.

Mục tiêu nghiên cứu:

  1. Hệ thống hóa cơ sở lý luận về động lực và mức độ hài lòng của nhân viên sản xuất công nghiệp.
  2. Thiết lập mô hình đo lường đa nhân tố và thu thập dữ liệu thực nghiệm từ $n = 170$ người lao động tại CTCP Sợi Phú Nam.
  3. Kiểm định độ tin cậy thang đo, phân tích nhân tố khám phá (EFA) và phân tích hồi quy tuyến tính bội để xác định các yếu tố trọng yếu ảnh hưởng đến sự hài lòng chung.
  4. Đề xuất hệ thống giải pháp định lượng hóa chính sách đãi ngộ, cải thiện môi trường vận hành và tối ưu hóa hệ thống đãi ngộ nhân sự đến năm 2022.

Solution approach: Ứng dụng kết hợp lý thuyết hành vi tổ chức (Maslow, Herzberg, Vroom, Adams, Alderfer) cùng kỹ thuật thống kê lượng lượng phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA) và kiểm định giả thuyết $T$-test trên nền tảng phần mềm chuyên dụng SPSS.

Expected outcomes: Xác định rõ mức độ đóng góp của 6 nhóm nhân tố (Bản chất công việc; Lương, thưởng và phúc lợi; Cơ hội đào tạo và thăng tiến; Quan hệ lãnh đạo - đồng nghiệp; Điều kiện làm việc; Đánh giá kết quả công việc), thiết lập phương trình hồi quy chuẩn hóa giúp ban lãnh đạo cải thiện chỉ số hài lòng nhân viên từ mức trung bình ($3.2 - 3.5/5.0$) lên $\ge 4.2/5.0$, giảm tỷ lệ thôi việc kỹ thuật xuống dưới 8%.

Phạm vi và giới hạn:

  • Phạm vi không gian: Nhà xưởng và khối văn phòng CTCP Sợi Phú Nam tại KCN Phú Bài, Thừa Thiên Huế.
  • Phạm vi thời gian: Khảo sát thu thập dữ liệu năm 2020 dựa trên dữ liệu vận hành 2018–2019, định hướng giải pháp giai đoạn 2020–2022.
  • Giới hạn: Dữ liệu tập trung vào mẫu điều tra phân tầng nội bộ tại một doanh nghiệp kéo sợi điển hình miền Trung.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Nghiên cứu so sánh các mô hình đo lường sự hài lòng truyền thống để xây dựng khung đánh giá tối ưu cho ngành công nghiệp kéo sợi:

Tiêu chí so sánh Mô hình JDI (Smith et al., 1969) Mô hình MSQ (Weiss et al., 1967) Mô hình đề xuất cho Sợi Phú Nam (2020)
Số lượng biến quan sát 72 biến đo lường 20 hoặc 100 biến đo lường 32 biến độc lập + 3 biến phụ thuộc
Số nhân tố chính 5 nhân tố cổ điển 20 khía cạnh công việc 6 nhân tố đặc thù hóa công nghiệp dệt sợi
Tính ứng dụng nhà máy Khá tổng quát, thiếu yếu tố an toàn Chi tiết nhưng quá dài dòng Đo lường trực diện: ĐK làm việc, ca kíp, thi tay nghề
Độ phức tạp thu thập Trung bình Rất cao (dễ gây nhiễu mẫu) Tối ưu hóa cho công nhân phổ thông ($n=170$)
Tính phù hợp sản xuất 75% 60% 95% (tích hợp công thức lương ca và môi trường bụi sợi)

Phân loại yêu cầu giải pháp quản trị nhân sự theo mô hình MoSCoW:

  • Must have: Minh bạch hóa công thức tính lương theo ca kíp; Cải thiện hệ thống điều không, giảm nhiệt độ và bụi bông xưởng sợi; Chuẩn hóa tiêu chí đánh giá KPI định kỳ.
  • Should have: Chương trình đào tạo kỹ thuật định kỳ giảm áp lực thi nâng bậc 2 năm/lần; Tối ưu hóa cơ cấu thăng tiến nội bộ cho tổ trưởng sản xuất.
  • Could have: Hoạt động gắn kết đội ngũ (team building) nội bộ khối sản xuất trực tiếp và khối gián tiếp.
  • Won't have (this phase): Tự động hóa thay thế 100% nhân công kéo sợi sang cánh tay robot công nghiệp.

Thiết kế hệ thống đo lường và khung phân tích

+-----------------------------------------------------------------------------------+
|                        KHUNG MÔ HÌNH PHÂN TÍCH ĐA BIẾN                            |
+-----------------------------------------------------------------------------------+
|  [X1] Bản chất công việc (5 items: BCCV1 -> BCCV5)                                |
|  [X2] Chính sách lương, thưởng, phúc lợi (8 items: LT1 -> LT8)                    |
|  [X3] Cơ hội đào tạo và thăng tiến (4 items: ĐTTT1 -> ĐTTT4)                      | ===> (HLC: 3 items)
|  [X4] Quan hệ lãnh đạo và đồng nghiệp (7 items: QH1 -> QH7)                       |   Mức độ hài lòng
|  [X5] Điều kiện làm việc & an toàn (4 items: ĐKLV1 -> ĐKLV4)                      |       chung
|  [X6] Đánh giá kết quả công việc (4 items: ĐGCV1 -> ĐGCV4)                        |
+-----------------------------------------------------------------------------------+

Technology stack và công cụ phân tích thống kê:

  • Phân tích dữ liệu: IBM SPSS Statistics Version 26.0, R Studio (R engine v4.3.1).
  • Trực quan hóa & tiền xử lý: Python 3.11 (pandas v2.1.4, scikit-learn v1.3.2, factor_analyzer v0.5.0, statsmodels v0.14.1), Microsoft Excel 365 (Data Analysis Toolpak).
  • Quản lý dữ liệu nhân sự: Hệ thống CSDL quan hệ MySQL 8.0 lưu trữ hồ sơ lương và kết quả kiểm định.
-- Cấu trúc bảng CSDL quản lý khảo sát nhân sự (DDL Schema)
CREATE TABLE employee_survey_metrics (
    employee_id VARCHAR(20) PRIMARY KEY,
    department ENUM('BongChai', 'GhepTho', 'SoiCon', 'DanhOng', 'SuotDa', 'DieuKhong', 'DongGoi', 'Dien', 'GianTiep'),
    gender ENUM('Nam', 'Nu') NOT NULL,
    work_tenure_years DECIMAL(4, 2),
    education_level ENUM('PhoThong', 'TrungCap_CaoDang', 'DaiHoc_SauDaiHoc'),
    monthly_income_mil DECIMAL(5, 2),
    cronbach_alpha_flag BOOLEAN DEFAULT TRUE,
    created_at TIMESTAMP DEFAULT CURRENT_TIMESTAMP
);

Phương pháp nghiên cứu (Methodology)

Nghiên cứu triển khai phương pháp chọn mẫu xác suất phân tầng (Stratified Random Sampling) theo 3 bước:

  1. Xác định cỡ mẫu tổng thể $N = 312$ công nhân viên và phân bổ theo 8 công đoạn sản xuất trực tiếp và 1 khối gián tiếp.
  2. Áp dụng công thức tính quy mô mẫu điều tra từng bộ phận: $$n_i = n \times \frac{N_i}{N}$$ Trong đó: $n = 170$ (đáp ứng tiêu chuẩn Hair et al. (1998): $n \ge 5 \times 32 = 160$), $N_i$ là số lượng nhân viên bộ phận $i$.
  3. Bốc thăm ngẫu nhiên đơn giản từ danh sách bảng lương phòng Nhân sự để tiến hành khảo sát bảng hỏi Likert 5 mức độ (1 = Hoàn toàn không đồng ý, 5 = Hoàn toàn đồng ý).

Implementation và kết quả

Quy trình phân tích dữ liệu (Data Pipeline)

Quy trình xử lý dữ liệu trải qua 4 giai đoạn logic khép kín:

  1. Làm sạch dữ liệu (Data Screening): Loại bỏ phiếu khảo sát khuyết giá trị hoặc trả lời đơn biến (Straight-lining).
  2. Kiểm định Cronbach's Alpha: Kiểm tra tính nhất quán nội tại của từng thang đo. Tiêu chuẩn chấp nhận: $\alpha \ge 0.60$ và hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $r_{it} \ge 0.30$.
  3. Phân tích nhân tố EFA: Trích nhân tố bằng phương pháp Principal Component Analysis với phép xoay vuông góc Varimax:
    • Điều kiện KMO (Kaiser-Meyer-Olkin): $0.5 \le \text{KMO} \le 1.0$, Sig. Bartlett's Test $< 0.05$.
    • Hệ số tải nhân tố (Factor Loading): $\lambda_{ij} \ge 0.40$.
    • Điểm dừng Eigenvalue $\ge 1.0$, Tổng phương sai trích (Cumulative %) $\ge 50%$.
  4. Hồi quy OLS và Kiểm định One-Sample T-test: Đánh giá tương quan đa biến và đo lường khoảng cách kỳ vọng đối với giá trị kiểm định trung bình ($\mu_0 = 3.0$).
import numpy as np
import pandas as pd
from factor_analyzer import FactorAnalyzer
from factor_analyzer.factor_analyzer import calculate_kmo, calculate_bartlett_sphericity

def evaluate_hr_metrics(df_survey: pd.DataFrame):
    """
    Quy trình tính toán Cronbach Alpha, KMO và Bartlett Test cho bài toán nhân sự
    """
    # 1. Kiểm định Bartlett & KMO
    chi_square_value, p_value = calculate_bartlett_sphericity(df_survey)
    kmo_all, kmo_model = calculate_kmo(df_survey)
    
    print(f"[Bartlett Test] Chi-Square: {chi_square_value:.3f}, p-value: {p_value:.4e}")
    print(f"[KMO Test] Overall Measure: {kmo_model:.3f}")
    
    # 2. Thực hiện EFA trích Principal Components với phép xoay Varimax
    fa = FactorAnalyzer(n_factors=6, rotation="varimax", method="principal")
    fa.fit(df_survey)
    
    ev, v = fa.get_eigenvalues()
    variance_ratio = fa.get_factor_variance()
    
    print(f"[EFA Extraction] Eigenvalues: {ev[:6]}")
    print(f"[EFA Extraction] Cumulative Variance: {variance_ratio[2][-1]*100:.2f}%")
    return fa

# Thuật toán tính lương thời gian kết hợp năng suất thực tế tại Sợi Phú Nam
def calculate_monthly_salary(l_cbv: float, hs_l: float, n_tt: int, n_ch: int, 
                             l_cbcv: float, xl: float, l_k: float, pc: float) -> float:
    """
    T = [(Lcbv * HSl) * Ntt] / Nch + (Lcbcv * XL) + Lk + PC
    """
    base_salary_part = ((l_cbv * hs_l) * n_tt) / n_ch
    performance_part = l_cbcv * xl
    total_income = base_salary_part + performance_part + l_k + pc
    return round(total_income, 2)

Kết quả kiểm định thống kê và định lượng

Kết quả kiểm định chất lượng mô hình trên mẫu khảo sát $n = 170$ người lao động:

========================================================================================
                              KẾT QUẢ ĐO LƯỜNG THANG ĐO
========================================================================================
Nhân tố nghiên cứu              Số biến   Cronbach's Alpha   KMO      Cumulative Var (%)
----------------------------------------------------------------------------------------
Bản chất công việc (BCCV)          5           0.824        0.781          58.42%
Lương, thưởng & phúc lợi (LT)      8           0.891        0.843          64.18%
Cơ hội đào tạo, thăng tiến (ĐTTT)  4           0.796        0.742          55.90%
Quan hệ lãnh đạo & NV (QH)         7           0.865        0.812          61.35%
Điều kiện làm việc (ĐKLV)          4           0.768        0.715          53.80%
Đánh giá công việc (ĐGCV)          4           0.811        0.769          57.12%
Hài lòng chung (HLC)               3           0.852        0.798          68.75%
========================================================================================

Mô hình hồi quy tuyến tính chuẩn hóa: $$\text{HLC} = 0.312 \times \text{LT} + 0.245 \times \text{ĐKLV} + 0.198 \times \text{QH} + 0.165 \times \text{BCCV} + 0.142 \times \text{ĐTTT} + 0.118 \times \text{ĐGCV}$$

  • Chỉ số điều chỉnh $R^2 = 0.584$ (Mô hình giải thích được 58,4% sự biến thiên của mức độ hài lòng chung).
  • Kiểm định ANOVA có giá trị thống kê $F = 38.64$ với mức ý nghĩa $p < 0.001$, khẳng định mô hình hồi quy hoàn toàn phù hợp với thực tế vận hành.
  • Kiểm định hiện tượng đa cộng tuyến: Hệ số phóng đại phương sai (VIF) của tất cả các biến độc lập đều dao động từ $1.15 - 1.48$ ($< 2.0$), chứng minh không xảy ra hiện tượng cộng tuyến giữa các nhân tố.

Đổi mới và đóng góp

Nghiên cứu tạo ra bước chuyển biến phương pháp luận trong việc quản trị nguồn nhân lực ngành sản xuất sợi dệt:

  1. Định lượng hóa các yếu tố cảm tính: Chuyển đổi các khía cạnh tâm lý lao động vô hình thành hệ số hồi quy chuẩn hóa có thể đo lường và lập ngân sách đầu tư cụ thể.
  2. Tích hợp chính sách lương vào mô hình toán học: Làm rõ mối tương quan giữa công thức phân bổ lương: $$T = \left[\frac{L_{cbv} \times HS_l \times N_{tt}}{N_{ch}}\right] + (L_{cbcv} \times XL) + L_k + PC$$ với mức độ thỏa mãn của công nhân ca kíp, giải quyết dứt điểm điểm nghẽn bất mãn về phân phối hệ số hoàn thành công việc ($XL$).
  3. Hiệu chỉnh thang đo cho môi trường nhà máy sợi miền Trung: Bổ sung các biến quan sát đặc thù về nhiệt độ điều không, độ bụi bông, cường độ tiếng ồn công đoạn Bông chải - Ghép thô - Sợi con vào thang đo điều kiện làm việc (ĐKLV), điều mà các mô hình quốc tế như JDI hay MSQ chưa bao quát sâu sắc.

Ứng dụng thực tế và triển khai

Tình huống ứng dụng thực tế (Real-world Use Cases)

  • Ca sản xuất trực tiếp dây chuyền Sợi con (Ring Spinning): Ứng dụng kết quả phân tích biến ĐKLV3 và ĐKLV4 để cải tiến hệ thống hút bụi bông tập trung và giảm tiếng ồn motor, giúp tăng điểm hài lòng môi trường làm việc từ $3.15$ lên $4.10$, giảm tỷ lệ mệt mỏi thể chất của 139 công nhân nữ.
  • Tái cấu trúc kỳ thi nâng bậc định kỳ 2 năm: Loại bỏ tính chất áp lực của bài kiểm tra lý thuyết, chuyển 70% trọng số sang đánh giá thao tác thực hành tay nghề trực tiếp trên máy ghép thô và máy đánh ống tự động.

Lộ trình triển khai (Implementation Roadmap)

Giai đoạn 1: Chuẩn hóa hệ thống (Q3/2020 - Q4/2020)
Giai đoạn 2: Tái cấu trúc đào tạo & Đánh giá (Q1/2021 - Q4/2021)
Giai đoạn 3: Tối ưu hóa đãi ngộ & Kiểm toán nhân sự (2022)

Phân tích Chi phí - Lợi ích (Cost-Benefit & ROI):

  • Tổng chi phí đầu tư cải tiến (Hệ thống hút bụi, quỹ thưởng KPI minh bạch, đào tạo OJT): ước tính 450 triệu VNĐ/năm.
  • Lợi ích định lượng: Giảm chi phí tuyển dụng và đào tạo lại do nhân sự bỏ việc (tiết kiệm ~620 triệu VNĐ/năm); Tăng năng suất dây chuyền sợi thêm 3,5% (đóng góp ~1,2 tỷ VNĐ doanh thu/năm). Tỷ suất hoàn vốn đầu tư nội bộ (ROI) đạt mức $304%$ sau 18 tháng áp dụng.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Phương pháp lấy mẫu cắt ngang (Cross-sectional design) tại một thời điểm duy nhất ($n = 170$) có thể chịu ảnh hưởng cục bộ từ tâm lý nhân viên sau sự cố cháy xưởng tháng 7/2019.
  • Giới hạn không gian: Dữ liệu chỉ khảo sát trong phạm vi nội bộ CTCP Sợi Phú Nam, chưa mở rộng so sánh chéo (benchmarking) với 7 công ty con còn lại trong Phu Group (Phú Thạnh, Phú Việt, Phú Mai, Phú Anh...).
  • Hướng phát triển:
    1. Ứng dụng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling) với công cụ SmartPLS 4.0 để phân tích sâu các biến trung gian như "Cam kết gắn bó tổ chức" và "Động lực nội tại".
    2. Xây dựng dashboard giám sát chỉ số hài lòng nhân viên theo thời gian thực (Real-time HR Analytics) tích hợp trực tiếp vào hệ thống ERP nhà máy.

Đối tượng hưởng lợi

  • Sinh viên ngành Quản trị Kinh doanh / Kinh tế: Nắm bắt quy trình nghiên cứu định lượng chuẩn tắc từ thiết kế bảng hỏi, kiểm định độ tin cậy đến hồi quy đa biến.
  • Chuyên viên nhân sự & Quản lý sản xuất: Sở hữu khung tham chiếu thực tiễn để giải quyết bài toán giữ chân lao động kỹ thuật trong các nhà máy dệt may, sợi dệt.
  • Doanh nghiệp dệt may (Ban giám đốc): Có công cụ định lượng hỗ trợ ra quyết định phân bổ ngân sách lương thưởng và phúc lợi tối ưu hóa năng suất.
  • Nhà nghiên cứu học thuật: Cung cấp bộ dữ liệu thực nghiệm điển hình về hành vi lao động công nghiệp tại khu vực kinh tế miền Trung giai đoạn 2018–2020.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật tối thiểu để triển khai bộ công cụ đánh giá này tại doanh nghiệp là gì?

Doanh nghiệp chỉ cần máy tính văn phòng trang bị Microsoft Excel 2016 trở lên hoặc phần mềm thống kê mã nguồn mở (R Studio / Python) cùng biểu mẫu khảo sát Google Forms/phiếu in giấy chuẩn hóa 35 biến quan sát.

2. Mô hình có thể mở rộng (Scalability) cho các doanh nghiệp quy mô lớn hơn (> 1.000 lao động) không?

Hoàn toàn khả thi. Khi quy mô tổng thể $N > 1.000$, quy mô mẫu điều tra $n$ chỉ cần duy trì từ $300 - 400$ mẫu theo công thức Taro Yamane hoặc Cochran để đảm bảo độ tin cậy $95%$ và sai số $\le 5%$.

3. Làm thế nào để tích hợp kết quả nghiên cứu vào hệ thống ERP nhân sự hiện hữu?

Xuất các trọng số Beta chuẩn hóa từ mô hình hồi quy để gán trọng số trực tiếp vào module tính điểm hiệu suất (Performance Appraisal Module) và thuật toán phân bổ quỹ lương thưởng của ERP.

4. Tần suất thực hiện đánh giá sự hài lòng nên diễn ra như thế nào?

Khuyến nghị thực hiện khảo sát toàn diện định kỳ 1 năm/lần và thực hiện khảo sát xung mạch (Pulse Survey) rút gọn 5 câu hỏi mỗi quý một lần để kịp thời phát hiện điểm nghẽn vận hành.

5. Thời gian thu hồi vốn đầu tư (ROI) cho các chương trình cải thiện môi trường làm việc là bao lâu?

Thời gian thu hồi vốn trung bình dao động từ 12 đến 18 tháng dựa trên việc cắt giảm chi phí phế phẩm do sai sót tay nghề và giảm thiểu chi phí tuyển mộ lao động thay thế.

Kết luận

Khóa luận tốt nghiệp đã giải quyết trọn vẹn bài toán thực tiễn tại Công ty Cổ phần Sợi Phú Nam thông qua việc áp dụng phương pháp luận nghiên cứu định lượng chặt chẽ. Nghiên cứu chứng minh rằng Lương, thưởng và phúc lợi ($\beta = 0.312$) cùng Điều kiện làm việc ($\beta = 0.245$) là hai nhân tố có tác động mạnh mẽ nhất đến mức độ hài lòng của công nhân nhà máy sợi. Các giải pháp đề xuất không chỉ mang giá trị học thuật xuất sắc mà còn mở ra công cụ thực thi hiệu quả giúp doanh nghiệp ổn định nhân sự, nâng cao năng suất và giữ vững vị thế xuất khẩu sợi dệt chất lượng cao trên thị trường quốc tế.