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, ngành Dệt may Việt Nam đối mặt với sự cạnh tranh gay gắt về nguồn nhân lực. Theo thống kê của Hiệp hội Dệt May Việt Nam (VITAS), tỷ lệ biến động nhân sự (turnover rate) trong ngành dao động trung bình từ 15% đến 25%/năm, gây thất thoát nghiêm trọng về chi phí tuyển dụng, đào tạo lại và gián đoạn dây chuyền công nghệ.

Công ty Cổ phần Đầu tư Dệt may Impulse Fashion Thiên An Phúc (IF THIENANPHUCGATEX) tại Cụm Công nghiệp An Hòa, TP. Huế sở hữu diện tích mặt bằng $22.800\text{ m}^2$, trang bị 14 chuyền may khép kín với quy mô nhân sự khoảng 600 lao động. Tuy nhiên, sau tái cấu trúc tổ chức và ảnh hưởng dịch tễ giai đoạn 2019–2021, doanh nghiệp đối mặt với bài toán nan giải: tần suất công nhân nghỉ việc gia tăng đột biến, nhiều vị trí chuyền may bị bỏ trống máy móc, gây lãng phí năng lực sản xuất định mức và giảm hiệu suất giao hàng xuất khẩu sang các thị trường quốc tế (Hàn Quốc, Nhật Bản, Đức, Pháp, Mỹ).

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

  1. Hệ thống hóa cơ sở lý luận: Tổng hợp các học thuyết quản trị nhân sự kinh điển (Tháp nhu cầu Maslow, Thuyết hai nhân tố Herzberg, Thuyết ERG Alderfer, Thuyết kỳ vọng Vroom, Thuyết công bằng Adams, Thuyết X-Y McGregor) để xây dựng khung phân tích chuẩn hóa cho ngành dệt may.
  2. Đánh giá thực trạng động lực lao động: Khảo sát và đo lường mức độ ảnh hưởng của các điều kiện làm việc, tiền lương, chính sách phúc lợi và phong cách lãnh đạo tại IF Thiên An Phúc.
  3. Mô hình hóa định lượng: Kiểm định độ tin cậy thang đo, phân tích nhân tố khám phá (EFA) và hồi quy tuyến tính bội để xác định chiều hướng, trọng số tác động ($\beta$) của từng nhóm nhân tố đến động lực làm việc ($Y$).
  4. Đề xuất giải pháp can thiệp: Thiết lập lộ trình tối ưu hóa chính sách đãi ngộ, cải thiện môi trường lao động và nâng cao năng suất tổng thể của 14 dây chuyền may.

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

  • Không gian: Toàn bộ các bộ phận chức năng, kho nguyên phụ liệu, xưởng cắt, xưởng may và hoàn thiện tại IF Thiên An Phúc.
  • Thời gian: Dữ liệu thứ cấp tổng hợp chu kỳ 2019–2021; dữ liệu sơ cấp điều tra từ tháng 10/2021 đến tháng 12/2021.
  • Giới hạn: Phương pháp lấy mẫu thuận tiện trong điều kiện kiểm soát dịch bệnh, tập trung vào lực lượng lao động trực tiếp và gián tiếp tại nhà máy.

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

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

Các mô hình đo lường động lực làm việc truyền thống thường áp dụng nguyên bản các khuôn mẫu từ phương Tây hoặc các ngành dịch vụ, dẫn đến sự sai lệch khi triển khai trong môi trường gia công may mặc thâm dụng lao động.

Tiêu chí Mô hình Kovach (1987) Mô hình Herzberg (1959) Mô hình Đề xuất cho IF Thiên An Phúc
Số lượng nhân tố 10 nhân tố rời rạc 2 nhóm nhân tố (Duy trì & Động viên) 8 nhân tố tích hợp đặc thù dệt may
Yếu tố đặc thù Thiếu liên kết nhóm & phúc lợi Tách biệt cứng nhắc hài lòng / bất mãn Tích hợp Văn hóa DN & Cân bằng cuộc sống
Môi trường ứng dụng Công nghiệp tổng hợp tại Mỹ Kỹ sư & Kế toán viên Nhà máy dệt may chuỗi cung ứng FOB/CMT
Khả năng lượng hóa Định tính/Xếp hạng cơ bản Phỏng vấn tình huống Hồi quy tuyến tính OLS & Kiểm định sai số

Ma trận ưu tiên yêu cầu giải pháp quản trị (MoSCoW)

  • Must have (Bắt buộc): Điều chỉnh chính sách đơn giá lương sản phẩm minh bạch; tái cơ cấu điều kiện chiếu sáng và thông gió xưởng may; tối ưu hóa quy trình đánh giá khen thưởng cuối năm.
  • Should have (Nên có): Xây dựng lộ trình đào tạo tay nghề bậc thợ định kỳ 6 tháng; cải thiện chế độ ăn ca tại căn tin $800\text{ m}^2$; hoàn thiện cơ chế đối thoại giữa Quản đốc và Công nhân.
  • Could have (Có thể có): Nâng cấp phần mềm điểm danh tích hợp mã QR tự động hóa; tổ chức hoạt động văn hóa - thể thao nội bộ theo quý.
  • Won't have (Chưa thực hiện): Tự động hóa thay thế hoàn toàn công nhân may bằng robot trong giai đoạn ngắn hạn 2022–2023.

Công nghệ và Phương pháp luận nghiên cứu

Xác định cỡ mẫu theo công thức Taro Yamane (1973)

Với quy mô tổng thể $N = 600$ cán bộ công nhân viên và mức sai số cho phép $e = 0{,}06$ (6%), dung lượng mẫu tối thiểu được xác định:

$$n = \frac{N}{1 + N \cdot e^2} = \frac{600}{1 + 600 \cdot (0{,}06)^2} = \frac{600}{1 + 600 \cdot 0{,}0036} = \frac{600}{3{,}16} \approx 189{,}87$$

Tác giả quyết định phát ra $n = 200$ phiếu khảo sát trực tiếp (thang đo Likert 5 mức độ: 1 = Hoàn toàn không đồng ý đến 5 = Hoàn toàn đồng ý) nhằm bù đắp rủi ro sai sót và tỷ lệ phản hồi không hợp lệ.

Stack công nghệ và công cụ thống kê

  • Phần mềm xử lý dữ liệu: IBM SPSS Statistics v26.0 (Phân tích Cronbach’s Alpha, EFA, Regression, ANOVA) và Microsoft Excel 365 (Làm sạch và mã hóa dữ liệu thô).
  • Công cụ thu thập dữ liệu an toàn: Hệ thống quét mã QR phân luồng vùng dịch tễ kiểm soát an toàn lao động trước khi tiếp cận phỏng vấn.

Implementation và kết quả

Pipeline xử lý và kiểm định dữ liệu

Quy trình xử lý dữ liệu định lượng tuân thủ nghiêm ngặt chuẩn nghiên cứu kinh tế lượng ứng dụng:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from factor_analyzer import FactorAnalyzer
from statsmodels.stats.outliers_influence import variance_inflation_factor

def quantitative_analysis_pipeline(survey_data_path: str):
    """
    Pipeline kiểm định kinh tế lượng cho dữ liệu khảo sát động lực lao động
    """
    # 1. Đọc dữ liệu đã mã hóa
    df = pd.read_csv(survey_data_path)
    independent_vars = ['CV', 'DK', 'TL', 'DT', 'LD', 'DN', 'VH', 'CB']
    dependent_var = 'DL'
    
    # 2. Kiểm định tương quan Pearson & Đa cộng tuyến (VIF)
    X = df[independent_vars]
    X_with_const = sm.add_constant(X)
    
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X_with_const.columns
    vif_data["VIF"] = [variance_inflation_factor(X_with_const.values, i) 
                        for i in range(X_with_const.shape[1])]
    
    # 3. Ước lượng mô hình Hồi quy Tuyến tính Bội (OLS)
    y = df[dependent_var]
    model = sm.OLS(y, X_with_const).fit()
    
    return {
        "summary": model.summary(),
        "vif": vif_data,
        "r_squared": model.rsquared,
        "f_pvalue": model.f_pvalue,
        "params": model.params
    }

Kết quả kiểm định thống kê

==============================================================================
                 KẾT QUẢ KIỂM ĐỊNH ĐỘ TIN CẬY & GIÁ TRỊ THANG ĐO
==============================================================================
1. Cronbach's Alpha:
   - Biến độc lập (34 quan sát): Cronbach's Alpha dao động 0.782 - 0.894 (> 0.70)
   - Biến phụ thuộc (3 quan sát): Cronbach's Alpha = 0.846 (> 0.80)
   - Tương quan biến - tổng (Corrected Item-Total Correlation) > 0.30

2. Phân tích Nhân tố Khám phá (EFA - Principal Axis Factoring, Varimax):
   - Hệ số KMO (Kaiser-Meyer-Olkin) = 0.832 (> 0.50, thích hợp cao)
   - Kiểm định Bartlett's: Chi-Square = 2418.52, Sig. = 0.000 (< 0.05)
   - Tổng phương sai trích (Total Variance Explained) = 63.45% (> 50%)
   - Giá trị Eigenvalue dừng ở mức 1.24 (> 1.0)
   - Tất cả Factor Loading > 0.50 (Không có biến bị loại)

3. Mô hình Hồi quy Tuyến tính Bội (OLS Model):
   - R-squared điều chỉnh (Adjusted R^2) = 0.582
   - Kiểm định F: F = 35.68, Sig. = 0.000 (Mô hình có ý nghĩa thống kê ở độ tin cậy 99%)
   - Hệ số Durbin-Watson = 1.912 (Nằm trong khoảng [1.5, 2.5] -> Không có tự tương quan)
   - Hệ số VIF của tất cả các biến < 2.0 (Không xảy ra hiện tượng đa cộng tuyến)
==============================================================================

Phương trình hồi quy chuẩn hóa

Dựa trên kết quả phân tích Coefficients, phương trình hồi quy tuyến tính biểu diễn mức độ ảnh hưởng đến động lực làm việc ($Y$) được lượng hóa như sau:

$$\hat{Y}{DL} = 0{,}312 \cdot X{TL} + 0{,}245 \cdot X_{LD} + 0{,}198 \cdot X_{DT} + 0{,}167 \cdot X_{CV} + 0{,}142 \cdot X_{DK} + 0{,}115 \cdot X_{CB} + 0{,}094 \cdot X_{DN}$$

(Trong đó: $X_{TL}$: Lương thưởng & Phúc lợi; $X_{LD}$: Phong cách lãnh đạo; $X_{DT}$: Đào tạo & Thăng tiến; $X_{CV}$: Đặc điểm công việc; $X_{DK}$: Điều kiện làm việc; $X_{CB}$: Cân bằng cuộc sống; $X_{DN}$: Mối quan hệ đồng nghiệp; tất cả đều có $p\text{-value} < 0{,}05$).

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

  1. Chuẩn hóa khung đánh giá nhân sự ngành Dệt may miền Trung: Bổ sung hai biến quan sát thực chứng là "Văn hóa doanh nghiệp thích ứng hậu dịch" và "Cân bằng cuộc sống - công việc đối với lao động nữ" (chiếm hơn 75% cơ cấu nhân sự dệt may), khắc phục nhược điểm thiếu tính ngữ cảnh của mô hình Kovach kinh điển.
  2. Quy trình kết hợp định tính - định lượng 2 giai đoạn: Phỏng vấn chuyên sâu 7 cán bộ quản lý cấp cao (Giám đốc sản xuất, Trưởng phòng Nhân sự, Kế hoạch, Kế toán) kết hợp thảo luận nhóm tiêu điểm 10 công nhân để tinh chỉnh bảng hỏi trước khi phân tích SPSS v26.
  3. Định lượng hóa các đòn bẩy phi tài chính: Chứng minh bằng thực nghiệm rằng nhân tố "Phong cách lãnh đạo" ($\beta = 0{,}245$) và "Đào tạo & Thăng tiến" ($\beta = 0{,}198$) có trọng số giải thích biến thiên động lực lên tới 44.3%, bác bỏ quan điểm truyền thống cho rằng công nhân may mặc chỉ làm việc vì thu nhập thuần túy.
Hạng mục so sánh Khảo sát định tính truyền thống Mô hình định lượng chuẩn hóa của Đề tài
Độ chính xác đo lường Cảm tính, phụ thuộc ý chí quản đốc Chỉ số $R^2 = 58{,}2%$, kiểm định F $p < 0{,}001$
Xác thực thang đo Không kiểm tra độ tin cậy Cronbach's Alpha $\ge 0{,}782$, KMO $= 0{,}832$
Phân hóa nhóm đối tượng Đồng nhất chính sách cho toàn xưởng Kiểm định ANOVA phân tách rõ thâm niên & thu nhập
Khả năng ứng dụng Giải pháp chung chung Phân bổ ngân sách theo trọng số $\beta$

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

Lộ trình triển khai 4 giai đoạn tại Doanh nghiệp

Phân tích Hiệu quả Đầu tư (Cost-Benefit & ROI Analysis)

  • Chi phí ước tính triển khai: 180.000.000 VNĐ (Bao gồm đào tạo quản lý cấp trung, cải thiện vi khí hậu xưởng may, quỹ khen thưởng sáng kiến chuyền may).
  • Lợi ích định lượng kỳ vọng:
    • Giảm tỷ lệ nghỉ việc từ $18%$ xuống dưới $8%$/năm (tiết kiệm ước tính ~320.000.000 VNĐ chi phí tuyển dụng, may đồng phục, đào tạo học việc cho 60 công nhân mới).
    • Tăng năng suất chuyền may trung bình $7{,}5%$ nhờ ổn định nhân lực, giảm thời gian ngừng máy (downtime).
    • ROI dự kiến: $\frac{320.000.000 - 180.000.000}{180.000.000} \times 100% = 77{,}7%$ ngay trong năm đầu tiên áp dụng.

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

Hạn chế kỹ thuật

  • Phương pháp chọn mẫu: Do rào cản giãn cách dịch tễ cục bộ, nghiên cứu sử dụng phương pháp lấy mẫu thuận tiện ($n = 200$), chưa áp dụng được phương pháp lấy mẫu ngẫu nhiên phân tầng theo từng chuyền may.
  • Thời gian quan sát: Nghiên cứu theo thiết kế cắt ngang (cross-sectional), chưa thu thập được dữ liệu chuỗi thời gian (longitudinal) để đánh giá biến thiên động lực theo mùa vụ đơn hàng dệt may (mùa cao điểm xuất khẩu Thu - Đông vs mùa thấp điểm).

Hướng nghiên cứu mở rộng

  • Ứ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ự hài lòng trong công việc" và "Lòng trung thành tổ chức".
  • Xây dựng Dashboard theo dõi chỉ số gắn kết nhân viên (eNPS - Employee Net Promoter Score) theo thời gian thực tích hợp với hệ thống ERP quản trị sản xuất dệt may.

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

  • Sinh viên chuyên ngành Quản trị Kinh doanh / Quản trị Nhân lực: Tài liệu tham khảo toàn diện về cách kết hợp lý thuyết động lực học cổ điển với các công cụ kinh tế lượng hiện đại.
  • Nhà quản lý Doanh nghiệp Dệt may: Bản thiết kế giải pháp can thiệp có thứ tự ưu tiên rõ ràng dựa trên hệ số tác động $\beta$, giúp phân bổ ngân sách nhân sự tối ưu.
  • Người lao động trực tiếp: Được làm việc trong môi trường thấu cảm, chế độ đãi ngộ công bằng và có lộ trình phát triển nghề nghiệp rõ ràng.

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

1. Cần chuẩn bị những điều kiện gì để triển khai hệ thống khảo sát định lượng nhân sự tương tự?

Doanh nghiệp cần hoàn thiện cơ sở dữ liệu nhân sự (mã nhân viên, thâm niên, xưởng sản xuất, bậc thợ), thành lập nhóm điều phối trung lập (không có sự giám sát trực tiếp của tổ trưởng khi công nhân điền phiếu) và sử dụng phần mềm phân tích như IBM SPSS Statistics hoặc Python/R để xử lý dữ liệu.

2. Mô hình có áp dụng được cho các nhà máy dệt may có quy mô trên 5.000 công nhân không?

Hoàn toàn khả thi. Khi quy mô tổng thể $N \ge 5.000$, công thức Taro Yamane hoặc Cochran sẽ yêu cầu dung lượng mẫu $n \approx 350 - 400$ quan sát. Khi đó, nên chuyển đổi phương pháp chọn mẫu sang kỹ thuật "Lấy mẫu ngẫu nhiên phân tầng" (Stratified Random Sampling) theo từng xưởng: Xưởng Cắt, Xưởng May, Xưởng Hoàn thiện.

3. Làm thế nào để giải quyết mâu thuẫn giữa nhân tố "Lương thưởng" và "Cân bằng cuộc sống" trong mùa cao điểm tăng ca?

Kết quả kiểm định hồi quy cho thấy cả $X_{TL}$ ($\beta = 0{,}312$) và $X_{CB}$ ($\beta = 0{,}115$) đều tác động cùng chiều. Trong mùa cao điểm xuất khẩu, doanh nghiệp cần áp dụng chế độ phụ cấp lũy tiến theo giờ làm thêm (Overtime incentive), bố trí bữa phụ giàu dinh dưỡng tại căn tin và sắp xếp ngày nghỉ bù luân phiên để tránh hiện tượng kiệt sức lao động (burnout).

4. Tại sao nhân tố "Văn hóa công ty" trong một số nghiên cứu có hệ số tác động thấp hoặc bị loại trong mô hình hồi quy?

Trong các doanh nghiệp sản xuất thuần túy gia công, nhận thức của công nhân trực tiếp về văn hóa doanh nghiệp thường bị đồng hóa với "Phong cách lãnh đạo của Quản đốc" và "Mối quan hệ đồng nghiệp". Khi phân tích EFA, các biến văn hóa thường có xu hướng gộp vào nhóm lãnh đạo hoặc không đạt độ phân biệt nếu truyền thông nội bộ chưa sâu.

5. Thời gian hoàn vốn (Payback Period) của các chương trình can thiệp nhân sự là bao lâu?

Dựa trên tính toán chi phí giữ chân so với chi phí tuyển dụng - đào tạo thay thế tại IF Thiên An Phúc, thời gian hoàn vốn đầu tư dao động từ 6 đến 9 tháng kể từ khi áp dụng đầy đủ quy chế lương thưởng minh bạch và chương trình đào tạo quản lý cấp trung.

Kết luận

Đề tài khóa luận tốt nghiệp "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 Đầu tư Dệt may Impulse Fashion Thiên An Phúc" của tác giả Nguyễn Hoàng Nhật đã giải quyết trọn vẹn bài toán thực tiễn của doanh nghiệp thông qua phương pháp luận nghiên cứu khoa học chuẩn mực.

Bằng việc kết hợp chặt chẽ giữa lý thuyết quản trị và phân tích kinh tế lượng (Cronbach's Alpha, EFA, OLS Regression, ANOVA), nghiên cứu đã chỉ rõ 7 nhóm nhân tố then chốt, trong đó Chính sách lương thưởng - phúc lợi ($\beta = 0{,}312$) và Phong cách lãnh đạo ($\beta = 0{,}245$) giữ vai trò quyết định. Các hàm ý quản trị và lộ trình triển khai đề xuất không chỉ trực tiếp tháo gỡ điểm nghẽn thiếu hụt lao động tại IF Thiên An Phúc mà còn đóng vai trò mô hình tham chiếu giá trị cho cộng đồng doanh nghiệp dệt may trong kỷ nguyên nâng cao năng suất dựa trên nguồn lực con người.