Giới thiệu dự án

Trong bối cảnh chuyển dịch mạnh mẽ sang nền kinh tế tri thức và công nghiệp hóa 4.0, nguồn nhân lực kỹ thuật cao trở thành yếu tố cốt lõi quyết định lợi thế cạnh tranh của các doanh nghiệp sản xuất và xây dựng hạ tầng. Theo số liệu thống kê ngành xây dựng - vật liệu tại khu vực miền Trung, mức độ luân chuyển nhân sự trung bình dao động từ 15% - 22%/năm, gây thất thoát chi phí tuyển dụng và đào tạo đáng kể. Công ty Cổ phần Phương Minh (FUMICO) – doanh nghiệp chuyên sản xuất cọc bê tông ly tâm, cơ khí đúc và xây lắp điện với quy mô tổng tài sản đạt 92,88 tỷ đồng và 383 lao động (năm 2019) – đang đối mặt với bài toán tối ưu hóa năng suất lao động trong bối cảnh tỷ lệ công nhân kỹ thuật chiếm tới 83,60% và lao động trực tiếp chiếm 76,11%.

Vấn đề then chốt đặt ra (Problem Statement) là doanh nghiệp chưa có công cụ định lượng khoa học để đo lường mức độ tác động của các chính sách quản trị nhân sự đến động lực làm việc thực tế của người lao động. Các quyết định về đãi ngộ, cải thiện môi trường làm việc hay kế hoạch đào tạo còn mang tính cảm tính, chưa dựa trên bằng chứng dữ liệu (data-driven decision making). Điều này dẫn tới nguy cơ suy giảm hiệu suất và gia tăng tỷ lệ bất mãn tiềm ẩn.

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Hệ thống hóa cơ sở lý luận về động lực và các mô hình hành vi tổ chức ứng dụng trong doanh nghiệp sản xuất - xây dựng.
  2. Thiết kế và kiểm định mô hình định lượng 5 nhóm yếu tố tác động đến động lực lao động thông qua thang đo Likert 5 mức độ với cỡ mẫu $n = 150$.
  3. Sử dụng các thuật toán phân tích thống kê (Cronbach's Alpha, EFA, OLS Multiple Linear Regression) để xác định trọng số ảnh hưởng của từng nhân tố.
  4. Đề xuất hệ thống giải pháp quản trị nhân sự (HRM) có khả năng can thiệp chính xác vào các yếu tố có hệ số hồi quy chuẩn hóa ($\beta$) cao nhất.

Giải pháp lựa chọn là xây dựng quy trình phân tích định lượng thực nghiệm (Empirical Quantitative Pipeline) kết hợp giữa lý thuyết tháp nhu cầu Maslow (1943), thuyết kỳ vọng Vroom (1964), thuyết công bằng Stacy Adams (1965) và thang đo JDI (Job Descriptive Index). Phạm vi nghiên cứu tập trung vào toàn bộ cán bộ, nhân viên, công nhân tại trụ sở và các phân xưởng sản xuất của Công ty Cổ phần Phương Minh (TX. Hương Trà, tỉnh Thừa Thiên Huế), với dữ liệu thứ cấp thu thập giai đoạn 2017–2019 và dữ liệu sơ cấp khảo sát từ tháng 10/2020 đến tháng 01/2021.

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

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

Trước khi nghiên cứu được triển khai, công tác đánh giá nhân sự tại Công ty Cổ phần Phương Minh chủ yếu dựa vào phản ánh định tính định kỳ, thiếu các chỉ số đo lường hiệu suất (KPIs) và thang đo tâm lý - hành vi chuẩn hóa.

Tiêu chí đánh giá Phương pháp đánh giá truyền thống Mô hình phân tích định lượng đề xuất (EFA & OLS)
Cơ sở dữ liệu Biên bản họp, nhận xét cảm tính cấp trên Bộ dữ liệu sơ cấp chuẩn hóa ($n=150$) qua thang đo Likert 5 điểm
Độ tin cậy Thấp, dễ phát sinh thiên lệch chủ quan Cao, kiểm định nghiêm ngặt qua Cronbach's Alpha ($\alpha \ge 0.6$)
Khả năng phân rã nhân tố Không phân tách được mức độ ảnh hưởng Tách bạch 5 nhóm độc lập bằng phép xoay Varimax trong EFA
Dự báo tác động Không thể lượng hóa Dự báo chính xác qua hệ số hồi quy $\beta$ và hệ số xác định $R^2$
Chi phí thực thi Thấp nhưng hiệu quả can thiệp kém Tối ưu hóa ngân sách HRM vào đúng nhân tố trọng yếu

Yêu cầu hệ thống đo lường được chuẩn hóa theo mô hình MoSCoW:

  • Must have: Kiểm định độ tin cậy thang đo ($\alpha \ge 0.6$, Corrected Item-Total Correlation $\ge 0.3$), kiểm định EFA ($KMO \ge 0.5$, Bartlett's Sig. $< 0.05$, Factor Loading $\ge 0.5$, Total Variance Explained $> 50%$).
  • Should have: Mô hình hồi quy tuyến tính OLS không vi phạm hiện tượng đa cộng tuyến ($VIF < 2$), tự tương quan ($1 < Durbin-Watson < 3$) và phương sai sai số thay đổi.
  • Could have: Phân tích sâu theo nhóm nhân khẩu học qua One-Sample T-test và Independent Samples T-test.
  • Won't have (this release): Mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) cho các biến trễ đa cấp.

Thiết kế hệ thống

Kiến trúc mô hình nghiên cứu bao gồm 5 biến độc lập cấu thành tác động lên 1 biến phụ thuộc:

Công nghệ và công cụ phân tích dữ liệu:

  • IBM SPSS Statistics 20.0: Xử lý dữ liệu định lượng, Cronbach's Alpha, EFA, OLS Regression, ANOVA, T-Test.
  • Microsoft Excel 2019: Làm sạch, tiền xử lý và mã hóa ma trận dữ liệu ($150 \times 23$).
  • Thang đo Likert 5 mức độ: $1 = \text{Rất không đồng ý}$ đến $5 = \text{Rất đồng ý}$.

Mô hình toán học hồi quy tuyến tính bội:

$$Y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \beta_3 X_3 + \beta_4 X_4 + \beta_5 X_5 + \epsilon$$

Trong đó:

  • $Y$: Động lực làm việc của người lao động (DLLV, 3 biến quan sát).
  • $X_1$: Điều kiện làm việc (DKLV, 4 biến quan sát: DKLV1 - DKLV4).
  • $X_2$: Đặc điểm công việc (DDCV, 4 biến quan sát: DDCV1 - DDCV4).
  • $X_3$: Lương, thưởng và phúc lợi (LTPL, 4 biến quan sát: LTPL1 - LTPL4).
  • $X_4$: Quan hệ với đồng nghiệp, cấp trên (QH, 4 biến quan sát: QH1 - QH4).
  • $X_5$: Đào tạo và thăng tiến (DTTT, 4 biến quan sát: DTTT1 - DTTT4).
  • $\beta_0$: Hằng số hồi quy; $\beta_i$: Hệ số hồi quy riêng phần; $\epsilon$: Sai số ngẫu nhiên.

Methodology

Nghiên cứu áp dụng quy trình 4 giai đoạn kết hợp phương pháp nghiên cứu hỗn hợp (Mixed Methods):

  1. Nghiên cứu sơ bộ định tính: Phỏng vấn sâu 10 chuyên gia và cán bộ quản lý phân xưởng nhằm điều chỉnh thang đo phù hợp với đặc thù sản xuất cọc bê tông ly tâm và xây lắp điện.
  2. Xác định cỡ mẫu: Áp dụng nguyên lý Hair et al. (1998) yêu cầu cỡ mẫu tối thiểu gấp 5 lần số biến quan sát ($n \ge 5 \times 20 = 100$) và Tabachnick & Fidell (2001) cho hồi quy tuyến tính ($n \ge 8m + 50 = 8 \times 5 + 50 = 90$). Cỡ mẫu thực tế được ấn định $N = 150$ phát ra, thu về 150 phiếu hợp lệ (tỷ lệ phản hồi 100%).
  3. Phương pháp chọn mẫu: Chọn mẫu thuận tiện phi xác suất (Convenience Sampling) kết hợp phân tầng theo các phân xưởng (Bê tông BTFMC, Xây lắp XLFMC, Cơ khí XCK, Đúc kim loại XĐKL, Đội ép cọc ECVC và khối Văn phòng).
  4. Kiểm soát rủi ro: Rủi ro sai lệch dữ liệu do tâm lý e ngại được triệt tiêu bằng cơ chế khảo sát ẩn danh (Anonymous Coding) và cam kết không sử dụng dữ liệu cho mục đích kỷ luật lao động.

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu được thực hiện tự động hóa thông qua cú pháp SPSS Syntax:

* 1. Kiem dinh do tin cay Cronbach's Alpha cho nhom Luong thuong va Phuc loi.
RELIABILITY
  /VARIABLES=LTPL1 LTPL2 LTPL3 LTPL4
  /SCALE('Luong Thuong Phuc Loi') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

* 2. Phan tich Nhan to Kham pha (EFA) cho cac bien doc lap.
FACTOR
  /VARIABLES DKLV1 DKLV2 DKLV3 DKLV4 DDCV1 DDCV2 DDCV3 DDCV4 
             LTPL1 LTPL2 LTPL3 LTPL4 QH1 QH2 QH3 QH4 DTTT1 DTTT2 DTTT3 DTTT4
  /MISSING LISTWISE 
  /ANALYSIS DKLV1 TO DTTT4
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORE
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PAF
  /ROTATION VARIMAX
  /METHOD=CORRELATION.

* 3. Phan tich Hoi quy Tuyen tinh Da bien (OLS Regression).
REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA COLLIN TOL
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN 
  /DEPENDENT DLLV
  /METHOD=ENTER DKLV DDCV LTPL QH DTTT
  /RESIDUALS DURBIN.

Testing và validation

Kết quả kiểm định thống kê trên mẫu khảo sát $n = 150$ chứng minh tính hợp lệ và độ tin cậy vượt chuẩn học thuật:

Thang đo / Khái niệm Số biến quan sát Cronbach's Alpha ($\alpha$) Hệ số tương quan biến - tổng nhỏ nhất Đánh giá độ tin cậy
Điều kiện làm việc (DKLV) 4 0.812 0.542 Tốt
Đặc điểm công việc (DDCV) 4 0.789 0.498 Chấp nhận được
Lương thưởng, phúc lợi (LTPL) 4 0.845 0.611 Tốt
Quan hệ đồng nghiệp, cấp trên (QH) 4 0.824 0.573 Tốt
Đào tạo và thăng tiến (DTTT) 4 0.867 0.638 Rất tốt
Động lực làm việc (DLLV - Phụ thuộc) 3 0.831 0.590 Tốt

Kết quả EFA các biến độc lập:

  • Hệ số Kaiser-Meyer-Olkin (KMO) đạt 0.828 ($0.5 < KMO < 1.0$).
  • Kiểm định Bartlett's Test: $\chi^2 = 1428.56$, $p\text{-value} = 0.000 < 0.05$, chứng minh các biến có tương quan tuyến tính chặt chẽ.
  • Rút trích được 5 nhân tố tại mức Eigenvalue = 1.214 (> 1.0).
  • Tổng phương sai trích (Cumulative Variance Explained) đạt 62.48% (> 50%), chứng minh 5 nhân tố giải thích được 62.48% độ biến thiên của tập dữ liệu.
  • Tất cả các hệ số tải nhân tố (Factor Loading) đều đạt từ 0.582 đến 0.841 (> 0.50), không có biến nào bị loại.

Kết quả đạt được

Kết quả ước lượng mô hình hồi quy đa biến (OLS Regression):

Biến độc lập Hệ số chưa chuẩn hóa ($B$) Sai số chuẩn ($SE$) Hệ số chuẩn hóa ($\beta$) Giá trị $t$ Mức ý nghĩa ($Sig.$) Hệ số phóng đại phương sai ($VIF$)
(Hằng số - Constant) 0.312 0.185 - 1.686 0.094 -
Điều kiện làm việc (DKLV) 0.176 0.048 0.182 3.667 0.000 1.245
Đặc điểm công việc (DDCV) 0.142 0.051 0.138 2.784 0.006 1.189
Lương, thưởng & phúc lợi (LTPL) 0.324 0.052 0.341 6.231 0.000 1.312
Quan hệ đồng nghiệp/cấp trên (QH) 0.189 0.046 0.198 4.108 0.000 1.201
Đào tạo và thăng tiến (DTTT) 0.228 0.047 0.235 4.851 0.000 1.278

Thống kê mức độ phù hợp mô hình:

  • Hệ số tương quan bội $R = 0.768$.
  • Hệ số xác định $R^2 = 0.590$; Hệ số $R^2$ điều chỉnh (Adjusted $R^2$) = 0.576.
  • Kiểm định ANOVA: $F(5, 144) = 41.439$ với $p = 0.000 < 0.001$, khẳng định mô hình hoàn toàn phù hợp với dữ liệu thực tế.
  • Tất cả $VIF < 1.4$ chứng minh hoàn toàn không xảy ra hiện tượng đa cộng tuyến.

Thứ tự mức độ tác động chuẩn hóa đến động lực làm việc: $$\text{Lương thưởng, phúc lợi } (\beta = 0.341) > \text{Đào tạo, thăng tiến } (\beta = 0.235) > \text{Quan hệ nội bộ } (\beta = 0.198) > \text{Điều kiện làm việc } (\beta = 0.182) > \text{Đặc điểm công việc } (\beta = 0.138)$$

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

  1. Chuẩn hóa thang đo đặc thù ngành sản xuất vật liệu & xây lắp: Khác biệt với mô hình của Boeve (2007) dành cho y tế hay Lưu Thị Bích Ngọc (2013) dành cho khách sạn, nghiên cứu này tích hợp các biến quan sát đặc thù về an toàn lao động công xưởng (DKLV4) và thiết bị bảo hộ hạng nặng phù hợp với cơ chế vận hành trạm trộn bê tông và máy ép cọc rô-bốt.
  2. Khám phá mối tương quan nghịch giữa đòn bẩy tài chính và chính sách đãi ngộ: Báo cáo tài chính 2017–2019 cho thấy nợ ngắn hạn chiếm tới 72,46 tỷ đồng trên tổng nguồn vốn 92,88 tỷ đồng (tỷ lệ nợ ~78%). Nghiên cứu chỉ ra rằng việc chiếm dụng vốn ngắn hạn để trang trải thanh khoản đã trực tiếp tạo áp lực chậm trễ tiền thưởng, làm suy giảm nhân tố có trọng số cao nhất là Lương thưởng ($\beta = 0.341$).
  3. Cung cấp mô hình định lượng cho khối doanh nghiệp SME: Xây dựng khung đánh giá mẫu với độ chuẩn xác $R^2 = 57.6%$, giúp doanh nghiệp vừa và nhỏ tự động hóa việc đo lường sức khỏe nguồn nhân lực định kỳ mà không cần thuê các đơn vị tư vấn đắt đỏ.

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

Kịch bản can thiệp quản trị nguồn nhân lực (HRM Playbook)

Dựa trên trọng số $\beta$, doanh nghiệp triển khai 3 gói can thiệp trọng điểm:

  • Gói 1: Tái cấu trúc chính sách Lương - Thưởng 3P (Position - Person - Performance):
    • Chuyển đổi từ trả lương theo thời gian đơn thuần sang kết hợp đơn giá sản phẩm hoàn thiện cho 291 lao động trực tiếp tại xưởng bê tông và đội ép cọc.
    • Thiết lập quỹ thưởng nóng theo tiến độ công trình và sáng kiến tiết kiệm vật tư định mức (xi măng, cốt thép).
  • Gói 2: Chuẩn hóa lộ trình Đào tạo - Thăng tiến (Career Path Framework):
    • Xây dựng ma trận kỹ năng (Skill Matrix) cho 318 công nhân kỹ thuật. Định kỳ 6 tháng tổ chức thi nâng bậc thợ gắn liền với mức tăng bậc lương tự động.
    • Tài trợ 100% chi phí chứng chỉ an toàn vận hành máy ép cọc thủy lực và đấu nối điện lưới 220kV.
  • Gói 3: Cải thiện Điều kiện An toàn & Môi trường làm việc:
    • Lắp đặt hệ thống hút lọc bụi xi măng tại xưởng bê tông và hệ thống thông gió xưởng đúc kim loại mạ kẽm.
    • Trang bị đầy đủ bảo hộ lao động đạt chuẩn QUACERT/ISO 9001:2008.

Lộ trình triển khai và Phân tích Chi phí - Lợi ích (ROI)

  • Chi phí triển khai ước tính: 280.000.000 VNĐ/năm (ngân sách đào tạo, nâng cấp bảo hộ, quỹ thưởng sáng kiến).
  • Lợi ích định lượng kỳ vọng:
    • Giảm tỷ lệ luân chuyển lao động kỹ thuật từ 18% xuống dưới 8%/năm, tiết kiệm 150.000.000 VNĐ chi phí tuyển dụng - đào tạo mới.
    • Tăng năng suất ép cọc và đúc cột bê tông bình quân 12,5%, tạo thêm giá trị gia tăng tương đương 620.000.000 VNĐ/năm.
    • Chỉ số ROI dự phóng:

$$\text{ROI} = \frac{770.000.000 - 280.000.000}{280.000.000} \times 100% = 175%$$

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

  • Hạn chế kỹ thuật: Phương pháp lấy mẫu thuận tiện ($n=150$) có thể tồn tại sai số đại diện nhất định đối với các đội thi công phân tán ngoài hiện trường xa. Phương pháp nghiên cứu cắt ngang (Cross-sectional) chỉ phản ánh trạng thái tâm lý tại một thời điểm nhất định.
  • Hướng phát triển đề xuất:
    • Mở rộng mô hình nghiên cứu theo chuỗi thời gian (Longitudinal Study) để theo dõi biến thiên động lực sau khi áp dụng chính sách 3P.
    • Nâng cấp mô hình phân tích sang SEM/PLS-SEM trên phần mềm SmartPLS 3.0 để đánh giá các biến trung gian như Sự gắn kết tổ chức (Organizational Commitment) và 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 ngành QTKD/HRM: Tài liệu tham khảo toàn diện về phương pháp luận nghiên cứu định lượng, cấu trúc xử lý SPSS Syntax và cách liên kết dữ liệu tài chính với hành vi nhân sự.
  • Nhà phân tích dữ liệu nhân sự (HR Data Analysts): Cung cấp khung biến quan sát chuẩn hóa Likert-5 và các ngưỡng kiểm định học thuật nghiêm ngặt.
  • Ban Lãnh đạo & Giám đốc Nhân sự các công ty sản xuất - xây dựng: Bản hướng dẫn chiến lược giúp phân bổ nguồn lực đãi ngộ chính xác dựa trên trọng số thực nghiệm, tối ưu hóa chi phí vận hành.
  • Các nhà nghiên cứu ứng dụng: Bộ khung mô hình mở có thể dễ dàng tái lập (reproducible) cho các ngành công nghiệp phụ trợ khác tại Việt Nam.

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

  1. Yêu cầu phần cứng và phần mềm tối thiểu để tái lập nghiên cứu này là gì?
    Hệ thống máy tính tối thiểu CPU Intel Core i3, 4GB RAM, cài đặt hệ điều hành Windows 7/10/11 hoặc macOS và phần mềm IBM SPSS Statistics phiên bản 20.0 trở lên cùng Microsoft Excel.

  2. Tại sao yếu tố Lương thưởng có hệ số tác động cao nhất ($\beta = 0.341$)?
    Do cơ cấu lao động của công ty có tới 83,6% là công nhân kỹ thuật và 76,11% lao động trực tiếp. Nhóm đối tượng này có nhu cầu sinh học và an toàn (theo tháp Maslow) rất lớn, phụ thuộc trực tiếp vào thu nhập thực tế hàng tháng để trang trải chi phí sinh hoạt gia đình.

  3. Mô hình có kiểm soát hiện tượng vi phạm giả định hồi quy không?
    Có. Mô hình đã kiểm soát đa cộng tuyến qua hệ số phóng đại phương sai ($VIF < 1.4$), tự tương quan qua hệ số Durbin-Watson, và kiểm định phương sai sai số đồng nhất đảm bảo tính ước lượng không chệch tốt nhất (BLUE).

  4. Doanh nghiệp quy mô dưới 100 nhân sự có thể áp dụng thang đo này không?
    Hoàn toàn có thể. Tuy nhiên cần điều chỉnh kích thước mẫu khảo sát tối thiểu $n \ge 50 - 80$ và chú ý điều chỉnh tiêu chuẩn Factor Loading trong EFA lên mức $\ge 0.70$ để đảm bảo ý nghĩa thực tiễn.

  5. Chi phí định kỳ để duy trì hệ thống đánh giá động lực định lượng này là bao nhiêu?
    Gần như bằng 0 đồng chi phí mua bản quyền mới nếu doanh nghiệp tự vận hành biểu mẫu khảo sát trực tuyến qua Google Forms kết hợp các script phân tích tự động bằng ngôn ngữ R hoặc Python (thư viện pandas, statsmodels, factor_analyzer).

Kết luận

Dự án nghiên cứu đã hoàn thành toàn diện việc xây dựng và kiểm định mô hình định lượng các yếu tố ảnh hưởng đến động lực làm việc tại Công ty Cổ phần Phương Minh. Bằng việc xử lý nghiêm ngặt tập dữ liệu $n = 150$ trên phần mềm SPSS 20.0, nghiên cứu đã chứng minh 5 nhân tố giải thích được 57,6% sự biến thiên của động lực làm việc, đồng thời xác lập vị trí ưu tiên: Lương thưởng ($\beta = 0.341$) và Đào tạo thăng tiến ($\beta = 0.235$). Đây là luận cứ khoa học thực nghiệm vững chắc giúp ban lãnh đạo chuyển đổi phương thức quản trị nhân sự từ cảm tính sang mô hình tối ưu hóa dựa trên dữ liệu, đảm bảo sự phát triển bền vững của doanh nghiệp trong giai đoạn mới. Các doanh nghiệp và nhà nghiên cứu có thể ứng dụng ngay bộ thang đo và quy trình phân tích này để nâng cao hiệu suất nguồn nhân lực trong tổ chức của mình.