Giới thiệu dự án

Trong bối cảnh công nghiệp dệt may Việt Nam đối mặt với sự cạnh tranh khốc liệt và biến động nhân sự lớn, việc duy trì và nâng cao hiệu suất lao động trở thành yếu tố sống còn của doanh nghiệp. Theo các khảo sát ngành dệt may, tỷ lệ luân chuyển lao động (turnover rate) trung bình tại các nhà máy sản xuất sợi thường dao động từ 15% đến 25%/năm, gây áp lực trực tiếp lên chi phí tuyển dụng và năng suất dây chuyền. Đề tài khóa luận nghiên cứu ứng dụng: "Phân tích các nhân tố ảnh hưởng đến động lực làm việc của người lao động tại Công ty Cổ phần Vinatex Phú Hưng" do tác giả Nguyễn Thị Như Ý thực hiện (GVHD: ThS. Hà Ngọc Thùy Liên, Trường Đại học Kinh tế – Đại học Huế) giải quyết bài toán cốt lõi về quản trị nguồn nhân lực (Human Resource Management - HRM) trong môi trường sản xuất sợi công nghiệp.

+-----------------------------------------------------------------------+
|                         BỐI CẢNH DOANH NGHIỆP                         |
|   Công ty CP Vinatex Phú Hưng (KCN Phú Bài, Hương Thủy, Thừa Thiên Huế)|
|   Quy mô nhân sự: N = 225 lao động | Ngành nghề cốt lõi: Sản xuất sợi |
+-----------------------------------------------------------------------+
                                   |
                                   v
+-----------------------------------------------------------------------+
|                           VẤN ĐỀ CẦN XỬ LÝ                            |
|  - Biến động động lực lao động theo ca kíp sản xuất nặng nhọc         |
|  - Chính sách đãi ngộ, tiền lương chưa tối ưu hóa theo hiệu suất       |
|  - Môi trường nhà xưởng công nghiệp (tiếng ồn, bụi bông, nhiệt độ)   |
+-----------------------------------------------------------------------+
                                   |
                                   v
+-----------------------------------------------------------------------+
|                         GIẢI PHÁP MÔ HÌNH HÓA                         |
|  Mô hình 7 nhân tố định lượng (34 biến quan sát) - Cronbach's Alpha   |
|  Phân tích nhân tố khám phá (EFA) & Hồi quy tuyến tính bội (OLS)      |
+-----------------------------------------------------------------------+

Vấn đề nghiên cứu thực tiễn

Công ty Cổ phần Vinatex Phú Hưng (đóng tại KCN Phú Bài, Thừa Thiên Huế) có quy mô $N = 225$ cán bộ công nhân viên, chuyên sản xuất sợi và nguyên phụ liệu dệt may. Đặc thù công việc đòi hỏi làm việc 3 ca liên tục, cường độ cao và chịu tác động của môi trường bụi xơ, tiếng ồn máy kéo sợi. Thực tế ghi nhận động lực lao động của công nhân trực tiếp có dấu hiệu suy giảm, biểu hiện qua tỷ lệ nghỉ việc cục bộ, năng suất chưa đạt mức tối đa và mức độ thỏa mãn với chính sách đãi ngộ chưa đồng đều.

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

  1. Hệ thống hóa cơ sở lý luận về động cơ, nhu cầu và động lực làm việc thông qua việc tích hợp Thuyết cấp bậc nhu cầu Maslow (1943), Thuyết công bằng Stacy Adams (1963) và Thuyết kỳ vọng Victor Vroom (1964).
  2. Đo lường và đánh giá thực trạng các yếu tố quản trị nhân lực ảnh hưởng đến động lực lao động tại doanh nghiệp giai đoạn 2016 – 2018.
  3. Xác định mức độ tác động ($Beta$) của 7 nhóm nhân tố đến động lực làm việc thông qua phương pháp định lượng hồi quy tuyến tính bội OLS (Ordinary Least Squares).
  4. Đề xuất hệ thống giải pháp can thiệp có tính khả thi cao nhằm cải thiện chính sách tiền lương, điều kiện làm việc và cơ chế quản lý.

Phạm vi và giới hạn nghiên cứu

  • Không gian: Công ty Cổ phần Vinatex Phú Hưng, KCN Phú Bài, Hương Thủy, Thừa Thiên Huế.
  • Thời gian: Dữ liệu thứ cấp tổng hợp giai đoạn 2016 – 2018; Dữ liệu sơ cấp khảo sát thực địa từ tháng 09/2019 đến tháng 12/2019.
  • Đối tượng: Toàn bộ công nhân trực tiếp và lao động gián tiếp đang làm việc tại doanh nghiệp.

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

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

Để lượng hóa hành vi tổ chức, nghiên cứu thực hiện rà soát các mô hình tạo động lực kinh điển và hiện đại trên thế giới cũng như tại Việt Nam, từ đó xác lập khoảng trống nghiên cứu (gap analysis).

Nghiên cứu / Mô hình Đối tượng & Phạm vi Ưu điểm Hạn chế
Kennett S. Kovach (1987) 10 nhân tố tạo động lực trong ngành công nghiệp tổng quát Toàn diện, bao quát cả động lực tài chính và phi tài chính Thiếu biến số đặc thù về văn hóa doanh nghiệp tại thị trường Á Đông
Marko Kukanja (2013) 191 nhân viên ngành dịch vụ du lịch (Slovenia) Xác định rõ vai trò trọng tâm của tiền lương và thời gian làm việc Đặc thù ngành dịch vụ khác biệt lớn so với dây chuyền sản xuất sợi
Bùi Thị Minh Thu & CS (2014) 223 công nhân ngành cơ khí lắp máy Việt Nam Kiểm định chặt chẽ EFA, phù hợp môi trường sản xuất kỹ thuật Chưa đi sâu vào biến động động lực theo ca kíp 3 ca/ngày
Mô hình đề xuất Vinatex Phú Hưng 180 lao động ngành sản xuất sợi dệt may Tích hợp 7 nhân tố toàn diện, thiết kế riêng cho công nhân dệt may Cần kiểm chứng định kỳ khi quy mô doanh nghiệp mở rộng

Áp dụng phương pháp phân tích yêu cầu MoSCoW trong thiết kế thang đo:

  • Must have: Thang đo Lương thưởng - Phúc lợi ($LTPL$), Điều kiện làm việc ($DKLV$), Đặc điểm công việc ($DDCV$).
  • Should have: Phong cách lãnh đạo ($PCLD$), Đào tạo & Thăng tiến ($DTTT$).
  • Could have: Quan hệ đồng nghiệp ($DN$), Văn hóa công ty ($VHCT$).
  • Won't have (kỳ này): Đo lường áp lực gia đình và biến động kinh tế vĩ mô ngoại cảnh.

Thiết kế hệ thống nghiên cứu định lượng

Mô hình nghiên cứu bao gồm 7 biến độc lập tác động đến 1 biến phụ thuộc thông qua 34 biến quan sát đo lường trên thang đo Likert 5 điểm (1: Rất không đồng ý $\rightarrow$ 5: Rất đồng ý).

+---------------------------+       +---------------------------+
|    Đặc điểm công việc     |       |    Điều kiện làm việc     |
|          (DDCV)           |       |          (DKLV)           |
+-------------+-------------+       +-------------+-------------+
              |                                   |
+-------------v-------------+       +-------------v-------------+
|  Lương, thưởng, phúc lợi  |       |    Đào tạo & thăng tiến   |
|          (LTPL)           +------>|          (DTTT)           |
+-------------+-------------+  (H)  +-------------+-------------+
              |               \                   |
              |                v                  |
              |     +-----------------------+     |
              +---->|  ĐỘNG LỰC LÀM VIỆC    |<----+
              |     |        (DLLV)         |     |
              |     +-----------------------+     |
              |                ^                  |
+-------------v-------------+ /     +-------------v-------------+
|    Phong cách lãnh đạo    +<      |         Đồng nghiệp       |
|          (PCLD)           |       |            (DN)           |
+-------------+-------------+       +-------------+-------------+
              |                                   |
+-------------v-----------------------------------v-------------+
|                     Văn hóa công ty (VHCT)                    |
+---------------------------------------------------------------+

Mô hình toán học hồi quy tuyến tính đa biến OLS được thiết lập:

$$DLLV = \beta_0 + \beta_1 \cdot DDCV + \beta_2 \cdot DKLV + \beta_3 \cdot LTPL + \beta_4 \cdot DTTT + \beta_5 \cdot PCLD + \beta_6 \cdot DN + \beta_7 \cdot VHCT + e_i$$

Trong đó:

  • $DLLV$: Biến phụ thuộc - Động lực làm việc của người lao động.
  • $\beta_0$: Hằng số tự do (Constant intercept).
  • $\beta_1 \dots \beta_7$: Hệ số hồi quy chuẩn hóa đại diện cho trọng số tác động của từng nhân tố.
  • $e_i$: Sai số ngẫu nhiên của mô hình.

Methodology (Phương pháp luận)

  1. Xác định cỡ mẫu nghiên cứu: Dựa trên công thức xác định cỡ mẫu của Taro Yamane (1973):

$$n = \frac{N}{1 + N \cdot (e)^2} = \frac{225}{1 + 225 \cdot (0.05)^2} \approx 144 \text{ mẫu}$$

Theo tiêu chuẩn phân tích nhân tố EFA của Hair et al. (1998) và Hoàng Trọng & Chu Nguyễn Mộng Ngọc (2008), cỡ mẫu tối thiểu cần gấp 5 lần số biến quan sát:

$$N_{\text{min}} \ge 34 \times 5 = 170 \text{ mẫu}$$

Để đảm bảo độ tin cậy và bù trừ bảng hỏi lỗi, nghiên cứu phát ra 190 phiếu, thu về 180 phiếu hợp lệ (đáp ứng trọn vẹn cả hai điều kiện).

  1. Quy trình phân tích dữ liệu 4 tầng:
    • Tầng 1: Thống kê mô tả cấu trúc nhân khẩu học (Demographics).
    • Tầng 2: Kiểm định độ tin cậy thang đo bằng Cronbach's Alpha ($\ge 0.60$, tương quan biến - tổng $\ge 0.30$).
    • Tầng 3: Phân tích nhân tố khám phá EFA (KMO $\ge 0.50$, Bartlett's Test $p < 0.05$, Factor Loading $\ge 0.50$, Phương sai trích $\ge 50%$, Eigenvalue $> 1$).
    • Tầng 4: Kiểm định giả định hồi quy (Đa cộng tuyến VIF, Tự tương quan Durbin-Watson) và kiểm định phi tham số / ANOVA, Independent-Sample T-Test.

Implementation và kết quả

Development process (Quy trình xử lý dữ liệu)

Toàn bộ dữ liệu sơ cấp được làm sạch, mã hóa và phân tích thông qua bộ công cụ IBM SPSS Statistics 22.0 kết hợp các thuật toán thống kê tự động hóa.

# Kịch bản phân tích hồi quy và kiểm định đa cộng tuyến (Tham chiếu thuật toán SPSS)
import pandas as pd
import numpy as np
import statsmodels.api as sm
from statsmodels.stats.outliers_influence import variance_inflation_factor

def run_hrm_econometric_analysis(dataset_path: str):
    # 1. Tải và xử lý cấu trúc dữ liệu khảo sát (N = 180)
    df = pd.read_csv(dataset_path)
    feature_cols = ['DDCV', 'DKLV', 'LTPL', 'DTTT', 'PCLD', 'DN', 'VHCT']
    target_col = 'DLLV'
    
    X = df[feature_cols]
    y = df[target_col]
    X_with_const = sm.add_constant(X)
    
    # 2. Ước lượng mô hình hồi quy tuyến tính OLS
    ols_model = sm.OLS(y, X_with_const).fit()
    
    # 3. Kiểm định hiện tượng đa cộng tuyến (VIF Calculation)
    vif_matrix = pd.DataFrame()
    vif_matrix["Factor"] = feature_cols
    vif_matrix["VIF"] = [
        variance_inflation_factor(X.values, i) for i in range(X.shape[1])
    ]
    
    # 4. Trích xuất chỉ số kiểm định mô hình
    print(f"R-squared: {ols_model.rsquared:.4f}")
    print(f"Adjusted R-squared: {ols_model.rsquared_adj:.4f}")
    print(f"F-statistic p-value: {ols_model.f_pvalue:.4e}")
    
    return ols_model.summary(), vif_matrix

* Kịch bản SPSS Syntax chạy kiểm định Cronbach's Alpha và EFA;
RELIABILITY
  /VARIABLES=LTPL1 LTPL2 LTPL3 LTPL4 LTPL5
  /SCALE('Thang do Luong Thuong Phuc Loi') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE
  /SUMMARY=TOTAL.

FACTOR
  /VARIABLES=DDCV1 TO VHCT5
  /MISSING MEANSUB
  /ANALYSIS DDCV1 TO VHCT5
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORE
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PC
  /ROTATION VARIMAX.

Testing và validation

1. Kiểm định độ tin cậy thang đo Cronbach's Alpha

Tất cả 7 nhóm nhân tố độc lập và 1 nhóm nhân tố phụ thuộc đều vượt ngưỡng tiêu chuẩn khắt khe ($\alpha > 0.60$, Corrected Item-Total Correlation $> 0.30$).

Nhóm nhân tố Ký hiệu Số biến quan sát Cronbach's Alpha Đánh giá chất lượng
Lương, thưởng & phúc lợi $LTPL$ 5 0.842 Thang đo rất tốt
Điều kiện làm việc $DKLV$ 5 0.789 Thang đo tốt
Đào tạo & thăng tiến $DTTT$ 4 0.765 Thang đo tốt
Phong cách lãnh đạo $PCLD$ 5 0.812 Thang đo rất tốt
Đặc điểm công việc $DDCV$ 4 0.734 Thang đo tốt
Quan hệ đồng nghiệp $DN$ 4 0.751 Thang đo tốt
Văn hóa công ty $VHCT$ 4 0.728 Thang đo tốt
Động lực làm việc (Phụ thuộc) $DLLV$ 3 0.805 Thang đo rất tốt

2. Phân tích nhân tố khám phá (EFA)

  • Hệ số KMO (Kaiser-Meyer-Olkin) $= 0.816 > 0.5$, khẳng định phân tích nhân tố là hoàn toàn thích hợp.
  • Kiểm định Bartlett's Test of Sphericity có mức ý nghĩa $Sig. = 0.000 < 0.05$, chứng minh các biến có tương quan tuyến tính chặt chẽ trong tổng thể.
  • Tổng phương sai trích (Cumulative Variance Explained) đạt 62.45% ($> 50%$), giá trị Eigenvalue tại nhân tố thứ 7 đạt $1.218 > 1$. Tất cả hệ số tải nhân tố (Factor Loading) đều đạt giá trị $\ge 0.55$.

Kết quả đạt được từ mô hình hồi quy

Kết quả ước lượng tham số mô hình OLS cho thấy giá trị $R^2 = 0.584$ và $R^2_{\text{Adj}} = 0.567$. Điều này chứng minh rằng 56.7% sự biến thiên của Động lực làm việc tại Vinatex Phú Hưng được giải thích bởi 7 nhân tố trong mô hình nghiên cứu.

Nhân tố độc lập Hệ số chưa chuẩn hóa ($B$) Sai số chuẩn ($Std. Error$) Hệ số chuẩn hóa ($\beta$) Mức ý nghĩa ($Sig.$) Hệ số phóng đại phương sai ($VIF$)
Hằng số ($Constant$) 0.412 0.185 - 0.027 -
Lương, thưởng & phúc lợi ($LTPL$) 0.315 0.045 0.342 0.000 1.341
Điều kiện làm việc ($DKLV$) 0.228 0.042 0.251 0.000 1.285
Phong cách lãnh đạo ($PCLD$) 0.174 0.039 0.188 0.001 1.412
Đào tạo & thăng tiến ($DTTT$) 0.142 0.038 0.155 0.003 1.298
Quan hệ đồng nghiệp ($DN$) 0.116 0.036 0.124 0.012 1.210
Đặc điểm công việc ($DDCV$) 0.098 0.035 0.105 0.035 1.189
Văn hóa công ty ($VHCT$) 0.082 0.034 0.089 0.048 1.156

Kiểm định hiện tượng đa cộng tuyến: Hệ số $VIF$ của tất cả các biến đều nằm trong khoảng $1.156 - 1.412 < 2.0$, chứng minh mô hình không vi phạm giả định đa cộng tuyến.

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

  1. Định lượng hóa đặc thù nhà xưởng sợi: Khác với các nghiên cứu tổng quát trong ngành dịch vụ hay hành chính, nghiên cứu này chứng minh thực nghiệm rằng trong môi trường công nghiệp nặng nhọc (sản xuất sợi), Điều kiện làm việc ($\beta = 0.251$) vươn lên thành nhân tố có sức ảnh hưởng lớn thứ hai, chỉ đứng sau Lương, thưởng & phúc lợi ($\beta = 0.342$).
  2. Thiết lập ma trận ưu tiên can thiệp nhân sự: Xác định thứ bậc tác động chuẩn hóa: $$\text{LTPL } (0.342) > \text{DKLV } (0.251) > \text{PCLD } (0.188) > \text{DTTT } (0.155) > \text{DN } (0.124) > \text{DDCV } (0.105) > \text{VHCT } (0.089)$$
  3. Phát hiện khác biệt nhân khẩu học qua One-Way ANOVA & T-Test: Kết quả kiểm định Levene ($p > 0.05$) và ANOVA cho thấy có sự khác biệt có ý nghĩa thống kê về mức độ thỏa mãn động lực làm việc theo nhóm thâm niên và thu nhập ($p < 0.05$), trong khi yếu tố giới tính không tạo ra sự phân hóa đáng kể.

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

Dựa trên kết quả mô hình định lượng, lộ trình can thiệp quản trị tại Công ty Cổ phần Vinatex Phú Hưng được hoạch định theo 3 gói giải pháp chiến lược:

+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH TRIỂN KHAI CAN THIỆP HRM                         |
+-----------------------------------------------------------------------------------+
  Giai đoạn 1 (Tháng 1 - 3): TỐI ƯU HÓA ĐIỀU KIỆN VẬT CHẤT & ĐÃI NGỘ
  - Cải tạo hệ thống lọc bụi xơ, giảm 15% tiếng ồn tại xưởng sợi con
  - Tái cấu trúc quy chế thưởng theo sản lượng (KPI-based bonus)
  ---------------------------------------------------------------------------------
  Giai đoạn 2 (Tháng 4 - 6): CHUẨN HÓA LÃNH ĐẠO & ĐÀO TẠO TẠI CHỖ
  - Tập huấn kỹ năng quản lý dân chủ cho quản đốc, tổ trưởng ca
  - Xây dựng lộ trình thăng tiến rõ ràng từ công nhân bậc 1 lên thợ bậc cao
  ---------------------------------------------------------------------------------
  Giai đoạn 3 (Tháng 7 - 12): GẮN KẾT VĂN HÓA VÀ ĐO LƯỜNG ĐỊNH KỲ
  - Thiết lập kênh tiếp nhận sáng kiến kỹ thuật Kaizen
  - Khảo sát lặp định kỳ (Pulse Survey) theo bộ chỉ số 34 biến quan sát

Dự toán tác động kinh tế và ROI

  • Giảm tỷ lệ luân chuyển lao động: Dự kiến giảm tỷ lệ thôi việc từ 18% xuống dưới 10%/năm, tiết kiệm ước tính 250 - 300 triệu VNĐ chi phí tuyển dụng và đào tạo lại mỗi năm.
  • Gia tăng năng suất chuyền: Tác động tích cực từ gói cải tiến môi trường làm việc và tiền lương giúp tăng sản lượng sợi bông thêm 4.5 - 6.0%/ca máy.

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

  • Hạn chế mẫu: Nghiên cứu áp dụng phương pháp chọn mẫu thuận tiện tại một thời điểm (cross-sectional design) tại một doanh nghiệp đơn lẻ ở Thừa Thiên Huế ($N=180$), hạn chế khả năng tổng quát hóa cho toàn ngành dệt may cả nước.
  • Mở rộng mô hình: Nghiên cứu tiếp theo nên ứng dụng mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) để phân tích các biến trung gian như Sự cam kết gắn bó (Organizational Commitment) hoặc Sự hài lòng công việc (Job Satisfaction).

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

  • Sinh viên & Học viên cao học: Cung cấp mẫu nghiên cứu thực nghiệm hoàn chỉnh từ khung lý thuyết, công thức chọn mẫu Yamane, mã lệnh phân tích định lượng SPSS/Python đến diễn giải kết quả.
  • Giám đốc Nhân sự & Quản đốc nhà máy: Bộ chỉ số 34 biến quan sát thực tế để kiểm toán sức khỏe nhân sự tại các nhà máy dệt may.
  • Nhà nghiên cứu HRM: Tài liệu tham khảo có độ tin cậy cao về hành vi lao động công nghiệp tại khu vực miền Trung Việt Nam.

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

  1. Tại sao nhân tố Điều kiện làm việc lại có trọng số tác động cao ($\beta = 0.251$)? Do ngành sản xuất sợi có tính chất vận hành máy liên tục, bụi bông và tiếng ồn cao. Việc đảm bảo thông gió, bảo hộ lao động và chế độ ca kíp ảnh hưởng trực tiếp đến sức khỏe thể chất của công nhân.
  2. Cỡ mẫu 180 có đảm bảo tính đại diện cho công ty 225 người không? Hoàn toàn đảm bảo. Theo công thức Taro Yamane với sai số $e=5%$, cỡ mẫu yêu cầu là 144. Mẫu 180 chiếm tới 80% toàn bộ lực lượng lao động của nhà máy.
  3. Mô hình có gặp hiện tượng tự tương quan hoặc đa cộng tuyến không? Không. Kiểm định $VIF$ tối đa chỉ đạt 1.412 ($< 2.0$), giá trị $Sig.$ kiểm định F đều đạt $0.000$, đảm bảo tính độc lập và vững của các ước lượng OLS.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Như Ý đã xây dựng thành công khung phân tích định lượng chuẩn mực về các nhân tố ảnh hưởng đến động lực làm việc của người lao động tại Công ty Cổ phần Vinatex Phú Hưng. Với mô hình giải thích được 56.7% phương sai biến thiên, đề tài khẳng định vai trò quyết định của chính sách Lương thưởng - Phúc lợi ($\beta = 0.342$) và Điều kiện làm việc ($\beta = 0.251$), cung cấp luận cứ khoa học vững chắc để ban lãnh đạo doanh nghiệp tối ưu hóa năng suất và phát triển bền vững nguồn nhân lực dệt may.